The Upper Indus Basin (UIB) covers an area of more than 200,000 km2 and has an elevation range from below 1000 to over 8000 m above sea level. Its water resources underpin Pakistan's food security and energy supply. Vertical and horizontal variations in key climate variables govern the runoff contributions of the UIB's various elevation zones and subcatchments. Remote sensing climatic data products from NASA's Moderate Resolution Imaging Spectrometer (MODIS) instrument platform provide an opportunity to develop a spatial characterization of the climatology of remote and rugged regions such as the UIB. Specifically, snow-covered area (SCA) and land surface temperature (LST) have been shown to provide good analogues, respectively, for precipitation and air temperature. As such, SCA and LST quantify regional variations in mass and energy inputs to runoff generation processes. Although the 10-year (2000–2010) MODIS observational record is not adequate to evaluate long-term trends, it does provide a consistent depiction of annual cycles and a preliminary assessment of interannual variability. This study presents a summary of the period means and interannual variability found in remotely sensed SCA and LST products for the UIB. It then provides an update of locally observed recent climate trends for the 1962 to 2007 period. Nonparametric trend tests are applied both to the local observations and to remote sensing records to assess patterns in recent variability. The climatic noise (intense variability) of the past decade, however, renders conclusions on nascent trends in SCA and LST premature. Finally, runoff sensitivity to temperature change—spatially applied as summer (JJA) nighttime 0°C LST isotherm migration—is assessed for a range of potential scenarios. Results indicate that changes in mean summer (JJA) runoff could range from −30 to 35% or more, depending on whether recent locally observed changes continue or scenarios derived from current regional climate model (RCM) simulations unfold.
The Indus River is the lifeblood of Pakistan, principally in supporting irrigated agriculture and hence food security (Archer et al 2010). Hydropower also underpins energy security (Tarbela Dam supplies approximately 20% of the national electricity demand). In addressing the need for improved knowledge management to support enhanced water productivity, Karki et al (2011) note, “Disruption in the hydrological regime can have serious impacts on people's lives in the basin …” Thus, as Karki et al point out, key priorities are (1) improved understanding of hydroclimatic variability and potential and (2) expanded use of available remote sensing data sets.
Recent studies have attempted to address the macroscale issues of future water availability in major Asian river basins through the use of remote sensing. For example, Immerzeel et al (2010) assessed precipitation using Tropical Rainfall Measurement Mission (TRMM) (Huffman 2007). However, TRMM data, derived from infrared and passive microwave-based rainfall estimation algorithms, show deficiencies when measuring orographic precipitation (Dinku 2008). Field observations of high-elevation precipitation in the Upper Indus Basin (UIB) show significant orographic dependency (Hewitt 2011). For example, Forsythe et al (2011) note that for tributary catchments with seasonal snowmelt-type hydrological regimes with minimal glacial contribution, TRMM precipitation estimates account for only ∼40 to 60% of observed river discharge (expressed as runoff depth in mm). Thus, care is needed when extrapolating conclusions from these products to variations in glacial mass balance and future trends in water availability. The approach presented here treats remotely sensed data products as quantitative indicators of catchment conditions rather than as absolute measures.
Study area and data limitations
Locally observed data from the Pakistan Meteorological Department (PMD) and the Pakistan Water and Power Development Authority (WAPDA) are of crucial value in understanding regional climate. Ground-based data have been shown to be representative of large-scale hydrological variability (Archer and Fowler 2008) but do not adequately define altitudinal variations. Local observations can also be used to validate remotely sensed data, which provide good spatial coverage, albeit with limited historical record lengths. Useful spatial data products for catchments dominated by runoff from snowmelt and glacial melt are estimates of snow-covered area (SCA) and land surface temperature (LST) as available from the Moderate Resolution Imaging Spectrometer (MODIS) instrument. SCA and LST data products are validated here by (1) direct comparison to point observations or estimates derived from them and (2) correlations of variability in spatially aggregated parameters to variability in point-based analogues.
Large-scale hydrology of the UIB
The scale of vertical range and variation within the UIB has few analogues elsewhere in the world. Elevation ranges from below 1000 m, where the river emerges on the plains at the 2 major controlling reservoirs of Tarbela and Mangla (Figure 1), to several mountain peaks above 8000 m. This vertical range, combined with the magnitude of the annual temperature cycle and the accumulation and ablation of precipitation as snow and ice, defines the hydrology of the basin.
Runoff in the UIB is primarily composed of meltwater from ephemeral snow and perennial glacial masses, with a smaller contribution from direct winter or monsoon rainfall from foothill catchments (Archer 2003). The relative contributions from these 3 sources of runoff (direct rainfall, seasonal snow, and perennial ice) give rise to 3 distinct hydrological regimes (Archer 2003). While the contribution of precipitation inputs to snowpack mass balance is immediately obvious, the energy controls on meltwater generation processes are less intuitive. Daily minimum temperatures are a key element, as nighttime refreezing disrupts meltwater processes (Bengtsson 1982). A key metric for assessing thermal constraints to runoff generation is thus the fraction of the catchment where monthly mean nighttime temperatures remain above 0°C, defined here as the “continuous melt area” (CMA). The particular hypsometry of the UIB—the high percentage of surface area in the elevation range where the average daily minimum temperature (Tmin) 0°C isotherm fluctuates during the summer—means that this parameter is highly sensitive to incremental temperature anomalies. This is illustrated in Figure 2.
The contrasting controls of temperature and precipitation on the timing and magnitude of the annual hydrological cycle are shown here for the River Astore, a catchment dominated by melt of seasonal snow (Figure 3A, B), and for the River Hunza, where runoff is mainly derived from glaciers and perennial snow (Figure 3C, D). For the River Astore (Figure 3A), melt runoff is initiated while the catchment is still warming and recedes before it cools, once the available snowpack has been exhausted. For the River Hunza (Figure 3C), meltwater generation is correlated with CMA (Tmin > 0°C): flows increase in late spring and remain high until temperatures fall in early autumn. Comparing the interannual variations between the 2 catchments, the years of maximum and minimum flow volume do not match because the controlling parameters (precipitation and temperature) are largely independent at this timescale.
Note that CMA (Tmin > 0°C) in Figure 3 is determined by linear extrapolation of the estimated lapse rate to calculate the elevation of the 0°C isotherm. This elevation is then cross-referenced against the catchment hypsometry (ie the vertical distribution of catchment area) to identify the percentage of the catchment below the 0°C isotherm.
Previous studies (Archer 2003, Archer and Fowler 2008, Forsythe et al 2011) have shown significant linear correlations between gauged river flows and observed climate in the UIB, notably, concurrent temperature or preceding season precipitation. These relationships differ between the 3 identified regimes. Examples of such relationships are shown by Forsythe et al (2011) for the Astore catchment and for the high-level Hunza catchment (their Figure 7; repeated here as Figure 3A, D). Strong relationships are found between the relative anomalies of summer (June to August) discharge and winter (October to March) precipitation for the Astore catchment, and summer discharge and concurrent temperature at Gilgit shown as CMA (Tmin > 0°C) for the Hunza catchment. These relationships indicate the sensitivity of UIB hydrology to both present climate variability and developing trends.
Archer and Fowler (2008) used measurements from long-record meteorological stations to identify relationships between local climate observations and regional hydrological responses in the UIB. Forsythe et al (2011) utilized remotely sensed spatial data products to extend this methodology. In place of air temperature and precipitation, LST and SCA were used as indicators of inputs of energy and mass to meltwater generation processes contributing the bulk of runoff in the UIB. We first summarize the results of their study before detailing updated recent trends in observed climate variables, using parametric and nonparametric trend testing, and then test their correlation with trends in LST and SCA.
The MODIS SCA and LST data products
The data products are derived from radiometric measurements made by the MODIS instrument flown aboard NASA's Terra platform. The algorithms for both SCA and LST depend upon “clear sky conditions.” This means that for a given pixel, in order to determine the snow cover state or calculate its LST, it must be cloud-free. For this reason, data used here are from consecutive 8-day temporal aggregates assimilating the available observations during the period. Given the cloud climatology of the Upper Indus, most satellite scenes are partially cloud covered. In general, pixels over valley bottoms will be cloud-free more often than pixels covering mountain tops. Thus, the count of observations contributing to 8-day aggregates will most often be greater for low elevation zones than higher ones.
SCA, as 8-day aggregates of maximum extent—note that maximum snow extent indicates that a pixel observed as both snow covered and snow free by daily observations during the 8-day period is classified as snow covered—at a 500-m horizontal resolution (product MOD10A2 [Hall et al 2001]), is used as an analogue for precipitation observed at local meteorological stations. Day and nighttime LST, in the form of 8-day aggregates (mean of available observations) at a 1-km horizontal resolution (product MOD11A2 [Wan 2008]), are used as analogues of locally observed average daily maximum (Tmax) and minimum (Tmin) 2-m air temperature. The MODIS data products are available from early 2000 to the present.
Forsythe et al (2011) assessed the LST and SCA data products with respect to air temperature observations at Astore for 2000 to 2007. A time series of average LST was extracted for the pixel overlying Astore and compared to Astore average temperature (Tavg). SCA was assessed by estimating snow cover as the fraction of the study area where mean daily air temperature is below freezing (0°C, as in Jain 2008). This was then extrapolated in the same way as for CMA, using the hypsometric data for the Astore catchment. Both LST and SCA were found to track well with respective local estimates (Figure 4; Figure 3 in Forsythe et al 2011) in timing and interannual variability, especially for annual maxima. However, there was a clear and consistent divergence in annual minima for both data products. Forsythe et al (2011) explained this divergence using snow cover dynamics. During the winter, snow cover is nearly complete and the snow pack surface temperature becomes appreciably colder than the overlying air temperature through the accumulation of thermal inertia, also known as “cold content” (Bras 1990). During the spring and summer, as the snow pack warms and then melts and SCA diminishes, the LST and air temperature converge. The divergence between remotely sensed SCA and the simple estimate based on local air temperature appears during the ablation phase when the reduction in SCA (from snowmelt) lags behind the air temperature.
Analysis of relationships between spatial data and point observations
To demonstrate that SCA and LST are equivalent counterparts to locally observed precipitation and temperature, specific indicators were chosen within the SCA and LST data products as the closest spatial approximation of the local indicators used to characterize UIB hydrology (see “Large-scale hydrology of the UIB” above).
At the individual pixel-scale, SCA is an effectively discrete measure of the surface state—ie snow covered (1) or snow-free (0). Therefore, an analogue for catchment winter precipitation mass input is not immediately obvious. For the accumulated winter precipitation (ONDJFM) indicator, March SCA was selected, as it is usually close to the annual maxima and represents the catchment snow state at the transition from accumulation to ablation phase. Fortunately, while not a perfect analogue for air temperature, LST tracks reliably enough to be used as derived. The time-varying behavior of these indicators, aggregated as spatial means for each catchment, was tested against local observations from three long-running stations using “standardized anomalies.” The standard anomaly time series were calculated individually for each indicator, both local and spatial, by subtracting the period means from each observation and then dividing by the period standard deviation.
The results of the correlation analyses (Pearson's “r,” 2-tailed significance) are shown in Tables S1 and S2 (Supplemental data; http://dx.doi.org/10.1659/MRD-JOURNAL-D-11-00027.S1) for the catchments in Figure 5. Due to the steep vertical precipitation gradients reported in the UIB (Hewitt 2011), March SCA was tested against both accumulated winter precipitation depth and the seasonal total number of wet days. As precipitation is highly variable in space and time, correlation would be expected to be limited. Statistical significance is also challenging given the short common observation record (and annual time-step comparison). However, strong correlations at statistically significant levels are found in several cases. These strong correlations demonstrate the March SCA is a valid analogue for the accumulated winter precipitation (ONDJFM) observed at the local meteorological stations.
For the LST to air temperature correlations, the greater sample size (required length divided by the shorter monthly time-step) yields more statistically significant relationships. The higher correlation coefficients for daytime/Tmax comparisons over nighttime/Tmin counterparts are also expected, as the same pattern is found when testing correlation between local stations. The primary finding here is that correlations for both “Tmax versus Tday” and “Tmin versus Tnight” are frequently greater than 0.5, thus demonstrating the validity of LST—even spatially aggregated over large tributary catchments—as an indicator for local observations of air temperature.
Correlation analysis was used to further test the coherence of the spatial data products within each catchment. Physically, snow cover and surface temperature are linked not only by their concurrent states but also by cumulative antecedent precipitation and persistent temperature anomalies. Even so, the simple comparison of concurrent anomalies in the correlation analysis results shown in Table S3 (Supplemental data; http://dx.doi.org/10.1659/MRD-JOURNAL-D-11-00027.S1) demonstrates that UIB subcatchments are quite responsive as spatial units. An example of these relationships, for the Gilgit catchment, is graphically represented in Figure 6.
Given that variations in temperature will often be due to large-scale systems taking more than a single day to pass through the study area, strong correlations would be expected between daytime and nighttime LST. Strong negative correlations between both LST measures and SCA would also be expected, given that SCA will decrease during warm conditions and persist or increase during cooler weather. The coherence of these results further strengthens the case for the use of these spatial data products as hydroclimatic indicators.
Determination of hydrological regime types and runoff sensitivities from SCA and LST
A core concern for improving the resilience of infrastructure and resource management systems to cope with climatic perturbations is the accurate assessment of the likely climatic variability (Wilby and Dessai 2010). An initial step in this assessment is quantification of present variability. Figure 7A, B depicts the annual cycles (period mean) and observed interannual variability (estimated 10th [c10] and 90th [c90] percentiles) of SCA and LST-derived continuous melt area, respectively, from the MODIS record for the large gauged tributaries of the UIB. As shown in “Large-scale hydrology of the UIB” above, variability in runoff from these tributaries is known to track locally observed variation in mass (precipitation) and energy (continuous melt area) inputs. As shown in this section, SCA and LST can provide viable spatial analogues for these local observations.
The differences between catchments in the shape of their annual SCA cycles, as shown in Figure 7A, align with established hydrological regime types. The catchments with high annual minima, Hunza and Shigar, are also those with high glacial melt contributions and minimal hydrological responsiveness to year-to-year variation in observed precipitation. Those with low annual minima, Gilgit and Astore, are snowmelt-dominated catchments dependent on annual replenishment of the snowpack from winter precipitation. Snowmelt regime catchments tend to have very little contributing area above 6000 m above sea level (masl). The relative variability in annual SCA minima (values for c10 compared to c90) is striking, particularly for the snowmelt-type catchments, Astore and Gilgit, where c90 is roughly twice c10 for the minimum month.
For water resources management purposes, either extreme in annual minima, c10 or c90, could indicate resource stress conditions. An extremely low annual minimum could indicate low preceding snowpack input (snowfall drought) and thus low catchment yield. Alternatively, a high annual minimum could indicate inadequate melt energy (cold summer) and thus inadequate ablation. Clearly the annual SCA minimum alone is an incomplete indicator. At the opposite end of the melt season, however, a lower than normal annual SCA maximum (in February or March) could provide an early warning of summer water shortfalls, giving local stakeholders time to prepare and enact contingency plans.
As for SCA, differences in the shape of annual cycles of CMA, as shown in Figure 7B, are indicative of differing hydrological regime types. While for snowmelt-type catchments (Gilgit and Astore) annual SCA minima showed the greatest relative variability, for glacial melt catchments (Hunza and Shigar) the relative variability annual maximum CMA is the key parameter in interannual climate variability. These glacial catchments have a low mean annual CMA maximum but a large relative variability (ratio of c90 to c10). Snowmelt-type catchments, in contrast, regularly attain very high CMA maxima because little of their area remains melt-free during the high summer. Another interesting element from Figure 7B is the consistency of the very large ranges in CMA during the early stages of the ablation phase (May and June); this can affect the timing of peak flows for both regime types as well as the availability of meltwater for irrigation during crucial crop growth phases. Some of this variability could be a result of apparent seasonal low-elevation nighttime LST inversions found when studying the “parent data product” (MOD11A2). At present it is unclear whether this phenomenon is a real climatic feature of the UIB, perhaps linked to local-scale cloud cover variations, or an artifact resulting from a particular element in the algorithm.
Locally observed recent climate trends
The multidecadal life span of most water resources' infrastructure necessitates design not only for present climate variability but also for probable future conditions. A simple analysis of annual mean temperature change is inadequate to capture key hydrological processes, especially in meltwater-type regimes. Determining the direction and rate of year-on-year change throughout the annual cycle is essential in assessing probable future water resource availability (Rees and Collins 2005). Fowler and Archer (2006) analyzed the available records of 6 UIB meteorological stations for the period 1961 to 1999. While they identified a general upward trend of 0.007°C in annual mean temperatures, similar to the global trend over the same period, notable differences appeared when the contributions from individual seasons and the diurnal components (Tmax and Tmin) were separated. Specifically, they detected not only rising winter temperatures, primarily Tmax, but also cooling summer temperatures, especially summer Tmin. These asymmetrical trends in Tmax and Tmin thus yield an increase in diurnal temperature range (DTR) in all seasons. This contradicts general global patterns whereby faster increases in Tmin than Tmax yield decreasing DTR (Karl et al 1993; Easterling et al 1997). The present analysis updates the ground-based assessment to 2007 for 3 key stations.
Table 1 shows the estimated linear trends rates based on observations from the updated period (1962 to 2007) for the selected stations (Gilgit, Skardu, and Astore). Table 1 shows that the trend sign and magnitude, for either Tmin or Tmax, for each station generally agree with the others. The updated analysis supports the previous findings of substantial interseasonal asymmetry, as for most seasons and stations trends in Tmax and Tmin have opposite signs and contrasting rates. Tmax is observed to be rising quickly in winter (DJF), spring (MAM), and autumn (SON) but shows little if any change in summer (JJA). Indeed Tmin, by contrast, is decreasing rapidly in summer (JJA) with some decreases also in autumn (SON). In winter (DJF), weak increases or no changes in Tmin are observed. Thus, this analysis shows the increasing trend in DTR found by Fowler and Archer (2006) to be continuing throughout the annual cycle. DTR increases have also been reported for parts of India (Kumar et al 1994).
Linear rate (°C/y) statistical measures for trends in daily minimum temperature for 3 stations in the UIB.a)
Nonparametric trend assessment
While linear estimation of observed change is a familiar and easily conceptualized measure of nonstationary systems, other methods of trend assessment that are less sensitive to outliers and less dependent on record length for test performance are available. Two popular methods are the “rank Spearman” and “Mann-Kendall” tests. These tests have been previously applied by Shrestha et al (1999) to assess temperature trends across Nepal by physiographic region.
While the Spearman and Mann-Kendall methods are statistically robust, they are still subject to the underlying climatic behavior of the chosen periods. A well-publicized example would be to considered global temperatures across a range of time periods; for example, 1961 to 2010 or 1998 to 2010. The underlying physical drivers of temperature change do not vary between the 2 time windows, but owing to the large range of underlying interannual variability, the magnitude and even signs of the estimated trend can be affected by what is otherwise an arbitrary choice of analysis window. To illustrate this with specific reference to the UIB, Table 2 presents the Spearman (rho) and Mann-Kendall (tau) measures of correlations of the observed series to time. Both Spearman's rho and Mann-Kendall's tau have the same bounds: −1 to +1. Respective estimates of statistical significance (P value) are also given. Two analysis windows are presented for each parameter: (1) summer (JJA) Tmin and (2) winter (ONDJFM) wet days. The longer analysis window corresponds to the entire available record length. The shorter window corresponds to part of the record overlapping with the MODIS record.
Nonparametric measures of trend detection for local observations in summer (JJA) Tmin and winter (ONDJFM) wet days.
The signs, magnitudes, and P values from the Spearman and Mann-Kendall tests agree well within each parameter–window pairing. The magnitude, P value, and even sign can change dramatically between the analysis windows, however. The most striking example here is the opposite sign but equivalent magnitude for Tmin at Gilgit when switching from long window (1962 to 2007) to short window (2000 to 2007). The weaker (larger) P values for the shorter analysis windows are a function primarily of the shorter record length but also of the considerable interannual variability during the past decade. The UIB is similar to many regions globally in this regard. This variability should be borne in mind when considering the results of the Spearman and Mann-Kendall tests for the spatial data products as shown in Table 3.
Nonparametric measures of trend detection for spatial data: summer (JJA) CMA and March SCA.
It could be tempting to overinterpret the results presented in Table 3, to note consistent if weak negative correlation measures across the catchments for summer CMA and similarly consistent if weak positive values for March SCA. Objective consideration of the reported P values, however, inescapably points to the absence of strong, statistically significant trends. The latter finding is in fact the important element—ie whatever spatially distributed changes are under way, it will take time before trends clearly emerge from the ambient noise of interannual variability.
Results and discussion
In the study of the UIB, spatial climate data derived from satellite imagery and point measurements from local weather stations are complementary tools for understanding the drivers of a complex and dynamic hydrological system. Spatial data products provide key insights into patterns linked to steep topographic gradients. Long-record data series provided by a network of local weather stations permit quantitative assessment of regional patterns of variability and change. Previous studies have established that predictive relationships exist between meteorological point observations and discharge in the large gauged tributary basins of the UIB. Forsythe et al (2011) demonstrated that remotely sensed spatial data products (MODIS SCA and LST) can provide adequate analogues for these point observations. Used together, the 2 types of data may improve assessments of the likely hydrological impacts of projected climate change. These impacts may well vary between UIB subcatchments based on their hydrological regime type (eg perennial ice or seasonal snow). Examination of the variability only in the record of the available data products (early 2000 to present) indicates considerable potential hydrological variability in key seasons for water resources management.
In catchments such as the Hunza, where glacial melt provides a significant contribution to runoff, summer ablation season energy inputs—indexed by temperature—govern runoff volumes (Archer 2003). Specifically, these catchments appear to be sensitive to Tmin. Observations from ground-based measurements seem to indicate a strong decreasing trend in summer Tmin. Propagation of this trend from valley stations upward through the catchment is thus of particular concern, as it may greatly inhibit summer runoff generation. Projections from the Intergovernmental Panel on Climate Change (Solomon et al 2007), however, suggest increasing temperatures in the UIB in all seasons over the next century. Thus, it is not yet clear whether the recent local summer (JJA) Tmin observations simply show a transitory phenomenon as the climate seeks a new equilibrium state or if they are evidence of a regional divergence from the global pattern due to realignment of large-scale circulation or relative changes in components of the local energy balance.
Given this uncertainty, it is prudent to examine the potential consequences from a range of potential temperature changes. These consequences are quantified in Table 4. The isotherm elevation changes are based on a lapse rate of 0.0064°C/m. The potential changes in relative runoff volume for the Hunza catchment are calculated using the coefficient (1.3816) shown in Figure 3D. As can be seen in Table 4, even moderate temperature change (±2°C) could result in changes of ±20% in mean summer runoff. It is important to understand that these shifts in isotherm elevation would also have important impacts on glacial dynamics in the UIB by shifting the equilibrium line altitude (ELA) and thus the accumulation area ratio (AAR).
Potential impacts for the Hunza catchment from summer (JJA) minimum/nighttime temperature change.
In light of the discrepancy between recently locally observed trends and projected future change, it would be extremely valuable to corroborate the local ground-based measurements with long-record spatial data products. To this end, parallel work is currently under way on development of multidecadal spatial data products—equivalent to MODIS SCA and LST—using imagery from the Advanced Very High Resolution Radiometer (AVHRR) instrument flown on several generations of U.S. National Oceanic and Atmospheric Agency (NOAA) satellites (Chokmani et al 2005, Khlopenkov and Trischenko 2006, Li and Becker 1993, Pinheiro et al 2006). At present AVHRR sensor data covering the MODIS record have been acquired and processed using the initial versions of the algorithms adapted from the methodologies of the authors cited. The validation process will thus include comparison with the two equivalent MODIS data products as well as evaluation against local observations using the analysis approach presented here.
If successful, construction of a multidecadal spatial record of SCA and LST would allow both direct estimation of recent trends and assessment of relationships between these regional patterns and larger-scale atmospheric phenomena. Long-record LST analyses aim to assess the multidecade evolution of 0°C isotherm elevations and thus the distribution of solid (snowfall) versus liquid (rainfall) precipitation. Long-record SCA data will be assessed as depletion versus available discharge. Extending the record of spatial data backward in time would be especially useful because this would greatly increase the overlap with available river gauging data. Combinations of scenarios within the observed record that are judged likely to become more frequent in the future may then be used as analogues for future hydrological conditions.
Nathan Forsythe was supported by a Graduate Research Fellowship award from the U.S. National Science Foundation as well as by project-based support from the British Council and the School of Civil Engineering & Geosciences of Newcastle University. Hayley Fowler was supported by NERC Postdoctoral Fellowship award NE/D009588/1 (2006–2010). Funds for travel between the UK and Pakistan and for capacity building activities were provided by the British Council through its PMI2Connect and INSPIRE research exchange award programs. We wish to thank the Pakistan Meteorological Department and the Water and Power Development Authority for the provision of climate and flow data.