From the 1990s climate change has decreased cool season catchment precipitation reducing river heights in Australia’s southern Murray-Darling Basin

The Murray-Darling Basin (MDB) is Australia’s major agricultural region. The southern MDB receives most of its annual catchment runoff during the cool season (April–September). Focusing on the Murrumbidgee River measurements at Wagga Wagga and further downstream at Hay, cool season river heights are available year to year. The 27-year period April–September Hay and Wagga Wagga river heights exhibit decreases between 1965 and 1991 and 1992–2018 not matched by declining April-September catchment rainfall. However, permutation tests of means and variances of late autumn (April–May) dam catchment precipitation and net inflows, produced p-values indicating a highly significant decline since the early 1990s. Consequently, dry catchments in late autumn, even with average cool season rainfall, have reduced dam inflows and decreased river heights downstream from Wagga Wagga, before water extraction for irrigation. It is concluded that lower April–September mean river heights at Wagga Wagga and decreased river height variability at Hay, since the mid-1990s, are due to combined lower April–May catchment precipitation and increased mean temperatures. Machine learning attribute detection revealed the southern MDB drivers as the southern annular mode (SAM), inter-decadal Pacific oscillation (IPO), Indian Ocean dipole (IOD) and global sea-surface temperature (GlobalSST). Continued catchment drying and warming will drastically reduce future water availability.


Drying trends in the Murray-Darling Basin
Aside from the 2017-2019 drought, the MDB previously experienced major multi-year droughts in 1997-2009 (Millennium Drought), 1938-1946 (World War II Drought) and 1895-1903. There are downward trends in rainfall since the 1990s in southern parts of eastern Australia, and in the higher rainfall areas of the MDB catchments that generally are located in the highlands of the Great Dividing Range 5 . These trends are concentrated in the cooler half of the year, from April to September 3 . A major influence on this drying trend has been the strengthening and extension of the subtropical high pressure ridge during winter, which shifts rainbearing weather systems southward 5,6 . Specifically, the lack of negative-phase Indian Ocean Dipole (IOD) events has been identified as a contributing factor to the drying and resulting droughts in southeast Australia since the 1990s 7 , where the influence of the IOD is greatest in June-October 8 . Other important climate drivers of cool season precipitation are the atmospheric-based southern annular mode (SAM), which in its positive phase is also a major influence on the drying trend, by its association with rain-bearing systems shifting southward, away from southern Australia 3,4 , and the basin-wide Pacific Ocean phenomenon, the inter-decadal Pacific oscillation (IPO). These and other possible climate drivers are assessed in the Results Section, using Machine Learning (ML) techniques to attribute identification.

Effect of regulation on river heights
Prior to the 1990s, the main regulation of the Murrumbidgee River upstream of Hay occurred over two periods during which large storages were built. These are the Burrinjuck Dam in 1910-1927, which started filling during the construction phase and Blowering Dam 1956-1961. Before the Burrinjuck Dam was built, annual river flows at Hay followed a pattern of variability and quantity similar to annual river flows measured at Wagga Wagga (Fig. 1). After the Burrinjuck Dam was built (1928), the water flows diverged in quantity although their annual variability patterns coincided 9 . After the Blowering Dam construction was completed in the 1960s, the relationship between river flows at Wagga Wagga and Hay were affected. There were similar patterns of variability but there was a substantial difference in the quantity of water reaching Hay from Wagga Wagga, compared with that before river regulation 9,10 . It was concluded that there was a 56% reduction in the lower Murrumbidgee River at Hay 11 , which occurred mostly after 1957 when there was a change point between low regulated flow and highly regulated flow. Even after the 1957 change point, river flows were seasonal, driven by reliable winter and spring rainfall, together with headwater ice and snow melt, and much of the flow reached the floodplain below Hay 11 . However, there has been a sharp drop to zero in JJAS Murrumbidgee River flooding events at Hay since 1991 (Fig. 2a), with just one minor flood in 1995 at Wagga Wagga and one major flood in 2016 (Fig. 2b), in addition to much reduced dam inflows since the mid-1990s (Fig. 2c,d).
In Fig. 2a there is a data overlap of approximately 13 years at Hay town gauge and Hay Private Irrigation District (PID) gauge for which the coefficient of determination (R 2 = 0.9235), obtained from linear regression of the weekly river height values (N = 668), explains 92% of the variance in the river heights between the two locations. The data quality is high as it takes into account missing values, outliers, continuity between the two gauge locations, and measurement practice 12 . The flooding level of 6.7 m at Hay is defined as the level above which water starts to spill over the riverbanks at Hay 12 and the minor flood level at the river gauge at Wagga Wagga is 7.3 m.
It is hypothesized here that from the early 1990s the main influence on decreased downstream river water availability is due to the reduction in dam net inflows from decreased precipitation. It is well-known that rainfall in southeast Australia, where the Murrumbidgee River catchment is located, usually is highest in winter and  www.nature.com/scientificreports/ has a minimum through summer into early autumn 5 . Hence the so-called 'winter cropping season' starts with planting in autumn and harvesting in spring-summer.

Monthly-seasonal maximum river flood heights at Hay 1874-2019
Most river height exceedances above 6.7 m at Hay in JJAS occur in late winter (July, August) and September (Table 1) after sufficient catchment wetting. Spring (SON) floods are the most common (Table 2) as river and  dam inflows are a result of both rainfall and snowmelt. However, since 1991 there has been just one April to  September flood at both Hay and Wagga Wagga resulting from Burrinjuck Dam overspill which was in September  2016 (Table 2). Extreme SON precipitation totals are associated with La Niña over eastern Australia 13 . In terms of river flood events at Hay and Wagga Wagga, and net inflows at the catchment dams of Burrinjuck and Blowering in southeast Australia, this relationship broke down, as shown by the presence only of neutral phases of ENSO in SON after the 1970s (Table 2) and in agreement with eastern Australian rainfall in general 13 . However, there was a strong relationship with ENSO for the 2010 flood in summer (DJF). In autumn (MAM) there are very few flood events (Table 3) because the catchment typically dries out sufficiently through summer to reduce river flow.

Precipitation time series in the Murrumbidgee river catchment area
Representing major dam storages, Burrinjuck Dam and Blowering Dam Figs. 3(a,b) each receives a mean JJAS precipitation of about 300 mm, whereas Tumut receives approximately 350 mm (Fig. 3c). There are fewer percentiles near the extremes from the 1990s, with no years above the 95th percentile apart from 2016 and a reduction in years below the 5th percentile. Statistical significance is discussed in the next section. Similarly, for April-May precipitation there are no years above the 90th percentile from the late 1990s for the three locations ( Fig. 3d-f).
All three locations, particularly Burrinjuck Dam and Blowering Dam, are exposed to topographically influenced precipitation from mid-latitude wet season (JJAS) cold frontal weather systems.
Monthly to seasonal precipitation typically exhibits large ranges due to its episodic nature, whereas seasonal means show little change. Therefore, an important part of this study is to analyse trends in both the means and variances of catchment rainfall and inflows, in addition to JJAS river heights at both Hay and Wagga Wagga. The focus is on the period after the 1960s, since there has been no new regulatory infrastructure upstream of Hay. The main focus is from the early to mid-1990s corresponding to the accelerated global warming signal 14 .
Considering changes in regulation of the Murrumbidgee River, box and whisker plots and permutation testing are used to analyse changes in both JJAS and April-May inflows at Burrinjuck/Blowering Dams, and precipitation at key rainfall catchment locations including Burrinjuck Dam, Blowering Dam and Tumut, in addition to river heights at both Hay and Wagga Wagga. The 27-year period chosen for the time series analysis is justified in the Methods Section. For the April-May and JJAS Burrinjuck Dam, Blowering Dam and Tumut combined precipitation, ML techniques are used to identify possible key climate drivers of Murrumbidgee River catchment precipitation. Wavelets then are used to provide some further understanding of possible climate drivers affecting precipitation and net inflows at those locations, and how their influence has changed over time.

Results
For each of the Burrinjuck Dam, Blowering Dam and Tumut time series there has been no statistically significant trend after bootstrapping and box plotting the precipitation into 27-year intervals, for 1965-1991 and 1992-2018. While their JJAS means appear to have decreased slightly (not shown), no locations have observed statistically significant changes (Table 4). For the Tumut JJAS data, precipitation was significantly more variable prior to 1965 (p-values = 0.045 and 0.0984; Table 4). However, for all three locations, their mean late autumn (April-May) precipitation shows a decline in both the mean and variance in the 27-year boxplots between 1965 and 1991 and 1992-2018 ( Fig. 4a-c). Statistically, there is a highly significant decrease in both their means (p-values < 0.01; Table 4), and their variances (p-values ~ 0.05; Table 4). It is noteworthy that while neither the mean nor variance of Tumut April-May precipitation between 1884 and 1910 and 1992-2018 is significant, their variances are similarly very low ( Fig. 4c-right panel), owing to the very low rainfall that occurred during the Federation Drought (1895-1903).
There is an apparent decrease in the mean and variance of JJAS Murrumbidgee River heights at Hay between the bootstrapped intervals 1938-1964 to 1965-1991 highlighted by their low p-values (Table 4; figure not shown), which corresponded to a change from low to highly regulated river flow 11 . However, the apparent mean and variance decrease is not statistically significant, whereas the variance between the 27-year bootstrapped intervals 1965-1991 to 1992-2018 is highly significant (p-value ~ 0.005; Table 4). The mean river height in this period for Hay (p = 0.236; Table 4) is not statistically significant most likely due to snow melt starting to reach Hay in September. Both Hay and Wagga Wagga river heights are significantly lower in both mean and variance from 1911 to 1937 to 1992-2018 owing to a high rainfall period 1911-1937 compared to the period 1992-2018 (p-values = 0.0068, 0.001, respectively for Hay; 0.0002, 0.013, respectively for Wagga Wagga). For Wagga Wagga the decrease in both the mean and variance from 1965 to 1991 to 1992-2018 is highly significant (p-values ~ 0.0068 and 0.001, respectively; Table 4). Importantly, the significant mean decrease at Wagga Wagga suggests that the cause is climate related with river water used for irrigation not a major contributor since almost all irrigation occurs downstream between Wagga Wagga and Hay. There is no significant change in October-March precipitation mean or variance from 1965 to 1991 to 1992-2018 for the three catchment rainfall locations (p-values not shown).
The mean JJAS precipitation wavelets for Burrinjuck Dam, Blowering Dam, Tumut and the mean JJAS weekly Murrumbidgee River heights at Hay (Fig. 5a-d), respectively, exhibit a significant ENSO-like periodicity of approximately 2-7 years, which mostly weakens in the three precipitation locations from the 1990s, and, in www.nature.com/scientificreports/ the case of the river heights at Hay, disappears after the 1990s. A possible reason for the weakening ENSO-like periodicity is that the nonlinear ENSO amplitude has weakened (less strong El Niños) in response to global warming 15,16 .
There is a clear periodicity of 2-7 years, suggestive of the ENSO time scale, at the three catchment precipitation locations (Fig. 5a-c). However, the ENSO signal disappears from the late 1990s for the downstream river heights at Wagga Wagga and Hay (Fig. 5d,e). This period corresponds to a lack of major precipitation-driven net inflows except for 2010-2012 and 2016.
Apart from ENSO there are many potential climate drivers that could influence precipitation in the Murrumbidgee catchment, and therefore river heights along the Murrumbidgee River. Potential attributes considered in this study include the IOD, global sea surface temperature anomalies (GlobalSSTA), global temperature anomalies (GlobalT), Niño3.4, the IPO, the SAM, the Southern Oscillation Index (SOI), and Tasman Sea sea surface temperature anomalies (TSSST). Additionally, two-way interaction terms between these predictors were considered as potential attributes, obtained through multiplication of two variables with each other (e.g., SAM*IPO). For April-May precipitation, the important attributes were assessed as GlobalT, SAM, SOI, GlobalT*SAM and GlobalT*IPO. For JJAS precipitation, the key attributes were IOD, SAM, SOI, IPO and GlobalT*GlobalSSTA.
There is a clear increasing trend in JJAS mean TMax and TMin at both Burrinjuck Dam and Cabramurra (located close to Blowering Dam), revealed by their 27-year interval boxplots ( Fig. 6a-d), and consistent with the known global warming signal since 1950 17 which has accelerated in the last 50 years 18,19 and particularly since the mid-1990s in Australia 13 . This is confirmed by the very high levels of significance (all p-values are < 0.05) for the difference in TMax and TMin means between the intervals 1965-1991 to 1992-2018 (Table 5).
Cabramurra (altitude 1482 m) shows a highly significant increase in TMin between the two bootstrapped intervals in contrast to a non-significant increase at Burrinjuck Dam (altitude 390 m). With such a large difference in altitude between the two locations, it is possible that other factors may be involved, including local differences in meteorological variables such as wind speed/wind direction and/or cloud cover.

Discussion
The JJAS river heights at Hay (Fig. 2a) have clearly reduced variability over the 27-year period 1965-1991 compared with 1992-2018 (p-value = 0.005, Table 4). Despite the apparent decrease in the mean and variance of JJAS Murrumbidgee River heights at Hay from the 1960s, as shown by the low p-values (0.309, 0.332), respectively, for the 27-year intervals 1938-1964 to 1965-1991 (Table 4), the decrease is not statistically significant. However, the decrease is consistent with the suggestion that a change point occurred from the late 1950s between unregulated to regulated flow at Hay 11 .
While the mean JJAS precipitation of the three catchment locations of Burrinjuck Dam, Blowering Dam and Tumut indicate a slight decrease in percentile extremes from the 1990s, with 2016 (due to September) only above the 95th percentile (Fig. 3), they exhibit no significant change in mean or variance based on the bootstrapped intervals JJAS 1964-1991 to 1992-2018. Consequently, the question arises of why the significant decrease in mean and variance of the JJAS Murrumbidgee River height at Wagga Wagga and in variance at Hay does not match a similar significant decrease in mean or variance of JJAS catchment rainfall. A rainfall decline in recent decades was found to be most pronounced in late autumn 20,21 and that without sufficient autumn rainfall to moisten catchments in southern Australia, follow-up rainfall in winter cannot be efficiently converted to runoff and catchment inflows 22 . There have been statistically significant decreases in April-May mean precipitation at the catchment locations of Blowering Dam, Burrinjuck Dam and Tumut from 1964 to 1991 to 1992-2018 and also for the mean inflows to the two Dams (Table 4). Furthermore, as a result of the Millennium Drought (1997-2009), modelling experiments indicate that, starting from very dry conditions, the run-off response to rainfall only will return to the normal pre-drought conditions after about 10-20 years of average rainfall 23 . Therefore, the significant decrease in variance of Murrumbidgee River heights at Hay and in mean and variance at Wagga Wagga, is most likely due to the April-May reduced dam inflows and precipitation, and from average JJAS catchment precipitation since 1991. Any role played by water extraction for irrigation between Wagga Wagga and Hay, where irrigation is concentrated, is likely to be small owing to the highly significant mean river height reduction at Wagga Wagga which is upstream from Hay. However, irrigation, and other water usage, is sourced from the dams, so there is a long-term impact of irrigation over the months preceding JJAS on flows at Wagga Wagga, due to the reduction in water stored in the upstream dams. The dams integrate the www.nature.com/scientificreports/ water extracted for irrigation and all other usage since the last spill event, and therefore the extractions over an extended period can have an impact on when the dam will fill, and hence on the flows downstream, including Wagga Wagga. The minimum water level in the dams, which typically occurs near the end of Autumn, is due to the reduction in inflows (impacted by climate change), and extractions from the dam. Coupled with the tendency for a slower fill rate due to reduced inflows, this results in fewer spill events. As a consequence, there is a change in the distribution between spill events and irrigation releases, changing the frequency distribution of flows. This will be particularly the case for the JJAS period, as a delayed dam filling will have a major impact in dam levels in that period. Before the 1990s the river at Wagga Wagga and Hay reached flood level height or close to flood level height regularly in JJAS from precipitation-driven inflows regardless of the amount of water that was extracted (see Fig. 2a,b). Since the mid-1990s less water reaching Wagga Wagga has significantly reduced the river height owing to significantly decreased April-May and JJAS precipitation-driven inflows at the upstream catchment dams of Burrinjuck and Blowering Dams and significantly decreased mean precipitation at Tumut, which also represents the catchment area of Blowering Dam (Fig. 2c,d; Table 4). Moreover, there has been no overallocation or hoarding of water found in the southern MDB 24 . In a different southern MDB catchment study of the Millennium Drought 1997-2008, factors for a disproportionate reduction in rainfall run-off were reduced mean annual rainfall, less interannual variability of rainfall, changed seasonality of rainfall and lastly increased potential evaporation 25 . However, the last two factors mentioned have since become well established www.nature.com/scientificreports/ in the last decade with reference to the work in this study. It was suggested that a rainfall reduction alone does not explain the observed inflow reduction trend 26 . Even after a major rain event, the soils are so dry that they absorb more water than before the rain event, and less reaches the dams and rivers than on a wet catchment. In the last three decades it is unknown what the effect on run-off into dams and the Murrumbidgee river has been in JJAS from major rain events because, apart from August-September 2016, there have been no major catchment net inflows since 1991 (Fig. 2c,d). There were significant precipitation-driven inflows during SON 2010 which led to flood level exceedances at Wagga Wagga and Hay in December 2010. In June and July 1991 there was a series of rain-producing cut-off low pressure systems over inland NSW and the adjacent coast influencing the catchment, interspersed with persistent, precipitation-producing frontal systems embedded in the westerly airflow during July and August. Rain producing inland cut-off low pressure systems over southeast Australia are the main influence on enhancing JJAS rainfall totals 8 . www.nature.com/scientificreports/   www.nature.com/scientificreports/ Decreased JJAS precipitation in continental southeast Australia has been evident for at least the last two decades, as anticipated by climate scientists. The naturally periodic La Niña phenomenon provided spring and summer precipitation during much of 2010 to 2012, which ended the Millennium Drought (1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)). The only other recent widespread significant rainfall in southeast continental Australia was in August-September 2016 due to a negative phase of the Indian Ocean Dipole (IOD). A negative IOD phase typically is associated with wetter than normal spring conditions for southeast Australia 7,8 .
Although the SAM is an atmospheric index with a time scale typically of a few weeks, an annual average SAM reconstruction shows that since the 1970s it is in its most positive state over at least the past 1000 years 27 . Prior to the 1990s soil wetness would have been in phase with the annual cycle of winter/spring peak rainfall, dry summer/early autumn and without a long term trend in SAM. However, because SAM has trended positive since the 1970s, the annual cycle of soil wetness of the MDB has been increasingly disrupted particularly since the Millennium Drought 23 and there is also a potential long-term impact from groundwater systems 28 . This is supported by the most recent available annual area-averaged actual evapotranspiration and soil moisture deciles in Fig. 7a,b. These figures show the anomalously dry MDB catchment area in the period 2018-2019. In the southeast corner of the MDB, actual evapotranspiration is below average and soil moisture is very much below average. Two Supplementary Tables showing historical April-May (S1) and JJAS (S2) maximum river heights above flood level at Hay (6.7 m) in IPO phases indicate, as expected, more in negative IPO phases than positive phases and importantly a dissociation with the IPO resulting from none in the most recent negative phase from 1998 and the preceding positive phase after the early 1990s. The implication is that accelerated global warming since the 1990s has overwhelmed the influence of negative IPO on precipitation.
The MDB plan, introduced from 2013 29 provided, for the first time, regulated allocations to environmental flows for ecosystem sustainability of rivers in southeast Australia such as the Murrumbidgee. However, the plan requires that each year on 1 July a fixed amount of water is locked in for future consumption, split three ways with the highest priority for human consumption and irrigation for permanent crops (e.g., fruit trees and nuts). The remaining allocations are split between non-permanent crops (e.g., cotton, rice) and environmental flow. A Table 4. P-values from permutation testing differences in interval means and variances for, April-May and JJAS precipitation for Blowering Dam, Burrinjuck Dam, and Tumut; inflows into Blowering Dam and Burrinjuck Dam; Murrumbidgee River heights at Blowering Dam, Burrinjuck Dam, Wagga Wagga and Hay. Significant values (p < 0.10) are italicized with an asterisk. Note that the p-value for each variance test is calculated after one sample has had bias correction in the mean. Note that the p-value for each variance test is calculated after one sample has had bias correction in the mean. Key points to note are highly significant April-May precipitation decreases from 1965 to 1991 to 1992-2018, highly significant recent JJAS net inflows into Blowering and Burrinjuck Dams, and some significant recent decreases in mean or variance of JJAS river heights at Wagga Wagga and Hay.

Catchment location Precipitation
Net Inflows (% Volume) Murrumbidgee River level  www.nature.com/scientificreports/ major issue is that the forecast net inflows upon which the allocations are based are the minimum inflows experienced in the 120 years up to the end of the twentieth century. However, as shown, even lower inflows have been experienced in the past two decades. It is not surprising that there is a significant decrease in the JJAS variance of the Murrumbidgee River height at Hay (p-value = 0.005) and both the JJAS mean and variance at Wagga Wagga from the periods 1965-1991 to 1992-2018 (mean p-value = 0.0044, variance p-value = 0.095; Table 4) since this period corresponds with the significantly reduced mean April-May catchment precipitation and mean April-May dam net inflows. The fact that there has been no significant change in the mean Murrumbidgee River height since 1991 is an indication that there has been a lack of major April-September rain events. The lack of significant catchment rainfall events from April to September is the reason for the reduction in the mean and variance of river heights at Wagga Wagga. Floods in April-May are rare along the Murrumbidgee River and the six years since 1874 in which April-May floods occurred at Hay prior to 1991 (Table 3), were dominated by precipitation that occurred as a result of mid-latitude interaction with either tropical or subtropical moisture, whereas the last flood that occurred in March 2012, was the result of a rain-producing tropical low pressure trough in the easterly wind regime that extended from northwest Australia to a low pressure centre in southern New South Wales near the Murrumbidgee catchment. Moreover, given the significant decline in April-May, Murrumbidgee catchment rainfall, JJAS run-off into the dams and Murrumbidgee River height at Wagga Wagga since 1991, the implication for water allocations of irrigated agriculture downstream from Wagga Wagga and for flood plain environmental flows required for sustainable wetlands downstream from Hay, will continue to be a problem.    www.nature.com/scientificreports/  www.nature.com/scientificreports/ April-May precipitation, decreased April-May and JJAS dam net inflows, and increased mean temperatures, which represents the accelerated global warming signal since the 1990s, all reduced catchment area run-off. Future work therefore is planned to address the high priority of searching for attributes, to assist in identifying and understanding the role of the key meteorological drivers, especially those related to the April-May decrease in catchment rainfall. Regardless of the mandated environmental flows in the last decade and the annually determined sustainable extraction limits on irrigation, water availability from the Murrumbidgee River at Wagga Wagga and Hay over recent decades continues to be affected by changes in catchment precipitation and run-off. In this study the decreased river heights at Wagga Wagga and Hay in the 27-year period, 1992-2018, occurred as a result of a change in seasonality of rainfall and increase in potential evaporation during the current accelerated period of global warming as described above, suggesting the need for a new review of water availability and sustainability in the Murrumbidgee River system and also of other river systems in the southern MDB.

Methods
The available data sets were Murrumbidgee River heights at Hay in the southern MDB, monthly precipitation, and maximum and minimum screen temperatures at four stations at Blowering Dam, Burrinjuck Dam, Tumut and Cabramurra. The aims of the methodology were to identify trends in the data sets and to use wavelet analysis to identify possible climate drivers of the trends. In addition, ML techniques used to assess the importance of a wide range of climate drivers related to catchment precipitation, are described briefly.
Trend analysis. The time series data were first plotted with their percentiles to obtain an overview of any trends which might be present in the data. Each data set was grouped into two equal 27-year periods, by taking account of the accelerated global warming signal from the mid-1990s and the archived length of each variable. Hence, the ending period was 1992-2018, and bootstrap resampling was applied with 5000 resamples. Box plots of the mean and variance of the bootstrapped intervals were created which provide a deeper understanding of any trends that might be present in the data. Permutation testing was applied with replacement to test for statistical significance in any apparent changes between two 27-year periods.
Wavelet analysis. Wavelet analysis 30,31 was applied to each time series to detect potential climate drivers such as the El-Niño Southern Oscillation. This approach provides both the local wavelet power spectrum (e.g., Fig. 6a, left panel) and the global power spectrum (e.g., Fig. 6a, right panel). The local wavelet power spectrum shows how the influence of climate drivers changes over time, while the global power spectrum provides an overview of which drivers are dominant in the time series. In this study, we used the Morlet wavelet as the mother wavelet.

Machine learning techniques.
Multiple linear regression, support vector regression (with both radial basis function and polynomial kernels), and random forests were used with either forward or backward selection through the space of potential attributes to determine the key attributes of precipitation in the Murrumbidgee catchment. These techniques were employed in a detailed study of the Sydney catchment 32 . By applying tenfold cross-validation on the data set, which covers the years 1965-2018, those attributes that appear in at least 50% of the folds across the eight training methods are considered key attributes of precipitation.

Data availability
The rainfall and temperature station datasets in this study are available online at the Australian Bureau of Meteorology at: http:// www. bom. gov. au/ clima te/ change/ index. shtml# tabs= Datas ets. The digitized weekly river height data at Hay are available at:-https:// doi. org/ 10. 5281/ zenodo. 37794 90. Data sets for the SAM, IOD, IPO, ENSO and Global SST time series data were obtained from the NOAA-ESRL Physical Sciences Laboratory, Boulder Colorado from their Web site at: http:// www. esrl. noaa. gov/ psd/ gcos_ wgsp/ Times eries/.