Change in beta diversity of riverine fish during and after supra-seasonal drought

A core theme in ecohydrology is understanding how hydrology affects spatial variation in the composition of species assemblages (i.e., beta diversity). However, most empirical evidence is from research in upland rivers spanning small spatial extents. Relatively little is known of the consequences of hydrological variation for beta diversity across multiple spatial scales in lowland rivers. We sought to examine how spatial variation in hydrology and fish beta diversity within and among rivers changed over time in response to intensification and cessation of hydrological drought. We used monitoring data of fish assemblages, coupled with hydrological and biophysical data, to test how spatial variation in hydrology and multiple components of fish beta diversity in lowland rivers of the Murray—Darling Basin (Australia) varied across spatial scales during contrasting hydrological phases. Spatial variation in hydrology among rivers declined with increasing duration of drought before increasing during a return to above-average flows. Spatial variation in hydrology within rivers did not show consistent changes between hydrological phases. Beta diversity among and within rivers showed variable, river-specific changes among hydrological phases for both incidence- and abundance-based components of assemblage composition. Inconsistent hydrology—beta diversity patterns found here suggest that mechanisms and outcomes of drought and flooding impacts to beta diversity are context-dependent and not broadly generalisable. Our findings indicate that hydrological fluctuations occurring in the Murray—Darling Basin in the period analysed here did not cause significant or consistent homogenisation or differentiation of freshwater fish assemblages.


Introduction
Ecological research is increasingly expanding to understanding and predicting the factors responsible for variation in species assemblages across multiple spatial scales (McGill et al. 2015;Pound et al. 2019;Chase et al. 2020). A fundamental theme in this expansion is spatial variation in the composition of assemblages or 'beta diversity' (Anderson et al. 2006;Socolar et al. 2016;Bush et al. 2016;Soininen et al. 2018). Beta diversity quantifies the similarity or difference in communities across scales and can be applied to different facets of biodiversity: taxonomic (emphasising species identity), functional (ecological traits), or phylogenetic (evolutionary history and potential) (Rocha et al. 2018;Carvalho et al. 2020). Our understanding of the factors affecting beta diversity across both terrestrial and aquatic ecosystems is largely based on fixed or 'static' differences in ecological variables across landscapes (e.g., Heino et al. 2015 for review). We know far less about how beta diversity varies over time time in response to temporally variable natural and anthropogenic factors (Cañedo-Argüelles et al. 2020;Crabot et al. 2020;Li et al. 2020;Lindholm et al. 2020).
Hydrological variability in space and time is a driver of beta diversity of freshwater, riparian, and floodplain biota via multiple ecohydrological mechanisms (Bogan et al. 2013;Rolls et al. 2018;Crabot et al. 2020). Variation in hydrological connectivity influences beta diversity by determining organism dispersal among patches of habitat, with floods frequently reducing beta diversity by increasing connectivity, and low water levels increasingly fragmenting patches of habitat leading to increased beta diversity as local assemblages diverge (Larned et al. 2010;de Macedo-Soares et al. 2010). Additionally, spatial variation in hydrology drives beta diversity because differences in hydrological regimes favour selection of contrasting species assemblages across space (Lytle and Poff 2004;Leigh and Sheldon 2009;Bogan et al. 2013). When spatial 'hydro-diversity' declines, such as during broad-scale hydrological drought, increasingly harsh environmental conditions trigger community disassembly (e.g., Chase 2007;O'Neill, 2016) by eliminating intolerant species from individual locations up to the extent of entire regions. Changes in the spatial variation in riverine hydrology over time are predicted to reduce or increase beta diversity ('biotic homogenisation' and 'biotic differentiation', respectively;Olden and Rooney 2006), possibly by changes in spatial 'hydro-diversity' and loss and recovery of species during hydrological fluctuations. Most empirical evidence of hydrology-beta diversity linkages is from studies spanning either headwater streams or single river systems (see Rolls et al. 2018). In comparison, we are unaware of empirical evidence of the effect of hydrology on beta diversity in a multi-scaled context across large, lowland river systems.
Climate fluctuations cause temporal variation in hydrology and provide opportunities to examine how freshwater beta diversity is driven by hydrology across multiple spatial scales. Water-resource development has been the predominantly perceived cause of altered spatial variation in hydrology; hydrological regimes can either increase or decrease in similarity among and within rivers because of dams and water diversions (Thoms and Parsons 2003;Poff et al. 2007;McManamay et al. 2012;Mackay et al. 2014;Rolls et al. 2021). However, hydrological fluctuations (drought and flooding) imposed by climate variability may also lead to periods of spatial hydrological homogenisation or differentiation, which in turn drive changes in riverine beta diversity through spatial variation in assembly-disassembly processes (e.g., Thomaz et al. 2007;Leigh and Sheldon, 2009;Larned et al., 2010). Assessing how climate fluctuations affect beta diversity of riverine assemblages would help to understand and predict hydrological scenarios that lead to biotic homogenisation or differentiation across spatial scales. Such evidence is necessary because change in beta diversity is a fundamental component of global biodiversity change (McGill et al. 2015), especially for freshwater fish (Villéger et al. 2011). Yet most empirical research on change in beta diversity in freshwaters has tested effects of non-native species invasions (e.g., Leprieur et al. 2008;Taylor 2010;Toussaint et al. 2014;Sommerwerk et al. 2017), with only six of 116 studies published to May 2020 testing the effects of temporal hydrological variability (Rolls, unpublished data).
One approach to test the effect of hydrology on beta diversity of river biota is to examine how spatial variation in hydrology differs among distinct hydrological phases, and the associated patterns of beta diversity across multiple spatial scales. Australia's Murray -Darling Basin spans multiple climatic zones contributing to spatial variability in hydrology among its major rivers. Inter-annual variation in rainfall and temperature in the basin is substantial, with prolonged periods of below-average annual rainfall and above-average annual temperature punctuated with wetter periods leading to distinct deviations in river discharge as low flows and floods, respectively. The fish fauna of the Murray -Darling Basin comprises species that have contrasting linkages between hydrology and population dynamics (e.g., Baumgartner et al. 2014). Studies of beta diversity of fish assemblages in the basin have revealed (i) substantial spatial variation in composition associated with variation in hydrology among rivers (e.g., Gehrke et al. 1995), (ii) that assemblage composition in individual rivers shifts between hydrological events such as floods (e.g., Rolls and Wilson 2010;Rayner et al. 2015;Cruz et al. 2020), and (iii) that spatial variation in assemblage composition within rivers differs in association with low-and high-flow events (e.g., Balcombe et al. 2006Balcombe et al. , 2011. Elsewhere (e.g., in the USA), declining surface water availability was associated with both differentiation and homogenisation of fish assemblages at different spatial scales (Taylor et al. 2019;Rogosch and Olden 2019). The scalespecific effects of temporal variation in hydrology on beta diversity have not yet been explored across large rivers such as those of the Murray -Darling Basin.
The objective of this study was to examine how beta diversity of lowland fish assemblages varied among multiple spatial scales over nine years of contrasting hydrological phases driven by climate fluctuations (Table 1). We first examined how spatial variation in hydrology differed among and within rivers of the Murray-Darling Basin during multiple phases of hydrological drought and flooding. Second, we tested how beta diversity of fish assemblages varied among hydrological phases to determine whether drought or post-drought periods were associated with apparent homogenisation or differentiation. We examined beta diversity in the context of replacement and richness difference components of total incidence-based beta diversity (Podani and Schmera 2011;Legendre 2014) in conjunction with temporal variation in abundancebased compositional heterogeneity. Examining these different components of beta diversity is informative because evidence is scant regarding whether temporal changes in beta diversity in response to hydrological variation manifest as changes in total beta diversity or replacement, richness-difference, or abundance components (Ruhí et al. 2017;Taylor et al., 2019;Crabot et al. 2020), particularly for fish assemblages. Finally, we tested how hydrological variability affected the relative contributions of antecedent hydrology, distances among pairs of sampling sites, and local environmental variables to each aspect of beta diversity to determine whether assemblage-environment relationships were variable or consistent across contrasting hydrological conditions.

Study area
The Murray -Darling Basin spans 1,061,000 km 2 of south-eastern Australia and is drained by 23 major rivers (Fig. 1). Most of the basin is low gradient, yet portions along its eastern boundary are undulating and steeper terrain. The climate varies across the basin, transitioning from hot to cool semi-arid (Köppen climate classification: Bsh, Bsk) from the north to the south in the western half, and from humid subtropical (Cfa) to oceanic (Cfb) in the eastern half, with small patches of warm-summer Mediterranean (Csb) along the southern perimeter. Consequently, annual rainfall increases from west to east, with the seasonality of rainfall switching from summer (November-March) dominated in the north to winter (May-October) dominated in the south. The Murray -Darling Basin is Australia's major agricultural region, producing 39% of the nation's total agricultural production (ABS 2008).
Natural hydrological variability within rivers of the Murray-Darling Basin is substantial (Kennard et al. 2010), with discharge seasonality switching from predominantly late summer in the northern portion (i.e., rivers draining into the Darling River) to late winter and early spring in the south (rivers draining into the Murray River) (MDBC 2008). Natural hydrological regimes have been altered substantially by the combined effects of large dams (generally located in upland and headwater regions) and both water extraction and irrigation diversion (particularly in lowland regions). Anthropogenic impacts on hydrological during a return to above-average flows. Spatial variation in hydrology within rivers did not show consistent changes between hydrological phases 2. Change in riverine fish beta diversity during increasing drought duration and initial post-drought phases PERMANOVA, PCoA, PERMDISP and group-boxplots of distances to group centroids tested using total and decomposed incidence-based dissimilarity and abundance-based dissimilarity matrices Beta diversity of fish assemblages among and within rivers differs among hydrological phases, with pronounced homogenisation of fish assemblages (i) within rivers occurring during a flood period following drought and (ii) among rivers during prolonged drought Fish beta diversity among and within rivers showed variable, river-specific changes among hydrological periods for all incidence-and abundance-based components of assemblage composition, with no consistent response occurring in all rivers 3. Assess the relative statistical associations of antecedent hydrology, spatial variables and local environmental variables with beta diversity of fish assemblages Variation partitioning and db-RDA using distance-based Moran's Eigenvector Maps, local environmental variables, and antecedent hydrology against each incidence-and abundance-based component of beta diversity for each hydrological phase The relative statistical contribution of antecedent hydrology to components of beta diversity would be greater during the prolonged-drought phase than during the initial-drought or post-drought phases Antecedent hydrology explained a small but significant component of fish beta diversity in the initial-drought and post-drought phases, yet it explained either none or very little of beta-diversity during prolonged drought regimes have occurred as decreased magnitude and frequency of flooding, altered discharge seasonality (particularly in the southern portion), and increased baseflow magnitude, overall leading to 'temporal flow homogenisation' at intra-annual and inter-annual time scales (Maheshwari et al. 1995;McMahon and Finlayson 2003;Humphries et al. 2008). While a detailed analysis of the effects of anthropogenic hydrological alteration on the magnitude of spatial variation in hydrological regimes within and among rivers of the Murray-Darling Basin has not been done, water-resource development has reduced longitudinal hydrological variability in some rivers in the basin's north (Thoms and Parsons 2003). During 2002-early 2010, the Murray-Darling Basin had a protracted period of below-average rainfall, leading to increasing hydrological drought ('Millennium Drought'), followed by a La-Niña event of substantial rainfall and flooding during late 2010-2012 (Fig. 2). We used a nine-year period spanning this prolonged drought and return to wetter conditions to examine how beta diversity of riverine fish assemblages varied over time. We distinguished hydrological phases using temporal analysis of the annual standardised discharge averaged across 198 stream flow records from 1960 to 2018. Hydrological phases were distinguished via trends in standardised discharge as 'initial drought ' (2004-2006), 'prolonged drought ' (2007-2009), and 'post-drought' (2010-2012) (Fig. 2).

Data sources
Fish: Fish were sampled in each of the 23 major river systems at three-year cycles during 2004-2013 as part of a "river health" monitoring program (Davies et al. 2010;MDBA 2012). Sites were selected with a stratified random design, whereby rivers were stratified into elevation zones (i.e., lowland, slope, upland, and montane), and the locations of sampling sites within each zone were randomly selected with GIS. The number of sites in each zone differed among rivers because of variation in river network length. At each site, fish were sampled with a standardised protocol of boat and/or backpack electrofishing combined with unbaited shrimp traps. If randomly selected sites were unable to sampled effectively at the time of sampling (e.g., reach was dry, access was unable Fig. 1 Location of sites (circles) where fish assemblages were sampled from eight major river systems of the southern and northern Murray-Darling Basin during November 2004-May 2013. Filled circles are sites that were paired with nearby hydrological gauging stations (see text for details) and used for modelling spatial, environmental, and hydrological drivers of beta diversity across seven of the selected study rivers. Unfilled circles indicate sites that were not close enough to be paired to hydrological gauging stations (and therefore excluded from beta diversity modelling), yet were included in broader analyses of temporal change in beta diversity within and among the eight study rivers to be gained), replacement sites were sampled (see Chessman 2018). Of the 23 river systems sampled, we selected the Border Rivers, Condamine, Darling, Gwydir, Lachlan, Macquarie, Murrumbidgee, and Mean annual discharge is presented as a proportion of the long-term mean annual discharge recorded between 1960 and 2018 for sites with complete gauging records. Hydrological phases were defined as 'initial drought ' (2004-2006), 'prolonged drought ' (2007-2009), and 'post-drought' (2010-2012), represented by the yellow, orange, and green bars respectively. The horizontal lines in each plot indicate mean values for each variable Namoi rivers as they span the major climate zones of the Basin and ≥ eight sites were sampled within the slope and lowland zones of each river during each hydrological phase. Of these eight selected rivers, the Lachlan and Murrumbidgee are part of the 'southern basin' and the Border Rivers, Condamine, Darling, Gwydir, Macquarie, and Namoi rivers represent the 'northern basin' (Fig. 1).
Hydrology: We sourced daily discharge data (ML day −1 ) from 1990 to 2013 from government water agencies (WaterNSW and the Queensland Department of Natural Resources, Mines and Energy) for all streamflow gauging stations that were located within 20 km of fish monitoring sites. Gaps in records (mean days across gauging records = 1.34%, maximum = 7.31% of all days during 1990-2014) were infilled by multiple regression from the nearest gauge record within each river using the River Analysis Package (Marsh 2004), providing the R 2 value was ≥ 0.80 (stream gauge data were discarded if there were gaps that could not be filled confidently with an R 2 < 0.80).
Spatial and local environmental variables: Spatial variation in riverine fish assemblages is influenced by a variety of spatial, environmental, and biotic factors (e.g., Jackson et al. 2001). We matched study sites to segments of the Australian national stream (GIS) layer using ArcMap, to extract values for environmental variables as the Australian Stream Environmental Attributes Database (Stein and Nevill, 2011). We selected a subset of these environmental variables for each sampling site to estimate the physical, climatic, and vegetation cover (reach elevation, stream order, proportion of natural forest within the stream segment, average annual mean air temperature, average coldest-month minimum temperature, and average hottest-month maximum temperature calculated during .

Data processing and variable calculation
Fish: We extracted all fish data produced from sampling in these eight river systems during 2004-2013 from the Murray -Darling Basin Authority Fish and Macroinvertebrate database (www. data. gov. au) to produce a sample × species (count) matrix. Two different sets of analyses were done to examine fish beta diversity patterns and beta diversityenvironmental relationships among hydrological phases. The first made use of all fish data from the eight study rivers recorded during the study period, and the second used a 'reduced' dataset where fish assemblage data were used only if fish monitoring sites were matched with hydrological gauging station data.

Hydrology
For each streamflow record, we calculated values of eight hydrological metrics for the three-year period (Online Appendix S1) antecedent to fish sampling. We chose this period to match the time interval between sampling events, and because we were specifically interested in how the magnitude of spatial variation in hydrology within and among rivers changed over time. Hydrological metrics reflected the mean duration, frequency, and timing of overbank flows (defined as the two-year annual return interval from the annual flood series) and low-flow events (less than the long-term 90th percentile flow), and the frequency and duration of in-channel flow pulses (defined as events > 0.5 the long-term mean daily flow). Collinear variables were identified if they had Variance Inflation Factors (VIFs) greater than five (5) using the 'vifcor' function in the USDM package. A single variable, 'duration of in-channel flow pulses', was highly correlated to other variables and was removed prior to all subsequent analyses.

Spatial variables
We used the origin-destination function in the network analysist tool of ArcMap 10.6 to calculate the network distance between pairs of sites based on the river network to create a pairwise distance matrix of watercourse distance between all pairs of study sites (reflecting the distances between pairs of sites relevant for dispersal by strictly aquatic organisms; Heino et al. 2017). From this distance matrix, spatial variables were created using distance-based Moran's Eigenvector Maps (dbMEM) , retaining only positive and significantly spatially correlated eigenvectors. No further manipulation of reach-scale environmental variables was needed prior to analyses.

Temporal changes in spatial variation in hydrology within and among rivers
Prior to analyses, we standardised each hydrological variable (mean = 0, SD = 1) to convert all variables to the same scale and prevent any variable having an overwhelming influence on patterns in the dataset by having a broader range of values than another. Dissimilarity in hydrological metrics for pairs of samples was quantified by Euclidean distance. We used a two-factor Permutational Analysis of Variance (PERMANOVA; Anderson 2001) with the 'adonis' function in vegan to test how spatial variation in hydrology among rivers ('River'; fixed) varied among hydrological phases ('Phase'; fixed). PERMANOVA tests for 'location' (i.e., centroids) differences among groups of samples. Therefore, we calculated the mean distance between group centroids for groups of sites from each river during each three-year hydrological phase (n = 24; eight rivers x three periods) using the 'meandist' function in the vegan package. Centroids of each group were displayed with Principal Coordinates Analysis (PCoA). Using the same Euclidean distance matrix, we then tested how within-river variation in hydrology varied among hydrological phases using tests of homogeneity of dispersions (PER-MDISP; Anderson 2006). Patterns of within-river variation in hydrology were displayed in box plots (x-axis: time, y-axis: mean pairwise distance among group centroids). PERMANOVA and PERMDISP analyses were done using the 'adonis' and 'betadisper' functions in the vegan package (Oksanen et al. 2019), and PCoA ordinations were produced with the ape package (Paradis et al. 2019) in R 4.0.0 (R Core Development Team 2020).

Temporal changes in beta diversity of fish assemblages
We produced four dissimilarity matrices from the site × species matrix, each quantifying different aspects of incidence and abundance composition to determine the underlying characteristics of beta diversity over time (sensu Anderson et al. 2011). Sørensen's distance (1) was used to calculate a 'broadsense' dissimilarity in incidence-based composition among samples, which was then decomposed into the replacement (2) and richness-difference (3) dissimilarities using the 'beta' function in the BAT package (Cardoso et al. 2015). A fourth dissimilarity matrix used was abundance-standardised assemblage dissimilarity (percentage-difference or 'Bray-Curtis' distance; Legendre and De Cáceres 2013). We retained all species not native to the Murray -Darling Basin or possibly translocated among study rivers, because non-native species may be important manifestations of changes in beta diversity across scales (Olden and Poff 2003).
For each dissimilarity matrix, we used a combination of Permutational Multivariate Analysis of Variance (PERMANOVA; Anderson 2001) and tests of homogeneity of multivariate dispersions (PER-MDISP; Anderson 2006) to determine how spatial variation in fish assemblage composition varied over time among and within river systems. This combination of statistical tests helps to determine the specific scale(s) of variation in beta diversity over time (see Anderson et al. 2011). We displayed patterns in beta diversity with both PCoA ordinations and graphing distances from group centroids for each group over time. To assist in interpreting changes in beta diversity among hydrological phases (both within and among rivers), we (i) plotted values for recorded sitelevel species richness (alpha diversity) as grouped boxplots, (ii) tested differences in alpha diversity using PERMANOVA based on a Euclidean distance, and (iii) presented summary data of gamma diversity (γ), mean alpha diversity (α), and Whittaker's beta diversity (as γ/ α) (see Supplementary Material S3).

Assessing the association of antecedent hydrology and spatial and environmental variables with beta diversity across hydrological periods
We used variation partitioning in constrained ordination (Borcard et al. 1992;Peres-Neto et al. 2006) to assess the relative importance of antecedent hydrology and spatial and local environmental variables in explaining variation in each of the four distance matrices of fish assemblages during each hydrological phase. Predictor variables were grouped into 'spatial', 'local environmental', and 'antecedent hydrology' sets. For each set of predictors, we used forward selection ('ordiR2step' in vegan) to determine which specific variables to retain in subsequent analyses (Blanchet et al. 2008).
Vol.: (0123456789) Variation partitioning identifies the unique ('pure') and combined contribution of predictor variables in explaining variation in assemblage composition (Borcard et al. 1992;Anderson and Gribble 1998). Variation partitioning and the significance of each pure fraction were assessed by distance-based redundancy analysis (dbRDA; Legendre and Anderson 1999) with the 'varpart' and 'capscale' functions in vegan (Oksanen et al. 2019). The statistical significance of each pure fraction was tested with the 'anova' function.

Spatial variation in hydrology across hydrological periods
Differences in hydrology among rivers were inconsistent among hydrological phases (significant River x Phase interaction; Table 2). Overall, spatial variation in hydrology among rivers was significantly lower during the prolonged-drought phase than during both initial-drought and post-drought phases (PERMDISP test of centroids P = 0.048), as evidenced by clustering of river centroids in the Principal Coordinates Analysis (Fig. 3a). Within-river hydrological variability was similar among the three hydrological phases for all rivers (Fig. 3b) except for the Condamine River, where spatial variation in hydrology was lowest during initial drought, increased during prolonged drought, and was highest after the drought.
Scale-dependent change in beta diversity of fish assemblages associated with hydrological phases Beta diversity among rivers was inconsistent across hydrological phases for each aspect of assemblage composition except the richness-difference component (significant interaction term for River × Hydrological phase; Table 3). Ordination of total incidence-based beta diversity (Sørensen) revealed that fish assemblages in the northern Murray-Darling Basin (Border Rivers, Condamine, Darling, Gwydir, Macquarie and Namoi rivers) were most dissimilar during initial drought, increasingly similar during prolonged drought, and most similar during the postdrought phase (Fig. 4). In contrast, mean assemblage composition of both southern (Lachlan and Murrumbidgee) rivers changed minimally between hydrological phases. Decomposition of overall incidence-based beta diversity into replacement and richness-difference components indicated that lower beta diversity during the post-drought phase was driven by replacement rather than richness differences (Fig. 4). This was further supported by the fact that the richness difference component of beta diversity differed consistently among rivers in all three hydrological phases (non-significant River x Phase interaction; Table 3). For abundance-based assemblage composition, both rivers in the southern basin showed little change in mean composition across hydrological phases, whereas average composition of the remaining rivers varied substantially among hydrological phases (Fig. 4). For each dissimilarity matrix, significant differences among rivers explained the largest proportion of variation (14.7-22.7%), and while the interaction of River × Hydrological phase was usually significant, it contributed only a moderate amount of explained variation (Table 3).
Overall tests of homogeneity of dispersions (PERMDISP) found that beta diversity within rivers differed significantly among combinations of rivers and hydrological phases for each aspect of assemblage composition (Table 4). For overall incidence-based beta diversity (Sørensen dissimilarity), post-hoc pairwise tests determined that within-river beta diversity in the Gwydir River was significantly lower (P = 0.006) during prolonged drought than during initial drought. In both the Condamine and Darling rivers, within-river beta diversity was significantly lower during the post-drought period ( Fig. 5; Table 2 Multivariate PERMANOVA tests (two-factor fixedeffects model) of differences in mean river-scale hydrology calculated using Euclidean distance of eight hydrological variables among rivers, hydrological phases, and their interaction Significant differences (P ≤ 0.05) are in bold, and variance components (VC) are expressed as the percentage of total variation explained by each level in the model post-hoc pairwise tests: Condamine P = 0.0008, Darling P = 0.024). Decomposing overall beta diversity into replacement and richness-difference components showed that these changes in within-river beta diversity among hydrological phases were driven primarily by replacement rather richness differences among sites (Fig. 5). Alpha diversity among hydrological phases showed inconsistent patterns among rivers, with the Lachlan, Macquarie, Murrumbidgee, and Namoi rivers showing no significant change in alpha diversity, the Border Rivers and Darling rivers having higher richness during the post-drought phase (compared to both drought-phases), the Condamine River showing lower richness post-drought, and the Gwydir River having lower richness during the initial drought phase (Supplementary Material S4, S5). Spatial variation in abundance-based composition within rivers generally did not differ significantly among hydrological phases, the only exception being the Murrumbidgee River where spatial variation in fish assemblages was significantly lower during the postdrought period (c.f. both the initial and prolongeddrought periods) ( Fig. 5; PERMDISP post-hoc pairwise tests P ≤ 0.001).

Assemblage-environment relationships during contrasting hydrological phases
Variation partitioning showed that the contribution of antecedent hydrology and spatial and environmental variables to beta diversity varied among hydrological phases. The relative contribution and significance of each predictor variable group also varied among beta-diversity components (Table 5). Antecedent hydrology explained a small but significant component of fish beta diversity in the initial-drought and post-drought phases, yet it explained either none or very little of beta diversity during prolonged drought. Across all three hydrological phases, hydrology was significantly associated with overall incidencebased beta diversity, yet it had no significant association with the richness-difference component of beta diversity (Table 5). Local environmental variables were significantly associated with most or all aspects of beta diversity across all hydrological phases. In contrast to environmental factors, spatial variables were not significantly associated with any aspect of incidence-based beta diversity during initial drought, yet were significantly associated with both incidence-based and abundance-based beta diversity during the prolonged-drought and post-drought phases.

Discussion
Beta diversity of fish assemblages in lowland rivers of the Murray -Darling Basin varied inconsistently with contrasting hydrological phases depending on the scale of analysis. Both incidence-based and abundance-based beta diversity among rivers were higher during both initial and prolonged drought phases and lower following drought in the northern Murray -Darling Basin. However, differences in assemblage composition between the two southern rivers remained similar through time. Changes in withinriver beta diversity were variable and unique to each river, and depended on whether analyses were based on incidence or abundance composition. These findings suggest that changes in spatial beta diversity associated with temporal variation in hydrology, manifesting as biotic homogenisation or biotic differentiation, are not broadly generalisable to fish assemblages of lowland rivers.

Change in spatial hydrological variability with climate fluctuations
Australia's Millennium drought was the most severe in terms of the magnitude and duration of reduced precipitation to have occurred during the period 1783-2009 (Gergis et al. 2012), contributing to substantial declines in river discharge in the Murray -Darling Basin (Leblanc et al. 2012). River-scale hydrology differed significantly between the postdrought phase and both initial-drought and prolongeddrought phases. Spatial variation in hydrology among rivers was lowest during prolonged drought compared to both initial drought and post-drought phases, indicating that increased duration of broad-spanning drought leads to loss of spatial hydro-diversity. Such reductions in hydro-diversity are expected because increasing duration of broad-scale drought increases the prevalence and simultaneous occurrence of low or zero-flow conditions (Leblanc et al. 2009;Van Loon 2015). We are unaware of other research examining changes in spatial variation in river hydrology across spatial scales linked with climate fluctuations using a Within-river (i.e. longitudinal) hydrological variability was similar among the three hydrological phases for all but one of the rivers where hydrological data were analysed. The exception was the Condamine River, where longitudinal hydrological Ordinations by Principal Coordinates Analysis (PCoA) illustrating differences in average fish assemblage composition among rivers and between hydrological phases based on Sørensen dissimilarity (and replacement and richness differ-ence components) for incidence data and percentage-difference (Bray-Curtis) for abundance data. The percentages of total variance explained by each of the first two PCoA axes are indicated in brackets Table 4 Results of tests of homogeneity of dispersions (PER-MDISP), based on a two-factor interaction (river x hydrological phase) examining change in within-river (i.e. among-site) variation in assemblage composition among hydrological phases Beta diversity aspect F-ratio P variability was lowest during initial drought, increased during prolonged drought and was highest during the post-drought phase. Under modified (regulated) conditions, the hydrological regime of the Condamine River has substantial longitudinal variation (Thoms and Parsons 2003), and our analysis indicates that the magnitude of longitudinal hydrological variability varies with hydrological drought. Increasing longitudinal variability between initial-drought and prolonged-drought phases may have been caused by increased irrigation pressure from agriculture as the drought progressed, leading to increased extraction of connected surface and subsurface water, particularly in lowland reaches (e.g. Wen et al. 2011). Further increases in longitudinal hydrological variability following drought may have been due to either geographical variation in runoff that is intercepted by offriver storages that account for most of the stored water Grouped box plots of distances from group (river × hydrological phase) medians of fish assemblages for eight large rivers of the Murray -Darling Basin among climate phases based on Sørensen dissimilarity (and replacement and richness components) for incidence data and percentagedifference (Bray-Curtis) for abundance data. Each boxplot displays the median (horizontal line), 25th and 75th percentile values and the whiskers (vertical lines) display 1.5 × the upper and lower interquartile range for each combination of river and hydrological phase; increasing distance from group median indicates increasing within-river beta diversity in the Condamine River system, or possible losses of surface water via groundwater recharge that are common in alluvial lowland rivers (sensu Larned et al. 2008). These mechanisms are also likely to operate in the remaining study rivers across the northern Murray -Darling Basin, yet the lack of a consistent effect may be due to differences in extraction demands for water masking effects of climate fluctuations.
Temporal changes in beta diversity associated with hydrological phases across spatial scales Significant changes in among-river beta diversity were associated with contrasting hydrological phases. As Australia's Millennium Drought progressed, seven of 39 fish species present in the Murray-Darling Basin had significant changes in their occurrence (decreased occurrence: two species, increased occurrence: five species) (Chessman 2013). Significant changes in abundance were also detected for five of these seven species (Chessman 2013). These changes in occurrence and abundance suggest that effects of hydrological variation would be reflected as changes in fish assemblage composition as either incidence-based or abundance-based beta diversity. In our study, beta diversity was lowest among the five northern study rivers during the post-drought phase and higher during both initial-drought and prolongeddrought phases. One hypothesis for this pattern is the role of spatial variation in community assemblydisassembly patterns during contrasting hydrological phases (e.g. inundation, discharge recession) (O'Neill 2016). During hydrological drought, loss of flowing water and habitat contraction likely vary substantially among rivers, meaning that assemblages occupying multiple rivers are at contrasting states of disassembly and therefore have increasing beta diversity. Flooding, in contrast, potentially leads to more pronounced consistency in environmental conditions among rivers (and reaches within rivers), likely contributing to lower beta diversity as multiple assemblages are in more synchronised phases of assembly. Overall, change in beta diversity associated with hydrological phases varied among different aspects of assemblage composition. Specifically, abundancebased composition showed a pronounced shift in ordination space between drought phases and the post-drought phase. In contrast, incidence-based Table 5 Summary results of variation partitioning analysis testing the relative importance of spatial, environmental, and antecedent hydrological variables in explaining spatial vari-ation in fish assemblage composition during initial-drought, prolonged-drought, and post-drought hydrological phases The strength of the overall model (adjusted R 2 ) and unique fractions explained by each variable component are presented Negative fractions are not tested and denoted by "-" composition showed minimal changes in overall ordination space among hydrological phases but showed a trend of declining beta diversity among rivers from initial-drought to post-drought. These findings highlight that hydrological variation likely affects different aspects of beta diversity in different ways and combined use of incidence and abundance-based composition is necessary to reveal distinct aspects of multi-species diversity patterns across landscapes (Taylor et al. 2019).
Within-river beta diversity did not change significantly among hydrological phases for most rivers. The lack of a consistent change in within-river beta diversity may have occurred because change in longitudinal hydrological variation was minimal, and therefore often insufficient to affect the distribution and abundance of fish species. Yet, changes in within-river beta diversity among hydrological phases did occur for individual rivers. For example, fish assemblages in the Gwydir River became increasingly homogeneous during prolonged drought and this homogenisation was predominantly driven by changes (i.e., reduction) in species replacement among samples. This decline of beta diversity due to reduced replacement during prolonged drought suggests that some species are 'filtered out' of the regional species pool (sensu Chase 2007). Two species (Australian smelt, freshwater catfish) showed marked declines in occurrence, and one species (bony bream) had marked increases in occurrence in the Gwydir River during prolonged drought compared to both initial-drought and post-drought phases (Online Appendix S2). Of these three species, bony bream is predicted to have low vulnerability to drought (Chessman 2013), and Australian smelt has either moderate vulnerability to drought (Chessman 2013) or intermediate resistance but high resilience to drought (Crook et al. 2010). Overall, the progressive filtering of species during increasingly harsh environmental conditions plays a major role in broader changes in beta diversity during hydrological fluctuations. In contrast to the findings from the Gwydir River, within-river beta diversity was significantly lower in the Condamine and Darling rivers during the post-drought phase than during both initial-and prolonged-drought phases. Lower beta diversity in these rivers during the post-drought phase may have been due to the low storage capacity in reservoirs in these two rivers compared to the other study rivers, possibly contributing to an inability to sustain base flows (and hence hydrological connectivity) during drought (Mallen-Cooper and Zampatti 2020). Increased hydrological connectivity enhances organism dispersal among patches during flooding episodes leading to reduced beta diversity (e.g. de Macedo-Soares et al. 2010;Cruz et al. 2020), and this process possibly explains the decline in within-river beta diversity during the post drought phase in the Condamine and Darling Rivers (Larned et al. 2010).
Previous analyses of effects of hydrological variability on fish beta diversity in the Murray -Darling Basin and elsewhere are based on assessments at the scale of individual rivers and frequently report reduced beta diversity associated with increased hydrological connectivity during flooding or loss of species during drought phases (Balcombe et al. 2006(Balcombe et al. , 2011Rolls and Wilson 2010;Rayner et al. 2015;Miyazono et al. 2015;Cruz et al. 2020). By examining beta diversity -hydrology associations among multiple lowland river systems, our work here indicates that outcomes of hydrological variability (drought, flooding) are variable across river basins. These inconsistent hydrology-beta diversity associations suggest that mechanisms and outcomes of drought and flooding for beta diversity are context-specific and not broadly generalisable. Such conclusions are further supported by the fact that the contribution of antecedent hydrology to assemblage variation relative to the contributions of spatial and environmental variables was highly variable among hydrological phases.
Limitations and caveats applicable to assessing hydrological consequences for beta diversity The absence of consistent significant change in beta diversity of fish assemblages associated with contrasting hydrological phases in the Murray -Darling Basin can possibly be explained by at least four general hypotheses worthy of consideration in future hydrology -beta diversity research. First, the hydrology of all the rivers in this study has long been altered to some degree by water-resource development prior to fish sampling (MDBC 2008). Such hydrological alterations may have reduced the abundance of, or eliminated, flow-dependent species (sensu Baumgartner et al. 2014) from the river systems studied here, with the remaining species being resistant to hydrological variability and therefore unresponsive to drought, in turn contributing to stable beta diversity. A second hypothesis (an alternative to the first) is that while hydrological fluctuations during the study period were substantial, the biota of the Murray -Darling Basin is likely adapted to even more major hydrological variability over evolutionary time scales that has occurred in the basin over millennia (Ho et al. 2015), and hence show limited response to drought events at the magnitude of the one examined here. Third, if randomly selected study sites could not be sampled in this monitoring program (e.g., inadequate depth for sampling or completely dry), they were replaced with sites in the same zone. Replacement of study sites during drought phases possibly dampened observed changes in the occurrence and abundance of fish species (and therefore beta diversity) by concentrating sampling on drought refuges (Elliot 2000;Magoulick and Kobza 2003;Rayner et al. 2009). The fourth hypothesis relates to the timing over which changes in beta diversity are apparent in response to hydrological events such as drought and flooding. Changes in beta diversity may have occurred primarily during the early stages of the Millennium Drought before SRA fish sampling began, or else beta diversity responses may be delayed and become evident over longer time periods. Although some sampling of fish assemblages was undertaken across the Murray -Darling Basin during the 1990s when rainfall, temperature, and discharge were near to average (Gehrke and Harris 2000), the resulting data are insufficient for comparison with those used in the present study. In contrast to responses during the initial phases of drought, changes in beta diversity driven by hydrological events (particularly flooding, yet also drought) may not be immediately evident but rather may be delayed. For example, beta diversity changes associated with post-drought hydrology in this study may be underestimated by monitoring being restricted to the first two-three years after the end of the drought. Understanding the complete linkages between hydrology and beta diversity therefore requires sufficiently long time-series hydrology and fish-assemblage data to explore how beta diversity at a specific point in time is associated with hydrological events occurring over various antecedent periods.

Conclusions and recommendations
Spatial variation in hydrological regimes is potentially altered over time with climate fluctuations, with possible outcomes for biodiversity across landscape scales. Our study found that changes in (i) the spatial variation in hydrology associated with climate fluctuations and (ii) corresponding spatial beta diversity were river specific (Table 1). Overall, these findings indicate that hydrological fluctuations occurring in the Murray -Darling Basin in the period analysed here did not cause consistent homogenisation or differentiation of freshwater fish assemblages. The increasing availability of broad-scale species-assemblage data combined with hydrological data spanning increasingly large extents in Australia and elsewhere provides ongoing opportunities to determine the ecohydrological mechanisms, hydrological scenarios, and scales where change in freshwater beta diversity is expected to occur.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.