USGS - science for a changing world
Header graphic spacer
Watershed Modeling
 
Home | People | Research Topics | Publications | Software | Training | Related Links | FAQs

 

 

A Modular Approach to Integrating Environmental Modeling and GIS

 

Table of Contents

Keywords: MMS, GIS Weasel, GIS, Modeling
Pre-process Component
Model Component
Post-process Component
Delineation Phase
Characterization Phase
Modification Phase
Parameterization Phase
References

 

G.H. Leavesley, R.J. Viger, S.M. Markstrom, and M.S. Brewer
US Geological Survey - Water Resources Division
MS 412, Box 25046 DFC, Denver CO, USA, 80225
email: george@usgs.gov
email: rviger@usgs.gov

 

Keywords: MMS, GIS Weasel, GIS, Modeling

ABSTRACT
The Modular Modeling System (MMS) is an integrated system of computer software that is being developed to provide the research and operational framework needed to support development, testing, and evaluation of physical process algorithms and to facilitate integration of user-selected sets of algorithms into operational environmental-process models. A geographic information system (GIS) interface, the GIS Weasel, has been developed to support MMS in model development, application, and analysis. The GIS Weasel permits application of a variety of GIS tools to delineate, characterize, and parameterize the topographic, hydrologic, and biologic features of a physical system for use in a variety of lumped and distributed parameter modeling approaches. The integration of the GIS Weasel and MMS provide a flexible framework in which to integrate and apply environmental models and analytical tools.

INTRODUCTION
The interdisciplinary nature and increasing complexity of environmental problems require the use of a variety of modeling approaches that can be coupled together to adequately address all elements of a given problem. Selection of models and model components to address these problems is difficult given the large number of available models and the potentially wide range of study objectives, data constraints, and spatial and temporal scales of application. To address the problems of model selection, application, and analysis, a set of modular modeling tools, termed the Modular Modeling System (MMS)1, has been developed. MMS uses a module library that contains compatible modules for simulating a variety of hydrologic and ecosystem processes. A model is created by selectively coupling the most appropriate modules from the library to create an "optimal" model, or models, for the desired application. Where existing physical process simulation modules are not appropriate, new modules can be developed and easily added to the system.

In addition to the issues of model creation and selection are the associated issues of study area characterization and the estimation of parameters for these models. Guidelines for characterization and parameterization are normally few and users commonly develop their own methods to analyze available data and derive initial parameter estimates. The subjective nature of many such methods makes it difficult to compare model applications among different users and to effectively transfer methods to other study areas. Geographic information systems (GIS) provide the capability to develop more objective and reproducible methods of characterization and parameter estimation. A GIS tool set, termed the GIS Weasel, has been developed to use a wide variety of methodologies in support of lumped and distributed parameter models. Parameter estimation procedures specific to MMS modules are also provided by the GIS Weasel.

The linking of modeling and GIS tools provides a common framework in which to focus multidisciplinary research and operational efforts to provide a flexible modeling system for developing, analyzing, and applying selected modeling approaches for a wide variety of environmental problems. The purpose of this paper is to provide an overview of MMS, the GIS Weasel, and the integration of these technologies.

MODULAR MODELING SYSTEM
The conceptual framework for MMS has three major components: pre-process, model, and post-process (Fig. 1). All three components retrieve and store data in an underlying data storage scheme. The storage scheme may be composed of one or more data structures ranging from simple flat files to SQL type databases. Data management interfaces (DMIs) handle the transfer and reformatting of information among system components and the data structures. DMIs are essential to the interaction of the three components and are the only elements that are data-structure specific. Thus, changing data storage structures requires only a change of the DMI.

The GIS Weasel supports MMS as a pre-processing component. However, it can also provide GIS post-processing capabilities for the visualization and analysis of spatial and temporal model output fields.

figure1

Figure 1. A schematic diagram showing the integration of the GIS Weasel with the
components of the Modular Modeling System (MMS).
Graphical user interfaces (GUIs) provide access to all the features of MMS and the GIS Weasel. The present framework has been developed for UNIX-based workstations and uses X-windows and Motif for the GUIs. The GUIs provide an interactive environment for users to access model components, apply selected options, and graphically display simulation and analysis results. The current GUIs are being rewritten in the JAVA programming language to enable MMS application to a wider range of computer platforms.

 

Pre-process Component

Pre-process component functions include all data preparation and analysis functions needed to meet the data and selected parameterization requirements of a user-created or user-selected model. Time-series data from existing databases as well as from field instrumentation are prepared for use in selected model applications. The spatial analysis needed to derive model parameters is accomplished using the GIS Weasel.

 

Model Component

The model component is the core of the system and includes the tools to selectively link process modules from the module library together into a model and to interact with this model to perform a variety of simulation and analysis tasks. The library can contain several modules for a given process, each representing an alternative conceptualization for simulating that process. The user, through an interactive model-building interface (Xmbuild), selects and links modules to create a specific model. Once a model has been built, it may be saved for future use without repeating the MBUILD step. This capability allows `canned' models to be provided to users.

 

Post-process Component

Tools to analyze model results are included in the post-processing component. These include a variety of visualization and statistical tools as well as the ability to interface with resource management models and decision support systems (DSS). Optimization (Opt) and sensitivity (Sens) tools are provided to analyze model parameters and evaluate the extent to which uncertainty in model parameters affects uncertainty in simulation results. A modified version of the National Weather Service's Extended Streamflow Prediction Program (ESP)2 provides forecasting capabilities using historic or synthesized meteorological data.

THE GIS WEASEL
The GIS Weasel has been designed to aid the preparation and analysis of spatial information for the estimation of parameters for lumped- and distributed-parameter physical-process models. It is composed of a combination of Arc/Info* GIS software, C language programs, and Unix shell scripts. All user interfaces are menu and map driven. This software focuses on the user's understanding of how to accommodate model assumptions regarding the treatment of variations in spatial attributes in the context of lumped- or distributed-parameter applications. The user is not required to have any knowledge of the command line operating procedures of Arc/Info.

The GIS Weasel provides tools to delineate, characterize, modify and parameterize an area and its associated "model response units" (MRUs) using a digital elevation model (DEM) and ancillary digital data (Fig. 2). MRUs are defined as areas delineated within a watershed, or region of interest, that reflect a model's treatment of spatially distributed characteristics, such as slope, aspect, elevation, soils, and vegetation. An MRU may represent an area whose character or composition is assumed to be homogenous with respect to one or more attributes or it may be heterogeneous with respect to all attributes.

In addition to an MRU map and estimates of model parameters for each MRU, the GIS Weasel also delineates the drainage network and computes the connectivity of MRUs with the drainage network. The location of data collection sites can be overlaid on the MRU map and the connectivity relations among MRUs, the drainage network, and the data sites also can be computed. The history of previously completed processing steps is also recorded to enable the user to reenter the computational process at any point and use an alternative computational procedure.

----------------------------------------
* The use of trade, product, industry, or firm names if for descriptive purpose only and does not imply endorsement by the U.S. Government.
----------------------------------------

 

Delineation Phase

In this phase, the DEM, existing in an Arc/Info GRID format, is used to derive additional topographic surfaces, including the area of application, drainage network, and preliminary MRU map. The GIS Weasel will fill depressions in the DEM in order to insure that all paths of hydrologic flow will eventually reach the edge of the DEM. Flow direction and accumulation surfaces are created from this input. The area of interest within the DEM is then delineated from either a set of automatically extracted watersheds, a user-selected point on the flow accumulation surface, or a pre-existing Arc/INFO coverage or GRID.

figure2

Next, the drainage network is extracted using a user-provided threshold value that represents the minimum upslope area to be associated with the initiation of a first-order channel3. The resulting drainage network is generated using the flow accumulation surface. The magnitude of the threshold effects the density of the drainage network by constraining the minimum area of a subwatershed. The user-defined first order drainage configuration may or may not reflect the actual configuration depending on the level of spatial aggregation desired. The interface allows the user to experiment with different thresholds until an acceptable value is found.

The last step in the delineation phase of the GIS Weasel is the selection of a methodology for sub-dividing the application area into MRUs. Currently, there are seven different methods available: 1-3) three types of irregular (traditional flow-plane based) polygons; 4) a fixed-interval lattice of square MRUs; 5) non-contiguous or "pixelated" MRUs (highlighting statistical similarities rather than spatial contiguity); 6) contributing areas associated with a coverage of points (e.g.: stream gages); and 7) use of a pre-existing Arc/Info coverage or GRID. A preliminary MRU map is manufactured based on the method selected.

The delineation phase is a highly automated sequence of processing steps with minimal user input. The characterization, modification, and parameterization phases, however, are more heuristic in nature and the options available are dependent on the type of delineation method selected. These three phases are not a fixed sequence of processing steps, but rather an open editing session where a user can iteratively evaluate, modify and parameterize the MRU map.

 

Characterization Phase

In this phase, a set of tools allows the user to explore the nature of the MRU map and any ancillary thematic information that may be associated with the area of application. Data derived from the original elevation grid (e.g.; slope, aspect) or other user-supplied grids of attribute data (e.g.; vegetation, soils) can be examined on the basis of individual or groups of MRUs, the entire watershed, or individual grid cells using various queries and displays of attribute distribution. User-defined boolean-based logical queries can be to used to display and evaluate various combinations of data attributes. Histograms of the distribution of selected attributes can displayed.

 

Modification Phase

This phase enables the user to make revisions to the MRU map based on information obtained during the characterization phase or on the basis of potential changes in the characteristics of the area of interest. An MRU map reflecting current conditions may be modified to reflect anticipated changes in land use as the result of anthropogenic or natural forces. MRUs can be grouped, divided, created or eliminated. The GIS Weasel provides version control and documentation to track modifications of MRU maps.

 

Parameterization Phase

This phase generates estimates of model parameters using MRU attributes and statistical measures of these attributes. Parameter estimates can be defined as measures of central tendency or measures of variance. Some parameters related to geographic and topographic characteristics, such as area, elevation, and slope, are common among many models and can be derived from the DEM used in the delineation phase. More module-specific parameters may require the use of multiple digital data coverages and the development of specific GIS procedures to compute parameter estimates. Developed procedures for selected modules in the MMS module library can be selected and applied using menu lists of modules and parameter options. Parameter estimates can be written directly into MMS parameter files for application in the selected model

INTEGRATED SYSTEMS APPLICATION
The integration of the GIS Weasel and MMS provide a flexible set of tools for application to a variety of environmental issues. One application of these systems is in the development of a water-resources management decision support system that address a number of environmental concerns, including water quantity, water quality, endangered species, and the optimization of operational objectives such as power generation and irrigation within existing environmental constraints4. The GIS Weasel provides the delineation, characterization, and parameterization facilities for a number of models in MMS. These include watershed models that provide streamflow to water storage reservoirs, channel flow and sediment transport models for routing water and sediment in selected channel reaches, and biological models that evaluate the effects of flow and sediment on stream biota.

A second application is the evaluation of the effects of alternative watershed and ecosystem management strategies5. Landscape-generation models provide spatial representations of alternative land-use patterns associated with alternative forest-management scenarios. Management scenarios might range from maximum timber yield to maximum habitat protection. Scenario outputs from the landscape-generation models can be input to the GIS Weasel and used to characterize and parameterize each scenario for application with a variety of watershed and ecosystem process models in MMS. The effects of these management strategies on issues such as streamflow, sediment, and selected chemical and biological processes can be evaluated.

While the GIS Weasel provides the mechanisms for delineating, characterizing, and estimating parameters for watersheds and MRUs, much work remains to be done on the development, testing, and identification of the specific parameterization procedures that are most appropriate for each of the environmental processes that may be simulated using MMS. The ability to recreate sets of spatial parameters for an application area is crucial to isolating and comparing the physical response predictions of different models. In addition, the development of consistent methodologies for extracting spatial information from different application areas enables the comparison of physical responses of the different areas, as well as the extension of calibrated model settings from application areas with observed data to application areas lacking relevant observations. As these procedures are defined for selected modules, they will be linked to the module. When a model is built or selected for execution, the parameterization procedures linked with each module will be incorporated in the GIS Weasel for automated execution if the needed spatial data bases are available.

Continued advances in physical and biological sciences, GIS technology, computer technology, and data resources will expand the need for a dynamic set of tools to apply these advances to a wide range of interdisciplinary research and operational environmental issues. The GIS Weasel and MMS provide a flexible framework in which to integrate and apply these tools.

 

REFERENCES


[1] G.H.Leavesley, P.J. Restrepo, L.G. Stannard, L.A.Frankowski, A.M. Sautins, MMS: A modeling framework for multidisciplinary research and operational applications, in GIS and Environmental Modeling: Progress and Research Issues, GIS World Books, Fort Collins, CO, p. 155 (1996).
[2] G.N. Day, Extended streamflow forecasting using NWSRFS. J. of Water Res. Planning and Manag., 111(2), p. 157 (1985).
[3] S.K. Jenson, J.O. Domingue, Extracting topographic structure from digital elevation data for geographic information system analysis, Photogrammetric Eng. and Remote Sens. 54, p. 1593 (1988).
[4] G.H. Leavesley, S.L. Markstrom, M.S. Brewer, R.J. Viger, The Modular Modeling System (MMS) -- The physical process modeling component of a database-centered decision support system for water and power management, Water, Air, and Soil Pollution 90, p. 303 (1996).
[5] G.H. Leavesley, G.E. Grant, S.L. Markstrom, R.J. Viger, M.S. Brewer, A modular modeling approach to watershed analysis and ecosystem management, in Proceedings of Watershed '96 A National Conference on Watershed Management, U.S. Government Printing Office, p. 794 (1996).

 

 

USA.gov logo   Take Pride in America 
button