Thirty-seven years of freshwater fish and crustacean presence/absence data were obtained from the NZFFD. That was all entries on flowing water dating from January 1970 to June 2007. This consisted of 22,546 sites, over a broad geographic coverage (Fig. 1).

Map showing the locations of the 22,546 river sites with records in the New Zealand Freshwater Fish Database from 1970 – 2007. These sites cover a broad geographic area throughout New Zealand.
The number of sites added to the database has increased over time but there were different patterns of increase related to land cover. The number of sites sampled in exotic forest, urban and scrub land cover showed a gradual increase, approximately doubling every decade. Sampling indigenous forest sites increased at a much greater rate but peaked in the 1990s, while pasture sites increased exponentially over the entire period. Between 2000 and 2007 the number of pasture sites sampled was more than all the other classes combined (Fig. 2).

Graph showing sites added to the New Zealand Freshwater Fish Database (NZFFD) over decades for 1970 to 2007 by the five land cover types:
Sites added to the NZFFD by land cover type, 1970-2007
| Land cover type | 1970-1979 | 1980-1989 | 1990-1999 | 2000-2007 |
|---|---|---|---|---|
| Exotic forest | 55 | 186 | 491 | 591 |
| Indigenous forest | 451 | 1584 | 2098 | 1401 |
| Scrub | 45 | 146 | 525 | 482 |
| Pasture | 384 | 2087 | 3075 | 4390 |
| Urban | 26 | 106 | 306 | 733 |
The average IBI score was significantly higher at indigenous forest and scrub sites than the other land cover classes, and tussock was significantly lower than all other land-cover classes (ANOVA F7,22538 = 247; P < 0.0001) (Table 1; Fig. 3). Pasture sites had the next lowest scores but were not significantly different from urban, exotic and unvegetated (bare land) sites.
| REC class | Pasture | Urban | Exotic forest | Bare land | Tussock | Indigenous forest | Scrub |
|---|---|---|---|---|---|---|---|
| Mean | 29.68 | 30.33 | 30.47 | 31.4 | 17.98 | 36.22 | 36.51 |
| Median | 32 | 34 | 32 | 40 | 0 | 38 | 40 |
| Standard deviation | 17.62 | 16.69 | 19.93 | 20.09 | 19.89 | 17.85 | 19.27 |
| Standard error | 0.18 | 0.49 | 0.55 | 0.88 | 0.37 | 0.24 | 0.56 |
| Number of sites | 9,932 | 1,167 | 1,319 | 522 | 2,806 | 5,530 | 1,194 |

Graph showing average Index of Biotic Integrity (IBI) score for all river sites monitored between 1970 and 2007, grouped by River Environment Classification land-cover type, with standard error (SE).
| Land cover type | Average IBI ± SE |
|---|---|
| Tussock | 17.98 ± 0.37 |
| Pasture | 29.68 ± 0.18 |
| Urban | 30.33 ± 0.49 |
| Exotic forest | 30.47 ± 0.55 |
| Bare | 31.40 ± 0.88 |
| Indigenous forest | 36.22 ± 0.24 |
| Scrub | 36.51 ± 0.56 |
The IBI scores for all sites show there has been a significant decline in indigenous freshwater biodiversity over the past 37 years. This decline was significant for both years and decades (Table 2) with the biggest reduction in the last decade (Fig. 4). To assess which of the land-cover classes contributed to this decline the different classes were analysed separately. Indigenous forest sites showed a significant increase for both years and decades, peaking in the 1990s (Fig. 5; Table 2). Pasture sites showed a significant decrease in IBI scores for both years and decades especially in the last decade (Fig. 6; Table 2). The sites in scrub land cover showed no significant trend over decades but did for years (Fig. 7; Table 2). Urban sites showed a significant decline in IBI scores over the 37 years for both years and decades (Fig. 8; Table 2). The exotic forest sites showed a dip in the 1990s but there was no significant linear trend for both years and decades (Fig. 9; Table 2). Tussock sites showed a significant decline in IBI scores for both years and decades (Fig 10; Table 2).
| REC land-use class | Direction of change | Number of sites | All years | Decades | ||||||
|---|---|---|---|---|---|---|---|---|---|---|
| F-value | P-value | F-value | P-value | |||||||
| All sites | negative | 22545 | 191.2 | 0.0001 | 223.7 | 0.0001 | ||||
| Pasture | negative | 9931 | 92.0 | 0.0001 | 118.4 | 0.0001 | ||||
| Tussock | negative | 2805 | 21.1 | 0.0001 | 38.83 | 0.0001 | ||||
| Indigenous | positive | 5529 | 41.5 | 0.0001 | 24.7 | 0.0001 | ||||
| Urban | negative | 1157 | 29.6 | 0.0001 | 19.9 | 0.001 | ||||
| Scrub | negative/ns | 1193 | 3.9 | 0.047 | 1.21 | 0.27 | ||||
| Exotic | ns | 1318 | 2.4 | 0.13 | 0.09 | 0.77 | ||||

Graph showing average Index of Biotic Integrity (IBI) score for all river sites monitored between 1970 and 2007 by decade, with standard error (SE).
| Decade | Number of sites | Average IBI score ± SE |
|---|---|---|
| 1970-1979 | 1142 | 32.86 ± 0.47 |
| 1980-1989 | 4598 | 32.11 ± 0.26 |
| 1990-1999 | 7415 | 31.64 ± 0.22 |
| 2000-2007 | 9390 | 27.91 ± 0.21 |

Graph showing average Index of Biotic Integrity (IBI) score for river sites with indigenous forest land-cover type monitored between 1970 and 2007 by decade, with standard error (SE). Land-cover type defined using River Environment Classification
| Decade | Number of sites | Average IBI score ± SE |
|---|---|---|
| 1970-1979 | 451 | 30.40 ± 0.81 |
| 1980-1989 | 1584 | 34.27 ± 0.49 |
| 1990-1999 | 2097 | 36.61 ± 0.44 |
| 2000-2007 | 1401 | 35.94 ± 0.51 |

Graph showing average Index of Biotic Integrity (IBI) score for river sites with pasture land-cover type monitored between 1970 and 2007 by decade, with standard error (SE). Land-cover type defined using River Environment Classification
| Decade | Number of sites | Average IBI score ± SE |
|---|---|---|
| 1970-1979 | 384 | 33.57 ± 0.76 |
| 1980-1989 | 2087 | 31.79 ± 0.34 |
| 1990-1999 | 3075 | 31.28 ± 0.30 |
| 2000-2007 | 4390 | 27.63 ± 0.28 |

Graph showing average Index of Biotic Integrity (IBI) score for river sites with scrub land-cover type monitored between 1970 and 2007 by decade, with standard error (SE). Land-cover type defined using River Environment Classification
| Decade | Number of sites | Average IBI score ± SE |
|---|---|---|
| 1970-1979 | 45 | 41.23 ± 1.98 |
| 1980-1989 | 146 | 38.92 ± 1.53 |
| 1990-1999 | 525 | 34.17 ± 0.91 |
| 2000-2007 | 482 | 36.58 ± 0.90 |

Graph showing average Index of Biotic Integrity (IBI) score for river sites with urban land-cover type monitored between 1970 and 2007 by decade, with standard error (SE). Land-cover type defined using River Environment Classification
| Decade | Number of sites | Average IBI score ± SE |
|---|---|---|
| 1970-1979 | 26 | 33.36 ± 3.24 |
| 1980-1989 | 106 | 33.96 ± 1.34 |
| 1990-1999 | 306 | 29.67 ± 0.89 |
| 2000-2007 | 733 | 26.74 ± 0.69 |

Graph showing average Index of Biotic Integrity (IBI) score for river sites with exotic forest land-cover type monitored between 1970 and 2007 by decade, with standard error (SE). Land-cover type defined using River Environment Classification
| Decade | Number of sites | Average IBI score ± SE |
|---|---|---|
| 1970-1979 | 55 | 35.41 ± 2.23 |
| 1980-1989 | 186 | 33.98 ± 1.30 |
| 1990-1999 | 491 | 26.62 ± 0.86 |
| 2000-2007 | 591 | 32.20 ± 0.83 |

Graph showing average Index of Biotic Integrity (IBI) score for river sites with tussock land-cover type monitored between 1970 and 2007 by decade, with standard error (SE). Land-cover type defined using River Environment Classification
| Decade | Number of sites | Average IBI score ± SE |
|---|---|---|
| 1970-1979 | 93 | 34.15 ± 1.57 |
| 1980-1989 | 429 | 21.72 ± 0.92 |
| 1990-1999 | 734 | 15.03 ± 0.69 |
| 2000-2007 | 1550 | 17.37 ± 0.52 |
A number of reaches in the database had been sampled more than once over the 37 years and these were analysed to look for changes. This could be done using the number of species as they were directly comparable, which negated the need to use the IBI. The number of sites that had more or fewer species in 2000 – 2007 compared to the 5-year periods between 1970 and 1999 were counted. There were consistently more sites sampled on the same reach that had less species than more species in the 2000s, when compared to the 5-year periods between 1975 and 1999 (Fig. 11). The only exception was the 1970 – 1975 period when more sites gained species than lost species over time. However, there were only 40 sites that were sampled in 1970 – 1975 that changed from the 2000s.

Graph showing the number of sites with fewer or more fish species found prior to 2000 compared with when the same stream reach was sampled between 2000 and 2007.
| Time period | Sites with fewer species in 2000-07 | Sites with more species in 2000-07 |
|---|---|---|
| 1970-1974 | 17 | 22 |
| 1975-1979 | 56 | 45 |
| 1980-1984 | 35 | 26 |
| 1985-1989 | 82 | 50 |
| 1990-1994 | 162 | 161 |
| 1995-1995 | 101 | 86 |
Using rank correlation analysis there were very few significant differences found in IBI scores over time for the stream reaches that had multiple sampling events. However, there were more significant positive than negative changes for native vegetation, no differences at exotic forest, but more negative correlations at pasture and urban sites (Fig. 12).

Graph showing significant changes in IBI scores over time for stream reaches with multiple sampling events by land cover type
| Land cover type | Number of sites with negative trends | Number of sites with positive trends |
|---|---|---|
| Pasture | 15 | 11 |
| Native | 4 | 10 |
| Exotic | 2 | 2 |
| Urban | 3 | 2 |