The methodology we adopt is based on the standard NRFA trend testing approach outlined by Hannaford et al. (2021). In brief, monotonic trends were assessed using the Mann-Kendall (MK) test (Mann 1945; Kendall 1975), a non-parametric rank-based approach that is widely supported for use in streamflow analysis (e.g. Hannaford & Marsh 2008; Murphy et al. 2013). The MK Z statistic (MKZ) follows the standard normal distribution with a mean of zero and a variance of one. A positive (or negative) value of MKZ indicates an increasing (or decreasing) trend. The probability of Type 1 errors set at the 5% significance level allowed the evaluation of statistical significance. A two-tailed MK test was chosen; hence, the null hypothesis of ‘no trend present’ (increasing or decreasing) is rejected when MKZ is outside ±1.96 using traditional statistical testing.
The analysis was undertaken on the full AMAX record available at each station (with the exclusion of rejected and unrepresentative periods). As is widely noted in the literature (e.g. MacDonald & Sangster 2017), decadal-scale variability results in flood-rich and flood-poor periods that can affect the robustness of trends based on the period of data that is selected, users should be aware of this and note that trends have been carried out on the available data and may not reflect longer-term trends (outside of the period of observed data).
The trend analysis is repeated with every release of the NRFA Peak Flow Dataset and run on stations that are ‘Suitable for Pooling’ or ‘Suitable for QMED’ .
The stations were selected based on the following missing data criteria to mitigate the impacts of inevitable gaps in data:
- No more than 10% of missing data.
- 27 or more years of data (≥27 AMAX).