Spatial overlap between predator and prey is a prerequisite for predation, but the degree of overlap is not necessarily proportional to prey consumption. This is because many of the behavioural processes that precede ingestion are non-linear and depend on local prey densities. In aquatic environments, predators and prey distribute not only across a surface, but also vertically in the water column, adding another dimension to the interaction. Integrating and simplifying behavioural processes across space and time can lead to systematic biases in our inference about interaction strength. To recognise situations when this may occur, we must first understand processes underlying variation in prey consumption by individuals. Here we analysed the diet of a major predator in the Barents Sea, the Atlantic cod (Gadus morhua), aiming to understand drivers of variation in cod’s feeding on its main prey capelin (Mallotus villosus). Cod and capelin only partly share habitats, as cod mainly reside near the seafloor and capelin inhabit the free water masses. We used data on stomach contents from ~2000 cod individuals and their surrounding environment collected over 12 years, testing hypotheses on biological and physical drivers of variation in cod’s consumption of capelin, using Generalized Additive Models. Specifically, effects of capelin abundance, capelin depth distribution, bottom depth, and cod abundance on capelin consumption were evaluated at a resolution scale of 2 km. We found no indication of food competition as cod abundance had no effect on capelin consumption. Capelin abundance had small effects on consumption, while capelin depth distribution was important. Cod fed more intensively on capelin when capelin came close to the seafloor, especially at shallow banks and bank edges. Spatial overlap as an indicator for interaction strength needs to be evaluated in three dimensions instead of the conventional two when species are partly separated in the water column.

The data was collected by Norwegian vessels participating in the Norwegian-Russian Barents Sea Ecosystem Survey in 2004-2015. The survey covers the entire ice-free Barents Sea in August-October each year and follows a regular grid design with sampling stations spaced approximately 65 km apart, collecting data on environmental conditions, species composition and abundance for several trophic levels. We selected stations from areas where cod and capelin overlap spatially (Fall, et al. 2018), i.e., stations were both species were observed. At each station, a CTD (Conductivity-Temperature-Depth) probe is lowered to measure depth-specific temperature, a Campelen 1800 demersal shrimp trawl is used for near-bottom sampling over a distance of 1.4 km (0.75 nautical miles, nmi), and a pelagic trawl (‘Harstad trawl’, Godø, et al. 1993) samples the upper approximately 60 m of the water column. Continuously during the survey, Simrad EK60 echo sounders with 18, 38, 120 and 200 kHz split beam transducers (on some vessels Simrad EK500 during the first years) record fish echoes along the survey tracks. The acoustic backscatter at 38 kHz is manually allocated to target groups based on species-specific acoustic signatures and the catch composition in pelagic and bottom trawls, then stored at a horizontal resolution of 1.9 km (1 nmi) and a vertical resolution of 10 m. One individual cod from each 5-cm length group present in the bottom trawl catch is randomly chosen for sampling of age, sex, mass, length, maturity stage, and stomach contents. The stomachs are frozen on board, and the contents identified to the lowest possible taxonomic level in a laboratory on land. Prey items are assigned a qualitative digestion stage from 1-5, where 1 corresponds to undigested prey and 5 to prey that is so heavily digested that it is unidentifiable by visual inspection. Prey categories are weighed (wet mass) and, if possible, counted and length measured. The dataset contains stomach data from cod in the size range 30-70 cm. Empty stomachs (9.2%) were excluded from the analyses, leaving 1944 stomachs from 455 stations across the 12 years.

The following covariates are included in the dataset: capelin density, capelin weighted median depth, cod density, cod length (for each individual within the station), and bottom depth (measured by the vessel-mounted echosounder at the start of trawling). The local capelin density (in units of Nautical Area Scattering Coefficient; m²/nmi², integrated throughout the water column) was taken from the acoustic segment of 1.9 km length that had the highest temporal overlap with each trawl haul. This included acoustic recordings from before, during and after trawling. For a more intuitive representation of prey density, we converted the acoustic values to number of individuals/km² based on the length distribution in the closest pelagic trawl hauls taken during the survey, and the relationship between capelin length and acoustic target strength (Gjøsæter, et al. 1998). For capelin weighted depth, we used the same acoustic segment resolved in 10 m depth channels to calculate the weighted median depth of the capelin backscatter at each trawl location. This depth was then standardised with respect to time of day at each trawl location by fitting a Generalised Additive model from which capelin depth was predicted for a fixed time of day (see Standardisation of capelin vertical distribution in the methods section of the article). The local cod density (number of individuals ≥ 30 cm/km²) was estimated from each demersal trawl haul using standard methods for cod swept area calculation in the Barents Sea, which assume that the effective sweep width of the trawl increases with cod size up to 62 cm (Mehl, et al. 2014).

Details on the dataset and variables

Date of data collection

2004-08-16 -- 2015-09-18

Geographic location of data collection

Barents Sea; longitude 18.41167 -- 46.57667 E, latitude 72.33167 -- 80.10333 N (decimal degrees)

Number of variables

42. Comma-separated.

Number of cases/rows

1944. Each line represents a cod individual.

Variable list

1. Name: year. Description: The year the cod was sampled. Used as a covariate in the analysis.
2. Name: date. Description: The date the cod was sampled.
3. Name: sampling.day. Description: The day of the year the cod was sampled, in relation to the first day of the year that was sampled in the dataset (i.e., August 8^th is day 0). Used as a covariate in the analysis.
4. Name: station.id. Description: A unique ID for the sampling station. Several cod (lines in the dataset) may be sampled at each station. Used as a covariate in the analysis.
5. Name: start_time. Description: start date and time for the trawl haul in which the cod was sampled. On the form year-month-day hour:min:sec.
6. Name: stop_time. Description: stop date and time for the trawl haul in which the cod was sampled. On the form year-month-day hour:min:sec.
7. Name: lon. Description: longitude of sampling station (decimal degrees).
8. Name: lat. Description: latitude of sampling station (decimal degrees).
9. Name: x. Description: easting of sampling station, expressed as distance in km from 35 degrees E (approximate center of the Barents Sea). Used as a covariate in the analysis.
10. Name: y. Description: northing of sampling station, expressed as distance in km from 75 degrees N (approximate center of the Barents Sea). Used as a covariate in the analysis.
11. Name: bottomdepthstart. Description: Bottom depth at the start of trawling.
12. Name: sunheight. Description: Angle of the sun relative to the horizon (angle = 0) at the time and location of sampling. Calculated from the day of the year, geographical location, and time of day.
13. Name: sunrise. Description: Angle of the sun at sunrise at the sampling location. Calculated from the day of the year and geographical location.
14. Name: corr.sun. Description: Solar elevation angle (sunheight) adjusted for the angle of sunrise, which changes during the survey. corr.sun = sunheight-sunrise.
15. Name: bottom_temp. Description: Temperature (degrees Celsius) measured near the seafloor in the CTD cast taken closest to the trawl station.
16. Name: dist.CTD.nmi. Description: distance from the trawl station to the nearest CTD cast, in nautical miles.
17. Name: cod.tot. Description: Swept area density of cod < 30 cm in the trawl haul (individuals/km²).
18. Name: best.ak. Description: Start date and time of the acoustic segment with which the trawl haul had the highest temporal overlap.
19. Name: capelin.sa. Description: acoustic backscatter (Nautical Area Scattering Coefficient, NASC; m²/nmi²) of capelin from the acoustic segment with highest temporal overlap with the trawl station (best.ak). Integrated in the water column.
20. Name: SA.per.ind. Description: Acoustic backscatter per capelin individual estimated from the capelin length distribution in the closest pelagic trawl haul and the relationship between length and acoustic target strength.
21. Name: cap.ind. Description: Capelin density in individuals/km² (capelin.sa divided by sa.per.ind). Integrated in the water column.
22. Name: capdepth.med. Description: Weighted median depth of the capelin acoustic backscatter, calculated from backscatter resolved in 10 m depth channels.
23. Name: corr.capdepth. Description: capelin weighted median depth adjusted so that it is not deeper than the bottom depth at the start of trawling. This can happen since the acoustic transect covers 1.9 km, over which the bottom depth may change.
24. Name: pred.capdepth. Description: Predicted weighted median depth of capelin acoustic backscatter at median solar elevation angle and survey day. See main article for details.
25. Name: se.pred.capdepth. Description: Standard error of the capelin depth prediction.
26. Name: fish.id. Description: Unique identifier for each cod individual.
27. Name: ind_weight_kg: Description: Weight of the cod individual in kilograms.
28. Name: ind_length_m: Description: Length of the cod individual in meters.
29. Name: age: Description: age of the cod individual in years estimated from otolith reading.
30. Name: sex: Description: Sex of the cod individual. 1: female, 2: male, NA: not determined.
31. Name: tot.weight.prey. Description: Total weight (wet mass, kilograms) of all prey items found in the cod stomach.
32. Name: n.prey. Description: Number of prey items found in the cod stomach.
33. Name: totcap. Description: Total weight (wet mass, kilograms) of capelin prey found in stomach.
34. Name: which.cap. Description: Whether the cod individual has fed on capelin only, other prey only, or capelin and other prey.
35. Name: wdigcap. Description: Mean digestion stage of capelin prey weighted by prey mass (in cases where the cod had capelin of different digestion stages in its stomach). The digestion stage ranges from 1 to 5, where 1 is newly eaten prey and 5 is prey that is completely digested and usually cannot be identified. Heavily digested capelin prey may sometimes be identified based on its particular smell.
36. Name: capelin.consumption.rate. Description: Estimated consumption rate of capelin based on a stomach evacuation model fitted to the mass capelin in the stomach (totcap), bottom temperature (bottom_temp), cod mass (ind_weight_kg*1000) and an evacuation constant specific to capelin prey (see main article).
37. Name: capelin.consump.scaled. Description: Capelin consumption rate per g of cod (capelin.consumption.rate divided by ind_weight_kg*1000). Response variable in the analysis.
38. Name: st.capind. Description: Capelin density (cap.ind) standardized to z-scores. Used as a covariate in the analysis.
39. Name: st.pred.capdepth. Description: Predicted capelin weighted median depth (pred.capdepth) standardized to z-scores. Used as a covariate in the analysis.
40. Name: st.codden. Description: Cod density (cod.tot) standardized to z-scores. Used as a covariate in the analysis.
41. Name: st.codlength. Description: Cod length (ind_length_m) standardized to z-scores. Used as a covariate in the analysis.
42. Name: st.depth. Description: Bottom depth (bottomdepthstart) standardized to z-scores. Used as a covariate in the analysis.