DOI: 10.6084/m9.figshare.25260232
title: Nicholson_etal_2024_LandscapesOfRisk

description: These data include observed pesticide residues in pollens gathered by honey bees (<i>Apis mellifera</i>), bumble bees (<i>Bombus</i><i> </i><i>terrestris</i>), and mason bees (<i>Osmia bicornis</i>) across 41 sites in southern Sweden from 2019 through 2021 (site network centroid: 55°43'09.2"N, 13°47'12.2"E). From these observed compound concentrations we calculated toxicity-weighted exposure (<i>TWE</i>) to estimate the direct pesticide risk to bees (Knapp et al. 2023), where the <i>TWE</i> for each compound <i>(TWE</i><sub><em>i</em></sub>) is the ratio of its exposure value (<i>c</i><sub><em>i</em></sub>) and its respective acute toxicity endpoint (LD<sub>50i</sub>). We summed TWE<sub>i</sub>s to calculate risk (see Materials and Methods).We analyzed the landscape surrounding our sites at three spatial scales according to the average foraging range for our three genera (<i>Osmia</i>: 500; <i>Bombus</i>: 1,500; and <i>Apis: </i>2,000 m). Landscape metrics were: proportion agricultural land, the proportion of pesticide-treated area, calculated landscape crop diversity, the interspersion and juxtaposition index, edge density, mean patch area (see Materials and Methods).Using a national data set of pesticide use and agricultural parcel data we calculated landscape-level metrics of pesticide use: total landscape pesticide load, total landscape toxic load (see Materials and Methods).We also used the spatialized pesticide use data as input a previously developed and validated mechanistic model (Lonsdorf et al., 2024) to predict bees' pesticide exposure and subsequent risk (i.e., summed TWE<sub>i</sub>s, as above).These data were used to compare the ability of the three classes of landscape-scale variables to predict this observed pesticide risk: 1) landscape composition and configuration metrics, 2) landscape load based on national pesticide use data, and 3) predictions from the bee pesticide exposure model. For the second and third class we distinguish between potential and realized load and risk (see Materials and Methods).The associated paper has been been published as an article in <i>Journal of Applied Ecology</i>.The columns included in this data are:<b>SiteID </b>- Field Identifier for the sites used in this analysis.<b>BeeSpecies </b>- the bee species from which pollen samples were collected. BB: <i>Bombus terrestris</i>; HB: <i>Apis melifera</i>; SB: <i>Osmia bicornis.</i><b>Year </b>- When the study took place.<b>Month </b>- Month number of when pollen samples were collected.<b>ObservedPesticideRisk </b>- pesticide risk (summed toxicity weighted concentrations) observed in pollen samples.<b>TotalLandscapeLoad_Realized </b>- The amount of different compounds applied in the landscape within the foraging range of the bee species from which the pollen sample was collected. This variable is calculated for all compounds that <i>were detected</i> in the sample (i.e., realized load). A log transformed and z-standardized version of this variable is also included (this is what was used in statistical models).<b>TotalLandscapeLoad_Potential </b>- The amount of different compounds applied in the landscape within the foraging range of the bee species from which the pollen sample was collected. This variable is calculated for all compounds that <i>were screened for</i> in the sample (i.e., potential load). A log transformed and z-standardized version of this variable is also included (this is what was used in statistical models).<b>TotalLandscapeToxicLoad_Realized </b>- The amount of different compounds applied in the landscape <i>weighting by their toxicity</i> within the foraging range of the bee species from which the pollen sample was collected. This variable is calculated for all compounds that <i>were detected</i> in the sample (i.e., realized load). A log transformed and z-standardized version of this variable is also included (this is what was used in statistical models).<b>TotalLandscapeToxicLoad_Potential </b>- The amount of different compounds applied in the landscape <i>weighting by their toxicity</i> within the foraging range of the bee species from which the pollen sample was collected. This variable is calculated for all compounds that <i>were screened for</i> in the sample (i.e., potential load). A log transformed and z-standardized version of this variable is also included (this is what was used in statistical models).<b>TotalLandscapeRisk_Realized </b>- The level of different compounds that bees are predicted to encounter based on a spatially explicit model of landscape exposure (SEMLE) for bees (Lonsdorf et al. 2024). These exposure predictions are weighted by their toxicity and summed (as above). This variable is calculated for all compounds that <i>were detected</i> in the sample (i.e., realized load). A log transformed and z-standardized version of this variable is also included (this is what was used in statistical models).<b>TotalLandscapeRisk_Potential </b>- The level of different compounds that bees are predicted to encounter based on a spatially explicit model of landscape exposure (SEMLE) for bees (Lonsdorf et al. 2024). These exposure predictions are weighted by their toxicity and summed (as above). This variable is calculated for all compounds that <i>were </i><i>s</i><i>creened for </i>in the sample (i.e., realized load). A log transformed and z-standardized version of this variable is also included (this is what was used in statistical models).<b>propAg </b>- The proportion of agricultural land within the foraging range of the bee species from which the pollen sample was collected. A z-standardized version of this variable is also included (this is what was used in statistical models).<b>propTreatedAg </b>- The proportion of pesticide-treated agricultural land within the foraging range of the bee species from which the pollen sample was collected. A z-standardized version of this variable is also included (this is what was used in statistical models).<b>CropDiversity </b>- Shannon index of crop types (richness) and their coverage (evenness) within the foraging range of the bee species from which the pollen sample was collected. A z-standardized version of this variable is also included (this is what was used in statistical models).<b>iji </b>- interspersion and juxtaposition index (McGarigal & Marks, 1995). A z-standardized version of this variable is also included (this is what was used in statistical models).References:Knapp, J. L., Nicholson, C. C., Jonsson, O., de Miranda, J. R., & Rundlöf, M. (2023). Ecological traits interact with landscape context to determine bees’ pesticide risk. <i>Nature Ecology & Evolution</i>, <i>7</i>(4), 547-556.Lonsdorf, E. V., Rundlöf, M., Nicholson, C. C., & Williams, N. M. (2024). A spatially explicit model of landscape pesticide exposure to bees: Development, exploration, and evaluation. <i>Science of the Total Environment</i>, <i>908</i>, 168146.McGarigal, K., & Marks, B. J. (1995). Spatial pattern analysis program for quantifying landscape structure. <i>Gen. Tech. Rep. PNW-GTR-351. US Department of Agriculture, Forest Service, Pacific Northwest Research Station</i>, 1-122.Nicholson, C. C., Lonsdorf, E. V., Andersson, G. K. S., Knapp, J., Svensson, G. P., Gönczi, M., Jonsson, O., de Miranda, J. R., Williams, N. M., & Rundlöf, M. (2024). Landscapes of risk: a comparative analysis of landscape metrics for the ecotoxicological assessment of pesticide risk to bees. <i>Journal of Applied Ecology. </i>doi: TBD