2.7.5.4 Receptor impact models


2.7.5.4.1 Pilliga riverine (upland and lowland)

Table 39 summarises the elicitation design matrix for the projected foliage cover of riparian trees in the riverine landscape classes in the Pilliga region. The first three design points – design point identifiers 1 to 3 – address the predicted variability (across the streams in the landscape class during the reference interval) in the zero-flow regime. Note that the design point identifiers are simply index variables that identify the row of the elicitation design matrix. They are included here to maintain an auditable path between analysis and reporting.

The first design points provide for an estimate of the uncertainty in mean projected foliage cover across the landscape class in the reference year 2012 (Yref). The remaining design points represent hydrological scenarios that span the uncertainty in the values of the hydrological response variables in the relevant time period of hydrological history associated with the short- (2042) and long- (2102) assessment years.

Table 39 Elicitation design matrix for the annual mean projected foliage cover of riparian trees along Pilliga riverine landscape classes that includes: yellow box, white cypress pine, Eucalyptus crebra, dirty gum, Blakely's red gum, Angophora floribunda, Eucalyptus fibrosa, fuzzy box in the Namoi subregion zone of potential hydrological change


Id

ZQD

dmaxRef

Yref

Year

tmaxRef

3

176.67

na

na

2012

0

1

0.00

na

na

2012

0

2

88.34

na

na

2012

0

26

83.40

56.80

0.12

2042

2103

28

0.00

0.20

0.20

2042

2027

48

178.00

0.20

0.20

2042

2103

15

178.00

8.16

0.12

2042

2065

36

178.00

56.80

0.20

2042

2027

52

0.00

56.80

0.20

2042

2103

86

83.40

8.16

0.20

2102

2027

61

0.00

56.80

0.12

2102

2027

73

0.00

0.20

0.12

2102

2103

57

178.00

0.20

0.12

2102

2027

92

83.40

0.20

0.20

2102

2065

108

178.00

56.80

0.20

2102

2103

Receptor impact modelling elicitation design matrix for annual mean projected foliage cover in a transect of 50 m length and 20 m width that extends from first bench (‘toe’) on both sides of the stream in permanent and temporary upland streams. Design points for Yref in the future (short- and long-assessment periods) are calculated during the receptor impact modelling elicitation workshop using elicited values for the receptor impact variable in the reference period. All other design points (with identifiers Id) are either default values or values determined by groundwater and surface water modelling. Hydrological response variables are as described in Table 37. na = not applicable

Data: Bioregional Assessment Programme (Dataset 4)

Design point identifiers 15 through to 108 (as listed in Table 39) represent combinations of the three hydrological response variables (dmaxRef, tmaxRef and ZQD), together with high and low values of Yref (see companion submethodology M08 (as listed in Table 1) for receptor impact modelling (Hosack et al., 2018)). The high and low values for Yref were calculated during the receptor impact modelling workshop following the experts’ response to the first design point, and then automatically included within the design for the elicitations at the subsequent design points.

The receptor impact modelling methodology allows for a very flexible class of statistical models to be fitted to the values of the receptor impact variables elicited from the experts at each of the design points (companion submethodology M08 (as listed in Table 1) for receptor impact modelling (Hosack et al., 2018)). The model fitted to the elicited values of mean foliage cover for the Pilliga riverine landscape classes is summarised in Figure 29 and Table 40. The fitted model takes the form:

eta equals h open parenthesis y close parenthesis equals beta subscript 0 end subscript x subscript 0 end subscript plus beta subscript f end subscript x subscript f end subscript plus beta subscript l end subscript x subscript l end subscript plus beta subscript r end subscript x subscript r end subscript plus sum from j equals 1 to 3 of beta subscript h subscript j end subscript end subscript x subscript h subscript j end subscript end subscript

(7)

where x subscript 0 end subscript is an intercept term (a vector of ones), x subscript f end subscript is a binary indicator variable scored 1 for the case of an assessment in the short- or long-assessment year, x subscript l end subscript is a binary indicator variable scored 1 for the case of an assessment in the long-assessment year, x subscript r end subscript is a continuous variable that represent the value of the receptor impact variable in the reference year (Yref, set to zero for the case of an assessment in the reference year), and x subscript h subscript j end subscript end subscript comma j equals 1 dot dot dot 4 are the (continuous or integer) values of the three hydrological response variables (dmaxRef, tmaxRef and ZQD). Note that the modelling framework provides for more complex models, including quadratic value of, and in interactions between, the hydrological response variables but in this instance the simple linear model (Equation 7) was identified as the most parsimonious representation of the experts’ responses.

The model estimation procedure adopts a Bayesian approach. The model coefficients (beta subscript 0 end subscript comma beta subscript f end subscript comma beta subscript l end subscript comma beta subscript r end subscript comma beta subscript h subscript j end subscript end subscript are assumed to follow a multivariate normal distribution. The Bayesian estimation procedure quantifies how compatible different values of the parameters of this distribution are with the data (the elicited expert opinion) under the model. The (marginal) mean and 80% central credible intervals[3] of the three hydrological response variable coefficients are summarised in partial regression plots in Figure 29, while Table 40 summarises the same information for all seven model coefficients.

The model indicates that the experts’ opinion provides strong evidence for Yref having a positive effect on average percent projected foliage cover. This suggests that given a set of hydrological response variable values in the future, a site with a higher foliage cover at the 2012 reference point is more likely to have a higher foliage cover in the future than a site with a lower foliage cover value at this time point. This reflects the lag in the response of foliage cover to changes in hydrological response variables that would be expected of mature trees with long life spans.

The model also indicates that the experts’ opinion provides some evidence for dmaxRef having a negative effect on average percent projected foliage cover (Figure 29). This suggests that percent projected foliage cover will decrease as groundwater drawdown increases due to coal resource development. The model predicts that (holding all other hydrological response variables constant at the mid-point of their elicitation range) the mean of the average percent projected foliage cover will drop from just under 12% without any change in groundwater level, to about 10% if the levels decrease by 50 m relative to the reference level in 2012. There is, however, considerable uncertainty in these predictions, with an 80% chance that the foliage cover will lie somewhere between approximately 8% and 16% on the short-assessment period, and somewhere between roughly 7% and 17% in the long-assessment period, with a 6 m drop in groundwater level. In relation to dmaxRef, Yrs2tmax is also found slightly significant, but with a positive effect. This suggests that a maximum drawdown happening very early would have a larger effect on foliage cover than a maximum drawdown reached after a long time.

The model indicates that the experts’ opinion provides strong evidence for ZQD (zero-flow days (averaged over 30 years) subsequently referred to in this Section as ‘zero-flow days’) having a negative effect on average percent projected foliage cover. The model predicts that (holding all other hydrological response variables constant at the mid-point of their elicitation range) the mean of the average percent projected foliage cover will decrease from just under 25% without any zero-flow days to about 10% if the number of zero-flow days increases to 180.

Finally, the model also indicates some diverging influence between short-term and long-term influence (holding hydrological response variables constant). For the short-assessment period, experts believe in a relative increase of foliage cover, while they are left uncertain about the long‑term assessment effect. An interpretation is that the effects of changes in hydrology are not immediate on foliage cover, and 2102 is a very long time into the future to make assessments without uncertainty.

Figure 29

Figure 29 (Top row) Predicted mean (black dot) and 80% central credible interval (grey polygon or line) of annual mean projected foliage cover along Pilliga riverine landscape classes under reference hydrological conditions. (Middle and bottom rows) Predicted future effect (mean = black line, 80% central credible interval = grey polygon) of each hydrological response variable on annual mean projected foliage cover along Pilliga riverine landscape classes, holding all other hydrological response variables constant at the midpoint of their elicitation range (during risk estimation all hydrological response variables vary simultaneously)

Dashed vertical lines show hydrological response variable range used in the elicitation. Hydrological response variables ZQD and dmaxRef are as defined in Table 37. Yrs2tmaxRef is the difference between tmaxRef and the assessment year that is relevant for the prediction (2012, 2042 or 2102). The numbers on the y-axis range from 0 to 1 as the receptor impact model was constructed using the proportion for the statistical modelling. They should be interpreted as a percent foliage cover ranging from 0 to 100%.

Data: Bioregional Assessment Programme (Dataset 4)

Table 40 Mean, 10th and 90th percentile of the coefficients of the receptor impact model for annual mean projected foliage cover along Pilliga riverine landscape classes


Mean

q10

q90

(Intercept)

–1.4

–1.69

–1.11

future1

1.47

0.335

2.6

long1

0.161

–0.332

0.653

Yref

0.938

0.238

1.64

dmaxRef

–0.00493

–0.0116

0.00171

Yrs2tmaxRef

0.0018

–0.00306

0.00666

ZQD

–0.00477

–0.00678

–0.00276

Future1 is the indicator variable for the short-assessment year (2042). Long1 is the indicator variable for the long-assessment year (2102). Yref quantifies the value of the receptor impact variable during the reference period. Hydrological response variables dmaxRef and ZQD are as defined in Table 37. Yrs2tmaxRef is the difference between tmaxRef and the assessment year that is relevant for the prediction (2012, 2042 or 2102).

Data: Bioregional Assessment Programme (Dataset 4)

Table 41 summarises the elicitation matrix for the average number of families of aquatic macroinvertebrates in instream pool habitats. The first six design points – design points 1 to 7 as shown – address the predicted variability (across the landscape class in the reference interval) in ZQD and ZME, capturing the lowest and highest predicted values together with two intermediate values. These design points provide for an estimate of the uncertainty in aquatic macroinvertebrate family abundance across the landscape class in the reference year 2012 (Yref).

Design points 21 to 213 inclusive (as listed in Table 41) represent scenarios that span the uncertainty in the predicted values of ZQD and ZME in the relevant time period of hydrological history associated with the short- (2042) and long- (2102) assessment years, combined with high and low values of Yref, as well as groundwater drawdown and time to drawdown combinations. Again, the high and low values for Yref were calculated during the receptor impact modelling workshop.

The fitted model for number of families of macroinvertebrates in instream pool habitats takes the form:

eta equals h open parenthesis y close parenthesis equals beta subscript 0 end subscript x subscript 0 end subscript plus beta subscript f end subscript x subscript f end subscript plus beta subscript l end subscript x subscript l end subscript plus beta subscript r end subscript x subscript r end subscript plus sum from j equals 1 to 4 of beta subscript h subscript j end subscript end subscript x subscript h subscript j end subscript end subscript

(8)

where x subscript 0 end subscript is an intercept term (a vector of ones), x subscript f end subscript is a binary indicator variable scored 1 for the case of an assessment in the short- or long-assessment year, x subscript l end subscript is a binary indicator variable scored 1 for the case of an assessment in the long-assessment year, x subscript r end subscript is a continuous variable that represent the value of the receptor impact variable in the reference year (Yref, set to zero for the case of an assessment in the reference year), and x subscript h subscript j end subscript end subscript comma j equals 1 dot dot dot 4 are the (continuous or integer) values of the four hydrological response variables (dmax, tmax, ZQD and ZME). Note that the modelling framework provides for more complex models, including quadratic value of, and in interactions between, the hydrological response variables but in this instance the simple linear model (Equation 8) was identified as the most parsimonious representation of the experts’ responses.

Table 41 Elicitation design matrix of the average number of families of aquatic macroinvertebrates in instream pool habitat along Pilliga riverine landscape classes in Namoi subregion zone of potential hydrological change


Id

ZQD

ZME

dmaxRef

Yref

Year

tmaxRef

2

88.33

0.32

na

na

2012

0

4

88.34

39.87

na

na

2012

0

5

176.66

0.94

na

na

2012

0

1

0.00

0.00

na

na

2012

0

3

88.34

19.93

na

na

2012

0

7

176.67

39.87

na

na

2012

0

59

177.99

0.95

56.80

13

2042

2103

29

0.00

0.00

8.16

13

2042

2065

70

178.00

39.60

0.20

19

2042

2027

87

177.99

0.95

0.20

19

2042

2065

42

178.00

39.60

56.80

13

2042

2065

41

83.40

19.30

0.20

13

2042

2027

111

83.40

39.60

0.20

19

2042

2103

67

83.40

19.30

0.20

19

2042

2027

21

83.40

0.30

0.20

13

2042

2027

213

177.99

0.95

0.20

19

2102

2065

132

83.40

39.60

0.20

13

2102

2027

210

178.00

39.60

56.80

19

2102

2027

139

83.40

39.60

8.16

13

2102

2027

212

83.40

0.30

0.20

19

2102

2065

190

0.00

0.00

0.20

19

2102

2027

170

83.40

0.30

0.20

13

2102

2103

131

178.00

19.30

0.20

13

2102

2027

Receptor impact modelling elicitation design matrix for average number of families of aquatic macroinvertebrates in instream pool habitat in permanent and temporary upland and lowland streams (Pilliga riverine landscape classes). Design points for Yref in the future (short- and long-assessment periods) are calculated during the receptor impact modelling elicitation workshop using elicited values for the receptor impact variable in the reference period. All other design points (with identifiers Id) are either default values or values determined by groundwater and surface water modelling. Hydrological response variables are as described in Table 37.

na = not applicable

Data: Bioregional Assessment Programme (Dataset 4)

The hydrological response variable in the macroinvertebrate model varies during the reference interval and the future interval. The model indicates that the experts’ elicited information supports the hypothesis that an increase in ZQD and/or ZME will have a slightly negative effect on the number of families of macroinvertebrates despite the experts being somewhat uncertain about its average value. The model suggests that it can vary across the landscape class from less than 15 to almost 20 under conditions of constant flow (ZQD = zero), holding all other covariates at their mid-values. As the number of zero-flow days increases, however, experts were of the opinion that the number of families would drop quite dramatically with values between 13 and 6 (under very intermittent flow conditions ZQD >150 days) (Figure 30).

There was very little evidence in the elicited data to suggest that this effect would be substantially different in the future assessment years. Again, this is indicated by the almost identical partial regression plots in the reference, short- and long-assessment years (Figure 30), and the relatively large negative 10th and positive 90th percentiles for the long and future coefficients in Table 42. The model does, however, suggest that the experts’ uncertainty increased a lot for predictions in the future assessment years relative to the reference year.

The best-fitting model is unable to eliminate the possibility that the average number of families of macroinvertebrates in the reference years has no influence on its number in the future years. This is indicated by the fact the model automatically dropped this variable from the model. This suggestion is consistent with the hypothesis that there is likely to be very little lag in the response of this short-lived species to changes in the hydrological response variables.

The model also indicates that the experts’ opinion provides strong evidence for dmaxRef having a negative effect on macroinvertebrates in instream pool habitats. This suggests that the number of families of macroinvertebrates will decrease as groundwater drawdown increases due to coal resource development. The model predicts that (holding all other hydrological response variables constant at the mid-point of their elicitation range) the number of families will drop from just under 12 without any change in groundwater level, to about 6 if the levels decrease by 55 m relative to the reference level in 2012. There is, however, considerable uncertainty in these predictions, with an 80% chance that the number of families will lie somewhere between approximately 1 and 12 on the short-assessment period and in the long-assessment period. In relation to the time of drawdown, Yrs2tmax is not found significant, meaning that the time at which the maximum drawdown will occur does not carry much influence over the number of families of aquatic macroinvertebrates.

Figure 30

Figure 30 (Top row) Predicted mean (black dot) and 80% central credible interval (grey polygon) of average number of families of aquatic macroinvertebrates in instream pool habitat in Pilliga riverine landscape classes under reference hydrological conditions. (Middle and bottom rows) Predicted future effect (mean = black line, 80% central credible interval = grey polygon) of each hydrological response variable on average number of families of aquatic macroinvertebrates in instream pool habitat in Pilliga riverine landscape classes, holding all other hydrological response variables constant at the midpoint of their elicitation range (during risk estimation all hydrological response variables vary simultaneously)

Dashed vertical lines show hydrological response variable range used in the elicitation. Hydrological response variables are as defined in Table 37. Yrs2tmaxRef is the difference between tmaxRef and the assessment year that is relevant for the prediction (2012, 2042 or 2102).

Data: Bioregional Assessment Programme (Dataset 4)

Table 42 Mean, 10th and 90th percentile of the coefficients of the receptor impact model for average number of families of aquatic macroinvertebrates in instream pool habitat in Pilliga riverine landscape classes


Mean

q10

q90

(Intercept)

2.87

2.66

3.08

future1

–0.0138

–0.198

0.17

long1

0.0151

–0.268

0.298

dmaxRef

–0.0163

–0.0285

–0.0041

Yrs2tmaxRef

0.000373

–0.00274

0.00349

ZQD

–0.00379

–0.0084

0.000817

ZME

–0.000412

–0.00186

0.00101

Future1 is the indicator variable for the short-assessment year (2042). Long1 is the indicator variable for the long-assessment year (2102). Hydrological response variables are as defined in Table 37. Yrs2tmaxRef is the difference between tmaxRef and the assessment year that is relevant for the prediction (2012, 2042 or 2102).

Data: Bioregional Assessment Programme (Dataset 4)

Last updated:
6 December 2018
Thumbnail of the Namoi subregion

Product Finalisation date

2018
PRODUCT CONTENTS

ASSESSMENT