L Leontodon hispidus L. - Malus sylvestris (L.) Mill.

L.1 Leontodon hispidus L.

Table L.1: Summary table of data for Leontodon hispidus L..
Effect of year on the probability of Leontodon hispidus 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 L.1: Effect of year on the probability of Leontodon hispidus 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 L.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 L.2: The same as L.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 Leontodon hispidus 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 L.3: Visualisation of the spatial smooth effect on the probability of Leontodon hispidus 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.

L.2 Leontodon saxatilis Lam.

Table L.2: Summary table of data for Leontodon saxatilis Lam..
Effect of year on the probability of Leontodon saxatilis Lam. 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 L.4: Effect of year on the probability of Leontodon saxatilis Lam. 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 L.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 L.5: The same as L.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 Leontodon saxatilis Lam. 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 L.6: Visualisation of the spatial smooth effect on the probability of Leontodon saxatilis Lam. 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.

L.3 Lepidium ruderale L.

Table L.3: Summary table of data for Lepidium ruderale L..
Effect of year on the probability of Lepidium ruderale 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 L.7: Effect of year on the probability of Lepidium ruderale 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 L.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 L.8: The same as L.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 Lepidium ruderale 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 L.9: Visualisation of the spatial smooth effect on the probability of Lepidium ruderale 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.

L.4 Lepidium virginicum L.

Table L.4: Summary table of data for Lepidium virginicum L..
Effect of year on the probability of Lepidium virginicum 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 L.10: Effect of year on the probability of Lepidium virginicum 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 L.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 L.11: The same as L.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 Lepidium virginicum 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 L.12: Visualisation of the spatial smooth effect on the probability of Lepidium virginicum 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.

L.5 Leucanthemum vulgare Lam.

Table L.5: Summary table of data for Leucanthemum vulgare Lam..
Effect of year on the probability of Leucanthemum vulgare Lam. 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 L.13: Effect of year on the probability of Leucanthemum vulgare Lam. 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 L.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 L.14: The same as L.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 Leucanthemum vulgare Lam. 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 L.15: Visualisation of the spatial smooth effect on the probability of Leucanthemum vulgare Lam. 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.

L.6 Ligustrum ovalifolium Hassk.

Table L.6: Summary table of data for Ligustrum ovalifolium Hassk..
Effect of year on the probability of Ligustrum ovalifolium Hassk. 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 L.16: Effect of year on the probability of Ligustrum ovalifolium Hassk. 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 L.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 L.17: The same as L.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 Ligustrum ovalifolium Hassk. 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 L.18: Visualisation of the spatial smooth effect on the probability of Ligustrum ovalifolium Hassk. 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.

L.7 Ligustrum vulgare L.

Table L.7: Summary table of data for Ligustrum vulgare L..
Effect of year on the probability of Ligustrum 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 L.19: Effect of year on the probability of Ligustrum 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.

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

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

L.8 Linaria vulgaris Mill.

Table L.8: Summary table of data for Linaria vulgaris Mill..
Effect of year on the probability of Linaria vulgaris 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 L.22: Effect of year on the probability of Linaria vulgaris 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 L.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 L.23: The same as L.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 Linaria vulgaris 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 L.24: Visualisation of the spatial smooth effect on the probability of Linaria vulgaris 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.

L.9 Linum usitatissimum L.

Table L.9: Summary table of data for Linum usitatissimum L..
Effect of year on the probability of Linum usitatissimum 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 L.25: Effect of year on the probability of Linum usitatissimum 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 L.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 L.26: The same as L.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 Linum usitatissimum 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 L.27: Visualisation of the spatial smooth effect on the probability of Linum usitatissimum 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.

L.10 Listera ovata (L.) R. Brown

Table L.10: Summary table of data for Listera ovata (L.) R. Brown.
Effect of year on the probability of Listera ovata (L.) R. Brown 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 L.28: Effect of year on the probability of Listera ovata (L.) R. Brown 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 L.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 L.29: The same as L.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 Listera ovata (L.) R. Brown 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 L.30: Visualisation of the spatial smooth effect on the probability of Listera ovata (L.) R. Brown 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.

L.11 Lobularia maritima (L.) Desv.

Table L.11: Summary table of data for Lobularia maritima (L.) Desv..
Effect of year on the probability of Lobularia maritima (L.) 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 L.31: Effect of year on the probability of Lobularia maritima (L.) 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.

The same as L.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 L.32: The same as L.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 Lobularia maritima (L.) 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.

Figure L.33: Visualisation of the spatial smooth effect on the probability of Lobularia maritima (L.) 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.

L.12 Lolium multiflorum Lam.

Table L.12: Summary table of data for Lolium multiflorum Lam..
Effect of year on the probability of Lolium multiflorum Lam. 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 L.34: Effect of year on the probability of Lolium multiflorum Lam. 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 L.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 L.35: The same as L.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 Lolium multiflorum Lam. 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 L.36: Visualisation of the spatial smooth effect on the probability of Lolium multiflorum Lam. 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.

L.13 Lolium perenne L.

Table L.13: Summary table of data for Lolium perenne L..
Effect of year on the probability of Lolium perenne 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 L.37: Effect of year on the probability of Lolium perenne 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 L.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 L.38: The same as L.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 Lolium perenne 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 L.39: Visualisation of the spatial smooth effect on the probability of Lolium perenne 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.

L.14 Lonicera periclymenum L.

Table L.14: Summary table of data for Lonicera periclymenum L..
Effect of year on the probability of Lonicera periclymenum 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 L.40: Effect of year on the probability of Lonicera periclymenum 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 L.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 L.41: The same as L.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 Lonicera periclymenum 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 L.42: Visualisation of the spatial smooth effect on the probability of Lonicera periclymenum 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.

L.15 Lotus corniculatus L.

Table L.15: Summary table of data for Lotus corniculatus L..
Effect of year on the probability of Lotus corniculatus 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 L.43: Effect of year on the probability of Lotus corniculatus 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 L.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 L.44: The same as L.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 Lotus corniculatus 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 L.45: Visualisation of the spatial smooth effect on the probability of Lotus corniculatus 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.

L.16 Lotus pedunculatus Cav.

Table L.16: Summary table of data for Lotus pedunculatus Cav..
Effect of year on the probability of Lotus pedunculatus 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 L.46: Effect of year on the probability of Lotus pedunculatus 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 L.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 L.47: The same as L.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 Lotus pedunculatus 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 L.48: Visualisation of the spatial smooth effect on the probability of Lotus pedunculatus 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.

L.17 Lunaria annua L.

Table L.17: Summary table of data for Lunaria annua L..
Effect of year on the probability of Lunaria annua 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 L.49: Effect of year on the probability of Lunaria annua 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 L.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 L.50: The same as L.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 Lunaria annua 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 L.51: Visualisation of the spatial smooth effect on the probability of Lunaria annua 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.

L.18 Luzula campestris (L.) DC.

Table L.18: Summary table of data for Luzula campestris (L.) DC..
Effect of year on the probability of Luzula campestris (L.) DC. 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 L.52: Effect of year on the probability of Luzula campestris (L.) DC. 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 L.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 L.53: The same as L.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 Luzula campestris (L.) DC. 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 L.54: Visualisation of the spatial smooth effect on the probability of Luzula campestris (L.) DC. 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.

L.19 Luzula multiflora (Ehrh.) Lej.

Table L.19: Summary table of data for Luzula multiflora (Ehrh.) Lej..
Effect of year on the probability of Luzula multiflora (Ehrh.) Lej. 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 L.55: Effect of year on the probability of Luzula multiflora (Ehrh.) Lej. 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 L.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 L.56: The same as L.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 Luzula multiflora (Ehrh.) Lej. 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 L.57: Visualisation of the spatial smooth effect on the probability of Luzula multiflora (Ehrh.) Lej. 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.

L.20 Luzula pilosa (L.) Willd.

Table L.20: Summary table of data for Luzula pilosa (L.) Willd..
Effect of year on the probability of Luzula pilosa (L.) Willd. 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 L.58: Effect of year on the probability of Luzula pilosa (L.) Willd. 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 L.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 L.59: The same as L.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 Luzula pilosa (L.) Willd. 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 L.60: Visualisation of the spatial smooth effect on the probability of Luzula pilosa (L.) Willd. 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.

L.21 Luzula sylvatica (Huds.) Gaudin

Table L.21: Summary table of data for Luzula sylvatica (Huds.) Gaudin.
Effect of year on the probability of Luzula sylvatica (Huds.) Gaudin 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 L.61: Effect of year on the probability of Luzula sylvatica (Huds.) Gaudin 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 L.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 L.62: The same as L.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 Luzula sylvatica (Huds.) Gaudin 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 L.63: Visualisation of the spatial smooth effect on the probability of Luzula sylvatica (Huds.) Gaudin 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.

L.22 Lychnis coronaria (L.) Desr.

Table L.22: Summary table of data for Lychnis coronaria (L.) Desr..
Effect of year on the probability of Lychnis coronaria (L.) Desr. 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 L.64: Effect of year on the probability of Lychnis coronaria (L.) Desr. 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 L.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 L.65: The same as L.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 Lychnis coronaria (L.) Desr. 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 L.66: Visualisation of the spatial smooth effect on the probability of Lychnis coronaria (L.) Desr. 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.

L.23 Lychnis flos-cuculi L.

Table L.23: Summary table of data for Lychnis flos-cuculi L..
Effect of year on the probability of Lychnis flos-cuculi 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 L.67: Effect of year on the probability of Lychnis flos-cuculi 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 L.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 L.68: The same as L.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 Lychnis flos-cuculi 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 L.69: Visualisation of the spatial smooth effect on the probability of Lychnis flos-cuculi 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.

L.24 Lycopus europaeus L.

Table L.24: Summary table of data for Lycopus europaeus L..
Effect of year on the probability of Lycopus 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 L.70: Effect of year on the probability of Lycopus 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 L.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 L.71: The same as L.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 Lycopus 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 L.72: Visualisation of the spatial smooth effect on the probability of Lycopus 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.

L.25 Lysimachia nemorum L.

Table L.25: Summary table of data for Lysimachia nemorum L..
Effect of year on the probability of Lysimachia nemorum 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 L.73: Effect of year on the probability of Lysimachia nemorum 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 L.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 L.74: The same as L.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 Lysimachia nemorum 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 L.75: Visualisation of the spatial smooth effect on the probability of Lysimachia nemorum 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.

L.26 Lysimachia nummularia L.

Table L.26: Summary table of data for Lysimachia nummularia L..
Effect of year on the probability of Lysimachia nummularia 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 L.76: Effect of year on the probability of Lysimachia nummularia 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 L.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 L.77: The same as L.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 Lysimachia nummularia 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 L.78: Visualisation of the spatial smooth effect on the probability of Lysimachia nummularia 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.

L.27 Lysimachia vulgaris L.

Table L.27: Summary table of data for Lysimachia vulgaris L..
Effect of year on the probability of Lysimachia vulgaris 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 L.79: Effect of year on the probability of Lysimachia vulgaris 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 L.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 L.80: The same as L.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 Lysimachia vulgaris 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 L.81: Visualisation of the spatial smooth effect on the probability of Lysimachia vulgaris 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.

L.28 Lythrum portula (L.) D.A. Webb

Table L.28: Summary table of data for Lythrum portula (L.) D.A. Webb.
Effect of year on the probability of Lythrum portula (L.) D.A. Webb 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 L.82: Effect of year on the probability of Lythrum portula (L.) D.A. Webb 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 L.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 L.83: The same as L.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 Lythrum portula (L.) D.A. Webb 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 L.84: Visualisation of the spatial smooth effect on the probability of Lythrum portula (L.) D.A. Webb 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.

L.29 Lythrum salicaria L.

Table L.29: Summary table of data for Lythrum salicaria L..
Effect of year on the probability of Lythrum salicaria 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 L.85: Effect of year on the probability of Lythrum salicaria 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 L.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 L.86: The same as L.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 Lythrum salicaria 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 L.87: Visualisation of the spatial smooth effect on the probability of Lythrum salicaria 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.

L.30 Mahonia aquifolium (Pursh) Nutt.

Table L.30: Summary table of data for Mahonia aquifolium (Pursh) Nutt..
Effect of year on the probability of Mahonia aquifolium (Pursh) Nutt. 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 L.88: Effect of year on the probability of Mahonia aquifolium (Pursh) Nutt. 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 L.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 L.89: The same as L.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 Mahonia aquifolium (Pursh) Nutt. 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 L.90: Visualisation of the spatial smooth effect on the probability of Mahonia aquifolium (Pursh) Nutt. 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.

L.31 Maianthemum bifolium (L.) F.W. Schmidt

Table L.31: Summary table of data for Maianthemum bifolium (L.) F.W. Schmidt.
Effect of year on the probability of Maianthemum bifolium (L.) F.W. Schmidt 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 L.91: Effect of year on the probability of Maianthemum bifolium (L.) F.W. Schmidt 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 L.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 L.92: The same as L.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 Maianthemum bifolium (L.) F.W. Schmidt 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 L.93: Visualisation of the spatial smooth effect on the probability of Maianthemum bifolium (L.) F.W. Schmidt 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.

L.32 Malus sylvestris (L.) Mill.

Table L.32: Summary table of data for Malus sylvestris (L.) Mill..
Effect of year on the probability of Malus sylvestris (L.) 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 L.94: Effect of year on the probability of Malus sylvestris (L.) 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 L.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 L.95: The same as L.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 Malus sylvestris (L.) 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 L.96: Visualisation of the spatial smooth effect on the probability of Malus sylvestris (L.) 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.