H Dryopteris dilatata (Hoffmann) A. Gray - Erysimum cheiranthoides L.

H.1 Dryopteris dilatata (Hoffmann) A. Gray

Table H.1: Summary table of data for Dryopteris dilatata (Hoffmann) A. Gray.

Figure H.1: Effect of year on the probability of Dryopteris dilatata (Hoffmann) A. Gray presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.2: The same as H.1, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.3: Visualisation of the spatial smooth effect on the probability of Dryopteris dilatata (Hoffmann) A. Gray presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.2 Dryopteris filix-mas (L.) Schott

Table H.2: Summary table of data for Dryopteris filix-mas (L.) Schott.

Figure H.4: Effect of year on the probability of Dryopteris filix-mas (L.) Schott presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.5: The same as H.4, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.6: Visualisation of the spatial smooth effect on the probability of Dryopteris filix-mas (L.) Schott presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.3 Duchesnea indica (Andrews) Focke

Table H.3: Summary table of data for Duchesnea indica (Andrews) Focke.

Figure H.7: Effect of year on the probability of Duchesnea indica (Andrews) Focke presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.8: The same as H.7, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.9: Visualisation of the spatial smooth effect on the probability of Duchesnea indica (Andrews) Focke presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.4 Echinochloa crus-galli (L.) Beauv.

Table H.4: Summary table of data for Echinochloa crus-galli (L.) Beauv..

Figure H.10: Effect of year on the probability of Echinochloa crus-galli (L.) Beauv. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.11: The same as H.10, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.12: Visualisation of the spatial smooth effect on the probability of Echinochloa crus-galli (L.) Beauv. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.5 Echium vulgare L.

Table H.5: Summary table of data for Echium vulgare L..

Figure H.13: Effect of year on the probability of Echium vulgare L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.14: The same as H.13, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.15: Visualisation of the spatial smooth effect on the probability of Echium vulgare L. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.6 Eleocharis multicaulis (Smith) Desv.

Table H.6: Summary table of data for Eleocharis multicaulis (Smith) Desv..

Figure H.16: Effect of year on the probability of Eleocharis multicaulis (Smith) Desv. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.17: The same as H.16, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.18: Visualisation of the spatial smooth effect on the probability of Eleocharis multicaulis (Smith) Desv. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.7 Eleocharis palustris (L.) Roem. et Schult.

Table H.7: Summary table of data for Eleocharis palustris (L.) Roem. et Schult..

Figure H.19: Effect of year on the probability of Eleocharis palustris (L.) Roem. et Schult. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.20: The same as H.19, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.21: Visualisation of the spatial smooth effect on the probability of Eleocharis palustris (L.) Roem. et Schult. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.8 Elodea canadensis Michaux

Table H.8: Summary table of data for Elodea canadensis Michaux.

Figure H.22: Effect of year on the probability of Elodea canadensis Michaux presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.23: The same as H.22, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.24: Visualisation of the spatial smooth effect on the probability of Elodea canadensis Michaux presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.9 Elodea nuttallii (Planch.) St John

Table H.9: Summary table of data for Elodea nuttallii (Planch.) St John.

Figure H.25: Effect of year on the probability of Elodea nuttallii (Planch.) St John presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.26: The same as H.25, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.27: Visualisation of the spatial smooth effect on the probability of Elodea nuttallii (Planch.) St John presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.10 Elymus athericus (Link) Kerguélen

Table H.10: Summary table of data for Elymus athericus (Link) Kerguélen.

Figure H.28: Effect of year on the probability of Elymus athericus (Link) Kerguélen presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.29: The same as H.28, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.30: Visualisation of the spatial smooth effect on the probability of Elymus athericus (Link) Kerguélen presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.11 Elymus repens (L.) Gould

Table H.11: Summary table of data for Elymus repens (L.) Gould.

Figure H.31: Effect of year on the probability of Elymus repens (L.) Gould presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.32: The same as H.31, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.33: Visualisation of the spatial smooth effect on the probability of Elymus repens (L.) Gould presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.12 Epilobium ciliatum Rafin.

Table H.12: Summary table of data for Epilobium ciliatum Rafin..

Figure H.34: Effect of year on the probability of Epilobium ciliatum Rafin. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.35: The same as H.34, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.36: Visualisation of the spatial smooth effect on the probability of Epilobium ciliatum Rafin. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.13 Epilobium hirsutum L.

Table H.13: Summary table of data for Epilobium hirsutum L..

Figure H.37: Effect of year on the probability of Epilobium hirsutum L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.38: The same as H.37, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.39: Visualisation of the spatial smooth effect on the probability of Epilobium hirsutum L. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.14 Epilobium montanum L.

Table H.14: Summary table of data for Epilobium montanum L..

Figure H.40: Effect of year on the probability of Epilobium montanum L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.41: The same as H.40, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.42: Visualisation of the spatial smooth effect on the probability of Epilobium montanum L. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.15 Epilobium palustre L.

Table H.15: Summary table of data for Epilobium palustre L..

Figure H.43: Effect of year on the probability of Epilobium palustre L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.44: The same as H.43, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

## Error in nullGrob(): could not find function "nullGrob"

H.16 Epilobium parviflorum Schreb.

Table H.16: Summary table of data for Epilobium parviflorum Schreb..

Figure H.45: Effect of year on the probability of Epilobium parviflorum Schreb. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.46: The same as H.45, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.47: Visualisation of the spatial smooth effect on the probability of Epilobium parviflorum Schreb. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.17 Epilobium roseum Schreb.

Table H.17: Summary table of data for Epilobium roseum Schreb..

Figure H.48: Effect of year on the probability of Epilobium roseum Schreb. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.49: The same as H.48, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.50: Visualisation of the spatial smooth effect on the probability of Epilobium roseum Schreb. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.18 Epilobium tetragonum L.

Table H.18: Summary table of data for Epilobium tetragonum L..

Figure H.51: Effect of year on the probability of Epilobium tetragonum L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.52: The same as H.51, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.53: Visualisation of the spatial smooth effect on the probability of Epilobium tetragonum L. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.19 Epipactis helleborine (L.) Crantz

Table H.19: Summary table of data for Epipactis helleborine (L.) Crantz.

Figure H.54: Effect of year on the probability of Epipactis helleborine (L.) Crantz presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.55: The same as H.54, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.56: Visualisation of the spatial smooth effect on the probability of Epipactis helleborine (L.) Crantz presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.20 Equisetum arvense L.

Table H.20: Summary table of data for Equisetum arvense L..

Figure H.57: Effect of year on the probability of Equisetum arvense L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.58: The same as H.57, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.59: Visualisation of the spatial smooth effect on the probability of Equisetum arvense L. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.21 Equisetum fluviatile L.

Table H.21: Summary table of data for Equisetum fluviatile L..

Figure H.60: Effect of year on the probability of Equisetum fluviatile L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.61: The same as H.60, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.62: Visualisation of the spatial smooth effect on the probability of Equisetum fluviatile L. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.22 Equisetum palustre L.

Table H.22: Summary table of data for Equisetum palustre L..

Figure H.63: Effect of year on the probability of Equisetum palustre L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.64: The same as H.63, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.65: Visualisation of the spatial smooth effect on the probability of Equisetum palustre L. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.23 Equisetum telmateia Ehrh.

Table H.23: Summary table of data for Equisetum telmateia Ehrh..

Figure H.66: Effect of year on the probability of Equisetum telmateia Ehrh. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.67: The same as H.66, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.68: Visualisation of the spatial smooth effect on the probability of Equisetum telmateia Ehrh. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.24 Eragrostis minor Host

Table H.24: Summary table of data for Eragrostis minor Host.

Figure H.69: Effect of year on the probability of Eragrostis minor Host presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.70: The same as H.69, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.71: Visualisation of the spatial smooth effect on the probability of Eragrostis minor Host presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.25 Eragrostis pilosa (L.) Beauv.

Table H.25: Summary table of data for Eragrostis pilosa (L.) Beauv..

Figure H.72: Effect of year on the probability of Eragrostis pilosa (L.) Beauv. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.73: The same as H.72, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.74: Visualisation of the spatial smooth effect on the probability of Eragrostis pilosa (L.) Beauv. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.26 Erica tetralix L.

Table H.26: Summary table of data for Erica tetralix L..

Figure H.75: Effect of year on the probability of Erica tetralix L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.76: The same as H.75, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.77: Visualisation of the spatial smooth effect on the probability of Erica tetralix L. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.27 Erigeron annuus (L.) Desf.

Table H.27: Summary table of data for Erigeron annuus (L.) Desf..

Figure H.78: Effect of year on the probability of Erigeron annuus (L.) Desf. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.79: The same as H.78, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.80: Visualisation of the spatial smooth effect on the probability of Erigeron annuus (L.) Desf. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.28 Conyza canadensis (L.) Cronq.

Table H.28: Summary table of data for Conyza canadensis (L.) Cronq..

Figure H.81: Effect of year on the probability of Conyza canadensis (L.) Cronq. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.82: The same as H.81, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.83: Visualisation of the spatial smooth effect on the probability of Conyza canadensis (L.) Cronq. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.29 Eriophorum polystachion L.

Table H.29: Summary table of data for Eriophorum polystachion L..

Figure H.84: Effect of year on the probability of Eriophorum polystachion L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.85: The same as H.84, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.86: Visualisation of the spatial smooth effect on the probability of Eriophorum polystachion L. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.30 Erodium cicutarium (L.) L’Hérit.

Table H.30: Summary table of data for Erodium cicutarium (L.) L’Hérit..

Figure H.87: Effect of year on the probability of Erodium cicutarium (L.) L’Hérit. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.88: The same as H.87, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.89: Visualisation of the spatial smooth effect on the probability of Erodium cicutarium (L.) L’Hérit. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.31 Erophila verna (L.) Chevall.

Table H.31: Summary table of data for Erophila verna (L.) Chevall..

Figure H.90: Effect of year on the probability of Erophila verna (L.) Chevall. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.91: The same as H.90, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.92: Visualisation of the spatial smooth effect on the probability of Erophila verna (L.) Chevall. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.

H.32 Erysimum cheiranthoides L.

Table H.32: Summary table of data for Erysimum cheiranthoides L..

Figure H.93: Effect of year on the probability of Erysimum cheiranthoides L. presence in 1 km x 1 km squares where the species has been observed at least once. The fitted line shows the sum of the overall mean (the intercept), a conditional effect of list-length equal to 130 and the year-smoother. The vertical dashed lines indicate the year(s) where the year-smoother is zero. The 95% confidence band is shown in grey (including the variability around the intercept and the smoother). a: 1950 - 2018, b: 1990 - 2018.

Figure H.94: The same as H.93, but the vertical axis is scaled to the range of the predicted values such that relative changes can be seen more easily. a: 1950 - 2018, b: 1990 - 2018.

Figure H.95: Visualisation of the spatial smooth effect on the probability of Erysimum cheiranthoides L. presence in 1 km x 1 km squares where the species has been observed at least once. The probabilities (values on the contour lines) are conditional on the final year of observation and a list-length equal to 130. The dashed contour line demarcates zones where the species is expected to be more prevalent (red shades) from zones where the species is less prevalent (blue shades). a: 1950 - 2018, b: 1990 - 2018.