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.     

REAL-TIME SNOW SIMULATION MODEL FOR THE MONONGAHELA RIVER BASIN1

ABSTRACT: The Pittsburgh District, U.S. Army Corps of Engineers, is responsible for operating two multipurpose reservoirs in the 7384 square mile (19198 square kilometer) Monongahela Basin. A third reservoir, presently under construction, will soon be operating. The real-time forecasting of runoff for operational purposes requires simulation of snow accumulation and snowmelt throughout the Basin during the winter season. This article describes capabilities of SNOSIM, a model being developed for performing such simulation. The application of this model as part of a comprehensive system of water control software, and some initial simulation results are presented.     

BASIN WATER MONITORING: HYDROLOGIC IMPACT OR LEGAL COMPLIANCE1

ABSTRACT: The Surface Mining Control and Reclamation Act of 1977 (PL 95–87) requires data collection prior to mining which will allow a determination to be made of the probable hydrologic consequences. Modeling strategies exist which allow for such determination, but selection of a specific model form will dictate the length of data record needed. The absence of a sound management strategy by the regulatory authority in‘Maryland does not result in valid data collection and suggests that the agency is concerned more with legal compliance than with integrated resource management.     

SHORT-TERM FORECASTING OF SNOWMELT RUNOFF USING ARMAX MODELS1

ABSTRACT: Time series models of the ARMAX class were investigated for use in forecasting daily riverflow resulting from combined snowmelt/rainfall. The Snowmelt Runoff Model (Martinec-Rango Model) is shown to have a form similar to the ARMAX model. The advantage of the ARMAX approach is that analytical model identification and parameter estimation techniques are available. In addition, previous forecast errors can be included to improve forecasts and confidence limits can be estimated for the forecasts. Diagnostic checks are available to determine if the model is performing properly. Finally, Kalman filtering can be used to allow the model parameters to vary continuously to reflect changing basin runoff conditions. The above advantages result in improved flow forecasts with fewer model parameters.     

A hierarchical approach to fisheries planning and modeling in the Columbia River Basin

The Columbia River Basin is the scene of a massive effort to restore populations of Pacific salmon (Oncorhynchus spp.) and steelhead (O. mykiss). Efficient restoration is confounded by a high level of complexity, competing sociopolitical goals and values, and uncertainty about key system properties. Simulation models and other tools of systems analysis are important to development of a comprehensive, regionally acceptable strategy. Hierarchy theory provides a useful paradigm for organized complexity within the Columbia Basin and the basis for a trilevel hierarchical structure for organizing and integrating models. Life-stage models compose the most basic simulation units at the lowest level in the proposed hierarchical modeling structure. Each life-stage model simulates a distinct period in the life cycle of anadromous salmonids. Population models at the intermediate level simulate the complete life cycles of salmon and steelhead populations. At the highest level in the hierarchy, interpopulation models simulate extensive, long-term processes that affect multiple species and stocks. A hierarchical system of models is preferable to a single model or to a group of models lacking formal structure. A principal advantage is that models have the correct spatial and temporal resolution for analyzing questions at different scales. A hierarchical structure also facilitates the flow of information among models, and aids in understanding the impacts of uncertainty. Constructing a hierarchy of models should involve both bottom-up and top-down perspectives that maintain logical consistency among models, while allowing unique model structures appropriate for each level in the hierarchy.     

DESIGN BASIS FOR EVERGLADES STORMWATER TREATMENT AREAS1

ABSTRACT: The State of Florida (1994) has adopted a plan for addressing Everglades eutrophication problems by reducing anthropogenic phosphorus loads. The plan involves implementation of Best Management Practices in agricultural watersheds and construction of regional treatment marshes (Stormwater Treatment Areas or STA's). This paper describes the development, testing, and application of a mass-balance model for sizing STA's to achieve treatment objectives. The model is calibrated and tested against peat and water-column data collected in Water Conservation Area-2A (WCA-2A), where phosphorus dynamics and eutrophication impacts have been intensively studied. The 26-year-average rate of phosphorus accretion in peat is shown to be proportional to average water-column phosphorus concentration, with a proportionality constant of 10.2 m/yr (90 percent Confidence Interval = 8.9 to 11.6 m/yr). Spatial and temporal variations in marsh water-column data suggest that drought-induced recycling of phosphorus was important during periods of low stage in WCA-2A. Maintaining wet conditions will be important to promote phosphorus removal in STA's. Sensitivity analysis of STA performance is conducted over the range of uncertainty in model parameter estimates to assess the adequacy of the model as a basis for STA design.     

NUMERICAL MODELING OF SALTWATER INTRUSION IN THE NORTHERN GUAM LENS1

ABSTRACT This paper deals with the application of a two-dimensional, saltwater intrusion model to the aquifer in Northern Guam. The model used finite element theory and the Galerkin, weighted-residual technique as its basis. The Northern Guam lens was discretized into 299 linear, triangular elements and 189 nodes. The model was calibrated using 1978 hydrologic data. The output of the model was compared with measured water levels in six observation wells. The calibrated values of permeability and porosity were then used to verify the model using 1979 data. A calibrated and verified model can be used to make an infinite variety of management and planning studies. In this study, three applications are provided that would be considered typical management runs. Steady state runs were made to compare the four conditions of no pumping, 1978 pumping levels, twice 1978 pumping levels, and five times 1978 pumping levels. The water levels due to these conditions are shown in plan and in cross sections of the aquifer. The effect of zero recharge to the aquifer is next demonstrated for the pumping levels existing during 1978. The final run shows how long the aquifer takes to reach steady state when the pumping rate is increased from the 1978 pumping level to twice that value. The program can be used for numerous other studies for management and planning purposes.     

VELOCITY EQUATION FOR WATER QUALITY MODELING IN GEORGIA1

ABSTRACT: Tabletop water quality modeling still plays an important role in the water pollution control activities of the Georgia Environmental Protection Division. Tabletop models are those developed with out the aid of extensive field data. One important component of GEORGIA DOSAG, our basic water quality model, is the equation used to predict flow through velocity. However, Georgia is characterized by wide physiographic diversity which reduces the effectiveness of uncalibrated velocity equations. Using 15 years of accumulated time-of-travel studies, a series of empirical velocity equations were developed and calibrated to various physiographic conditions in Georgia. Equations are available for each major soil province and for three stream flow ranges within each province - Q<100 cfs, 100<Q<1000 cfs, and Q>1000 cfs. Now, in the absence of extensive field data, we have data based velocity equations which can be tailored to each site under study.     

A dynamic model for intertemporal allocation of old-growth forests in the Pacific Northwest

Across the globe, continued policy debates regarding the management of old-growth forests center around the difficult task of balancing economic and ecological considerations. Though the forests of the Pacific Northwest United States are among the most studied old-growth ecosystems, ecological and economic analyses have yielded public land management directives that remain controversial. Specifically, the recently adopted Northwest Forest Plan lacks explicit goals for maintaining intergenerational equity for the use of forest resources and the diversity of old-growth ecosystems. Unlike previous studies which rely on monetary quantification of costs and benefits, this study develops and applies a conceptual framework for evaluating socially optimal Pacific Northwest old-growth forest utilization strategies. Conditions for the optimal management of old-growth forests are derived using dynamic programming. The objective function synthesizes relevant biological and economic attributes of the old-growth allocation problem. Results in the form of extraction paths are compared given social pressure for consumptive and non-consumptive benefits, as well as different planning horizons, rates of social time preference, and environmental variance. Lengthening the planning horizon results in a vast divergence of optimal policies in the absence of discounting. Extraction rates appear to approach zero as the planning horizon approaches infinity. While higher rates of social time preference increase the rate of extraction, forest stocks remaining at the terminal time period equal levels remaining with a lower discount rate. Increasing environmental variance results in a higher level of stock remaining at the terminal time period. This analysis, while specific to the old-growth controversy of the Pacific Northwest, does provide general guidelines for addressing similar problems of multiple uses of natural areas, particularly where such uses are mutually incompatible, or where one use may be irreversibly destructive to another.     

ARCGIS‐SWAT: A GEODATA MODEL AND GIS INTERFACE FOR SWAT1

This paper presents ArcGIS‐SWAT, a geodata model and geographic information system (GIS) interface for the Soil and Water Assessment Tool (SWAT). The ArcGIS‐SWAT data model is a system of geodatabases that store SWAT geographic, numeric, and text input data and results in an organized fashion. Thus, it is proposed that a single and comprehensive geodatabase be used as the repository of a SWAT simulation. The ArcGIS‐SWAT interface uses programming objects that conform to the Component Object Model (COM) design standard, which facilitate the use of functionality of other Windows‐based applications within ArcGIS‐SWAT. In particular, the use of MS Excel and MATLAB functionality for data analysis and visualization of results is demonstrated. Likewise, it is proposed to conduct hydrologic model integration through the sharing of information with a not‐model‐specific hub data model where information common to different models can be stored and from which it can be retrieved. As an example, it is demonstrated how the Hydrologic Modeling System (HMS) ‐ a computer application for flood analysis ‐ can use information originally developed by ArcGIS‐SWAT for SWAT. The application of ArcGIS‐SWAT to the Seco Creek watershed in Texas is presented.