- Research article
- Open Access
Model-supported estimation of mortality rates in Baltic cod (Gadus morhua callarias L.) larvae: the varying impact of 'critical periods'
© Voss et al; licensee BioMed Central Ltd. 2001
Received: 20 September 2001
Accepted: 19 November 2001
Published: 19 November 2001
Changes in the survival-rate during the larval phase may strongly influence the recruitment level in marine fish species. During the larval phase different 'critical periods' are discussed, e.g. the hatching period and the first-feeding period. No such information was available for the Baltic cod stock, a commercially important stock showing reproduction failure during the last years. We calculated field-based mortality rates for larval Baltic cod during these phases using basin-wide abundance estimates from two consecutive surveys. Survey information was corrected by three dimensional hydrodynamic model runs.
The corrections applied for transport were of variable impact, depending on the prevailing circulation patterns. Especially at high wind forcing scenarios, abundance estimates have the potential to be biased without accounting for transport processes. In May 1988 mortality between hatch and first feeding amounted to approximately 20% per day. Mortality rates during the onset of feeding were considerably lower with only 7% per day. In August 1991 the situation was vice versa: Extremely low mortality rates of 0.08% per day were calculated between hatch and first feeding, while the period between the onset of feeding to the state of an established feeder was more critical with mortality rates of 22% per day.
Mortality rates during the different proposed 'critical periods' were found to be highly variable. Survival rates of Baltic cod are not only influenced by a single 'critical period', but can be limited at different points during the larval phase, depending on several biotic and abiotic factors.
Due to the missing data on larval mortality in the field, it was not possible to investigate the influence of different processes on larval survival rates. Moreover, it was not possible to test the applicability of general ideas of larval survival for the Baltic Sea, as data on stage specific mortality rates are missing. Especially the concept of a 'critical period' connected to the onset of feeding [11, 12], but also the impact of hatching mortality have to be considered.
In this approach we use basin-wide abundance estimates from two consecutive surveys to calculate mortality rates. Abundance estimates are corrected for transport processes using the initial horizontal distribution of larvae and applying a 3 dimensional hydrodynamic model. Mortality rates within two periods, from hatch (egg stage IV) to first feeding (larval stage 5–7) and from onset of feeding to the state of an established feeder (larval stage 8) were calculated. We compared two sampling periods, May 1988 and August 1991.
The resulting mortality estimates are discussed with respect to available information on possible sources of mortality, as these are suggesting a high variability in the impact of 'critical periods' on larval survival.
Initial horizontal distributions
Hydrodynamic model runs
Cod larval mortality rates:
May 1988 Time interval between surveys: 5.66 days;
August 1991 Time interval between surveys: 5.92 days;
age [days] *
The applied correction for transport processes had a varying impact. Especially for high speed, non-circular transport regimes these advection processes have to be taken into account, although the principle differences in mortality rates can also be seen without a corresponding correction (see Table 1). It is well known that the flow dynamics in the Baltic Sea are highly complex. They are mainly determined by the ephemeral character of wind stress, the baroclinic mass field and the complicated bottom topography . However, the Baltic Sea Model used for the simulations has generally proven to be applicable for simulating egg and larval drift processes [8, 14]. To explicitly test the applicability of the model for predicting the horizontal distribution patterns of ichthyoplankton during the investigated time periods, separate model runs were performed during an earlier study . The simulated results of this study were in good accordance with field observations, confirming the applicability of the Baltic Sea Model. This good accordance could be reached, although, as in the present study, simplistic vertical distribution patterns had to be used. The vertical position of the eggs and larvae does have an influence on their horizontal transport. Unfortunately, no direct measurements for the time periods under investigation are available. While stage-specific vertical distribution patterns of cod eggs are well investigated [6, 7, 15], corresponding information for cod larvae is rather scarce. The available information on stage-specific vertical distribution of cod larvae in the Baltic [5, 16] can only be used to extract rough distribution patterns. If older cod larvae perform diurnal vertical migrations, as observed in other cod stocks  can so far not be answered. As a first approximation of the vertical distribution vertically resolving samplings 4 weeks apart from the periods under investigation were used. It is known that the vertical distribution patterns can change over the spawning period, but strong changes within 4 weeks are not to be expected. Furthermore, a sensitivity analysis conducted for the transport processes in the Bornholm Basin  revealed no significant differences in the drift patterns of cod larvae within a vertical range of 12 m. As a conclusion, we are quite confident, that the modelled development of the horizontal distribution patterns are, despite the various shortcomings, close to reality.
The importance of corrections for transport processes will even increase for planned studies on sprat larval mortality, as these are distributed shallower in the water column [15, 19] and therefore are subject to higher current speeds.
There is no indication for a strong small-scale patchiness of cod eggs or larvae, which would not be resolved by the station grid and thereby considerably bias the abundance estimates. As well the results from the 15 years ichthyoplankton time series as results from high-resolution samplers like the Longhurst-Hardy Plankton Recorder (Raab, pers. comm., Inst. Marine Sciences, Kiel) and the Ichthyoplankton Recorder (Mees, pers. comm., Inst. Marine Sciences, Kiel) rather led to the assumption, that the distribution and abundance of ichthyoplankton is well covered by the station grid in use. The uncertainties in the abundance estimates arising from the non-synoptic sampling procedure as well as the spatial sampling resolution can so far not be quantified, but a theoretical study dealing with this problem has been conducted by Voss and Hinrichsen (in prep.) showing mean errors in abundance estimates of 10–20%.
Even after accounting for these uncertainties, the results presented here for cod larvae suggest a high variability in the importance of different 'critical periods' for recruitment. In May 1988 the period between hatch and the onset of feeding was characterised by high mortalities (~20 % day-1), whereas during the same period in August 1991 very low mortality rates were observed. Different processes have been proposed to influence survival on this stage. Viable hatch requires a minimum oxygen level of 2 ml/l . To account for year to year variations in the hydrographic environment occupied by the eggs, the so called Reproductive Volume (RV) has been defined [21, 22]. The RV defines the volume of water suitable of successful development of cod eggs by threshold levels in temperature (>1.5°C), salinity (>11 psu) and oxygen content (> 2 ml/l) as experimental studies have shown that below these limits development of cod eggs ceases before hatch [23, 24, 20]. In August 1991 the RV in the Bornholm Basin was 127 km3, substantially higher than in May 1988 with only 93 km3. Unfortunately, no field sampling for the vertical distribution of eggs was performed during the two periods, but these results as well as modelling results  suggest less favourable hydrographic conditions for the late egg and early larval stages in May 1988. Even if hatching occurs, low oxygen concentrations negatively impact larval activity . Under such circumstances hatched cod larvae might sink to deeper, less oxygenated water layers, as specific gravity increases after hatch and the larvae are not able to counterbalance the sinking rate. Another potentially important source of mortality is predation by the clupeids sprat and herring. Performed stomach content analysis of sprat and herring indicate considerable predation pressure on all egg stages in May 1988 , while for August 1991 the consumption of cod eggs by the clupeids was low . This may have the potential to explain the strong differences in mortality rates during the period from hatch to first feeding encountered in this study.
The second investigated 'critical period' ranges from the onset of feeding (larval stage 5) to the state of an established feeder (larval stage 8). During this phase the larvae have to migrate vertically into upper water layers with sufficient light conditions and higher prey concentrations for successful feeding . Predation on these larvae is of less importance, as the vertical overlap between prey and predator is rather limited, with the clupeids feeding within and below the permanent halocline [26, 27]. Correspondingly, in May 1988 almost no larvae were identified in the diet of herring and sprat. In a time series principally covering the spawning seasons from 1988 to 1996, only in August 1991 relatively high numbers of cod larvae were found in the stomachs of herring , indicating an increased predation pressure. Another possible source of variation in survival are differences in prey availability for the larvae. Unfortunately, there are no data on abundance and distribution of suitable prey organisms, like copepod nauplii and small copepodite stages [28–30], for these time periods. The only indication available is that the total meso-zooplankton biomass in the Bornholm Basin was higher in August 1991 than in May 1988 . However, it has been suggested, that prey abundance in the Bornholm Basin is always sufficient to ensure feeding success of cod larvae . Therefore the potential influence of changing prey availability remains to be clarified.
The ability of the larvae to migrate vertically and to establish as feeding larvae might be influenced by their energy reserves, i.e. the size and quality of their yolk sac. Egg size is a function of fish size and the batch number spawned . In the beginning and at the end of the spawning season the egg size is smallest. While it was main spawning time in May 1988, August 1991 corresponds to the late spawning season . There are no direct egg diameter measurements available, but due to the relative spawning time egg diameter and therefore yolk reserves should have been higher in May 1988, favouring survival of the larvae.
The different mortality rates found for the two investigated periods can not easily be attributed as inter-annual or seasonal differences. Due to the complex abiotic and biotic relationships in the Baltic Sea, it is rather a combination of both effects. However, the two time periods, May 1988 and August 1991 revealed the peak larval abundance in each year, so that changes in mortality rates probably had a comparable impact on the resulting recruitment.
Baltic cod recruitment decreased strongly in the beginning of the 1980ies. Since 1985 it fluctuated around a rather low level compared to the late 1970's and early 1980's. Both investigation dates exhibited outstanding high abundance values for cod larvae in the last 15 years . While the spawning stock biomass as well as the potential egg production were higher in 1988 than in 1991 , recruitment in 1988 reached only an intermediate level (compared to the 'low recruitment phase' since 1985). The 1991 year-class was above average, but by far also not reaching the values from the beginning of the 1980ies. It becomes obvious that also after reaching the state of an established feeder processes must act as bottlenecks of survival, limiting year-class strength. Besides others, transport processes have been shown to influence survival probability [35, 8].
Transport processes have the potential to strongly influence abundance estimates, especially for cod larvae. These processes have therefore to be accounted for when calculating mortality rates.
Mortality rates during the different proposed 'critical periods' were found to be highly variable. Survival rates of Baltic cod are not only influenced by a single 'critical period', but can be limited at different points during the larval phase, depending on several biotic and abiotic factors. It becomes further obvious that also after reaching the state of an established feeder processes must act as bottlenecks of survival, limiting year-class strength. Field based mortality estimates, as presented here, may help to disentangle the influence of the various stage-specific processes under consideration and to better understand the mechanisms regulating year-class strength.
Materials and Methods
- 25% in 27–30 m, 50% in 30–33 m and 25% in 33–36 m for feeding cod larvae (stages > 5)
- 25% in 66–69 m, 50% in 69–72 m and 25% in 72–75 m for cod eggs and non-feeding larvae
Estimation of mortality rates
Mortality rates were estimated by comparing the daily production of a given developmental stage at the first sampling date with that of a corresponding stage at the second sampling date (cohort method). Temperature-dependent incubation time of eggs were calculated  with an estimated ambient mean temperature . As the consecutive surveys were approximately 6 days apart, the developmental stages had to be chosen accordingly. Mean stage-specific ages as well as exact time intervals between successive samplings are given in Tab. 1. As the time interval between the sampling dates did not match exactly the difference in age between developmental stages, the daily production observed at the second sampling date was adjusted . For the young (hatching) larvae mortality rates were calculated from egg stage IV (oldest egg stage) to larval stage L5-7 (first feeding). For the older larvae larval stages 5–7 were compared with larval stage 8 (established feeding).
While the station grids covered the area > 60 m depth, all calculations were based on the mean abundance (n/m2) of the developmental stages inside the 80 m depth contour line. This was done, as abundance values outside the observation area were unknown and therefore no reliable exchange rates due to transport could be calculated. By restricting the calculations to the area > 80 m depth, reliable exchange rates between the areas 60–80 m depth and > 80 m depth could be derived from a hydrodynamic model.
The abundance values were then corrected for the amount of transport gains or losses as obtained by the model runs.
Correction of mortality rates by application of a hydrodynamic model
As the two surveys were ~6 days apart, changes in abundance due to transport processes had to be taken into account. Therefore, the Bornholm Basin was divided in 3 sub-areas: >80 m, 60–80 m and <60 m. While the initial stage specific abundance is known from the survey results for the areas >60 m depth, the abundance for the area <60 m depth was set to be 0. The initial horizontal distributions as obtained by the first surveys were entered into the model as a tracer and the exchange rates between the areas were calculated over time. By this, the abundance found during the second survey in the area > 80 m depth could be corrected for transport gains or losses. The numerical simulations of the Bornholm Basin circulation were performed by application of a three-dimensional (3-D) eddy resolving baroclinic model of the Baltic Sea . The Baltic Sea Model is based on the free surface Bryan-Cox-Semtner model  which is a special version of the Cox numerical ocean general circulation model [41–43]. The Baltic Sea model comprises the whole Baltic with a horizontal resolution of 5 km and 41 vertical levels specified. For the region of the Bornholm Basin this results in a vertical resolution of 3 m layers. A grid size of 5 km and a time step of 5 min were chosen. Within the Bornholm Basin the model was initialized with three-dimensional hydrographic data (temperature and salinity) obtained during the research surveys. No hydrographic data were taken outside the observational area for running the model and in order to overcome this lack of data, the general features of the Baltic were utilized by incorporation of hydrographic characteristics typical for these regions and time periods obtained from previous model runs. For each depth level of the model, observational data were interpolated onto the model grid by objective analysis . The model was forced for all simulations with actual wind data for the entire Baltic provided by the SMHI (Swedish Meteorological and Hydrological Institute). In order to adapt the initial fields to the model dynamics and to the prescribed mass field outside the Bornholm Basin the model was allowed to spin-up for 4 days without external forcing prior to incorporation of tracers representing the horizontal distribution of cod eggs and larvae. After this initialization period, the forcing was switched on and the model was run for a period of 10 days.
The two periods (19–28 May 88 and 11–20 Aug 91) exhibit contrasting wind forcing scenarios: In May 1988 wind forcing was relatively low in magnitude and variable in direction. In contrast, the simulation for August 1991 was affected by relatively strong wind forcing of mainly western directions.
To visualize the persistent circulation patterns averaged model results for different depth layers in the Bornholm Basin will be presented. However, averaged currents give no information about their variability. Thus, we calculated the stability which is defined as the ratio of the averaged vectorial velocity and the averaged arithmetic velocity.
We thank Lonny Hansen of the Danish Meteorological Institute for making available the wind data from the weather station at Christiansø. The work was carried out as part of the BASYS project and of the Baltic STORE project and was funded by a research grant from the European Union (FAIR CT 98–3959).
- Dickson R, Brander K: Effects of a changing windfield on cod stocks of the North Atlantic. Fish. Oceanogr. 1993, 2: 124-153.View ArticleGoogle Scholar
- ICES: Report of the Baltic Fisheries Assessment Working Group. ICES CM. ACFM: 18-Google Scholar
- Bagge O, Thurow F: The Baltic cod stock, fluctuations and the possible causes. ICES mar. Sci. Symp. 1994, 198: 254-268.Google Scholar
- Wieland K, Jarre-Teichmann A: Prediction of vertical distribution and ambient development temperature of Baltic cod (Gadus morhua L.) eggs. Fish. Oceanogr. 1997, 6(3): 172-183. 10.1046/j.1365-2419.1997.00038.x.View ArticleGoogle Scholar
- Grønkjær P, Wieland K: Ontogenetic and environmental effects on vertical distribution of cod larvae in the Bornholm Basin, Baltic Sea. Mar. Ecol. Prog. Ser. 1997, 154: 91-105.View ArticleGoogle Scholar
- Wieland K: Einfluβ der Hydrographie auf die Vertikalverteilung und Sterblichkeit der Eier des Ostseedorsches (Gadus morhua callarias) im Bornholmbecken, südliche zentrale Ostsee. Ber. Inst. f. Meeresk. Kiel. 1995, 266: 114-Google Scholar
- Wieland K, Hinrichsen HH, Grønkjær P: Stage-specific mortality of Baltic cod (Gadus morhua L.) eggs. J. Appl. Ichthyol. 2000, 16: 266-272. 10.1046/j.1439-0426.2000.00241.x.View ArticleGoogle Scholar
- Voss R, Hinrichsen HH, St John MA: Variations in the drift of larval cod (Gadus morhua L.) in the Baltic Sea: combining field observations and modelling. Fish. Oceanogr. 1999, 8: 199-211. 10.1046/j.1365-2419.1999.00106.x.View ArticleGoogle Scholar
- Fortier IC, Leggett WC: A drift study of larval fish survival. Mar. Ecol. Prog. Ser. 1985, 25: 245-257.View ArticleGoogle Scholar
- Heath MR, MacLachlan P: Dispersion and mortality of yolk-sac herring (Clupea harengus L.) larvae from a spawning ground to the west of the Outer Hebrides. J. Plankton Res. 1987, 9: 613-630.View ArticleGoogle Scholar
- Hjort J: Fluctuations in the year classes of important food fishes. J. Cons. Perm. int. Explor. Mer. 1914, 164: 73-76.Google Scholar
- Cushing DH: Plankton production and year class strength in fish populations: an update of the match/mismatch hypothesis. Ad. Mar. Biol. 1990, 26: 249-293.View ArticleGoogle Scholar
- Lehmann A: A three-dimensional baroclinic eddy-resolving model of the Baltic Sea. Tellus,. 1995, 47A: 1013-1031.View ArticleGoogle Scholar
- Hinrichsen HH, Lehmann A, St John MA, Brügge B: Modeling the cod larvae drift in the Bornholm Basin in summer 1994. Cont. Shelf Res. 1997, 17: 1765-1784. 10.1016/S0278-4343(97)00045-9.View ArticleGoogle Scholar
- Wieland K, Zuzarte F: Vertical distribution of cod and sprat eggs and larvae in the Bornholm Basin (Baltic Sea) 1987–1990. ICES CM. 1991, J: 37-Google Scholar
- Grønkjær P, Clemmensen C, St John MA: Nutritional condition and vertical distribution of Baltic cod larvae. J. Fish Biol. 1997, 51: 352-369.View ArticleGoogle Scholar
- Lough RG, Potter DC: Vertical distribution patterns and diel migrations of larval and juvenile haddock Melanogrammus aeglefinus and Atlantic cod Gadus morhua on Georges Bank. Fish. Bull US. 1993, 91: 281-303.Google Scholar
- Voss R: Horizontale Verteilung und Drift von Dorschlarven in der südlichen zentralen Ostsee in Abhängigkeit vom mesoskaligen Strömungssystem. Diploma Thesis University Kiel. 1996, 112-Google Scholar
- Makarchouk A, Hinrichsen HH: The vertical distribution of ichthyoplankton in relation to the hydrographic conditions in the Eastern Baltic. ICES CM. 1998, R: 11-Google Scholar
- Wieland K, Waller U, Schnack D: Development of Baltic cod eggs at different levels of temperature and oxygen content. Dana. 1994, 10: 163-177.Google Scholar
- M Plikshs, M Kalejs, G Graumann: The influence of environmental conditions and spawning stock size on the year-class strength of eastern Baltic cod. ICES CM. 1993, J: 22-Google Scholar
- MacKenzie BR, St John MA, Wieland K: Eastern Baltic cod: perspectives from existing data on processes affecting growth and survival of eggs and larvae. Mar. Ecol. Prog. Ser. 1996, 134: 265-281.View ArticleGoogle Scholar
- Thompson BM, Riley JD: Egg and larval developmental studies in the North Sea cod (Gadus morhua L.). Rapp. P.-v. Réun. Cons. int. Explor. Mer. 1981, 178: 553-559.Google Scholar
- Nissling A, Westin L: Egg mortality and hatching rate of Baltic cod (Gadus morhua) in different salinities. Mar. Biol. 1991, 111: 29-32.View ArticleGoogle Scholar
- Rohlf N: Verhaltensänderungen der Larven des Ostseedorsches (Gadus morhua) während der Dottersackphase. Ber Inst. Meeresk. Kiel. 1999, 312: 54-Google Scholar
- Köster FW: Der Einfluβ von Bruträubern auf die Sterblichkeit früher Jugendstadien des Dorsches (Gadus morhua) und der Sprotte (Sprattus sprattus) in der zentralen Ostsee. Ber. Inst. f. Meeresk., Kiel. 1994, 261: 286-Google Scholar
- Köster FW, Möllmann C: Trophodynamic control on recruitment success in Baltic cod?. ICES J. Mar. Sci. 2001, 57: 310-323. 10.1006/jmsc.1999.0528.View ArticleGoogle Scholar
- Last JM: The food af three species of gadoid larvae in the eastern English Channel and southern North Sea. Mar. Biol.,. 1978, 48: 377-386.View ArticleGoogle Scholar
- van der Meeren T, Næss T: How does cod (Gadus morhua) cope with variability in feeding conditions during early life stages?. Mar. Biol. 1993, 116: 637-647.View ArticleGoogle Scholar
- Fossum P, Ellertsen B: Gut content analysis of first-feeding cod larvae (Gadus morhua L.) sampled at Lofoten, Norway, 1979–1986. ICES mar. Sci. Symp. 1994, 198: 430-437.Google Scholar
- HELCOM: Third periodic Assessment of the State of the Marine Environment of the Baltic Sea, 1989–93. Baltic Sea Environ. Proc. 64B,. 1996Google Scholar
- Krajewska-Soltys A, Lingkowski TB: Densities of potential prey for cod larvae in deep-water basins of the southern Baltic. ICES CM,. 1994, J: 17-Google Scholar
- Vallin L, Nissling A: Maternal effects on egg size and egg buoyancy of Baltic cod, Gadus morhua. Implications for stock structure effects on recruitment. Fish. Res. 2000, 49: 21-37. 10.1016/S0165-7836(00)00194-6.View ArticleGoogle Scholar
- Köster FW, Möllmann C, Neuenfeldt S, Plikshs M, Voss R: Developing Baltic cod recruitment models I. Resolving spatial and temporal dynamics of spawning stock and recruitment for cod, herring, and sprat. Can. J. Fish. Aquat. Sci.,. 2001, 58: 1516-1533. 10.1139/cjfas-58-8-1516.View ArticleGoogle Scholar
- Hinrichsen HH, St John MA, Aro E, Grønkjær P, Voss R: Testing the larval drift hypothesis in the Baltic Sea: Retention vs. dispersion due to the influence of the wind driven circulation. ICES Mar. Sci. Symp. 2001,Google Scholar
- Westernhagen Hv: Erbrütung der Eier von Dorsch (Gadus morhua), Flunder (Pleuronectes flesus) und der Scholle(Pleuronectes platessa) unter kombinierten Temperatur- und Salzgehaltsbedingungen. Helgoländer wiss. Meeresunters. 1970, 21: 21-102.View ArticleGoogle Scholar
- Wieland K: Distribution and mortality of cod eggs in the Bornholm Basin (Baltic Sea) during two patch studies in 1986. Kieler Meeresforsch. Sonderh. 1988, 6: 331-340.Google Scholar
- Wieland K, Köster FW: Size and visibility of Baltic cod eggs with reference to size-selective and stage-dependent predation mortality. J. Appl. Ichthyol. 1996, 12: 83-89.View ArticleGoogle Scholar
- Fossum P: A staging system for larval cod. Fisk. Dir. Skr. Ser. Hav. Unders. 1986, 18: 69-76.Google Scholar
- Killworth PD, Stainforth D, Webbs DJ, Paterson SM: The development of a free-surface Bryan-Cox-Semtner ocean model. J. Phys. Oceanogr. 1991, 21: 1333-1348. 10.1175/1520-0485(1991)021<1333:TDOAFS>2.0.CO;2.View ArticleGoogle Scholar
- Bryan K: A numerical method for the study of the circulation of the world ocean. J. Phys. Oceanogr. 1969, 15: 1312-1324.Google Scholar
- Semtner AJ: A general circulation model for the World Ocean. UCLA Dept. of Meteorology Tech. Rep. 1974, 8: 99-Google Scholar
- Cox MD: A primitive equation 3-dimensional model of the ocean. GFDL/Princeton University. GFDL Ocean Group Tech. Rep. 1984, 1: 144.Google Scholar
- Hiller W, Kaese RH: Objective analysis of hydrographic data from mesoscale surveys. Ber. Inst. Meeresk. Kiel. 1983, 116: 78-Google Scholar
- G Neumann, WJ Pierson: Principles of Physical Oceanography. London, Prentice-Hall Int. Inc. 1967Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.