Inferences from bowhead whale ovarian and pregnancy

SC/56/BRG8
Inferences from bowhead whale ovarian and pregnancy data: age
estimates, length at sexual maturity and ovulation rates
J.C. GEORGE1, E. FOLLMANN2, J. ZEH3, M. SOUSA2, R. TARPLEY4, AND R. SUYDAM1
1. Department of Wildlife Management, Barrow, AK
2. University of Alaska Fairbanks, Fairbanks, AK
3. University of Washington, Department of Statistics, Box 354322, Seattle, WA, 98195-4322
4. School of Veterinary Medicine, St. George's University, Grenada, West Indies
Contact email: [email protected]
ABSTRACT
We used life history and reproductive data (number of corpora in both ovaries) to estimate mean ovulation rate,
length at sexual maturity, number of ova per estrous cycle, and age. For this analysis, we define ovulation as the
number of ‘readable’ corpora formed per year in both ovaries. We used data archive bh.corpora which included
data for 86 female bowhead whales, 75 with both ovaries examined. Of these 75, 35 were immature (0 corpora
albicantia (CA), 0 corpora lutea (CL)) and 40 were mature (at least 1 CA or CL). Corpora counted in a single ovary
ranged from 0 to 22, and 1 to 41 in both ovaries for mature females with both. Using logistic regression and
including data from bh.reproduction, we estimated length at maturity (females) to be 13.41 m. While no ovaries in
our sample had more than 1 CL, diameters of CAs in Tarpley and Hillmann (1999) suggest multiple ovulation is
possible, and was estimated to be ~1.13 corpora/estrous based on their data. To estimate ovulation rate we used data
from Barrow and Kaktovik where most whales are carefully examined for pregnancy and maturity. These examined
whales have known probability of being mature: 0 for immature and 1 for mature females. Since not all females are
examined for maturity, we estimated the probability that the unexamined females were mature as a function of body
length using data from bh.individual_whale and the logistic regression results. The denominator of the ovulation rate
estimate is the sum of the probabilities of being mature and the numerator is the number of females either pregnant
or having a CL. The estimated ovulation rate was 0.368/yr for an ovulation interval of 2.7yr. Estimates of ages for
the 40 mature females (with both ovaries available) were computed from CA counts, using the estimates of age at
sexual maturity, number of ova released per estrous and ovulation rate. Corpora age estimates are consistent with
those of George et al. (1999) and Rosa et al. (2004) using the aspartic acid racemization (AAR) technique. A
regression of AAR age on ‘corpora age’ results in an R2 value of 0.66. A problem with corpora age estimates is, of
course, that they cannot be calculated for males or immature whales. In contrast to AAR age estimates in George et
al. (1999) and Rosa et al. (2004), which lack very old females, these data suggest that some females may live nearly
as long as males. Of interest is whale 92B2 which was not used in this analysis since only one ovary was recovered,
but her single ovary held 20 corpora. This whale also had a stone implement lodged in the dorsal blubber, suggesting
she was very old (George et al., 1999). Estimated standard errors of the corpora ages are lower than those for AAR
ages. An attractive aspect of the corpora ageing technique is that corpora counts are relatively easy to do and there
is relatively little uncertainty in them. Our results suggest that CA counts give a reasonable estimate of age, past
pregnancies, and an upper bound for lifetime calf production.
KEYWORDS: BOWHEAD WHALE, REPRODUCTION, AGE ESTIMATION, OVULATION, CORPORA
1
SC/56/BRG8
INTRODUCTION
The female reproductive system of placental mammals is a highly derived character that is unique to Class
Mammalia. Many other vertebrates are oviparous (lay eggs) but only Prototherian mammals (e.g., the duck billed
platypus) do this following internal fertilization. Young mammals develop within the uterus, are nourished via the
placenta and are born in various states of development. In the case of cetaceans, calves are anatomically and
physiologically prepared to cope with the rigors of life in ice-covered seas (precocial). However, most neonate
mammals are altritial requiring considerable parental care (Vaughan et al., 2000).
In adult mammals, ova develop within ovarian follicles, and are released when the follicle ruptures during
estrus. The follicles then transform into a corpus luteum (CL). If pregnancy occurs, the CL maintains the
pregnancy through hormone secretions, principally progesterone. Following pregnancy, or if pregnancy fails, the
CL regresses into a corpus albicans (CA). The corpus albicans is composed of fibrous scar tissue (Tarpley and
Hillmann, 1999) and in some mammalian orders, such as Cetacea, these appear to be retained in the ovary leaving a
“record” of ovulation.
The reproductive cycle of the bowhead is unclear, however, some aspects are known. Based on fetus
recoveries, conception likely occurs in March with birth about 14 months later (Reese et al., 2001). A single calf is
born the following spring (late April to early June) with a mode in late May. The calf is presumably nursed for less
than a year based on examinations of whales believed to be a year old (ingutuks). Relatively short lactation (~9
months) is characteristic of mysticete whales (Costa and Williams, 1999). Calving intervals are believed to be
about 3-4 years in length (Miller et al., 1992; Rugh et al., 1992; Koski et al., 1993).
If corpora are persistent in the ovary, then they offer a record of ovulation for a female. Early researchers
examining large samples of whales taken commercially suggested they were, in fact, persistent (Tarpley and
Hillmann, 1999). However, additional problems of interpretation exist. Questions that remain unresolved for
whales are 1) whether ovulation is induced or spontaneous, 2) the number of ova released during an estrus cycle, and
3) how many ovulations do not lead to fertilization and/or implantation. Uterine scars from pregnancies do not form
in cetaceans because of the epitheliochorial type of placental attachment (Vaughn et al., 2000).
Tarpley and Hillmann (1999) carefully assessed the diameter, regression, and persistence of 257 corpora.
CAs ranged in diameter from 6.3 cm to 0.3 cm. Histograms of CA diameter indicated a normal distribution with a
rather abrupt size break at about 0.5 cm. This is either evidence of maximum regression or sudden absorption, and
conclusions from Tarpley and Hillmann (1999) were tentative. However, based on a qualitative analysis of the
distribution of 41 CAs from whale 89B3, a very large female, they suggested that corpora were persistent (at least
for this animal).
Nerini et al. (1984) conducted some of the early investigations of corpora counts in bowhead whale
ovaries. Not surprisingly, she did not find a strong correlation between whale length and corpora number, since the
variation in age at length for mature bowheads is considerable.
To estimate age from corpora counts, several types of data are required and assumptions must be made.
These include:
i.
Corpora are persistent, or disappear at some known rate. In this analysis we assume they are
persistent.
ii.
Multiple ovulations during estrous must be quantified. An estimate of how often estrous occurs must
be obtained.
iii.
The average age at sexual maturity must be estimated.
We discuss these assumptions and data requirements in the text at some detail.
Based on the above discussion, the basic objectives of this paper are to:
i)
assemble all corpora count data into a database by whale identification number
ii)
test the possibility of using corpora counts to estimate the age of bowhead whales.
METHODS
We used new and published reproductive and life history data to estimate several statistics. These include: the mean
ovulation rate, the number of corpora in both ovaries, length at sexual maturity, age at sexual maturity and number
of ova per ovulation event. Throughout this paper, the terms ‘mature’ and ‘maturity’ refer to sexual maturity.
Bowheads likely reach physical maturity at a greater length and older age than sexual maturity (George et al., 1999).
For this analysis, we define ovulation rate as the number of ‘readable’ corpora formed per year in both ovaries.
2
SC/56/BRG8
Corpora counts
We followed methods of Tarpley and Hillmann (1999) to prepare and ‘read’ the ovaries. After fixing in 10%
formalin, ovaries were sliced in 5 mm sections and the numbers of CLs and CAs were counted. The approximate
diameter of each CL was noted as well. Only whales for which both ovaries were examined were used to estimate
age, but the presence of a CL or CA in a single ovary defined a female as mature.
Length and age at sexual maturity
We used binary logistic regression to estimate the length at maturity. The data archives bh.corpora and
bh.reproduction provided data on females known to be either immature or mature on the basis of corpora counts,
ovary size or pregnancy. We scored whales as mature (1) based on the presence of corpora or a fetus or immature
(0) if these were absent. One whale (79WW1) was scored as mature in the absence of corpora data because her
ovaries were extremely large (13kg). If only one ovary was examined and no corpora were found, we did not
classify the whale as immature and include it in the analysis unless it was very small; e.g. 84S2 (14.4m) was not
classified, but 02B5 (8.5m) was treated as immature. The logistic regression model is:
⎛ p ⎞
⎟⎟ = β 0 + β 1 × BL
log⎜⎜
⎝1− p ⎠
where log is the natural logarithm, p is the probability that a female bowhead with body length BL is sexually
mature, and β0 and β1 are the parameters to be estimated. Given the parameters, the length LSM at which the
probability is 0.5 that the whale is mature can easily be estimated. Since p = 0.5, p/(1 – p) = 1 and log(1) = 0.
Therefore
LSM =
− β0
β1
.
The estimated probability that a female bowhead with body length BL is mature is
Pr(mature|BL) = exp( β 0 + β 1 × BL ) / [1 + exp( β 0 + β 1 × BL )].
To determine ASM, the mean age at sexual maturity, we compute the mean of ages estimated from aspartic
acid racemization data in the bh.aspartic_age data archive (George et al., 1999; Rosa et al., 2004). We include
females with lengths at which, according to Pr(mature|BL), there is a 2% or better chance that sexual maturity is
reached. However, we omit whales in this size range if they are known from corpora counts to be well past the age
at sexual maturity, and we adjust ages of whales with only one CA by subtracting the ovulation interval. The
ovulation interval (yr) is computed as 1/OR, where OR is the “ovulation rate” defined earlier. Note that our use of
the term ovulation rate differs slightly from other studies, and we assume that ovulations (and pregnancies) result in
a readable CA. We estimate V(ASM), the variance of ASM, by the sample variance of the ages from which ASM was
estimated.
Ovulation rate
A crude estimate of ovulation rate (yr-1) (OR) can be obtained by dividing the number of females with a CL or fetus
by the number of mature females. We used whales from fall harvests at Barrow and Kaktovik for this calculation
since whales from these villages were most consistently examined by biologists. In addition, there is a higher
probability of missing a small spring fetus than a larger fall fetus. Since not all females are examined for maturity,
we estimated Pr(mature|BL) as described in the previous section and used it in the following manner:
OR =
∑ pf
∑ mf
where the sums are over harvested potentially mature females; pf = 1 if the harvested female has a CL or fetus, pf =
0 otherwise; and mf = the probability that the harvested female is mature. We treat females 12m or more in length
as potentially mature; the shortest female found to date with a CL or fetus was 12.6m long (whale 99B6). The
3
SC/56/BRG8
longest female verified to be immature was 14.2m long. We therefore set mf = 1 for whales >15m long, as well as
for those known to be mature from corpora or pregnancy data. For females 12-15m long not known to be mature, mf
= Pr(mature|BL). Since OR can be viewed as the probability that a mature whale ovulates in a given year,
Pr(ovulation) = OR × Pr(mature), so this way of calculating mf gives the correct weight to potentially mature whales
whose status is uncertain, all of whom contribute pf = 0 to the numerator. We estimate the variance of OR by
V (OR) = OR × (1 − OR )
∑ mf .
Number of ova per estrous
Plots of CA diameters in Tarpley and Hillmann (1999) were examined to estimate the mean number of ova per
estrous cycle. If the largest and presumably most recent corpora from the same ovary had the same diameter, it was
assumed that they represented multiple ova from the same estrous n event. If there was a unique largest corpora, we
assumed that it represented a single ovulation for that year. In either case, if the mean number of ova released, AR,
is greater than one, then corpora counts must be adjusted downwards accordingly as described below. The variance
of AR was estimated as the square of the standard error of the mean.
Ages estimated from corpora for mature whales
The estimated age of the whale is
age = ASM + TCA×1/(AR×OR)
where ASM = mean age at sexual maturity, TCA = total CAs counted from both ovaries, AR = mean number of ova
released per ovulation event, and OR = ovulation rate. If corpora are counted without error, an estimated variance
of the age estimate (derived by the delta method) is:
V(age) = V(ASM) + [TCA/(AR×PR)2]2 × [PR2×V(AR) + AR2×V(PR) – V(AR)×V(PR)].
We use only CA in computing the age estimate. This is because we define maturity by the presence of
corpora or a fetus. So the year a female attains sexual maturity, she should have exactly 1 CL and/or early (not
term) fetus, and her age is estimated by ASM (Fig. 1). Now suppose she has 1 CL and 1 or 2 CAs. Then she has just
1 or 2 ovulation intervals to add to ASM. And so on. Now suppose she has no CL, but 1 CA. Then she reached
sexual maturity some fraction of an ovulation interval ago, so the first CA , if there is no CL, should only contribute
a fraction of an ovulation interval. However, to simplify, we give all CAs the same weight. As we will see below,
the uncertainty about the interval represented by the first CA is a minor source of uncertainty compared to other
sources accounted for in V(age).
We also assume that the corpora can be counted with negligible error. This assumption is discussed later.
Comparisons with aspartic acid age estimates
Some of the whales examined had age independently estimated using the aspartic acid racemization (AAR)
technique (George et al., 1999; Rosa et al., 2004). These age estimates were regressed against the number of
corpora as a means of estimating ovulation rates. We also compared AAR ages with the ages estimated from
corpora counts.
RESULTS and DISCUSSION
Corpora counts
The data archive bh.corpora included data on 86 female bowhead whales. Corpora counted in a single ovary ranged
from 0 to 22. There was no significant difference in number of corpora in the right versus left ovary among the
mature whales (P-value = 0.3, paired t-test), but the magnitude of the difference ranged from 0 to 9. Therefore,
except for the length at maturity analysis, we restricted our analyses to the 75 whales that had both ovaries examined
for corpora. Of these 75, 35 were immature (0 CA, 0 CL) and 40 were mature (at least 1 CA or CL).
The distributions of both CA and CA+CL for the mature whales with both ovaries examined were highly
skewed. Although 2 whales with only one ovary examined had counts of 20 and 22 CA respectively, only 4 of the
40 whales with both ovaries examined had 18 or more total corpora. Basic statistics are shown in Table 1 and the
distribution of total corpora in Fig. 2.
None of the ovaries had more than one CL.
4
SC/56/BRG8
Length and age at maturity
From the 86 females in the bh.corpora archive and 37 potentially pregnant females from the bh.reproduction
archive, we assembled a data set with 98 females known to be either mature or immature. The smallest known
mature female was whale 99B6, a 12.6m female early in her pregnancy. The largest female known to be immature
was whale 83B2, 14.2m. The 35 immatures included 6 whales between 13 and 13.9 m, all the rest were < 12.3 m.
The results of the logistic regression were β0= -32.840, β1 = 2.448 (Table 2). The estimated length at which
the probability is 0.5 that the whale is mature is LSM = 13.41m. Pr(mature|BL) in the range of lengths in which over
95% of the whales likely reach maturity ranges from 0.03 at BL = 12m to 0.98 at BL = 15m (Fig 3).
The increase in Pr(mature|BL) over Pr(mature|BL-0.1m) is 0.016 or more from BL = 12.4 through BL =
14.5 and is <0.015 outside that range, suggesting that around 2% or more of females reach maturity with each
increase of 0.1m in length within that range. This range, not surprisingly, corresponds well to the range within
which the smallest known mature and largest immature female fall. To estimate age at sexual maturity (ASM), we
therefore use ages determined by the aspartic acid racemization technique for females in this size range.
The only such whales in the bh.aspartic_age archive are 96B4, 96B6, 96B10, 99B6, 99B18 and 00B2. Of
these whales, we know from corpora counts that 96B4 (10CA) and 00B2 (5 CA) are well past ASM and therefore
should be omitted. The remaining whales can be used: 99B6 can reasonably be assumed to be newly mature, 99B18
has a CA and a CL so is at least one ovulation interval past sexual maturity, and the remaining two whales are of
unknown maturity status. Ages computed from the two D/L ratio measurements available for 99B18 differ as much
from each other as they do from the estimates for the other whales because of D/L ratio measurement error (Rosa et
al., 2004). We therefore use both age estimates for this whale along with one age for each of the three other whales
in estimating ASM and V(ASM). The mean of the ages is ASM = 24.4 with sample variance V(ASM) = 26.7,
corresponding to SD(ASM) = 5.2. Another method of estimating ASM (and its SD) is baleen ages (Lubetkin, et al
2004).
Estimated LSM = 13.41m is quite consistent with other estimates. Koski et al. (1993) estimated female
length at sexual maturity as 13-13.5m from photogrammetric data. Because of stretching when whales are landed
(George et al., 2004), we might expect estimates from harvest data to be longer. On the other hand,
photogrammetric estimates depend on seeing a cow with a calf, and estimates based on corpora counts do not
require a successful pregnancy, only an ovulation. It may be that miscarriages and/or neonatal mortality are higher
among the smallest/youngest mature females; hence smaller whales may be detected as mature in the harvest sample
than in the photogrammetry sample.
Ovulation rate and mean number of ova per ovulation event
In the bh.individual_whale archive, there were 39 females > 12m, harvested on known dates in the fall (August
through October) at Barrow or Kaktovik between 1976 and 2003. We restricted attention to those years because
only in those years was it reasonable to assume that all females landed were carefully examined for evidence of
ovulation and/or pregnancy. Two of these whales were known from corpora data to be immature and were not
considered further. Of the remaining 37, 21 were known to be mature from corpora or pregnancy data, 3 were
assumed mature because their lengths were >15m, and the remaining 13 had mf set to Pr(mature|BL) in the
denominator of OR. The number of mature females in the denominator was estimated to be 28.54. Among the 21
known mature whales, 10 were known to be pregnant and/or had a CL observed and 1 could not be conclusively
identified as pregnant or with a CL. In the numerator of OR, that whale, 76B20, was assigned pf = 0.5 instead of pf
= 1. The resulting estimated ovulation rate was OR = 10.5/28.54 = 0.368 with estimated V(PR) = 0.0081,
corresponding to SE(PR) = 0.090. The ovulation interval was estimated to be 1/0.368 = 2.7yr.
Based on examination of plots in Tarpley and Hillmann (1999), of the largest and probably most recent
corpora from the same ovary, multiple ovulations may have occurred in 4 of 38 ovulation events (38 ovaries). The
maximum was 3 CAs for whale 82WW1. The mean number of ova released per event was AR = 1.13. If this is
representative of multiple ovulations, then corpora counts must be adjusted downwards accordingly as described
below. The variance of AR was estimated to be V(AR) = 0.00451, corresponding to standard error SE(AR) = 0.067.
In some cases, CAs from different were close in size, so we may have underestimated OR to some extent.
Ovulation rates estimated from landed whales have several potential biases. Biases could result from: i)
failure to detect CL or early pregnancies, ii) difficulties in determining if females are mature, and iii) failure to
account for mature females that were NOT harvested because they were accompanied by a calf. These biases can be
minimized considerably if only fall animals are used – but not eliminated. In fall, fetus size has consistently been
about 1.5m which is relatively easy for a biologist or hunter to detect. Also, mothers do not always swim closely
5
SC/56/BRG8
with their calves or may have lost their calf, and, at times, may not be avoided by hunters. There are a couple of
examples in the fall hunt (Tarpley and Hillmann, 1999) in which lactating females were taken accidentally. On the
other hand, some calves appear to remain with their mothers for approximately 10 months and so, some avoidance
of mothers with calves is likely during the fall harvest as well. Also, mothers with calves avoid nearshore waters
less than 20 m deep (Koski and Miller, 2004) and that habitat preference may introduce a slight negative bias to the
fall harvest in the proportion of females with calves. This might produce a slight upward bias in the
ovulation/pregnancy rate estimate.
Corpora counts seemed relatively low in our dataset, with a median of 9 per animal. It could be that most
female bowheads in the BCBS population have had less than 10 calves in their lifetime. However, four females had
over 20 corpora suggesting that either some females are quite old, that fecundity is highly disproportionate between
females, or both.
There were no whales observed with more than one corpus luteum. This suggests either that the release of
multiple ova is extremely rare or that only one develops into a CL during pregnancy. Nerini et al., (1984) reported
that one bowhead had two active CLs, and we know of one observation by Eskimo hunters of twin calves; but no
bowheads have been observed with twin calves in any aerial or ice-based survey. Still, the possibility of multiple
ovulations within a year exists. This might happen if, say, an early pregnancy fails.
We have assumed that the corpora can be counted with negligible error. However, a small experiment with
two people counting the corpora from the same whale independently should be conducted to verify this.
Corpora age estimates
Estimates of ages for the 40 mature females were computed from CA counts, using the estimates given above for
ASM, AR and OR. They are listed in Table 3 with their standard errors and shown against body length in Fig. 4.
The corpora age estimates are consistent with those of George et al. (1999) and Rosa et al. (2004) using the
AAR technique. A problem with corpora age estimates is, of course, that they cannot be calculated for immature
whales. Since ASM was estimated from AAR ages, we expect positive correlation between corpora and AAR ages
among the younger females. However, the oldest female ages estimated from corpora data were >100yr, as were
AAR ages for the oldest males in Rosa et al. (2004). Whale 92B2 was not used in this analysis since only one
ovary was recovered, but her single ovary had 20 corpora, so the total could well be ~40. This whale also had a
stone implement lodged in the dorsal blubber, suggesting she was very old (George et al., 1999). In contrast to
AAR age estimates in George et al. (1999) and Rosa et al. (2004), which lack very old females, the corpora data
suggest that some females may live nearly as long as males.
Our age estimates are sensitive to the estimated ovulation rate and interval, which are intern sensitive to the
length at sexual maturity (LSM). Estimates of LSM will become more precise as more intermediated sized whales
are examined for maturity.
Correlations with aspartic acid ages
Results of a regression of AAR age on ‘corpora age’ for eight females with both are shown in Fig 5. This depends
on which age is used for whale 99B7. Using the L2 value of 27yr from Appendix A of Rosa et al. (2004), the R2
value is 0.66 between the two ageing techniques. Using an age of 71yr, the L1 value, gives an R2 of 0.21. By either
measure, our results suggest that corpora counts are useful for ageing bowhead whales. It should be noted that the
SEs are lower for the ovulation ages than AAR ages.
An attractive aspect of this ageing technique is that corpora counts are relatively easy to do and there is
very little uncertainty in the corpora counts. The calving interval, while still uncertain, may not introduce error of
the magnitude noted in the AAR method.
Regressing corpora counts on aspartic acid age estimates suggests an ovulation rate of 0.311/yr, which,
considering the large errors in the AAR ages, agrees reasonably well with our more direct estimate of 0.37.
Do corpora counts represent pregnancies? Most likely. Based on several lines of evidence: our results,
the polyandrous mating system of Balaenids, and the long periods between reproductive events, it seems likely that a
receptive female would become pregnant during estrous. Some females will undoubtedly fail to become pregnant;
however, the CA counts probably give a reasonable estimate of past pregnancies.
Do pregnancies represent calf production? Some portion of these pregnancies will undoubtedly fail;
nonetheless, CA counts give an upper bound for lifetime calf production. Furthermore, because the number of
corpora is relatively small, it makes sense that many of these must be successful births or the BCBS stock would not
be increasing at its current estimated rate.
6
SC/56/BRG8
ACKNOWLEDGEMENTS
We thank the Alaska Eskimo Whaling Commission (AEWC) and the Barrow Whaling Captains’ Association for
their confidence, guidance and support of our research. The postmortem examinations that we conduct on landed
whales often interfere with the butchering process. The whale hunters’ patience with us has been enormous. We
appreciate the field assistance of Dave Ramey, Victoria Woshner, Cheryl Rosa, Todd O’Hara, Matt Irinaga, Paul
Nader, Benny Akootchook, Cyd Hanns and many others; and the lab assistance of Thorsten Bentzen and Lara Dehn.
We recognize Mayor George N. Ahmaogak, Sr., and Dr. Thomas Albert, former Senior Scientist; Taqulik Hepa,
Deputy Director; and Charles D. N. Brower, Director of the North Slope Borough Department of Wildlife
Management for their steadfast support. We gratefully acknowledge funding provided by the North Slope Borough
Department of Wildlife Management and National Oceanic and Atmospheric Administration (through the AEWC).
Lisa Delarosa, Dolores Vinas, and April Brower provided office assistance.
REFERENCES
Costa, D.P, and Williams, T.M. 1999. Marine Mammal Energetics. Biology of Marine Mammals. Smithsonian.
George, J.C., Bada, J., Zeh, J., Scott, L., Brown, S., O’Hara, T. and Suydam, R. 1999. Age and growth estimates of
bowhead whales (Balaena mysticetus) via aspartic acid racemization. Canadian Journal of Zoology. 77:571-80.
George, C., Koski, B., Rugh, D. 2004b. Body stretching of bowhead whales during hauling and butchering during
the subsistence hunt. SC/56/BRG9.
Koski, W.R. and G.W. Miller. 2004. Habitat selection by different size classes of bowhead whales in the central
Beaufort Sea during late summer and autumn. (unpublished manuscript).
Koski, W.R, Davis, R.A., Miller, G.W. and Withrow, D.E. 1993. Reproduction. pp. 239-69. In: J.J. Burns, J. J.
Montague and C.J. Cowles (eds.) The Bowhead Whale. Special publication No. 2 of the Society of Marine
Mammalogy. i-xxxvi + 787pp.
Lubetkin, S., Zeh, J., Rosa, C., and George, J.C. 2004. Deriving von Bertalanffy age-length relationships for
bowhead whales (Balaena mysticetus) using a synthesis of age estimation techniques. SC/56/BRG3.
Miller, G. W., Davis, R.A., Koski, W.R., Crone, M.J., Rugh, D.J., Withrow, D.E. and Fraker, M.A. 1992. Calving
intervals of bowhead whales: an analysis of photographic data. Rep. int. Whal. Commn 42:501-6.
Nerini, M.K., Braham, H.W., Marquette, W.M. and Rugh, D.J. 1984. Life history of the bowhead whale, Balaena
mysticetus (Mammalia: Cetacea). J. Zool., London 204:443-68.
Reese, C.S., Calvin, J.A., George, J.C. and Tarpley, R.J. 2001. Estimation of fetal growth and gestation in bowhead
whales. J. Am. Stat. Assoc. 96:915-23.
Rosa, C., George, J.C., Zeh, J., O’Hara, T. Botta, O. and Bada, J. 2004. Update on age estimation of bowhead
whales using aspartic acid racemization. Paper SC/56/BRG6.
Rugh, D.J., Miller, G.W., Withrow, D.E. and Koski, W.R. 1992. Calving intervals of bowhead whales established
through photographic identifications. J. Mamm. 73 (3):487-490.
Tarpley, R.J. and Hillmann, D.J. 1999. Observations on ovary morphology, fetal size and functional correlates in
the bowhead whale Balaena mysticetus. Final Report to the Department of Wildlife Management, North Slope
Borough, Box 69, Barrow, AK.
Vaughan, T.A., Ryan J.M., Czaplewski, N.J. 2000. Saunders College Publishing, Orlando, FL. 565 pp. +
Appendices.
7
SC/56/BRG8
Table 1
Basic corpora count statistics for mature whales with both ovaries examined, N= 40. The mode for Total CA is NA
because counts of 1, 5, 6, 9 and 12 were each seen in 4 whales, other counts in fewer whales.
Statistic
Mean
SD
Minimum
Mode
25th percentile
Median
75th percentile
Maximum
Total CA
10.2
8.7
1
NA
5
9
12
41
CA+CL
10.7
8.7
2
2
5
9
12.5
41
Table 2
Parameters from logistic regression of maturity status on body length (BL) to estimate Pr(mature|BL).
Parameter
β0
β1
Value
-32.840
2.448
SE
9.404
0.683
T-value
-3.492
3.584
8
SC/56/BRG8
Table 3
Estimated “ovulation” ages for mature female bowhead whales using corpora data. Data are sorted by body length.
*99B7 age estimates from two D/L determinations were 27 and 71. WhaleID = the whale identification number,
TotCA = total CAs counted; TotCL = total CL counted, BL = body length (m), baleen = baleen length (cm), OA =
ovulation age, OA-se = standard error of the ovulation ages, AARage = aspartic acid racemization age (Rosa et al.,
2004), AARse = standard error of AAR age.
WhaleID TotCA TotCL BL baleen OA
99B18
1
1
13
342
27
97B8
2
0
13.6 260
29
93B17
5
0
13.6
36
01B17
2
0
13.9 255
29
92B7
1
1
14.2
27
81S2
5
1
14.2 344
36
76B20
6
0
14.3
39
96B4
10
0
14.4 300
48
00B2
5
0
14.5 273
36
92B10
5
0
14.6 320
36
00B3
11
1
14.6
51
92B3
15
0
14.6
60
89B2
6
1
14.7 280
39
99B16
1
1
14.8
27
88KK1
6
1
14.9 297
39
90B4
10
1
14.9 274
48
92B9
9
0
15
46
85WW2
12
1
15.1 320
53
95B8
1
1
15.2
27
99B7*
4
1
15.4
34
00B4
9
0
15.4
46
93B20
3
1
15.5
32
88G1
9
1
15.7 301
46
87B5
11
1
15.7 300
51
87B6
12
0
15.7 330
53
85WW1
6
0
16.2 355
39
81KK3
9
0
16.2 340
46
92B4
12
0
16.2 360
53
03B9
25
1
16.4 344
85
82WW2
7
0
16.5 280
41
81WW3
27
1
16.5
89
97B10
13
1
16.7 320
56
02B2
17
0
16.7 352
65
89B3
41
0
16.9 405
123
86KK2
12
0
17.2 356
53
81KK1
8
0
17.4 362
44
82WW1
16
0
17.6
63
86WW2
4
1
17.7 345
34
00B5
16
1
18.9 331
63
02B3
34
0
19.2 409
106
OA-se AARage
5
29
5
6
5
5
6
6
8
43
6
26
6
8
41
10
6
5
39
6
8
8
9
5
6
27
8
52
5
8
8
9
6
8
9
16
7
17
9
11
25
9
7
11
6
11
66
21
AARse
12
13
8
13
12
12
14
15
9
SC/56/BRG8
Year
0 Month
Type J F M
1
A MJ J A S O N D J
2
F M A M
J J A S O N D J
3
F M A MJ J A S O N D J
4
F M
A MJ J A S O N D J
5
F M
prep OC pregnancy
CL forms
3 yr prep 0C pregnancy
CL forms
B
lactation
CA transition
rest
CA regresses
4 yr prep 0C pregnancy
CL forms
B
lactation
CA transition
rest
CA regresses (in size)
B
A MJ J A S ON D J
F MA M
lactation
prep OC pregnancy
CL forms
B
Fig. 1. Diagram of possible reproductive cycles for bowhead whales. Rows 1 and 2 show details for possible 3year and 4-year calving intervals, respectively.
14
12
10
8
6
4
Std. Dev = 8.80
2
Mean = 10.7
N = 39.00
0
0.0
5.0
10.0
15.0
20.0
25.0
30.0
35.0
40.0
TOTCORP
Fig. 2. Distribution of total corpora for 40 bowhead whales in which both ovaries were examined. Note that nearly
all females had less than 18 corpora but 4 (10%) had 25 or more.
10
SC/56/BRG8
AA
AA
A
A
AAA
A
AA
A
A
A A
A
AA
A
A
A
A
A
A
A
A
A
A
Predicted probability
1.00000
AA
A
0.75000
A
A
A
0.50000
A
A
A
0.25000
A
A
A
A
A AA
A
A
A
A
A
AA
A
AA
0.00000
7.5
A
A
AA
A
10.0
12.5
15.0
17.5
bodylen
Fig. 3. Plot showing the results of the logistic regression of bowhead body length and maturity.
20
19
Body Length (m)
18
17
16
15
14
13
12
20
40
60
80
100
120
140
Ovulation Age (yr)
Fig. 4. Plot of ovulation age by length. These data suggest that some females grow quite old (> 100 yr) which was
only the case for males reported in George et al. (1999) and Rosa et al. (2004).
11
SC/56/BRG8
80
70
y = 0.8828x + 3.7389
R2 = 0.6645
AAR age (yr)
60
50
40
30
20
10
0
0
10
20
30
40
50
60
70
80
Ovulation Age (yr)
Fig. 5. Plot of aspartic acid racemization age estimates and corpora counts. Note that the lower age estimate (27 yr)
for whale 99B7 was used in this plot. With this small sample, there appears to be good agreement between these
two independent ageing methods.
12