# Q Reseda luteola L. - Rumex sanguineus L.

## Q.1Reseda luteola L.

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

## Q.2Reseda lutea L.

(ref:distinct-squares-Reselute.1) Summary table of data for Reseda lutea L..

(#tab:distinct-squares-Reselute.1)(ref:distinct-squares-Reselute.1)

(ref:p-year-Reselute.1) Effect of year on the probability of Reseda lutea 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. (#fig:p-year-Reselute.1)(ref:p-year-Reselute.1)

(ref:p-year-freey-Reselute.1) The same as @ref(fig:p-year-Reselute.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. (#fig:p-year-freey-Reselute.1)(ref:p-year-freey-Reselute.1)

(ref:p-space-Reselute.1) Visualisation of the spatial smooth effect on the probability of Reseda lutea 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. (#fig:p-space-Reselute.1)(ref:p-space-Reselute.1)

## Q.3Frangula alnus Mill.

Table Q.2: Summary table of data for Frangula alnus Mill.. Figure Q.4: Effect of year on the probability of Frangula alnus 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 Q.5: The same as Q.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 Q.6: Visualisation of the spatial smooth effect on the probability of Frangula alnus 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.

## Q.4Rhinanthus angustifolius C.C. Gmel.

Table Q.3: Summary table of data for Rhinanthus angustifolius C.C. Gmel.. Figure Q.7: Effect of year on the probability of Rhinanthus angustifolius C.C. Gmel. 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 Q.8: The same as Q.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 Q.9: Visualisation of the spatial smooth effect on the probability of Rhinanthus angustifolius C.C. Gmel. 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.

## Q.5Rhododendron ponticum L.

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

## Q.6Rhus hirta (L.) Sudworth

Table Q.5: Summary table of data for Rhus hirta (L.) Sudworth. Figure Q.13: Effect of year on the probability of Rhus hirta (L.) Sudworth 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 Q.14: The same as Q.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 Q.15: Visualisation of the spatial smooth effect on the probability of Rhus hirta (L.) Sudworth 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.

## Q.7Ribes nigrum L.

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

## Q.8Ribes rubrum L.

Table Q.7: Summary table of data for Ribes rubrum L.. Figure Q.19: Effect of year on the probability of Ribes rubrum 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 Q.20: The same as Q.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 Q.21: Visualisation of the spatial smooth effect on the probability of Ribes rubrum 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.

## Q.9Ribes uva-crispa L.

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

## Q.10Robinia pseudoacacia L.

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

## Q.11Rorippa amphibia (L.) Besser

Table Q.10: Summary table of data for Rorippa amphibia (L.) Besser. Figure Q.28: Effect of year on the probability of Rorippa amphibia (L.) Besser 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 Q.29: The same as Q.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 Q.30: Visualisation of the spatial smooth effect on the probability of Rorippa amphibia (L.) Besser 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.

## Q.12Nasturtium microphyllum (Boenningh.) Reichenb.

Table Q.11: Summary table of data for Nasturtium microphyllum (Boenningh.) Reichenb.. Figure Q.31: Effect of year on the probability of Nasturtium microphyllum (Boenningh.) Reichenb. 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 Q.32: The same as Q.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 Q.33: Visualisation of the spatial smooth effect on the probability of Nasturtium microphyllum (Boenningh.) Reichenb. 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.

## Q.13Nasturtium officinale R. Brown

Table Q.12: Summary table of data for Nasturtium officinale R. Brown. Figure Q.34: Effect of year on the probability of Nasturtium officinale 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 Q.35: The same as Q.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 Q.36: Visualisation of the spatial smooth effect on the probability of Nasturtium officinale 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.

## Q.14Rorippa palustris (L.) Besser

Table Q.13: Summary table of data for Rorippa palustris (L.) Besser. Figure Q.37: Effect of year on the probability of Rorippa palustris (L.) Besser 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 Q.38: The same as Q.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 Q.39: Visualisation of the spatial smooth effect on the probability of Rorippa palustris (L.) Besser 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.

## Q.15Rorippa sylvestris (L.) Besser

Table Q.14: Summary table of data for Rorippa sylvestris (L.) Besser. Figure Q.40: Effect of year on the probability of Rorippa sylvestris (L.) Besser 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 Q.41: The same as Q.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 Q.42: Visualisation of the spatial smooth effect on the probability of Rorippa sylvestris (L.) Besser 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.

## Q.16Rosa arvensis Huds.

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

## Q.17Rosa canina L.

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

## Q.18Rosa rubiginosa L.

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

## Q.19Rosa rugosa Thunb.

Table Q.18: Summary table of data for Rosa rugosa Thunb.. Figure Q.52: Effect of year on the probability of Rosa rugosa Thunb. 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 Q.53: The same as Q.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 Q.54: Visualisation of the spatial smooth effect on the probability of Rosa rugosa Thunb. 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.

## Q.20Rubus caesius L.

Table Q.19: Summary table of data for Rubus caesius L.. Figure Q.55: Effect of year on the probability of Rubus caesius 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 Q.56: The same as Q.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 Q.57: Visualisation of the spatial smooth effect on the probability of Rubus caesius 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.

## Q.21Rubus idaeus L.

Table Q.20: Summary table of data for Rubus idaeus L.. Figure Q.58: Effect of year on the probability of Rubus idaeus 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 Q.59: The same as Q.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 Q.60: Visualisation of the spatial smooth effect on the probability of Rubus idaeus 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.

## Q.22Rubus laciniatus Willd.

Table Q.21: Summary table of data for Rubus laciniatus Willd.. Figure Q.61: Effect of year on the probability of Rubus laciniatus 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 Q.62: The same as Q.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 Q.63: Visualisation of the spatial smooth effect on the probability of Rubus laciniatus 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.

## Q.23Rubus nessensis W. Hall

Table Q.22: Summary table of data for Rubus nessensis W. Hall. Figure Q.64: Effect of year on the probability of Rubus nessensis W. Hall 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 Q.65: The same as Q.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 Q.66: Visualisation of the spatial smooth effect on the probability of Rubus nessensis W. Hall 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.

## Q.24Rumex acetosa L.

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

## Q.25Rumex acetosella L.

(ref:distinct-squares-Rumeacet.1) Summary table of data for Rumex acetosella L..

(#tab:distinct-squares-Rumeacet.1)(ref:distinct-squares-Rumeacet.1)

(ref:p-year-Rumeacet.1) Effect of year on the probability of Rumex acetosella 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. (#fig:p-year-Rumeacet.1)(ref:p-year-Rumeacet.1)

(ref:p-year-freey-Rumeacet.1) The same as @ref(fig:p-year-Rumeacet.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. (#fig:p-year-freey-Rumeacet.1)(ref:p-year-freey-Rumeacet.1)

(ref:p-space-Rumeacet.1) Visualisation of the spatial smooth effect on the probability of Rumex acetosella 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. (#fig:p-space-Rumeacet.1)(ref:p-space-Rumeacet.1)

## Q.26Rumex conglomeratus Murray

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

## Q.27Rumex crispus L.

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

## Q.28Rumex hydrolapathum Huds.

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

## Q.29Rumex maritimus L.

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

## Q.30Rumex obtusifolius L.

Table Q.28: Summary table of data for Rumex obtusifolius L.. Figure Q.82: Effect of year on the probability of Rumex obtusifolius 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 Q.83: The same as Q.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 Q.84: Visualisation of the spatial smooth effect on the probability of Rumex obtusifolius 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.

## Q.31Rumex palustris Smith

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

## Q.32Rumex sanguineus L.

Table Q.30: Summary table of data for Rumex sanguineus L.. Figure Q.88: Effect of year on the probability of Rumex sanguineus 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 Q.89: The same as Q.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 Q.90: Visualisation of the spatial smooth effect on the probability of Rumex sanguineus 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.