Precipitation Extremes and Flood Frequency in a Changing Climate in Southeastern Virginia

Despite the advances in climate change modeling, extreme events pose a challenge to develop approaches that are relevant for urban stormwater infrastructure designs and best management practices. The study first investigates the statistical methods applied to the land‐based daily precipitation series acquired from the Global Historical Climatology Network‐Daily (GHCN‐D). Additional analysis was carried out on the simulated Multivariate Adaptive Constructed Analogs (MACA)‐based downscaled daily extreme precipitation of 15 General Circulation Models and Weather Research and Forecasting‐based hourly extreme precipitation of North American Regional Reanalysis to discern the return period of 24‐hr and 48‐hr events. We infer that the GHCN‐D and MACA‐based precipitation reveals increasing trends in annual and seasonal extreme daily precipitation. Both BCC‐CSM1‐1‐m and GFDL‐ESM2M models revealed that the magnitude and frequency of extreme precipitation events are projected to increase between 2016 and 2099. We conclude that the future scenarios show an increase in magnitudes of extreme precipitation up to three times across southeastern Virginia resulting in increased discharge rates at selected gauge locations. The depth‐duration‐frequency curve predicted an increase of 2–3 times in 24‐ and 48‐h precipitation intensity, higher peaks, and indicated an increase of up to 50% in flood magnitude in future scenarios.     

Trends in Raindrop Kinetic Energy with Modeled Climate Warming in the Lake Tahoe Basin

Two means by which climate change may increase surface soil erosion in mountainous terrain are: (1) increasing the frequency of extreme rainfall events and (2) decreasing the duration of snow cover on bare soil. We used output from four general circulation models (GCMs) and two greenhouse gas trajectories to produce a suite of hydrologic variables at a daily time‐step for historic and projected 21st Century conditions. We statistically disaggregated the daily rainfall to hourly, using hourly rainfall from a network of nine weather stations in the Tahoe Basin, and filtered out rain falling on a snowpack. We applied published equations to convert hourly intensity to raindrop kinetic energy (KE) for each day and grid cell in the Basin, averaged across grid cells, and created time series of total annual and maximum annual hourly kinetic energy (TKE and MKE) on snow‐free ground. Using the Generalized Extreme Value distribution, we calculated the significance of long‐term trends in KE on snow‐free ground, and estimated energy levels for return periods of 2, 20, and 100 years. We then detrended the snowpack data and compared the resulting trends in KE with the trends resulting from changes in both rainfall energy and snowpack under two GCMs. Principal findings include (1) upward trends in MKE, (2) stronger upward trends in TKE; and (3) an effect of increasing rainfall intensities on KE in some cases, and a strong effect of reduced snowpack in all cases examined.     

Evaluating and Extending CLIGEN Precipitation Generation for the Loess Plateau of China1

Abstract: Climate generator (CLIGEN) is widely used in the United States to generate long‐term climate scenarios for use with agricultural systems models. Its applicability needs to be evaluated for use in a new region or climate. The objectives were to: (1) evaluate the reproducibility of the latest version of CLIGEN v5.22564 in generating daily, monthly, and yearly precipitation depths at 12 stations, as well as storm patterns including storm duration (D), relative peak intensity (ip), and peak intensity (rp) at 10 stations dispersed across the Loess Plateau and (2) test whether an exponential distribution for generating D and a distribution‐free approach for inducing desired rank correlation between precipitation depth and D can improve storm pattern generations. Mean absolute relative errors (MAREs) for simulating daily, monthly, annual, and annual maximum daily precipitation depth across all 12 stations were 3.5, 1.7, 1.7, and 5.0% for the mean and 5.0, 4.5, 13.0, and 13.6% for the standard deviations (SD), respectively. The model reproduced the distributions of monthly and annual precipitation depths well (p > 0.3), but the distribution of daily precipitation depth was less well produced. The first‐order, two‐state Markov chain algorithm was adequate for generating precipitation occurrence for the Loess Plateau of China; however, it underpredicted the longest dry periods. The CLIGEN‐generated storm patterns poorly. It underpredicted mean and SD of D for storms ≥10 mm by ?60.4 and ?72.6%, respectively. Compared with D, ip, and rp were slightly better reproduced. The MAREs of mean and SD were 21.0 and 52.1% for ip, and 31.2 and 55.2% for rp, respectively. When an exponential distribution was used to generate D, MAREs were reduced to 2.6% for the mean and 7.8% for the SD. However, ip estimation became much worse with MAREs being 128.9% for the mean and 241.1% for the SD. Overall, storm pattern generation needs improvement. For better storm pattern generation for the region, precipitation depth, D, and rp may be generated correlatively using Copula methods.     

On the Role of Spatial,Temporal, and Climatic Forces on Stream Sediment Loading from Rural and Urban Ecosystems

In this study, we characterize the greatest sediment loading events by their sediment delivery behavior; dominant climate, watershed, and antecedent conditions; and their seasonal distribution for rural and urban land uses. The study area is Paradise Creek Watershed, a mixed land use watershed in northern Idaho dominated by saturation excess processes in the upstream rural area and infiltration excess in the downstream urban area. We analyzed 12 years of continuous streamflow, precipitation, and watershed data at two monitoring stations. We identified 137 sediment loading events in the upstream rural section of the watershed and 191 events in the downstream urban section. During the majority of these events conditions were transport limited and the sediment flush occurred early in the event, generally in the first 20% of elapsed event time. Statistical analysis including two dozen explanatory variables showed peak discharge, event duration, and antecedent baseflow explained most of the variation in event sediment load at both stations and for the watershed as a whole (R² = 0.73‐0.78). In the rural area, saturated soils combined with spring snowmelt in March led to the greatest loading events. The urban area load contribution peaked in January, which could be a re‐suspension of streambed sediments from the previous water year. Throughout the study period, one event contributed, on average, 33% of the annual sediment load but only accounted for 2% of the time in a year.     

Bankfull Regional Curves for the Alleghany Plateau/Valley and Ridge,Piedmont, and Coastal Plain Regions of Maryland

Regional curves are empirical relationships that can help identify the bankfull stage in ungaged watersheds and aid in designing the riffle dimension in stream restoration projects. Bankfull regional curves were developed from gage stations with drainage areas less than 102 mi² (264.2 km²) for the Alleghany Plateau/Valley and Ridge (AP/VR), Piedmont, and Coastal Plain regions of Maryland. The AP/VR regions were combined into one region for this project. These curves relate bankfull discharge, cross‐sectional area, width, and mean depth to drainage area within the same hydro‐physiographic region (region with similar rainfall/runoff relationship). The bankfull discharge curve for the Coastal Plain region was further subdivided into the Western Coastal Plain (WCP) and Eastern Coastal Plain (ECP) region due to differences in topography and runoff. Results show that the Maryland Piedmont yields the highest bankfull discharge rate per unit drainage area, followed by the AP/VR, WCP, and ECP. Likewise, the Coastal Plain and AP/VR streams have less bankfull cross‐sectional area per unit drainage area than the Piedmont. The average bankfull discharge return interval across the three hydro‐physiographic regions was 1.4 years. The Maryland regional curves were compared to other curves in the eastern United States. The average bankfull discharge return interval for the other studies ranged from 1.1 to 1.8 years.     

Long‐Term Trends in Streamflow and Precipitation in Northwest California and Southwest Oregon, 1953‐2012

Using nonparametric Mann‐Kendall tests, we assessed long‐term (1953‐2012) trends in streamflow and precipitation in Northern California and Southern Oregon at 26 sites regulated by dams and 41 “unregulated” sites. Few (9%) sites had significant decreasing trends in annual precipitation, but September precipitation declined at 70% of sites. Site characteristics such as runoff type (groundwater, snow, or rain) and dam regulation influenced streamflow trends. Decreasing streamflow trends outnumbered increasing trends for most months except at regulated sites for May‐September. Summer (July‐September) streamflow declined at many sites, including 73% of unregulated sites in September. Applying a LOESS regression model of antecedent precipitation vs. average monthly streamflow, we evaluated the underlying streamflow trend caused by factors other than precipitation. Decreasing trends in precipitation‐adjusted streamflow substantially outnumbered increasing trends for most months. As with streamflow, groundwater‐dominated sites had a greater percent of declining trends in precipitation‐adjusted streamflow than other runoff types. The most pristine surface‐runoff‐dominated watersheds within the study area showed no decreases in precipitation‐adjusted streamflow during the summer months. These results suggest that streamflow decreases at other sites were likely due to more increased human withdrawals and vegetation changes than to climate factors other than precipitation quantity.     

Groundwater Level Trends and Drivers in Two Northern New England Glacial Aquifers

We evaluated long‐term trends and predictors of groundwater levels by month from two well‐studied northern New England forested headwater glacial aquifers: Sleepers River, Vermont, 44 wells, 1992‐2013; and Hubbard Brook, New Hampshire, 15 wells, 1979‐2004. Based on Kendall Tau tests with Sen slope determination, a surprising number of well‐month combinations had negative trends (decreasing water levels) over the respective periods. Sleepers River had slightly more positive than negative trends overall, but among the significant trends (p < 0.1), negative trends dominated 67 to 40. At Hubbard Brook, negative trends outnumbered positive trends by a nearly 2:1 margin and all seven of the significant trends were negative. The negative trends occurred despite generally increasing trends in monthly and annual precipitation. This counterintuitive pattern may be a result of increased precipitation intensity causing higher runoff at the expense of recharge, such that evapotranspiration demand draws down groundwater storage. We evaluated predictors of month‐end water levels by multiple regression of 18 variables related to climate, streamflow, snowpack, and prior month water level. Monthly flow and prior month water level were the two strongest predictors for most months at both sites. The predictive power and ready availability of streamflow data can be exploited as a proxy to extend limited groundwater level records over longer time periods.     

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.     

SPATIAL VARIABILITY IN WATER‐BALANCE MODEL PERFORMANCE IN THE CONTERMINOUS UNITED STATES1

ABSTRACT: A monthly water‐balance (WB) model was tested in 44 river basins from diverse physiographic and climatic regions across the conterminous United States (U.S.). The WB model includes the concepts of climatic water supply and climatic water demand, seasonality in climatic water supply and demand, and soil‐moisture storage. Exhaustive search techniques were employed to determine the optimal set of precipitation and temperature stations, and the optimal set of WB model parameters to use for each basin. It was found that the WB model worked best for basins with: (1) a mean elevation less than 450 meters or greater than 2000 meters, and/or (2) monthly runoff that is greater than 5 millimeters (mm) more than 80 percent of the time. In a separate analysis, a multiple linear regression (MLR) was computed using the adjusted R‐square values obtained by comparing measured and estimated monthly runoff of the original 44 river basins as the dependent variable, and combinations of various independent variables [streamflow gauge latitude, longitude, and elevation; basin area, the long‐term mean and standard deviation of annual precipitation; temperature and runoff; and low‐flow statistics (i.e., the percentage of months with monthly runoff that is less than 5 mm)]. Results from the MLR study showed that the reliability of a WB model for application in a specific region can be estimated from mean basin elevation and the percentage of months with gauged runoff less than 5 mm. The MLR equations were subsequently used to estimate adjusted R‐square values for 1,646 gauging stations across the conterminous U.S. Results of this study indicate that WB models can be used reliably to estimate monthly runoff in the eastern U.S., mountainous areas of the western U.S., and the Pacific Northwest. Applications of monthly WB models in the central U.S. can lead to uncertain estimates of runoff.     

Evidence for Changing Flood Risk in New England Since the Late 20th Century1

Abstract: Long‐term flow records for watersheds with minimal human influence have shown trends in recent decades toward increasing streamflow at regional and national scales, especially for low flow quantiles like the annual minimum and annual median flows. Trends for high flow quantiles are less clear, despite recent research showing increased precipitation in the conterminous United States over the last century that has been brought about primarily by an increased frequency and intensity of events in the upper 10th percentile of the daily precipitation distribution – particularly in the Northeast. This study investigates trends in 28 long‐term annual flood series for New England watersheds with dominantly natural streamflow. The flood series are an average of 75 years in length and are continuous through 2006. Twenty‐five series show upward trends via the nonparametric Mann‐Kendall test, 40% (10) of which are statistically significant (p < 0.1). Moreover, an average standardized departures series for 23 of the study gages indicates that increasing flood magnitudes in New England occurred as a step change around 1970. The timing of this is broadly synchronous with a phase change in the low frequency variability of the North Atlantic Oscillation, a prominent upper atmospheric circulation pattern that is known to effect climate variability along the United States east coast. Identifiable hydroclimatic shifts should be considered when the affected flow records are used for flood frequency analyses. Special treatment of the flood series can improve the analyses and provide better estimates of flood magnitudes and frequencies under the prevailing hydroclimatic condition.     

Probabilistic Streamflow Generation Model for Data Sparse Arid Watersheds1

Abstract: The authors present a model that generates streamflow for ephemeral arid streams. The model consists of a stochastic hourly precipitation point process model and a conceptual model that transforms precipitation into flow. It was applied to the Santa Cruz River at the border crossing from Mexico into Southern Arizona. The model was constructed for four different seasons and three categories of inter‐annual variability for the wet seasons of summer and winter. The drainage area is ungauged and precipitation information was inferred from a precipitation gauge downstream. The precipitation gauge record was evaluated against simulated precipitation from a mesoscale numerical weather prediction model, and was found to be the representative of the regional precipitation variability. The flow generation was found to reproduce the variability in the observed record at the daily, seasonal and annual time scales, and it is suitable for use in planning studies for the study site.     

Statistical Models to Predict and Assess Spatial and Temporal Low‐Flow Variability in New England Rivers and Streams

In the northern hemisphere, summer low flows are a key attribute defining both quantity and quality of aquatic habitat. I developed one set of models for New England streams/rivers predicting July/August median flows averaged across 1985–2015 as a function of weather, slope, % imperviousness, watershed storage, glacial geology, and soils. These models performed better than most United States Geological Survey models for summer flows developed at a statewide scale. I developed a second set of models predicting interannual differences in summer flows as a function of differences in air temperature, precipitation, the North Atlantic Oscillation (NAO) index, and lagged NAO. Use of difference equations eliminated the need for transformations and accounted for serial autocorrelations at lag 1. The models were used in sequence to estimate time series for monthly low flows and for two derived flow metrics (tenth percentile [Q10] and minimum 3‐in‐5 year average flows). The first metric is commonly used in assessing risk to low‐flow conditions over time, while the second has been correlated with increased probability of localized extinctions for brook trout. The flow metrics showed increasing trends across most of New England for 1985–2015. However, application of summer flow models with average and extreme climate projections to the Taunton River, Massachusetts, a sensitive watershed undergoing rapid development, projected that low‐flow metrics will decrease over the next 50 years.     

Seasonal and Regional Patterns in Performance for a Baltic Sea Drainage Basin Hydrologic Model

This study evaluates the ability of the Catchment SIMulation (CSIM) hydrologic model to describe seasonal and regional variations in river discharge over the entire Baltic Sea drainage basin (BSDB) based on 31 years of monthly simulation from 1970 through 2000. To date, the model has been successfully applied to simulate annual fluxes of water from the catchments draining into the Baltic Sea. Here, we consider spatiotemporal bias in the distribution of monthly modeling errors across the BSDB since it could potentially reduce the fidelity of predictions and negatively affect the design and implementation of land‐management strategies. Within the period considered, the CSIM model accurately reproduced the annual flows across the BSDB; however, it tended to underpredict the proportion of discharge during high‐flow periods (i.e., spring months) and overpredict during the summer low flow periods. While the general overpredictions during summer periods are spread across all the subbasins of the BSDB, the underprediction during spring periods is seen largely in the northern regions. By implementing a genetic algorithm calibration procedure and/or seasonal parameterization of subsurface water flows for a subset of the catchments modeled, we demonstrate that it is possible to improve the model performance albeit at the cost of increased parameterization and potential loss of parsimony.     

Effect of Climate Warming and Drought on Urban Stream Temperatures in a Semiarid,Seasonal System

Stream temperatures are key indicators for aquatic ecosystem health, and are of particular concern in highly seasonal, water‐limited regions such as California that provide sensitive habitat for cold‐water species. Yet in many of these critical regions, the combined impacts of a warmer climate and urbanization on stream temperatures have not been systematically studied. We examined recent changes in air temperature and precipitation, including during the recent extreme drought, and compared the stream temperature responses of urban and nonurban streams under four climatic conditions and the 2008–2018 period. Metrics included changes in the magnitude and timing of stream temperatures, and the frequency of exceedance of ecologically relevant thresholds. Our results showed that minimum and average daily air temperatures in the region have increased by >1°C over the past 20 years, warming both urban and nonurban streams. Stream temperatures under drought warmed most (1°C–2°C) in late spring and early fall, effectively lengthening the summer warm season. The frequency of occurrence of periods of elevated stream temperatures was greater during warm climate conditions for both urban and nonurban streams, but urban streams experienced extreme conditions 1.5–2 times as often as nonurban streams. Our findings underscore that systematically monitoring and managing urban stream temperatures under climate change and drought is critically needed for seasonal, water‐limited urban systems.     

PRECIPITATION CHANGES FROM 1956 TO 1996 ON THE WALNUT GULCH EXPERIMENTAL WATERSHED1

ABSTRACT: The climate of Southern Arizona is dominated by summer precipitation, which accounts for over 60 percent of the annual total. Summer and non‐summer precipitation data from the USDA‐ARS Walnut Gulch Experimental Watershed are analyzed to identify trends in precipitation characteristics from 1956 to 1996. During this period, annual precipitation increased. The annual precipitation increase can be attributed to an increase in precipitation during non‐summer months, and is paralleled by an increase in the proportion of annual precipitation contributed during non‐summer months. This finding is consistent with previously reported increases in non‐summer precipitation in the southwestern United States. Detailed event data were analyzed to provide insight into the characteristics of precipitation events during this time period. Precipitation event data were characterized based on the number of events, event precipitation amount, 30‐minute event intensity, and event duration. The trend in non‐summer precipitation appears to be a result of increased event frequency since the number of events increased during nonsummer months, although the average amount per event, average event intensity, and average event duration did not. During the summer “monsoon” season, the frequency of recorded precipitation events increased but the average precipitation amount per event decreased. Knowledge of precipitation trends and the characteristics of events that make up a precipitation time series is a critical first step in understanding and managing water resources in semiarid ecosystems.     

Daily Updating of Operational Statistical Seasonal Water Supply Forecasts for the western U.S.1

Abstract: Official seasonal water supply outlooks for the western United States are typically produced once per month from January through June. The Natural Resources Conservation Service has developed a new outlook product that allows the automated production and delivery of this type of forecast year‐round and with a daily update frequency. Daily snow water equivalent and water year‐to‐date precipitation data from multiple SNOTEL stations are combined using a statistical forecasting technique (“Z‐Score Regression”) to predict seasonal streamflow volume. The skill of these forecasts vs. lead‐time is comparable to the official published outlooks. The new product matches the intra‐monthly trends in the official forecasts until the target period is partly in the past, when the official forecasts begin to use information about observed streamflows to date. Geographically, the patterns of skill also match the official outlooks, with highest skill in Idaho and southern Colorado and lowest skill in the Colorado Front Range, eastern New Mexico, and eastern Montana. The direct and frequent delivery of objective guidance to users is a significant new development in the operational hydrologic seasonal forecasting community.     

Temporal Changes in Streamflow and Attribution of Changes to Climate and Landuse in Wisconsin Watersheds

Previous historic trends analyses on 21st Century hydrologic data in the United States generally focus on annual flow statistics and have continued to use USGS hydro‐climatic data network (HCDN) stations, although post‐1988 diversions and runoff regulations are not reflected in the HCDN. Using a more recent dataset, Geospatial Attributes of Gages for Evaluating Streamflow, version II (GAGES II), compiled by Falcone (2012), which includes more watersheds with reference conditions, a comprehensive analysis of changes in seasonal, and annual streamflow in Wisconsin watersheds is demonstrated. Given the pronounced influence of seasonal hydrology in Wisconsin watersheds, the objective of this study is to elucidate the nature of temporal (annual, seasonal, and monthly) changes in runoff. Considerable temporal and regional variability was found in annual and seasonal streamflow changes between the two historic periods 1951‐1980 and 1981‐2010 considered in the study. For example, the northern watersheds show relatively small changes in streamflow discharge ranging from ?6.0 to 4.2%, while the southern watersheds show relatively large increases in streamflow discharge ranging from 13.1 to 18.2%. To apportion streamflow changes to climate and nonclimatic factors, a method based on potential evapotranspiration changes is demonstrated. Results show that nonclimatic factors account for more than 60% of changes in annual runoff in Wisconsin watersheds considered in the study.     

Spatial and Temporal Patterns of Low Streamflow and Precipitation Changes in the Chesapeake Bay Watershed

Spatial and temporal patterns in low streamflows were investigated for 183 streamgages located in the Chesapeake Bay Watershed for the period 1939–2013. Metrics that represent different aspects of the frequency and magnitude of low streamflows were examined for trends: (1) the annual time series of seven‐day average minimum streamflow, (2) the scaled average deficit at or below the 2% mean daily streamflow value relative to a base period between 1939 and 1970, and (3) the annual number of days below the 2% threshold. Trends in these statistics showed spatial cohesion, with increasing low streamflow volume at streamgages located in the northern uplands of the Chesapeake Bay Watershed and decreasing low streamflow volume at streamgages in the southern part of the watershed. For a small subset of streamgages (12%), conflicting trend patterns were observed between the seven‐day average minimum streamflow and the below‐threshold time series and these appear to be related to upstream diversions or the influence of reservoir‐influenced streamflows in their contributing watersheds. Using multivariate classification techniques, mean annual precipitation and fraction of precipitation falling as snow appear to be broad controls of increasing and decreasing low‐flow trends. Further investigation of seasonal precipitation patterns shows summer rainfall patterns, driven by the Atlantic Multidecadal Oscillation, as the main driver of low streamflows in the Chesapeake Bay Watershed.     

Forecasting the Probability of Future Groundwater Levels Declining Below Specified Low Thresholds in the Conterminous U.S.

We present a logistic regression approach for forecasting the probability of future groundwater levels declining or maintaining below specific groundwater‐level thresholds. We tested our approach on 102 groundwater wells in different climatic regions and aquifers of the United States that are part of the U.S. Geological Survey Groundwater Climate Response Network. We evaluated the importance of current groundwater levels, precipitation, streamflow, seasonal variability, Palmer Drought Severity Index, and atmosphere/ocean indices for developing the logistic regression equations. Several diagnostics of model fit were used to evaluate the regression equations, including testing of autocorrelation of residuals, goodness‐of‐fit metrics, and bootstrap validation testing. The probabilistic predictions were most successful at wells with high persistence (low month‐to‐month variability) in their groundwater records and at wells where the groundwater level remained below the defined low threshold for sustained periods (generally three months or longer). The model fit was weakest at wells with strong seasonal variability in levels and with shorter duration low‐threshold events. We identified challenges in deriving probabilistic‐forecasting models and possible approaches for addressing those challenges.     

Effects of Urbanization on Flow Duration and Stream Flashiness: A Case Study of Puget Sound Streams,Western Washington,USA

The overall influence of urbanization on how flows of different frequency might change over time, while important in hydrologic design, remains imprecisely known. In this study, we investigate the effects of urbanization on flow duration curves (FDCs) and flow variability through a case study of eight watersheds that underwent different amounts of growth, in the Puget Sound region in Western Washington State, United States. We computed annual FDCs from flow records spanning 1960‐2010 and, after accounting for the effects of precipitation, we conducted statistical trend analyses on flow metrics to quantify how key FDC percentiles changed with time in response to urbanization. In the urban watersheds, the entire FDC tended to increase in magnitude of flow, especially the 95th‐99th percentile of the daily mean flow series, which increased by an average of 43%. Stream flashiness in urban watersheds was found to increase by an average of 70%. The increases in FDC magnitude and flashiness in urbanizing watersheds are most likely a result of increasing watershed imperviousness and altered hydrologic routing. Rural watersheds were found to have decreasing FDC magnitude over the same time period, which is possibly due to anthropogenic extractions of groundwater, and increasing stream flashiness, which is likely a result of reductions in base flow and increasing precipitation intensity and variability.