Study System
This study was conducted in two different fjords in western Norway, Nordfjord and Hardangerfjord (Figure 1). With lengths of approximately 110 and 160 km respectively, Nordfjord and Hardangerfjord are among the longest fjords in Norway. Each provides a relatively complex migration route for the post-smolts to navigate, with widths ranging from 2-7 km and many branching fjord arms and inlets. Both fjords are in areas of Norway where the sea lice-induced mortality is estimated to be high (20).
Smolt individuals from two rivers in each fjord were tagged in both 2018 and 2019: Stryn and Eid in Nordfjord, and Eio and Granvin in Hardangerfjord. All of these rivers have their outlets in the inner fjord, such that smolts emigrating from these rivers must swim more than 50-120 km to reach the open sea. Similarly, each of these rivers have comparable hydrodynamics as they are all relatively short and drain from a large lake surrounded by mountains.
Due to the influx of freshwater from precipitation and snowmelt from the surrounding mountains during springtime, a brackish layer of highly variable extent is created in the upper water column of the fjords. Interannual and geographical variation in timing and magnitude of freshwater input, alongside variable hydrographic conditions in the marine environment, render this variation large at the relevant timescales of post-smolt migration. As smolts from these rivers generally do not begin their migration before the spring snowmelt has begun (21), post-smolts can use this brackish layer to graduate their acclimation to seawater (see Supplementary Materials).
In Hardangerfjord and neighboring Bjørnafjord, a total of 106 Thelma Biotel receivers (TBR700) were deployed. In Nordfjord, 71 Innovasea (VR2W) receivers were deployed (Figure 1). All receivers were moored and attached to ropes with the hydrophones oriented in a downward position, with buoys ensuring that the receivers were positioned 3-4 meters below the surface. Within each fjord, the receivers were grouped into zones A-D such that zone A consisted of the area surrounding each estuary, zone D consisted of the outer fjord near to the open sea, and zones B and C were intermediary zones between zone A and D (Figure 1). These are the same groupings used for the survival analysis in Bjerck et al., 2021.
Fish Sampling and Tagging
In April of each year, pre-smolts were captured using DC electrofishing (Ing. Paulsen, Norway, FA4 ,1600V, 80Hz) in each river. The fish were kept in 60 L holding tanks with flow of river water for 24 hours prior to tagging. These fish were then tagged with ThelmaBiotel acoustic tags (D-LP7), using the procedure described in Bjerck et al. (2021). These tags have a weight in air of 2.0 grams and a length of 21.5 millimeters and are designed to transmit both the ID of the tag and the depth of the tag with a resolution of 0.2 meters every 30-90 seconds. The pressure sensors within these tags have much higher internal resolution, but the acoustic protocol restricts the tags to transmit data as a single byte in order to minimize the amount of information being sent per transmission. This means that the maximum depth they are able to record is 51 meters. The tags transmitted signals with signal strength 139 dB Re 1 μPa @ 1 m and they were set to automatically deactivate after 155-200 days. This signal strength corresponds with a detection range on the order of 200 meters, though this is known to vary substantially through time and space with changing ocean conditions (22,23). The average length of smolts tagged with depth transmitters was 14.7 cm (SD= 1.1 cm). The tagging protocol was approved by Norwegian Authorities for animal welfare (FOTS IDs: 12002 and 15471).
Quality Control
As these tags measure depth by measuring water pressure, air pressure was controlled for by retrieving hourly weather data from the Norwegian Meteorological Institute (seklima.met.no). In Nordfjord, weather data from the Sandane Airport (SN58100) was used and, in Hardanger, weather data from Kvamsøy (SN50070) was used. There was little variation between these two stations despite the 160 kilometers separating them such that greater spatial resolution was not considered necessary. Further, variation in the calibration of the depth sensor was corrected for by retrieving the factory test value of each tag from Thelma Biotel and the air pressure from the closest weather station to the factory (Selbu II, SN68290).
The spatio-temporal migration trajectories of each individual were inspected visually in order to identify false detections and mortalities/tag losses. Detections occurring at unlikely or impossible locations in relation to the rest of an individual’s trajectory were removed. Mortality was identified based on depth data, as predated fish would often first exhibit erratic depth movements and then stop at a constant depth (or varying according to tidal cycles). Efforts were made to remove detections occurring after mortality or tag loss from the analysis.
Further, only smolts that were detected as successful migrants (i.e., smolts that were detected in the outer reaches of the fjord (zone D in Figure 1)) and which had at least 10 detections in the fjord were included in the analysis. This was a conservative approach to ensure that the data used in the quantitative analyses actually represented migrating salmon smolts rather than movements of predators (see e.g. Daniels et al., 2019). We assumed that it is unlikely that a predator of a tagged smolt will continue to exhibit a migratory trajectory similar to a post-smolt with the tag within its stomach. The use of acoustic tags designed to detect predation has revealed that this can happen within freshwater (Lennox et al. 2021), but this has yet to be documented within the fjord environment. Acoustic tagging of brown trout (Salmo trutta) from the rivers of Granvin and Eio, one of the primary potential predators of migrating smolts in these populations, showed that the majority of those that migrated did not venture out in the outer fjord (25), with similar results in Sognefjorden, a fjord in between the fjords studied here (26).
Generalized Linear Mixed Effects Modelling
In order to investigate to what extent variables of interest accounted for variation in depth use, a generalized linear mixed effects model with log-link as the link function was fitted to the data with ID as a random intercept effect using the function glmer in the R package lme4 (27). Candidate models reflecting hypotheses pertinent to the study objectives were subjected to model selection by using the Akaike Information Criterion (AIC) aiming at finding the model(s) that most efficiently explained the variance in depth use. Models attaining DAIC-values <2 were considered to have substantial support in the data (28). Variables fitted in the model included fjord zone, fjord, river of origin, waterway distance to the river mouth of origin, day of year, the position of the sun with respect to the horizon, and a parameterized version of the position of the sun.
In order to test if depth use was correlated with light conditions (sensu Davidsen et al., 2008), the sun’s vertical position was used as a predictor. The position of the sun with respect to the horizon in degrees for a given time and position was determined through the use of the R package suncalc (29). This was then parameterized with the function:

such that the parameterized values flattened out when darkness and true daylight arrived but changed continuously during dusk and dawn (see Supplementary Materials). The parameterized values therefore reflect more of a day/night switch than the raw values which change continuously through the day and night. This approach was used as, due to the high latitude of the study system, the angle at which the sun moves with respect to the horizon is very acute, leading to prolonged dusk/dawn periods. Additionally, later in the season, true darkness becomes an impossibility at these latitudes as the sun reaches its nadir just below the horizon.