USE OF A GEOGRAPHIC INFORMATION SYSTEM TO ASSEMBLE INPUT-DATA SETS FOR A FINITE-DIFFERENCE MODEL OF GROUND-WATER FLOW1

ABSTRACT: A Geographic Information System (GIS) was used to expedite assembly of input-data sets for a model of ground-water flow in the Middle Patuxent River basin in Howard County, Maryland. The model grid was developed with GIS and used to select attributes from GIS data coverages. These attributes were then output from GIS into the input-data sets, which included model boundaries, regolith thickness, bedrock-surface altitude, stream locations, stream length, stream-bottom altitude, and trans-missivity. The ability to change large sets of spatial data quickly and accurately with GIS enhances the model-calibration process. (KEY TERMS: GIS; ground water; Piedmont; Maryland; MODFLOW.)     

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.     

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.     

CONTINENTAL SCALE SIMULATION OF THE HYDROLOGIC BALANCE1

This paper describes the application of a continuous daily water balance model called SWAT (Soil and Water Assessment Tool) for the conterminous U.S. The local water balance is represented by four control volumes; (1) snow, (2) soil profile, (3) shallow aquifer, and (4) deep aquifer. The components of the water balance are simulated using “storage” models and readily available input parameters. All the required databases (soils, landuse, and topography) were assembled for the conterminous U.S. at 1:250,000 scale. A GIS interface was utilized to automate the assembly of the model input files from map layers and relational databases. The hydrologic balance for each soil association polygon (78,863 nationwide) was simulated without calibration for 20 years using dominant soil and land use properties. The model was validated by comparing simulated average annual runoff with long term average annual runoff from USGS stream gage records. Results indicate over 45 percent of the modeled U.S. are within 50 mm of measured, and 18 percent are within 10 mm without calibration. The model tended to under predict runoff in mountain areas due to lack of climate stations at high elevations. Given the limitations of the study, (i.e., spatial resolution of the data bases and model simplicity), the results show that the large scale hydrologic balance can be realistically simulated using a continuous water balance model.     

ANALYSIS AND MODELING OF LONG-TERM STREAM TEMPERATURES ON THE STEAMBOAT CREEK BASIN,OREGON: IMPLICATIONS FOR LAND USE AND FISH HABITAT1

ABSTRACT: Steamboat Creek basin is an important source of timber and provides crucial spawning and rearing habitat for anadromous steelhead trout (Oncorhynchus mykiss). Because stream temperatures are near the upper limit of tolerance for the survival of juvenile steelhead, the possible long-term effect of clear-cut logging on stream temperatures was assessed. Twenty-year (1969–1989) records of summer stream temperature and flow from four tributaries and two reaches of Steamboat Creek and Boulder Creek (a nearby unlogged watershed) were analyzed. Logging records for the Steamboat Creek basin and air temperature records also were used in the analysis. A time-series model of the components of stream temperature (seasonal cycle of solar radiation, air temperature, streamflow, an autoregressive term of order 1, and a linear trend variable) was fitted to the water-temperature data. The linear trend variable was significant in all the fitted models except Bend Creek (a tributary fed by cool ground-water discharge) and Boulder Creek. Because no trends in either climate (i.e., air temperature) or streamflow were found in the data, the trend variable was associated with the pre-1969 loss and subsequent regrowth of riparian vegetation and shading canopies.     

WATER QUALITY MODEL OUTPUT UNCERTAINTY AS AFFECTED BY SPATIAL RESOLUTION OF INPUT DATA1

ABSTRACT: Resolution of the input GIS data used to parameterize distributed‐parameter hydrologic/water quality models may affect uncertainty in model outputs and impact the subsequent application of model results in watershed management. In this study we evaluated the impact of varying spatial resolutions of DEM, land use, and soil data (30 × 30 m, 100 × 100 m, 150 × 150 m, 200 × 200 m, 300 × 300 m, 500 × 500 m, and 1,000 × 1,000 m) on the uncertainty of SWAT predicted flow, sediment, NO₃‐N, and TP transport. Inputs included measured hydrologic, meteorological, and watershed characteristics as well as water quality data from the Moores Creek watershed in Washington County, Arkansas. The SWAT model output was most affected by input DEM data resolution. A coarser DEM data resolution resulted in decreased representation of watershed area and slope and increased slope length. Distribution of pasture, forest, and urban areas within the watershed was significantly affected at coarser resolution of land use and resulted in significant uncertainty in predicted sediment, NO₃‐N, and TP output. Soils data resolution had no significant effect on flow and NO₃‐N predictions; however, sediment was overpredicted by 26 percent, and TP was underpredicted by 26 percent at 1,000 m resolution. This may be due to change in relative distribution of various hydrologic soils groups (HSGs) in the watershed. Minimum resolution for input GIS data to achieve less than 10 percent model output error depended upon the output variable of interest. For flow, sediment, NO₃‐N, and TP predictions, minimum DEM data resolution should range from 30 to 300 m, whereas minimum land use and soils data resolution should range from 300 to 500 m.     

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.     

DEVELOPMENT OF A WATER SUPPLY PROTECTION MODEL IN A GIS1

ABSTRACT: As part of a larger model to identify lands suitable for acquisition, a water supply protection model was developed using the Southwest Florida Water Management District's GIS. Several hydrologic and hydrogeologic data layers were overlaid to develop maps showing ground-water supply suitability, protection areas for surface-water supply, protection areas for major public supply wells, susceptibility to ground-water contamination, and recharge to the Floridan aquifer. These intermediate layers were combined into a final map to prioritize protection areas for water supply.     

CALCULATION OF NONLINEAR CONFIDENCE AND PREDICTION INTERVALS FOR GROUND-WATER FLOW MODELS1

ABSTRACT: A method is derived to efficiently compute nonlinear confidence and prediction intervals on any function of parameters derived as output from a mathematical model of a physical system. The method is applied to the problem of obtaining confidence and prediction intervals for manually-calibrated ground-water flow models. To obtain confidence and prediction intervals resulting from uncertainties in parameters, the calibrated model and information on extreme ranges and ordering of the model parameters within one or more independent groups are required. If random errors in the dependent variable are present in addition to uncertainties in parameters, then calculation of prediction intervals also requires information on the extreme range of error expected. A simple Monte Carlo method is used to compute the quantiles necessary to establish probability levels for the confidence and prediction intervals. Application of the method to a hypothetical example showed that includsion of random errors in the dependent variable in addition to uncertainties in parameters can considerably widen the prediction intervals.     

NUMERICAL MODEL OF A TRACER TEST ON THE SANTA CLARA RWER,VENTURA COUNTY,CALIFORNIA1

ABSTRACT: To better understand the flow processes, solute-trans. port processes, and ground-water/surface-water interactions on the Santa Clara River in Ventura County, California, a 24-hour fluorescent-dye tracer study was performed under steady-state flow conditions on a 45-km reach of the river. The study reach includes perennial (uppermost and lowermost) subreaches and ephemeral subreaches of the lower Piru Creek and the middle Santa Clara River. The tracer-test data were used to calibrate a one-dimensional flow model (DAFLOW) and a solute-transport model (BLTM). The dye-arrival times at each sample location were simulated by calibrating the velocity parameters in DAFLOW. The simulations of dye transport indicated that (1) ground-water recharge explains the loss of mass in the ephemeral middle subreaches, and (2) ground-water recharge does not explain the loss of mass in the perennial uppermost and lowermost subreaches. The observed tracer curves in the perennial subreaches were indicative of sorptive dye losses, transient storage, and (or) photodecay - these phenomena were simulated using a linear decay term. However, analysis of the linear decay terms indicated that photodecay was not a dominant source of dye loss.     

INTEGRATION OF A BASIN-SCALE WATER QUALITY MODEL WITH GIS1

ABSTRACT: Geographic Information Systems (GIS) have been successfully integrated with distributed parameter, single-event, water quality models such as AGNPS (AGricultural NonPoint Source) and ANSWERS (Areal Nonpoint Source Watershed Environmental Response Simulation). These linkages proved to be an effective way to collect, manipulate, visualize, and analyze the input and output date of water quality models. However, for continuous-time, basin large-scale water quality models, collecting and manipulating the input data are more time-consuming and cumbersome due to the method of disaggregation (subdivisions are based on topographic boundaries). SWAT (Soil and Water Assessment Tool), a basin-scale water quality model, was integrated with a GIS to extract input data for modeling a basin. This paper discusses the detailed development of the integration of the SWAT water quality model with GRASS (Geographic Resources Analysis Support System) GIS, along with an application and advantages. The integrated system was applied to simulated a 114 sq. km upper portion of the Seco Creek Basin by subdividing it into 37 subbasins. The average monthly predicted streamflw is in agreement with measured monthly streamflw values.     

THE DEVELOPMENT AND IMPLEMENTATION OF A GIS-BASED CONTAMINANT SOURCE INVENTORY OVER THE SPOKANE AQUIFER,SPOKANE COUNTY WASHINGTON1

ABSTRACT: U.S. Environmental Protection Agency (EPA) statutes require identification of all potential sources of contamination within a wellhead protection area. All wells over Spokane County's aquifer are included in one wellhead protection zone called the Aquifer Sensitive Area (ASA). A GIS-based Contaminant Source Inventory (CSI) was developed for the ASA. Datasets listing businesses and agencies within the ASA were imported into the GIS from state, county, city, and local agencies. These datasets were selected, joined, and sorted using GIS relational database capabilities into one ASA “business master file.” Map files were projected and transformed into common coordinates. Next, business sites within the master file were spatially related by address to the digital map files. Likely Critical Materials Users (CMU) were identified by sorting on selected standard Industrial Codes (SIC). Additional files of CMUs were imported into the Contaminant Source Inventory. GIS queries were performed to locate specific materials, quantities, and storage facilities, and to analyze CMU activity within selected buffer zones. This project demonstrated the usefulness of GIS technology in the development, management, maintenance, and analysis of vast quantities of data associated with a local wellhead protection pro. gram.     

EXPOSURE ASSESSMENT OF POPULATIONS USING ENVIRONMENTAL MODELING,DEMOGRAPHIC ANALYSIS,AND GIS1

ABSTRACT: Simulation of ground-water flow and fate of contaminants in the subsurface environment constitutes a major phase of most environmental assessment and site remediation studies. These simulation studies yield information on spatial and temporal distributions of contaminants in the subsurface media. An important use of this information is to conduct exposure assessment studies. Spatial and temporal distributions of both chemical concentrations and exposed populations render this integrated exposure analysis task rather difficult. Geographic Information Systems (GIS), on the other hand, provide a platform in which layered, spatially distributed databases can be manipulated with ease, thereby simplifying exposure analysis tasks significantly. In this paper, we describe procedures that combine the simulation models and demographic databases under a GIS platform to automate the exposure assessment phase of a typical health assessment study. Procedures developed herein significantly simplify the post-processing phase of the analysis, and render the overall task more ‘user friendly.’ A site-specific application is included as a demonstration of the proposed process.     

Development and Integration of a Groundwater Simulation Model to an Open Geographic Information System

This study presents the implementation of a spatial decision support system (SDSS) named ARENA. The program has been developed based on object‐oriented concepts using the Java programming language. The SDSS is made up of a groundwater simulation tool coupled to an open geographic information system (open GIS). Both the open GIS and groundwater simulation package share common spatial and nonspatial entities during the modeling process. A dedicated interface provides direct access to the GIS data without the need of external files. The finite element method is used to solve the partial differential equation that governs groundwater flow. The system implementation is presented by depicting the main classes and coupling procedures. A study case demonstrates the applicability of the simulation tool.     

THE DEVELOPMENT OF THE NILE HYDROMETEOROLOGICAL FORECAST SYSTEM1

ABSTRACT: The National Oceanic and Atmospheric Administration is developing a river forecast system for the Nile River in Egypt. The river forecast system operates on scientific work stations using hydrometeorological models and software to predict inflows into the high Aswan Dam and forecast flow hydrographs at selected gaging locations above the dam The Nile Forecasting System (NFS) utilizes satellite imagery from the METEOSAT satellite as the input to the forecast system. Satellite imagery is used to estimate precipitation over the Blue Nile Basin using five different techniques. Observed precipitation data and climatic statistics are used to improve precipitation estimation. Precipitation data for grid locations are input to a distributed water balance model, a hill slope routing model, and a channel routing model. A customized Geographic Information System (GIS) was developed to show political boundaries, rivers, terrain elevation, and gaging network. The GIS was used to develop hydrologic parameters for the basin and is used for multiple display features.     

WATER QUALITY MODELING IN THE UPPER REACHES OF THE CHOWAN RIVER1

ABSTRACT: The Chowan River system consists of three rivers in southeast Virginia that form two confluences before flowing into Albermarle Sound in North Carolina. A computer program was written to simulate flows through the river system to determine flow rates, velocities, and depths. The output of the flow program was input into a second program that calculated the concentrations of BOD₅, COD, DO, and four nitrogen parameters (organic, ammonia, nitrite-nitrate and algal-N). Measured field data were used to calibrate the model. The effect of reducing the concentration of nutrients from overland runoff on algal concentrations at the mouth of the river was studied. The program was also run to simulate the water quality of the watershed in a primitive condition, in which the watershed was assumed to consist only of forests. The results of the computer program indicate that the major changes in the water quality of the river are simulated satisfactorily. The program can be used to assess the impact of any management scheme to improve water quality.     

AN AUTO-CAD-BASED WATERSHED INFORMATION SYSTEM FOR THE HYDROLOGIC MODEL HEC-11

ABSTRACT: A micro computer based Watershed Information System (W.LS.) is developed to assist in the preparation of input files for the hydrologic simulation model HEC-1. This system consists of three phases. Phase I utilizes the capabilities of AutoCAD version 9 and three programs, BASINS, PLANES, and CHANNELS, to extract, organize, and display watershed data. Phase II uses the program CN to calculate some HEC-1 parameter values. Phase II utilizes the program HECUPDATE to create HEC-1 input files. The system input includes topographic, soils, land use, watershed geometry data, and a skeletal HEC-1 input file. Output from the system includes a summary User Reference File, a Soils File, a Land Use File, a Watershed Geometry File, a Curve Number File, and a HEC-1 input file, which is ready to run. The W.I.S. has been applied to Macks Creek Watershed in southwest Idaho.     

SENSITWITY ANALYSIS OF SIMULATED CONTAMINATED SEDIMENT TRANSPORT1

ABSTRACT: A simulation analysis of contaminated sediment transport involves model selection, data collection, model calibration and verification, and evaluation of uncertainty in the results. Sensitivity analyses provide information to address these issues at several stages of the investigation. A sensitivity analysis of simulated contaminated sediment transport is used to identify the most sensitive output variables and the parameters most responsible for the output variable sensitivity. The output variables included are streamflow and the flux of sediment and Cs¹³⁷. The sensitivities of these variables are measured at the field and intermediate scales, for flood and normal flow conditions, using the HSPF computer model. A sensitivity index was used to summarize and compare the results of a large number of output variables and parameters. An extensive database was developed to calibrate the model and conduct the sensitivity analysis on a 6.2 mi² catchment in eastern Tennessee. The fluxes of sediment and Cs¹³⁷ were more sensitive than streamflow to changes in parameters for both flood and normal flow conditions. The relative significance of specific parameters on output variable sensitivity varied according to the type of flow condition and the location in the catchment. An implications section illustrates how sensitivity analysis results can help with model selection, planning data collection, calibration, and uncertainty analysis.     

USING A GIS FOR ESTIMATING INPUT PARAMETERS IN URBAN STORMWATER QUALITY MODELING1

ABSTRACT: Geographic Information Systems (GIS) are being used increasingly as a method of preparing, analyzing, and displaying data for watershed analysis and modeling. Although GIS technology is a powerful tool for integrating and analyzing watershed characteristics, the initial preparation of the necessary database is often a time consuming and costly endeavor. This demonstration project assesses the viability of creating a cost-effective spatial database for urban stormwater modeling from existing digital and hard-copy data sources. The GIS was used to provide input parameters to the Source Loading and Management Model (SLANM), an empirical urban stormwater quality model. Land use characteristics, drainage boundaries, and soils information were geocoded and referenced to a base data layer consisting of transportation features. GIS overlay and data manipulation capabilities were utilized to preprocess the input data for the model. Model output was analyzed through postprocessing by GIS, and results were compared to a similar recent modeling study of the same watershed. The project, undertaken for a small urban watershed located in Plymouth, Minnesota, successfully demonstrates that the use of GIS in stormwater management can allow even small communities to reap the benefits of stormwater quality modeling.     

HYDROGEOLOGY AND PATHS OF FLOW IN THE CARBONATE BEDROCK AQUIFER,NORTHWESTERN INDIANA1

ABSTRACT: The U.S. Geological Survey (USGS) is assessing the ground-water resources of the carbonate bedrock aquifers in Indiana and Ohio as part of their Regional Aquifer Systems Analysis program. Part of this assessment includes the determination of unknown aspects of the hydraulic characteristics, boundaries, and flow paths of the carbonate aquifer. To accomplish this, the USGS drilled three wells through the carbonate aquifer near the Kankakee River in northwestern Indiana. Geophysical logs were used to help determine depths and thicknesses for testing and to help describe geology at the three wells. Packer tests were used to determine direction of ground-water flow and to provide data for an analysis of the distribution of transmissivity in the carbonate aquifer. Transmissivity of the carbonates is associated with two physical characteristics of the rocks: fractures and interconnected porosity. Almost all of the transmissivity is derived from horizontal fracturing; however, only a few of the fractures present in the carbonate are transmissive. Some transmissivity is associated with a zone of fossiliferous, vuggy dolomite, which yields water from the rock matrix. Most of the transmissivity is associated with large fractures and solution crevices in the upper 30 feet of the bedrock; less transmissivity is associated with the deeper vuggy reef material, even where extensively fractured. Transmissivity of individual fractures and fossiliferous zones ranges from 300 to 27,000 feet squared per day. The aquifer bottom is defined by a lack of transmissive fractures and an increased shale content near the contact of the Silurian and Ordovician sections. Water-level data from the three wells indicate that flow is horizontal at well site 1 north of the Kankakee River, upward at well site 2 near the river, and downward at well site 3 south of the river. Most of the flow occurs in the upper part of the carbonate bedrock where fracturing and solution-enlarged crevices are most developed. Water levels indicate the the Kankakee River is a hydrologic boundary for the regional carbonate bedrock aquifer.