Influence of Historical Climate Patterns on Streamflow and Water Demand in Wales, UK

: Ensuring reliable drinking water supplies is anticipated to be a key future challenge facing water service providers due to fluctuations in rainfall patterns and water demand caused by climate change. This study investigates historical trends and relationships between precipitation, air temperature and streamflow in five catchments in Wales, before correlating these with actual total abstraction data provided by the water company, to give insight into the supply-demand balance. Changes in seasonal and annual averages, as well as extreme events, are assessed for a 34-year period (1982–2015) and a breakpoint analysis is performed to better understand how climate has already changed and what this might mean for the future of water supply. Results show a north-south divide in changes in extreme temperature and streamflow; a strong warming trend in autumn average temperatures across Wales (Sen’s slope range: 0.38–0.41, p <0.05), but little change in precipitation. Abstraction, as a proxy for overall water demand, is shown to be positively correlated to temperature (Spearman’s ρ value range: 0.094–0.403, p <0.01; Pearson’s r value range 0.073–0.369, p <0.01) in four of five catchments. Our study provides new insight into the relationship between abstraction volume and hydroclimatic factors and highlights the need for catchment-scale water resource planning that accounts for hydroclimatic variations over small spatial distances, as these nuances can be vital. and summer one day minimum flow volume. The Conwy shows an increase in winter minimum flow volume, along with an increase in the maximum summer flow, while the Clwyd displays an increase in annual maximum flow volume.


Introduction
Water service providers face a vast array of challenges and uncertainties when planning their future operations and services. Brown et al. [1] compiled a list of 94 priority research questions for the UK water sector, in which the impact of climate change on water quantity was ranked as the second most important question. Their work also highlighted the need to better understand the drivers of water demand, both domestic and commercial, in order to improve future demand forecasting. Previous work on the characterization of hydroclimatic trends in the UK suggests that precipitation and streamflows have become more seasonal, a pattern that is expected to continue [2][3][4]. A study by Christiersen et al. [5] showed, using UK Climate Predictions 2009 (UKCP09) data, that by the late 2020s, increases in winter precipitation levels are likely to be more prominent in northern and western parts of the UK, while decreases in summer flows will be seen more generally across the whole country. Similarly, a hydrological modelling study by Prudhomme et al. [6] showed that summer precipitation and streamflows will decrease across the UK by varying amounts for the period 2040-2069; whereas future winter precipitation and streamflows showed an upward trend, especially for Wales. Extreme precipitation events are also projected to become more seasonal in the UK, with longer duration and more intense rainfall events in winter becoming more common [3]. Mayes [7] suggested that these anticipated changes will not be uniform across the UK, and current rainfall gradients are likely to be accentuated, i.e., the south getting dryer in summer and the north getting wetter in winter.
Regional scale understanding of water resource provision in the UK is particularly important because water supply is under the control of individual water companies that serve separate regions of varying sizes, populations, and physical characteristics. For instance, in south-east England, which is already a water stressed area, studies on future hydroclimatic trends suggest that summer streamflow levels will continue longer into autumn, with overall summer flow levels declining also. Furthermore, winter streamflows will increase and continue longer into spring, leading to accentuated seasonality in terms of season longevity and flow volumes [8][9][10][11][12]. However, an increase in winter precipitation will do little to combat summer shortages if no further storage capacity is developed soon [13]. Borgomeo et al. [14] suggested that, due to the combined effect of climate change and significant predicted population growth, the London water supply zone urgently required both supply and demand-side interventions if the current standard of water provision is to continue. In Scotland, although winter precipitation is predicted to increase in the future, a lower percentage of it will fall as snow [15,16]. This will make catchments more responsive to winter precipitation and increase the pressure on water managers to deal with larger discharge events [17]. For summer precipitation and streamflows, Blenkinsop & Fowler [18] noted that Scotland has a limited amount of groundwater storage capacity, which heightens the drought risk from any reduction in non-winter precipitation. In Wales, studies suggest that winter and summer season characteristics, e.g. wet winters and dry summers, will be exacerbated, especially in winter [19][20][21].
In this study, we use Wales as a case study region, a country often viewed as abundant in water resources, receiving some of the highest average annual rainfall totals in the UK [22], but which in reality does have zones of water deficit [23]. Wales is also important due to its role as an exporter of water to major metropolitan areas in England. Dŵr Cymru Welsh Water (DCWW), the major water service provider for Wales, has over 20 bulk water trades, the largest of which supplies 360 million liters per year to Severn Trent Water for distribution around Birmingham [23,24].
Past studies on water resources in Wales have predominantly been conducted either as part of UK-wide research [2][3][4]19,[25][26][27], or with a focus on the combined England and Wales region [11,28]. When focusing specifically on the area covering Wales, Fowler & Wilby [19] projected a much larger magnitude of increase in winter flows from the 1960-1990 average to thirty-year averages centered on 2025, 2055 and 2085, than the corresponding decreases in summer precipitation. Dixon et al. [20] showed a significant upward trend across 56 Welsh and West Midlands catchments between 1962-2001 for winter high flow values, but no significant changes in the mean annual values. Conflictingly, Macdonald et al. [29], demonstrated that during the period 1973-2002, there was no significant change in the seasonality of rainfall across 30 catchments in Wales, which they proved to have a significant link to streamflows; however, they did show that the frequency of occurrence of extreme precipitation events in Wales had increased during the study period. These contrasting results highlight the need for careful consideration when selecting study period length and timeframe due to the potential impact on trend analysis results and projections.
Up until recently, supply-side measures to tackle water scarcity and to manage water resources have traditionally been the main path towards a reliable sustainable future water network. However, it has increasingly been recognized over the past decade that demand-side interventions should also play a role as an adaptation measure [30]. For this to be a viable option, further work is needed to understand the relationship between prevailing and antecedent weather conditions, and demand for water in the UK. Several studies have looked at the general interplay between the two in the UK [31][32][33][34], and abroad [35][36][37]. This interplay has been under research since at least as far back as the 1990s, with Herrington [32] stating that at the time, up to 40% of total consumption in summer can be due to garden watering, which is obviously highly affected by the prevailing weather conditions. Goodchild [38] used summer daily domestic water demand data (55% of UK piped water supply at the time) from 41 domestic properties and daily meteorological data to develop a demand prediction model. The model included ten weather variables to account for current and antecedent conditions; evapotranspiration, days since rain, and temperature were all important functions. This modelling work projected a 2.1% increase in average summer 7-day household demand by the 2020s. More recently, Parker & Wilby [39] reviewed domestic water demand in the UK and noted the lack of studies on weather and climate. It is also important to look not only at domestic demand, but also industrial, agricultural and non-revenue water use (e.g. leakage), as they all influence the long-term sustainability of water supply.
In this study, we look at the implications of past trends in hydroclimatic data on two of the problems identified by Brown et al. [1]: (1) impact of climate variations on water availability and (2) understanding the factors affecting water demand. The first problem has been addressed by assessing trends in seasonal and annual average climate and streamflow data as well extreme event frequency and magnitude, using Mann-Kendall trend analysis and breakpoint analysis. We have addressed the second problem by investigating the historical links between hydroclimatic factors and total water demand, using actual abstraction data provided by DCWW as a proxy for demand. To our knowledge, this is one of the first studies conducted independently of a water service provider to utilize actual abstraction data provided by a water service provider in this manner in the United Kingdom. The reliance on Wales for water supply in other regions, combined with a potentially inaccurate assumption of national water abundance in Wales, makes the region a crucial area of study in terms of water management and water supply availability. Our research has therefore been undertaken in order to provide information for future water resource planning and policy decisions, as well as future research. It is hoped that this will be achieved by providing evidence of the longterm trends and links between prevailing weather and flow conditions, as well as total demand for water in Wales.

Study Catchments
Our study focusses on the rivers Clwyd and Conwy in the north, the Dyfi in the west, and the Teifi and Tywi in the south of Wales ( Figure 1). We selected these catchments due to them being among the largest systems within the region, encompassing a range of Land Use/Land Cover (LULC), and exhibiting a variety of different catchment characteristics. Furthermore, the catchments are also mostly encompassed within DCWW's water supply zone (Figure 1), and all have multiple surface water abstraction locations licensed to DCWW for use for public water supply. Brief catchments descriptions have been provided below and key comparative details about each are given in Table 1. Figures for LULC in the catchments were calculated from 2012 CORINE Landcover data (Copernicus Land Monitoring Data), while the catchment size, longest stream length, elevation and slope values were calculated using the 5 m Digital Elevation Model (DEM) data provided by Ordnance Survey.
In terms of LULC, the Clwyd, Teifi, and Tywi are dominated by agricultural land use ( Figure 1); while the largest proportion of LULC in Conwy and Dyfi is scrubland (mainly in the form of moors and heathland in Conwy, and natural grasslands in Dyfi). Agricultural land in all study catchments is mostly pastureland, however, Clwyd does have a larger proportion of its total LULC as arable land (10.4%). Forested land in all catchments is mainly coniferous, except in Conwy where there is an equal proportion of coniferous and mixed forest (5.1% LULC each) and 3% broadleaf forest. The bedrock geology is predominantly mudstone, siltstone and sandstone in all five catchments, ranging from about 60% in Conwy and Clwyd, 73% in Tywi, and over 85% in Dyfi and Teifi. The catchments range in area covered from 1363 km 2 at the Tywi to less than half of that for the Conwy (564 km 2 ). The steepest catchment is the Dyfi, with an average slope of 14.2 degrees, being largely in the south of the Snowdonia region, while the Teifi in the south is the least steep catchment, with an average slope of 6.6 degrees.

Data
Hydrological and meteorological data for all five study catchments was obtained for the longest complete period possible for all datasets, that being the 34-year period from 1 October 1981 to 30 September 2015 (except flow data for the Dyfi, which was only available to 5 May 2014). Daily precipitation and air temperature data were obtained from the Centre for Ecology and Hydrology's (CEH) Climate, Hydrology and Ecology research Support System (CHESS) dataset. The mean value of all of the 1 km grids contained within each catchment was taken for each day for both precipitation and temperature, giving a daily average value for each variable for each catchment. Daily streamflow data for the catchments was obtained from CEH's National River Flow Archive dataset; gauging station locations are shown in Figure 1.
Daily water abstraction volumes were obtained from DCWW, which is the main water supply company for this region. However, these data were only available for a five-year period from 1 January 2012 to 30 September 2016. Therefore, comparison of the water abstraction data with the hydroclimatic data could only be completed for the overlapping dates between the datasets, 1 January 2012 to 30 September 2015. The daily water abstraction data was provided for 22 abstraction locations within the five study catchments and has been used as a proxy for overall demand across the network for all consumers, as well as leakage.

Trend Analysis
We selected the nonparametric Mann-Kendall (MK) test [40,41] to detect any consistent trends in the hydroclimatic data over time. The decision to use a nonparametric test was taken due to the nature of both climate and hydrology data generally being non-normally distributed and displaying seasonality, which goes against the assumption of constant distribution [42]. Moreover, hydrology data often displays auto-correlation, making it unsuitable for parametric testing [43]. The MK test maintains the sequential order of the data and calculates Kendall's tau (a measure of association between two samples) between each value and all values proceeding it, to test for a monotonic increase or decrease relationship in the data over time [44]. MK has also been successfully applied to similar hydroclimatic data in various locations globally [45][46][47][48] so was deemed suitable for use in this study. Details on the specific equations relating to the MK test can be found in Jaiswal et al. [49].
Long-term average trends in the datasets were examined by analyzing seasonal and annual averages for each catchment; where winter is December to February; spring is March to May; summer is June to August; autumn is September to November; and where, for example, the 1982 hydrological year runs from 1 October 1981 to 30 September 1982. For the abstraction data, we have taken the sum of abstractions in each catchment, with trend analysis being performed on this value, in order that the results are comparable to the hydroclimatic factors analyzed. We also examined trends in the frequency and magnitude of extreme events relating to hydroclimatic factors. This was achieved firstly by analysis of trends in maximum and minimum temperatures; maximum one day precipitation total and cumulative rainfall totals; and maximum and minimum average one day streamflow volumes. Secondly, we undertook an analysis of "events over threshold", in order to establish whether there are generally more or less extreme weather events per year/season over the study period. This was completed by taking the 5th and 95th percentile values of a whole dataset (seasonally and annually) and analyzing the number of times in each year and season that that value is surpassed (95th percentile) or not reached (5th percentile) for temperature, and surpassing the 95th percentile only for precipitation. The MK test was applied with Sen's slope estimator [50] in order to estimate the size and direction of trends in the data; Hamed and Rao's method of auto-correlation correction was also applied in order to remove any apparent trend which the data exhibits with itself over time [43].
In order to detect any sudden changes in the hydroclimatic data, we carried out a breakpoint analysis on the seasonal/annual average and extreme event data. This analysis was completed using two well-established methods, the Standard Normal Homogeneity Test (SNHT) method [51,52] and the Pettitt method [53]. These two methods were used to ensure the widest possible detection of breaks in the data; with the SNHT being more accurate at detecting breaks at the start and end of time series, and the Pettitt method being more reliable in the middle [54,55]. If annual values of the variable being tested are identically distributed and independent, both methods accept the null hypothesis. The alternative hypothesis (that there is a change point in the series) will be accepted however, if a shift in the value of the mean has occurred [54]. Both methods are also location specific, so will identify the year at which the mean changed within the time series. For further information on the equations that drive the two methods, we suggest Hänsel et al. [54] and Jaiswal et al. [49].

Correlation Analysis of Hydroclimatic and Abstraction Data
In order to investigate the relationship between the three hydroclimatic factors and water abstraction we calculated both Pearson's & Spearman's rank correlation coefficients. The analysis was performed for each catchment with daily and monthly average temperature, precipitation and streamflow all being compared separately to daily and monthly total abstraction volume. We calculated Pearson's coefficient values (r) using Equation (1): where ℎ refers to the daily average hydroclimatic variables studied, refers to daily total abstraction volume, and ℎ & represent the mean of the entire respective datasets. The nonparametric Spearman's correlation coefficient (ρ) was calculated using Equation (2), a modified version of Pearson's correlation coefficient which calculates correlation between ranks, as opposed to raw data: where (ℎ ) refers to the rank of the daily average hydroclimatic variables studied, ( ) refers to the rank of the daily total abstraction volume, and (ℎ ) & ( ) represent the mean rank of the entire respective datasets.

Results
Our results indicate a north-south divide when looking at extreme temperature and streamflow changes. Additionally, a strong warming trend in average autumn temperatures across Wales is observed, however little change in precipitation is seen. Abstraction, as a proxy for overall water demand, is strongly positively correlated to temperature in all catchments and negatively with streamflow and precipitation in all catchments except the Dyfi. Below we provide a detailed breakdown and explanation of these results.

Precipitation
We found no significant trends during MK analysis or breakpoint analysis of annual and seasonal precipitation averages in any of the five catchments for the study period; but Sen's slope indicator values do suggest a slight decrease in average spring precipitation across all catchments. In addition, only two significant trends were present in terms of extreme precipitation events; the first is a decrease in total cumulative spring precipitation in the Clwyd catchment throughout the period (Tau = −0.184, p = 0.040; Sen's Slope = −1.196). The second is a decrease in the number of extreme wet autumn days in the Teifi catchment (Tau = −0.593, p = 0.003; Sen's Slope = −0.138). Although both of these trends are mirrored in the other catchments, those trends are not statistically significant, highlighting the need to study individual catchments to account for varying characteristics. Additionally, the Pettitt breakpoint analysis showed a marked increase in winter cumulative precipitation in the Conwy after 1989 (p = 0.040), the mean for the post 1989 period increasing by 13.4% compared to the pre-1989 period.

Temperature
Unlike precipitation, air temperature data does display significant trends in Wales over the course of the study period. In all catchments, we observe a warming trend in average autumn temperature ( Figure 2); the rate of increase is marginally higher in the two north Wales catchments compared to the two in the south of the country ( Table 2). These findings correlate with both the Pettitt and SNHT breakpoint analysis in all catchments, which shows a step increase in autumn temperatures in 1994, with the percentage change in pre-and post-1994 mean temperatures being larger in the north (Table 3). We found no significant trends in the MK analysis of annual, winter, spring or summer datasets, however when looking at the breakpoint analysis, further changes are seen. In the two most northerly catchments, winter temperatures show a break and increase under SNHT analysis in 1987 (p = 0.042 & 0.048 for the Clwyd and Conwy respectively). All catchments also show a breakpoint increase in average spring temperatures in 1987 and annual average temperatures in 1988, both under SNHT analysis.  When looking at extreme temperature events, the three most southerly catchments, the Dyfi, Teifi and Tywi, display both a decrease in the number of hottest days annually (days above the whole dataset 95th percentile; Figure 3), as well as decrease in the number of coldest summer days (days below the summer dataset 5th percentile; Figure 4) across the period; suggesting a narrowing of temperature ranges, especially in summer. The Conwy catchment in the north also displays the latter trend of fewer of the coldest summer days (Table 4). This narrowing of temperature ranges has however, not been abrupt enough to cause a breakpoint in the data, with no significant changes seen.

Streamflow
We observe only one statistically significant trend in seasonal and annual average streamflow data, this being in annual flows in the Teifi; here a significant increase in average annual flow ( Figure  5) is detected over the study period (Tau = 0.240; Sen's slope = 0.204; p = 0.044). The Teifi also shows the only observed breakpoint in the average flow data, that being a step increase in winter flow in 2013 (p = 0.019) the mean for the post 2013 period increasing by 63.0% compared to the pre-2013 period. We found various trends in terms of extreme events in streamflow in Wales (Table 5); once again the three most southerly catchments show similar trends, in this instance, an increase in volume of the annual and summer one day minimum flow volume. The Conwy shows an increase in winter minimum flow volume, along with an increase in the maximum summer flow, while the Clwyd displays an increase in annual maximum flow volume.

Correlation Analysis of Hydroclimatic and Abstraction Data
When looking at the daily actual abstraction data (1 January 2012 to 30 September 2015) under both correlation coefficient tests, all catchments display a positive relationship between air temperature and volume of water abstracted (Table 6), except for the Dyfi, which shows the reverse trend. All catchments, except the Dyfi, also show a statistically significant negative correlation between precipitation and volume of water abstracted, and consequently a negative relationship between streamflow and abstraction volume ( Table 6). In order to investigate the situation for Wales as a whole, we have also included in Table 6 the correlation results for total daily and monthly abstraction volume in all five catchments with average daily/monthly temperature, total daily/monthly rainfall volume, and total daily/monthly streamflow volume, for all five catchments combined. The correlations seen in this dataset are consistent with those seen in most individual catchments. In addition, all catchments also display a positive relationship between streamflow and precipitation, and a negative one between streamflow and air temperature. Table 6. Pearson's and Spearman's rank correlation coefficient results for daily and monthly abstraction data with hydroclimatic factors (2012-2015). Statistics relating to a combined analysis of total daily/monthly abstraction in all catchments with average daily/monthly temperature, total streamflow volume and total precipitation volume across all catchments is also provide. Values for r and ρ underlined are significant at p <0.05, and those in bold are significant at p <0.01.

River
Hydroclimatic  Figure 6 shows the aforementioned combined data at a daily time-step from all catchments which has been normalized to each factor's maximum and minimum dataset value. Clear relationships between the hydroclimatic factors and actual abstraction data are observed. When looking at temperature, a positive correlation can be seen in the scatter plot, while the time series plot clearly also shows a large amount of consistency between the average temperature and total water abstraction volume. Arguably the clearest trend seen in the scatter plots is total streamflow volume and total actual abstraction volume, showing a negative correlation. This can also be seen in the time series plot that clearly shows drops in abstraction volume at times of largest flow volumes. Generally, the weakest correlations shown in Table 6 relate to precipitation and actual abstraction, the same is true for the graphs presented in Figure 6. Although peaks in precipitation do tend to coincide with a drop in abstraction, there is much more noise in the precipitation dataset, which makes these relationships more difficult to pinpoint. The correlations and relationships seen in the daily data are made more evident when looking at the combined normalized monthly data, as shown in Figure 7. Clear negative trends can be seen between total abstraction volume with both total precipitation and total streamflow volume, when looking at scatter and time series line graphs. Clear peaks in abstraction can be seen at time of lowest streamflow and precipitation levels, and vice versa. Temperature and total abstraction are shown to also be broadly in-line when looking at the time series plot, with a clear positive trend shown when looking at the scatter data.

Hydroclimatic Trends
The results of the trend, breakpoint, and correlation analyses, show a selection of spatially varying changes across the five catchments over the 34-year study period. Wales has a maritime climate which is strongly influenced by the North Atlantic Oscillation (NAO), a major cause of atmospheric circulation variability in the north Atlantic region [20,56]. This is important to consider when analyzing changes that have occurred across all catchments. NAO displays year-to-year and longer-term variability, therefore breakpoints in temperature, precipitation, and streamflow data that occur across all five catchments and likely to be rooted in changes in NAO. However, when explaining changes that affect a subset of the catchments, it is likely that catchment characteristics in terms of topography, LULC and surrounding landscapes have had an impact by interacting with these large scale meteorological processes [57]. The breakpoint analysis is a prime example of this interplay, a combination of NAO, climate change and other non-climate/weather related changes at a catchment and national level are likely to have impacted on the climatology breaks shown. No single factor accounts for all of breaks shown across the catchments. While a change in NAO is likely to have been a key driver, other factors such as LULC change and climate variations may well have had an impact on spreading the breakpoints out over the range of approximately 7 years shown (1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994).
When looking at both temperature and streamflows, a north-south divide in results can often be seen. Taking the extreme events results for streamflows as an example, it can be seen that in the two northerly catchments, the Clwyd and Conwy, the largest annual one-day discharge events have become larger over the study period. Meanwhile the two most southerly catchments, the Teifi and Tywi, have seen the smallest one-day flow values becoming larger, in particular when looking annually and at summer flows. These results are consistent with the findings of Dixon et al. [20] who showed that in the period 1977-2001, summer minimum flow volume values had increased significantly in south Wales catchments, but not the north-mainly due to the more mountainous terrain of north Wales (Figure 1). Furthermore, increases in the lowest summer flow values were also seen in Osborn and Hulme [58], and were linked to increased light summer rainfall over the period. The Dyfi catchment in mid-Wales shows a combination of these two patterns with maximum and minimum one day flows increasing in the summer and annually, showing a shift to generally wetter summers over the period, again corresponding to the reasoning and results put forward in Dixon et al. [20] for mid-Wales. These results are also consistent with the mean catchment slope data shown earlier (Table 2) with the Dyfi and Conwy in particular being the steepest catchments, and arguably the flashiest, of the five studied. This factor may well contribute to the increase in summer maximum flows observed over the study period.
It is interesting to note that the aforementioned changes in streamflows occur despite a lack of corresponding statistically significant change in annual or seasonal average precipitation, an observation also made by Macdonald et al. [29]. This mismatch of significant trends between precipitation volumes and streamflows is surprising, as the two factors are correlated significantly in all catchments, suggesting that flows in the catchments are highly sensitive to very small changes in precipitation. Similarly, the only significant trend seen in extreme precipitation events, a decreasing trend in cumulative spring rainfall in the Clwyd, does not relate to a corresponding change in average or extreme spring streamflows during the study period. Although not analyzed in this study, it is also possible that changing LULC during the study period could have had an effect on streamflows, independent of precipitation. Changes in agricultural land use in particular could have a large impact in the study catchments, affecting processes such as water infiltration and runoff [59]. Between the early 1980s and the late 2000s there was an increase in the amount of agricultural land classified as permanent grassland and rough grazing, despite overall agricultural land area remaining relatively stable [60]. From the late 2000s to present, the total amount of agricultural land has grown to its highest level since World War II, with the vast majority of this being for grazing [60]. Furthermore, forest cover has been continually increased throughout the study period, with greater planting taking place at the start of the study period and the rate steadily declining throughout [61]. Forested land also impacts on the connectivity of river flows and precipitation, with processes such as evapotranspiration and interception storage diverting water from, or delayed water reaching rivers [59]. This could once again explain some of the mismatch between significant precipitation and streamflow trends.
Changes in seasonal and annual averages, as well as extreme events are most observed in temperature, for the factors studied. Once again, a north-south divide can be seen in parts of these results, with the three southernmost catchments displaying a decrease in both the annual number of days that are hotter than the dataset 95 th percentile, and the number of summer days that are colder than the dataset 5 th percentile. These two changes suggest a narrowing of the temperature ranges in these catchments, with less extreme hot days annually and less extreme cold days in the summer.
These findings are in line with wider UK research such as Dessai & Hulme [62] who showed, via the Central England Temperature record, that annual average summer temperatures increased between 1960-2007 when compared to a 1961-1990 mean. Nesbitt et al. [63] also presented a trend for summer days becoming warmer over the period 1954 to 2012 for south-east and south-central UK, as did Luterbacher et al. [64] for the period 1977 to 2003. Annual NAO index has been largely negative on average throughout the study period, especially so in the latter two thirds (Figure 8). Negative NAO has been linked to colder maximum temperatures [56], while climate change has been causing increasing summer temperatures on average. These two factors combined could be contributing to the narrowing temperature range observed, with NAO bringing down the maximum, and climate change bringing up the minimum. Furthermore, in all catchments, an increase in average autumn temperatures across the study period has been observed, with this being slightly more pronounced in the two north Wales catchments. All of the observed trends also fit with general UK observations of a greater degree of warming in the south of the country than the north, and the exacerbation of the temperature gradient between them. Furthermore, the north-south divide also fits when considering the surrounding geography of the catchments, the mountainous Snowdonia region lies just to the south of the two most northerly catchments, heavily influencing the climate here and causing local variation in the weather that is brought in from the Atlantic [20]. Our study period has been highlighted in yellow; updated and adapted from Osborn [65].

Weather and Demand for Water
When looking at the correlations between daily and monthly total abstraction data and hydroclimatic factors, it is clear that both average temperature and precipitation volume play a crucial role in influencing the total amount of water abstracted; volumes abstracted increasing in higher temperatures, and decreasing on wetter days. Given that the data takes account for total abstraction volume, these relationships are not surprising, as it is likely that on hotter days there will be both more domestic demand (more showers, car washing, garden watering etc.), as well as higher demand from agriculture, in particular for water for livestock in the catchments studied. The opposite is true for wetter days, which reduce the agricultural demand from public water supplies and reduce domestic demand, especially for water use on external areas of a property. These overall relationships do fit with other work that has sought to quantify the relationship between climate variables and domestic demand, both in the UK and abroad. Slavíková et al. [66] showed that in the Czech Republic, air temperature relative to the season average, accounted for the most variability in residential water consumption. Similarly in the UK, Goodchild [38] showed that domestic water demand in 41 houses in an Essex case study correlated with an R 2 value of 0.44 with days when maximum temperatures were over 25 °C. Additionally, total daily sunshine hours were shown to have a stronger correlation with demand, having an R 2 value of 0.53 [38]. Work by Xenochristou et al. [67] also found that sunshine hours and air temperature were the most influential weather variables on domestic demand, along with humidity.
However, using total abstraction volumes as a proxy for demand does present some challenges. For example, factors such as the usage of water internally, within water treatment works, as well as network water leakage, may mask or alter observed trends in the dataset. In addition, it is difficult to apportion the total abstraction volume to different user groups, such as domestic, industrial and agricultural on a daily basis. Furthermore, different user groups may have different relationships with weather conditions, for example Massoud et al. [68] showed for California's Central Valley, that while agricultural demand increased in dryer years, precipitation volume had little impact on urban water demand. Additionally, we recognize that the comparison period between the hydroclimatic factors and abstraction data is relatively short in this study, at four years. This does not give the opportunity to investigate relationships such as the effect of prolonged drought on water use, or other longer terms patterns. Nevertheless, the overall abstraction data does give a baseline relationship to work with when considering the impact of future climate change on total water demand.

Study Implications
Our study has shown that the climate of Wales has changed since the early 1980s, and that this will have contributed to both the supply and demand of water in the region. On the basis of the results found, it is clear that these changes have been more keenly felt in terms of average air temperatures than precipitation volumes; this applies to both annual and seasonal averages as well as extreme events. When looking at the impact that these climatic changes have had on streamflows it can be seen that seasonal and annual average flows have remained largely unchanged, instead it is extreme flow events that have been more greatly affected. These changes occur with a north-south divide, with the largest annual flow events getting larger in the north, and the south becoming less dry in the summer when looking at the very lowest flows. These changes could have important impacts if they are continued into the future, with implications not only for water supply, but also in terms of water resource management, to prevent flooding and other related natural disasters. Furthermore, large industrial users as well as applications such as hydroelectric developments could be impacted by changes in river regime. In particular the viability of some small scale hydroelectric installations could be called in to question in some areas due to changing streamflow characteristics and flow duration curves, further emphasizing the importance of a solid understanding of the relationship between, and emerging trends in, hydroclimatic factors.
Climate change induced alterations in the future timing, quantity and quality of water available for supply, as well as policy relating to adaptation and management methods to cope with these predicted changes, needs to be further researched for Wales. In particular, hydrological modelling studies comparing current baseline streamflow (such as that presented in this paper) to future streamflow under various climate change scenarios, could prove to be particularly useful. This work is crucial to better inform future water supply-demand dynamic assessments, water resource management, and adaption planning. This is especially true if projected increases in the reoccurrence, duration, and intensity of extreme events under future climate change, as suggested in other research [69][70][71][72][73], are correct. This suggested research would also go some way to addressing some of the aforementioned priority research questions laid out in Brown et al. [1]; but must however also keep in mind the cost and practicality of adaptation measures, in order to ensure both a continued unbroken water supply service, and affordable water for all. More broadly, this research has shown the need to research and understand historical trends and future projections of hydroclimatic factors at local, catchment levels. This is clear when looking at the presented differences in observed trends seen over small spatial distances, due to changing land characteristics, and it is these nuances that are vital to incorporate in to future planning for any industry that relies on surface water abstraction.

Conclusions
This research has highlighted the potential for water scarcity problems even in a relatively waterrich region such as Wales. For example, with observed trends such as warmer average autumn temperatures providing for potentially greater water use in the season, the pressures on summer water supply could in the future extend further into the autumn. Although potentially increased demand could be countered by a trend of the largest discharge events becoming larger in the north of Wales, and summers becoming less dry in the south of Wales, any increase in flow is of little use if the capacity to store this additional water is not sufficient to make use of it.
Finally, we suggest that further research should focus on how future climate change will affect the relationship between weather factors, streamflow, and water demand, both in Wales and globally. For example, research concerning "trigger temperatures" for significant increases in water use, or the effect of long-term higher than averages temperatures on water demand, would aid understanding of the finer detail of the dynamic between hydroclimatic factors and total water abstracted. We also hope that this paper will set a frame onto which future climate change research focusing on surface waters, and the future provision of water services can be built; being one of the first steps in securing the long-term sustainability of water supply services in the region and further afield. Funding: This study was conducted as part of the Dŵr Uisce project, which aims at improving the long-term sustainability of water supply, treatment and end-use in Ireland and Wales. The project has been supported by the European Regional Development Fund (ERDF) through the Interreg Ireland-Wales Co-operation Programme 2014-2020. The work has been used as a case study for Dŵr Uisce Work Package 7 "Climate Change", and forms part of a study assessing the impact of future climate change on water supply in Wales.