Abstract

A snow data assimilation study was undertaken in which real data were used to update a conceptual model, SNOW-17. The aim of this study is to improve the model’s estimate of snow water equivalent (SWE) by merging the uncertainties associated with meteorological forcing data and SWE observations within the model. This is done with a view to aiding the estimation of snowpack initial conditions for the ultimate objective of streamflow forecasting via a distributed hydrologic model. To provide a test of this methodology, the authors performed experiments at 53 stations in Colorado. In each case the situation of an unobserved location is mimicked, using the data at any given station only for validation; essentially, these are withholding experiments. Both ensembles of model forcing data and assimilated data were derived via interpolation and stochastic modeling of data from surrounding sources. Through a process of cross validation the error for the ensemble of model forcing data and assimilated observations is explicitly estimated. An ensemble square root Kalman filter is applied to perform assimilation on a 5-day cycle. Improvements in the resulting SWE are most evident during the early accumulation season and late melt period. However, the large temporal correlation inherent in a snowpack results in a less than optimal assimilation and the increased skill is marginal. Once this temporal persistence is removed from both model and assimilated observations during the update cycle, a result is produced that is, within the limits of available information, consistently superior to either the model or interpolated observations.

Abstract

This paper describes a flexible method to generate ensemble gridded fields of precipitation in complex terrain. The method is based on locally weighted regression, in which spatial attributes from station locations are used as explanatory variables to predict spatial variability in precipitation. For each time step, regression models are used to estimate the conditional cumulative distribution function (cdf) of precipitation at each grid cell (conditional on daily precipitation totals from a sparse station network), and ensembles are generated by using realizations from correlated random fields to extract values from the gridded precipitation cdfs. Daily high-resolution precipitation ensembles are generated for a 300 km × 300 km section of western Colorado (dx = 2 km) for the period 1980–2003. The ensemble precipitation grids reproduce the climatological precipitation gradients and observed spatial correlation structure. Probabilistic verification shows that the precipitation estimates are reliable, in the sense that there is close agreement between the frequency of occurrence of specific precipitation events in different probability categories and the probability that is estimated from the ensemble. The probabilistic estimates have good discrimination in the sense that the estimated probabilities differ significantly between cases when specific precipitation events occur and when they do not. The method may be improved by merging the gauge-based precipitation ensembles with remotely sensed precipitation estimates from ground-based radar and satellites, or with precipitation and wind fields from numerical weather prediction models. The stochastic modeling framework developed in this study is flexible and can easily accommodate additional modifications and improvements.

Abstract

The timing of snowmelt runoff (SMR) for 84 rivers in the western United States is examined to understand the character of SMR variability and the climate processes that may be driving changes in SMR timing. Results indicate that the timing of SMR for many rivers in the western United States has shifted to earlier in the snowmelt season. This shift occurred as a step change during the mid-1980s in conjunction with a step increase in spring and early-summer atmospheric pressures and temperatures over the western United States. The cause of the step change has not yet been determined.

Abstract

This paper examines an archive containing over 40 years of 8-day atmospheric forecasts over the contiguous United States from the NCEP reanalysis project to assess the possibilities for using medium-range numerical weather prediction model output for predictions of streamflow. This analysis shows the biases in the NCEP forecasts to be quite extreme. In many regions, systematic precipitation biases exceed 100% of the mean, with temperature biases exceeding 3°C. In some locations, biases are even higher. The accuracy of NCEP precipitation and 2-m maximum temperature forecasts is computed by interpolating the NCEP model output for each forecast day to the location of each station in the NWS cooperative network and computing the correlation with station observations. Results show that the accuracy of the NCEP forecasts is rather low in many areas of the country. Most apparent is the generally low skill in precipitation forecasts (particularly in July) and low skill in temperature forecasts in the western United States, the eastern seaboard, and the southern tier of states. These results outline a clear need for additional processing of the NCEP Medium-Range Forecast Model (MRF) output before it is used for hydrologic predictions.

Techniques of model output statistics (MOS) are used in this paper to downscale the NCEP forecasts to station locations. Forecasted atmospheric variables (e.g., total column precipitable water, 2-m air temperature) are used as predictors in a forward screening multiple linear regression model to improve forecasts of precipitation and temperature for stations in the National Weather Service cooperative network. This procedure effectively removes all systematic biases in the raw NCEP precipitation and temperature forecasts. MOS guidance also results in substantial improvements in the accuracy of maximum and minimum temperature forecasts throughout the country. For precipitation, forecast improvements were less impressive. MOS guidance increases the accuracy of precipitation forecasts over the northeastern United States, but overall, the accuracy of MOS-based precipitation forecasts is slightly lower than the raw NCEP forecasts.

Four basins in the United States were chosen as case studies to evaluate the value of MRF output for predictions of streamflow. Streamflow forecasts using MRF output were generated for one rainfall-dominated basin (Alapaha River at Statenville, Georgia) and three snowmelt-dominated basins (Animas River at Durango, Colorado; East Fork of the Carson River near Gardnerville, Nevada; and Cle Elum River near Roslyn, Washington). Hydrologic model output forced with measured-station data were used as “truth” to focus attention on the hydrologic effects of errors in the MRF forecasts. Eight-day streamflow forecasts produced using the MOS-corrected MRF output as input (MOS) were compared with those produced using the climatic Ensemble Streamflow Prediction (ESP) technique. MOS-based streamflow forecasts showed increased skill in the snowmelt-dominated river basins, where daily variations in streamflow are strongly forced by temperature. In contrast, the skill of MOS forecasts in the rainfall-dominated basin (the Alapaha River) were equivalent to the skill of the ESP forecasts. Further improvements in streamflow forecasts require more accurate local-scale forecasts of precipitation and temperature, more accurate specification of basin initial conditions, and more accurate model simulations of streamflow.

Abstract

At least four different modeling studies indicate that variability in snow cover over Asia may modulate atmospheric circulation over the North Pacific Ocean during winter. Here, satellite data on snow extent for east Asia for 1971–95 along with atmospheric fields from the National Centers for Environmental Prediction–National Center for Atmospheric Research reanalysis are used to examine whether the circulation signals seen in model results are actually observed in nature. Anomalies in snow extent over east Asia exhibit a distinct lack of persistence. This suggests that understanding the effects of east Asian snow cover is more germane for short- to medium-range weather forecasting applications than for problems on longer timescales. While it is impossible to attribute cause and effect in the empirical study, analyses of composite fields demonstrate relationships between snow cover extremes and atmospheric circulation downstream remarkably similar to those identified in model results. Positive snow cover extremes in midwinter are associated with a small decrease in air temperatures over the transient snow regions, a stronger east Asian jet, and negative geopotential height anomalies over the North Pacific Ocean. Opposing responses are observed for negative snow cover extremes. Diagnosis of storm track feedbacks shows that the action of high-frequency eddies does not reinforce circulation anomalies in positive snow cover extremes. However, in negative snow cover extremes, there are significant decreases in high-frequency eddy activity over the central North Pacific Ocean, and a corresponding decrease in the mean cyclonic effect of these eddies on the geopotential tendency, contributing to observed positive height anomalies over the North Pacific Ocean. The circulation signals over the North Pacific Ocean are much more pronounced in midwinter (January–February) than in the transitional seasons (November–December and March–April).

Abstract

This paper provides a detailed description of the relationship between spring snow mass in the mountain areas of the western United States and summertime precipitation in the southwestern United States associated with the North American monsoon system and examines the hypothesis that antecedent spring snow mass can modulate monsoon rains through effects on land surface energy balance. Analysis of spring snow water equivalent (SWE) and July–August (JA) precipitation for the period of 1948–97 confirms the inverse snow–monsoon relationship noted in previous studies. Examination of regional difference in SWE–JA precipitation associations shows that although JA precipitation in New Mexico is significantly correlated with SWE over much larger areas than in Arizona, the overall strength of the correlations are just as strong in Arizona as in New Mexico. Results from this study also illustrate that the snow–monsoon relationship is unstable over time. In New Mexico, the relationship is strongest during 1965–92 and is weaker outside that period. By contrast, Arizona shows strongest snow–monsoon associations before 1970. The temporal coincidence between stronger snow–monsoon associations over Arizona and weaker snow–monsoon associations over New Mexico (and vice versa) suggests a common forcing mechanism and that the variations in the strength of snow–monsoon associations are more than just climate noise. There is a need to understand how other factors modulate monsoonal rainfall before realistic predictions of summertime precipitation in the Southwest can be made.

Abstract

This study introduces medium-range meteorological ensemble inputs of temperature and precipitation into the Ensemble Streamflow Prediction component of the National Weather Service River Forecast System (NWSRFS). The Climate Diagnostics Center (CDC) produced a reforecast archive of model forecast runs from a dynamically frozen version of the Medium-Range Forecast (MRF) model. This archive was used to derive statistical relationships between MRF variables and historical basin-average precipitation and temperatures. The latter are used to feed the Ensemble Streamflow Prediction (ESP) component of the NWSRFS. Two sets of ESP reforecasts were produced: A control run based on historically observed temperature and precipitation and an experimental run based on MRF-derived temperature and precipitation. This study found the MRF reforecasts to be generally superior to the control reforecasts, although there were situations when the downscaled MRF output actually degraded the forecast. Forecast improvements were most pronounced during the rising limb of the hydrograph—at this time accurate temperature forecasts improve predictions of the rate of snowmelt. Further improvements in streamflow forecasts at short forecast lead times may be possible by incorporating output from high-resolution regional atmospheric models into the NWSRFS.

Abstract

This study compares methods to incorporate climate information into the National Weather Service River Forecast System (NWSRFS). Three small-to-medium river subbasins following roughly along a longitude in the Colorado River basin with different El Niño–Southern Oscillation signals were chosen as test basins. Historical ensemble forecasts of the spring runoff for each basin were generated using modeled hydrologic states and historical precipitation and temperature observations using the Ensemble Streamflow Prediction (ESP) component of the NWSRFS.

Two general methods for using a climate index (e.g., Niño-3.4) are presented. The first method, post-ESP, uses the climate index to weight ensemble members from ESP. Four different post-ESP weighting schemes are presented. The second method, preadjustment, uses the climate index to modify the temperature and precipitation ensembles used in ESP. Two preadjustment methods are presented. This study shows the distance-sensitive nearest-neighbor post-ESP to be superior to the other post-ESP weighting schemes. Further, for the basins studied, forecasts based on post-ESP techniques outperformed those based on preadjustment techniques.

Abstract

This paper presents a set of analytical tools to evaluate the performance of three land surface models (LSMs) that are used in global climate models (GCMs). Predictions of the fluxes of sensible heat, latent heat, and net CO₂ exchange obtained using process-based LSMs are benchmarked against two statistical models that only use incoming solar radiation, air temperature, and specific humidity as inputs to predict the fluxes. Both are then compared to measured fluxes at several flux stations located on three continents. Parameter sets used for the LSMs include default values used in GCMs for the plant functional type and soil type surrounding each flux station, locally calibrated values, and ensemble sets encompassing combinations of parameters within their respective uncertainty ranges. Performance of the LSMs is found to be generally inferior to that of the statistical models across a wide variety of performance metrics, suggesting that the LSMs underutilize the meteorological information used in their inputs and that model complexity may be hindering accurate prediction. The authors show that model evaluation is purpose specific; good performance in one metric does not guarantee good performance in others. Self-organizing maps are used to divide meteorological “‘forcing space” into distinct regions as a mechanism to identify the conditions under which model bias is greatest. These new techniques will aid modelers to identify the areas of model structure responsible for poor performance.

Abstract

Analysis is performed on the spatiotemporal attributes of North American monsoon system (NAMS) rainfall in the southwestern United States. Trends in the timing and amount of monsoon rainfall for the period 1948–2004 are examined. The timing of the monsoon cycle is tracked by identifying the Julian day when the 10th, 25th, 50th, 75th, and 90th percentiles of the seasonal rainfall total have accumulated. Trends are assessed using the robust Spearman rank correlation analysis and the Kendall–Theil slope estimator. Principal component analysis is used to extract the dominant spatial patterns and these are correlated with antecedent land–ocean–atmosphere variables. Results show a significant delay in the beginning, peak, and closing stages of the monsoon in recent decades. The results also show a decrease in rainfall during July and a corresponding increase in rainfall during August and September. Relating these attributes of the summer rainfall to antecedent winter–spring land and ocean conditions leads to the proposal of the following hypothesis: warmer tropical Pacific sea surface temperatures (SSTs) and cooler northern Pacific SSTs in the antecedent winter–spring leads to wetter than normal conditions over the desert Southwest (and drier than normal conditions over the Pacific Northwest). This enhanced antecedent wetness delays the seasonal heating of the North American continent that is necessary to establish the monsoonal land–ocean temperature gradient. The delay in seasonal warming in turn delays the monsoon initiation, thus reducing rainfall during the typical early monsoon period (July) and increasing rainfall during the later months of the monsoon season (August and September). While the rainfall during the early monsoon appears to be most modulated by antecedent winter–spring Pacific SST patterns, the rainfall in the later part of the monsoon seems to be driven largely by the near-term SST conditions surrounding the monsoon region along the coast of California and the Gulf of California. The role of antecedent land and ocean conditions in modulating the following summer monsoon appears to be quite significant. This enhances the prospects for long-lead forecasts of monsoon rainfall over the southwestern United States, which could have significant implications for water resources planning and management in this water-scarce region.