In the November 2023 release of P2R Projector two new project types (Gully and Streambank) were developed. These project types are adapted from the Gully and Streambank Toolbox methods, developed for Reef Trust by CSIRO, The Australian National University, Griffith University and the Queensland Department of Agriculture and Fisheries, for calculating sediment savings from erosion control activities. For more detailed information on the application of these methods, see the Gully and Streambank Toolbox Technical Guide.

Gully and Streambank project explanatory notes on this page:

Gully erosion

Gully and Streambank Toolbox automatic data retrieval

November 2023 release

To support the use of P2R Projector Gully and Streambank Toolbox for ‘prospecting’ projects, some required data values are retrieved from publicly available API sources to i) prevent calculations from failing and ii) provide initial erosion estimates that are representative of expected values. These surrogate data values are provided as a convenience and can be overridden with locally derived data in all cases when a greater focus on accuracy and comparability are required.

Climate Correction Metrics

Metrics of long term and recent period daily rainfall averages are source from the SILO Data Drill API, using a point (latitude, longitude) that is calculated as the geometric centre of the geometric feature. The Climate Correction Factor is the Long Term Daily Rainfall Average divided by the Recent Period Daily Rainfall Average.

Soil Order

A representative Australian Soil Order is estimated using the Modelled Australian Soil Classification Map  WMS Service.

Trials querying the WMS by polygon and polygon and polyline features revealed that the ArcGIS REST Service ultimately performed a point-based query using a geometric centre. Until that behaviour is changed, the Toolbox API call will use a simple point (latitude, longitude) that is calculated as the geometric centre of the geometric feature. The Soil Order returned is mapped from the numeric according to its numeric pixel value.

Numeric pixel valueSoil order
1Vertosol
2Sodosol
3Dermosol
4Chromosol
5Ferrosol
6Kurosol
7Tenosol
8Kandosol
9Hydrosol
10Podosol
11Rudosol
12Calcarosol
13Organosol
14Anthroposol

Soil Parameters

Soil parameter values for Silt, Clay, Bulk Density, Total Nitrogen and Total Phosphorus are retrieved from the data layers made available as the Soil and Landscape Grid Of Australia. The layers are queried using the geometric centre of the geometric feature, calling the ‘extractSLGAdata’ service provided by the API. Data values from the soil layer with upper depth 30cm and lower depth 60cm are used to provide surrogate data values.

Streambank ‘Catchment modelled method’ computation

November 2023 release

To support the ‘Catchment modelled method’ of determining streambank erosion rate, parameters and outputs from the Great Barrier Reef (GBR) Dynamic SedNet models (Report Card 2019, RC10) were harvested and processed.

The primary data element to be calculated for Toolbox use is an estimate of an annual average mass of fine sediment (kg) eroded per unit of representative stream length (km) in each GBR Dynamic SedNet subcatchment. Using the RC10 prediction of annual average streambank erosion mass (fine sediment, kg per year) it is possible to directly calculate the required streambank erosion rate using the model’s internally stored stream length parameter. However, it is acknowledged that the density of the stream network represented by the GBR Dynamic SedNet models is much sparser than the existing stream network providing streambank sediments throughout much of the GBR catchments. A search for a more representative stream network data set led to the identification of the HydroRIVERS layer (part of the HydroSHEDS collection).

The HydroRIVERS data set was intersected with the modelling subcatchments from the GBR Dynamic SedNet models (RC10), with the length of intersected river segments calculated. On average, the length of intersected HydroRIVERS stream segments is 3.3 times greater than the relevant GBR Dynamic SedNet stream length.

Where the GBR Dynamic SedNet fine sediment erosion mass was greater than zero, an initial fine sediment erosion rate (kg/km) was calculated using the derived HydroRIVERS stream lengths (for convenience, erosion rates for particulate nitrogen and particulate phosphorus were calculated in the same manner). Recognising that the custodians of P2R Projector wanted no instances of zero streambank erosion represented (preferring a low rate represented) in the Toolbox, each of the 6 GBR NRM Regions (Cape York, Wet Tropics, Burdekin, Mackay Whitsundays, Fitzroy, Burnett Mary) had the initial fine sediment rates analysed for a suitable minimum value to be applied.

Using Microsoft Excel, the fine sediment erosion representing the 15th Percentile was selected as the minimum value to apply in low and zero erosion modelling subcatchments for each NRM Region. Similarly, to provide an upper cap of fine sediment erosion rates in P2R Projector Toolbox, the value representing the 97th Percentile was selected. It is these fine sediment erosion rates (kg/km, constrained by the 15th and 97th Percentile of regional rates) that support the P2R Projector implementation of the Toolbox Catchment Modelling Method. The same process and percentile constraints have been applied to the erosion rates of particulate nitrogen and particulate phosphorus.

To compliment the infilled streambank erosion rates, any GBR Dynamic SedNet subcatchments with a zero or missing stream length have had a default of 1 km provided for Toolbox use.