Risks are ranked from 1 (lowest risk) to 10 (highest risk).
The risk scores were assigned based on expert opinion which was supported by a review of the literature on known mechanisms of microbial transport and a review of all modelling studies published in NZ. We thus use multiple lines of evidence to support the matrix.
Risk scores are applied on the basis of land cover, soils, altitude and climate.
The aim of this work was to develop a simple tool to help a regional council make decisions at a freshwater management unit (FMU) scale.
Note: The risk score map is made with “current landuse” which relates to 2018 year data. But users should be able to change land use and see the effect on the risk matrix.
Date: August 2023 Version: v1
Contact: Richard Muirhead, Sandy Elliott, Ton Snelder, Annette Semadeni-Davies
100m (1ha); coarsest datasets are climate normals (500m)
Evaluation method (Validation)
Evaluation result (Numeric)
Evaluation result (Categorical)
Uncertainty data format (Numeric)
Uncertainty data format (Categorical)
The risk matrix mapping was done using a bespoke Python script that reads in pre-processed rasterised (1 ha resolution) spatial data, mostly sourced from publicly available special datasets, classifies that data into the classes used in the matrix, determines the risk score for each 1ha cell according to the matrix using raster calculations and writes out a georeferenced tiff file of the risk scores. The datasets and classes are listed below, the data sources are in the data links table.
Land use was was derived from LCDB5 and Agribase and relates to the year 2018 . There are five broad land use classes in the matrix, pasture, arable land, horticulture, urban land use and low emission land uses that includes native and exotic forestry, tussock, bare rock and ungrazed pasture. Elevation
There are two elevation classes in the matrix, lowland (<= 350 m.a.s.l) and hill to high country (>= 350 m.a.s.l). The layer was derived from the 15 m Digital Elevation Model (DEM) grid created by the University of Otago upscaled to 1ha resolution.
Wetness is an index derived from three data sets, irrigated land, rainfall and potential evapotranspiration (PET).
Irrigated land was taken from the MfE 2017 irrigation shapefile, this is a Boolean dataset where a value of 1 represents irrigated land and a value of 0 no irrigation. The shapefile was rasterised for mapping.
Mean annual Rainfall and PET were taken from the NIWA’s climate normal grids (1991-2020). The unit for both datasets is mm depth. The grids have a 500m resolution and were downscaled to 1 ha.
The wetness index is assigned as follows:
all irrigated land
non-irrigated land with a mean annual rainfall <= 1700 m
non-irrigated land with a mean annual rainfall < 1700 m where the effective rainfall (rain – PET) < 1100 mm
non-irrigated land with a mean annual rainfall < 1700 m where the effective rainfall <= 1100 mm