The identification of hydrological indices for the characterization of macroinvertebrate community response to flow regime variability

Abstract

The importance of flow regime variability for maintaining ecological functioning and integrity of river ecosystems has been firmly established in both natural and anthropogenically modified systems. River flow regimes across lowland catchments in eastern England are examined using 47 variables, including those derived using the Indicators of Hydrologic Alteration (IHA) software. A principal component analysis method was used to identify redundant hydrological variables and those that best characterized the hydrological series (1986–2005). A small number of variables (<6) characterized up to 95% of the statistical variability in the flow series. The hydrological processes and conditions that the variables represent were found to be significant in structuring the in-stream macroinvertebrate community Lotic-invertebrate Index for Flow Evaluation (LIFE) scores at both the family and species levels. However, hydrological variables only account for a relatively small proportion of the total ecological variability (typically <10%). The research indicates that a range of other factors, including channel morphology and anthropogenic modification of in-stream habitats, structure riverine macroinvertebrate communities in addition to hydrology. These factors need to be considered in future environmental flow studies to enable the characterization of baseline/reference conditions for management and restoration purposes.

Editor Z.W. Kundzewicz; Guest editor M. Acreman

Citation Worrall, T.P., Dunbar, M.J., Extence, C.A., Laizé, C.L.R., Monk, W.A., and Wood, P.J., 2014. The identification of hydrological indices for the characterization of macroinvertebrate community response to flow regime variability. Hydrological Sciences Journal, 59 (3–4), 645–658.     

The response of flow duration curves to afforestation

The hydrologic effect of replacing pasture or other short crops with trees is reasonably well understood on a mean annual basis. The impact on flow regime, as described by the annual flow duration curve (FDC) is less certain. A method to assess the impact of plantation establishment on FDCs was developed. The starting point for the analyses was the assumption that rainfall and vegetation age are the principal drivers of evapotranspiration. A key objective was to remove the variability in the rainfall signal, leaving changes in streamflow solely attributable to the evapotranspiration of the plantation. A method was developed to (1) fit a model to the observed annual time series of FDC percentiles; i.e. 10th percentile for each year of record with annual rainfall and plantation age as parameters, (2) replace the annual rainfall variation with the long term mean to obtain climate adjusted FDCs, and (3) quantify changes in FDC percentiles as plantations age. Data from 10 catchments from Australia, South Africa and New Zealand were used. The model was able to represent flow variation for the majority of percentiles at eight of the 10 catchments, particularly for the 10–50th percentiles. The adjusted FDCs revealed variable patterns in flow reductions with two types of responses (groups) being identified. Group 1 catchments show a substantial increase in the number of zero flow days, with low flows being more affected than high flows. Group 2 catchments show a more uniform reduction in flows across all percentiles. The differences may be partly explained by storage characteristics. The modelled flow reductions were in accord with published results of paired catchment experiments. An additional analysis was performed to characterise the impact of afforestation on the number of zero flow days (N_zero) for the catchments in group 1. This model performed particularly well, and when adjusted for climate, indicated a significant increase in N_zero. The zero flow day method could be used to determine change in the occurrence of any given flow in response to afforestation. The methods used in this study proved satisfactory in removing the rainfall variability, and have added useful insight into the hydrologic impacts of plantation establishment. This approach provides a methodology for understanding catchment response to afforestation, where paired catchment data is not available.     

Channel movement and erosion response to rainfall variability in southeast Australia

Intrinsic and extrinsic forces on the catchment and stream channel network drive morphological change. Separating individual forcings is difficult given the complexity of such nonlinear systems. Here a modelling approach is used to investigate the sensitivity of channel position and movement under a series of realistic rainfall scenarios for a catchment in southeastern Australia. The results demonstrate the sensitivity of the catchment to different rainfall patterns and how relatively small changes in rainfall can lead to much larger sediment outputs revealing sensitivity to subtle changes in climate. Channel movement occurs as an avulsion. This is the first time such a process has been observed and modelled in an ephemeral stream environment and demonstrates fluvial geomorphic change at human time scales. Human intervention by rock lining channels was demonstrated to prevent the movement of the main channel. Overall the CAESAR landscape evolution and erosion model used in this study is able to replicate both erosion rates and the variation in past channel movement. The modelling suggests that any landscape change is based on both internal and external forcing and that landscape history also plays a significant role. Here, we demonstrate the potential to quantify many of the nonlinearities and thresholds in soil‐mantled catchments using a landscape evolution model. Copyright © 2011 John Wiley & Sons, Ltd.     

Subtidal flow and its variability at the entrance to a subtropical lagoon

Spatial and temporal variability of the subtidal exchange flow at West Pass, an inlet at the entrance to a subtropical lagoon (St. Andrew Bay, Florida), was studied using moored and towed current velocity profiles and hydrographic data. Towed and hydrographic measurements were captured over one diurnal tidal cycle to determine intratidal and spatial changes in flow. Hydrographic profiles over the tidal cycle showed that tidal straining modified density stratification asymmetrically, thus setting up the observed mean flow within the inlet. During the towed survey, the inlet's mean flow had a two-layer exchange structure that was moderately frictional and weakly influenced by Coriolis accelerations. Moored current profiles revealed the additional contribution to the dynamics from centrifugal accelerations. Along channel residual flows changed between unidirectional and exchange flow, depending on the forcing from the along-estuary wind stress and, to a lesser extent, the spring–neap tidal cycle. Increases in vertical shear in the along channel subtidal flow coincided with neap tides and rain pulses. Lateral subtidal flows showed the influence on the dynamics of centrifugal accelerations through a well-developed two-layer structure modulated in magnitude by the spring–neap tidal cycle.     

Vegetation response to rainfall seasonality and interannual variability in tropical dry forests

Projected changes in rainfall seasonality and interannual variability are expected to have severe impacts on arid and semi‐arid tropical vegetation, which is characterized by a fine‐tuned adaptation to extreme rainfall seasonality. To study the response of these ecosystems and the related changes in hydrological processes to changes in the amount and seasonality of rainfall, we focused on the caatinga biome, the typical seasonally dry forest in semi‐arid Northeast Brazil. We selected four sites across a gradient of rainfall amount and seasonality and analysed daily rainfall and biweekly Normalized Difference Vegetation Index (NDVI) data for hydrological years 2000 to 2014. Rainfall seasonal and interannual statistics were characterized by recently proposed metrics describing duration, timing and intensity of the wet season and compared to similar metrics of NDVI time series. The results show that the caatinga tends to have a more stable response with longer and less variable growing seasons (3.1 ± 0.1 months) compared to the duration wet seasons (2.0 ± 0.5 months). The ecosystem ability to buffer the interannual variability of rainfall is also evidenced by the stability in the timing of the growing season compared to the wet season, which results in variable delays (ranging from 0 to 2 months) between the peak of the rainfall season and the production of leaves by the ecosystem. The analyses show that the shape and size of the related hysteresis loops in the rainfall–NDVI relations are linked to the buffering effects of soil moisture and plant growth dynamics. Finally, model projections of vegetation response to different rainfall scenarios reveal the existence of a maximum in ecosystem productivity at intermediate levels of rainfall seasonality, suggesting a possible trade‐off in the effects of intensity (i.e. amount) and duration of the wet season on vegetation growth and related soil moisture dynamics and transpiration rates. Copyright © 2016 John Wiley & Sons, Ltd.     

Subtidal variability of transverse flow off a cape

Time series of velocity profiles at two Chesapeake Bay entrance sites were used to characterize the subtidal variability of transverse flows off a cape. A shallow sampling site was located near Cape Henry over 6 m of water and separated from a deep site, 20 m deep, by a distance of 4 km. The velocity profiles showed that wind-induced subtidal variations in general masked curvature effects (centrifugal accelerations) that may produce secondary circulation associated with tidal flow around a cape. Such secondary circulation, consisting of flow away from the cape at surface and toward the cape at depth, was observed only during periods of weak winds. Most of the time, transverse flows were unidirectional throughout the water column and moved in opposite directions at the two sites examined. This caused convergence of transverse flow between the two sites under the influence of northerly winds and divergence of transverse flow with southwesterly winds. In addition to unidirectional and curvature-induced secondary flows, other modes of subtidal variability consisted of (1) two-layered responses with surface flow toward the cape, and (2) three-layered responses. These two- and three-layered structures were observed more frequently at the deep site, where greater stratification is expected, than at the shallow site.Responsible Editor: Iris Grabemann     

Regional response of dengue fever epidemics to interannual variation and related climate variability

Stochastic Environmental Research and Risk Assessment - Dengue is a major international public health concern and one of the most important vector-borne diseases. The purpose of this article is to...     

Stochastic analysis of spatial variability in unconfined groundwater flow

In this paper, spatial variability in steady one-dimensional unconfined groundwater flow in heterogeneous formations is investigated. An approach to deriving the variance of the hydraulic head is developed using the nonlinear filter theory. The nonlinear governing equation describing the one-dimensional unconfined groundwater flow is decomposed into three linear partial differential equations using the perturbation method. The linear and quadratic frequency response functions are obtained from the first- and second-order perturbation equations using the spectral method. Furthermore, under the assumption of the exponential covariance function of log hydraulic conductivity, the analytical solutions of both the spectrum and the variance of the hydraulic head produced from the linear system are derived. The results show that the variance derived herein is less than that of Gelhar (1977). The reason is that the log transmissivity is linearized in Gelhars work. In addition, the analytical solutions of both the spectrum and the variance of the hydraulic head produced from the quadratic system are derived as well. It is found that the correlation scale and the trend in mean of log hydraulic conductivity are important to the dimensionless variance ratio.     

Interannual variability of Gulf Stream warm-core rings in response to the North Atlantic Oscillation

Analysis of a quality-controlled database of Gulf Stream warm-core rings (WCRs) between 75° and 50°W during 1978–1999 demonstrates a significant correlation between WCR occurrences and variations in large-scale atmospheric forcing related to the state of the North Atlantic Oscillation (NAO). The mechanisms for linking the NAO with the rate of WCR occurrences are two-fold: (1) the influence of the NAO on Gulf Stream (GS) position, which could affect the interaction of the Gulf Stream with the New England Seamounts chain and thus allow for a higher/lower number of WCR occurrences; (2) the NAO-induced eddy kinetic energy (EKE) variability in the Gulf Stream region (GSR), which is indicative of the baroclinic instability processes necessary for WCR formation. Variability in GS movement is studied by analyzing annual mean positions of the Gulf Stream North Wall obtained from satellite-derived sea surface temperature (SST) frontal charts. Response of GSR EKE to fluctuations in the state of the NAO is examined with a numerical simulation of the North Atlantic basin from 1980–1999. The North Atlantic basin is simulated using a 1/6°-resolution eddy-resolving Regional Ocean Modeling System (ROMS) model that spins up with Southampton Oceanography Center (SOC) ocean-atmosphere atlas-derived atmospheric forcing fields. Model-derived EKE estimates are observed to be in good agreement with TOPEX/Poseidon altimeter-based EKE estimates as well as with results from other modeling studies for the North Atlantic basin. We suggest that lateral movement of the GS may not be the primary mechanism causing variation in the rate of WCR occurrences, because GS position is observed to respond at a lag of one year, whereas annual rates of WCR occurrences respond at 0-year lag to the NAO. Based on results from numerical simulations of the North Atlantic basin, adjustment to NAO-induced wind forcing is seen to impact the GSR EKE intensity and possibly the related baroclinic instability structure of the GS at 0-year lag. These results suggest that NAO-induced interannual variability in GSR EKE is the most likely mechanism affecting WCR occurrences. Numerical simulations show that high (low) phases in the state of the NAO exhibit higher (lower) EKE in the GSR, providing a greater (lesser) source of baroclinic instability to the GS front, possibly resulting in higher (lower) occurrences of WCRs.     

Intra-catchment variability in solute response an east devon example

Modifications to the land use of a small catchment, including clearance, cultivation and an increase in drainage density, are shown to have influenced a number of solute characteristics. Contrasts are drawn between three source areas; an unmodified subcatchment draining an area of mixed deciduous woodland, and a tile drain system and drainage ditch receiving water from cultivated areas. Not only are the modifications identified in terms of mean solute concentrations, but it is also shown that a change in land use alters the interaction between independent hydrological and meteorological variables, and solute response. A tentative explanation is provided for long term solute behaviour in relation to a number of biological geochemical and hydrological controls.     

Spatial patterns of simulated transpiration response to climate variability in a snow dominated mountain ecosystem

Transpiration is an important component of soil water storage and stream‐flow and is linked with ecosystem productivity, species distribution, and ecosystem health. In mountain environments, complex topography creates heterogeneity in key controls on transpiration as well as logistical challenges for collecting representative measurements. In these settings, ecosystem models can be used to account for variation in space and time of the dominant controls on transpiration and provide estimates of transpiration patterns and their sensitivity to climate variability and change. The Regional Hydro‐Ecological Simulation System (RHESSys) model was used to assess elevational differences in sensitivity of transpiration rates to the spatiotemporal variability of climate variables across the Upper Merced River watershed, Yosemite Valley, California, USA. At the basin scale, predicted annual transpiration was lowest in driest and wettest years, and greatest in moderate precipitation years (R² = 0·32 and 0·29, based on polynomial regression of maximum snow depth and annual precipitation, respectively). At finer spatial scales, responsiveness of transpiration rates to climate differed along an elevational gradient. Low elevations (1200–1800 m) showed little interannual variation in transpiration due to topographically controlled high soil moistures along the river corridor. Annual conifer stand transpiration at intermediate elevations (1800–2150 m) responded more strongly to precipitation, resulting in a unimodal relationship between transpiration and precipitation where highest transpiration occurred during moderate precipitation levels, regardless of annual air temperatures. Higher elevations (2150–2600 m) maintained this trend, but air temperature sensitivities were greater. At these elevations, snowfall provides enough moisture for growth, and increased temperatures influenced transpiration. Transpiration at the highest elevations (2600–4000 m) showed strong sensitivity to air temperature, little sensitivity to precipitation. Model results suggest elevational differences in vegetation water use and sensitivity to climate were significant and will likely play a key role in controlling responses and vulnerability of Sierra Nevada ecosystems to climate change. Copyright © 2008 John Wiley & Sons, Ltd.     

Photosynthetic pigment changes and adaptations in biofilms in response to flow intermittency

Among the environmental factors affecting benthic algae and cyanobacteria in streams, the one often producing the largest effects is flow intermittency. This study aimed to characterize the responses of algal assemblages to flow intermittency in a Mediterranean intermittent stream during the drying, non-flow (112 days), and rewetting phases. Algae growing in the epilithic, epipsammic and hyporheic streambed compartments were analyzed for pigment composition, and for the existence of structural changes in cells. Chlorophyll-a concentrations decreased between 60 to 90 % during the non-flow phase, indicating low resistance of algal assemblages to desiccation. In contrast, fast recoveries of Chlorophyll-a when flow resumed indicated high resilience. Pigment composition revealed that the epilithic algal assemblage was considerably different than the epipsammic and hyporheic ones. These differences were mainly attributed to the physical conditions prevailing on each streambed compartment that allowed the growth of different algal assemblages. During the non-flow phase, the synthesis of protective carotenoids (i.e. echinenone and scytonemin) and the occurrence of cell resistance structures (i.e. enlarged membrane thickness and resistant spores) enhanced resistance of the epilithic biofilm. The resistance observed in the epilithic biofilm might also be related to the tightly adhered growth-form of algae on this substratum. Main results suggest that algal assemblages in the epilithic compartment, which were the most exposed to desiccation, were structurally and functionally better adapted to flow interruption than those colonizing other streambed compartments, and that this compartment plays a crucial role in maintaining ecosystem functions under varying flow periods.     

SPM response to tide and river flow in the hyper-turbid Ems River

In this paper, we analyse the behaviour of fine sediments in the hyper-turbid Lower Ems River, with focus on the river’s upper reaches, a stretch of about 25 km up-estuary of Terborg. Our analysis is based on long records of suspended particulate matter (SPM) from optical backscatter (OBS) measurements close to the bed at seven stations along the river, records of salinity and water level measurements at these stations, acoustic measurements on the vertical mud structure just up-estuary of Terborg and oxygen profiles in the lower 3 m of the water column close to Leerort and Terborg. Further, we use cross-sectionally averaged velocities computed with a calibrated numerical model. Distinction is made between four timescales, i.e. the semi-diurnal tidal timescale, the spring–neap tidal timescale, a timescale around an isolated peak in river flow (i.e. about 3 weeks) and a seasonal timescale. The data suggest that a pool of fluid/soft mud is present in these upper reaches, from up-estuary of Papenburg to a bit down-estuary of Terborg. Between Terborg and Gandersum, SPM values drop rapidly but remain high at a few gram per litre. The pool of fluid/soft mud is entrained/mobilized at the onset of flood, yielding SPM values of many tens gram per litre. This suspension is transported up-estuary with the flood. Around high water slack, part of the suspension settles, being remixed during ebb, while migrating down-estuary, but likely not much further than Terborg. Around low water slack, a large fraction of the sediment settles, reforming the pool of fluid mud. The rapid entrainment from the fluid mud layer after low water slack is only possible when the peak flood velocity exceeds a critical value of around 1 m/s, i.e. when the stratified water column seems to become internally supercritical. If the peak flood velocity does not reach this critical value, f.i. during neap tide, fluid mud is not entrained up to the OBS sensors. Thus, it is not classical tidal asymmetry, but the peak flood velocity itself which governs the hyper-turbid state in the Lower Ems River. The crucial role of river flow and river floods is in reducing these peak flood velocities. During elongated periods of high river flow, in e.g. wintertime, SPM concentrations reduce, and the soft mud deposits consolidate and possibly become locally armoured as well by sand washed in from the river. We have no observations that sediments are washed out of the hyper-turbid zone. Down-estuary of Terborg, where SPM values do not reach hyper-turbid conditions, the SPM dynamics are governed by classical tidal asymmetry and estuarine circulation. Hence, nowhere in the river, sediments are flushed from the upper reaches of the river into the Ems-Dollard estuary during high river flow events. However, exchange of sediment between river and estuary should occur because of tide-induced dispersion.     

Evaluating Duero's basin (Spain) response to the NAO phases: spatial and seasonal variability

In this study, we analysed the influence of the North Atlantic Oscillation (NAO) on the hydrological response of the Duero River (central Spain) and its tributaries. We identified the positive and negative phases of the winter NAO for the period 1961–2006 and related precipitation and river discharge anomalies. Significant differences in precipitation and river discharge were found between the positive and negative NAO phases with negative anomalies (dry conditions) during positive NAO periods, and positive anomalies (wet conditions) during negative NAO periods. Marked differences were also found in the precipitation and river discharge response times to the NAO episodes, with an immediate response for precipitation but a lag and more sustained response for river discharges. Differing patterns were also identified in the response of river discharge to the NAO throughout the Duero basin. The physical characteristics of watersheds (including area, altitude and permeability explained most of the differences in the timing and magnitude of anomalies in the river discharge in response to the NAO. The findings highlight the great variability in the hydrological response of rivers to the NAO episodes. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of stiffness variability on the response of isolated structures

Results are presented of an investigation, the objective of which was to determine the relationship between the stiffness variability of the bearings of an isolation system and the response variability of the structure. The system is modeled as a rigid, rectangular structure that is free to translate and rotate. The isolation system consists of N isolation bearings arranged in a rectangular pattern, each with a stiffness k_i that is an independent, normally distributed, random variable. Response spectrum analysis is used to obtain the analytical solution for the structure response. Approximate closed‐form expressions are obtained for the variance of the centreline displacement, rotation, corner displacement and base shear, that are in terms of the variability of the isolator stiffness, aspect ratio of the structure, and the number and layout of isolation bearings. Results show that the standard deviation of the centreline displacement and base shear decrease with increasing number of isolation bearings, and are independent of the aspect ratio and layout of isolators, and in all cases are less than 1/4 the standard deviation of the isolator stiffness. The standard deviation of the corner displacement is a function of all of the system parameters, and is bounded below by the standard deviation of the centreline displacement and above by the standard deviation of a bar aligned perpendicular to the direction of ground motion with m isolation bearings distributed along the length. The approximate expressions are shown to be in good agreement with the results of Monte Carlo simulations. The results should be of use to designers of isolated structures and manufacturers of isolation systems, in assessing the significance of stiffness variability on the response of the isolated structure. Copyright © 2000 John Wiley & Sons, Ltd.     

Spatial variability of the poroelastic response of the heterogeneous medium

In this work, a stochastic methodology is applied to analyze the variability of the poroelastic response of the heterogeneous medium at the field scale. To solve the problem analytically, we restrict our attention to the one-dimensional models, where fluid flow as well as deformation occurs in one direction only under a constant applied stress. Assuming statistic homogeneity, the closed-form solutions that describe the variability of fluid pressure head, and a solid's strain and displacement are developed using a spectral approach based on Fourier–Stieltjes representations for the perturbed quantities. The influence of the correlation length of the log hydraulic conductivity on these results is investigated. It is found that the variances of the solid's strain and displacement increase with the correlation length of the log hydraulic conductivity, while the correlation length of the log hydraulic conductivity plays the role in reducing the variability of the specific discharge.     

Overland flow generation in chaparral ecosystems: temporal and spatial variability

Fire is an important and natural process in the lifecycle of chaparral systems, removing old growth and recycling nutrients. Recent catastrophic wildfires in southern California chaparral have heightened concerns about increased runoff and nutrient export. The goal of this study was to improve understanding of how overland flow is generated in unburned and post‐fire chaparral watersheds. Samples of overland flow were collected from burned and unburned watersheds after rainfall events and multiple regression analysis was used to examine the influence of individual storm characteristics and system moisture on overland flow volume. The results indicate that variation in overland flow generation in the unburned watershed is best explained by storm size, while overland flow in the burned watershed was positively related to storm size and time between storms. These findings suggest that the burned system had decreased infiltration rates and increased soil water repellency. In contrast, there is a statistically significant negative relationship between overland flow 1 year after a fire against different system and precipitation factors revealed a negative correlation with drying period and a positive relationship with rainfall intensity, a combination that suggests reduced repellency. Copyright © 2009 John Wiley & Sons, Ltd.     

ENSO and the natural variability in the flow of tropical rivers

This paper examines the relationship between the annual discharges of the Amazon, Congo, Paran á, and Nile rivers and the sea surface temperature (SST) anomalies of the eastern and central equatorial Pacific Ocean, an index of El Niño-Southern Oscillation (ENSO). Since river systems are comprehensive integrators of rainfall over large areas, accurate characterization of the flow regimes in major rivers will increase our understanding of large-scale global atmospheric dynamics. Results of this study reveal that the annual discharges of two large equatorial tropical rivers, the Amazon and the Congo, are weakly and negatively correlated with the equatorial Pacific SST anomalies with 10% of the variance in annual discharge explained by ENSO. Two smaller subtropical rivers, the Nile and the Paraná, show a correlation that is stronger by about a factor of 2. The Nile discharge is negatively correlated with the SST anomaly, whereas the Paraná river discharge shows a positive relation. The tendency for reduced rainfall/discharge over large tropical convection zones in the ENSO warm phase is attributed to global scale subsidence associated with major upwelling in the eastern Pacific Ocean.     

Sediment particle selectivity and its response to overland flow hydraulics within grass strips

Particle selectivity plays an important role in clarifying sediment transport processes in vegetative filter strips (VFS). 10-m long grass strips at slopes of 5° and 15° were subjected to a series of sediment-laden inflows experiments with different particle sizes to investigate the sediment transport and its response to overland flow hydraulics. The inflow sediments came from local soil, river-bed sand, and mixed, with median particle size d₅₀ of 39.9, 207.9 and 77.4 μm, respectively. Three independent repeated experiments were carried for each treatment. The results show that when the sediment trapping lasted for a certain length of time, the re-entrainment of some small-sized particles was greater than the deposition; that is, net loss occurred, which was not erosion of the original soil. Net loss of particles is mainly determined by the particle diameter. The coarser the inflow sediment particles and/or the steeper the slope, the coarser the particles can be net lost. Deposited sediment causes the VFS bed surface to become smooth and hydraulic resistance decrease exponentially. Unit stream power P is more suitable than shear stress τ of overland flow to be used to describe the process of sediment particle transport in VFS. The relationship between P and d₅₀ of outflow sediment is very consistent with the form of power function with a constant term. These results are helpful to understand the physical process of sediment transport on vegetation hillslopes.     

Precipitation variability and its relationship to hydrologic variability in Alberta

By applying wavelet‐based empirical orthogonal function (WEOF) analysis to gridded precipitation (P) and empirical orthogonal function (EOF) analysis to gridded air temperature (T), potential evapotranspiration (PET), net precipitation (P‐PET) and runoff (Q), this paper examines the spatial, temporal and frequency patterns of Alberta's climate variability. It was found that only WEOF‐based precipitation patterns, possibly modulated by El Nino Southern Oscillation (ENSO) and Pacific Decadal Oscillation(PDO), delineated Alberta into four major regions which geographically represent northern Alberta Boreal forests, southern Alberta grasslands and Aspen Parklands and the Rocky Mountains and Foothills. The leading mode of wavelet‐based precipitation variability WPC1 showed that between 1900 and 2000, a wet climate dominated northern Alberta with significant 4–8, 11 and 25‐year periodic cycles, while the second mode WPC2 showed that between 1960 and 2000, southern Alberta grasslands were characterized by decreasing precipitation, dominated by 11‐year cycles, and the last two modes WPC3 and WPC4 were characterized by 4–7 and 25‐year cycles and both delineated regions where moisture from the Pacific Ocean penetrated the Rocky Mountains, accounted for much of the sub‐alpine climate. These results show that WEOF is superior to EOF in delineating Alberta precipitation variability to sub‐regions that more closely agree with its eco‐climate regions. Further, it was found that while WPC2 could not explain runoff variations in southern Alberta, WPC1, WPC3 and WPC4 accounted for runoff variability in their respective sub‐regions. Copyright © 2009 John Wiley & Sons, Ltd.