Generation of Ensemble Streamflow Forecasts Using an Enhanced Version of the Snowmelt Runoff Model1

Harshburger, Brian J., Von P. Walden, Karen S. Humes, Brandon C. Moore, Troy R. Blandford, and Albert Rango, 2012. Generation of Ensemble Streamflow Forecasts Using an Enhanced Version of the Snowmelt Runoff Model. Journal of the American Water Resources Association (JAWRA) 48(4): 643‐655. DOI: 10.1111/j.1752‐1688.2012.00642.x Abstract: As water demand increases in the western United States, so does the need for accurate streamflow forecasts. We describe a method for generating ensemble streamflow forecasts (1‐15 days) using an enhanced version of the snowmelt runoff model (SRM). Forecasts are produced for three snowmelt‐dominated basins in Idaho. Model inputs are derived from meteorological forecasts, snow cover imagery, and surface observations from Snowpack Telemetry stations. The model performed well at lead times up to 7 days, but has significant predictability out to 15 days. The timing of peak flow and the streamflow volume are captured well by the model, but the peak‐flow value is typically low. The model performance was assessed by computing the coefficient of determination (R²), percentage of volume difference (Dv%), and a skill score that quantifies the usefulness of the forecasts relative to climatology. The average R² value for the mean ensemble is >0.8 for all three basins for lead times up to seven days. The Dv% is fairly unbiased (within ±10%) out to seven days in two of the basins, but the model underpredicts Dv% in the third. The average skill scores for all basins are >0.6 for lead times up to seven days, indicating that the ensemble model outperforms climatology. These results validate the usefulness of the ensemble forecasting approach for basins of this type, suggesting that the ensemble version of SRM might be applied successfully to other basins in the Intermountain West.     

AI-based runoff simulation based on remote sensing observations: A case study of two river basins in the United States and Canada

Data-driven techniques are used extensively for hydrologic time-series prediction. We created various data-driven models (DDMs) based on machine learning: long short-term memory (LSTM), support vector regression (SVR), extreme learning machines, and an artificial neural network with backpropagation, to define the optimal approach to predicting streamflow time series in the Carson River (California, USA) and Montmorency (Canada) catchments. The moderate resolution imaging spectroradiometer (MODIS) snow-coverage dataset was applied to improve the streamflow estimate. In addition to the DDMs, the conceptual snowmelt runoff model was applied to simulate and forecast daily streamflow. The four main predictor variables, namely snow-coverage (S-C), precipitation (P), maximum temperature (T_max), and minimum temperature (T_min), and their corresponding values for each river basin, were obtained from National Climatic Data Center and National Snow and Ice Data Center to develop the model. The most relevant predictor variable was chosen using the support vector machine-recursive feature elimination feature selection approach. The results show that incorporating the MODIS snow-coverage dataset improves the models' prediction accuracies in the snowmelt-dominated basin. SVR and LSTM exhibited the best performances (root mean square error = 8.63 and 9.80) using monthly and daily snowmelt time series, respectively. In summary, machine learning is a reliable method to forecast runoff as it can be employed in global climate forecasts that require high-volume data processing.     

Improved Bankfull Discharge Prediction Using 2‐Year Recurrence‐Period Discharge1

Abstract:  Knowledge of bankfull discharge (Q_bf) is essential for planners, engineers, geomorphologists, environmentalists, agricultural interests, developments situated on flood prone lands, surface mining and reclamation activities, and others interested in floods and flooding. In conjunction with estimating Q_bf, regionalized bankfull hydraulic geometry relationships, which relate Q_bf and associated channel dimensions (i.e., width, depth, and cross‐section area) to drainage basin area (A_da), are often used. This study seeks to improve upon the common practice of predicting Q_bf using A_da exclusively. Specifically, we hypothesize that predictions of Q_bf can be improved by including estimates of the 2‐year recurrence‐period discharge (Q₂) in regression models for predicting Q_bf. For testing this hypothesis, we used Q_bf estimates from 30 reports containing data for streams that span 34 hydrologic regions in 16 states. Corresponding values of Q₂ and A_da were compiled from flood‐frequency reports and other sources. By comparing statistical measures (i.e., root mean squared error, coefficient of determination, and Akaike’s information criterion), we determined that predicting Q_bf from Q₂ rather than A_da yields consistently better estimates of Q_bf. Other principal findings are (1) data are needed for at least 12 sites in a region for reliable hydraulic geometry model selection and (2) an approximate range of values for Q_bf/Q₂ is 0.10‐3.0.     

Empirical Estimation of Stream Discharge Using Channel Geometry in Low‐Gradient,Sand‐Bed Streams of the Southeastern Plains

Manning's equation is used widely to predict stream discharge (Q) from hydraulic variables when logistics constrain empirical measurements of in‐bank flow events. Uncertainty in Manning's roughness (n_M) is the major source of error in natural channels, and sand‐bed streams pose difficulties because flow resistance is affected by flow‐dependent bed configuration. Our study was designed to develop and validate models for estimating Q from channel geometry easily derived from cross‐sectional surveys and available GIS data. A database was compiled consisting of 484 Q measurements from 75 sand‐bed streams in Alabama, Georgia, South Carolina, North Carolina (Southeastern Plains), and Florida (Southern Coastal Plain), with six New Zealand streams included to develop statistical models to predict Q from hydraulic variables. Model error characteristics were estimated with leave‐one‐site‐out jackknifing. Independent data of 317 Q measurements from 55 Southeastern Plains streams indicated the model (Q = A_cR_H^0.6906S^0.1216; where A_c is the channel area, R_H is the hydraulic radius, and S is the bed slope) best predicted Q, based on Akaike's information criterion and root mean square error. Models also were developed from smaller Q range subsets to explore if subsets increased predictive ability, but error fit statistics suggested that these were not reasonable alternatives to the above equation. Thus, we recommend the above equation for predicting in‐bank Q of unbraided, sandy streams of the Southeastern Plains.     

Decoupling Streamflow Responses to Climate Variability and Land Use/Cover Changes in a Watershed in Northern China

Restored annual streamflow (Q_r) and measured daily streamflow of the Chaohe watershed located in northern China and associated long‐term climate and land use/cover data were used to explore the effects of land use/cover change and climate variability on the streamflow during 1961‐2009. There were no significant changes in annual precipitation (P) and potential evapotranspiration, whereas Q_r decreased significantly by 0.81 mm/yr (p < 0.001) over the study period with a change point in 1999. We used 1961‐1998 as the baseline period (BP) and 1999‐2009 the change period (CP). The mean Q_r during the CP decreased by 39.4 mm compared with that in the BP. From 1979 to 2009, the grassland area declined by 69.6%, and the forest and shrublands increased by 105.4 and 73.1%, respectively. The land use/cover change and climate variability contributed for 58.4 and 41.6% reduction in mean annual Q_r, respectively. Compared with the BP, median and high flows in the CP decreased by 38.8 and up to 75.5%, respectively. The study concludes that large‐scale ecological restoration and watershed management in northern China has greatly decreased water yield and reduced high flows due to the improved land cover by afforestation leading to higher water loss through evapotranspiration. At a large watershed scale, land use/cover change could play as much of an important role as climate variability on water resources.     

Improved Bankfull Channel Geometry Prediction Using Two‐Year Return‐Period Discharge1

He, Laien and Gregory V. Wilkerson, 2011. Improved Bankfull Channel Geometry Prediction Using Two‐Year Return‐Period Discharge. Journal of the American Water Resources Association (JAWRA) 47(6):1298–1316. DOI: 10.1111/j.1752‐1688.2011.00567.x Abstract:  Bankfull discharge (Q_bf) and bankfull channel geometry (i.e., width, W_bf; mean depth, D_bf; and cross‐section area, A_bf) are important design parameters in stream restoration, habitat creation, mined land reclamation, and related projects. The selection of values for these parameters is facilitated by regional curves (regression models in which Q_bf, W_bf, D_bf, and A_bf are predicted as a function of drainage area, A_da). This paper explores the potential for the two‐year return‐period discharge (Q₂) to improve predictions of W_bf, D_bf, and A_bf. Improved predictions are expected because Q₂ estimates integrate the effects of basin drainage area, climate, and geology. For conducting this study, 29 datasets (each representing one hydrologic region) spanning 14 states in the United States were analyzed. We assessed the utility of using Q₂ by comparing statistical measures of regression model performance (e.g., coefficient of determination and Akaike’s information criterion). Compared to using A_da, Q₂ is shown to be a “clearly superior” predictor of W_bf, D_bf, and A_bf, respectively, for 21, 13, and 25% of the datasets. By contrast, A_da yielded a clearly superior model for predicting W_bf, D_bf, and A_bf, respectively, for 0, 0, and 14% of the datasets. Our conclusion is that it alongside with developing conventional regional curves using A_da it is prudent to develop regional curves that use Q₂ as an independent variable because in some cases the resulting model will be superior.     

Performance of the National Water Model in Iowa Using Independent Observations

This paper explores the performance of the analysis‐and‐assimilation configuration of the National Water Model (NWM) v1.0 in Iowa. The NWM assimilates streamflow observations from the United States Geological Survey (USGS), which increases the performance but also limits the available data for model evaluation. In this study, Iowa Flood Center Bridge Sensors (IFCBS) data provided an independent nonassimilated dataset for evaluation analyses. The authors compared NWM outputs for the period between May 2016 and April 2017, with two datasets: USGS streamflow and velocity observations; Stage and streamflow data from IFCBS. The distribution of Spearman rank correlation (r_s), Nash–Sutcliffe efficiency (E), and Kling–Gupta efficiency (KGE) provided quantification of model performance. We found the performance was linked with the spatial scale of the basins. Analysis at USGS gauges showed the strongest performance in large (>10,000 km²) basins (r_s = 0.9, E = 0.9, KGE = 0.8), with some decrease at small (<1,000 km²) basins (r_s = 0.6, E = ?0.25, KGE = ?0.2). Analysis with independent IFCBS observations was used to report performance at large basins (r_s = 0.6, KGE = 0.1) and small basins (r_s = 0.2, KGE = ?0.4). Data assimilation improves simulations at downstream basins. We found differences in the characterization of the model and observed data flow velocity distributions. The authors recommend checking the connection of USGS gauges and NHDPlus reaches for selected locations where performance is weak.     

BASE FLOW RECESSION RATES,LOW FLOWS,AND HYDROLOGIC FEATURES OF SMALL WATERSHEDS IN PENNSYLVANIA,USA1

This paper examines the relationships between measurable watershed hydrologic features, base flow recession rates, and the Q_7,10 low flow statistic (the annual minimum seven‐day average streamflow occurring once every 10 years on average). Base flow recession constants were determined by analyzing hydrograph recession data from 24 small (>130 km²), unregulated watersheds across five major physiographic provinces of Pennsylvania, providing a highly variable dataset. Geomorphic, hydrogeologic, and land use parameters were determined for each watershed. The base flow recession constant was found to be most strongly correlated to drainage density, geologic index, and ruggedness number (watershed slope); however, these three parameters are intercorrelated. Multiple regression models were developed for predicting the recession rate, and it was found that only two parameters, drainage density and hydrologic soil group, were required to obtain good estimates of the recession constant. Equations were also developed to relate the recession rates to Q_7,10 per unit area, and to the Q_7,10/Q₅₀ ratio. Using these equations, estimates of base flow recession rates, Q_7,10, and streamflow reduction under drought conditions can be made for small, ungaged basins across a wide range of physiography.     

APPLICATION OF THE SNOWMELT RUNOFF MODEL USING MULTIPLE-PARAMETER LANDSCAPE ZONES ON THE TOWANDA CREEK BASIN,PENNSYLVANIA1

ABSTRACT: The Snowmelt Runoff Model (SRM) is designed to compute daily stream discharge using satellite snow cover data for a basin divided into elevation zones. For the Towanda Creek basin, a Pennsylvania watershed with relatively little relief, analysis of snow cover images revealed that both elevation and land use affected snow accumulation and melt on the landscape. The distribution of slope and aspect on the watershed was also considered; however, these landscape features were not well correlated with the available snow cover data. SRM streamflow predictions for 1990, 1993 and 1994 snowmelt seasons for the Towanda Creek basin using a combination of elevation and land use zones yielded more precise streamflow estimates than the use of standard elevation zones alone. The use of multiple-parameter zones worked best in non-rain-on-snow conditions such as in 1990 and 1994 seasons where melt was primarily driven by differences in solar radiation. For seasons with major rain-on-snow events such as 1993, only modest improvements were shown since melt was dominated by rainfall energy inputs, condensation and sensible heat convection. Availability of GIS coverages containing satellite snow cover data and other landscape attributes should permit similar reformulation of multiple-parameter watershed zones and improved SRM streamflow predictions on other basins.     

EQUILIBRIUM EVAPORATION AND POSITWE EVAPORATION-SURFACE TEMPERATURE RELATIONSHIPS ACROSS A GRASSLAND1

Sensible (H) and latent (LE) heat fluxes, soil moisture (SM) and surface temperatures (T_s) were analyzed from seven sites at FIFE to evaluate relationships among the spatial variability of evaporative fraction, EF, SM, and the diurnal surface temperature range (T_dr). Intersite correlations between EF and T_dr were significantly negative for regional average soil moisture SM_r < 20 percent, insignificant for 20 < SM_r < 27 percent, and slightly positive for SM_r > 27 percent. Statistical analysis of the pooled correlation coefficient between EF and T_dr for SM_r < 20 percent indicates that it is less than zero at a very high level of significance, while the pooled correlation coefficient for regional SM_r > 27 percent is greater than zero at the 10 percent level. The positive EF:T_dr correlations are attributed to increased surface vapor pressure at warmer sites under nearly potential conditions. These results suggest that to characterize the spatial variability of the energy budget partitioning, a variable representing the thermal response of the site should be included. An important application of these findings relates to modeling the subgrid variability of a region by subdividing the region into a few classes within which surface variables and parameters are assumed invariant. The thermal response of the surface should be included as a variable in defining these classes.     

Observed Changes in Climate and Streamflow in the Upper Rio Grande Basin

Observed streamflow and climate data are used to test the hypothesis that climate change is already affecting Rio Grande streamflow volume derived from snowmelt runoff in ways consistent with model‐based projections of 21st‐Century streamflow. Annual and monthly changes in streamflow volume and surface climate variables on the Upper Rio Grande, near its headwaters in southern Colorado, are assessed for water years 1958–2015. Results indicate winter and spring season temperatures in the basin have increased significantly, April 1 snow water equivalent (SWE) has decreased by approximately 25%, and streamflow has declined slightly in the April–July snowmelt runoff season. Small increases in precipitation have reduced the impact of declining snowpack on trends in streamflow. Changes in the snowpack–runoff relationship are noticeable in hydrographs of mean monthly streamflow, but are most apparent in the changing ratios of precipitation (rain + snow, and SWE) to streamflow and in the declining fraction of runoff attributable to snowpack or winter precipitation. The observed changes provide observational confirmation for model projections of decreasing runoff attributable to snowpack, and demonstrate the decreasing utility of snowpack for predicting subsequent streamflow on a seasonal basis in the Upper Rio Grande Basin.     

Hydrologic Calibration and Validation of SWAT in a Snow‐Dominated Rocky Mountain Watershed,Montana, U.S.A.1

Abstract: The Soil and Water Assessment Tool (SWAT) has been applied successfully in temperate environments but little is known about its performance in the snow‐dominated, forested, mountainous watersheds that provide much of the water supply in western North America. To address this knowledge gap, we configured SWAT to simulate the streamflow of Tenderfoot Creek (TCSWAT). Located in central Montana, TCSWAT represents a high‐elevation watershed with ～85% coniferous forest cover where more than 70% of the annual precipitation falls as snow, and runoff comes primarily from spring snowmelt. Model calibration using four years of measured daily streamflow, temperature, and precipitation data resulted in a relative error (RE) of 2% for annual water yield estimates, and mean paired deviations (Dv) of 36 and 31% and Nash‐Sutcliffe (NS) efficiencies of 0.90 and 0.86 for monthly and daily streamflow, respectively. Model validation was conducted using an additional four years of data and the performance was similar to the calibration period, with RE of 4% for annual water yields, Dv of 43% and 32%, and NS efficiencies of 0.90 and 0.76 for monthly and daily streamflow, respectively. An objective, regression‐based model invalidation procedure also indicated that the model was validated for the overall simulation period. Seasonally, SWAT performed well during the spring and early summer snowmelt runoff period, but was a poor predictor of late summer and winter base flow. The calibrated model was most sensitive to snowmelt parameters, followed in decreasing order of influence by the surface runoff lag, ground water, soil, and SCS Curve Number parameter sets. Model sensitivity to the surface runoff lag parameter reflected the influence of frozen soils on runoff processes. Results indicated that SWAT can provide reasonable predictions of annual, monthly, and daily streamflow from forested montane watersheds, but further model refinements could improve representation of snowmelt runoff processes and performance during the base flow period in this environment.     

Long‐Term Streamflow Response to Climatic Variability in the Loess Plateau,China1

Abstract: The Loess Plateau region in northwestern China has experienced severe water resource shortages due to the combined impacts of climate and land use changes and water resource exploitation during the past decades. This study was designed to examine the impacts of climatic variability on streamflow characteristics of a 12‐km² watershed near Tianshui City, Gansu Province in northwestern China. Statistic analytical methods including Kendall’s trend test and stepwise regression were used to detect trends in relationship between observed streamflow and climatic variables. Sensitivity analysis based on an evapotranspiration model was used to detect quantitative hydrologic sensitivity to climatic variability. We found that precipitation (P), potential evapotranspiration (PET) and streamflow (Q) were not statistically significantly different (p > 0.05) over the study period between 1982 and 2003. Stepwise regression and sensitivity analysis all indicated that P was more influential than PET in affecting annual streamflow, but the similar relationship existed at the monthly scale. The sensitivity of streamflow response to variations of P and PET increased slightly with the increase in watershed dryness (PET/P) as well as the increase in runoff ratio (Q/P). This study concluded that future changes in climate, precipitation in particular, will significantly impact water resources in the Loess Plateau region an area that is already experiencing a decreasing trend in water yield.     

Characterizing Runoff and Water Yield for Headwater Catchments in the Southern Sierra Nevada

In a Mediterranean climate where much of the precipitation falls during winter, snowpacks serve as the primary source of dry season runoff. Increased warming has led to significant changes in hydrology of the western United States. An important question in this context is how to best manage forested catchments for water and other ecosystem services? Answering this basic question requires detailed understanding of hydrologic functioning of these catchments. Here, we depict the differences in hydrologic response of 10 catchments. Size of the study catchments ranges from 50 to 475 ha, and they span between 1,782 and 2,373 m elevation in the rain‐snow transitional zone. Mean annual streamflow ranged from 281 to 408 mm in the low elevation Providence and 436 to 656 mm in the high elevation Bull catchments, resulting in a 49 mm streamflow increase per 100 m (R² = 0.79) elevation gain, despite similar precipitation across the 10 catchments. Although high elevation Bull catchments received significantly more precipitation as snow and thus experienced a delayed melt, this increase in streamflow with elevation was mainly due to a reduction in evapotranspiration (ET) with elevation (45 mm/100 m, R² = 0.65). The reduction in ET was attributed to decline in vegetation density, growing season, and atmospheric demand with increasing elevation. These findings suggest changes in streamflow in response to climate warming may likely depend on how vegetation responds to those changes in climate.     

Streamflow Forecasting Using Singular Value Decomposition and Support Vector Machine for the Upper Rio Grande River Basin

The current study improves streamflow forecast lead‐time by coupling climate information in a data‐driven modeling framework. The spatial–temporal correlation between streamflow and oceanic–atmospheric variability represented by sea surface temperature (SST), 500‐mbar geopotential height (Z₅₀₀), 500‐mbar specific humidity (SH₅₀₀), and 500‐mbar east–west wind (U₅₀₀) of the Pacific and the Atlantic Ocean is obtained through singular value decomposition (SVD). SVD significant regions are weighted using a nonparametric method and utilized as input in a support vector machine (SVM) framework. The Upper Rio Grande River Basin (URGRB) is selected to test the applicability of the proposed model for the period of 1965–2014. The April–August streamflow volume is forecasted using previous year climate variability, creating a lagged relationship of 1–13 months. SVD results showed the streamflow variability was better explained by SST and U₅₀₀ as compared to Z₅₀₀ and SH₅₀₀. The SVM model showed satisfactory forecasting ability with best results achieved using a one‐month lead to forecast the following four‐month period. Overall, the SVM results showed excellent predictive ability with average correlation coefficient of 0.89 and Nash–Sutcliffe efficiency of 0.79. This study contributes toward identifying new SVD significant regions and improving streamflow forecast lead‐time of the URGRB.     

Effects of Land Use and Climate Change on Stream Temperature I: Daily Flow and Stream Temperature Projections

Freshwater mussels (order Unionida) are a highly imperiled group of organisms that are at risk from rising stream temperatures (T). There is a need to understand the potential effects of land use (LU) and climate change (CC) on stream T and have a measure of uncertainty. We used available downscaled climate projections and LU change simulations to simulate the potential effects on average daily stream T from 2020 to 2060. Monte Carlo simulations were run, and a novel technique to analyze results was used to assess changes in hydrologic and stream T response. Simulations of daily mean T were used as input to our stochastic hourly T model. CC effects were on average two orders of magnitude greater than LU impacts on mean daily stream T. LU change affected stream T primarily in headwater streams, on average up to 2.1°C over short durations, and projected CC affected stream T, on average 2.1‐3.3°C by 2060. Daily mean flow and T ratios from Monte Carlo simulations indicated greater variance in the response of streamflow (up to 55%) to LU change than in the response of stream T (up to 9%), and greater variance in headwater stream segments compared to higher order stream segments for both streamflow and T response. Simulations indicated that combined effects of climate and LU change were not additive, suggesting a complex interaction and that forecasting long‐term stream T response requires simulating CC and LU change simultaneously.     

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.     

Analysis of Appropriate Timescales for Water Diversion Permits in Brazil

Agreement on the criteria for granting the right to use water resources between governing bodies represents a significant advance in the process of sharing water use. To aid water resource management agencies in optimizing water use, the impact of using different criteria for permitting water use in the Paracatu river basin, Brazil, was evaluated in this study. The streamflow criteria corresponding to 30 % of the annual Q_7,10 (used by the governing body of Minas Gerais), 70 % of the annual Q₉₅ (used by the governing body of the union), 30 % of the monthly Q_7,10, and 70 % of the monthly Q₉₅ were evaluated. The use of criteria based on the monthly streamflow allows for better management of water use because it allows for greater utilization of this resource in times when there is high water availability and imposes a more realistic restriction during critical periods. Substitution of the annual Q_7,10 for the monthly Q_7,10 significantly increases the streamflow permitted in some months, for example, from December to May. Use of the criterion of 70 % of the annual Q₉₅ involves a high risk of drought in critical months, while the criterion of 70 % of the monthly Q₉₅ minimizes this risk.     

Comparison of the Performance of Statistical Models in Forecasting Monthly Total Dissolved Solids in the Rio Grande1

Abudu, S., J.P. King, Z. Sheng, 2011. Comparison of the Performance of Statistical Models in Forecasting Monthly Total Dissolved Solids in the Rio Grande. Journal of the American Water Resources Association (JAWRA) 48(1): 10‐23. DOI: 10.1111/j.1752‐1688.2011.00587.x Abstract: This paper presents the application of autoregressive integrated moving average (ARIMA), transfer function‐noise (TFN), and artificial neural networks (ANNs) modeling approaches in forecasting monthly total dissolved solids (TDS) of water in the Rio Grande at El Paso, Texas. Predictability analysis was performed between the precipitation, temperature, streamflow rates at the site, releases from upstream reservoirs, and monthly TDS using cross‐correlation statistical tests. The chi‐square test results indicated that the average monthly temperature and precipitation did not show significant predictability on monthly TDS series. The performances of one‐ to three‐month‐ahead model forecasts for the testing period of 1984‐1994 showed that the TFN model that incorporated the streamflow rates at the site and Caballo Reservoir release improved monthly TDS forecasts slightly better than the ARIMA models. Except for one‐month‐ahead forecasts, the ANN models using the streamflow rates at the site as inputs resulted in no significant improvements over the TFN models at two‐month‐ahead and three‐month‐ahead forecasts. For three‐month‐ahead forecasts, the simple ARIMA showed similar performance compared to all other models. The results of this study suggested that simple deseasonalized ARIMA models could be used in one‐ to three‐month‐ahead TDS forecasting at the study site with a simple, explicit model structure and similar model performance as the TFN and ANN models for better water management in the Basin.     

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.