Vertical density distributions of fish: a balance

ICES Journal of Marine Science, 59: 679–710. 2002
doi:10.1006/jmsc.2002.1249, available online at http://www.idealibrary.com on
Vertical density distributions of fish: a balance between
environmental and physiological limitation
Boonchai K. Stensholt, Asgeir Aglen, Sigbjørn Mehl,
and Eivind Stensholt
Stensholt, B. K., Aglen, A., Mehl, S., and Stensholt, E. 2002. Vertical density
distributions of fish: a balance between environmental and physiological limitation. –
ICES Journal of Marine Science, 59: 679–710.
Data (trawl, acoustic, CTD) from scientific surveys along the Norwegian coast, in the
North Sea, in the Barents Sea, west of the British Isles, and in the Irminger Sea are
used. The vertical density distributions of blue whiting, cod, haddock, redfish, saithe,
capelin, and herring are described in relation to environmental conditions and
physiological limitations. The first four surveys mainly cover banks and shelf areas
shallower than 500 m. The last two surveys, aimed at blue whiting and redfish, mainly
cover shelf edge and deep-sea areas with depths from 200 to 1300 m and from 440 to
3000 m. In regard to cod some information from data-storage tags is used.
For physoclists the relative vertical profile of each acoustic sample i.e. acoustic-area
backscattering coefficient (sA), is expressed in terms of the relative pressure reduction
level from seabed up to surface. Thus relative vertical profiles with different bottom
depths are normalized and are made compatible for a discussion in terms of the free
vertical range (FVR). This restriction to rapid vertical movement is evident in the
physoclist species studied. For samples in the shelf area, the profiles show that blue
whiting, haddock, saithe, cod, and redfish are mainly distributed within the bottom
half of the water column. Some fish adapt to pelagic living especially in areas with high
acoustic density and where the bottom is deep. Here a pelagically living fish is defined
as an individual fish having a current free vertical range that does not include the
seabed.
For demersal fish, day and night relative vertical profiles are corrected for unequal
day and night losses in the bottom acoustic dead zone, which is the zone near the
seabed where echoes from fish cannot be discriminated from the sea bottom echo. Day
and night samples are separated by the sun’s passing 5 below the horizon. In most
years evidence of diurnal vertical migration is found for all investigated species. In
many cases of demersal fish there is a higher relative acoustic density (sA-values) in the
mid-range of the bottom half of the water column in the daytime as opposed to the
night-time. At night there is a degree of separation, one group of fish descends to
aggregate near the seabed and another ascends. Inter-annual variations in the diel
movement from different parts of the stock are discussed in relation to the inter-annual
variations in age composition of the stock.
2002 International Council for the Exploration of the Sea. Published by Elsevier Science Ltd.
All rights reserved.
Keywords: area backscattering coefficient (sA), acoustic bottom dead zone loss,
demersal, free vertical range, fish vertical distribution, hydrostatic pressure, pelagic,
physostomes, physoclists, swimbladder.
Received 29 August 2001; accepted 8 April 2002.
Boonchai K. Stensholt, Asgeir Aglen, and Sigbjørn Mehl: Institute of Marine Research,
PO Box 1870, N-5817, Bergen, Norway; tel: 47 55 23 8667; fax: 47 55 23 8687; e-mail:
[email protected]. Eivind Stensholt: The Norwegian School of Economics and Business
Administration, Helleveien 30, N-5045 Bergen, Norway.
Introduction
In marine teleosts with a swimbladder, that feature
occupies 3.1–5.7% of the body volume (Alexander, 1966;
Marshall, 1972). One of its functions is to be a buoyancy
1054–3139/02/080679+32 $35.00/0
organ (Blaxter and Tytler, 1978). However, its physiological limitation in adaptation to pressure changes
imposes restrictions on vertical migration (Harden
Jones, 1949, 1951, 1952; Brawn, 1962; Kutchai and
Steen, 1971; Blaxter and Tytler, 1972, 1978; Ross, 1979;
2002 International Council for the Exploration of the Sea. Published by Elsevier Science Ltd. All rights reserved.
680
B. K. Stensholt et al.
Harden Jones and Scholes, 1981, 1985; Arnold and
Greer Walker, 1992). There are metabolic costs of
swimbladder adaptation and of compensatory swimming when the buoyancy is not neutral (Alexander,
1970, 1972; Blaxter and Tytler, 1978).
During pressure reduction fish with open swimbladders (physostomes, e.g. herring and capelin), may
release gas through a pneumatic duct. Gas secretion is
known to be slow and insignificant in physostomes.
Brawn’s experiment indicates that herring take in air at
the surface. The amount of air determines the level of
neutral buoyancy, which may be no more than 50 to
100 m. In experiments physostomes respond to rapid
pressure increases by ascending (Brawn, 1962; Blaxter
and Tytler, 1978; Sundnes and Bratland, 1972).
Fish with closed swimbladders (physoclists, e.g. blue
whiting, cod, haddock, redfish and saithe) produce gas
by secretion in response to under-buoyancy and remove
gas by resorption in response to over-buoyancy but
these processes take time (Blaxter and Tytler, 1978;
Harden Jones and Scholes, 1985). At any point in time
the gas content determines a depth level of neutral
buoyancy. For rapid vertical movement around that
depth is a zone of free vertical range (FVR). Within this
range the fish can swim freely and compensate for the
deviation from neutral buoyancy by means of swimming
movements (Harden Jones, 1949, 1951, 1952; Harden
Jones and Scholes, 1981, 1985; Tytler and Blaxter, 1973,
1977; Blaxter and Tytler, 1972, 1978; Arnold and Greer
Walker, 1992). By ascending rapidly too far above the
neutral buoyancy plane, a physoclist gets too overbuoyant, risking an uncontrollable rise to the surface or
a swimbladder rupture (Harden Jones, 1951, 1952;
Tytler and Blaxter, 1973; Harden Jones and Scholes,
1985), while physostomes get rid of surplus gas
through their duct (Brawn, 1962). On the other hand,
during deep dives physoclists have an advantage over
physostomes because their gas secretion reduces underbuoyancy. Thus physoclists and physostomes are
naturally disposed towards different patterns of vertical
distribution and vertical migration.
Given the current gas content of the swimbladder, the
current free vertical range (FVR) allows rapid vertical
movement between two planes, at M and m metres, the
limiting depth levels for under- and over-buoyancy. As
one atmospheric pressure is about the pressure of a 10 m
water column, the pressure ratio (m+10)/(M+10) can be
considered a ‘‘theoretical’’ physiological limit to pressure
changes within a free vertical range. Experiments
(Harden Jones and Scholes, 1985) have determined that
the ratio is ca. 50% in the case of cod. Blaxter and Tytler
(1972) concluded that the sensitivity to pressure decrease
was about the same for cod and saithe, and that haddock
was slightly more sensitive. We expect that saithe,
haddock, redfish and blue whiting have a free vertical
range which is not very different from that of cod.
Several studies conclude that a physoclist, which
undertakes extensive diel vertical migrations, generally
maintains neutral buoyancy only near the top of its daily
vertical range and otherwise is negatively buoyant.
Descending some distance below the current free vertical
range has no serious consequences however. It may
adapt at the bottom by gas secretion and, perhaps with
an extra effort, ascend into its FVR (Alexander, 1966;
Harden Jones and Scholes, 1981, 1985; Arnold and
Greer Walker, 1992; Rose and Porter, 1996; Righton
et al., 2001).
Because of gas loss by passive diffusion due to
permeability of the swimbladder wall, there may be a
maximal depth where neutral buoyancy maintenance is
possible (Kutchai and Steen, 1971; Ross, 1976, 1979;
Blaxter and Tytler, 1978).
Temperature and salinity impose environmental
barriers to fish distribution. In cod a physiological
adaptation process starts when it enters waters colder
than 2C that enables it to live in colder temperatures. It
also tolerates cold shocks from 4C to 0C in normal
seawater. High temperatures (above 10C) with low
salinity may be lethal for large cod (Woodhead and
Woodhead 1959, 1965; Harden Jones and Scholes,
1974). Björnsson et al. (2001) found that large Icelandic
cod (d2 kg) cannot tolerate prolonged exposure to 16C
with salinity 32‰. The optimal temperature for cod’s
growth found in laboratory experiments is higher than
its natural ambient temperature (Björnsson et al., 2001).
In the natural environment preferred temperature is just
one of many concerns, like availability of food and
avoidance of predation (Kristiansen et al., 2001).
Time series of depth and temperature from data
storage tags (DST) attached to adult cod show interplay
between the individual cod’s migration behaviour and
season, depth, temperature, tidal current and its physiological limitation (Stensholt, 2001; Righton et al. 2001).
Combined information from different sources helps the
interpretation of conventional fisheries survey data
(Righton et al., 2001).
In the present study information from acoustic
records, i.e. the area backscattering coefficient (sA) as
defined in MacLennan et al. (2001), trawl sampling, and
CTD, are combined to investigate the distribution of
cod, haddock, saithe, redfish, blue whiting, herring, and
capelin in relation to environmental and physiological
factors.
Besides the swimbladder physiology described above,
daylight, food availability, predator avoidance, tidal
currents, temperature, salinity, depth, and vessel noise
are among the factors that influence vertical migration
(Neilson and Perry, 1990; Ona and Godø, 1990; Godø
et al., 1999; Aglen, 1994). The observed large-scale
vertical distribution is expressed in physiological terms,
thus all samples with various depths are normalized,
and on that basis the influence of some environmental
Vertical density distributions of fish
factors is discussed. The structure of the spatial density
distribution of fish in relation to environment and
physiological limitation can be established. This can be
used in the estimation of total stock abundance.
Vertical migration behaviour affects the accuracy of
acoustic stock estimates of demersal physoclists, mainly
through the loss of acoustic fish information in the
bottom dead zone, i.e. fish very close to the bottom
cannot be distinguished from the bottom itself
(Aglen, 1994; Ona and Mitson, 1996; Aglen et al., 1999).
Moreover, the extent of losses due to negative buoyancy
near the bottom, body orientation, and avoidance
reactions to vessel passage among cod, haddock and
saithe is still unclear (Ona and Godø, 1990; Aglen, 1994;
Godø et al., 1999; Totland, pers. comm.). The accuracy
can be improved when the amount that is lost varies
according to predictable environmental conditions, e.g.
daylight and tidal current. How should one interpret
vertical profiles, which are based on the observed
acoustic abundance of the variable proportion of fish
above the dead zone? The difference between the losses
in acoustic density for day and night is estimated and the
vertical profiles are adjusted accordingly, under the
assumption that the dead zone is the main reason
for day–night differences in acoustic observations for
demersal fish during winter in the Barents Sea (Aglen
et al., 1999). In the case of cod the results are discussed
in the light of information from data-storage tags, which
provides details of the daily vertical ranges and the time
spent at each depth level.
Materials and methods
The Institute of Marine Research, Bergen, performs
several annual acoustic surveys for fish stock assessment.
In this study observations of acoustic backscattering
(sA-values) from several acoustic surveys aimed at
demersal fish and blue whiting have been used. Data
from the demersal fish surveys contain acoustic values
that are allocated to the demersal species cod, haddock,
saithe and redfish, and, in addition, values allocated to
the pelagic species blue whiting, capelin, and herring. In
the demersal fish surveys the latter are non-target species
and not completely covered.
The data considered in this paper are from annual
surveys in the main. Firstly those for saithe, haddock,
blue whiting, and herring from the saithe surveys along
the Norwegian coast from 62N to 72N during
October–November, 1992–2000 (Nedreaas, 1995;
Korsbrekke and Mehl, 2000) and from demersal fish
surveys in the North Sea from 55N to 62N during
October–November, 1992–1998 (Anon., 1992–1996) and
during August–September in 1999–2000 (Anon., 1999)
(Figure 1). Secondly those for cod, haddock, redfish,
capelin, and herring from demersal fish surveys in the
681
Barents Sea, one during February–March 1996–2002
(Mehl and Nakken, 1996; Mehl, 1997, 1998, 1999;
Aglen, et al., 2001; Aglen, 2001b, 2002) (Figure 2) and
one during July–August 1995–2000 (Aglen and Nakken,
1996; Aglen, 1999, 2001a; Michalsen et al., 2000)
(Figure 3). Thirdly from surveys of the spawning stock
of blue whiting found west of the British Isles during
March–April 1995–1996 (Monstad et al., 1995, 1996)
and finally a survey of oceanic redfish in the Irminger
Sea, southeast of Greenland, during June–July 2001
(Dalen and Nedreaas, 2001).
The saithe survey is a pure acoustic survey, which
means that the decision regarding when and where to
make trawl hauls (both demersal and pelagic) is largely
based on the acoustic observations. This survey covers
areas with highly variable topography, and the survey
design therefore varies between regions (Figure 1). The
other demersal fish surveys are combined bottom-trawl
and acoustic surveys, where there is a fixed grid of
predetermined bottom-trawl positions used both for a
swept-area estimate and for the interpretation of the
acoustic data. Some additional hauls (both demersal and
pelagic) are made for identifying acoustic recordings and
these are only used for the acoustic estimate. These
surveys mainly follow parallel transects with regular
spacing. In all the surveys CTD measurements (vertical
profiles of salinity and temperature) are made at nearly
all trawl positions. The exceptions to this are the Barents
Sea summer surveys in 1996 and 1997, and the Barents
Sea winter survey in 2000, where parts of the area were
covered by hired fishing vessels not equipped for CTD
measurements. Typical survey-line transects and trawl
stations for each of the demersal fish surveys are shown
in bathymetric maps in Figures 1–3. They cover mainly
the banks and shelf areas of depth <600 m. The blue
whiting survey is a combined pelagic trawl, bottom
trawl, and acoustic survey (Monstad et al., 1995, 1996),
which cover mainly the banks and shelf edge of depth
from 200 to 1400 m. The oceanic redfish (Sebastes
mentella) survey is a pelagic trawl and acoustic survey of
depths from 500 to 3000 m. The sA-values are split into
the redfish within and above the deep scattering layer
and, at the same time, above 500 m, and the rest of the
redfish because redfish biologists suspect that there are
two separate stocks. The sA-values within and below the
deep scattering layer are underestimated by 34% and
55% respectively according to data from a deep towed
vehicle (Magnússon et al., 1996; Dalen and Nedreaas,
2001).
In order to investigate the relationship between
fish and environmental variables satisfactorily the
respective data must be located at the same locations
over the entire studied area. The values at un-sampled
locations of temperature, salinity, sea bottom depth,
acoustic sA-values, and density of the studied species,
are estimated on the basis of observed values from
682
B. K. Stensholt et al.
69
71
67
5
65
10
63
Lo
15
61
ng
itu
de
E
20
59
25
0m
100
30
way
Nor
Norwegian Trench
20
0m
1996
Figure 1. 1996 acoustic surveys in the North Sea south of 62N, and along the Norwegian coast north of 62N. Trawl stations are
marked with . Depth isolines are in metres.
neighbouring locations by application of geostatistical
methods (Cressie, 1991; Stensholt and Sunnanå, 1996).
They are available in the software ISATIS (ISATIS,
1997), e.g. kriging with a moving window.
Acoustic data and vertical density profile
Area backscattering (sA-values) were sampled with echointegrators (Knudsen, 1990; Foote et al., 1991). During
the surveys the acoustic data were scrutinized and the
sA-values were allocated to a number of categories
(species or groups of species). For each survey series
there is a list of standard categories for the target
species and then several non-target species are defined
as separate categories usually. Cod, haddock, redfish,
herring, and capelin are defined as separate categories in
the Barents Sea surveys, while values representing saithe
are merged with several non-target species in the
category ‘‘other demersal’’. Saithe, haddock and blue
whiting are defined as separate categories in the coastal
saithe survey, the North Sea survey, and the blue
whiting survey.
During daily scrutinizing onboard the survey vessel,
the post-processing system (BEI, Knudsen, 1990)
presents the echogram data (from surface to bottom) in
blocks of 5 nautical miles (nmi) horizontal. Within such
a block the user can define a number of sub-blocks
(regular or irregular depth intervals or rectangular
boxes). Within each sub-block the operator partitions
the sA-values (averaged over 5 nmi) into the separate
categories by using the trawl catch species composition
of neighbouring stations and subjective judgment based
on experience. Thus for each acoustic sample within the
same sub-block all species have the same vertical sA
density profile (in relative terms). When the scrutinized
data are stored in the database single values are
averaged within blocks of 1 nmi (0.1 nmi for blue
whiting surveys) horizontal resolution and 10 m (25 m
for oceanic redfish) vertical resolution, relating to the sea
surface. In addition, the values from bottom up to 10 m
above bottom are stored with 1 m vertical resolution.
These profiles are used for investigating the vertical
distributions of fish in relation to environmental factors
and physiological limitations.
Cumulative relative frequency
To investigate the vertical distribution of a physostome
species or physoclists with pelagic distribution we
selected the acoustic samples with reasonably large total
Vertical density distributions of fish
683
10
15
20
25
n
e
erg
30
sb
pit
Lon
gitu
35
S
de E
40
45
50
55
80
10
0
50 0
0
20
00
Bear
B
Island Svalbard
0
20
100
200
Great B
78
300
200
300
76
74
va
10
No
0
Norway
ya
Ze
ml
ya
0
20
72
Ru
ssia
70
Summer 1996
Figure 2. 1996 winter survey in the Barents Sea. Trawl stations are marked with . Depth isolines are in metres.
aggregated sA-values over the water column, as in
Stensholt and Nakken (2001). To normalize the large
spatial variation in absolute values, for each sample s the
relative cumulative acoustic abundance C(s,v) of fish is
calculated in 10 m depth intervals from the surface down
to 10 · v metres, v=1, 2, . . ., D, with bottom depth in the
interval (10 · (D-1), 10 · D]. For each depth step (v=1, 2,
. . .) the calculations of C(s,v) need samples with bottom
depth at least 10 v metres, and thus the number of
samples is reduced with increasing v.
To correct for unequal bottom depth, the water
column is partitioned according to fractions of the
bottom depth, and the relative sA-values are accumulated from the surface down to the chosen fractions of
the bottom depth.
For each level of depth (Figure 16f) or relative depth
(Figure 17) the distribution over all samples of the
relative cumulative values are presented as box-plots.
The boxes indicate the inter-quartile ranges and a line
inside each box marks the median.
Relative pressure reduction
For physoclistous fish the extent of rapid vertical movement is restricted because of their closed swimbladder.
At any depth level, the change of swimbladder volume
and change of buoyancy status depend on the relative
change of pressure. For cod adapted to neutral buoyancy at depth d, the free vertical range (FVR) was
determined experimentally as the range between 25%
pressure reduction (at depth m) and 50% pressure
increase (at depth M) when compared to the pressure at
depth d (Harden Jones and Scholes, 1985). Thus the
extent of FVR increases with depth. With depths in
metres, the pressure ratios are
(m+10) · (d+10) 1 =0.75, (M+10) (d+10) 1 =1.5, and
(M+10) · (m+10) 1 =2
Thus m+10=0.5 · (M+10), i.e. the pressure is reduced
by 50% for fish ascending from M to m. A full free
vertical range for cod with the deeper limit at M (>10)
metres stretches up the water column to m=M/2–5 m. In
particular, when M is the bottom depth, the ‘‘bottom
restricted free vertical range’’ (FVRbot) stretches from
the seabed up to 5 m above the mid-water level. For
other physoclistous species the free vertical range may
have a different size, i.e. from m to M where
M+10=K · (m+10), and where K may be different from
2, but for the species considered here the results below
indicate that it is not much different.
The free vertical range of physoclistous fish adapted
to pelagic living, e.g. schooling fish such as blue whiting
in the deep sea, may be calculated using m equal to the
684
B. K. Stensholt et al.
10
15
20
25
en
erg
30
sb
pit
Lon
gitu
35
S
de E
40
45
50
55
80
100
200
00
Great B
10
0
50 0
0
20
78
200
300
76
00
74
va
10
No
0
Norway
ya
Ze
ml
ya
2
72
Ru
ssia
70
Summer 1996
Figure 3. 1996 summer survey in the Barents Sea. Trawl stations are marked with . Depth isolines are in metres. N shows the
release site of tagged cod.
minimum depth of the distribution layer or M equal to
the maximum depth of the distribution layer.
The acoustic density coefficient (sA-value) depends
on the fish buoyancy status. Moreover the acoustic
sampling unit represents a rather short time interval,
typically about 6 min per nautical mile along the
transect line. Therefore, a description of the vertical
distribution of acoustic density should be related to the
physiological limitation permitting rapid pressure
change and, consequently, rapid vertical migration for
individual physoclistous fish.
In the case of demersal physoclists, it is most useful to
investigate the relative vertical profiles of sA-values in
terms of relative pressure reduction with reference to
bottom pressure. Here the bottom-restricted, free vertical range (FVRbot) is defined as the FVR with its deepest
end (M) at the seabed, and its extension depends on the
fish physiology. It is used as a yardstick for measuring
the height above the seabed for investigation of the
acoustic vertical profiles. Since m+10=(M+10)K 1, a
fish that moves up one FVR experiences a relative
pressure reduction of 1K 1 (K=2 for cod).
The height above seabed is partitioned into steps
of the relative pressure reduction level (RPRL) with
reference to the bottom pressure. The samples with
various depth ranges due to unequal bottom depths are
transformed to have ranges from RPRL=0, at the
seabed with depth M, to M · (M+10) 1, at the sea
surface. In this study the relative pressure reduction at
the sea surface can be considered to be approximately 1
in all cases. Consequently, all samples with various
depth ranges have the bottom restricted FVR stretching
from RPRL=0 to RPRL=1K 1 (RPRL=0.5 in the
case of cod). This normalization allows the comparison
of relative profiles from locations with different bottom
depths. Then the factors affecting the vertical distribution can be investigated. The graphs show the extent of
the distribution range and the vertical distribution profile within each distribution range in terms of an FVR
and in relation to the bottom.
In the case of a pelagic distribution layer of a
physoclist species over great bottom depth, the reference
is to the deepest boundary of the distribution range
instead of bottom depth. The graphs show the extent of
the distribution range and the vertical distribution profile within each distribution range in terms of an FVR,
which is independent of the depth level.
Relative profiles
To normalize for the large variation in acoustic
abundance, the profile of each sample is converted into a
Vertical density distributions of fish
685
100
y75
80
y50
Box marks interquartile range,
25–75, of the distribution.
Line inside the box marks median
60
y25
40
20
cod
0
0.0
Seabed
x
0.2
0.4
0.6
0.8
1.0
RPRL
Figure 4. The horizontal axis shows the water column partitioned from seabed to surface according to the relative pressure
reduction level (RPRL) with reference to the pressure at the seabed. For each acoustic sample consisting of area backscattering
coefficients (sA-values), the relative acoustic density is accumulated in a stepwise manner from the seabed up to each RPRL. At
RPRL=x, the distribution over the entire sample set of the relative cumulative values up to RPRL=x is presented as a thin box,
showing the inter-quartile range and median. The relative cumulative distribution up to RPRL=x has median=y50, lower
quartile=y25, and upper quartile=y75, as shown on the vertical axis. Reference line indicates RPRL=0.5.
relative profile with the observed total sA-value in the
water column as 100%. Accumulating the relative
sA-value stepwise from the bottom up along the RPRL
axis, it grows from 0 to 100% as RPRL increases from 0
to 1.
In most of the cases studied here the fish are demersal,
and the increments of the relative cumulative vertical
profiles start from the seabed. The range of RPRL,
where the relative cumulative value grows from 0% until
it reaches 100% for most of the samples, may indicate
the presence of a physiological restriction in rapid
vertical movement due to the relative change of pressure
(Figure 4). This range of physiological restriction may
correspond to the free vertical range for the species, if
no other environmental factors influence the vertical
distribution.
In the case of pelagic distribution layers, the fish
vertical distribution is off the bottom entirely, as for
redfish and blue whiting shoals in the deep sea. Thus the
relative cumulative vertical profiles should be started
from the deepest end of each sample distribution range
(Figures 9e,f and 16a–e). When the profiles are described
in terms of the relative pressure reduction with respect to
the pressure at seabed the vertical extent of the school
stretches from RPRL=XM upward to RPRL=Xm, as
in the samples from 1995, 1996, 1998 of Figure 15f. In
order to compare the difference Xm XM to the FVR in
terms of relative pressure reduction, it must be converted
and expressed in terms of relative pressure reduction
from the deeper to the upper end of the distribution
range as a ratio
(Xm XM) · (1XM) 1 =
[(1XM)(1Xm)] · (1XM) 1 =
1[(1Xm) · (1XM) 1].
The last form corresponds to and may be compared to
1K 1 above.
Adaptation to pelagic life
For a physoclist, a rapid ascent too far above the upper
limit of its free vertical range is much more hazardous
686
B. K. Stensholt et al.
than a rapid descent below the lower limit (Tytler and
Blaxter, 1973). The latter results in temporary, excessive
under-buoyancy (Harden Jones, 1951, 1952; Harden
Jones and Scholes, 1981, 1985). When the observed
relative cumulative profile reaches 100% beyond one
bottom restricted FVR unit, there are samples with fish
detected more than one such unit above the sea bottom.
Based on the assumptions that the FVR restriction to
rapid vertical movement is effective, and that such a fish
has not ascended above its current FVR, we conclude
that the fish have adapted to pelagic living in the sense
that their current free vertical range does not include the
seabed.
Fish that are detected within the bottom-restricted
FVR may well have an individual, current FVR that
does not include the seabed. This means that they just
happen to be in the part of their current FVR that
partially overlaps with the bottom-restricted FVR, or
that they have accepted excessive negative buoyancy. It
is impossible to recognize from the acoustic signals if
there are fish adapted to pelagic living among those
detected within the bottom-restricted FVR. Therefore it
is necessary to use the very strict criterion that some fish
must have been detected above one bottom-restricted
FVR before concluding that part of the stock has
adapted to pelagic living.
If the interval [m, M] is the bottom-restricted FVR,
and there are fish with individual current free vertical
range [c, C], where c<m<C<M, then some of them may
happen to be detected in the depth interval [m, C], and
others, with excessive negative buoyancy, in [C, M].
Only fish that are detected in [c, m] satisfy the criterion.
If a fish detected in [m, M] actually has adapted to
pelagic living, this cannot be verified from the observed
samples.
Distribution of relative cumulative sA-values
In Figure 4 the water column is partitioned into RPRL
steps of 0.05 from ‘‘bottom’’ to surface, shown here
along the horizontal axis. The word ‘‘bottom’’ here
means the seabed in case of demersal distribution and
the deepest end of each sample distribution range for a
clear pelagic distribution. At each discrete RPRL, the
box-plot shows the distribution of the observed relative
cumulative sA-values. In each acoustic sample the relative sA-value is accumulated from bottom up to
RPRL=x. The distribution of these relative cumulative
values up to each discrete RPRL over the entire sample
set is presented as a box showing the inter-quartile range
with a line inside each box marking the median. The
accumulation levels, in percent, appear on the vertical
axis. The distribution of cumulative relative sA-values
up to RPRL=x has median=y50, lower quartile=y25,
and upper quartile=y75. For example, the lower
quartile=y25 means that 25% of the samples have the
values below, and 75% of the samples have the values
above y25.
The line joining the medians for all RPRL is used to
represent the relative cumulative acoustic density from
bottom up to each RPRL. The inter-quartile range
represents the variation of the distribution. If an RPRL
has large inter-quartile range, there is large variability
among the samples of the relative cumulative acoustic
density up to that level.
The FVR unit for each species is then estimated from
the median, i.e. an RPRL range where the median starts
to grow from 0% until it approaches 100%, and with the
quartiles indicating the uncertainty.
As explained above, a sample with fish above the
bottom-restricted FVR indicates that some of the fish
have adjusted to pelagic living. When the lower quartile
is clearly below 100% at the RPRL corresponding to the
upper limit of the bottom-restricted FVR (RPRL=0.5 in
the case of cod), at least 25% of the samples record fish
adapted to pelagic living. Of course this interpretation of
the figures builds on the assumption that the defined free
vertical range is an effective restriction for rapid vertical
movement.
The day-night time is calculated according to the sun’s
angle with the horizon, for winter surveys daytime starts
and ends when the sun is 5 below horizon. Each
acoustic sample is accordingly classified as a day- or
night-sample, and the vertical profiles may be presented
for day and night separately and compared.
The acoustic bottom dead zone and the day/night
cumulative relative frequencies
The observed total acoustic sA-value in the water
column may differ significantly between day and night.
Such differences could be caused by diurnal changes due
to fish behaviour and buoyancy status, combined with
technical limitations of the hydro-acoustic method to
detect fish (Harden Jones and Scholes, 1981, 1985;
Aglen, 1994). Thus the basis for calculating relative
frequencies changes correspondingly. Based on studies
of demersal species during winter in the Barents Sea
(Aglen et al., 1999), it is reasonable to assume that the
main reason for the difference in day and night loss
is migration of fish in and out of the near bottom
acoustic dead zone. An adjustment for the difference in
the day and night loss is based on this assumption and
that the total amount of fish is equal during day and
night.
Only when the daytime and night-time vertical profiles
are based on the same proportion of the fish present in
the water column, will the comparison of day/night
profiles give evidence of fish vertical movements. Thus
the observed relative frequencies have been adjusted
using the relative change in observed total acoustic
sA-value between day and night.
Vertical density distributions of fish
Starting with one acoustic sample, let A(h) be the
observed cumulative acoustic density from seabed up to
RPRL=h. With bottom depth b metre,
h=1[(d+10) · (b+10) 1]=(bd) · (b+10) 1
at depth d metres. At the surface, d=0, so
h=b · (b+10) 1, and N=A[b · (b+10) 1] is the
observed total acoustic sA-value.
The observed relative cumulative acoustic density is
F(h)=A(h) · N 1. Let X be the loss due to the acoustic
dead-zone. If the equipment had been able to detect fish
at the bottom, A(h) would have been replaced by
X+A(h), the total acoustic sA-value would have been
X+N, and the relative cumulative density at level h
would have been G(h)=[X+A(h)] · (X+N) 1. The relation between the observed F and the ‘‘correct’’ G is
therefore expressed by
G(h)=X · (X+N) 1 +[N · (X+N) 1] · F(h)
=X · (X+N) 1 +{1[X · (X+N) 1]} · F(h)
(1)
The intercept G(0)=X · (X+N) 1 is the relative loss
in the dead-zone, and N · (X+N) 1 is the reduction
factor for the slope. Increased loss X means higher
intercept and lower slope for G. In particular,
F(0)=0%, G(0)=X · (X+N) 1,
F(b · (b+10) 1)=G(b · (b+10) 1)=100%,
G(h)=F(h) · N · (X+N) 1,
1G(h)=[1F(h)] · N · (X+N) 1.
For each sample, s, one could draw a similar curve.
The observed total acoustic sA-value is N(s), and the
unknown relative loss is X(s) · [X(s)+N(s)] 1, which
could have been used to correct the curve. All sample
curves would have been shifted upward according to
Equation (1). This method can be applied to correct for
the unequal losses at the bottom due to different
reasons, e.g. day-night, or tidal current, or before and
after the passage of a vessel.
We use the median curve as a ‘‘virtual sample’’ to
represent the data. Let virtual dead-zone losses be Xd
and Xn, and observed total acoustic sA-value be Nd and
Nn for day and night, respectively. The total amount of
fish at any location does not change between day and
night, hence Xd +Nd =Xn +Nn. For Nd and Nn respectively, we may use the median of the observed day and
night total acoustic sA-value. For the purpose of day and
night comparison, we only need the difference
Xn Xd =Nd Nn.
Typically the daytime and night-time median curves
are concave, reach 100%, cross each other at RPRL=hc
(Figure 5), and in an interval around hc the daytime
curve is steeper than the night-time curve. There is then
a higher relative density of fish at day than at night in
this interval.
687
Before correction, the day and night cumulative frequencies are based on different total observed acoustic
densities. In order to draw conclusions about vertical
migration, one must consider the correction for difference in dead-zone loss, so that both curves are based on
the same total acoustic sA-value. Then one median curve
is shifted accordingly.
If Xn Xd =0, no correction is needed, and hc remains
the same.
If Xn Xd =Nd Nn >0, the correction shifts the
night-time
curve
upwards
to
intercept
at
(Nd Nn) · (Nd) 1. This intercept indicates the daynight difference in the relative acoustic density that is
lost due to fish in the dead-zone.
With X=Nd Nn and N=Nn, Equation (1) can be
rewritten in terms of the observed values as
G(h)=(Nd Nn) · (Nd) 1 +
{1[(Nd Nn) · (Nd) 1]} · F(h)
=F(h)+[(Nd Nn) · (Nd) 1] · [1F(h)].
Now, in the typical case, daytime accumulation
reaches 100% at a level RPRL=h* where the night value
is <100%. In other words, the median curves of day/
night relative cumulative acoustic sA-values indicate that
RPRL=h* is the upper limit of the daytime vertical
distribution, but that among night samples a significant
proportion of fish were still found above RPRL=h*. It
is then geometrically clear that the corrected night curve
still crosses the day curve, say at RPRL=hk, where
hc <hk <h* (Figure 5). Hence the median curves show
that in fact there must be net evening migrations both
upwards and downwards away from the level
RPRL=hk.
If Xn Xd =Nd Nn <0, there must be a similar
upward shift of the daytime curve (Figure 7). For small
shifts the intersections at 0 and hc are replaced by two
intersections between 0 and hc, and for large shifts there
is no intersection.
The normalizations for unequal bottom depths and
unequal acoustic densities together with corrections of
unequal bottom dead zone loss enable us to compare
vertical profiles and discuss vertical migration from the
overall line transect survey data. In particular the influence of some environmental factors on the vertical
migration can be discussed, e.g. daylight, bottom depth,
current and fish density.
Results
Vertical distribution of fish over different seasons
and geographical areas
The distribution of relative cumulative acoustic density
(sA-values) is calculated for all years combined as well as
annually for each physoclist and physostome of each
688
B. K. Stensholt et al.
h*
hk
100
80
night
60
40
20
day
haddock
50–600 m
0
0.0
hc
0.2
0.4
0.6
0.8
1.0
Figure 5. The median curves of day/night relative cumulative acoustic profiles indicate RPRL=h* as the upper limit of the daytime
vertical distribution, but among night time samples a significant proportion of fish are still found above RPRL=h*. The adjusted
night curve (broken line) crosses the day curve at RPRL=hk where hc <hk <h*. The higher slope of the daytime curve around
RPRL=hk shows higher density at daytime than night-time. Hence the median curves show that there must be net evening
migration both upwards and downwards away from the level RPRL=hk.
survey. The results are shown in Figures 6–17. For all
physoclists the relative cumulative vertical profiles are
started from the seabed, except for redfish when it is
pelagic (Figure 9e,f) and for blue whiting in the deep sea
(Figure 16a–e). In these cases the profiles are started
from the deepest boundary of each sample distribution
range. The figures are presented for each survey and
species, and within each species the acoustic samples are
grouped according to day-night or to density level and
bottom depth range, as indicated at the bottom right
corner of the figures. The species are cod and haddock
from winter surveys in the Barents Sea (Figures 6–8);
redfish from winter surveys in the Barents Sea and
Irminger Sea (Figure 9); cod, haddock, and redfish from
summer surveys in the Barents Sea (Figures 10–12);
saithe, haddock, and blue whiting from surveys along
the Norwegian coast (Figures 13, 14); saithe and blue
whiting from surveys in the North Sea (Figure 15); blue
whiting from surveys west of the British Isles (Figure 16)
and capelin and herring from winter surveys in the
Barents Sea (Figure 17).
The results for demersal fish surveys reflect the vertical
distribution of fish found in the shelf sea (Figures 1–3).
There are data sets where a few samples have very high
density and account for a high proportion of the
sampled fish. They have special profiles and this is
not reflected in the median of the entire collection
of samples. Therefore, such samples are investigated
separately.
The annual median curves of the relative cumulative
sA-values show the annual variation in vertical distribution (e.g. Figure 6b). The reference line at RPRL=0.5
indicates one bottom restricted FVR for cod. For other
species, the RPRL where the median gets close to 100%
and the inter-quartile range approaches 0 may indicate
the extent of their free vertical ranges. They seem to be
of roughly the same size as for cod (Figures 6a,
8a–13a,c,e, 9f, 15a,c, 16d,e,f). Below the reference line
Vertical density distributions of fish
100
100
(a)
80
80
60
60
40
40
20
20
689
(b)
sA>10
cod
cod
0
0
0.0
100
0.2
0.4
0.6
0.8
1.0
0.0
100
(c)
80
0.2
0.4
0.6
0.8
1.0
(d)
day
80
60
60
night
40
40
20
20
night
sA>1
cod
cod
50–600 m
0
0
0.0
100
0.2
0.4
0.6
0.8
1.0
0.0
100
(e)
80
80
60
60
40
40
sA>1
20
0.2
0.4
0.6
0.8
(f)
day
night
sA>10
20
cod
cod
50–200 m
0
0.0
0.2
0.4
0.6
0.8
1.0
1.0
200–600 m
0
0.0
0.2
0.4
0.6
0.8
1.0
Figure 6. Distribution of cod relative cumulative vertical profiles from the seabed to the surface in terms of RPRL, from winter
surveys in the Barents Sea, 1996–2001. Median and inter-quartile range (represented by thin boxes, left box for day, right box for
night) of the distribution, (a) all data; (b) annual median curves; (c) day- and night samples; (d) day–night median curves with
upward adjusted night median curve (broken line); (e) samples with bottom depth <200 m; (f) samples with bottom depth
200–600 m, more cod adapt to pelagic living.
690
B. K. Stensholt et al.
100
80
night
60
40
day
cod 1996
50–600 m
20
0
0.0
0.2
0.4
0.6
0.8
1.0
Figure 7. Day and night median curves of the distributions of
cod vertical profiles of relative cumulative acoustic sA-values in
relation to relative pressure reduction level, RPRL, with
upward adjusted day median curve (broken line). Samples are
from the winter survey in the Barents Sea, 1996. Reference line
indicates RPRL=0.5.
and close to the bottom the inter-quartile range is
relatively large, indicating a large variation of the relative cumulative values among different acoustic samples.
This means that within one FVR the fish are free to
distribute along the water column in a way suited to
local conditions.
Pelagic living for some fish in the samples is indicated
when the relative cumulative sA-values show there are
fish higher than RPRL=0.5, i.e. the lower quartile
reaches 100% beyond the reference line (Figures 6c,f, 8c,
9c, 10e, 11e, 13c, 15c,e). In such cases more than 25% of
the samples have some fish adapted to pelagic living.
This is observed especially in samples from areas of high
acoustic abundance and where the seabed is deeper than
200 m (Figures 6f, 8c, 9c, 10e, 11e, 14b,d). The most
notable species may be oceanic redfish and blue whiting
(Figures 15e,f, 16e,f), which are pelagic at great depth.
For surveys in the Barents Sea, the distribution at
each RPRL has larger inter-quartile range in summer
than winter (Figures 6a, 8a–12a). The distribution in
summer 1999 shows more pelagic living for cod and
haddock than in most years. The depth distribution
100
100
(a)
(b)
80
80
60
60
40
40
20
20
haddock
haddock
0
0
0.0
0.2
0.4
0.6
0.8
1.0
0.0
100
100
(c)
80
80
60
60
40
40
0.2
0.4
0.6
day
1.0
(d)
night
20
20
0.8
haddock
50–600 m
haddock
0
0
0.0
0.2
0.4
0.6
0.8
1.0
0.0
0.2
0.4
0.6
0.8
1.0
Figure 8. Distribution of haddock relative cumulative vertical profiles from the seabed to the surface in terms of RPRL, from
winter surveys in the Barents Sea, 1996–2001. Median and inter-quartile range (represented by thin boxes, left box for day, right
box for night) of the distribution (a) all data; (b) annual median curves; (c) day- and night samples; (d) day–night median curves
with upward adjusted night median curve (broken line).
Vertical density distributions of fish
100
691
100
(a)
(b)
80
80
60
60
40
40
20
20
sA>1
redfish
0
0.0
0.2
0.4
0.6
redfish
0.8
100
1.0
(c)
0
0.0
0.2
0.4
0.6
0.8
1.0
100
(d)
day
80
80
60
60
night
40
40
20
20
sA>1
redfish
0
0.0
0.2
0.4
0.6
redfish
50–1000 m
0.8
1.0
0
0.0
0.2
0.4
0.6
1.0
100
100
(e)
80
60
60
40
40
sA>1
20
(f)
day
80
night
20
sA>1
oceanic redfish
oceanic redfish
450–3000 m
450–3000 m
0
0.0
0.8
0.2
0.4
0.6
0.8
1.0
0
0.0
0.2
0.4
0.6
0.8
1.0
Figure 9. Distribution of redfish relative cumulative vertical profiles from the seabed to the surface in terms of RPRL, from winter
surveys in the Barents Sea, 1996–2001. Median and inter-quartile range (represented by thin boxes, left box for day, right box for
night) of the distribution (a) all data; (b) annual median curves; (c) day- and night samples; (d) day–night median curves with
upward adjusted night median curve (broken line); relative cumulative profiles from the deepest of the distribution layer to surface
of samples from Irminger Sea, 2001 (e) entire redfish distribution range; (f) only above and within deep scattering layer.
692
B. K. Stensholt et al.
100
100
(a)
(b)
80
80
60
60
40
40
20
20
sA>1
cod
cod
0
0.0
0.2
0.4
0.6
0.8
1.0
100
0
0.0
0.2
0.4
0.6
0.8
100
(c)
(d)
80
80
60
60
40
40
20
20
cod
sA>1
cod
50–200 m
0
0.0
0.2
0.4
0.6
50–200 m
0.8
1.0
100
0
0.0
0.2
0.4
0.6
0.8
1.0
100
(e)
(f)
80
80
60
60
40
40
20
20
cod
sA>10
cod
200–600 m
200–600 m
0
0.0
1.0
0.2
0.4
0.6
0.8
1.0
0
0.0
0.2
0.4
0.6
0.8
1.0
Figure 10. Distribution of cod relative cumulative vertical profiles from the seabed to the surface in terms of RPRL, from summer
surveys in the Barents Sea, 1995–2000. Median and inter-quartile range (represented by thin boxes) of the distribution (a) all data;
(b) annual median curves; (c) and (d) samples with bottom depth<200 m; (e) and (f) samples with bottom depth 200–600 m,
showing more cod adapt to pelagic living.
of individual cod from data storage tags has larger
variation in summer than in winter (Tables 1, 2), in
agreement with the above results.
Proportion of fish in the bottom channel
For the physoclist species, the ratio of the observed
sA-value from the bottom channel (c10m above seabed)
Vertical density distributions of fish
100
693
100
(a)
(b)
80
80
60
60
40
40
sA>1
haddock
20
0
0.0
0.2
0.4
0.6
haddock
20
0.8
1.0
100
0
0.0
0.2
0.4
0.6
0.8
100
(c)
(d)
80
80
60
60
40
40
haddock
20
0
0.0
1.0
0.4
0.6
haddock
20
50–200 m
0.2
sA>1
0.8
1.0
100
0
0.0
100
(e)
80
50–200 m
0.2
0.4
40
40
20
20
1.0
(f)
99
97
60
0.8
2000
80
60
0.6
96
98
sA>10
haddock
haddock
200–600 m
200–600 m
0
0.0
0.2
0.4
0.6
0.8
1.0
0
0.0
0.2
0.4
0.6
0.8
1.0
Figure 11. Distribution of haddock relative cumulative vertical profiles from the seabed to the surface in terms of RPRL, from
summer surveys in the Barents Sea, 1995–2000. Median and inter-quartile range (represented by thin boxes) of the distribution (a)
all data; (b) annual median curves; (c) and (d) samples with bottom depth<200 m; (e) and (f) samples with bottom depth
200–600 m, showing more haddock adapt to pelagic living.
to the observed total sA-value was calculated in each
sample and the distribution of the ratio investigated.
The median of the distribution is shown in Tables 3–4. A
general feature is that in deeper waters fish distributions spread relatively higher up into the water column.
Table 3 shows that cod and haddock have a higher
proportion in the bottom channel during night than
day but this is based on the observed acoustic values
without adjustment for unequal losses during day and
night.
694
B. K. Stensholt et al.
100
100
(a)
(b)
80
80
60
60
40
40
sA>1
20
sA>1
20
redfish
redfish
0
0.0
0.2
0.4
0.6
0.8
1.0
0
0.0
0.2
0.4
0.6
0.8
1.0
100
100
(c)
(d)
80
80
60
60
40
40
20
sA>1
sA>1
redfish
50–200 m
redfish
50–200 m
20
0
0.0
0.2
0.4
0.6
0.8
1.0
100
0
0.0
0.2
0.4
0.6
0.8
1.0
100
(f)
(e)
80
80
60
60
40
40
sA>10
sA>10
redfish
20
redfish
20
200–1000 m
200–1000 m
0
0.0
0.2
0.4
0.6
0.8
1.0
0
0.0
0.2
0.4
0.6
0.8
1.0
Figure 12. Distribution of redfish relative cumulative vertical profiles from the seabed to the surface in terms of RPRL, from
summer surveys in the Barents Sea, 1995–2000. Median and inter-quartile range (represented by thin boxes) of the distribution (a)
all data; (b) annual median curves; (c) and (d) samples with bottom depth <200 m; (e) and (f) samples with bottom depth
200–1000 m.
Table 2, based on depth records from data storage
tags, shows that individual cod spends more time in
the lower part of the daily depth range (within 10 m
above the daily maximum depth) during December–
January than during other periods. When cod is in the
‘‘lower part’’ it is likely that for some of the time it
would have been detected by acoustics and would have
contributed to the sA-value in the bottom channel,
and at other times it is in the acoustic bottom dead
zone.
Vertical density distributions of fish
100
695
100
(a)
(b)
80
80
60
60
40
40
sA>1
20
0
0.0
20
saithe
30–600 m
0.2
0.4
0.6
0.8
1.0
100
saithe
30–600 m
0
0.0
0.2
0.4
0.6
0.8
1.0
100
(c)
(d)
80
80
60
60
40
40
day
sA>1
20
20
haddock
30–600 m
0
0.0
0.2
0.4
0.6
0.8
1.0
100
haddock
30–600 m
0
0.0
0.2
0.4
0.6
0.8
1.0
100
(e)
(f)
day
80
80
night
60
night
60
40
40
sA>1
20
20
blue whiting
blue whiting
30–600 m
30–600 m
0
0.0
0.2
0.4
0.6
0.8
1.0
0
0.0
0.2
0.4
0.6
0.8
1.0
Figure 13. The distribution of relative cumulative vertical profiles from the seabed to the surface in terms of RPRL, from surveys
along the Norwegian coast 1992–2000. Day and night medians and inter-quartile ranges (represented by thin boxes, left box for
day, right box for night), from combined data for (a) saithe; (c) haddock; (e) blue whiting. Day–night median curves with upward
adjusted night median curve (broken line) for (b) saithe; (d) haddock; (f) blue whiting.
696
B. K. Stensholt et al.
100
100
(a)
(b)
80
80
60
60
40
40
sA>1
20
0
0.0
saithe
30–200 m
0.2
0.4
0.6
sA>5
20
0.8
1.0
100
0
0.0
saithe
30–600 m
0.2
0.4
0.6
0.8
100
(c)
(d)
80
80
60
60
40
40
sA>1
20
0
0.2
0.4
0.6
sA>1
20
haddock
30–200 m
0.0
0.8
1.0
100
0
0.0
haddock
30–600 m
0.2
0.4
0.6
0.8
1.0
100
(e)
(f)
80
80
60
60
40
40
sA>1
20
sA>10
20
blue whiting
blue whiting
300–600 m
200–300 m
0
0.0
1.0
0.2
0.4
0.6
0.8
1.0
0
0.0
0.2
0.4
0.6
0.8
1.0
Figure 14. The distribution of relative cumulative vertical profiles from the seabed to the surface in terms of RPRL, from surveys
along the Norwegian coast 1992–2000. Annual median curves for samples with bottom depth <200 m of (a) saithe; (c) haddock;
(e) blue whiting; and for samples with bottom depth 300–600 m of (b) saithe; (d) haddock; (f) blue whiting. The level of sA-value
indicated at bottom right corner.
Vertical density distributions of fish
697
100
100
(a)
(b)
80
80
60
60
40
40
sA>1
20
0
0.0
sA>1
20
saithe
30–600 m
saithe
30–600 m
0
0.2
0.4
0.6
0.8
1.0
100
0.0
0.2
0.4
0.6
0.8
1.0
100
(c)
(d)
80
80
60
60
1999
1998
40
40
sA>1
20
0
0.0
20
blue whiting
30–600 m
0.2
0.4
0.6
0.8
1.0
100
0
0.0
blue whiting
30–600 m
0.2
0.4
0.6
0.8
100
(e)
(f)
80
80
60
60
40
40
94
98
99
96
20
blue whiting
95
97
sA>1
20
300–600 m
0
0.0
1.0
0.2
0.4
0.6
0.8
1.0
0
0.0
blue whiting
300–600 m
0.2
0.4
0.6
0.8
1.0
Figure 15. The distribution of relative cumulative vertical profiles from the seabed to the surface in terms of RPRL, from surveys
in the North Sea 1992–2000. Median and inter-quartile range (represented by thin boxes) of (a) saithe; (c) and (e) blue whiting.
Annual median curves of (b) saithe; (d) and (f) blue whiting. Acoustic sA-values level, species and bottom depth ranges are shown
at the bottom right corner.
698
B. K. Stensholt et al.
100
100
(a)
(b)
80
80
60
60
night
40
40
1995
20
1995
sA>10
20
sA>10
300–600 m
0
0.0
0.2
0.4
0.6
0.8
600–1500 m
1.0
0
0.0
0.2
0.4
0.6
0.8
1.0
100
100
(c)
(d)
80
80
60
60
40
40
night
1996
20
1996
sA>10
night
20
sA>10
200–300 m
0
0.0
0.2
0.4
0.6
0.8
300–600 m
1.0
0
0.0
0.2
0.4
0.6
0.8
1.0
100
100
(e)
(f)
80
80
60
60
night
40
40
1996
20
sA>10
night
1996
sA>10
20
600–1500 m
600–1500 m
0
0
0.0
0.2
0.4
0.6
0.8
1.0
0
200
400
600
800
Figure 16. Median (full or broken line join day or night median) and inter-quartile range (represented by thin boxes, left box for
day, right box for night) of blue whiting from the distribution of relative cumulative vertical profiles from the deepest end of each
sample distribution layer to surface in terms of RPRL. Samples are from surveys west of British Isles. Bottom depth and year are
as specified at bottom right corner. The pelagic distribution layers are shown in (f), with accumulation from surface down in
metres.
Vertical density distributions of fish
100
100
(a)
80
80
60
60
40
40
20
699
(d)
20
capelin
0
0.0
herring
0
0.2
0.4
0.6
0.8
1.0
100
0.0
100
(b)
80
80
60
60
40
40
0.2
0.4
0.6
(e)
capelin
herring
0
0
100
100
(c)
(f)
80
80
60
60
day
night
40
20
night
20
capelin
0
0.0
1.0
20
20
40
0.8
0.2
0.4
0.6
0.8
herring
1.0
0
0.0
0.2
0.4
0.6
0.8
1.0
Figure 17. The distribution of relative cumulative vertical profiles from the surface to the bottom in terms of relative depth (0.0 at
surface, 1.0 at bottom), from winter surveys in the Barents Sea, 1996–2001. Median and inter-quartile range (represented by thin
boxes, left box for day, right box for night) of the distribution for (a) capelin; (d) herring; annual median curves (b) capelin;
(e) herring; day and night for (c) capelin; (f) herring.
Dead-zone loss during day and night
The typical case in the Barents Sea winter surveys
1996–2002 was that the observed acoustic abundance
during the daytime was higher than at night, Nd >Nn, i.e.
that night-time loss in the acoustic dead-zone exceeds
daytime loss. The night-time median curve is adjusted by
an upward shift. The values for Nd and Nn are the
medians of the total acoustic sA-value in the water
column over the day and night samples, respectively.
Only the samples with total sA-value greater than 1 are
used here.
700
B. K. Stensholt et al.
Table 1. Temperature and depth distribution, in winter (January–March) and summer (July–August). From 1996–1997 time series
of 6 data-storage tags attached to North East Arctic cod released in the Barents Sea.
Tag no.
Temperature
117
131
191
204
206
246
Depth
117
131
191
204
206
246
Summer
Winter
Summer
Winter
Summer
Winter
Summer
Winter
Summer
Winter
Summer
Winter
Summer
Winter
Summer
Winter
Summer
Winter
Summer
Winter
Summer
Winter
Summer
Winter
Mean
Std
Min
Max
3.1
3.9
2.2
2.4
3.9
4.0
2.6
4.0
3.1
2.4
2.4
4.2
2.0
0.8
1.5
1.4
1.2
0.3
1.5
0.4
1.2
0.9
1.5
0.3
1
1.5
1.3
0.1
1.3
3.4
0.7
3.0
0.3
0.1
1.0
3.5
6.4
5.2
5.0
5.0
9.2
4.7
4.9
4.9
6.3
3.6
5.3
5.3
237
286
208
261
209
318
236
309
130
201
235
311
78
75
69
43
96
40
84
31
92
94
81
11
0
0
60
81
1
174
2.4
193
0
49
3
241
For the six winter surveys and the three species
cod, haddock, and redfish, the dead zone corrections,
i.e. the relative difference in losses at the bottom
(Nd Nn) · (Nd) 1 if Nd >Nn (night correction), or
(Nn Nd) · (Nn) 1 if Nd <Nn (day correction), are given
in Table 5. The material may be insufficient to allow the
conclusion that 1996, with its large day corrections, is an
exceptional year (Figure 7). One possibility is that there
is a separation of both the cod and the haddock stocks
into groups with very different diurnal vertical migration
patterns and that switches between Nd <Nn and Nd >Nn
from one year to another is explained by relative
changes in the abundance of such groups. The relative
475
475
406
371
410
404
389
376
353
353
370
331
P5
0.6
2.5
0.1
0.2
2.1
3.6
0.4
3.4
1.0
0.6
0.2
3.5
113
176
92
198
52
235
103
256
20
72
92
297
Q1
1.3
3.4
0.9
1.0
3.4
3.8
1.3
3.8
2.1
1.7
0.9
3.9
174
271
151
237
136
303
150
288
55
86
166
303
Median
3.9
4.2
2.1
2.9
3.9
4.0
3
4.0
3.3
2.6
2.7
4.1
243
296
218
255
224
321
272
308
99
207
264
312
Q3
4.7
4.5
3.6
3.5
4.3
4.2
3.8
4.4
3.9
3.1
3.7
4.4
294
322
257
293
281
344
305
334
207
293
300
320
P95
5.4
4.8
4.1
4.3
6.3
4.6
4.4
4.6
5.3
3.3
4.4
4.6
353
379
309
335
351
377
344
352
300
327
329
328
cumulative acoustic frequency figures for day and night
(Figures 6d, 8d, 9d) indicate that such separation exists,
although the acoustic data give no information on which
parts of the stock ascend or descend. The pattern of
vertical migration in 1996 corresponds to the high
proportion of the one-year-old fish (Table 6), in agreement with Aglen et al. (1999). Unless the species and age
allocation of acoustic samples can be made according to
age and species composition in smaller depth intervals,
such profiles may also occur in areas where different
species and age groups have opposite vertical migration.
Diel pattern in vertical distribution
Table 2. Percentage of time spent within 10 m above the daily
maximal depth from depth records during April–June (4–6),
July–August (7–8), September–November (9–11), December–
January (12–1), and February–March (2–3). Depth time series
from April 1996 to March 1997 are from 6 data-storage tags
attached to North East Arctic cod released in the Barents Sea.
Tag number
4–6
7–8
9–11
12–1
2–3
117
131
191
204
206
246
43
46
45
46
35
67
41
33
45
45
48
37
71
37
33
57
75
42
76
65
50
70
73
74
29
51
57
31
44
56
In most cases, when the samples from winter surveys
(cod, haddock, and redfish 1996–2001) are classified as
day or night, the relative cumulative distribution of
sA-values near the sea bottom shows lower median
values at daytime than at night-time. With accumulation
further up the water column, the day median increases
faster than the night median, and the order is reversed
around mid-range in the lower half of the water column
(Figures 6d, 7d, 8d). This reversal appears as the day
median curve’s crossing and passing over the night
median curve, with a steeper slope of the day median
curve than the night median curve. When the crossing
remains after the correction for unequal dead zone loss,
it shows that during night-time a higher proportion of
Vertical density distributions of fish
701
Table 3. Cod and haddock: medians of the distributions of the relative acoustic sA-values from bottom channel (10 m above
seabed) in percent of total sA-values. Acoustic samples are from summer and winter surveys in the Barents Sea 1995–2002 and are
classified according to bottom depth interval, day and night.
Summer survey
Winter survey
50–600 m
Cod
Cod
Year
50–200 m
200–600 m
1995
1996
1997
1998
1999
2000
2001
2002
97.3
82.7
93.3
99.9
74.7
100.0
26.3
23.8
9.7
15.7
11.7
52.3
Haddock
Haddock
50–200 m
200–600 m
50–200 m
200–600 m
Day
Night
Day
Night
55.6
67.0
55.5
63.6
59.6
52.9
72.9
31.5
36.0
24.5
24.7
19.1
39.3
37.2
72.7
58.2
69.1
57.6
49.5
42.8
39.9
25.5
22.4
23.2
18.0
13.1
37.5
16.0
29.4
38.4
20.5
22.5
13.7
22.5
26.2
37.0
39.3
30.2
31.5
29.0
52.6
46.9
21.9
18.8
18.2
10.4
13.6
25.6
12.7
37.5
28.7
30.5
29.2
19.9
49.5
21.4
fish is distributed near the bottom while another higher
proportion also distributes further up into the water
column than during daytime. During daytime there is a
higher proportion of fish around mid-range of the
bottom half of the water column. This vertical movement and replacement may cause unequal amounts of
acoustic information to be lost in the dead zone at
night-time and daytime. Moreover, the slope of the
median curve indicates the density level, and thus
doubling the slope, e.g. from 2 to 4 (i.e. a small angle at
the crossing), means a doubling of the density.
Table 5 shows that for cod and haddock there is a
higher loss during night than day in five out of seven
years, and for redfish night loss is the higher in six out of
seven years. After correction, the median profiles for
‘‘day’’ and ‘‘night’’ are based on equal dead zone losses.
Table 4. Median and lower quartile (P25) of the relative
acoustic sA-values from bottom channel (10 m above seabed) in
percent of total sA-values. Blue whiting from saithe survey
along the Norwegian coast north of 62N 1992–1997. Saithe
from North Sea survey south of 62N 1992–2000. Classification
by bottom depth.
Saithe survey
Blue whiting
50–200 m
300–600 m
Cod
North Sea survey
Saithe
30–600 m
Year
Median
Median
P25
Median
1992
1993
1994
1995
1996
1997
1998
1999
2000
100.0
68.3
68.7
100.0
15.7
100.0
2.5
13.2
35.1
14.9
11.2
9.5
41.3
100.0
100.0
100.0
46.9
100.0
63.9
100.0
97.2
81.2
100.0
100.0
100.0
86.3
100.0
100.0
100.0
100.0
The acoustic data show the net effect if there are
migrations in opposite directions but do not show to
what extent migrations with opposite direction cancel
out.
Table 5. The relative difference between day and night average
acoustic sA-values in percent of the larger of these two. The day
column shows the correction percentage when day loss>night
loss, and similarly for the night column.
Cod
Year
Day
1996
1997
1998
1999
2000
2001
2002
24.3
Haddock
Night
24.7
2.3
Day
Night
Redfish
Day
21.8
4.8
28.8
31.7
21.1
38.1
30.8
4.7
11.4
63.4
38.6
Night
20.8
31.6
54.7
7.0
29.3
39.9
9.7
Table 6. NE-Arctic cod: age composition based on
acoustic abundance indices in percent (Aglen, 2002) from the
Norwegian winter survey in the Barents Sea.
Year
1995
19961
19971*
19981*
1999
2000
2001
2002
1
Age 1
Age 2
Age 3
Age 4
Age 5
Age 6+
74.0
77.7
67.4
71.4
36.7
15.6
55.5
3.2
8.3
11.0
17.9
13.3
31.2
22.5
5.6
39.1
4.7
3.0
7.8
9.0
15.4
24.9
12.2
10.0
4.5
2.3
2.1
3.8
9.9
16.1
15.1
19.3
5.9
2.8
2.0
0.9
4.6
14.4
6.8
17.0
2.5
3.2
2.7
1.6
2.2
6.4
4.7
11.4
Survey covered a larger area.
*Adjusted indices, survey covered only Norwegian zone.
702
B. K. Stensholt et al.
The eastern part of the Barents Sea, which is an
important feeding area for young cod, was not surveyed
in 1997–1998. In 1996, 2000, and 2001 there are sufficient samples east of 3530E for analysis, and separate
investigation shows a relative day correction of 36% for
1996 and relative night corrections of 60% in 2000 and
57% in 2001. The results from the areas west of 3530E
are similar to the overall results in Table 5.
Day and night samples along the Norwegian coast
and from west of the British Isles are situated far apart
in terms of geographic locations, so only cautious conclusions can be drawn. The diel vertical migration of
saithe and haddock is not immediately clear from the
combined data 1992–2000 (Figure 13a,c). Their diel
vertical migration patterns appear more clearly when the
relative cumulative vertical profiles for night-time have
been adjusted upward by correcting for the estimated
losses in the acoustic dead-zone (Figure 13b). Yearly
data show some evidence of diel vertical migration but
the direction of migration changes from year to year.
Diel vertical migration of blue whiting is seen all years
both in samples from the Norwegian coast (Figure 13e,f)
and from the North Sea survey, when the RPRL is
calculated based on the bottom restricted FVR.
West of the British Isles the RPRL is calculated based
on the deepest boundary of each distribution range. The
1995 samples show diel behaviour in areas of depth
<600 m (Figure 16a). For depth >600 m, the day
samples are mainly from deeper areas than the
night samples; however, there is no clear diel pattern
(Figure 16b). The 1996 samples show clear diel pattern,
with higher daytime concentration at the deeper end of
the vertical distribution range (Figure 16c,d,e). Samples
from areas deeper than 600 m have the night-time distribution about 50–100 m higher up than the daytime
distribution (Figure 16f).
For oceanic redfish, the RPRL is calculated from the
deep end of the distribution range. In the relative
cumulative vertical profiles of total sA-values, i.e. above,
within, and below the deep scattering layer (Figure 9e),
the median mark is higher during the day than at night
from just above the deepest end of the distribution
range. Further up, but within RPRL=0.5, the night
median is higher, but beyond RPRL=0.5 the day
median is again higher, and the night distribution
stretches higher in the water column. This may be an
effect of how the samples are gathered, as biologists
suspect that there are two separate stocks of redfish
(Magnússon et al., 1996; Dalen and Nedreaas, 2001).
When only fish from above and within the deep scattering layer are used, then each daytime sample is mainly
distributed within one (cod) FVR. However, night
samples stretch higher up (Figure 9f).
Diel vertical migration is also found for capelin in all
winter surveys. The median indicates that capelin are
distributed throughout the entire water column during
night-time and are concentrated closer to the bottom
during the day (Figure 17c). The large inter-quartile
range may indicate that the capelin’s vertical profile
varies according to local conditions. Similar diel behaviour is found for herring in the years for which we have
a significant amount of data (1999–2000). In 1998 and
2000 most herring samples were distributed in the
bottom half of the water column (Figure 17d,e,f).
Some results species by species
Cod
Data from both winter and summer surveys
(Figures 6a,b, 10a,b) show that the vertical distribution
is in agreement with a theoretical free vertical range
(FVR) as determined by Harden Jones and Scholes
(1985). The samples with bottom depth less than 200 m
have relative vertical profiles within the bottomrestricted FVR (Figures 6e, 10c,d). With bottom depth
more than 200 m combined with high acoustic abundance, more cod adapt to pelagic living in both winter
and summer survey samples (Figures 6f, 10e,f and Table
3). Both acoustic samples and data-storage tag records
show larger variation in depth distribution during summer than in winter (Figures 6a, 10a and Tables 1, 2).
Haddock
Data from winter and summer surveys show a similar
vertical profile to that of cod. In data from all four
surveys haddock mainly stay within one ‘‘cod’’
bottom restricted FVR, but more often adapt to pelagic
living (Figures 8a,b, 11a,b, 13c,d, and 14c,d). Summer
samples have larger inter-quartile range than winter
samples (Figures 8a and 11a). This suggests a greater
variation in the relative vertical profiles in summer.
The vertical distribution in terms of pressure change
in comparison to the bottom restricted FVR seems
to depend on bottom depth and the size of acoustic
abundance (Figures 11c,d,e,f, 14c,d). Haddock samples
from areas deeper then 200 m or with high acoustic
abundance show clearly that some of the fish adapt to
pelagic living (Figures 11e,f, 14c,d, and Table 3). In the
North Sea surveys and summer surveys, however, they
spread more upward into the water column than in
winter-survey samples. The summer 1998 median curve
in Figure 11f shows that most haddock samples had
more than half of the relative acoustic density well above
the bottom restricted FVR (RPRL=0.5). The two selection criteria together were satisfied by samples in one
stretch of line transects along the northern coast of
Norway. Most individual samples were similar to the
median curve, with two separate high concentrations of
acoustic abundance, one close to the bottom and the
other in the upper half of the water column.
Vertical density distributions of fish
Redfish
In general redfish are found in regions deeper than
100 m. The vertical profile shows that redfish mainly
stay well within the 50% RPRL with respect to the
bottom pressure in both summer and winter (Figures 9
and 12). Even with depths greater than 300 m only small
parts of the samples go into the upper half of the water
column. Redfish is not a target species and none of the
surveys cover its distribution areas completely.
Saithe
Along the Norwegian coast, saithe mainly stay in the
bottom half of the water column for all bottom depths
(Figures 13a,b, 14a,b), but as depth increases offshore
they spread relatively higher into the water column
(Figure 14b). For depths less than 200 m the saithe stay
close to seabed (Figure 14a), but for depth greater than
300 m (Figure 14b) the median of relative cumulative
frequency in some years is well below 100%, at
RPRL=0.5.
In the North Sea survey, saithe are found mainly in
the north and northwest part of the survey area, west of
the Norwegian Trench. Saithe in the North Sea mainly
stay within the bottom channel (10 m above seabed)
(Figure 15a,b) and Table 4.
Blue Whiting
Along the Norwegian coast in areas deeper than 300 m
blue whiting usually have a relatively high acoustic
abundance and the distribution stretches over more than
one cod bottom-restricted FVR in some years (Figure
14f). In the total samples (Figure 13e,f) and samples in
areas shallower than 300 m (Figure 14e and Table 4),
however, the relative cumulative sA-values mainly
reaches 100% well below the RPRL=0.5, with reference
to the bottom pressure.
Blue whiting samples in the North Sea are few in
number and there are only night samples in 1998. The
median curves of most years, except 1999, have relative
cumulative acoustic density profiles within or stretching
just above RPRL=0.5, with reference to the bottom
pressure (Figure 15c,d). For 1999 (August–September)
the relative cumulative frequency indicates a mixture of
pelagic and demersal living. Samples from deep areas
(>300 m) are mainly from the area north of Shetland.
The annual median curves show that effective accumulation starts higher in the water column and reaches
100% at a level higher than RPRL=0.5 (Figure 15e,f),
indicating pelagic living in mid-water. The reason for the
S-shape of the median curve may be that a large part of
the echo abundance comes from a distribution range
stretching over roughly one FVR in mid-water. In
samples from 1995, 1996, 1998, the vertical extent of a
distribution range has a relative pressure reduction at
the upper end, with reference to the pressure at the
deeper end at approximately 0.5.
703
For samples of spawning blue whiting from surveys
west of the British Isles, the relative cumulative vertical
profiles of sA-values show that most blue whiting adapt
to pelagic living. The main distribution range gradually
leaves the seabed for samples with depths greater than
300 m and is all well above seabed at depths greater than
600 m (Figure 16f). Thus the relative vertical profiles are
accumulated from, and the RPRL=0 refers to, the
deepest end of each sample distribution range. The
distribution of both years’ profiles shows that the
main concentration of the stock is within RPRL=0.5
(Figure 16).
Physostome species
Capelin and herring have open swimbladders and thus
the relative vertical profiles of acoustic density are
accumulated from surface down to bottom. The depth
axis is partitioned as a relative proportion of depth from
sea surface to bottom.
Capelin
Data from winter surveys show that the capelin’s
relative vertical profile varies greatly over the survey
area (Figure 17a,b,c). High acoustic abundance is found
along the Polar Front in mid-water and mainly in areas
with a bottom depth more than 200 m. There are only a
few samples during summer and all of the fish are then
found close to the surface.
Herring
Samples from the North Sea show that herring are
mainly found in the upper half of the water column. The
vertical profiles of samples along the Norwegian coast
vary greatly. The influx of herring into the Barents Sea is
variable. In most years the data are insufficient for
conclusions to be made about diel vertical migration. In
winter, herring in the Barents Sea are mainly distributed
in the bottom half of the water column (Figure 17d,e,f),
and samples from year 2000 have diel vertical migration
with night-time ascents.
Temperature and salinity
Temperature, salinity, and sA-values from the winter
and summer surveys in the Barents Sea have been
interpolated onto the same locations. The distribution of
temperature and salinity, weighted by sA-values for each
survey and species, are presented in Tables 7–8.
Summer and winter the vertical profiles for cod,
haddock, and redfish show that in most cases the higher
concentration is below 100 m where the vertical temperature gradient is low. The survey coverage of the
eastern parts varies, and in 1997–1998 the survey vessels
did not enter the Russian zone. This may affect the range
704
B. K. Stensholt et al.
Table 7. Temperature and salinity weighted by acoustic sA-values for cod, haddock, and redfish.
Summer surveys in the Barents Sea, 1996–2000. All years and annual average. s.d. is standard
deviation; P1, P5, . . . are the percentiles.
Species
Mean
s.d.
P75
Median
P25
P5
3.30
5.77
4.39
4.32
9.07
3.06
3.40
3.11
2.10
8.56
1.99
3.18
2.22
1.26
7.38
0.21
1.93
1.54
0.11
3.80
3.6
6.5
4.8
7.0
5.6
3.3
3.9
3.0
4.0
3.6
3.16
2.40
2.27
3.32
2.59
2.7
1.0
1.1
2.5
1.1
1.0
0.3
0.5
1.0
0.2
0.4
0.5
1.1
0.4
0.8
Haddock summer temperature, Barents Sea
1996
3.30
0.97
6.6
6.2
1997
5.69
0.78
8.1
7.6
1998
5.93
0.61
8.8
8.1
1999
5.26
0.92
8.7
8.3
2000
5.19
0.60
8.3
7.2
3.3
6.6
7.0
7.2
6.1
3.2
5.8
6.1
4.9
5.4
3.2
5.0
5.1
3.6
4.5
1.8
3.0
2.8
2.3
2.3
1.0
2.1
1.8
0.9
0.8
Redfish summer temperature, Barents Sea
1996
3.81
0.88
6.5
1997
2.99
0.36
6.2
1998
3.58
0.53
6.2
1999
3.54
0.70
6.6
2000
3.10
0.69
6.2
6.4
5.2
6.1
6.5
5.4
5.5
3.7
5.2
4.1
4.0
3.9
2.8
3.1
3.1
3.0
1.9
2.4
2.5
2.6
1.9
1.2
1.0
1.2
1.6
1.6
1.0
0.1
0.6
0.8
1.6
35.02
35.01
35.11
35.05
35.04
34.94
34.89
35.05
35.01
34.92
34.87
34.83
35.02
34.97
34.66
34.76
34.55
35.00
34.90
34.38
34.39
34.22
34.92
34.71
34.11
34.22
34.13
34.71
34.35
33.90
Summer temperature 1996–2000,
Cod
2.76
1.45
Haddock
4.33
1.03
Redfish
3.41
0.68
Capelin
2.80
2.57
Herring
7.98
2.01
P99
P95
Barents Sea
7.13
5.20
8.34
7.58
6.49
6.23
8.09
6.74
10.39
9.73
Cod summer temperature, Barents Sea
1996
2.86
1.51
5.4
1997
2.62
1.76
7.8
1998
2.10
1.32
6.2
1999
3.36
1.52
8.2
2000
2.45
0.87
6.4
Summer salinity 1996–2000, Barents Sea
Cod
34.82
0.20
35.06
Haddock
34.70
0.15
35.05
Redfish
35.02
0.03
35.12
Capelin
34.93
0.17
35.07
Herring
34.64
0.36
35.07
of temperature distribution and also the relationship
between temperature and fish distribution. Most of the
cod, haddock and redfish are found in the range 2–6C
but a significant proportion of cod is in subzero temperatures where no haddock and very little redfish are
found. The subzero records for cod stem from its high
concentration along the Polar Front. For all species the
temperature variance is smaller in winter than in
summer, which may simply reflect the fact that fish
spread more out in summer. For cod this is in agreement
with the distribution of temperature and depth in
storage-tag data (Table 1).
The variations in salinity are, to our knowledge, too
small to make a significant physiological difference
(Tables 7, 8). The particularly low salinity variation for
redfish may reflect that it does not migrate as much
between different water masses as the other species. The
average temperature where the fish were found can vary
from year to year, the clearest difference is seen between
the median winter temperatures of the water masses for
redfish between 1996 and 2000.
P1
0.57
0.98
0.86
0.78
2.71
Discussion
Vertical distribution and relative pressure
reduction
The rapid vertical movements of a physoclist are limited
by the relative change in the ambient pressure because of
swimbladder physiology. Therefore, in the case of
demersal physoclists, it is meaningful to investigate the
relative vertical profiles of sA-values in terms of relative
pressure reduction with reference to bottom pressure
and to measure their extent in FVR units. This allows
the comparison of relative profiles from locations with
different bottom depth. Then the factors affecting the
vertical distribution can be investigated.
Each acoustic sample is a ‘‘snapshot’’ of the fish
distribution in the water column expressed as an integral
over 1 nmi, each fish staying in its own individually
adaptable range. The relative cumulative frequency of
the observed sA-values in each water column represents
the profile of the sample. If all individual FVR include
Vertical density distributions of fish
705
Table 8. Temperature and salinity weighted by acoustic sA-values for cod, haddock, and redfish.
Winter surveys in the Barents Sea, 1996–2000. All years and annual average. s.d. is the standard
deviation; P1, P5, . . . are the percentiles.
Species
P95
P75
Median
P25
P5
P1
Winter temperature 1996–2000, Barents Sea
Cod
3.10
1.28
6.8
5.8
Haddock
4.12
0.63
6.8
5.9
Redfish
3.66
0.62
6.8
6.1
Capelin
1.25
1.51
5.6
4.0
Herring
3.23
0.75
6.8
5.7
4.3
4.6
4.7
2.3
4.1
3.17
4.09
3.70
1.46
3.32
1.9
3.6
2.4
0.02
2.2
0.1
2.5
1.8
1.2
1.2
0.5
1.8
1.3
1.3
0.4
Cod winter temperature, Barents
1996
2.33
1.11
1997
2.97
1.00
1998
2.97
1.37
1999
3.13
1.18
2000
4.10
1.33
5.3
4.5
5.3
6.1
6.6
3.3
4.0
4.3
4.5
5.2
2.3
3.4
3.1
3.2
4.2
1.3
2.0
1.8
1.8
3.0
0.3
0.4
0.2
0.3
1.5
0.5
0.4
0.4
1.1
0.03
Haddock winter temperature, Barents Sea
1996
3.81
0.64
6.2
1997
4.00
0.49
5.7
1998
4.21
0.43
6.4
1999
4.17
0.60
6.9
2000
4.54
0.78
7.0
5.9
4.6
5.8
6.4
6.6
4.4
4.4
4.7
5.0
5.2
3.6
4.0
4.2
4.2
4.5
2.9
3.7
3.7
3.3
3.7
2.3
3.0
2.6
2.2
2.8
1.6
2.1
1.7
1.5
1.9
Redfish winter temperature, Barents Sea
1996
3.23
0.60
6.2
1997
3.28
0.59
5.5
1998
3.85
0.73
6.5
1999
4.19
0.37
6.9
2000
4.87
0.47
7.0
5.6
5.0
6.3
6.4
6.8
4.2
4.0
4.7
5.1
5.9
2.72
3.58
4.04
4.32
5.08
2.1
2.3
2.7
3.3
3.8
1.7
1.6
1.8
1.6
2.4
1.0
1.3
1.4
1.3
2.3
35.06
35.06
35.08
35.04
35.07
35.00
35.00
35.03
35.01
34.93
34.96
34.96
35.00
34.94
34.73
34.89
34.83
34.96
34.87
34.63
34.64
34.56
34.76
34.56
34.38
34.41
34.24
34.34
34.50
34.21
Winter salinity
Cod
Haddock
Redfish
Capelin
Herring
Mean
s.d.
P99
Sea
6.0
5.4
6.3
6.9
7.0
1996–2000, Barents Sea
34.92
0.10
35.10
34.90
0.11
35.10
34.97
0.06
35.11
34.92
0.12
35.05
34.76
0.11
35.15
the seabed, the distribution of the relative cumulative
sA-values, with its lower quartile, will reach 100% within
one bottom-restricted FVR, i.e. at the RPRL where the
median and quartile curves approach 100%. When the
distribution of the relative cumulative frequency reaches
100% well beyond one bottom-restricted FVR, e.g. when
the inter-quartile range remains large at an RPRL where
the median curve approaches 100%, it indicates that
some fraction of the stock has adapted to pelagic living.
The free vertical range and the extent of the
vertical distribution
The cod’s vertical distribution is related to the physiological restrictions of the individual cod’s vertical movements. This was recognized by Ellis (1956), who
reported from acoustic observations at the Skolpen
Bank in the Barents Sea, September 1955, that cod
shoals were found to be well within the theoretical limits
imposed by the swimbladder. The cod’s FVR was determined experimentally by Harden Jones and Scholes
(1981, 1985) as a depth interval with a pressure ratio of
0.5 between the upper and the lower end. It is used here
as a unit for comparison although other physoclists may
have a somewhat different FVR. The restriction to an
FVR is not absolute, as cod, saithe and haddock tolerate
60–70% pressure reduction from the pressure at neutral
buoyancy level before the swimbladder bursts (Tytler
and Blaxter, 1973).
The relative cumulative sA-values often reach 100% at
RPRL=0.5, which indicates that the experimentally
determined FVR is in fact a reasonable description of
the cod’s natural behaviour. The restriction to change of
pressure during rapid vertical movement, imposed by
swimbladder physiology, is also evident in the analysis
of acoustic samples of other physoclists. Their free
vertical ranges seem to be of roughly the same size as for
cod. This restriction is shown both in physoclists living
as demersals (Figures 4–14) and in pelagic distribution
layers (Figures 9e,f, 15f, 16).
Consider the deepest end of the vertical distribution
range of a sample, at depth M, to be the deepest end of
an FVR. Then the neutral buoyancy level, at depth d, is
the level where the relative pressure reduction is 1/3, and
706
B. K. Stensholt et al.
the upper end of FVR, at depth m, is the level where the
relative pressure reduction is 1/2. In most cases the
relative cumulative values are approximately 80% at
RPRL=1/3 (Figures 4–16). If the bottom-restricted
FVR is their natural FVR these results indicate that
demersal physoclistous fish spend most of their time
being under-buoyant. This supports the hypothesis of
Harden Jones and Scholes (1981), that if the overall
strategy is to conserve energy, such behaviour would be
adopted by demersal physoclists, while pelagic physoclists might maintain neutral buoyancy throughout
the day. In the case of blue whiting shoals, they appear
to be mainly within one FVR whether they distribute above seabed or as pelagic shoals in mid-water
(Figures 15d,f, 16).
which was mainly along the Polar Front and around
Svalbard, including a steep shelf edge. Acoustic signals
show demersal fish shoals that are distributed at the
seabed on a steep shelf and then extending outwards in
mid-water over the edge.
There is great diversity of prey species for cod,
haddock and redfish during the summer surveys, while
capelin is the main prey during the winter. The variation
in spatial distribution and behaviour among different
prey species may result in a large variation of the vertical
distribution of the predators among the different
samples collected from different locations (Neilson and
Perry, 1990). The depth distribution of individual cod
from data-storage tags has larger variation in summer
than in winter (Table 1, Stensholt, 2001) in agreement
with these results.
The distribution of vertical profiles
The FVR unit expresses the capacity of the physoclists
to perform rapid vertical moves but the fish may adapt
to make use of this capacity to a variable degree. This
degree is reflected in the vertical profile of the relative
cumulative acoustic density. Investigating such profiles
may reveal factors that influence the vertical migration
behaviour.
The median profile represents most of the samples.
Small samples of fish may have vertical profiles that are
modified due to special factors, e.g. the prey species’
vertical distribution, depth, density, and current
(Neilson and Perry, 1990). Thus the median of the entire
set of samples may not be a good representative of these
special cases. In these cases the analysis of such samples
should be carried out separately. Our results indicate
that great bottom depth, high density, and strong
currents affect the vertical distribution in all studied
physoclist species.
Because of the survey design only one sample is
available in each geographic location. Classification of
the locations for uniform environmental conditions,
which may affect the vertical profiles, could improve the
analysis. On the other hand such grouping of the data
means a reduced number of samples within each group
and so over-refinement could also reduce reliability.
Strong local environmental gradients in the saithe
surveys along the Norwegian coast lead to some
uncertainty with respect to the geographical separation
of day-night effects from environmental effects. The
larger variation in the extent of vertical distribution
observed during summer surveys than winter surveys
could be caused partly by the variation over different
years in the summer survey coverage, partly by fish
migration due to climate change (Midttun, 1965),
and partly by the vertical distribution being location
dependent. The distribution in summer 1999 shows more
pelagic living for cod and haddock than in most years.
This may be due to the survey line-transect coverage,
Evidence for pelagic living
The relative sA-values from the bottom channel (10 m
above seabed) generally decrease with increasing bottom
depth (Tables 3–4). The relative cumulative profiles
(Figures 6, 8–15) show more detail. In particular, the
lower quartile sometimes reaches 100% beyond one
bottom-restricted FVR, which is evidence that there are
fish above the bottom-restricted FVR in 25% of the
samples.
Neilson and Perry (1990) summarized from several
reports evidence that predators’ migration behaviour is
modified or governed according to food availability and
behaviour of their main prey species. Physoclists must be
adapted to pelagic living in order to forage more than
one bottom-restricted FVR above sea bottom. Such
adaptation means that the fish have an individual free
vertical range that does not include the seabed. After a
descent to the bottom these fish will be excessively
negatively buoyant staying below the lower free vertical
range limit. They compensate by gas secretion but the
process consumes time and energy (Harden Jones, 1951;
Blaxter and Tytler, 1972; Harden Jones and Scholes,
1985). Generally, with higher fish density and deeper
water, a larger fraction of the fish is found to adapt to
pelagic living. With increased density there is increased
competition for food, and this may cause more individuals to search for prey relatively higher in the water
column.
With increasing bottom depth there is an increasing
number of FVR-units from the seabed to the surface.
For physoclist fish adapted to stay near the bottom this
means that the upper layers become less accessible via
rapid ascents. With increased bottom depth fish adapted
to pelagic living can more easily reach the upper layers.
There are clear differences between the species. Among
the physoclists studied here blue whiting is most often
seen to be pelagic, especially in areas deeper than 300 m,
while redfish almost always stay within the bottom
Vertical density distributions of fish
half of the water column. The observations that the
physoclists in the demersal fish surveys stay within one
bottom-restricted free vertical range, may come from the
fact that these surveys are mainly in shelf areas of depths
<500 m. Moreover, the allocation of acoustic density to
species is based mainly on bottom trawl catches. Therefore the species composition found is not considered
reliable in the pelagic channel. In the surveys west of the
British Isles, blue whiting is the target species and the
allocation of acoustic values in the pelagic zone is much
better supported by pelagic trawl catches (Monstad et al.
1995, 1996). The pelagic distribution layers of blue
whiting are mostly within the limit of RPRL=0.5 (one
FVR for cod), when the RPRL is calculated based on
the deepest of the distribution layers. Similar results are
found for oceanic redfish. Observations from other
surveys report that occasionally cod, redfish and blue
whiting are found to distribute in or above mid-water
over great depths (Harden Jones and Scholes, 1985;
Magnússon et al., 1996; Monstad et al. 1995, 1996).
Records from data storage tags attached on cod
released in the Barents Sea (Stensholt, 2001) and the
North Sea (Righton et al., 2001) show that occasionally
individual cod have their daily range larger than the free
vertical range. In these exceptional cases the cod adapts
to pelagic living and seems to accept excessive negative
buoyancy at the bottom part of the daily range. Righton
et al. (2001) found that this occurs more often at shallow
depths. This is natural as in shallower waters the fish has
a shorter interval of excessive negative buoyancy to pass.
Even though the tags from the Barents Sea show that
this occasionally happens also during the periods of
summer and winter surveys, the acoustic profiles
demonstrate that the cod are mainly distributed well
within the bottom-restricted FVR in shallow depths.
Such a large daily range may be too rare to change the
general picture of the distribution of relative cumulative
acoustic profiles. However, in shallow areas (<200 m),
the profile of relative acoustic density may be distorted
or altered by avoidance reactions among fish in the
upper layers due to the approaching survey vessel (Ona
and Godø, 1990). Studies of the cod, haddock, and
saithe show that avoidance reactions are not always
evident and it is still unclear whether the avoidance
reaction is predominantly vertical or horizontal (Aglen,
1994; Godø et al., 1999; Totland, pers. comm.). Recent
studies have documented cases where no avoidance
could be detected (Fernandes et al., 2000a). These
observations were made with RV ‘‘Scotia’’ which is far
less noisy than the vessel used in our study (Fernandes
et al., 2000b).
Tidal currents
Most saithe found in the North Sea stay, for the most
part, within 10 m from the bottom (Figure 15a,b) and
707
Table 4. They are found mainly in the northern and
northeastern part of the survey area, along the west side
of the Norwegian Trench. Saithe along the Norwegian
coast north of 62N, are clearly distributed higher in the
water column (Figures 13a,b, 14a,b). The differing
behaviour in the two survey areas is most likely due to a
difference in the environmental conditions with the
strong tidal currents in the North Sea as a possible
cause. A similar observation was made in the cod’s
acoustic abundance on the Svalbard bank. Analysis of
storage-tag data shows examples of cod with reduced
vertical activity in connection with tidal currents
(Stensholt, 2001), and the saithe may stay close to the
bottom for the same reason.
Age composition, vertical movement, and
dead-zone loss
In most winter surveys the vertical profiles show that
during night-time a higher proportion of fish is distributed both near the bottom and further up into the
water column respectively than during daytime (Table 3,
Figures 5, 6d, 8d, 9d). During daytime there is a higher
proportion of fish around mid-range of the bottom half
of the water column. The degree of separation and
replacement of one part with another part, which results
in the observed day-night vertical profiles, may depend
on the size composition of the stock.
Several studies indicate that large fish ascend during
daytime and descend to the bottom during night-time,
while the opposite holds for small fish (Beamish, 1966;
Aglen et al., 1999; Korsbrekke and Nakken, 1999;
Einarsson, 2001; Stensholt and Nakken, 2001). The
storage-tag data also confirm a daytime ascent of
50–100 m for large cod during winter (Stensholt, 2001).
In 1996 there is higher relative cumulative acoustic
density near the bottom during day than night and
higher relative cumulative acoustic density at night
compared with day from the upper part of the vertical
distribution range (Figure 7). This may be due to the
very high proportion of one- and two-year-old cod in the
stock composition (Table 6). Smaller fish may migrate to
a depth level with a lower concentration of larger fish to
avoid predation. In 1996 the survey transect has the
most extensive coverage in the eastern part of the
Barents Sea, the main living area for young cod. In
1997–1999 there are no surveys east of 40E. In 1992–
1997 the 0-group index was relatively high for cod and
intermediate for haddock (Anon., 1998).
Table 5 shows very high night-time corrections in
2001. A possible cause is that too much of the acoustic
abundance from cod, haddock, and redfish was ascribed
to blue whiting, which were unusually abundant in the
survey area due, it is thought, to relatively high temperatures. There is a high proportion of one-year-old cod in
2001 and, if the proportion that ascends at night has
708
B. K. Stensholt et al.
been recorded as blue whiting, the acoustic loss at night
may have been correspondingly over-estimated, resulting in too high correction of the night-time acoustic
dead zone loss.
Logarithm of pressure ratio
With bottom depth D metres, a depth of x metres is
converted to RPRL=(Dx) · (D+10) 1. One may
regard the FVR as a logarithmic yardstick and convert a
depth of x metres to the logarithm of the pressure
ratio, i.e. log[(10+D) · (10+x) 1]. For small Dx,
which is usual for demersal fish, the RPRL is roughly
proportional to the difference of the logarithms:
ln[10+D]ln[10+x]= ln[1(Dx) · (D+10) 1]]
(Dx) · (D+10) 1.
Thus it makes no practical difference whether depths
are converted to RPRL or to logarithms of pressure
ratio. However, one FVR unit for cod corresponds
to x=0.5 · D5 and the pressure ratio is
(10+D) · (10+x) 1 =2, and the logarithm base 2 of the
pressure ratio is 1.
Conclusions
In order to compare the vertical profiles from the
overall acoustic line-transect survey data, large spatial
variations in absolute values among samples are
normalized using the relative acoustic sA-values. The
samples with unequal bottom depths are normalized by
expressing the relative vertical profiles of each acoustic
sample in terms of the relative pressure reduction level
(RPRL) from the seabed, or from the deepest end of the
vertical distribution range, up to the surface. The
normalized profile is suitable for a discussion of the
physiological limitation to change of pressure during
rapid vertical movement, the free vertical range (FVR),
among demersal physoclists. The observed acoustic vertical distribution of cod is generally within one FVR,
RPRL ranging from 0 to 0.5. Other physoclists show
similar distributions, which indicate that they have
roughly the same FVR as cod. For blue whiting, the
extent of observed pelagic distribution layer is mainly
within one cod FVR.
Taking into account the physiological limitation, the
fish vertical distribution is shown to change with
changes in factors such as age, fish density, prey species’
vertical migration patterns, day–night, bottom depth,
pressure change, strong current, and vessel noise. In
deep water or areas of high fish density, cod, haddock,
saithe, and blue whiting adapt to pelagic living in the
sense that their FVR does not include the seabed. In
areas with strong current or in shallow waters the fish
distribution is closer to the seabed and the vertical
distribution range is mainly less than the FVR.
Due to unequal day and night losses of acoustic
sA-values the profiles have been adjusted, and a method
for adjustment is suggested. The adjusted relative
acoustic vertical profiles for day–night suggest a net
vertical movement, while different parts of the stock may
have different diel migration. Discussion in terms of
stock composition, in the case of cod, indicates that this
is linked to changes in the age composition. When
vertical profiles are expressed in terms of FVR, then
trawl, acoustic and storage-tag information can be combined. A structure of spatial density distribution of fish
in relation to the environment and physiological limitation can be established. This information can be used
to improve the accuracy of stock estimates.
Acknowledgements
Special thanks to colleagues at the Institute of Marine
Research, Bergen, Norway, who gave excellent help on
different aspects: Kaare Hansen for organizing the
acoustic data, Jari Jacobsen and Bjørn Helge Bjørnestad
at the IT department, Odd Nakken for comment and
suggestion, Egil Ona for information on swimbladder
physiology, Adnan Ajiad, Tore Jacobsen, Terje
Monstad, Kjell H. Nedreaas, Odd Smestad, Atle
Totland. Thanks are also offered to editor John
Ramster and referees for their interest and comments.
References
Aglen, A. 1994. Sources of error in acoustic estimation of fish
abundance. In Marine Fish behaviour in Capture and
Abundance Estimation, pp. 107–133. Ed. by A. Fernö, and
S. Olsen. Fishing News Books, Oxford, UK.
Aglen, A. 1999. Report on demersal fish surveys in the Barents
Sea and Svalbard area during summer/autumn 1996 and
1997. Fisken og Havet, nr.7, 46 pp.
Aglen, A. 2001a. Report on demersal fish surveys in the Barents
Sea and Svalbard area during summer/autumn 1998 and
1999.
http://ressurs.imr.no/bunnfisk/raporter/sommertokt98.htm. 41 pp.
Aglen, A. 2001b. Results from the joint Russian-Norwegian
demersal fish survey in the Barents Sea, Winter 2001. Working Document AFWG, Bergen, 24 April–3 May 2001, 29 pp.
Aglen, A. 2002. Results from the joint Norwegian-Russian
demersal fish survey in the Barents Sea, February 2002.
Working Document AFWG, Copenhagen, 16–26 April 2002,
29 pp.
Aglen, A., Drevetnyak, K., Jakobsen, T., Korsbrekke, K.,
Lepesevich, Y., Mehl, S., Nakken, O., and Nedreaas, K. H.
2001. Investigations on demersal fish in the Barents Sea
winter 2000. Detailed report. IMR-PINRO Joint Report
Series no. 5, Bergen, Murmansk. 74 pp.
Aglen, A., Engås, A., Huse, I., Michalsen, K., and Stensholt, B.
1999. How vertical fish distribution may affect survey results.
ICES Journal of marine Science, 56: 345–360.
Vertical density distributions of fish
Aglen, A., and Nakken, O. 1996. Bunnfiskundersøkelser i
Barentshavet og ved Svalbard August 1995. In Norwegian.
Fisken og Havet, nr. 30, 33 pp.
Alexander, R. McN. 1966. Physical aspects of swimbladder
function. Biological reviews of the Cambridge Philosophical
Society, 41: 141–176.
Alexander, R. McN. 1970. Swimbladder gas secretion and
energy expenditure in vertically migrating fishes. In Proceedings of an International Symposium on Biological Sound
Scattering in the Ocean, pp. 75–86. Ed. by G. B. Farquhar.
US. Gov. Print. Off., Washington, D.C.
Alexander, R. McN. 1972. The energetics of vertical migration
by fishes. Symposia of the Society for Experimental Biology,
26: 273–294.
Anon. ICES Preliminary Report of the International Saithe
Survey in October–November 1992–1996.
Anon. Preliminary report of the international 0-group
fish survey in the Barents Sea and adjacent waters in
August–September 1998. 26 pp.
Anon. 1999. Internal cruise report ‘‘R/V Michael Sars’’ survey
No. 1999110. In Norwegian.
Arnold, G. P., and Greer Walker, M. 1992. Vertical movements
of cod (Gadus morhua L.) in the open sea and the hydrostatic
function of the swimbladder. ICES Journal of marine
Science, 49: 357–372.
Beamish, F. W. H. 1966. Vertical migration by demersal fish in
the Northwest Atlantic. Journal of Fisheries Research Board,
Canada, 23: 109–139.
Björnsson, B., Steinarsson, A., and Oddgeirsson, M. 2001.
Optimal temperature for growth and feed conversion of
immature cod (Gadus morhua L.). ICES Journal of Marine
Science, 8: 29–38.
Blaxter, J. H. S., and Tytler, P. 1972. Pressure discrimination in
Teleost fish: The effects of Pressure on organisms. Symposia
of the society for experimental biology, 26. The University
Press, Cambridge. 417–443.
Blaxter, J. H. S., and Tytler, P. 1978. Physiology and function
of the swimbladder. Advances in Comparative Physiology
and Biochemistry, 7: 311–367.
Brawn, V. M. 1962. Physical properties and hydrostatic function of the swimbladder of herring (Clupea harengus L.).
Journal of Fisheries Research Board, Canada, 19(4):
635–656.
Cressie, N. A. C. 1991. Statistics for Spatial Data. John Wiley
& Sons. 900 pp.
Dalen, J., and Nedreaas, K. H. 2001. Acoustic abundance
estimation and target strength observations of redfish
(Sebastes Mentella) in the Irminger Sea. Institute of Marine
Research, Bergen, Norway. 13 pp.
Einarsson, H. A. 2001. Length related diurnal vertical migration of cod (Gadus morhua L.), haddock (Melanogrammus
aeglefinus L.) and redfish (Sebastes spp.) in the Barents Sea.
Dept. of Fisheries and Marine Biology, University of Bergen.
Ellis, G. H. 1956. Observations on the shoaling behaviour of
cod (Gadus Gallarias) in deep water relative to daylight.
Journal of the Marine Biological Association of the UK, 35:
415–417.
Fernandes, P. G., Brierley, A. S., Simmonds, E. J., Millard,
N. W., McPhail, S. D., Armstrong, F., Stevemson, P., and
Sqires, M. 2000a. Oceanography: Fish do not avoid survey
vessels. Nature, 404(6773): 35–36.
Fernandes, P. G., Brierley, A. S., Simmonds, E. J, Millard,
N. W., McPhail, S. D., Armstrong, F., Stevemson, P., and
Sqires, M. 2000b. Oceanography: Fish do not avoid survey
vessels. Nature, 404(6801): 152.
Foote, K. G., Knudsen, H. P., Korneliussen, R. I., Nordbö,
P. E., and Röang, K. 1991. Postprocessing system for echo
709
sounder data. The Journal of the Acoustical Society of
America, 90: 37–47.
Godø, O. R., Somerton, D., and Totland, A. 1999. Fish
behaviour during sampling as observed from free floating
buoys – Application for bottom trawl survey assessment.
ICES CM 1999/J: 10.
Harden Jones, F. R. 1949. The teleostean swimbladder and the
vertical migration. Nature, 164: 847.
Harden Jones, F. R. 1951. The swimbladder and the vertical
movements of Teleostean fishes, I. Physical factors. Journal
of Experimental Biology, 28: 553–566.
Harden Jones, F. R. 1952. The swimbladder and the vertical
movements of Teleostean fishes, II. The restriction to rapid
and slow movements. Journal of Experimental Biology, 29:
94–109.
Harden Jones, F. R., and Scholes, P. 1974. The effect of
low temperature on cod, Gadus morhua. Journal du
Conseil International pour l’Exploration de la Mer, 35(3):
258–271.
Harden Jones, F., and Scholes, P. 1981. The swimbladder,
vertical movements, and the target strength of fish. In
Meeting on hydroacoustical methods for the estimation of
marine fish populations, 25–29 June 1979. II: Contributed
papers, discussion and comments. The Charles Stark
Draper Laboratory, Inc., Cambridge, Massachusetts, USA.
964 pp.
Harden Jones, F. R., and Scholes, P. 1985. Gas secretion and
resorption in the swimbladder of cod Gadus morhua. Journal
of Comparative Physiology, 155b: 319–331.
ISATIS. 1997. Transvalor and Geovariance Avon France.
Knudsen, H. P. 1990. The Bergen Echo Integrator: An introduction. Journal du Conseil International pour l’Exploration
de la Mer, 47: 167–174.
Korsbrekke, K., and Mehl, S. 2000. Abundance of saithe
Finnmiark-Møre autumn 2000. Int. rep. Institute of Marine
Research, Bergen, Norway, 20pp. In Norwegian with English
legends.
Korsbrekke, K., and Nakken, O. 1999. Length and species
dependent diurnal variation of catch rates in the Norwegian
Barents Sea bottom-trawl surveys. ICES Journal of Marine
Science, 56: 284–291.
Kristiansen, T. S., Michalsen, K., Jacobsen, J. A., and Huse, I.
2001. Optimal selection of temperature areas by juvenile cod
(Gadus morhua L.) in the Barents Sea modeled by dynamic
optimization. ICES Journal of Marine Science, 58: 172–182.
Kutchai, H., and Steen, J. B. 1971. The permeability of the
swimbladder. Comparative Biochemistry and Physiology,
39A: 119–123.
MacLennan, D. N., Fernandes, P. G., and Dalen, J. 2001. A
consistent approach to definitions and symbols in fisheries
acoustics. ICES Journal of Marine Science (in press).
Magnússon, J., Magnússon, J. V., Sigurðsson, Þ., Reynisson,
P., Hammer, C., Bethke, E., Pedchenko, A., Gavrilov, E.,
Melnikov, S., Antsilerov, M., and Kiseleva, V. 1996. Report
on the Joint Icelandic/German/Russian Survey on Oceanic
Redfish in the Irminger Sea and Adjacent Waters in June/
July 1996. ICES CM 1996/G: 8 Ref. H.
Marshall, N. B. 1972. Swimbladder organization and depth
ranges of deep-sea teleosts. Symposia of the Society for
Experimental Biology, 26: 261–272.
Mehl, S. 1997. Botnfiskundersøkingar I Barentshavet (Norsk
sone) vinteren 1997. Fisken og Havet, nr. 11. 72 pp.
Mehl, S. 1998. Botnfiskundersøkingar I Barentshavet (redusert
område) vinteren 1998. Fisken og Havet, nr. 7. 69 pp.
Mehl, S. 1999. Botnfiskundersøkingar I Barentshavet vinteren
1999. Fisken og Havet, nr. 13, 70 pp.
Mehl, S., and Nakken, O. 1996. Botnfiskundersøkingar I
Barentshavet vinteren 1996. Fisken og Havet, nr. 11. 68 pp.
710
B. K. Stensholt et al.
Michalsen, K., Aglen, A., and Pennington, M. 2000. Report on
demersal fish surveys in the Barents Sea and Svalbard area
during summer/autumn 2000. Internal report Institute of
Marine Research, Bergen, Norway. 41 pp.
Midttun, L. S. 1965. The relation between temperature
conditions and fish distribution in the Southeastern Barents
Sea. ICNAF Special Publication, 6: 213–220.
Monstad, T., Belikov, S. V., Shamrai, E. A., and McFadzen,
I. R. B. 1995. Investigations on Blue Whiting in the area west
of the British Isles, spring 1995. ICES CM 1995/H: 7 Pelagic
fish committee, Ref. B, 21 pp.
Monstad, T., Belikov, S. V., and Shamrai, E. A. 1996. Report
of the joint Norwegian-Russian acoustic survey on Blue
Whiting during spring 1996. ICES CM 1996/H: 12 Pelagic
fish committee, 23 pp.
Nedreaas, K. H. 1995. Annual acoustic surveys on 2-5 year old
North-East Arctic saithe in October. ICES Saithe study
group, Aberdeen 30 May–2 June 1995, 8 pp.
Neilson, J. D., and Perry, R.I. 1990. Diel vertical migrations of
marine fishes: an obligate or facultative process? Advances in
Marine Biology, 26: 115–168.
Ona, E., and Godø, O. R. 1990. Fish reaction to trawling
noise: the significance for trawl sampling. Rapports et
Procès-Verbaux des Réunions du Conseil International pour
l’Exploration de la Mer, 189: 159–166.
Ona, E., and Mitson, R. B. 1996. Acoustic sampling and signal
processing near the seabed: the deadzone revisited. ICES
Journal of Marine Science, 53: 677–690.
Righton, D., Metcalfe, J. D., and Arnold, G. P. 2001. Vertical
reality: utilizing knowledge of cod behaviour to interpret
survey results. ICES Annual Science Conference 2001, CM
2001/Q: 20, 15 pp.
Rose, G. A., and Porter, D. R. 1996. Target-strength studies on
Atlantic cod (Gadus morhua) in Newfoundland waters.
Simmond, E. J., and MacLennan, D. N. (eds.) Proceedings of
an ICES International Symposium on Fisheries and Plankton Acoustics, Aberdeen, Scotland, 12–16 June 1995. ICES
Journal of Marine Science, 53(2): 259–265.
Ross, L. G. 1976. The permeability to oxygen of the swimbladder of the meso-pelagic fish Ceratoscopelus maderensis.
Marine Biology (Berl.), 37: 83–87.
Ross, L. G. 1979. The permeability to oxygen and the guanin
content of the swimbladder of a physoclist fish, Pollachius
viren. Journal of the Marine Biological Association of the
UK, 59: 437–441.
Stensholt, B. K. 2001. Cod migration patterns in relation to
temperature: analysis of storage tag data. ICES Journal of
Marine Science, 58: 770–793.
Stensholt, B., and Nakken, O. 2001. Environmental factors,
spatial density and size distributions of 0-group fish. In
Proceeding of the symposium on Spatial Processes and
Management of Marine Populations, Oct 27–30, 1999,
pp. 395–413. Ed. by G. H. Kruse et al. Anchorage Alaska,
USA.
Stensholt, B. K., and Sunnanå, K. 1996. Spatial distributions of
variables in marine environ-mental and fisheries research:
Geostatistics and autocorrelated environmental and fisheries
data. ICES Annual Science Conference, Reykjavik, Iceland
27 Sept–1 Oct 1996. 27 pp.
Sundnes, G., and Bratland, P. 1972. Notes on the gas content
and neutral buoyancy in physostome fish. Fiskeridirektoratets skrifter, Serie Havundersøkelser, 16: 89–97.
Tytler, P., and Blaxter, J. H. S. 1973. Adaptation by cod and
saithe to pressure changes. Netherlands Journal of Sea
Research, 7: 31–45.
Tytler, P., and Blaxter, J. H. S. 1977. The effect of swimbladder
deflation on pressure sensitivity in the saithe Pollachius
Virens. Journal of the Marine Biological Association of the
UK, 57: 1057–1064.
Woodhead, P. M. J., and Woodhead, A. D. 1959. The effects of
low temperatures on the physiology and distribution of the
cod, Gadus morhua L., in the Barents sea. Proceeding of the
Zoological Society of London: 181–199.
Woodhead, P. M. J., and Woodhead, A. D. 1965. Seasonal
changes in the physiology of the Barents Sea cod, Gadus
morhua L., in relation to its environmient. II. Physiological
reactions to low temperature. ICNAF Special Publications,
6: 717–734.