Streamflow Response Research Articles

Hydrological impacts of climate and land‐use changes in a mountain watershed: uncertainty estimation based on model comparison

AbstractThis paper evaluates the response of streamflow and other components of the water balance to changes in climate and land use in a Pyrenean watershed. It further provides a measure of uncertainty in water resources forecasts by comparing the performance of two hydrological models: the Soil and Water Assessment Tool (SWAT) and Regional Hydro‐ecological Simulation System (RHESSys). Regional climate model outputs for the 2021–2050 time frame and hypothetical (but plausible) land‐use scenarios considering re‐vegetation and wildfire processes were used as inputs to the models. Results indicate an overall decrease in river flows (up to 30%, depending on the combination of scenarios) when the scenarios are considered, except for the post‐fire vegetation scenario, in which streamflows are simulated to increase (between 2% and 10%). However, the magnitude of these projections varies between the two models used, as SWAT tends to produce larger hydrological changes under climate change scenarios and RHESSys shows more sensitivity to changes in land cover. The final prediction will therefore depend largely on the combination of the land‐use and climate scenarios and on the model used. Copyright © 2014 John Wiley & Sons, Ltd.

Seeing the climate through the trees: observing climate and forestry impacts on streamflow using a 60‐year record

AbstractPaired watershed experiments involving the removal or manipulation of forest cover in one of the watersheds have been conducted for more than a century to quantify the impact of forestry operations on streamflow. Because climate variability is expected to be large, forestry treatment effects would be undetectable without the treatment–control comparison. New understanding of climate variability provides an opportunity to examine whether climate variability interacts with forestry treatments, in a predictable manner. Here, we use data from the H. J. Andrews Experimental Forest, Oregon, USA, to examine the impact of the El Niño‐Southern Oscillation on streamflow linked to forest harvesting. Our results show that the contrast between El Niño and La Niña events is so large that, whatever the state of the treated watershed in terms of regrowth of the forest canopy, extreme climatic variability related to El Niño‐Southern Oscillation remains the more dominant driver of streamflow response at this location. Improvements in forecasting interannual variation in climate might be used to minimize the impact of forestry treatments on streamflow by avoiding initial operations in La Niña years. Copyright © 2014 John Wiley & Sons, Ltd.

Hydrological Impacts of Climate Change Simulated by HIMS Models in the Luanhe River Basin, North China

This paper applied a HIMS (hydroinformatic modeling system) model to simulate streamflow in the Luanhe River Basin. This model was compared with SIMHYD and XAJ models for eight sub-basins of the Luanhe River. The results showed HIMS model performed better than SIMHYD and XAJ models, in these areas. We then investigated the streamflow response to climate changes in the different sub-basins. Twenty hypothetical climate change scenarios (perturbed temperatures and precipitation) were used to test the sensitivity of HIMS model simulated annual and mean monthly streamflow. Our results demonstrated that: (i) the annual streamflow was positively related to precipitation, and there was a negative relationship between streamflow and temperature for all the eight sub-basins; (ii) in all sub-basins, the relationship of annual streamflow change to precipitation change was highly non-linear, but the relationship of annual streamflow change with temperature change was approximately linear; (iii) the annual streamflow response to precipitation change was more sensitive when increasing than decreasing; (iv) the annual streamflow response to climate change was more sensitive in the Xingzhouhe River sub-basin, followed by the Wuliehe River sub-basin, and the Sahe River sub-basin was least sensitive; (iv) there were few differences in inner-streamflow response to climate change in the Laoniuhe, Yimatuhe, and Yixunhe Rivers. But for other rivers, when the temperature changed, larger streamflow differences happened in winter and summer; when the precipitation decreased or was unchanged, the larger differences happened in winter months, and when the precipitation increased, larger differences happened in winter and summer.

Statistical emulation of streamflow projections from a distributed hydrological model: Application to CMIP3 and CMIP5 climate projections for British Columbia, Canada

AbstractA recent hydrological impacts study in British Columbia, Canada, used an ensemble of 23 climate change simulations to assess potential future changes in streamflow. These Coupled Model Intercomparison Project Phase 3 (CMIP3) simulations were statistically downscaled and used to drive the Variable Infiltration Capacity (VIC) hydrology model over several watersheds. Due to computational restrictions, the 23 member VIC ensemble is a subset of the full 136 member CMIP3 archive. Extending the VIC ensemble to cover the full range of uncertainty represented by CMIP3, and incorporating the latest generation CMIP5 ensembles, poses a considerable computing challenge. Thus, we extend the VIC ensemble using a computationally efficient statistical emulation model, which approximates the combined output of the two‐step process of statistical downscaling and hydrologic modeling, trained with the 23 member VIC ensemble. Regularized multiple linear regression links projected changes in monthly temperature and precipitation with projected changes in monthly streamflow over the Fraser and Peace River watersheds. Following validation, the statistical emulator is forced with the full suite of CMIP3 and CMIP5 climate change projections. The 23 member VIC ensemble has a smaller spread than the full ensemble; however, both ensembles provide the same consensus estimate of monthly streamflow change. Qualitatively, CMIP5 shows a similar streamflow response as CMIP3 for snow‐dominated hydrologic regimes. However, by end‐century, the CMIP5 worst‐case RCP8.5 has a larger impact than CMIP3 A2. This work also underscores the advantage of using emulation to rapidly identify those future extreme projections that may merit further study using more computationally demanding process‐based methods.

Analyzing effects of climate change on streamflow in a glacier mountain catchment using an ARMA model

Impacts of the climate change on streamflow have been of great concern over recent decades. The streamflow from a glacier mountain in an arid or semiarid region, in particular, plays a vital role in ecological, social and economic developments of the region. In this study, a long-term climate and runoff data were collected from an upstream glacier mountain area in the Ürümqi River basin, China. An autoregressive-moving-average model (ARMA) was used to quantitatively analyze the influence of the air temperature and the precipitation on the streamflow that originated from mountain glaciers. We used the differences between two consecutive years in the monthly air temperature and the monthly precipitation over a period of 48 years (from January 1959 to December 2006) as our climate change indexes. Similarly, the corresponding differences of monthly runoff were treated in the same way to derive the streamflow response indicator time series. These three new time series (temperature, precipitation, and runoff) were then transformed by removing their seasonal trends to construct three stationary time series. The ARMA model was subsequently applied to analyze the transformed time series. Results of this analysis indicated that the runoff was related to the temperature and the precipitation at any given time and that the precipitation was more important than the temperature in controlling the streamflow. The runoff in the upstream of the Ürümqi River increased approximately 1 m3/s every 10 years due to the climate change.

Seasonal changes in runoff generation in a small forested mountain catchment

AbstractThis study aimed to investigate the seasonal variability of runoff generation processes, the sources of stream water, and the controls on the contribution of event water to streamflow for a small forested catchment in the Italian pre‐Alps. Hydrometric, isotopic, and electrical conductivity data collected between August 2012 and August 2013 revealed a marked seasonal variability in runoff responses. Noticeable differences in runoff coefficients and hydrological dynamics between summer and fall/spring rainfall events were related to antecedent moisture conditions and event size. Two‐component and three‐component hydrograph separation and end‐member mixing analysis showed an increase in event water contributions to streamflow with event size and average rainfall intensity. Event water fractions were larger during dry conditions in the summer, suggesting that stormflow generation in the summer consisted predominantly of direct channel precipitation and some saturated overland flow from the riparian zone. On the contrary, groundwater and hillslope soil water contributions dominated the streamflow response during wet conditions in fall. Seasonal differences were also noted between event water fractions computed based on isotopic and electrical conductivity data, likely because of the dilution effect during the wetter months. Copyright © 2014 John Wiley & Sons, Ltd.

Streamflow sensitivity to drought scenarios in catchments with different geology

AbstractStreamflow drought response depends to a large degree on groundwater recharge. To better predict and preempt streamflow droughts, the relationship between recharge deficit and streamflow response needs to be studied more systematically. We present a combined soil‐vegetation‐atmosphere transfer and conceptual groundwater model experiment that applies a novel set of recharge scenarios preceding drought events in humid‐temperate catchments with different dominant aquifer types. The recharge scenarios are based on the permutation of historical time series and on modified time series matching extreme (50 year) drought events. Karstic and fractured aquifers show short‐term sensitivity to drought with an event‐specific relationship between recharge and streamflow response. Porous and complex aquifers show long‐term sensitivity and a more catchment‐controlled propagation of drought. The sensitivity of drought deficit and recovery time correlates with the water age distribution in baseflow, which is trackable in the groundwater model and a characteristic that should be exploited to improve streamflow drought prediction.

Climate change poised to threaten hydrologic connectivity and endemic fishes in dryland streams.

Protecting hydrologic connectivity of freshwater ecosystems is fundamental to ensuring species persistence, ecosystem integrity, and human well-being. More frequent and severe droughts associated with climate change are poised to significantly alter flow intermittence patterns and hydrologic connectivity in dryland streams of the American Southwest, with deleterious effects on highly endangered fishes. By integrating local-scale hydrologic modeling with emerging approaches in landscape ecology, we quantify fine-resolution, watershed-scale changes in habitat size, spacing, and connectance under forecasted climate change in the Verde River Basin, United States. Model simulations project annual zero-flow day frequency to increase by 27% by midcentury, with differential seasonal consequences on continuity (temporal continuity at discrete locations) and connectivity (spatial continuity within the network). A 17% increase in the frequency of stream drying events is expected throughout the network with associated increases in the duration of these events. Flowing portions of the river network will diminish between 8% and 20% in spring and early summer and become increasingly isolated by more frequent and longer stretches of dry channel fragments, thus limiting the opportunity for native fishes to access spawning habitats and seasonally available refuges. Model predictions suggest that midcentury and late century climate will reduce network-wide hydrologic connectivity for native fishes by 6-9% over the course of a year and up to 12-18% during spring spawning months. Our work quantifies climate-induced shifts in stream drying and connectivity across a large river network and demonstrates their implications for the persistence of a globally endemic fish fauna.

Refining understanding of hydrological connectivity in a boreal catchment

AbstractLiterature has long documented how streamflow response in boreal hillslopes and catchments is influenced by storage capacities, thresholds, and landscape and runoff pathway heterogeneity. More recently, the influence these traits have on streamflow has been interpreted through the concept of hydrological connectivity. However, the nature of hydrological connectivity in boreal catchments has only begun to be described. The focus of hydrological connectivity studies at the catchment scale has been on discerning from which areas does runoff originate and when, but not necessarily discovering the source waters of this streamflow. This has been the realm of studies investigating residence time and runoff pathways. This article summarizes an investigation in a 155‐km2 catchment in Canada's Northwest Territories that applied both hydrometric and geochemical methods to measure streamflow response, storage state, connectivity and source waters. The goal of this research was to determine if a catchment scale metric of connectivity could be sensitive to changes in source waters and runoff pathways in a boreal landscape, so as to improve the predictive role of connectivity metrics. Nine runoff events from three water years were evaluated. There were distinct patterns of hydrological connectivity, which included, among others, a non‐linear relationship with the runoff ratio and hysteresis with streamflow. The results indicate that one metric of connectivity alone could not encapsulate connectivity extent and quality, identify water sources and be used successfully to predict runoff response. Multiple metrics were needed that further encapsulated the role of precipitation and hydrological processes on both structural and dynamic connectivity. © 2014 Her Majesty the Queen in Right of Canada. Hydrological Processes. © 2014 John Wiley & Sons, Ltd.

Paradox of Peak Flows in a Changing Climate

Paradox of Peak Flows in a Changing Climate

Is there a superior conceptual groundwater model structure for baseflow simulation?

AbstractPrevious work has shown that streamflow response during baseflow conditions is a function of storage, but also that this functional relationship varies among seasons and catchments. Traditionally, hydrological models incorporate conceptual groundwater models consisting of linear or non‐linear storage–outflow functions. Identification of the right model structure and model parameterization however is challenging. The aim of this paper is to systematically test different model structures in a set of catchments where different aquifer types govern baseflow generation processes. Nine different two‐parameter conceptual groundwater models are applied with multi‐objective calibration to transform two different groundwater recharge series derived from a soil‐atmosphere‐vegetation transfer model into baseflow separated from streamflow data. The relative performance differences of the model structures allow to systematically improve the understanding of baseflow generation processes and to identify most appropriate model structures for different aquifer types. We found more versatile and more aquifer‐specific optimal model structures and elucidate the role of interflow, flow paths, recharge regimes and partially contributing storages. Aquifer‐specific recommendations of storage models were found for fractured and karstic aquifers, whereas large storage capacities blur the identification of superior model structures for complex and porous aquifers. A model performance matrix is presented, which highlights the joint effects of different recharge inputs, calibration criteria, model structures and aquifer types. The matrix is a guidance to improve groundwater model structures towards their representation of the dominant baseflow generation processes of specific aquifer types. Copyright © 2014 John Wiley & Sons, Ltd.

Identifying dominant controls on hydrologic parameter transfer from gauged to ungauged catchments – A comparative hydrology approach

Summary Daily streamflow information is critical for solving various hydrologic problems, though observations of continuous streamflow for model calibration are available at only a small fraction of the world’s rivers. One approach to estimate daily streamflow at an ungauged location is to transfer rainfall–runoff model parameters calibrated at a gauged (donor) catchment to an ungauged (receiver) catchment of interest. Central to this approach is the selection of a hydrologically similar donor. No single metric or set of metrics of hydrologic similarity have been demonstrated to consistently select a suitable donor catchment. We design an experiment to diagnose the dominant controls on successful hydrologic model parameter transfer. We calibrate a lumped rainfall–runoff model to 83 stream gauges across the United States. All locations are USGS reference gauges with minimal human influence. Parameter sets from the calibrated models are then transferred to each of the other catchments and the performance of the transferred parameters is assessed. This transfer experiment is carried out both at the scale of the entire US and then for six geographic regions. We use classification and regression tree (CART) analysis to determine the relationship between catchment similarity and performance of transferred parameters. Similarity is defined using physical/climatic catchment characteristics, as well as streamflow response characteristics (signatures such as baseflow index and runoff ratio). Across the entire US, successful parameter transfer is governed by similarity in elevation and climate, and high similarity in streamflow signatures. Controls vary for different geographic regions though. Geology followed by drainage, topography and climate constitute the dominant similarity metrics in forested eastern mountains and plateaus, whereas agricultural land use relates most strongly with successful parameter transfer in the humid plains.

Increased evaporation following widespread tree mortality limits streamflow response

AbstractA North American epidemic of mountain pine beetle (MPB) has disturbed over 5 million ha of forest containing headwater catchments crucial to water resources. However, there are limited observations of MPB effects on partitioning of precipitation between vapor loss and streamflow, and to our knowledge these fluxes have not been observed simultaneously following disturbance. We combined eddy covariance vapor loss (V), catchment streamflow (Q), and stable isotope indicators of evaporation (E) to quantify hydrologic partitioning over 3 years in MPB‐impacted and control sites. Annual control V was conservative, varying only from 573 to 623 mm, while MPB site V varied more widely from 570 to 700 mm. During wet periods, MPB site V was greater than control V in spite of similar above‐canopy potential evapotranspiration (PET). During a wet year, annual MPB V was greater and annual Q was lower as compared to an average year, while in a dry year, essentially all water was partitioned to V. Ratios of 2H and 18O in stream and soil water showed no kinetic evaporation at the control site, while MPB isotope ratios fell below the local meteoric water line, indicating greater E and snowpack sublimation (Ss) counteracted reductions in transpiration (T) and sublimation of canopy‐intercepted snow (Sc). Increased E was possibly driven by reduced canopy shading of shortwave radiation, which averaged 21 W m−2 during summer under control forest as compared to 66 W m−2 under MPB forest. These results show that abiotic vapor losses may limit widely expected streamflow increases.

Streamflow response to regional climate model output in the mountainous watershed: a case study from the Swiss Alps

Regional climate model (RCM) outputs are often used in hydrological modeling, in particular for streamflow forecasting. The heterogeneity of the meteorological variables such as precipitation, temperature, wind speed and solar radiation often limits the ability of the hydrological model performance. This paper assessed the sensitivity of RCM outputs from the PRUDENCE project and their performance in reproducing the streamflow. The soil and water assessment tool was used to simulate the streamflow of the Rhone River watershed located in the southwestern part of Switzerland, with the climate variables obtained from four RCMs. We analyzed the difference in magnitude of precipitation, maximum and minimum air temperature, and wind speed with respect to the observed values from the meteorological stations. In addition, we also focused on the impact of the grid resolution on model performance, by analyzing grids with resolutions of 50 × 50 and 25 × 25 km2. The variability of the meteorological inputs from various RCMs is quite severe in the studied watershed. Among the four different RCMs, the Danish Meteorological Institute provided the best performance when simulating runoff. We found that temperature lapse rate is significantly important in the mountainous snow and glacier dominated watershed as compared to other variables like precipitation, and wind speed for hydrological performance. Therefore, emphasis should be given to minimum and maximum temperature in the bias correction studies for downscaling climatic data for impact modeling in the mountainous snow and glacier dominated complex watersheds.

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

AbstractRestored annual streamflow (Qr) 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 Qr 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 Qr 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 Qr, 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.

Impact of modellers' decisions on hydrological a priori predictions

Abstract. In practice, the catchment hydrologist is often confronted with the task of predicting discharge without having the needed records for calibration. Here, we report the discharge predictions of 10 modellers – using the model of their choice – for the man-made Chicken Creek catchment (6 ha, northeast Germany, Gerwin et al., 2009b) and we analyse how well they improved their prediction in three steps based on adding information prior to each following step. The modellers predicted the catchment's hydrological response in its initial phase without having access to the observed records. They used conceptually different physically based models and their modelling experience differed largely. Hence, they encountered two problems: (i) to simulate discharge for an ungauged catchment and (ii) using models that were developed for catchments, which are not in a state of landscape transformation. The prediction exercise was organized in three steps: (1) for the first prediction the modellers received a basic data set describing the catchment to a degree somewhat more complete than usually available for a priori predictions of ungauged catchments; they did not obtain information on stream flow, soil moisture, nor groundwater response and had therefore to guess the initial conditions; (2) before the second prediction they inspected the catchment on-site and discussed their first prediction attempt; (3) for their third prediction they were offered additional data by charging them pro forma with the costs for obtaining this additional information. Holländer et al. (2009) discussed the range of predictions obtained in step (1). Here, we detail the modeller's assumptions and decisions in accounting for the various processes. We document the prediction progress as well as the learning process resulting from the availability of added information. For the second and third steps, the progress in prediction quality is evaluated in relation to individual modelling experience and costs of added information. In this qualitative analysis of a statistically small number of predictions we learned (i) that soft information such as the modeller's system understanding is as important as the model itself (hard information), (ii) that the sequence of modelling steps matters (field visit, interactions between differently experienced experts, choice of model, selection of available data, and methods for parameter guessing), and (iii) that added process understanding can be as efficient as adding data for improving parameters needed to satisfy model requirements.

Streamflow and hydrogen ion interrelationships identified using data-based mechanistic modelling of high frequency observations through contiguous storms

With the aim of quantifying the purely hydrological control on fast water quality dynamics, a modelling approach was used to identify the structure (and dynamic response characteristics or DRCs) of the relationship between rainfall and hydrogen ion (H+) load, with reference to rainfall to streamflow response. Unlike most hydrochemistry studies, the method used makes no a priori assumptions about the complexity of the dynamics (e.g., number of flow-paths), but instead uses objective statistical methods to define these (together with uncertainty analysis). The robust models identified are based on continuous-time transfer functions and demonstrate high simulation efficiency with a constrained uncertainty allowing hydrological interpretation of dominant flow-paths and behaviour of H+ load in four upland headwaters. Identified models demonstrated that the short-term dynamics in H+ concentration were closely associated with the streamflow response, suggesting a dominant hydrological control. The second-order structure identified for the rainfall to streamflow response was also seen as the optimal model for rainfall to H+ load, even given the very dynamic concentration response, possibly indicating the same two flow-paths being responsible for both integrated responses.

Modeling the effect of glacier recession on streamflow response using a coupled glacio-hydrological model

Abstract. We describe an integrated spatially distributed hydrologic and glacier dynamic model, and use it to investigate the effect of glacier recession on streamflow variations for the upper Bow River basin, a tributary of the South Saskatchewan River, Alberta, Canada. Several recent studies have suggested that observed decreases in summer flows in the South Saskatchewan River are partly due to the retreat of glaciers in the river's headwaters. Modeling the effect of glacier changes on streamflow response in river basins such as the South Saskatchewan is complicated due to the inability of most existing physically based distributed hydrologic models to represent glacier dynamics. We compare predicted variations in glacier extent, snow water equivalent (SWE), and streamflow discharge with satellite estimates of glacier area and terminus position, observed glacier mass balance, observed streamflow and snow water-equivalent measurements, respectively over the period of 1980–2007. Observations of multiple hydroclimatic variables compare well with those simulated with the coupled hydrology-glacier model. Our results suggest that, on average, the glacier melt contribution to the Bow River flow upstream of Lake Louise is approximately 22% in summer. For warm and dry years, however, the glacier melt contribution can be as large as 47% in August, whereas for cold years, it can be as small as 15% and the timing of the glacier melt signature can be delayed by a month. The development of this modeling approach sets the stage for future predictions of the influence of warming climate on streamflow in partially glacierized watersheds.

Hydrologic response to valley‐scale structure in alpine headwaters

AbstractFew systematic studies of valley‐scale geomorphic drivers of streamflow regimes in complex alpine headwaters have compared response between catchments. As a result, little guidance is available for regional‐scale hydrological research and monitoring efforts that include assessments of ecosystem function. Physical parameters such as slope, elevation range, drainage area and bedrock geology are often used to stratify differences in streamflow response between sampling sites within an ecoregion. However, these metrics do not take into account geomorphic controls on streamflow specific to glaciated mountain headwaters. The coarse‐grained nature of depositional features in alpine catchments suggests that these landforms have little water storage capacity because hillslope runoff moves rapidly just beneath the rock mantle before emerging in fluvial networks. However, recent studies show that a range of depositional features, including talus slopes, protalus ramparts and ‘rock‐ice’ features may have more storage capacity than previously thought.To better evaluate potential differences in streamflow response among basins with extensive coarse depositional features and those without, we examined the relationships between streamflow discharge, stable isotopes, water temperature and the amplitude of the diurnal signal at five basin outlets. We also quantified the percentages of colluvial channel length measured along the stepped longitudinal profile. Colluvial channels, characterized by the presence of surficial, coarse‐grained depositional features, presented sediment‐rich, transport‐limited morphologies that appeared to have a cumulative effect on the timing and volume of flow downstream. Measurements taken from colluvial channels flowing through depositional landforms showed median recession constants (Kr) of 0.9–0.95, δ18O values of ≥−14.5 and summer diurnal amplitudes ≤0.8 as compared with more typical surface water recession constant values of 0.7, δ18O ≤ −13.5 and diurnal amplitudes >2.0. Our results demonstrated strong associations between the percentage of colluvial channel length within a catchment and moderated streamflow regimes, water temperatures, diurnal signals and depleted δ18O related to groundwater influx. Copyright © 2014 John Wiley & Sons, Ltd.

Do catchment characteristics explain differences in coherence and trends in hydroclimatic behaviour in an upland region?

To explore the link between catchment characteristics and sensitivity to environmental change, longer-term data from 21 contrasting mesoscale (67–690 km2) catchments in northern Scotland were analysed to examine the influence of hydroclimatic drivers on hydrological response. Cluster analysis was used to classify the catchments into four distinct groups, with topography and annual precipitation being the two most significant differentiating factors. Using 12 years of data from the 21 sites, inter-catchment coherence of intra- and inter-annual variation of precipitation, air temperatures and stream discharge were assessed both within and between clusters. Whilst catchments in clusters characterised by lower elevations exhibited coherence in hydroclimatic drivers and streamflow response, the higher altitude catchments did not. Annual trends were evident for widespread increases in temperature and also increasing values of precipitation in the higher elevation catchments. Seasonal trends indicate a lack of consistent change during winter, but most other seasons exhibit an increase in temperatures, with isolated increases in measures of precipitation and discharge particularly during the transitional seasons of spring and autumn. These trends ultimately could not be related with catchment typology where increasing upland area results in a breakdown of hydroclimatic coherence between sites.