Discharge and Nutrient Transport between Lakes in a Hydrologically Complex Area of Voyageurs National Park,Minnesota, 2010‐2012

An acoustic Doppler velocity meter (ADVM) was deployed in the narrows between Namakan and Kabetogama Lakes in Voyageurs National Park, Minnesota, from November 3, 2010, through October 3, 2012. The ADVM can account for wind, seiche, and changing flow direction in hydrologically complex areas. The objectives were to (1) estimate discharge and document the direction of water flow, (2) assess whether specific conductance can be used to determine flow direction, and (3) document nutrient and chlorophyll a concentrations at the narrows. The discharge direction through the narrows was seasonal. Water generally flowed out of Kabetogama Lake and into Namakan Lake throughout the ice‐covered season. During spring, water flow was generally from Namakan Lake to Kabetogama Lake. During the summer and fall, the water flowed in both directions, affected in part by wind. Water flowed into Namakan Lake 70% of water year 2011 and 56% of water year 2012. Nutrient and chlorophyll a concentrations were highest during the summer months when water‐flow direction was unpredictable. The use of an ADVM was effective for assessing flow direction and provided flow direction under ice. The results indicated the eutrophic Kabetogama Lake may have a negative effect on the more pristine Namakan Lake. The results also provide data on the effects of the current water‐level management plan and may help determine if adjustments are necessary to help protect the aquatic ecosystem of Voyageurs National Park.     

Assessment of Future Climate Change Impact on Water Quality of Chungju Lake,South Korea,Using WASP Coupled with SWAT

This study is to evaluate the future potential impact of climate change on the water quality of Chungju Lake using the Water Quality Analysis Simulation Program (WASP). The lake has a storage capacity of 2.75 Gm³, maximum water surface of 65.7 km², and forest‐dominant watershed of 6,642 km². The impact on the lake from the watershed was evaluated by the Soil and Water Assessment Tool (SWAT). The WASP and SWAT were calibrated and validated using the monthly water temperatures from 1998 to 2003, lake water quality data (dissolved oxygen, total nitrogen [T‐N], total phosphorus [T‐P], and chlorophyll‐a [chl‐a]) and daily dam inflow, and monthly stream water quality (sediment, T‐N, and T‐P) data. For the future climate change scenario, the MIROC3.2 HiRes A1B was downscaled for 2020s, 2050s, and 2080s using the Change Factor statistical method. The 2080s temperature and precipitation showed an increase of +4.8°C and +34.4%, respectively, based on a 2000 baseline. For the 2080s watershed T‐N and T‐P loads of up to +87.3 and +19.6%, the 2080s lake T‐N and T‐P concentrations were projected to be 4.00 and 0.030 mg/l from 2.60 and 0.016 mg/l in 2000, respectively. The 2080s chl‐a concentration in the epilimnion and the maximum were 13.97 and 52.45 μg/l compared to 8.64 and 33.48 μg/l in 2000, respectively. The results show that the Chungju Lake will change from its mesotrophic state of 2000 to a eutrophic state by T‐P in the 2020s and by chl‐a in the 2080s. Editor's note: This paper is part of a featured series on Korean Hydrology. The series addresses the need for a new paradigm of river and watershed management for Korea due to climate and land use changes.     

SIMPLE TROPHIC STATE CLASSIFICATION OF THE CANADIAN NEARSHORE WATERS OF THE GREAT LAKES1

ABSTRACT: Trophic classification of the Canadian nearshore waters of the Great Lakes is attempted using summer, surface water quality data for the early 1970's. A generalized Composite Trophic Index is developed using paired linear relationships for total phosphorus, chlorophyll a, and Secchi depth data for 66 defined nearshore regions. The chlorophyll a and total phosphorus relationship indicates that the nearshore waters contain a low chlorophyll a concentration for a given total phosphorus concentration than observed for the open waters of the Great Lakes or for smaller Canadian lakes. The most eutrophic nearshore regions occur in areas of relatively restricted circulation and/or high nutrient loadings. These include the Bay of Quinte, Toronto and Hamilton harbours, and portions of Lake We's Western Basin. Lakes Huron and Superior are generally oligotrophic, except for some embayments. Although nearshore water quality is highly variable, this apprach represents a reasonable compromise with respect to analytical complexity. The Composite Trophic Index removes biases introduced through the use of a single trophic state indicator and uniquely describes the nearshore water quality in terms generally comparable to other water bodies.     

Integrating Seasonal Information on Nutrients and Benthic Algal Biomass into Stream Water Quality Monitoring

Benthic chlorophyll a (BChl a) and environmental factors that influence algal biomass were measured monthly from February through October in 22 streams from three agricultural regions of the United States. At‐site maximum BChl a ranged from 14 to 406 mg/m² and generally varied with dissolved inorganic nitrogen (DIN): 8 out of 9 sites with at‐site median DIN >0.5 mg/L had maximum BChl a >100 mg/m². BChl a accrued and persisted at levels within 50% of at‐site maximum for only one to three months. No dominant seasonal pattern for algal biomass accrual was observed in any region. A linear model with DIN, water surface gradient, and velocity accounted for most of the cross‐site variation in maximum chlorophyll a (adjusted R² = 0.7), but was no better than a single value of DIN = 0.5 mg/L for distinguishing between low and high‐biomass sites. Studies of nutrient enrichment require multiple samples to estimate algal biomass with sufficient precision given the magnitude of temporal variability of algal biomass. An effective strategy for regional stream assessment of nutrient enrichment could be based on a relation between maximum BChl a and DIN based on repeat sampling at sites selected to represent a gradient in nutrients and application of the relation to a larger number of sites with synoptic nutrient information.     

Mesoscale survey of western and northwestern Irish lakes – Spatial and aestival patterns in trophic status and phytoplankton community structure

Inland water bodies are considered as integrated parts of the landscape and the monitoring of water quality and aquatic resources need to be addressed on a regional basis for optimal assessment and management. In this study, a simple stratified sampling scheme was applied to a mesoscale survey of western and northwestern Irish lakes, which was carried out to identify, based on the distribution patterns of phytoplankton biomass, potential associations between lake trophic state and land cover attributes. Phytoplankton community analysis was also performed to determine whether taxa associations reflected meteorology-linked aestival succession or specific spatial distributions. The assessment was based on the typology of hydrogeomorphological and land cover attributes of river catchments through ArcGIS analysis. Sampling was carried out in 50 lakes and during a 15-week period in summer 2009. Results showed a general longitudinal gradient in the trophic status of the lakes sampled, with a greater frequency of mesotrophic lakes in the eastern part of the study area where land cover is dominated by agricultural surfaces. Significant relationships (p < 0.010) were found between chlorophyll-a concentration and the proportion of river catchment surface covered by agriculture land and wetlands, findings which might be considered further as proxies for developing an eutrophication risk index. Multivariate analysis of phytoplankton community data clustered the sampled lakes into three assemblages, with ordination along axis 1 being significantly correlated to time and temperature (p < 0.006). There was greater frequency of occurrence of diatoms in lakes from cluster III (Kruskal–Wallis, p < 0.05, H = 6.34, df = 2, n = 49), concomitant to lower chlorophyll-a concentrations, lake surface temperatures and Secchi depths, reflecting meteorological conditions dominated by precipitations. Those results support the potential of mesoscale surveys to assess water quality variables and detect environmental patterns at regional scales.     

Environmental Control and Limnological Impacts of a Large Recurrent Spring Bloom in Lake Washington,USA

A series of statistical analyses were used to identify temporal and spatial patterns in the phytoplankton and nutrient dynamics of Lake Washington, an mesotrophic lake in Washington State (USA). These analyses were based on fortnightly or monthly samples of water temperature, Secchi transparency, ammonium (NH₄), nitrate (NO₃), inorganic phosphorus (IP), total nitrogen (TN), total phosphorus (TP), dissolved oxygen (DO), pH and chlorophyll a (chl a) collected during 1995–2000 from 12 stations. Lake Washington has a very consistent and pronounced annual spring diatom bloom which occurs from March to May. During this bloom, epilimnetic chl a concentrations peak on average at 10 μg/L, which is 3 times higher than chl a concentrations typically seen during summer stratified conditions. The spring bloom on average comprised 62% diatoms, 21% chlorophytes and 8% cyanobacteria. During summer stratification, diatoms comprised 26% of the phytoplankton community, chlorophytes 37% and cyanobacteria 25%. Cryptophytes comprised approximately 8% of the community throughout the year. Overall, 6 phytoplankton genera (i.e., Aulacoseira, Fragilaria, Cryptomonas, Asterionella, Stephanodiscus, and Ankistrodesmus) cumulatively accounted for over 50% of the community. These analyses also suggest that the phytoplankton community strongly influences the seasonality of NO₃, IP, DO, pH and water clarity. According to a MANOVA, seasonal fluctuations explained 40% of the total variability for the major parameters, spatial heterogeneity explained 10% of variability, and the seasonal-spatial interaction explained 10% of variability. Distinctive patterns were identified between offshore and inshore sampling stations. The results of our analyses also suggest that spatial variability was substantial, but much smaller than temporal variability.     

Estimating Sediment and Nutrient Loads in Four Western Lake Superior Streams

Total suspended solids (TSS) and total phosphorus (TP) have been shown to be strongly correlated with turbidity in watersheds. High‐frequency in situ turbidity can provide estimates of these potential pollutants over a wide range of hydrologic conditions. Concentrations and loads were estimated in four western Lake Superior trout streams from 2005 to 2010 using regression models relating continuous turbidity data to grab sample measures of TSS and TP during differing flow regimes. TSS loads estimated using the turbidity surrogate were compared with those made using FLUX software, a standard assessment technique based on discharge and grab sampling for TSS. More traditional rating curve methodology was not suitable because of the high variability in the particulates vs. discharge relationship. Stream‐specific turbidity and TSS data were strongly correlated (r² = 0.5 to 0.8; p < 0.05) and less so for TP (r² = 0.3 to 0.7; p < 0.05). Near‐continuous turbidity monitoring (every 15 min) provided a good method for estimating both TSS and TP concentration, providing information when manual sample collection was unlikely, and allowing for detailed analyses of short‐term responses of flashy Lake Superior tributaries to highly variable weather and hydrologic conditions while the FLUX model typically resulted in load estimates greater than those determined using the turbidity surrogate, with 17/23 stream years having greater FLUX estimates for TSS and 18/23 for TP.     

A CARLSON-TYPE TROPHIC STATE INDEX FOR NITROGEN IN FLORIDA LAKES1

: A data base consisting of predominantly nitrogen limited Florida lakes from the National Eutrophication Survey (NES) was used to develop a trophic state index based on total nitrogen concentration. This index was compared with Carlson's (1977) index based on total phosphorus concentration, and the lesser of the two values for each lake was averaged with indices based on Secchi disk transparency and chlorophyll a concentration to assess the trophic state of the 40 Florida NES lakes.     

Nitrate Dynamics in Two Streams Impacted by Wastewater Treatment Plant Discharge: Point Sources or Sinks?

We examined nitrate processing in headwater stream reaches downstream of two wastewater treatment plant outfalls during low streamflow. Our objectives were to quantify nitrate mass flux before and after effluent discharge and to use field and laboratory techniques to assess the mechanism of nitrate uptake. Microcosm experiments were utilized to determine the location of nitrate processing, and molecular biomarkers were used to detect and quantify microbial denitrification. At one site, downstream nitrate mass flux was significantly (p = 0.01) lower than sum of upstream and wastewater effluent fluxes, indicating rapid stream assimilation of incoming nitrate in the vicinity of the point source. Microcosm experiments supported the theory that nitrate processing occurs in sediments. Molecular assays for denitrifcation‐associated functional genes nosZ, nirS, and nirK, provided evidence that effluent contained enriched denitrifying communities relative to ambient stream water. Nitrate loss at the site with greater uptake was correlated with sulfate loss (p < 0.01; r² = 0.86), suggesting a possible link between sulfate reducing bacteria and denitrifying bacterial communities. Results suggest there is an opportunity to better understand nitrate dynamics in cases where point sources may act as point sinks under specific sets of conditions.     

Resolution and Analysis of Spatial Variations and Patterns in an Urban Lake with Rapid Profiling Instrumentation

Rapid response vertical profiling instrumentation was used to document spatial variability and patterns in a small urban lake, Onondaga Lake, associated with multiple drivers. Paired profiles of temperature, specific conductance (SC), turbidity (T_n), fluorometric chlorophyll a (Chl_f), and nitrate nitrogen (NO₃^?) were collected at >30 fixed locations (a “gridding”) weekly, over the spring to fall interval of several years. These gridding data are analyzed (1) to characterize phytoplankton (Chl_f) patchiness in the lake's upper waters, (2) to establish the representativeness of a single long‐term site for monitoring lake‐wide conditions, and (3) to resolve spatial patterns of multiple tracers imparted by buoyancy effects of inflows. Multiple buoyancy signatures were resolved, including overflows from less dense inflows, and interflows to metalimnetic depths and underflows to the bottom from the plunging of more dense inputs. Three different metrics had utility as tracers in depicting the buoyancy signatures as follows: (1) SC, for salinity‐enriched tributaries and the more dilute river that receives the lake's outflow, (2) T_n, for the tributaries during runoff events, and (3) NO₃^?, for the effluent of a domestic waste treatment facility and from the addition of NO₃^? solution to control methyl mercury. The plunging inflow phenomenon, which frequently prevailed, has important management implications.     

Hydrogeomorphic Reference Condition and Its Relationship with Macroinvertebrate Assemblages in Southeastern U.S. Sand Hills Streams

Defining stream reference conditions is integral to providing benchmarks to ecological perturbation. We quantified channel geometry, hydrologic and environmental variables, and macroinvertebrates in 62 low‐gradient, SE United States (U.S.) Sand Hills (Level IV ecoregion) sand‐bed streams. To identify hydrogeomorphic reference condition (HGM), we clustered channel geometry deviation from expectations given watershed area (Aws), resulting in two HGM groups discriminated by area at the top of bank (Atob) residuals <0.6 m² and >0.6 m² predicted to be HGM reference/nonreference streams, respectively. Two independent partial least squares discriminate analyses used (1) hydrologic/environmental variables and (2) macroinvertebrate mean trait values (mT) on 10 reference/nonreference stream pairs of similar Aws for classification validation. Nonreference streams had flashier hydrographs and altered flow magnitudes, lower organic matter, coarser substrate, higher pH/specific conductivity compared with reference streams. Macroinvertebrate assemblages corresponded to HGM groupings, with mT indicative of multivoltinism, collector‐gatherer functional feeding groups, fast current‐preference taxa, and lower Ephemeroptera, Plecoptera, and Trichoptera richness and biotic integrity in nonreference streams. HGM classifications in Sand Hills, sand‐bed streams were determined from channel geometry. This easily implemented classification is indicative of contemporary hydrologic disturbance resulting in contrasting macroinvertebrate assemblages.     

RELATIONSHIPS BETWEEN WATER QUALITY AND PHOSPHORUS CONCENTRATIONS FOR PUGET SOUND REGION LAKES1

Regression relationships were developed between summer mean total phosphorus (P) concentrations in near-surface water and both chlorophyll a concentrations and Secchi disc transparency for Puget Sound region lakes. Total P concentrations in the lakes studied ranged from 7 to 66 μ/L. The relationship between total P and chlorophyll a, based on data from 69 lakes, explained 57 percent of the variance in chlorophyll a. Predicted chlorophyll a concentrations and 95 percent confidence intervals ranged from 1 ⁺³_-0.5μg/L for 7 μg/L P to about ⁺³⁵_-10μ/L for 66 μ/L P. The relationship between total P and Secchi disc, based on data from 71 lakes, explained 53 percent of the variance in Secchi disc. Predicted Secchi disc transparencies and 95 percent confidence intervals ranged from 5.5 ^+5.5_-3.0 m for 7 μ/L P to 1.4 ^+1.5_-0.7 m for 66 μ/L P.     

OPTICAL HETEROGENEITY IN GREEN BAY1

Differences in light penetration and light attenuating components and processes are documented along 112 km of the major (NE/SW) axis of Green Bay (Lake Michigan) during a three-day cruise (May 25–27, 1982). Measurements included diffuse attenuation of downwelling irradiance (k_d), Secchi disk transparency (SD), phytoplankton pigments (chlorophyll a and phaeophytin), turbidity (T_n), and dissolved color (absorbance). The relative importance of absorption and scattering to attenuation was calculated from paired measurements of k_d and SD. Absorption and scattering coefficients were calculated; value estimates were supported by a strong linear relationship between the scattering coefficient (b) and T_n (b = 0.99 *T_n; r²= 097, n = 12). Light attenuation characteristics, including the extent of light penetration and the magnitudes and relative importance of absorption, scattering and individual attenuating components, were found to be heterogenous in space. This heterogeneity is due to the characteristics and positions of entry of fluvial discharges to the bay as they influence levels of dissolved color (Gelbstoff), phytoplankton standing crop, and inorganic particulates. Identification of key processes regulating light penetration and their potential for response to pollution control measures can aid in the development of a water quality management plan for Green Bay.     

Climate and Water Conflicts Coevolution from Tropical Development and Hydro‐Climatic Perspectives: A Case Study of Costa Rica

Costa Rica is a nation with a vast wealth of water resources; however, recently the country has faced water conflicts (WC) due to social, economic, legal, and political impediments in response to limited water availability during El Niño events and inefficient use of its water resources. This study presents a spatial distribution and temporal analysis of WC in Costa Rica from 2005 to 2015. In total, 719 WC were analyzed of which 54% were among private individuals and government. The largest urban areas and the Grande de Tárcoles Basin were identified as the main “hot spot” for the conflicts. WC were mainly caused by spills of wastewater, water pollution, water shortage, infrastructure damage, and flooding, and can be predicted using a multiple linear model including the population size and the number of hydro‐meteorological events (HME) (R² = 0.77). The identified HME also coevolved significantly with the changes in precipitation regimes (r = 0.67, p = 0.021). Our results suggest that there is a need to recognize that water infrastructure longevity across the country concatenates and amplifies WC, mainly in the most populated area located in the Central Valley. Implications of our findings include the need for truly integrated water resources management plans that include, for example, WC as indicators of hydro‐climatic changing conditions and water supply and sanitation infrastructure status.     

Benthic Algal (Periphyton) Growth Rates in Response to Nitrogen and Phosphorus: Parameter Estimation for Water Quality Models

Nitrogen (N) and phosphorus (P) are significant pollutants that can stimulate nuisance blooms of algae. Water quality models (e.g., Water Quality Simulation Program, CE‐QUAL‐R1, CE‐QUAL‐ICM, QUAL2k) are valuable and widely used management tools for algal accrual due to excess nutrients in the presence of other limiting factors. These models utilize the Monod and Droop equations to associate algal growth rate with dissolved nutrient concentration and intracellular nutrient content. Having accurate parameter values is essential to model performance; however, published values for model parameterization are limited, particularly for benthic (periphyton) algae. We conducted a 10‐day mesocosm experiment and measured diatom‐dominated periphyton biomass accrual through time as chlorophyll a (chl a) and ash‐free dry mass (AFDM) in response to additions of N (range 5–11,995 µg nitrate as nitrogen [NO₃‐N]/L) and P (range 0.89–59.51 µg soluble reactive phosphorus/L). Resulting half‐saturation coefficients and growth rates are similar to other published values, but minimum nutrient quotas are higher than those previously reported. Saturation concentration for N ranged from 150 to 2,450 µg NO₃‐N/L based on chl a and from 8.5 to 60 µg NO₃‐N/L when based on AFDM. Similarly, the saturation concentration for P ranged from 12 to 29 µg‐P/L based on chl a, and from 2.5 to 6.1 µg‐P/L based on AFDM. These saturation concentrations provide an upper limit for streams where diatom growth can be expected to respond to nutrient levels and a benchmark for reducing nutrient concentrations to a point where benthic algal growth will be limited.     

Integrating Temporal Inequality into Conservation Planning to Improve Practice Design and Efficacy

In contrast to spatial inequality, there are currently no methods for leveraging information on temporal inequality to improve conservation efficacy. The objective of this study was to use Lorenz curves to quantify temporal inequality in surface runoff and tile drainage, identify controls on nutrient loading in these flowpaths, and develop design flows for structural conservation practices. Surface runoff (n = 94 site‐years) and tile drainage (n = 90 site‐years) were monitored on 40 fields in Ohio. Results showed, on average, 80% of nitrate‐nitrogen, soluble reactive phosphorus (P), and total P loads occurred between 7 and 12 days per year in surface runoff and between 32 and 58 days per year in tile drainage. Similar temporal inequality between discharge and load provided evidence that loading was transport‐limited and highlighted the critical role hydrologic connectivity plays in nutrient delivery from tile‐drained fields. Design flow criterion for sizing structural practices based on load reduction goals was developed by combining Lorenz curves and flow duration curves. Comparing temporal inequality between fields and the Maumee River, the largest tributary to the western Lake Erie Basin, revealed challenges associated with achieving watershed load reduction goals with field‐scale conservation. In‐field (i.e., improved nutrient and water management), edge‐of‐field (i.e., structural practices), and instream practices will all be required to meet nutrient reduction goals from tile‐drained watersheds.     

Decreased Runoff Response to Precipitation,Little Missouri River Basin,Northern Great Plains,USA

High variability in precipitation and streamflow in the semiarid northern Great Plains causes large uncertainty in water availability. This uncertainty is compounded by potential effects of future climate change. We examined historical variability in annual and growing season precipitation, temperature, and streamflow within the Little Missouri River Basin and identified differences in the runoff response to precipitation for the period 1976‐2012 compared to 1939‐1975 (n = 37 years in both cases). Computed mean values for the second half of the record showed little change (<5%) in annual or growing season precipitation, but average annual runoff at the basin outlet decreased by 22%, with 66% of the reduction in flow occurring during the growing season. Our results show a statistically significant (p < 0.10) 27% decrease in the annual runoff response to precipitation (runoff ratio). Surface‐water withdrawals for various uses appear to account for <12% of the reduction in average annual flow volume, and we found no published or reported evidence of substantial flow reduction caused by groundwater pumping in this basin. Results of our analysis suggest that increases in monthly average maximum and minimum temperatures, including >1°C increases in January through March, are the dominant driver of the observed decrease in runoff response to precipitation in the Little Missouri River Basin.     

PRISM vs. CFSR Precipitation Data Effects on Calibration and Validation of SWAT Models

Precipitation is one of the most important drivers in watershed models. Our objective was to compare two sources of interpolated precipitation data in terms of their effect on calibration and validation of two Soil and Water Assessment Tool (SWAT) models. One model was a suburban watershed in metropolitan Atlanta, Georgia. The precipitation sources were Parameter‐elevation Relationships on Independent Slopes Model (PRISM) data on a 4‐km grid and climate forecast system reanalysis (CFSR) data on a 38‐km grid. The PRISM data resulted in a better fit to the calibration data (Nash Sutcliffe efficiency [NSE] = 0.64, Kling‐Gupta efficiency [KGE] = 0.74, p‐factor = 0.84, and r‐factor = 0.43) than the CFSR data (NSE = 0.47, KGE = 0.53, p‐factor = 0.67, and r‐factor = 0.39). Validation results were similar. Sensitive parameters were similar in both the PRISM and CFSR models, but fitted values indicated more rapid groundwater flow to the streams with the PRISM data. The same comparison was made in the Big Creek watershed located approximately 1,000 km away, in central Louisiana. Results were similar with a more responsive groundwater system indicating PRISM data may produce better predictions of streamflow because of a more accurate estimate of rainfall within a watershed or because of a denser grid. Our study implies PRISM is providing a better estimate than CFSR of precipitation within a watershed when rain gauge data are not available, resulting in more accurate simulations of streamflows at the watershed outlet. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Modeling Discharge Rates Using a Coupled Modeled Approach for the Merced River in Yosemite National Park

This study describes the application of the NASA version of the Carnegie‐Ames‐Stanford Approach (CASA) ecosystem model coupled with a surface hydrologic routing scheme previously called the Hydrological Routing Algorithm (HYDRA) to model monthly discharge rates from 2000 to 2007 on the Merced River drainage in Yosemite National Park, California. To assess CASA‐HYDRA's capability to estimate actual water flows in extreme precipitation years, the focus of this study is the 2007 water year, which was very dry, and the 2005 water year, which was a moderately wet year in the historical record. Prior to comparisons to gauge records, CASA‐HYDRA snowmelt algorithms were modified with equations from the U.S. Department of Agriculture Snowmelt‐Runoff Model (SRM), which has been designed to predict daily streamflow in mountain basins where snowmelt is a major runoff factor. Results show that model predictions closely matched monthly flow rates at the Pohono Bridge gauge station (USGS#11266500), with R² = 0.67 and Nash‐Sutcliffe (E) = 0.65. By subdividing the upper Merced River basin into subbasins with high spatial resolution in the gridded modeling approach, we were able to determine which biophysical characteristics in the Sierra differed to the largest degree in extreme low‐flow and high‐flow years. Average elevation and snowpack accumulation were found to be the most important explanatory variables to understand subbasin contributions to monthly discharge rates.