Leveraging multiple data sources to assess competition between introduced wild pigs and native deer
Data files
Oct 01, 2025 version files 697.16 MB
-
Pig_and_Deer_Competition.zip
697.14 MB
-
README.md
20.64 KB
Abstract
Introduced species have diverse impacts on native wildlife and ecosystems. Negative effects of introduced species on native species through competition or amensalism are frequently hypothesized, but challenging to test in the field, particularly for large vertebrates. Wild pigs (Sus scrofa) are a widely distributed introduced species that are hypothesized to compete with native wildlife. We used multiple data sources to evaluate the strength of negative effects of wild pigs on native black-tailed deer (Odocoileus hemionus columbianus) by looking at multiple niche dimensions and fitness consequences for deer. We combined isotopic niche data with camera surveys that provided density, spatial overlap, and temporal overlap of deer and wild pigs. We also used cameras to gather data on deer body condition and reproduction to use as fitness metrics. We evaluated niche overlap between deer and wild pigs and evaluated if deer fitness was related to the density of wild pigs. Our results illustrate overlap in space, time, and isotopic niche between wild pigs and deer. Despite this overlap, we observed no relationship between wild pig density and deer body condition or reproduction, suggesting minimal effects of competition or amensalism on deer. Our results confirm the resource overlap that has been observed between deer and wild pigs in other systems, but contradict the perception that wild pigs are negatively affecting deer through resource competition. We demonstrate that the effects of wild pigs on deer may be weaker than has previously been thought, depending on the resource availability. Our approach also provides a framework for going beyond measuring niche overlap to provide a more thorough understanding of competition between large vertebrates in a field setting.
Metadata for the manuscript 'Leveraging multiple data sources to assess competition between introduced wild pigs and native deer' published in Ecosphere by MA Parsons and JK Young.
The data and R code provided here can be used to recreate all results from this manuscript. The data files included several raster layers of topographic and landcover variables, photo data from a camera grid, cougar kill site information, and isotope data for deer and pigs.
The R scripts conduct all analyses including estimating density of deer and pigs, calculating dietary overlap via isotopes, and estimating what landscape variables are related to deer body condition and reproduction.
DATA FILES
ProcessedData
The datafiles include:
- Raster layers of landcover and topographic data
- Camera station and photo data from a camera grid used to monitor deer and pigs
- Cougar kill site locations and data
- Isotope data
GIS Data
Vegetation
This raster layer is a 30x30m resolution raster of the log transformed distance to agricultural cover. This layer was derived from the 2021 National Landcover Database.
This raster layer is a 30x30m resolution raster that is 1/0 indicating whether a cell occurs on land or over the ocean. This was used for mapping and masking other data to avoid calculating animal density at sea. This was derived from the 2021 National Landcover Database by summing all terrestrial landcover types and reclassifying any terrestrial landcover as 1.
This raster layer is a 30x30m resolution raster of the log transformed distance to forest cover. This layer was derived from the 2021 National Landcover Database (https://www.mrlc.gov/data; classes 41, 42, and 43).
This raster layer is a 30x30m resolution raster of the log transformed distance to herbaceous cover. This layer was derived from the 2021 National Landcover Database (classes 71, 72, 73, 74, 81, 82).
This raster layer is a 30x30m resolution raster of the log transformed distance to riparian cover. This layer was derived from the 2021 National Landcover Database (classes 11, 12, 90, 95).
This raster layer is a 30x30m resolution raster of the log transformed distance to shrub cover. This layer was derived from the 2021 National Landcover Database (classes 51, 52).
These rasters are 30mx30m resolution rasters of the average seasonal NDVI value for four seasons: summer (May-July), fall (August-October), winter (November-January), and spring (February-April). These data are derived from MODIS 16-day data at the 250m scale.
This raster layer is a 30x30m resolution raster of the proportion of forest cover within a 17x17 cell moving window (510m) derived from the 2021 National Landcover Database.
This raster layer is a 30x30m resolution raster of the proportion of herbaceous cover within a 17x17 cell moving window (510m) derived from the 2021 National Landcover Database.
This raster layer is a 30x30m resolution raster of the proportion of riparian cover within a 17x17 cell moving window (510m) derived from the 2021 National Landcover Database.
This raster layer is a 30x30m resolution raster of the proportion of shrub cover within a 17x17 cell moving window (510m) derived from the 2021 National Landcover Database.
Topographic
This raster layer is a 30x30m resolution raster of elevation. This layer was derived from the USGS 1/3 arc second digital elevation model (https://apps.nationalmap.gov/downloader/)
This raster layer is a 30x30m resolution raster of slope calculated from resampled_DEM.tiff using the terra package in R.
This raster layer is a 30x30m resolution raster of the terrain ruggedness index calculated from resampled_DEM.tiff using the terra package in R.
This raster layer is a 30x30m resolution raster of the topographic position index calculated from resampled_DEM.tiff using the terra package in R.
This raster layer is a 30x30m resolution raster of the terrain ruggedness index averaged over a 17x17 cell moving window (510 m).
This raster layer is a 30x30m resolution raster of the topographic position index averaged over a 17x17 cell moving window (510m).
Animal Density Layers
These tiff files are 30x30m rasters of predicted deer and pig density for the study area based on the selected random encounter and staying time models developed in the manuscript. There are three files for each species:
- *_average.tif is the estimated species density ignoring season-specific covariates in the model.
- *_Summer.tif is the estimated species density incorporating summer-specific covariates in the model.
- *_Winter.tif is the estimated species density incorporating winter-specific covariates in the model.
These tiff files are 30x30m rasters of predicted deer and pig density averaged over a 17x17cell moving window (510m). These are derived from the DeerDensityPred*.tif and PigDensityPred*.tif files. There are three files for each species:
- *_average_289cell.tif is the estimated species density ignoring season-specific covariates in the model.
- *_summer_289cell.tif is the estimated species density incorporating summer-specific covariates in the model.
- *_winter_289cell.tif is the estimated species density incorporating winter-specific covariates in the model.
Other GIS Layers
This raster layer is a 30mx30m resolution raster that defines the subset of areas withing the study area that are closed to hunting.
Shapefiles
Shapefiles are a common format for vector-based geographic information system (GIS) data. They can be opened and used in any GIS software and in R or Python. A shapefile consists of multiple file types beyond the .shp (specifically, .cpg, .dbf, .prj, .sbn, and .sbx). The user only interacts directly with the .shp file but the other files need to be in the same directory.
These three shapefiles are the boundary and rivers of the Fort Hunter Liggett and are used for mapping and figure making.
Raw Data
Camera Data
This csv file contains the camera station locations and deployment dates. Each row represents a camera deployment and columns provide information on the deployment. There are 26 columns:
- Station ID - the name of the camera location, each location was sampled multiple times
- CSY - the camera-season-year combination of the deployment
- CSY2 - the camera-season-year combination plus a unique identifier when there were multiple deployments at a location within a season. This occurred when cameras were moved or replace because of technical issues. This is a unique identifier for each deployment.
- UTM_X - the UTM easting of the camera deployment
- UTM_Y - the UTM northing of the camera deployment
- Habitat - a broad habitat category for the camera location
- Distance - the distance of the 100% detection zone in meters. All detection zones used a 30 degree field of view
- Set - was this location part of the "even" set of cameras or "odd" set of cameras during our camera rotations
- Season - season of the camera deployment
- Year - year of the camera deployment
- Set Date - the date the camera was deployed
- Pull - the date the camera was removed from the location
- Deployment_days - the total length of the deployment
- Camera Type - The type of camera that was used: Browning (Brown), Bushnell (Bush), Reconyx. Old (O), New (N), UCD (from project partner)
- Problem*_From - The start date of a problem with the camera that prevented detection (e.g., card filled, batteries died, animal damaged/displaced). Max of 4 problems for a single deployment. NAs indicate lack of a problem period
- Problem*_to - The end date of a problem with the camera. NAs indicate lack of a problem period
- Problem*_days - The number of days the problem persisted
This csv file contains the information of visits to cameras by deer and pigs and is used to run the Random Encounter Staying Time model. Each row represents a visit by a single animal to the camera. There are 10 columns:
- CSY - the camera-season-year combination of the camera deployment where the visit occurred
- CameraID - the station ID of the camera where the visit occurred
- Date - the date of the visit
- Time - the time the visit began
- season - the season of the visit
- year - the year of the visit
- Species - the species for each visit. either deer or pig
- Start - the start time of the visit in seconds (origin 1/1/1970)
- Stop - the stop time of the visit in seconds (origin 1/1/1970)
- visit_length - the duration of the visit in seconds
This csv file contains information on deer sex, age, and body condition. Each row is data from a single photograph and there are 21 columns:
- File - the name of the photo file
- Camera - the camera station where the photo was taken
- Season - the season when the photo was taken
- Year - the year when the photo was taken
- RelativePath - the relative path to access the photo
- Folder - the working folder during photo processing
- Date - the date the photo was taken
- Time - the time the photo was taken
- ImageQuality - autogenerated value of image quality produced by the software
- DeleteFlag - autogenerated value by the software used for processing photos
- Doe_count - the number of adult female deer in the photo
- Buck_count - the number of adult female deer in the photo
- Fawn_count - the number of deer fawns in the photo
- Comments - any comments about the photo. NAs indicate no comments were recorded.
- Animal_*_Condition - the body condition score (0-4) of each individual adult deer in the photo. Followed methods from Smiley et al. (2020) https://doi.org/10.1002/wsb.1070. NAs indicate that there was no animal to score (i.e., there were fewer individuals than the column suggests (e.g., NA in Animal_2_Condition means there was only one animal, NA in Animal_5_Condition means there were fewer than 5 animals))
- Confidence_Condition - the confidence score (0-3) for the body condition score of the deer in the photo
Cougar kill site information
We are unable to share these raw data because of their sensitive nature. We uploaded a censored version of this data that excludes GPS locations (cougar_clusters_noGPS.csv)
This csv file includes information on cougar kill site investigations. Each row represents a cluster of cougar GPS points that was investigated because it met a criteria to be a possible kill site. Each column provides information on the cluster. There are 62 columns.
- cluster_id - the name of the GPS cluster investigated
- cougar_id - the ID of the individual cougar associated with the cluster
- type - what type of cluster was investigated. The majority are "current" clusters indicated that they were investigated soon after formation. Other options include April (clusters investigated long after formation), day (presumed daybed locations), Den (confirmed den locations), and trapping (area where we deployed bait to trap animals). Only 'current' clusters were used in analysis.
- collar_serial - the serial number of the cougar's GPS collar
- n_points - the number of points in the GPS cluster
- radius - the radius of the GPS cluster
- duration_h - the duration between the first and last point at the GPS cluster
- form_date - the date of the first point at the GPS cluster
- centroid_lat - the latitude of the centroid of the GPS cluster
- centroid lon - the longitude of the centroid of the GPS cluster
- inv_date - the date we investigated the GPS cluster
- personnel - initials of the investigators
- time_200m - the time of day we reached 200m from the cluster centroid
- habitat - broad habitat category that the cluster was in
- habitat2 - broad secondary habitat category that the cluster was in. NAs indicate that there was only one dominant habitat type (habitat column)
- canopy - whether the site had open, moderate, or closed canopy cover
- carcass_found - whether or not a carcass was found at the site
- time_found - the time of day the carcass was found. NA indicates that no carcass was found
- carcass_lat - the latitude the carcass was located at. NA indicates that no carcass was found
- carcass_lon - the latitude the carcass was located at. NA indicates that no carcass was found
- species - the species of the located carcass. NA indicates that no carcass was found
- how_id - what features were used to identify the carcass. NA indicates that no carcass was found
- sex - the sex of the carcass. F = female, M = male, UNK = unknown. NA indicates that no carcass was found
- how_sex - what features were used to determine the sex of the carcass. NA indicates that no carcass was found
- age - the age category of the animal that was killed. Neonate (<6 mo), Fawn/Calf (6 - 18 mo), yearling (18 - 30 mo), adult (>30 mo). NA indicates that no carcass was found
- how_age - what features were used to determine the age of the carcass. NA indicates that no carcass was found
- est_age - the estimated age of the carcass in months/years. NA indicates that no carcass was found
- cacehd - whether or not the carcass was cached currently, previously, or not at all. NA indicates that no carcass was found
- n_prey - the number of prey items at the cluster. NA indicates that no carcass was found
- n_piles - the number of separate piles of remains and evidence. NA indicates that no carcass was found
- condition - the condition of the carcass (intact, dismantled, dispersed). NA indicates that no carcass was found
- parts - a description of what parts of the carcass were located. NA indicates that no carcass was found
- bite_wounds - whether bite wounds were located. NA indicates that no carcass was found
- hemorrhaging - whether hemorrhaging was located and if so the location of hemorrhaging. NA indicates that no carcass was found
- fed_on - whether the carcass had been fed on. NA indicates that no carcass was found
- per_rem - the rough percentage of biomass remaining on the carcass. NA indicates that no carcass was found
- feeding_use - a description of the observed feeding sign on the carcass. NA indicates that no carcass was found
- marrow - the color and texture of the bone marrow from long bones. NA indicates that no carcass was found or marrow was not assessed
- drag_trail - whether a drag trail was located. NA indicates that no carcass was found
- length - the length of the drag trail. NA indicates that no drag trail was found
- kill_site - whether a kill site was located. NA indicates that no carcass was found
- kill_lat - the latitude of the kill site. NA indicates that no kill site was found
- kill_lon - the longitude of the kill site. NA indicates that no kill was found
- evidence_kill_site - the evidence used to identify the kill site. NA indicates that no carcass was found
- consistent_date - whether the condition of the carcass is consistent with the cluster formation date. NA indicates that no carcass was found
- confidence_cougar_kill - confidence level in whether a cougar killed or scavenged the prey. positive, probable, unknown, or scavenge. NA indicates that no carcass was found
- other_sign - sign of other species that was observed at the carcass. NA indicates that no carcass was found
- scavenged - whether or not the carcass was scavenged by other species. Y, N, or Likely. NA indicates that no carcass was found
- scavened_by - what other species scavenged the carcass. NA indicates that no evidence of scavenging was found
- kill_2 - whether there was a second kill at the cluster. NA indicates that no carcass was found
- species - the species of the second prey item. NA indicates that no second carcass was found
- cached - whether the second prey item was cached. NA indicates that no second carcass was found
- notes - any notes about the cluster investigation. NA indicates that no notes were recorded
- cache_photo - whether photos were taken of the cache
- carcass_photo - whether photos were taken of the carcass
- bite_photo - whether photos were taken of the bite wounds
- hem_photo - whether photos were taken of the hemorrhaging
- drag_photo - whether photos were taken of the drag trail
- kill_photo - whether photos were taken of the kill site
- coug_sign_photo - whether photos were taken of cougar sign
- scav_sign_photo - whether photos were taken of scavenger sign
- time_leave - the time of day we left the GPS cluster
Isotopes
This csv file contains the raw isotope data for deer and pig hair samples. Each row represents an individual hair sample. There are 13 columns.
- Sample ID - the ID of the hair sample
- Species - the species of the sample. pig = SUSC and deer = ODHE
- Tray # - the tray number the sample was in for analysis
- TrayID - the samples location within the tray
- Weight (mg) - the mass of the sample in mg
- d15N - the nitrogen isotope value in parts per mille
- d13C - the carbon isotope value in parts per mille
- %N (wt.) - the nitrogen percent mass of the sample
- %C (wt.) - the carbon percent mass of the sample
- C:N - the carbon to nitrogen ratio
- ClusterID - the location each sample came from. Samples that begin with F/M were collected from live captured pigs. All others were collected from cougar-killed prey
- Season - the season each sample was collected
- Year - the year each sample was collected
CODE FILES
- Extracting camera covariates. Also used later to make layers of predicted deer and pig density
- Run the random encounter staying time model for deer and pigs. Separate file for each species
- Analyze deer body condition and fawn to doe ratios to evaluate if there is evidence of competition with pigs
- Analyze the isotope data
This file uses the camera station data and GIS data to extract covariates for each camera location
It creates and saves a new file (CameraCovariates_NDVI_289_log.csv) that is needed by the random encounter staying time model
Then, starting on Line 174, this file creates rasters of predicted deer and pig density needed for the deer body condition models
The second half of this script cannot be run until the REST model has been run and output saved.
These files use the camera station data, camera covariate data (created from previous script), and deer and pig visit data.
It processes data and runs the random encounter staying time model for the species.
It then runs a cross validation of the model to test model performance
Note that these models can take several hours to run depending on the number of iterations, chains, and computer specs
The beginning of this code is also used to estimate deer and pig temporal overlap
This file test is deer fawn to doe ratios or body condition are related to pig density
It uses the deer body condition data from cameras, GIS data, the camera station data, and kill site data (not necessary for the manuscript)
It processes the data and incorporates covariates from the GIS data
It then runs generalized linear models to test if deer fawn to doe ratios or body condition are related to landscape covariates, seasonal covariates, or estimated pig density
This file estimates the isotopic niche overlap of deer and pigs
It uses the raw isotope data, creates visualizations of the overlap, and estimates overlap using Bayesian standard ellipses for each species
