DECISION SUPPORT SYSTEM FOR MANAGING GROUND WATER RESOURCES IN THE CHOUSHUI RIVER ALLUVIAL IN TAIWAN1

ABSTRACT: Ground water is a vital water resource in the Choushui River alluvial fan in Taiwan. A significantly increased demand for water, resulting from rapid economic development, has led to large scale ground water extraction. Overdraft of ground water has considerably lowered the ground water level, and caused seawater intrusion, land subsidence, and other environmental damage. Sound ground water management thus is essential. This study presents a decision support system (DSS) for managing ground water resources in the Choushui River alluvial fan. This DSS integrates geographic information, ground water simulation, and expert systems. The geographic information system effectively analyzes and displays the spatially varied data and interfaces with the ground water simulation system to compute the dynamic behavior of ground water flow and solute transport in the aquifer. Meanwhile, a ground water model, MODFLOW‐96, is used to determine the permissible yield in the Choushui River alluvial fan. Additionally, an expert system of DSS employs the determined aquifer permissible yield to assist local government agencies in issuing water rights permits and managing ground water resources in the Choushui River alluvial fan.     

Nitrogen Subsidies from Hillslope Alder Stands to Streamside Wetlands and Headwater Streams,Kenai Peninsula,Alaska

We examined nitrogen transport and wetland primary production along hydrologic flow paths that link nitrogen‐fixing alder (Alnus spp.) stands to downslope wetlands and streams in the Kenai Lowlands, Alaska. We expected that nitrate concentrations in surface water and groundwater would be higher on flow paths below alder. We further expected that nitrate concentrations would be higher in surface water and groundwater at the base of short flow paths with alder and that streamside wetlands at the base of alder‐near flow paths would be less nitrogen limited than wetlands at the base of long flow paths with alder. Our results showed that groundwater nitrate‐N concentrations were significantly higher at alder‐near sites than at no‐alder sites, but did not differ significantly between alder‐far sites and no‐alder sites or between alder‐far sites and alder‐near sites. A survey of ¹⁵N stable isotope signatures in soils and foliage in alder‐near and no‐alder flow paths indicated the alder‐derived nitrogen evident in soils below alder is quickly integrated downslope. Additionally, there was a significant difference in the relative increase in plant biomass after nitrogen fertilization, with the greatest increase occurring in the no‐alder sites. This study demonstrates that streamside wetlands and streams are connected to the surrounding landscapes through hydrologic flow paths, and flow paths with alder stands are potential “hot spots” for nitrogen subsidies at the hillslope scale.     

ISOTOPE HYDROLOGY DYNAMICS OF RIVERINE WETLANDS IN THE KANKAKEE WATERSHED,INDIANA1

ABSTRACT: Wetland restoration activities may disturb shallow ground‐water flow dynamics. There may be unintentional sources of water flowing into a constructed wetland that could compromise the long‐term viability of a wetland function. Measurement of naturally‐occurring isotopes in the hydrosphere can provide an indication of provenance, flow paths or components, and residence times or ages of wetland ground‐water flow systems. Hydraulic head measurements may not provide sufficient detail of shallow flow disturbances and can be complemented by analyzing isotopes in waters flowing through the wetland. Two north‐central Indiana wetlands in the Kankakee watershed are being studied to determine the adequacy of wetland restoration activities. The native LaSalle wetland and the restored Hog Marsh wetland have contrasting ground‐water flow regimes. The conservative water isotopes ¹⁸O, ²H, and ³H, and selected solute isotopes ¹³C, ¹⁴C, ¹⁵N, ³⁴S, ⁸⁷Sr, and 206–208Pb, demonstrate the complexity of ground‐water flow in Hog Marsh compared to the established flow regime at the LaSalle wetland.     

Elemental contaminants in groundwater: A study of trace metals from residential area in the vicinity of an industrial area in Lagos, Nigeria

Groundwater (well water) from a residential area within the vicinity of an industrial estate in Lagos, Nigeria were sampled and analysed by Flame Atomic Absorption Spectroscopy for their heavy metals content. This was with a view of assessing the quality of the water, which was being used for domestic activities, especially, drinking usually without treatment. Total trace metal determination by mineral acid digestion of water samples was applied. This method proved to be better than an extractive concentration technique in the quality assurance protocols with the recovery range being 90.7 ± 0.006–97.6 ± 0.003%. Mean concentration of trace metals in water samples ranged from Fe: 0.05–0.47 mg l⁻¹; Al: 0.1–1.54 mg l⁻¹; Cu: 0.14–1.39 mg l⁻¹; Zn: 0.04–0.43 mg l⁻¹; Cd: trace–0.02 mg l⁻¹; Pb: trace–0.03 mg l⁻¹, Mn: 0.01–0.18 mg l⁻¹ and Ni: 0.02–0.11 mg l⁻¹. Physical parameters of water samples examined were within the drinking water safety limits except for conductivity. Results generally indicate the presence of heavy metal constituents in groundwater samples. Detection of metals such as cadmium and lead which have serious health implications above WHO and USEPA limits in drinking water gives cause for concern.     

Surface water pollution in three urban territories of Nepal, India, and Bangladesh

In South Asian countries such as Nepal, India, and Bangladesh, pollution of rivers is more severe and critical near urban stretches due to huge amounts of pollution load discharged by urban activities. The Bagmati River in the Kathmandu valley, the Yamuna River at Delhi, and peripheral rivers (mainly Buriganga River) of Dhaka suffer from severe pollution these days. The observed dry season average of biochemical oxygen demand (BOD) in all these rivers is in the range of 20–30 mg/liter and total coliform are as high as 10⁴–10⁵ MPN/100 ml. Per capita pollution load discharge of urban areas has been estimated to be about 31, 19, and 25 gBOD/capita/day in Bagmati, Yamuna, and the rivers of Dhaka, respectively. Regression analysis reveals pollution loads steadily increasing nearly in step with the trend in urbanization. The dissolved oxygen (DO) level of the Bagmati and Buriganga rivers is declining at an average annual rate of nearly 0.3 mg/liter/year. Unplanned urbanization and industrialization occurring in these cities may be largely responsible for this grave situation. Inadequate sewerage, on-site sanitation, and wastewater treatment facilities in one hand, and lack of effective pollution control measures and their strict enforcement on the other are the major causes of rampant discharge of pollutants in the aquatic systems.     

Estimating the Unknown Components of Nutrient Mass Balances for Forestry Plantations in Mine Rehabilitation, Upper Hunter Valley, New South Wales, Australia

Commercial forestry plantations as a postmining land use in the Upper Hunter Valley of New South Wales, Australia are restricted by both the poor nutrient availability of mining substrates and low regional rainfall. An experiment was conducted to investigate whether municipal waste products and saline groundwater from coal mining operations could improve early tree growth without impacting on the environment through salt accumulation and/or nutrient enrichment and changes in groundwater quality. Potential impacts were investigated by quantifying the nutrient cycling dynamics within the plantation using an input–output mass balance approach for exchangeable calcium (Ca²⁺), exchangeable magnesium (Mg²⁺), exchangeable potassium (K⁺), exchangeable sodium (Na⁺), nitrogen (N), and phosphorus (P). Measured inputs to and outputs from the available nutrient pool in the 0–30 cm of the overburden subsystem were used to estimate the net effect of unmeasured inputs and outputs (termed “residuals”). Residual values in the mass balance of the irrigated treatments demonstrated large leaching losses of exchangeable Ca, Mg, K, and Na. Between 96% and 103% of Na applied in saline mine-water irrigation was leached below the 0–30-cm soil profile zone. The fate of these salts beyond 30 cm is unknown, but results suggest that irrigation with saline mine water had minimal impact on the substrate to 30 cm over the first 2 years since plantation establishment. Accumulations of N and P were detected for the substrate amendments, suggesting that organic amendments (particularly compost) retained the applied nutrients with very little associated losses, particularly through leaching.     

Basic hydrologic studies for assessing impacts of flow diversions on riparian vegetation: Examples from streams of the eastern Sierra Nevada,California, USA

As the number of proposals to divert streamflow for power production has increased in recent years, interest has grown in predicting the impacts of flow reductions on riparian vegetation. Because the extent and density of riparian vegetation depend largely on local geomorphic and hydrologic setting, site-specific geomorphic and hydrologic information is needed. This article describes methods for collecting relevant hydrologic data, and reports the results of such studies on seven stream reaches proposed for hydroelectric development in the eastern Sierra Nevada, California, USA. The methods described are: (a) preparing geomorphic maps from aerial photographs, (b) using well level records to evaluate the influence of streamflow on the riparian water table, (c) taking synoptic flow measurements to identify gaining and losing reaches, and (d) analyzing flow records from an upstream-downstream pair of gages to document seasonal variations in downstream flow losses. In the eastern Sierra Nevada, the geomorphic influences on hydrology and riparian vegetation were pronounced. For example, in a large, U-shaped glacial valley, the width of the riparian strip was highly variable along the study reach and was related to geomorphic controls, whereas the study reaches on alluvial fan deposits had relatively uniform geomorphology and riparian strip width. Flow losses of 20% were typical over reaches on alluvial fans. In a mountain valley, however, one stream gained up to 275% from geomorphically controlled groundwater contributions.     

Parameter Uncertainty Reduction for SWAT Using Grace,Streamflow, and Groundwater Table Data for Lower Missouri River Basin1

Abstract: This study incorporates the newly available Gravity Recovery and Climate Experiment (GRACE) water storage data and water table data from well logs to reduce parameter uncertainty in Soil and Water Assessment Tool (SWAT) calibration using a SUFI2 (sequential uncertainty fitting) framework for the Lower Missouri River Basin. Model evaluations are performed in multiple stages using a multiobjective function consisting of multisite streamflow and GRACE water storage data as well as a groundwater component. Results show that (1) a model calibrated with both streamflow and GRACE data simultaneously can maintain the water balance for the whole basin, but may improperly partition surface flow and base flow. Additional inclusion of the groundwater constraint can significantly improve the model performance in groundwater hydrological processes. In our case, the estimation of specific yield of shallow aquifers has been increased to 10^?2 from previous much underestimated level (<10^?3). (2) The daily streamflow data are needed to confine the parameters related to water flow in channels such as the Manning’s coefficient, which are less sensitive to the monthly simulations. (3) Parameters are nonuniformly sensitive for different goal variables, and thus, proper specification of a prior distribution of parameters may be the key factor for global optimization algorithms to obtain stable and realistic model performance.     

Assessment of Water Supply as an Ecosystem Service in a Rural-Urban Watershed in Southwestern Mexico City

Studies from the ecosystem services perspective can provide a useful framework because they allow us to fully examine the benefits that humans obtain from socio-ecological systems. Mexico City, the second largest city in the world, has faced severe problems related to water shortages, which have worsened due to increasing population. Demand for space has forced changes in land cover, including covering areas that are essential for groundwater recharge. The city has 880 km² of forest areas that are crucial for the water supply. The Magdalena River Watershed was chosen as a model because it is a well-preserved zone within Mexico City and it provides water for the population. The general aim of this study was to assess the ecosystem service of the water supply in the Magdalena River Watershed by determining its water balance (SWAT model) and the number of beneficiaries of the ecosystem services. The results showed that the watershed provides 18.4 hm³ of water per year. Baseflow was dominant, with a contribution of 85%, while surface runoff only accounted for 15%. The zone provides drinking water to 78,476 inhabitants and could supply 153,203 potential beneficiaries. This work provides an example for understanding how ecosystem processes determine the provision of ecosystem services and benefits to the population in a rural–urban watershed in Mexico City.     

Modeling channel management impacts on river migration: a case study of Woodson bridge state recreation area,Sacramento river,California, USA

Understanding how hydraulic factors control alluvial river meander migration can help resource managers evaluate the long-term effects of floodplain management and bank stabilization measures. Using a numerical model based on the mechanics of flow and sediment transport in curved river channels, we predict 50 years of channel migration and suggest the planning and ecological implications of that migration for a 6.4-km reach (river miles 218–222) of the Sacramento River near the Woodson Bridge State Recreation Area, California, USA. Using four different channel management scenarios, our channel migration simulations suggest that: (1) channel stabilization alters the future channel planform locally and downstream from the stabilization; (2) rock revetment currently on the bank upstream from the Woodson Bridge recreation area causes more erosion of the channel bank at the recreation area than if the revetment were not present; (3) relocating the channel to the west and allowing subsequent unconstrained river migration relieves the erosion pressure in the Woodson Bridge area; (4) the subsequent migration reworks (erodes along one river bank and replaces new floodplain along the other) 26.5 ha of land; and (5) the river will rework between 8.5 and 48.5 ha of land in the study reach (over the course of 50 years), depending on the bank stabilization plan used. The reworking of floodplain lands is an important riparian ecosystem function that maintains habitat heterogeneity, an essential factor for the long-term survival of several threatened and endangered animal species in the Sacramento River area.     

High-Frequency Zone of River Desiccation Disasters in China and Influencing Factors

In recent years, the desiccation of the lower Yellow River has raised concerns in the government, public, and scientific community in China. Long-term and widespread desiccation of rivers is a disaster with many adverse environmental implications. It has been found in this study that there exists a high-frequency zone of river desiccation disasters at 34°–42°N in the North China Plain. The hazardous environment is characterized by semiarid climate, widely distributed thick loess in the basin and a “hanging river bed” in the plain as well as unfavorable man–water–land coupling relationships. In this setting, the sharply increased water diversion by man since the late 1950s led to the occurrence of river desiccation disasters in the lower reaches of the river in this area.     

Physiological and Morphological Response Patterns of <Emphasis Type="Italic">Populus deltoides</Emphasis> to Alluvial Groundwater Pumping

We examined the physiological and morphological response patterns of plains cottonwood [Populus deltoides subsp. monilifera (Aiton) Eck.] to acute water stress imposed by groundwater pumping. Between 3 and 27 July 1996, four large pumps were used to withdraw alluvial groundwater from a cottonwood forest along the South Platte River, near Denver, Colorado, USA. The study was designed as a stand-level, split-plot experiment with factorial treatments including two soil types (a gravel soil and a loam topsoil over gravel), two water table drawdown depths (∼0.5 m and >1.0 m), and one water table control (no drawdown) per soil type. Measurements of water table depth, soil water potential (Ψ_s), predawn and midday shoot water potential (Ψ_pd and Ψ_md), and D/H (deuterium/hydrogen) ratios of different water sources were made in each of six 600-m² plots prior to, during, and immediately following pumping. Two additional plots were established and measured to examine the extent to which surface irrigation could be used to mitigate the effects of deep drawdown on P. deltoides for each soil type. Recovery of tree water status following pumping was evaluated by measuring stomatal conductance (g _s ) and xylem water potential (Ψ_xp) on approximately hourly time steps from before dawn to mid-afternoon on 11 August 1996 in watered and unwatered, deep-drawdown plots on gravel soils. P. deltoides responded to abrupt alluvial water table decline with decreased shoot water potential followed by leaf mortality. Ψ_pd and percent leaf loss were significantly related to the magnitude of water table declines. The onset and course of these responses were influenced by short-term variability in surface and ground water levels, acting in concert with physiological and morphological adjustments. Decreases in Ψ_pd corresponded with increases in Ψ_md, suggesting shoot water status improved in response to stomatal closure and crown dieback. Crown dieback caused by xylem cavitation likely occurred when Ψ_pd reached −0.4 to −0.8 MPa. The application of surface irrigation allowed trees to maintain favorable water status with little or no apparent cavitation, even in deep-drawdown plots. Two weeks after the partial canopy dieback and cessation of pumping, g _s and Ψ_xp measurements indicated that water stress persisted in unwatered P. deltoides in deep-drawdown plots.     

Macroinvertebrate Responses to Constructed Riffles in the Cache River,Illinois, USA

Stream restoration practices are becoming increasingly common, but biological assessments of these improvements are still limited. Rock weirs, a type of constructed riffle, were implemented in the upper Cache River in southern Illinois, USA, in 2001 and 2003–2004 to control channel incision and protect high quality riparian wetlands as part of an extensive watershed-level restoration. Construction of the rock weirs provided an opportunity to examine biological responses to a common in-stream restoration technique. We compared macroinvertebrate assemblages on previously constructed rock weirs and newly constructed weirs to those on snags and scoured clay streambed, the two dominant substrates in the unrestored reaches of the river. We quantitatively sampled macroinvertebrates on these substrates on seven occasions during 2003 and 2004. Ephemeroptera, Plecoptera, and Trichoptera (EPT) biomass and aquatic insect biomass were significantly higher on rock weirs than the streambed for most sample periods. Snags supported intermediate EPT and aquatic insect biomass compared to rock weirs and the streambed. Nonmetric multidimensional scaling (NMDS) ordinations for 2003 and 2004 revealed distinct assemblage groups for rock weirs, snags, and the streambed. Analysis of similarity supported visual interpretation of NMDS plots. All pair-wise substrate comparisons differed significantly, except recently constructed weirs versus older weirs. Results indicate positive responses by macroinvertebrate assemblages to in-stream restoration in the Cache River. Moreover, these responses were not evident with more common measures of total density, biomass, and diversity.     

Comparison of Aquifer Sustainability Under Groundwater Administrations in Oklahoma and Texas1

Mittelstet, Aaron R., Michael D. Smolen, Garey A. Fox, and Damian C. Adams, 2011. Comparison of Aquifer Sustainability Under Groundwater Administrations in Oklahoma and Texas. Journal of the American Water Resources Association (JAWRA) 1‐8. DOI: 10.1111/j.1752‐1688.2011.00524.x Abstract: We compared two approaches to administration of groundwater law on a hydrologic model of the North Canadian River, an alluvial aquifer in northwestern Oklahoma. Oklahoma limits pumping rates to retain 50% aquifer saturated thickness after 20 years of groundwater use. The Texas Panhandle Groundwater Conservation District’s (GCD) rules limit pumping to a rate that consumes no more than 50% of saturated thickness in 50 years, with reevaluation and readjustment of permits every 5 years. Using a hydrologic model (MODFLOW), we simulated river‐groundwater interaction and aquifer dynamics under increasing levels of “development” (i.e., increasing groundwater withdrawals). Oklahoma’s approach initially would limit groundwater extraction more than the GCD approach, but the GCD approach would be more protective in the long run. Under Oklahoma rules more than half of aquifer storage would be depleted when development reaches 65%. Reevaluation of permits under the Texas Panhandle GCD approach would severely limit pumping as the 50% level is approached. Both Oklahoma and Texas Panhandle GCD approaches would deplete alluvial base flow at approximately 10% development. Results suggest periodic review of permits could protect aquifer storage and river base flow. Modeling total aquifer storage is more sensitive to recharge rate and aquifer hydraulic conductivity than to specific yield, while river leakage is most sensitive to aquifer hydraulic conductivity followed by specific yield.     

Assessment of the Hydrogeochemistry and Groundwater Quality of the Tarim River Basin in an Extreme Arid Region, NW China

The concentrations of the major and trace elements in the groundwater of the Tarim River Basin (TRB), the largest inland river basin of China, were analyzed before and during rainy seasons to determine the hydrogeochemistry and to assess the groundwater quality for irrigation and drinking purposes. The groundwater within the TRB was slightly alkaline and characterized by high ionic concentrations. The groundwater in the northern sub-basin was fresh water with a Ca²⁺–HCO₃ ^? water type, whereas the groundwater in the southern and central sub-basins was brackish with a Na⁺–Cl^? water type. Evaporite dissolution and carbonate weathering were the primary and secondary sources of solutes in the groundwater within the basin, whereas silicate weathering played a minor role. The sodium adsorption ratio (SAR), water quality index (WQI), and sodium percentage (%Na) indicated that the groundwater in the northern sub-basin was suitable for irrigation and drinking, but that in the southern and central sub-basins was not suitable. The groundwater quality was slightly better in the wet season than in the dry season. The groundwater could be used for drinking after treatment for B³⁺, F^?, and SO₄ ^2? and for irrigation after control of the sodium and salinity hazards. Considering the high corrosivity ratio of the groundwater in this area, noncorrosive pipes should be used for the groundwater supply. For sustainable development, integrated management of the surface water and the groundwater is needed in the future.     

Occurrence of arsenic in ground water in the Choushui River alluvial fan, Taiwan

An investigation of shallow ground water quality revealed that high arsenic (As) concentrations were found in both aquifers and aquitards in the southern Choushui River alluvial fan of Taiwan. A total of 655 geological core samples from 13 drilling wells were collected and analyzed. High As contents were found primarily in aquitards, to a maximum of 590 mg/kg. The contents were correlated with the locations of the marine sequences. Additionally, strong correlations among the As concentrations of core samples, the clay, and the geological age of the Holocene transgression were identified. Most of the As in ground water originated from the aquitard of the marine sequence. The high As content in marine formations with high clay contents may be attributable to the bioaccumulation of As in the sea organisms, which accrued and were deposited in the formation. A preliminary geogenic model of the origin of the high As concentration in the shallow sedimentary basin of the Choushui River alluvial fan of Taiwan is proposed.     

Signal processing and wavelet analysis of simultaneously registered blood pressure and laser Doppler flow signals during extremely low frequency electromagnetic field exposure in humans in vivo

Regulation of microcirculation and other physiological processes have strong non-linear character and involves complex of different processes, every process with own hierarchy in time and different frequencies. Traditional Fourier analysis does not provide sufficient power and resolution to elucidate characteristic of low vasomotor frequencies. Therefore, we apply a Time–frequency (wavelet) analysis on the signal obtained by Laser Doppler flow meter (LDF) at 25 healthy volunteers, exposed at the same time to low frequency electromagnetic fields, used for physiotherapy. Signal processing include Matlab based algorithms for digital signal processing (DSP) and Matlab Spectral analysis toolbox of simultaneous registered variations in Blood Pressure (BP), Laser Doppler Flow (LDF), and Intravital microscopy (IVM). It provides useful information about regulatory mechanisms and vegetative nervous system regulation of peripheral blood flow. Continuous changes in blood pressure variations and perfusion of extremities were measured prior and after 10, 20, and 30 min ELF-EMF (10, 16, 20, and 30 mT), exposure. After wavelet analysis of the blood flow signals and vasomotion changes signals, several frequency bands were distinguished: 0.0095–0.02 Hz; (α), 0.02–0.06 Hz; (β), 0.06–0.15 Hz; (γ), 0.15–0.4 Hz; (δ), and 0.4–1 Hz; (θ) for LDF data and 0.0095–0.4 Hz; (α), 0.4–0.75 Hz; (β), 0.75–0.9 Hz; (γ), 0.9–1.2 Hz; (δ), and 1.2–2 Hz; (θ) for IVM data. In this study, overlapping of some frequency bands between IVM and LDF data were found. Overlapping of the frequency bands has two ways of interpretation, one related with similarity of the structures and tissues and other related with output of ELF-EMF stimulation. We used also correlation and cross-correlation analysis to compare non-invasive (BP measurements and LDF) data, with invasive intravital microscopy (IVM) data (obtained on animals in vivo), during ELF-EMF stimulation. IVM data were used as a reference value, for certain information of possible mechanisms of biological response at the tissue and blood vessel level after ELF-EMF exposure with frequency in the range from 10 to 50 Hz and magnetic flux density of 20 mT. Comparative analysis of IVM and LDF, frequency bands show that they have statistical significant changes after ELF-EMF stimulation. Five subintervals were confirmed (α-, β-, γ-, δ-, and θ). The findings indicate that local ELF-EMF exposure at the constant temperature of the media increases skin blood flow at the upper extremities which have a contribution to the α-frequency band at IVM.     

The Spatio-Temporal Dynamic Pattern of Rural Residential Land in China in the 1990s Using Landsat TM Images and GIS

Through interpreting Landsat TM images, this study analyzes the spatial distribution of rural settlements in China in 2000. It calculates rural residential land percentage for every 1-km² cell. The entire country is divided into 33 regions to investigate the spatio-temporal dynamic patterns of rural residential land during the 1990s. According to the remote sensing survey, the rural residential land increased by 7.88 × 10⁵ ha in the 1990s. The increment of rural residential land was 0.55 million ha in 1990–1995 and 0.23 million ha in 1995–2000. In 1990–1995, rural residential land increased dramatically in the eastern regions such as the Yangtze River Delta, Pearl River Delta, and North China Plain, accounting for 80.80% of the national growth; the expansion in the western regions was much more moderate. In 1995–2000, the expansion of rural residential land in eastern regions slowed, accounting for only 58.54% of the increase at the national level, whereas the expansion in the western regions accelerated. Rapid rural residential development resulted from increasing home construction and the limited control on rural land. The great regional disparity reflected the regional economic development and land-use policy change. Our finding shows that nearly 60% of the rural residential area came from cropland.     

Improvement in HSPF's Low‐Flow Predictions by Implementation of a Power Law Groundwater Storage‐Discharge Relationship

We have enhanced the ability of a widely used watershed model, Hydrologic Simulation Program — FORTRAN (HSPF), to predict low flows by reconfiguring the algorithm that simulates groundwater discharge. During dry weather periods, flow in most streams consists primarily of base flow, that is, groundwater discharged from underlying aquifers. In this study, HSPF's groundwater storage‐discharge relationship is changed from a linear to a more general nonlinear relationship which takes the form of a power law. The nonlinear algorithm is capable of simulating streamflow recession curves that have been found in some studies to better match observed dry weather hydrographs. The altered version of HSPF is implemented in the Chesapeake Bay Program's Phase 5 Model, an HSPF‐based model that simulates nutrient and sediment loads to the Chesapeake Bay, and is tested in the upper Potomac River basin, a 29,950 km² drainage area that is part of the Bay watershed. The nonlinear relationship improved median Nash‐Sutcliffe efficiencies for log daily flows at the model's 45 calibration points. Mean absolute percent error on low‐flow days dropped in five major Potomac River tributaries by up to 12 percentage points, and in the Potomac River itself by four percentage points, where low‐flow days were defined as days when observed flows were in the lowest 5th percentile range. Percent bias on low‐flow days improved by eight percentage points in the Potomac River, from ?11 to ?3%.