I Eupatorium cannabinum L. - Geum urbanum L.

I.1 Eupatorium cannabinum L.

Table I.1: Summary table of data for Eupatorium cannabinum L..
Effect of year on the probability of Eupatorium cannabinum 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 I.1: Effect of year on the probability of Eupatorium cannabinum 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.

The same as I.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 I.2: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Eupatorium cannabinum 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.

Figure I.3: Visualisation of the spatial smooth effect on the probability of Eupatorium cannabinum 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.

I.2 Euphorbia helioscopia L.

Table I.2: Summary table of data for Euphorbia helioscopia L..
Effect of year on the probability of Euphorbia helioscopia 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 I.4: Effect of year on the probability of Euphorbia helioscopia 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.

The same as I.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 I.5: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Euphorbia helioscopia 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.

Figure I.6: Visualisation of the spatial smooth effect on the probability of Euphorbia helioscopia 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.

I.3 Euphorbia lathyris L.

Table I.3: Summary table of data for Euphorbia lathyris L..
Effect of year on the probability of Euphorbia lathyris 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 I.7: Effect of year on the probability of Euphorbia lathyris 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.

The same as I.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 I.8: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Euphorbia lathyris 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.

Figure I.9: Visualisation of the spatial smooth effect on the probability of Euphorbia lathyris 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.

I.4 Euphorbia peplus L.

Table I.4: Summary table of data for Euphorbia peplus L..
Effect of year on the probability of Euphorbia peplus 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 I.10: Effect of year on the probability of Euphorbia peplus 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.

The same as I.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 I.11: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Euphorbia peplus 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.

Figure I.12: Visualisation of the spatial smooth effect on the probability of Euphorbia peplus 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.

I.5 Euonymus europaeus L.

Table I.5: Summary table of data for Euonymus europaeus L..
Effect of year on the probability of Euonymus europaeus 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 I.13: Effect of year on the probability of Euonymus europaeus 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.

The same as I.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 I.14: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Euonymus europaeus 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.

Figure I.15: Visualisation of the spatial smooth effect on the probability of Euonymus europaeus 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.

I.6 Fagus sylvatica L.

Table I.6: Summary table of data for Fagus sylvatica L..
Effect of year on the probability of Fagus sylvatica 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 I.16: Effect of year on the probability of Fagus sylvatica 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.

The same as I.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 I.17: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Fagus sylvatica 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.

Figure I.18: Visualisation of the spatial smooth effect on the probability of Fagus sylvatica 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.

I.7 Festuca arundinacea Schreb.

Table I.7: Summary table of data for Festuca arundinacea Schreb..
Effect of year on the probability of Festuca arundinacea 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 I.19: Effect of year on the probability of Festuca arundinacea 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.

The same as I.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 I.20: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Festuca arundinacea 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.

Figure I.21: Visualisation of the spatial smooth effect on the probability of Festuca arundinacea 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.

I.8 Festuca gigantea (L.) Vill.

Table I.8: Summary table of data for Festuca gigantea (L.) Vill..
Effect of year on the probability of Festuca gigantea (L.) Vill. 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 I.22: Effect of year on the probability of Festuca gigantea (L.) Vill. 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.

The same as I.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 I.23: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Festuca gigantea (L.) Vill. 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.

Figure I.24: Visualisation of the spatial smooth effect on the probability of Festuca gigantea (L.) Vill. 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.

I.9 Festuca pratensis Huds.

Table I.9: Summary table of data for Festuca pratensis Huds..
Effect of year on the probability of Festuca pratensis Huds. 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 I.25: Effect of year on the probability of Festuca pratensis Huds. 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.

The same as I.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 I.26: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Festuca pratensis Huds. 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.

Figure I.27: Visualisation of the spatial smooth effect on the probability of Festuca pratensis Huds. 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.

I.10 Filago minima (Smith) Pers.

Table I.10: Summary table of data for Filago minima (Smith) Pers..
Effect of year on the probability of Filago minima (Smith) Pers. 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 I.28: Effect of year on the probability of Filago minima (Smith) Pers. 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.

The same as I.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 I.29: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Filago minima (Smith) Pers. 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.

Figure I.30: Visualisation of the spatial smooth effect on the probability of Filago minima (Smith) Pers. 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.

I.11 Filipendula ulmaria (L.) Maxim.

Table I.11: Summary table of data for Filipendula ulmaria (L.) Maxim..
Effect of year on the probability of Filipendula ulmaria (L.) Maxim. 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 I.31: Effect of year on the probability of Filipendula ulmaria (L.) Maxim. 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.

The same as I.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 I.32: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Filipendula ulmaria (L.) Maxim. 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.

Figure I.33: Visualisation of the spatial smooth effect on the probability of Filipendula ulmaria (L.) Maxim. 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.

I.12 Foeniculum vulgare Mill.

Table I.12: Summary table of data for Foeniculum vulgare Mill..
Effect of year on the probability of Foeniculum vulgare Mill. 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 I.34: Effect of year on the probability of Foeniculum vulgare Mill. 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.

The same as I.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 I.35: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Foeniculum vulgare Mill. 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.

Figure I.36: Visualisation of the spatial smooth effect on the probability of Foeniculum vulgare Mill. 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.

I.13 Fragaria vesca L.

Table I.13: Summary table of data for Fragaria vesca L..
Effect of year on the probability of Fragaria vesca 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 I.37: Effect of year on the probability of Fragaria vesca 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.

The same as I.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 I.38: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Fragaria vesca 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.

Figure I.39: Visualisation of the spatial smooth effect on the probability of Fragaria vesca 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.

I.14 Fraxinus excelsior L.

Table I.14: Summary table of data for Fraxinus excelsior L..
Effect of year on the probability of Fraxinus excelsior 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 I.40: Effect of year on the probability of Fraxinus excelsior 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.

The same as I.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 I.41: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Fraxinus excelsior 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.

Figure I.42: Visualisation of the spatial smooth effect on the probability of Fraxinus excelsior 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.

I.15 Fumaria officinalis L.

Table I.15: Summary table of data for Fumaria officinalis L..
Effect of year on the probability of Fumaria officinalis 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 I.43: Effect of year on the probability of Fumaria officinalis 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.

The same as I.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.

Figure I.44: The same as I.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.

Visualisation of the spatial smooth effect on the probability of Fumaria officinalis 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.

Figure I.45: Visualisation of the spatial smooth effect on the probability of Fumaria officinalis 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.

I.16 Galanthus nivalis L.

Table I.16: Summary table of data for Galanthus nivalis L..
Effect of year on the probability of Galanthus nivalis 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 I.46: Effect of year on the probability of Galanthus nivalis 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.

The same as I.46, 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 I.47: The same as I.46, 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.

Visualisation of the spatial smooth effect on the probability of Galanthus nivalis 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.

Figure I.48: Visualisation of the spatial smooth effect on the probability of Galanthus nivalis 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.

I.17 Galinsoga parviflora Cav.

Table I.17: Summary table of data for Galinsoga parviflora Cav..
Effect of year on the probability of Galinsoga parviflora Cav. 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 I.49: Effect of year on the probability of Galinsoga parviflora Cav. 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.

The same as I.49, 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 I.50: The same as I.49, 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.

Visualisation of the spatial smooth effect on the probability of Galinsoga parviflora Cav. 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.

Figure I.51: Visualisation of the spatial smooth effect on the probability of Galinsoga parviflora Cav. 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.

I.18 Galinsoga quadriradiata Ruiz et Pav.

Table I.18: Summary table of data for Galinsoga quadriradiata Ruiz et Pav..
Effect of year on the probability of Galinsoga quadriradiata Ruiz et Pav. 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 I.52: Effect of year on the probability of Galinsoga quadriradiata Ruiz et Pav. 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.

The same as I.52, 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 I.53: The same as I.52, 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.

Visualisation of the spatial smooth effect on the probability of Galinsoga quadriradiata Ruiz et Pav. 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.

Figure I.54: Visualisation of the spatial smooth effect on the probability of Galinsoga quadriradiata Ruiz et Pav. 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.

I.19 Galium aparine L.

Table I.19: Summary table of data for Galium aparine L..
Effect of year on the probability of Galium aparine 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 I.55: Effect of year on the probability of Galium aparine 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.

The same as I.55, 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 I.56: The same as I.55, 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.

Visualisation of the spatial smooth effect on the probability of Galium aparine 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.

Figure I.57: Visualisation of the spatial smooth effect on the probability of Galium aparine 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.

I.20 Galium mollugo L.

Table I.20: Summary table of data for Galium mollugo L..
Effect of year on the probability of Galium mollugo 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 I.58: Effect of year on the probability of Galium mollugo 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.

The same as I.58, 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 I.59: The same as I.58, 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.

Visualisation of the spatial smooth effect on the probability of Galium mollugo 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.

Figure I.60: Visualisation of the spatial smooth effect on the probability of Galium mollugo 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.

I.21 Galium palustre L.

Table I.21: Summary table of data for Galium palustre L..
Effect of year on the probability of Galium 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 I.61: Effect of year on the probability of Galium 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.

The same as I.61, 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 I.62: The same as I.61, 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.

Visualisation of the spatial smooth effect on the probability of Galium 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.

Figure I.63: Visualisation of the spatial smooth effect on the probability of Galium 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.

I.22 Galium saxatile L.

Table I.22: Summary table of data for Galium saxatile L..
Effect of year on the probability of Galium saxatile 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 I.64: Effect of year on the probability of Galium saxatile 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.

The same as I.64, 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 I.65: The same as I.64, 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.

Visualisation of the spatial smooth effect on the probability of Galium saxatile 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.

Figure I.66: Visualisation of the spatial smooth effect on the probability of Galium saxatile 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.

I.23 Galium uliginosum L.

Table I.23: Summary table of data for Galium uliginosum L..
Effect of year on the probability of Galium uliginosum 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 I.67: Effect of year on the probability of Galium uliginosum 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.

The same as I.67, 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 I.68: The same as I.67, 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.

Visualisation of the spatial smooth effect on the probability of Galium uliginosum 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.

Figure I.69: Visualisation of the spatial smooth effect on the probability of Galium uliginosum 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.

I.24 Galium verum L.

Table I.24: Summary table of data for Galium verum L..
Effect of year on the probability of Galium verum 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 I.70: Effect of year on the probability of Galium verum 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.

The same as I.70, 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 I.71: The same as I.70, 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.

Visualisation of the spatial smooth effect on the probability of Galium verum 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.

Figure I.72: Visualisation of the spatial smooth effect on the probability of Galium verum 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.

I.25 Genista anglica L.

Table I.25: Summary table of data for Genista anglica L..
Effect of year on the probability of Genista anglica 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 I.73: Effect of year on the probability of Genista anglica 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.

The same as I.73, 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 I.74: The same as I.73, 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.

Visualisation of the spatial smooth effect on the probability of Genista anglica 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.

Figure I.75: Visualisation of the spatial smooth effect on the probability of Genista anglica 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.

I.26 Geranium dissectum L.

Table I.26: Summary table of data for Geranium dissectum L..
Effect of year on the probability of Geranium dissectum 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 I.76: Effect of year on the probability of Geranium dissectum 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.

The same as I.76, 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 I.77: The same as I.76, 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.

Visualisation of the spatial smooth effect on the probability of Geranium dissectum 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.

Figure I.78: Visualisation of the spatial smooth effect on the probability of Geranium dissectum 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.

I.27 Geranium molle L.

Table I.27: Summary table of data for Geranium molle L..
Effect of year on the probability of Geranium molle 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 I.79: Effect of year on the probability of Geranium molle 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.

The same as I.79, 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 I.80: The same as I.79, 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.

Visualisation of the spatial smooth effect on the probability of Geranium molle 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.

Figure I.81: Visualisation of the spatial smooth effect on the probability of Geranium molle 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.

I.28 Geranium purpureum Vill.

Table I.28: Summary table of data for Geranium purpureum Vill..
Effect of year on the probability of Geranium purpureum Vill. 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 I.82: Effect of year on the probability of Geranium purpureum Vill. 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.

The same as I.82, 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 I.83: The same as I.82, 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.

Visualisation of the spatial smooth effect on the probability of Geranium purpureum Vill. 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.

Figure I.84: Visualisation of the spatial smooth effect on the probability of Geranium purpureum Vill. 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.

I.29 Geranium pusillum L.

Table I.29: Summary table of data for Geranium pusillum L..
Effect of year on the probability of Geranium pusillum 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 I.85: Effect of year on the probability of Geranium pusillum 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.

The same as I.85, 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 I.86: The same as I.85, 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.

Visualisation of the spatial smooth effect on the probability of Geranium pusillum 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.

Figure I.87: Visualisation of the spatial smooth effect on the probability of Geranium pusillum 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.

I.30 Geranium pyrenaicum Burm. f.

Table I.30: Summary table of data for Geranium pyrenaicum Burm. f..
Effect of year on the probability of Geranium pyrenaicum Burm. f. 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 I.88: Effect of year on the probability of Geranium pyrenaicum Burm. f. 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.

The same as I.88, 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 I.89: The same as I.88, 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.

Visualisation of the spatial smooth effect on the probability of Geranium pyrenaicum Burm. f. 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.

Figure I.90: Visualisation of the spatial smooth effect on the probability of Geranium pyrenaicum Burm. f. 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.

I.31 Geranium robertianum L.

Table I.31: Summary table of data for Geranium robertianum L..
Effect of year on the probability of Geranium robertianum 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 I.91: Effect of year on the probability of Geranium robertianum 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.

The same as I.91, 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 I.92: The same as I.91, 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.

Visualisation of the spatial smooth effect on the probability of Geranium robertianum 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.

Figure I.93: Visualisation of the spatial smooth effect on the probability of Geranium robertianum 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.

I.32 Geum urbanum L.

Table I.32: Summary table of data for Geum urbanum L..
Effect of year on the probability of Geum urbanum 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 I.94: Effect of year on the probability of Geum urbanum 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.

The same as I.94, 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 I.95: The same as I.94, 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.

Visualisation of the spatial smooth effect on the probability of Geum urbanum 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.

Figure I.96: Visualisation of the spatial smooth effect on the probability of Geum urbanum 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.