README These data contain observations of cabbage stem flea beetle abundance, associated crop damage and levels of larval infestation, and arthropod trapping data for 14 oilseed rape fields in Lower Saxony, Germany. All observations were made between September and November 2021. Data were analysed in response to a range of calculated landscape metrics. Data are associated with the bioRxiv article: "Sharing the burden: Cabbage stem flea beetle pest pressure and crop damage are lower in rapeseed fields surrounded by other rapeseed crops" By D J Leybourne, A M C Pahl, P Melloh, and E A Martin. To protect privacy of the participating farmers we are unable to publicly share any geodata associated with the project. However, this information can be requested by contacting the corresponding author. Associated datafiles: 1) "Rapeseed_Site_Info.csv" - Dataset containing agronomic and other informaiton for each field site 2) "Pest_Pressure.csv" - Dataset containing the pest pressure data for each field site. Column description: field (field ID), Transect (Edge or Crop), Week (calander number for 2021), Beetle_Pressure (beetle pressure, total number of cabbage stem flea beetle, for the given week), and Pressure_Cummulative (cumulative number of cabbage stem flea beetle) 3) "Adult_Damage.csv" - Dataset containing the leaf damage observations for each field site Column description: field (field ID), Transect (Edge or Crop), Assessment_Area (Quadrat), Assessed_Plant (The plant being examined, 1-20), Assessment_Round (1-2), Assessment_Month (Sep or Oct), and Percent_Damage (The percent of cabbage stem flea beetle damage observed on each plant). 4) "Larva_Load.csv" - Dataset containing the larvae infestaion (larvae load) data for each field site Column description: field (field ID), Transect (Edge or Crop),Assessment_Area (Quadrat), Assessed_Plant (The plant being examined, 1-5), Plant_Diameter (the diamter of the plant stem, cm), Larvae_Stage_1 (The number of first-instar larvae), Larvae_Stage_2 (The number of second-instar larvae), Larvae_Stage_3 (The number of third-instar larvae), Larvae_Stage_4 (The number of fourth-instar larvae), Larvae_Total (the total number of all larvae) 5) "Insect_Traps.csv" - Dataset containing the data detailing the insect trapping observations Column description: field (field ID), Transect (Edge or Crop), Trap_Type (Associated trap, pan or pitfall), Assessment_Area (Quadrat), Assessment_Round (1-3), Assessment_Month (Sep, Oct, Nov). All other columns contain abundance information for the stated taxa: Arthropod Order (Column F-Q); total abundance (Column R); Family for Coleoptera (S-AA), Araneae (AB-AG), and Hymenoptera (AH-AO) 6) "Landscape_500m.csv" - Dataset detailing the landscape data, characterised at the 500 m spatial scale 7) "Landscape_750m.csv" - Dataset detailing the landscape data, characterised at the 750 m spatial scale 8) "Landscape_1000m.csv" - Dataset detailing the landscape data, characterised at the 1000 m spatial scale Column descriptions for all Landscape datasets: Field_ID (Field ID), Columns then detail a series of landscape features for several crops. See descriptions below. 9) "Landscape_500m_Radii_Correlations.csv" - Dataset detailing landscape correlations at the 500 m spatial scale 10) "Landscape_750m_Radii_Correlations.csv" - Dataset detailing landscape correlations at the 750 m spatial scale 11) "Landscape_1000m_Radii_Correlations.csv" - Dataset detailing landscape correlations at the 1000 m spatial scale Column descriptions for all Landscape correlation datasets: Field_ID (Field ID), RP_2021 (proportion of oilseed rape grown in previous year), CD_2021 (crop heterogeneity in previous year), RP_2022 (proportion of oilseed rape grown in experimental year), CD_2022 (crop heterogeneity in experimental year), ChangeRaps (Ratio change in rapeseed proportion between years), SNH (proportion of seminatural habitat in the landscape) Codes and abbreviations used in the datafiles: - In all landscape datafiles: 1) OSR = Oilseed rape 2) Adj_OSR = Presence/absence of an adjacent oilseed rape field 3) No_ = Number 4) Prop_ = Proportion of associated landscape feature in the assessed area 5) Config_ = Configuration of the associated landscape feature in the assessed area (i.e., number of fields) 6) Y2021_ = Harvest year 2021 ("Previous season") 7) Y2022_ = Harvest year 2022 ("Experimental season") 8) CropDiv21 = Total crop heterogeneity in the previous season 9) CropDiv22 = Total crop heterogeneity in the experimental season 10) TotalConfig21 = Total configuration of all crop fields in the assessed area in the previous season 11) TotalConfig22 = Total configuration of all crop fields in the assessed are in the experimental season 12) ChangeRapsRatio = The ratio change in oilseed rape production between the two seasons 13) veg02_f and veg03_f = Total proportion of forest, woodland, heath, moor, swamp, and uncultivated land (extracted from Lower Saxpny ATKIS database) 14) SNH = Semi-natural habitat (The sum of veg02_f and veg03_f) - In datafile "Insect_Traps.csv": 1) Col = Coleoptera 2) Ara = Aranae 3) Hym = Hymenoptera Associated R Script: 1) "Landscape_Beetle_Analysis_Code_Final.R" - Script containing the code used for data analysis and creation of graphs Note: Session Info for R/R Studio used for analysis. R version 4.2.0 (2022-04-22 ucrt) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 10 x64 (build 19045) attached packages: [1] GGally_2.1.2 piecewiseSEM_2.3.0 DHARMa_0.4.6 MuMIn_1.47.5 ggpubr_0.6.0 vegan_2.6-4 [7] lattice_0.21-8 permute_0.9-7 car_3.1-2 carData_3.0-5 performance_0.10.4 glmmTMB_1.1.7 [13] lme4_1.1-34 Matrix_1.6-1 lubridate_1.9.2 forcats_1.0.0 stringr_1.5.0 dplyr_1.1.2 [19] purrr_1.0.2 readr_2.1.4 tidyr_1.3.0 tibble_3.2.1 ggplot2_3.4.3 tidyverse_2.0.0 loaded via a namespace (and not attached): [1] TH.data_1.1-2 minqa_1.2.5 colorspace_2.1-0 ggsignif_0.6.4 ellipsis_0.3.2 class_7.3-22 [7] estimability_1.4.1 rstudioapi_0.15.0 proxy_0.4-27 farver_2.1.1 fansi_1.0.4 mvtnorm_1.2-2 [13] codetools_0.2-19 splines_4.2.0 doParallel_1.0.17 jsonlite_1.8.7 nloptr_2.0.3 pbkrtest_0.5.2 [19] broom_1.0.5 cluster_2.1.4 shiny_1.7.5 DiagrammeR_1.0.10 compiler_4.2.0 emmeans_1.8.8 [25] backports_1.4.1 fastmap_1.1.1 cli_3.6.1 later_1.3.1 prettyunits_1.1.1 visNetwork_2.1.2 [31] htmltools_0.5.6 tools_4.2.0 igraph_1.5.1 coda_0.19-4 gtable_0.3.4 glue_1.6.2 [37] Rcpp_1.0.11 vctrs_0.6.3 nlme_3.1-163 iterators_1.0.14 insight_0.19.3 rbibutils_2.2.15 [43] timechange_0.2.0 mime_0.12 lifecycle_1.0.3 gap.datasets_0.0.6 rstatix_0.7.2 MASS_7.3-60 [49] zoo_1.8-12 scales_1.2.1 ragg_1.2.5 hms_1.1.3 promises_1.2.1 parallel_4.2.0 [55] sandwich_3.0-2 TMB_1.9.6 RColorBrewer_1.1-3 gridExtra_2.3 reshape_0.8.9 stringi_1.7.12 [61] gap_1.5-3 foreach_1.5.2 e1071_1.7-13 boot_1.3-28.1 qgam_1.3.4 systemfonts_1.0.4 [67] Rdpack_2.5 rlang_1.1.1 pkgconfig_2.0.3 sf_1.0-14 labeling_0.4.2 htmlwidgets_1.6.2 [73] cowplot_1.1.1 tidyselect_1.2.0 plyr_1.8.8 magrittr_2.0.3 R6_2.5.1 generics_0.1.3 [79] multcomp_1.4-25 DBI_1.1.3 pillar_1.9.0 withr_2.5.0 mgcv_1.9-0 units_0.8-3 [85] survival_3.5-7 datawizard_0.8.0 abind_1.4-7 crayon_1.5.2 KernSmooth_2.23-22 utf8_1.2.3 [91] tzdb_0.4.0 progress_1.2.2 grid_4.2.0 digest_0.6.33 classInt_0.4-9 xtable_1.8-6 [97] httpuv_1.6.11 numDeriv_2020.2-1 textshaping_0.3.6 stats4_4.2.0 munsell_0.5.0