THE ILWAS MODEL: FORMULATION AND APPLICATION STEVEN A. GHERINI, LINGFUNG MOK, ROBERT J. M. HUDSON, GEORGE F. DAVIS Tetra Tech., Inc., 3746Mt. Diablo Boulevard, Suite 300, Lafayette, CA 94549, U.S.A. CARL W. CHEN Systech., Inc., 3744 Mr. Diablo Boulevard, Suite 101, Lafayette, CA 94549, U.S.A. and ROBERT A. GOLDSTEIN Electric Power Research Institute, 3412 Hillview Avenue P.O. Box 10412, Palo Alto, CA 94303, U.S.A. (Received November 1, 1984; revised May 14, 1985) Abstract. The Integrated Lake-Watershed Acidification Study (ILWAS) model was developed to predict changes in surface water acidity given changes in the acidity of precipitation and dry deposition. The model routes precipitation through the forest canopy, soil horizons, streams and lakes using mass balance concepts and equations which relate flow to hydraulic gradients. The physical-chemical processes which change the acid-base characteristics of the water are simulated by rate (kinetic) and equilibrium expressions and include mass transfers between gas, liquid and solid phases. The aqueous constituents simulated include: pH, alkalinity, the major cations (Ca 2+ , Mg 2+ , K + , Na +, and NH4 + ) and anions (SO4 z- , NO3, CI-, F-), monomeric A1 and its inorganic and organic complexes, organic acid analogues and dissolved inorganic carbon (CT). Since free hydrogen ion (H + ) (and hence pH) is not conserved, its concentration is derived from the solution alkalinity and the total concentrations of inorganic C, organic acid, and monomeric AI. The ILWAS model has been used to predict changes in the acidity of Woods Lake (typical lake pH 4.5 to 5.0) and Panther Lake (typical lake pH 6 to 7) given reductions in total atmospheric S loads. The two basins are located within 30 km of each other in the Adirondack Mountains and receive similar acidic deposition. The response to a halving in the total atmospheric S load was basin-specific: In Panther Lake, little pH change occurred even 12 yr after the load reduction; in Woods Lake, the change was considerably larger. Hypothesis testing with the model has shown that the routing of water through soils (shallow versus deep flow) largely determines the extent to which incident precipitation is neutralized. Analysis of the two lake basins using the model and field data showed the watersheds to be net suppliers of base to the through-flowing water, although internal watershed production of strong acidity did occur. This internal production of acidity was approximately two-thirds the amount of the atmospheric load. 1. Introduction The ILWAS model was developed to predict changes in surface water acidity given changes in deposition acidity incident to forested ecosystems. Specifically, the model was developed to predict surface water H + and A1 concentrations because of their importance to fish. As conceptualization work progressed, it became obvious that the model must simulate both the routing of water to the lakes and the major alkalinity producing and consuming reactions occurring along the flow pathways. The result was a unified theory of lake acidification quantified by equations and algorithms embedded in the model code. Water, Air, and Soil Pollution 26 (1985) 425-459. 0049-6979/85.15. © 1985 by D. Reidel Publishing Company.