New publication by GEMM Lab reveals sub-population health differences in gray whales 

Dr. KC Bierlich, Postdoctoral Scholar, OSU Department of Fisheries, Wildlife, & Conservation Sciences, Geospatial Ecology of Marine Megafauna (GEMM) Lab

In a previous blog, I discussed the importance of incorporating measurement uncertainty in drone-based photogrammetry, as drones with different sensors, focal length lenses, and altimeters will have varying levels of measurement accuracy. In my last blog, I discussed how to incorporate photogrammetric uncertainty when combining multiple measurements to estimate body condition of baleen whales. In this blog, I will highlight our recent publication in Frontiers in Marine Science (https://doi.org/10.3389/fmars.2022.867258) led by GEMM Lab’s Dr. Leigh TorresClara Bird, and myself that used these methods in a collaborative study using imagery from four different drones to compare gray whale body condition on their breeding and feeding grounds (Torres et al., 2022).

Most Eastern North Pacific (ENP) gray whales migrate to their summer foraging grounds in Alaska and the Arctic, where they target benthic amphipods as prey. A subgroup of gray whales (~230 individuals) called the Pacific Coast Feeding Group (PCFG), instead truncates their migration and forages along the coastal habitats between Northern California and British Columbia, Canada (Fig. 1). Evidence from a recent study lead by GEMM Lab’s Lisa Hildebrand (see this blog) found that the caloric content of prey in the PCFG range is of equal or higher value than the main amphipod prey in the Arctic/sub-Arctic regions (Hildebrand et al., 2021). This implies that greater prey density and/or lower energetic costs of foraging in the Arctic/sub-Arctic may explain the greater number of whales foraging in that region compared to the PCFG range. Both groups of gray whales spend the winter months on their breeding and calving grounds in Baja California, Mexico. 

Figure 1. The GEMM Lab field team following a Pacific Coast Feeding Group (PCFG) gray whale swimming in a kelp bed along the Oregon Coast during the summer foraging season. 

In January 2019 an Unusual Mortality Event (UME) was declared for gray whales due to the elevated numbers of stranded gray whales between Mexico and the Arctic regions of Alaska. Most of the stranded whales were emaciated, indicating that reduced nutrition and starvation may have been the causal factor of death. It is estimated that the population dropped from ~27,000 individuals in 2016 to ~21,000 in 2020 (Stewart & Weller, 2021).

During this UME period, between 2017-2019, the GEMM Lab was using drones to monitor the body condition of PCFG gray whales on their Oregon coastal feeding grounds (Fig. 1), while Christiansen and colleagues (2020) was using drones to monitor gray whales on their breeding grounds in San Ignacio Lagoon (SIL) in Baja California, Mexico. We teamed up with Christiansen and colleagues to compare the body condition of gray whales in these two different areas leading up to the UME. Comparing the body condition between these two populations could help inform which population was most effected by the UME.

The combined datasets consisted of four different drones used, thus different levels of photogrammetric uncertainty to consider. The GEMM Lab collected data using a DJI Phantom 3 Pro, DJI Phantom 4, and DJI Phantom 4 Pro, while Christiansen et al., (2020) used a DJI Inspire 1 Pro. By using the methodological approach described in my previous blog (here, also see Bierlich et al., 2021a for more details), we quantified photogrammetric uncertainty specific to each drone, allowing cross-comparison between these datasets. We also used Body Area Index (BAI), which is a standardized relative measure of body condition developed by the GEMM Lab (Burnett et al., 2018) that has low uncertainty with high precision, making it easier to detect smaller changes between individuals (see blog here, Bierlich et al., 2021b). 

While both PCFG and ENP gray whales visit San Ignacio Lagoon in the winter, we assume that the photogrammetry data collected in the lagoon is mostly of ENP whales based on their considerably higher population abundance. We also assume that gray whales incur low energetic cost during migration, as gray whale oxygen consumption rates and derived metabolic rates are much lower during migration than on foraging grounds (Sumich, 1983). 

Interestingly, we found that gray whale body condition on their wintering grounds in San Ignacio Lagoon deteriorated across the study years leading up to the UME (2017-2019), while the body condition of PCFG whales on their foraging grounds in Oregon concurrently increased. These contrasting trajectories in body condition between ENP and PCFG whales implies that dynamic oceanographic processes may be contributing to temporal variability of prey available in the Arctic/sub-Arctic and PCFG range. In other words, environmental conditions that control prey availability for gray whales are different in the two areas. For the ENP population, this declining nutritive gain may be associated with environmental changes in the Arctic/sub-Arctic region that impacted the predictability and availability of prey. For the PCFG population, the increase in body condition across years may reflect recovery of the NE Pacific Ocean from the marine heatwave event in 2014-2016 (referred to as “The Blob”) that resulted with a period of low prey availability. These findings also indicate that the ENP population was primarily impacted in the die-off from the UME. 

Surprisingly, the body condition of PCFG gray whales in Oregon was regularly and significantly lower than whales in San Ignacio Lagoon (Fig. 2). To further investigate this potential intrinsic difference in body condition between PCFG and ENP whales, we compared opportunistic photographs of gray whales feeding in the Northeastern Chukchi Sea (NCS) in the Arctic collected from airplane surveys. We found that the body condition of PCFG gray whales was significantly lower than whales in the NCS, further supporting our finding that PCFG whales overall have lower body condition than ENP whales that feed in the Arctic (Fig. 3). 

Figure 2. Boxplots showing the distribution of Body Area Index (BAI) values for gray whales imaged by drones in San Ignacio Lagoon (SIL), Mexico and Oregon, USA. The data is grouped by phenology group: End of summer feeding season (departure Oregon vs. arrival SIL) and End of wintering season (arrival Oregon vs. departure SIL). The group median (horizontal line), interquartile range (IQR, box), maximum and minimum 1.5*IQR (vertical lines), and outliers (dots) are depicted in the boxplots. The overlaid points represent the mean of the posterior predictive distribution for BAI of an individual and the bars represents the uncertainty (upper and lower bounds of the 95% HPD interval). Note how PCFG whales at then end of the feeding season (dark green) typically have lower body condition (as BAI) compared to ENP whales at the end of the feeding season when they arrive to SIL after migration (light brown).
Figure 3. Boxplots showing the distribution of Body Area Index (BAI) values of gray whales from opportunistic images collected from a plane in Northeaster Chukchi Sea (NCS) and from drones collected by the GEMM Lab in Oregon. The boxplots display the group median (horizontal line), interquartile range (IQR box), maximum and minimum 1.5*IQR (vertical lines), and outlies (dots). The overlaid points are the BAI values from each image. Note the significantly lower BAI of PCFG whales on Oregon feeding grounds compared to whales feeding in the Arctic region of the NCS.

This difference in body condition between PCFG and ENP gray whales raises some really interesting and prudent questions. Does the lower body condition of PCFG whales make them less resilient to changes in prey availability compared to ENP whales, and thus more vulnerable to climate change? If so, could this influence the reproductive capacity of PCFG whales? Or, are whales that recruit into the PCFG adapted to a smaller morphology, perhaps due to their specialized foraging tactics, which may be genetically inherited and enables them to survive with reduced energy stores?

These questions are on our minds here at the GEMM Lab as we prepare for our seventh consecutive field season using drones to collect data on PCFG gray whale body condition. As discussed in a previous blog by Dr. Alejandro Fernandez Ajo, we are combining our sightings history of individual whales, fecal hormone analyses, and photogrammetry-based body condition to better understand gray whales’ reproductive biology and help determine what the consequences are for these PCFG whales with lower body condition.

Did you enjoy this blog? Want to learn more about marine life, research, and conservation? Subscribe to our blog and get a weekly message when we post a new blog. Just add your name and email into the subscribe box below.

Loading

References

Bierlich, K. C., Hewitt, J., Bird, C. N., Schick, R. S., Friedlaender, A., Torres, L. G., … & Johnston, D. W. (2021). Comparing Uncertainty Associated With 1-, 2-, and 3D Aerial Photogrammetry-Based Body Condition Measurements of Baleen Whales. Frontiers in Marine Science, 1729.

Bierlich, K. C., Schick, R. S., Hewitt, J., Dale, J., Goldbogen, J. A., Friedlaender, A.S., et al. (2021b). Bayesian Approach for Predicting Photogrammetric Uncertainty in Morphometric Measurements Derived From Drones. Mar. Ecol. Prog. Ser. 673, 193–210. doi: 10.3354/meps13814

Burnett, J. D., Lemos, L., Barlow, D., Wing, M. G., Chandler, T., & Torres, L. G. (2018). Estimating morphometric attributes of baleen whales with photogrammetry from small UASs: A case study with blue and gray whales. Marine Mammal Science35(1), 108–139.

Christiansen, F., Rodrı́guez-González, F., Martı́nez-Aguilar, S., Urbán, J., Swartz, S., Warick, H., et al. (2021). Poor Body Condition Associated With an Unusual Mortality Event in Gray Whales. Mar. Ecol. Prog. Ser. 658, 237–252. doi:10.3354/meps13585

Hildebrand, L., Bernard, K. S., and Torres, L. G. (2021). Do Gray Whales Count Calories? Comparing Energetic Values of Gray Whale Prey Across Two Different Feeding Grounds in the Eastern North Pacific. Front. Mar. Sci. 8. doi: 10.3389/fmars.2021.683634

Stewart, J. D., and Weller, D. (2021). Abundance of Eastern North Pacific Gray Whales 2019/2020 (San Diego, CA: NOAA/NMFS)

Sumich, J. L. (1983). Swimming Velocities, Breathing Patterns, and Estimated Costs of Locomotion in Migrating Gray Whales, Eschrichtius Robustus. Can. J. Zoology. 61, 647–652. doi: 10.1139/z83-086

Torres, L.G., Bird, C., Rodrigues-Gonzáles, F., Christiansen F., Bejder, L., Lemos, L., Urbán Ramírez, J., Swartz, S., Willoughby, A., Hewitt., J., Bierlich, K.C. (2022). Range-wide comparison of gray whale body condition reveals contrasting sub-population health characteristics and vulnerability to environmental change. Frontiers in Marine Science. 9:867258. https://doi.org/10.3389/fmars.2022.867258

Dive into Oregon’s underwater forests

By Lisa Hildebrand, PhD student, OSU Department of Fisheries, Wildlife, & Conservation Sciences, Geospatial Ecology of Marine Megafauna Lab

When I was younger, I aspired to be a marine mammal biologist. I thought it was purely about knowing as much about marine mammal species as possible. However, over time and with experience in this field, I have realized that in order to understand a species, you need to have a holistic understanding of its prey, habitat, and environment. When I first applied to be advised by Leigh in the GEMM Lab, I had no idea how much of my time I would spend looking at tiny zooplankton under a microscope, thinking about the different benefits of different habitat types, or reading about oceanographic processes. But these things have been incredibly vital to my research to date and as a result, I now refer to myself as a marine ecologist. This holistic understanding that I am gaining will only grow throughout my PhD as I am broadly looking at the habitat use, site fidelity, and population dynamics of the Pacific Coast Feeding Group (PCFG) of gray whales for my thesis research. 

The PCFG display many foraging tactics and occupy several habitat types along the Oregon coast while they spend their summer feeding seasons here (Torres et al. 2018). Here, I will focus on one of these habitats: kelp. When you hear the word kelp, you probably conjure an image of long, thick stalks that reach from the ocean floor to the surface, with billowing fronds waving around (Figure 1a). However, this type is only one of three basic morphologies (Filbee-Dexter & Scheibling 2014) and it is called canopy kelp, which often forms extensive forests. The other two morphologies are stipitate and prostrate kelps. The former forms midwater stands (Figure 1b) while the latter forms low-lying kelp beds (Figure 1c). All three of these morphologies exist on the Oregon coast and create a mosaic of understory and canopy kelp patches that dot our coastline.

Figure 1. Examples of the three different kelp morphologies. a: bull kelp (Nereocystis luetkeana) is a type of canopy kelp and the dominant kelp on the Oregon coast (Source: Oregon Coast Aquarium); b: sea palm (Postelsia palmaeformis) is a type of stipitate kelp that forms mid-water stands (Source: Oregon Conservation Strategy); c: sea cabbage (Saccharina sessilis) is a type of prostrate kelp that is stipeless and forms low-lying kelp beds (Source: Central Coast Biodiversity).

One of the most magnificent things about kelp is that it is not just a species itself, but it provides critical habitat, refuge, and food resources to a myriad of other species due to its high rates of primary production (Dayton 1985). Kelp is often referred to as a foundation species due to all of these critical services it provides. In Oregon, many species of rockfish, which are important commercial and recreational fisheries, use kelp as habitat throughout their life cycle, including as nursery grounds. Lingcod, another widely fished species, forages amongst kelp. A large number of macroinvertebrates can be found in Oregon kelp forests, including anemones, limpets, snails, sea urchins, sea stars, and abalone, to name a fraction of them. 

Kelps grow best in cold, nutrient-rich waters (Tegner et al. 1996) and their growth and distribution patterns are highly naturally variable on both temporal and spatial scales (Krumhansl et al. 2016). However, warm water, low nutrient or light conditions, intensive grazing by herbivores, and severe storm activity can lead to the erosion and defoliation of kelp beds (Krumhansl et al. 2016). While these events can occur naturally in cyclical patterns, the frequency of several of these events has increased in recent years, as a result of climate change and anthropogenic impacts. For example, Dawn’s blog discussed increasing marine heatwaves that represent an influx of warm water for a prolonged period of time. In fact, kelps can be useful sentinels of change as they tend to be highly responsive to changes in environmental conditions (e.g., Rogers-Bennet & Catton 2019) and their nearshore, coastal location directly exposes them to human activities, such as pollution, harvesting, and fishing (Bennett et al. 2016).

Due to its foundational role, changes or impacts to kelp can reverberate throughout the ecosystem and negatively affect many other species. As mentioned previously, kelp is naturally highly variable, and like many other ecological processes, undergoes boom and bust cycles. For over four decades, dense, productive kelp forests have been shown to transition to sea urchin barrens, and back again, in natural cycles (Sala et al. 1998; Pinnegar et al. 2000; Steneck et al. 2002; Figure 2). These transitions are called phase shifts. In a healthy, balanced kelp forest, sea urchins typically passively feed on detrital plant matter, such as broken off pieces of kelp fronds that fall to the seafloor. A phase shift occurs when the grazing intensity of sea urchins increases, resulting in them actively feeding on kelp stalks and fronds to a point where the kelp in an area can become greatly reduced, creating an urchin barren. Sea urchin grazing intensity can change for a number of reasons, including reduction in sea urchin predators (e.g., sea otters, sunflower sea stars) or poor kelp recruitment events (e.g., due to warm water temperature). Regardless of the reason, the phases tend to transition back and forth over time. However, there is concern that sea urchin barrens may become an alternative stable state of the subtidal ecosystem from which kelp in an area cannot recover (Filbee-Dexter & Scheibling 2014). 

Figure 2. Screenshots from GoPro videos from 2016 (left) and 2018 (right) at the same kayak sampling station in Port Orford showing the difference between a dense kelp forest and what appears to be an urchin barren. (Source: GEMM Lab).

For example, in 2014, bull kelp canopy cover in northern California was reduced by >90% and has not shown signs of recovery since (Rogers-Bennet & Catton 2019; Figure 3). This massive decline was attributed to two major events: 1) the onset of sea star wasting disease (SSWD) in 2013 and 2) the “warm blob” of 2014-2016. SSWD affected over 20 sea star species along the coast from Mexico to Alaska, with the predatory sunflower sea star, which consumes purple sea urchins, most affected, including population declines of 80-100% along the coast (Harvell et al. 2019). Following this SSWD outbreak, the “warm blob”, which was an extreme marine heatwave in the Pacific Ocean, caused ocean temperatures to spike. These two events allowed purple sea urchin populations to grow unchecked by their predators, and created nutrient-poor and warm water conditions, which limited kelp growth and productivity. Intense grazing on bull kelp by growing urchin populations resulted in the >90% reduction in bull kelp canopy cover and has left behind widespread urchin barrens instead (Rogers-Bennet & Catton 2019). Consequently, there have been ecological and economic impacts on the ecosystem and communities in northern California. Without bull kelp, red abalone and red sea urchin populations starved, leading to a subsequent loss of the recreational red abalone (estimated value of $44 million/year) and commercial red urchin fisheries in northern California (Rogers-Bennet & Catton 2019).

Figure 3. Surface kelp canopy area pre- and post-impact from sites in Sonoma and Mendocino counties, northern California from aerial surveys (2008, 2014-2016). Figure and figure caption taken from Rogers-Bennett & Catton (2019).

As I mentioned earlier, while phase shifts between kelp forests and urchin barrens are common cycles, the intensity of the events described above in northern California are an example of sea urchin barrens potentially becoming a stable state of the subtidal ecosystem (Filbee-Dexter & Scheibling 2014). Given that marine heatwaves are only expected to increase in intensity and frequency in the future (Frölicher et al. 2018), the events documented in northern California may not be an isolated incidence. 

Considering that parts of the Oregon coast, particularly the southern portion, are very similar to northern California biogeographically, and that it was not exempt from the “warm blob”, similar changes in kelp forests may be occurring along our coast. There are many individuals and groups that are actively working on this issue to examine potential impacts to kelp and the species that depend on the services it provides. For more information, check out the Oregon Kelp Alliance

Figure 4. A gray whale surfaces in a large kelp bed during a foraging bout along the Oregon coast. (Source: GEMM Lab).

So, what does all of this information have to do with gray whales? Given their affinity for kelp habitats (Figure 4) and their zooplankton prey that aggregates there, changes to kelp ecosystems may affect gray whale health and ecology. This aspect of the complex kelp trophic web has not been examined to date; thus one of my PhD chapters focuses on the response of gray whales to changing kelp ecosystems along the southern Oregon coast. To do this, I am examining 6 years of data collected during the TOPAZ/JASPER project in Port Orford, to look at the relationships between kelp health, sea urchin density, zooplankton abundance, and gray whale foraging effort over space and time. Documenting impacts of changing kelp forests on gray whales is important to assist management efforts as healthy and abundant kelp seems critical in providing ample food opportunities for these iconic Pacific Northwest marine predators.

Did you enjoy this blog? Want to learn more about marine life, research, and conservation? Subscribe to our blog and get a weekly alert when we make a new post! Just add your name into the subscribe box on the left panel.

References

Bennett S, et al. The ‘Great Southern Reef’: Social, ecological and economic value of Australia’s neglected kelp forests. Marine and Freshwater Research 67:47-56.

Dayton PK (1985) Ecology of kelp communities. Annual Review of Ecology and Systematics 16:215-245.

Filbee-Dexter K, Scheibling RE (2014) Sea uechin barrens as alternative stable states of collapsed kelp ecosystems. Marine Ecology Progress Series 495:1-25.

Frölicher TL, Fischer EM, Gruber N (2018) Marine heatwaves under global warming. Nature 560:360-364.

Harvell CD, et al. (2019) Disease epidemic and a marine heat wave are associated with the continental-scale collapse of a pivotal predator (Pycnopodia helianthoides). Science Advances 5(1) doi:10.1126/sciadv.aau7042.

Krumhansl KA, et al. (2016) Global patterns of kelp forest change over the past half-century. Proceedings of the National Academy of Sciences of the United States of America 113(48):13785-13790.

Pinnegar JK, et al. (2000) Trophic cascades in benthic marine ecosystems: lessons for fisheries and protected-area management. Environmental Conservation 27:179-200.

Rogers-Bennett L, Catton CA (2019) Marine heat wave and multiple stressors tip bull kelp forest to sea urchin barrens. Scientific Reports 9:15050.

Sala E, Boudouresque CF, Harmelin-Vivien M (1998) Fishing, trophic cascades and the structure of algal assemblages; evaluation of an old but untested paradigm. Oikos 82:425-439.

Steneck RS, et al. (2002) Kelp forest ecosystems: biodiversity, stability, resilience and future. Environmental Conservation 29:436-459.

Tegner MJ, Dayton PK, Edwards PB, Riser KL (1996) Is there evidence for the long-term climatic change in southern California kelp forests? California Cooperative Oceanic Fisheries Investigations Report 37:111-126.

Torres LG, Nieukirk SL, Lemos L, Chandler TE (2018) Drone up! Quantifying whale behavior from a new perspective improves observational capacity. Frontiers in Marine Science doi:10.3389/fmars.2019.00319.

A pregnancy test for whales?! Why and how?

Dr. Alejandro A. Fernández Ajó, Postdoctoral Scholar, Marine Mammal Institute – OSU Department of Fisheries, Wildlife, & Conservation Sciences, Geospatial Ecology of Marine Megafauna (GEMM) Lab.

I often receive two reactions when asked what I am currently working on; one is “Wow! That is a very cool job, it must be amazing to work with such incredible animals!”, the other is “How do you do that and why is that important?”. So, today I decided to blog about some of the reasons why it is important to develop a pregnancy test for gray whales and how we are doing this.

In a previous blogpost, I described the many ways in which whales play critical roles in sustaining marine ecosystem. Briefly, whales can enhance marine productivity by vertically and horizontally mixing of ocean waters, promoting primary production, and mitigating climate change by sequestering carbon with their large biomass and long life-span (1-3). Even after they die, their carcasses can contribute to biodiversity creating new habitat on the seafloor (4). But, over several decades, the whaling industry drastically removed whales around the globe, with some species and populations depleted to near extinction (5). Consequently, these depleted whale populations now play a diminished role in ocean ecosystem processes and their recovery is currently challenged by an increasing number of modern anthropogenic impacts. Hence, working towards whale conservation is essential for keeping a healthy marine ecosystem.

Working and designing effective strategies for conservation biology often involves gaining knowledge regarding the reproductive parameters of individual animals in wild populations. This information is critical for understanding population trends and the underlaying mechanisms that affect animal welfare and their potential for recovery. However, getting such information from free-living whales can be challenging (see Hunt et al. 2013). While we know that whales typically have long life-spans, lengthy generation times, extended parental care, and high survival rates, detailed knowledge on the life history and general reproductive biology of free-ranging whales is limited for the majority of the whale populations. In fact, much of what we do know about whale reproduction is derived from whaling records. Only recently, conservation physiology approaches (see our previous post here) have contributed alternative and non-invasive methods for monitoring key physiological processes that can help monitor a whale’s reproductive biology and determine reproductive parameters such as sexual maturity and pregnancy (6-9).

In this clip you can see an example of a fecal sample collection from a gray whale off the Oregon coast. We can look at hormones in the fecal samples which are useful indicators for endocrine assessments of free-swimming whales. Fecal sample and footage filmed under NOAA/NMFS permit #16111.

Gray whales (Eschrichtius robustus) in the Eastern North Pacific (ENP) typically undertake annual migrations between their lower latitude breeding grounds in the coastal waters of the Baja California Peninsula, Mexico, and the foraging grounds located on the Bering and Chukchi Seas (10). However, among the ENP whales a distinct subgroup of about 230 whales shorten their migration to feed in the coastal waters of Northern California, Oregon, and southeastern Alaska (11). This group of whales is known as the gray whale Pacific Coast Feeding Group (PCFG).

Since 2016, the GEMM Lab has monitored individual gray whales within the PCFG off the Oregon coast (check the GRANITE project). Gray whales have a distinct mottled skin; and each individual whale presents a unique pigmentation pattern that allows for the individual identification of whales. We can identify who is who among the whales who visit the Oregon coast. In this way, we can keep a detailed record of re-sightings of known individuals (visit our new web site to know more about the lives of individual whales that visit the Oregon coast).  We have high individual re-sighting rates, so this unique opportunity helps us keep a long-term data series for individual whales to monitor their health, body condition, and reproductive status over time, and thus further develop and advance our non-invasive study methods.

We are combining behavioral and feeding ecology with drone photogrammetry and endocrinology of the same individual whales to help us understand the relationships between natural and anthropogenic drivers with biological parameters. In this way, following individual whales, we are developing sensitive biomarkers to monitor and infer about the population health, population trends, and identify stressors that impact their recovery and welfare. In particular, we are now working to develop a noninvasive approach to detect pregnancy in gray whales based on fecal hormone analyses.

In this picture you can see “Rose”, a gray whale calf, on top of her mother “Scarlett”. Scarlett is one of the most recognizable whales from the PCFG, due to a large scar on the right side of her back (not visible in this picture). She has been observed along the Pacific NW coast since 1996, so she is at least 26 years old today. We know 3 of her calves. Following individual whales like Scarlett is helping us to better understand the gray whale reproductive biology. Photo by Alejandro Fernandez Ajo taken under NOAA/NMFS permit #21678.

In marine mammals, the progesterone hormone is secreted in the ovaries during the estrous cycle and gestation, and is the predominant hormone responsible for sustaining pregnancy (12). As the hormones are cleared from the blood into the gut, they are metabolized and eventually excreted in feces; fecal samples represent a cumulative and integrated concentration of hormone metabolites (13-14), which are useful indicators for endocrine assessments of free-swimming whales. Several studies show that changes in hormone concentration correlate in meaningful ways with exposure to stressors (15-16) and changes in reproductive status (17-19). We are using our long data series of fecal hormones and individual life histories to advance our understanding on the gray whales’ reproductive biology. We are close to developing a technique that will allow us to detect pregnancy in whales based in fecal hormones analyses and photogrammetry. Stay tuned for results from this pregnancy test!

Did you enjoy this blog? Want to learn more about marine life, research, and conservation? Subscribe to our blog and get a weekly alert when we make a new post! Just add your name into the subscribe box on the left panel.

References:

1- Pershing AJ, Christensen LB, Record NR, Sherwood GD, Stetson PB (2010) The impact of whaling on the ocean carbon cycle: Why bigger was better. PLoS ONE 5(8): e12444.

2- Roman J and McCarthy JJ. 2010. The whale pump: marine mammals enhance primary productivity in a coastal basin. PLoS ONE. 5(10): e13255.

3- Morissette L, Kaschner K, and Gerber LR. 2010. “Whales eat fish”? Demystifying the myth in the Caribbean marine ecosystem. Fish Fish 11: 388–404.

4- Smith CR, Roman J, Nation JB. A metapopulation model for whale-fall specialists: The largest whales are essential to prevent species extinctions. J. Mar. Res. 77, 283–302 (2019).

5- Branch TA, Williams TM. Legacy of industrial whaling. Whales. Whal. Ocean Ecosyst. 2006, 262–278 (2006).

6- Kellar NM, Keliher J, Trego ML, Catelani KN, Hanns C, George JC, et al. Variation of bowhead whale progesterone concentrations across demographic groups and sample matrices. Endanger Species Res 2013; 22:61–72. https://doi.org/10.3354/esr00537.

7- Pallin L, Robbins J, Kellar N, Berube M, Friedlaender A. Validation of a blubber-based endocrine pregnancy test for humpback whales. Conserv Physiol 2018;6:1 11. https://doi.org/10.1093/conphys/coy031PMID:29942518.

8-Hunt KE, Robbins J, Buck CL, Bérubé M, Rolland RM (2019) Evaluation of fecal hormones for noninvasive research on reproduction and stress in humpback whales (Megaptera novaeangliae). Gen Comp Endocrinol 280: 24–34.

9-Melica, V., Atkinson, S., Calambokidis, J., Lang, A., Scordino, J., & Mueter, F. (2021). Application of endocrine biomarkers to update information on reproductive physiology in gray whale (Eschrichtius robustus). Plos one, 16(8), e0255368.

10-Swartz SL. Gray Whale. In: Wursig B, Thewissen JGM, Kovacs KM, editors. Encyclopedia of Marine Mammals (Third Edition). Elsevier;2018,p. 422–8.https://doi.org/10.1016/B978-0-12-804327-1.00140–0.

11-Calambokidis J, Darling JD, Deecke V, Gearin P, Gosho M, Megill W, et al. Abundance, range and movements of a feeding aggregation of gray whales (Eschrichtius robustus) from California to south-eastern Alaska in 1998. J Cetacean Res Manag 2002;4:267–76.

12- Bronson, F. H. (1989). Mammalian reproductive biology. University of Chicago Press.

13-Wasser SK, Hunt KE, Brown JL, Cooper K, Crockett CM, Bechert U, Millspaugh JJ, Larson S, Monfort SL (2000) A generalized fecal glucocorticoid assay for use in a diverse array of nondomestic mammalian and avian species. Gen Comp Endocrinol120:260–275.

14- Hunt, K.E., Rolland, R.M., Kraus, S.D., Wasser, S.K., 2006. Analysis of fecal glucocorticoids in the North Atlantic right whale (Eubalaena glacialis). Gen. Comp. Endocrinol. 148, 260–272. https://doi.org/10.1016/j.ygcen.2006.03.01215.

15- Lemos, L.S., Olsen, A., Smith, A., Burnett, J.D., Chandler, T.E., Larson, S., Hunt, K.E., Torres, L.G., 2021. Stressed and slim or relaxed and chubby? A simultaneous assessment of gray whale body condition and hormone variability. Mar. Mammal Sci. 1–11. https://doi.org/10.1111/mms.12877

16- Rolland, R., McLellan, W., Moore, M., Harms, C., Burgess, E., Hunt, K., 2017. Fecal glucocorticoids and anthropogenic injury and mortality in North Atlantic right whales Eubalaena glacialis. Endanger. Species Res. 34, 417–429. https://doi.org/10.3354/esr00866.

17-Rolland, R.M., Hunt, K.E., Kraus, S.D., Wasser, S.K., 2005. Assessing reproductive status of right whales (Eubalaena glacialis) using fecal hormone metabolites. Gen. Comp. Endocrinol. 142, 308–317. https://doi.org/10.1016/j.ygcen.2005.02.002

18- Valenzuela Molina M, Atkinson S, Mashburn K, Gendron D, Brownell RL. Fecal steroid hormones reveal reproductive state in female blue whales sampled in the Gulf of California, Mexico. Gen Comp Endocrinol 2018;261:127–35.https://doi.org/10.1016/j.ygcen.2018.02.015 PMID:29476760.

19- Hunt, K. E., Robbins, J., Buck, C. L., Bérubé, M., & Rolland, R. M. (2019). Evaluation of fecal hormones for noninvasive research on reproduction and stress in humpback whales (Megaptera novaeangliae). General and Comparative Endocrinology, 280, 24-34.

What drives individual specialization?

Clara Bird, PhD Student, OSU Department of Fisheries, Wildlife, and Conservation Sciences, Geospatial Ecology of Marine Megafauna Lab

When I wrote my first blog on individual specialization well over a year ago, I just skimmed the surface of the literature on this topic and only started to recognize the importance of studying individual specialization. The question, “is there individual specialization in the PCFG of gray whales?” is the focus of my first thesis chapter and the results will affect all my subsequent work. Therefore, the literature and concepts of individual specialization are a focus of my literature review and studies.

In my previous blog I focused on common characteristics of individuals that are similarly specialized as an underlying driver of individual specialization. While these characteristics (often attributable to age, sex, or physical traits) are important to consider, I’ve learned that the list of drivers of individual specialization is long and that many variables are dynamic. Of all the drivers I’ve learned about, competition is among the most common.

Competition is a major driver of individual specialization, and a common driver of competition is resource availability. When resource availability decreases, whether caused by increasing population density or changing environmental conditions, competition for that resource increases. As competition increases, individuals have a choice. They can choose to engage in competition, either by racing, fighting, or sharing [1], which can be costly, or they can diffuse the competition by focusing on a different resource.  This second approach would be considered niche partitioning in the prey dimension. Niche partitioning is a way for individuals to share ecological space by using different resources. Essentially, individuals can share habitat without having to engage in direct competition by pursuing different prey types [2]. 

This switch to different prey types can change the degree of individual specialization present in the population (Figure 1). But the direction of the change is not constant. If all individuals were pursuing the same prey type under low competition conditions but then switched to different alternate prey types under high competition, then individual specialization would increase (Figure 1a). This direction has been observed across a range of species including sharks [3], otters [4]–[7], dolphins [8], [9], stickleback fish [10], [11], largemouth bass [12], banded mongoose [13], fur seals [14], and baleen whales [15].

However, if individuals were pursuing different prey types under low competition conditions (maybe because of underlying differences such as age or sex) but then switched to the same alternate prey types under high competition, diet overlap would increase, and individual specialization would decrease (Figure 1b). Furthermore, an individual might not switch to an entirely new prey type but instead add prey items to its diet [16]. This diet expansion under competition would also decrease individual specialization. Fewer studies have reported this direction but it’s been found in the common bumblebee [17] and in several neotropical vertebrate species [18], [19].

Figure `1. Figure 3 from Araújo et al. 2011 [20]. Illustration of how ecological mechanisms may affect the degree of individual specialization. Arrows linking resources to individual consumers indicate resource consumption (relative thickness indicates proportional contribution). 
Horizontal arrows indicate the sign (positive or negative) of the effect on the degree of individual specialization. (a) Consumers share the same preferred resource (dark gray tangle) but have different alternative resources (white and light gray triangles). As the preferred resource becomes scarce, consumers switch to different alternatives, increasing the degree of individual specialization. (b) Alternatively, consumers have distinct preferred resources, so that as resources become scarce, individuals converge to the alternative resource (dark gray triangle), reducing diet variation.

Interestingly, its hypothesized that individual specialization driven by competition is one of the factors that facilitates the formation and existence of stable groups [21]. For example, a study on resident female dolphins in Sarasota Bay, FL, USA found that females with high spatial overlap used distinct foraging specializations [8](Figure 2). This study illustrates how partitioning prey enabled spatial and social coexistence. A study on banded mongooses reached a similar conclusion [13]. They found that specialization was highest in the biggest groups (with the most competition) and not explained by sex, age, or other inherent differences. They hypothesized that specialization increasing with competition reduced conflict and allowed the groups to remain stable. This study also highlighted the role of learning to determine an individual’s specialization.

Figure 2. A bottlenose dolphin.
Source: https://sarasotadolphin.org

Learning drives the distribution of knowledge throughout a population, which can lead to either specialization or generalization. ‘One-to-one’ learning, where one individual learns from one demonstrator, tends to promote individual specialization [21]. This form of transmission drives specialization because the individuals who learn the specialization tend to then carry on using, and eventually teaching, that specialization [6]. A common example of ‘one-to-one’ learning is vertical transmission from parent to offspring. It has been shown to transmit specializations in dolphins [22] and otters [6]. ‘One-to-one’ learning can occur outside of parent-offspring pairs; non-parent-offspring ‘one-to-one’ learning has been shown to drive specialization in banded mongooses [13](Figure 3).

However, other forms of social learning can promote more generalized foraging strategies. ‘Many-to-one’ or ‘one-to-many’ learning  can reduce the presence of specialization in species [13], [21] as can the presence of conformity in a group [23], [24].

Figure 3. A group of banded mongooses.
Source: http://socialisresearch.org/about-the-banded-mongoose-project/

The multiple drivers of specialization and their dynamic quality means that it is important to contextualize specialization. For example, a study on four species of neotropical frogs found varying degrees of specialization across multiple populations of each species [18]. The degree of specialization was dependent on a variety of drivers including predation and both intra- and inter-specific competition. Notably, the direction of the relationship between degree of specialization and each driver was species specific. This study highlights that one species may not always be more specialized than another, but that a populations’ specialization is context dependent.

Therefore, it is important to not only be aware of the degree of specialization present in a population, but to also understand its dynamics and drivers. These relationships can then be used to understand how, and why, a population may react to competition from other species, predators, and changes in resource availability [20].  A population’s specialization can also affect the specialization of other populations and community dynamics [25], therefore, it’s important to consider and study individual specialization on both the population and community level. I am excited to start using our incredible six-year dataset to start investigating these questions for PCFG gray whales, so stay tuned for results!

Did you enjoy this blog? Want to learn more about marine life, research, and conservation? Subscribe to our blog and get a weekly alert when we make a new post! Just add your name into the subscribe box on the left panel.  

References

[1]       M. Taborsky, M. A. Cant, and J. Komdeur, The Evolution of Social Behaviour. Cambridge: Cambridge University Press, 2021. doi: 10.1017/9780511894794.

[2]       E. R. Pianka, “Niche Overlap and Diffuse Competition,” vol. 71, no. 5, pp. 2141–2145, 1974.

[3]       P. Matich et al., “Ecological niche partitioning within a large predator guild in a nutrient-limited estuary,” Limnol. Oceanogr., vol. 62, no. 3, pp. 934–953, 2017, doi: https://doi.org/10.1002/lno.10477.

[4]       S. D. Newsome et al., “The interaction of intraspecific competition and habitat on individual diet specialization: a near range-wide examination of sea otters,” Oecologia, vol. 178, no. 1, pp. 45–59, May 2015, doi: 10.1007/s00442-015-3223-8.

[5]       M. T. Tinker, G. Bentall, and J. A. Estes, “Food limitation leads to behavioral diversification and dietary specialization in sea otters,” Proc. Natl. Acad. Sci., vol. 105, no. 2, pp. 560–565, Jan. 2008, doi: 10.1073/pnas.0709263105.

[6]       M. T. Tinker, M. Mangel, and J. A. Estes, “Learning to be different: acquired skills, social learning, frequency dependence, and environmental variation can cause behaviourally mediated foraging specializations,” Evol. Ecol. Res., vol. 11, pp. 841–869, 2009.

[7]       M. T. Tinker et al., “Structure and mechanism of diet specialisation: testing models of individual variation in resource use with sea otters,” Ecol. Lett., vol. 15, no. 5, pp. 475–483, 2012, doi: 10.1111/j.1461-0248.2012.01760.x.

[8]       S. Rossman et al., “Foraging habits in a generalist predator: Sex and age influence habitat selection and resource use among bottlenose dolphins (Tursiops truncatus),” Mar. Mammal Sci., vol. 31, no. 1, pp. 155–168, 2015, doi: https://doi.org/10.1111/mms.12143.

[9]       L. G. Torres, “A kaleidoscope of mammal , bird and fish : habitat use patterns of top predators and their prey in Florida Bay,” vol. 375, pp. 289–304, 2009, doi: 10.3354/meps07743.

[10]     M. S. Araújo et al., “Network Analysis Reveals Contrasting Effects of Intraspecific Competition on Individual Vs. Population Diets,” Ecology, vol. 89, no. 7, pp. 1981–1993, 2008, doi: 10.1890/07-0630.1.

[11]     R. Svanbäck and D. I. Bolnick, “Intraspecific competition drives increased resource use diversity within a natural population,” Proc. R. Soc. B Biol. Sci., vol. 274, no. 1611, pp. 839–844, Mar. 2007, doi: 10.1098/rspb.2006.0198.

[12]     D. E. Schindler, J. R. Hodgson, and J. F. Kitchell, “Density-dependent changes in individual foraging specialization of largemouth bass,” Oecologia, vol. 110, no. 4, pp. 592–600, May 1997, doi: 10.1007/s004420050200.

[13]     C. E. Sheppard et al., “Intragroup competition predicts individual foraging specialisation in a group-living mammal,” Ecol. Lett., vol. 21, no. 5, pp. 665–673, 2018, doi: 10.1111/ele.12933.

[14]     L. Kernaléguen, J. P. Y. Arnould, C. Guinet, and Y. Cherel, “Determinants of individual foraging specialization in large marine vertebrates, the Antarctic and subantarctic fur seals,” J. Anim. Ecol., vol. 84, no. 4, pp. 1081–1091, 2015, doi: 10.1111/1365-2656.12347.

[15]     E. M. Keen and K. M. Qualls, “Respiratory behaviors in sympatric rorqual whales: the influence of prey depth and implications for temporal access to prey,” J. Mammal., vol. 99, no. 1, pp. 27–40, Feb. 2018, doi: 10.1093/jmammal/gyx170.

[16]     R. H. MacArthur and E. R. Pianka, “On Optimal Use of a Patchy Environment,” Am. Nat., vol. 100, no. 916, pp. 603–609, 1966, doi: 10.1086/282454.

[17]     C. Fontaine, C. L. Collin, and I. Dajoz, “Generalist foraging of pollinators: diet expansion at high density,” J. Ecol., vol. 96, no. 5, pp. 1002–1010, 2008, doi: 10.1111/j.1365-2745.2008.01405.x.

[18]     R. Costa-Pereira, V. H. W. Rudolf, F. L. Souza, and M. S. Araújo, “Drivers of individual niche variation in coexisting species,” J. Anim. Ecol., vol. 87, no. 5, pp. 1452–1464, 2018, doi: 10.1111/1365-2656.12879.

[19]     M. M. Pires, P. R. Guimarães Jr, M. S. Araújo, A. A. Giaretta, J. C. L. Costa, and S. F. dos Reis, “The nested assembly of individual-resource networks,” J. Anim. Ecol., vol. 80, no. 4, pp. 896–903, 2011, doi: 10.1111/j.1365-2656.2011.01818.x.

[20]     M. S. Araújo, D. I. Bolnick, and C. A. Layman, “The ecological causes of individual specialisation,”Ecol. Lett., vol. 14, no. 9, pp. 948–958, 2011, doi: https://doi.org/10.1111/j.1461-0248.2011.01662.x.

[21]     C. E. Sheppard, R. Heaphy, M. A. Cant, and H. H. Marshall, “Individual foraging specialization in group-living species,” Anim. Behav., vol. 182, pp. 285–294, Dec. 2021, doi: 10.1016/j.anbehav.2021.10.011.

[22]     S. Wild, S. J. Allen, M. Krützen, S. L. King, L. Gerber, and W. J. E. Hoppitt, “Multi-network-based diffusion analysis reveals vertical cultural transmission of sponge tool use within dolphin matrilines,” Biol. Lett., vol. 15, no. 7, p. 20190227, Jul. 2019, doi: 10.1098/rsbl.2019.0227.

[23]     L. M. Aplin, D. R. Farine, J. Morand-Ferron, A. Cockburn, A. Thornton, and B. C. Sheldon, “Experimentally induced innovations lead to persistent culture via conformity in wild birds,” Nature, vol. 518, no. 7540, pp. 538–541, Feb. 2015, doi: 10.1038/nature13998.

[24]     E. Van de Waal, C. Borgeaud, and A. Whiten, “Potent Social Learning and Conformity Shape a Wild Primate’s Foraging Decisions,” Science, Apr. 2013, doi: 10.1126/science.1232769.

[25]     D. I. Bolnick et al., “Why intraspecific trait variation matters in community ecology,” Trends Ecol. Evol., vol. 26, no. 4, pp. 183–192, Apr. 2011, doi: 10.1016/j.tree.2011.01.009.

The many dimensions of a fat whale: Using drones to measure the body condition of baleen whales 

Dr. KC Bierlich, Postdoctoral Scholar, OSU Department of Fisheries, Wildlife, & Conservation Sciences, Geospatial Ecology of Marine Megafauna (GEMM) Lab

In my last blog, I discussed how to obtain morphological measurements from drone-based imagery of whales and the importance of calculating and considering uncertainty, as different drone platforms have varying levels of measurement uncertainty. But how does uncertainty scale and propagate when multiple measurements are combined, such as when measuring body condition of the whole animal? In this blog, I will discuss the different methods used for measuring body condition of baleen whales from drone-based imagery and how uncertainty differs between these metrics.

Body condition is defined as the energy stored in the body as a result of feeding and is assumed to indicate an animal’s overall health, as it reflects the balance between energy intake and investment toward growth, maintenance and reproduction (Peig and Green, 2009). Thus, body condition reflects the foraging success of an individual, as well as the potential for reproductive output and the quality of habitat. For example, female North American brown bears (Ursus arctos) in high quality habitats were in better body condition, produced larger litter sizes, and lived in greater population densities compared to females in lower quality habitats (Hilderbrand et al., 1999). As Dawn Barlow and Will Kennerley discussed in their recent blog, baleen whales are top predators and serve as ecosystem sentinels that shed light not only on the health of their population, but on the health of their ecosystem. As ocean climate conditions continue to change, monitoring the body condition of baleen whales is important to provide insight on how their population and ecosystem is responding. 

As discussed in a previous blog, drones serve as a valuable tool for obtaining morphological measurements of baleen whales to estimate their body condition. Images are imported into photogrammetry software, such as MorphoMetriX (Torres and Bierlich, 2020), to measure the total length of an individual and that is then divided into perpendicular width segments (i.e., in 5 or 10% increments) down the body (Fig. 1). These total length and width measurements are then used to estimate body condition in either 1-, 2-, or 3-dimensions: a single width (1D), a projected dorsal surface area (2D), or a body volume measure (3D). These 1D, 2D, and 3D measurements of body condition can then be standardized by total length to produce a relative measure of an individual’s body condition to compare among individuals and populations. 

Figure 1. An example of a Pacific Coast Feeding Group (PCFG) gray whale measured in MorphoMetriX (Torres & Bierlich, 2020).

While several different studies have used each of these dimensions to assess whale body condition, it is unclear how these measurements compare amongst each other. Importantly, it is also unclear how measurement uncertainty scales across these multiple dimensions and influences inference, which can lead to misinterpretation of data. For example, the surface area and volume of two geometrically similar bodies of different sizes are not related to their linear dimensions in the same ratio, but rather to the second and third power, respectively (i.e., x2 vs. x3).  Similarly, uncertainty should not be expected to scale linearly across 1D, 2D, and 3D body condition measurements. 

The second chapter of my dissertation, which was recently published in Frontiers in Marine Science and includes Clara Bird and Leigh Torres as co-authors, compared the uncertainty associated with 1D, 2D, and 3D drone-based body condition measurements in three baleen whale species with different ranges in body sizes: blue, humpback, and Antarctic minke whales (Figure 2) (Bierlich et al., 2021). We used the same Bayesian model discussed in my last blog, to incorporate uncertainty associated with each 1D, 2D, and 3D estimate of body condition. 

Figure 2. An example of total length and perpendicular width (in 5% increments of total length) measurements of an individual blue, humpback and Antarctic minke whale. Each image measured using MorphoMetriX (Torres and Bierlich, 2020). 

We found that uncertainty does not scale linearly across multi-dimensional measurements, with 2D and 3D uncertainty increasing by a factor of 1.45 and 1.76 compared to 1D, respectively. This result means that there is an added cost of increased uncertainty when utilizing a multidimensional body condition measurement. Our finding is important to help researchers decide which body condition measurement best suits their scientific question,  particularly when using a drone platform that is susceptible to greater error – as discussed in my previous blog. However, a 1D measurement only relies on a single width measurement, which may be excluding other regions of an individual’s body condition that is important for energy storage. In these situations, a 2D or 3D measure may be more appropriate.

We found that when comparing relative measures of body condition (standardized by total length of the individual), each standardized metric was highly correlated with one another. This finding suggests that 1D, 2D, and 3D metrics will draw similar relative predictions of body condition for individuals, allowing researchers to be confident they will draw similar conclusions relating to the body condition of individuals, regardless of which standardized metric they use. However, when comparing the precision of each of these metrics, the body area index (BAI) – a 2D standardized metric – displayed the highest level of precision. This result highlights how BAI can advantageously detect small changes in body condition, which is useful for comparing individuals or even tracking the same individual over time.

BAI was developed by the GEMM Lab (Burnett et al., 2018) and was designed to be similar to body mass index (BMI) in humans [BMI = mass (kg)/(height (m))2], where BAI uses the calculated surface area as a surrogate for body mass. In humans, a healthy BMI range is generally considered 18.5–24.9, below 18.5 is considered underweight, above 24.9 is considered overweight, and above 30 is considered obese (Flegal et al., 2012). Identifying a healthy range in BAI for baleen whales is challenged by a limited knowledge of what a “healthy” body condition range is for a whale. We found strong evidence that a healthy range of BAI is species-specific, as each species displayed a distinctive range in BAI: blue whales: 11–16; AMW: 17–24; humpback whales: 23–32; humpback whale calves: 23–28 (Fig. 3). These differences in BAI ranges likely reflect differences in the body shape of each species (Fig. 4). For example, humpbacks have the widest range of BAI compared to these other two species, which was also reflected in their larger variation in perpendicular widths (Figs. 2-4). Thus, it seems that BAI offers conditionally “scalefree” comparisons between species, yet it is unreasonable to set a single, all-whale BAI threshold to determine “healthy” versus “unhealthy” body condition.  Collecting a large sample of body condition measurements across many individuals and demographic units over space and time with information on vital rates (e.g., reproductive capacity) will help elucidate a healthy BAI range for each species.

Figure 3. Body area index (BAI) for each species. AMW = Antarctic minke whale.  Figure from Bierlich et al. (2021).
Figure 4. A) Absolute widths (m) and B) relative widths, standardized by total length (TL) to help elucidate the different body shapes of Antarctic minke whales (AMW; n = 40), blue whales (n = 32), humpback whales (n = 40), and humpback whale calves (n = 15). Note how the peak in body width occurs at a different percent body width between species, demonstrating the natural variation in body shape between baleen whales. Figure from Bierlich et al. (2021).

Over the past six years, the GEMM Lab has been collecting drone images of Pacific Coast Feeding Group (PCFG) gray whales off the coast of Oregon to measure their BAI (see GRANITE Project blog). Many of the individuals we encounter are seen across years and throughout the foraging season, providing an opportunity to evaluate how an individual’s BAI is influenced by environmental variation, stress levels, maturity, and reproduction. These data will in turn help determine what the healthy range in BAI for gray whales is. For example, linking BAI to pregnancy – whether a whale is currently pregnant or becomes pregnant the following season – will help determine what BAI is needed to support calf production. We are currently analyzing hundreds of body condition measurements from 2016 – 2021, so stay tuned for upcoming results!

References

Bierlich, K. C., Hewitt, J., Bird, C. N., Schick, R. S., Friedlaender, A., Torres, L. G., … & Johnston, D. W. (2021). Comparing Uncertainty Associated With 1-, 2-, and 3D Aerial Photogrammetry-Based Body Condition Measurements of Baleen Whales. Frontiers in Marine Science, 1729.

Burnett, J. D., Lemos, L., Barlow, D., Wing, M. G., Chandler, T., & Torres, L. G. (2018). Estimating morphometric attributes of baleen whales with photogrammetry from small UASs: A case study with blue and gray whales. Marine Mammal Science35(1), 108–139.

Flegal, K. M., Carroll, M. D., Kit, B. K., & Ogden, C. L. (2012). Prevalence of Obesity and Trends in the Distribution of Body Mass Index Among US Adults, 1999-2010. JAMA307(5), 491. https://doi.org/10.1001/jama.2012.39

Hilderbrand, G. V, Schwartz, C. C., Robbins, C. T., Jacoby, M. E., Hanley, T. A., Arthur, S. M., & Servheen, C. (1999). The importance of meat, particularly salmon, to body size, population productivity, and conservation of North American brown bears. Canadian Journal of Zoology77(1), 132–138.

Peig, J., & Green, A. J. (2009). New perspectives for estimating body condition from mass/length data: the scaled mass index as an alternative method. Oikos118(12), 1883–1891.

Torres, W., & Bierlich, K. C. (2020). MorphoMetriX: a photogrammetric measurement GUI for morphometric analysis of megafauna. Journal of Open Source Software5(45), 1825–1826.

Wavelet analysis to describe biological cycles and signals of non-stationarity

By Allison Dawn, GEMM Lab Master’s student, OSU Department of Fisheries, Wildlife, and Conservation Sciences, Geospatial Ecology of Marine Megafauna Lab 

During my second term of graduate school, I have been preparing to write my research proposal. The last two months have been an inspiring process of deep literature dives and brainstorming sessions with my mentors. As I discussed in my last blog, I am interested in questions related to pattern and scale (fine vs. mesoscale) in the context of the Pacific Coast Feeding Group (PCFG) of gray whales, their zooplankton prey, and local environmental variables.

My work currently involves exploring which scales of pattern are important in these trophic relationships and whether the dominant scale of a pattern changes over time or space. I have researched which analysis tools would be most appropriate to analyze ecological time series data, like the impressive long-term dataset the GEMM lab has collected in Port Orford as part of the TOPAZ  project, where we have monitored the abundance of whales and zooplankton, as well as environmental variables since 2016. 

A useful analytical tool that I have come across in my recent coursework and literature review is called wavelet analysis. Importantly, wavelet analysis can handle non-stationarity and edge detection in time series data. Non-stationarity is when a dataset’s mean and/or variance can change over time or space, and edge detection is the identification of the change location (in time or space). For example, it is not just the cycles or “ups and downs” of zooplankton abundance I am interested in, but when in time or where in space these cycles of “ups and downs” might change in relation to what their previous values, or distances between values, were. Simply stated, non-stationarity is when what once was normal is no longer normal. Wavelet analysis has been applied across a broad range of fields, such as environmental engineering (Salas et al. 2020), climate science (Slater et al. 2021), and bio-acoustics (Buchan et al. 2021). It can be applied to any time series dataset that might violate the traditional statistical assumption of stationarity. 

In a recent review of climate science methodology, Slater et al. (2021) outlined the possible behavior of time series data. Using theoretical plots, the authors show that data can a) have the same mean and variance over time, or b) have non-stationarity that can be broken into three major groups – trend, step change, or shifts in variance. Figure 1 further demonstrates the difference between stationary vs. non-stationary data in relation to a given variable of interest over time. 

Figure 1. Plots showing the possible magnitude of a given variable across a time series: a) Stationary behavior, b) Non-stationary trend, step-change, and a shift in variance. [Taken from Slater et. al (2021)].

Traditional correlation statistics assumes stationarity, but it has been shown that ecological time series are often non-stationary at certain scales (Cazelles & Hales, 2006). In fact, ecological data rarely meets the requirements of a controlled experiment that traditional statistics require. This non-stationarity of ecological data means that while widely-used methods like generalized linear models and analyses of variances (ANOVAs) can be helpful to assess correlation, they are not always sufficient on their own to describe the complex natural phenomena ecologists seek to explain. Non-stationarity occurs frequently in ecological time series, so it is appropriate to consider analysis tools that will allow us to detect edges to further investigate the cause.

Wavelet analysis can also be conducted across a time series of multiple response variables to assess if these variables share high common power (correlation). When data is combined in this way it is called a cross-wavelet analysis. An interesting paper used cross-wavelet analysis to assess the seasonal response of zooplankton life history in relation to climate warming (Winder et. al 2009). Results from their cross-wavelet analysis showed that warming temperatures over the past two decades increased the voltinism (number of broods per year) of copepods. The authors show that where once annual recruitment followed a fairly stationary pattern, climate warming has contributed to a much more stochastic pattern of zooplankton abundance. From these results, the authors contribute to the hypothesis that climate change has had a temporal impact on zooplankton population dynamics, and recruitment has increasingly drifted out of phase from the original annual cycles. 

Figure 2. Cross-wavelet spectrum for immature and adult Leptodiaptomus ashlandi for 1965 through either 2000 or 2005. Plots show a) immatures and temperature, b) adults and temperature, c) immatures and phytoplankton, and d) adults and phytoplankton. Arrows indicate phase between combined time series. 0 degrees is in-phase and 180 degrees is anti-phase. Black contour lines show “cone of influence” or the 95% significance level, every value within the cone is considered significant. Left axis shows the temporal period, and the color legend shows wavelet frequency power, with low frequencies in blue and high frequencies in red. Plots show strong covariation of high common power at the 12-month period until the 1980s. This pattern is especially evident in plot c) and d). [Taken from (Winder et. al 2009)].

While wavelet and cross-wavelet analyses should not be the only tool used to explore data, due to its limitations with significance testing, it is still worth implementing to gain a better understanding of how time series variables relate to each other over multiple spatial and/or temporal scales. It is often helpful to combine multiple methods of analysis to get a larger sense of patterns in the data, especially in spatio-temporal research.

When conducting research within the context of climate change, where the concentration of CO2 in ppm in the atmosphere is a non-stationary time series itself (Figure 3), it is important to consider how our datasets might be impacted by climate change and wavelet analysis can help identify the scales of change. 

Figure 3. Plot showing the historic fluctuations of CO^2 and the recent deviation from normal levels. Source: https://globalclimate.ucr.edu/resources.html

When considering our ecological time series of data in Port Orford, we want to evaluate how changing ocean conditions may be related to data trends. For example, has the annual mean or variance of zooplankton abundance changed over time, and where has that change occurred in time or space? These changes might have occurred at different scales and might be invisible at other scales. I am eager to see if wavelet analysis can detect these sorts of changes in the abundance of zooplankton across our time series of data, particularly during the seasons of intense heat waves or upwelling. 

Did you enjoy this blog? Want to learn more about marine life, research and conservation? Subscribe to our blog and get a weekly email when we make a new post! Just add your name into the subscribe box on the left panel.

References

Buchan, S. J., Pérez-Santos, I., Narváez, D., Castro, L., Stafford, K. M., Baumgartner, M. F., … & Neira, S. (2021). Intraseasonal variation in southeast Pacific blue whale acoustic presence, zooplankton backscatter, and oceanographic variables on a feeding ground in Northern Chilean Patagonia. Progress in Oceanography, 199, 102709.

Cazelles, B., & Hales, S. (2006). Infectious diseases, climate influences, and nonstationarity. PLoS Medicine, 3(8), e328.

Salas, J. D., Anderson, M. L., Papalexiou, S. M., & Frances, F. (2020). PMP and climate variability and change: a review. Journal of Hydrologic Engineering, 25(12), 03120002.

Slater, L. J., Anderson, B., Buechel, M., Dadson, S., Han, S., Harrigan, S., … & Wilby, R. L. (2021). Nonstationary weather and water extremes: a review of methods for their detection, attribution, and management. Hydrology and Earth System Sciences, 25(7), 3897-3935.

Winder, M., Schindler, D. E., Essington, T. E., & Litt, A. H. (2009). Disrupted seasonal clockwork in the population dynamics of a freshwater copepod by climate warming. Limnology and Oceanography, 54(6part2), 2493-2505.

Introducing IndividuWhale!

By Lisa Hildebrand, PhD student, OSU Department of Fisheries, Wildlife, & Conservation Sciences, Geospatial Ecology of Marine Megafauna Lab

If you are an avid reader of our blog, you probably know quite a bit about gray whales, specifically the Pacific Coast Feeding Group (PCFG) of gray whales. Of the just over 50 GEMM Lab blogs written in 2021, 43% of them were about PCFG gray whales (or at least mentioned gray whales in some way). I guess this statistic is not too surprising when you consider that six of the 10 GEMM Lab members conduct gray whale-related research. You might think that we would have reached our annual limit of online gray whale content with that many blogs featuring these gentle giants, but you would in fact be wrong…

At the end of 2021, we launched a brand new website all about gray whales called IndividuWhale! It features stories of some of the Oregon coast’s most iconic gray whales, as well as information about how we study them, stressors they experience in our waters, and even a game to test your gray whale identification skills. IndividuWhale is a true labor of love that took over a year to create and that we are extremely proud to share with you today. Before I tell you more about the website, I want to take a moment to give a huge shout out to Erik Urdahl who was instrumental in getting this website off the ground and making it as interactive and beautiful as it is – hurrah Erik!

Equal‘s right side with visible boat propeller scars. Source: GEMM Lab.

Like us humans, gray whales have individual personalities and stories. They experience life-altering events, go through periods of stress, must provide for their offspring, and can behave differently to one another. Since Leigh & co. have been conducting in-depth research about PCFG gray whales in Oregon waters since 2016, we have been able to document several fascinating stories and events that these individuals have experienced. Take Equal, for example, a male whale that is at least 7 years old. The GEMM Lab observed Equal on consecutive days in June 2018, where on the first day he looked healthy and normal, but on the second day had fresh boat propeller-like scars on his back. Not only did we document these scars in photographs, but we were also able to collect a fecal sample from Equal the day we observed him with these scars. After analyzing his fecal sample for stress hormones, we discovered that Equal had very high stress levels compared to previous samples collected – unsurprising seeing as he had been hit by a boat! While this event was certainly sad for Equal (although don’t worry – we have seen him many times again in the years after this event looking healthy & normal once again), it was a very fortuitous occurrence for us since we were able to “validate” our stress hormone data relative to the value from Equal when he was clearly stressed out. Find out more about Equal as well as seven other gray whales here!

You might be wondering, how we knew that the whale with the boat propeller scar was Equal and how we recognize him again years after the incident. Gray whales have unique pigmentation patterns on their bodies and flukes that allow us to re-identify individuals between years and distinguish them from one another. Additionally, scars, such as those that Equal now carries on his back, can also be useful in telling whales apart. Therefore, we take photographs of every whale we see to match markings and identify whales. This process is called photo ID. Some individuals can have very distinctive markings, such as Roller Skate who has two big white dots on her right side, while others can look more inconspicuous, like Clouds. Therefore, conducting photo ID requires a lot of attention to detail and perseverance. To learn more about the different features we use to identify individuals, check out the “Studying Whales With Photographs” page. Do you think you have what it takes to tell individuals apart? Then try your luck at our photo ID game after!

Test your photo ID skills in our whale match game!

Unfortunately, these whales do not live in a pristine environment, as is evidenced by Equal’s story. Their main objective during the summer when in Oregon waters is to gain weight (energy stores) by consuming large amounts of prey, which is made more difficult by a number of stressors, including potential fishery entanglements, ocean noise, vessel traffic, and habitat changes. We describe these four stressors on the IndividuWhale website since we are trying to assess the impacts of them on gray whales through our research, however they are certainly not the only stressors that these whales experience. Little is known about the level at which these stressors may have a negative impact on the whales, and how whales react when they experience some of these stressors in concert. The answers to these questions are difficult to tease apart but the GRANITE project is aiming to do so through a framework called Population Consequences of Multiple Stressors (read more about it here). This approach requires a lot of data on a lot of individuals in a population and as you can see from the IndividuWhale website, we are slowly starting to get there! 2022 will certainly bring many more gray whale-themed blogs to this website, and you can share in our journey of learning about the lives of PCFG gray whales by exploring the IndividuWhale website (https://www.individuwhale.com).

Drones with lasers: almost as cool as “sharks with laser beams attached to their heads”

Dr. KC Bierlich, Postdoctoral Scholar, OSU Department of Fisheries, Wildlife, & Conservation Sciences, Geospatial Ecology of Marine Megafauna (GEMM) Lab

The recent advancement in drones (or unoccupied aircraft systems, UAS) has greatly enhanced opportunities for scientists across a broad range of disciplines to collect high-resolution aerial imagery. Wildlife researchers in particular have utilized this technology to study large elusive animals, such as whales, to observe their behavior (see Clara Bird’s blog) and obtain morphological measurements via photogrammetry (see previous blog for a brief history on photogrammetry and drones). However, obtaining useful measurement data is not as easy as flying the drone and pressing record. For this blog, I will provide a brief overview on the basics of using photogrammetry to extract morphological measurements from images collected with drones, as well as the associated uncertainty from using different drone platforms. 

During my PhD at Duke University, I co-developed an open-source photogrammetry software called MorphoMetriX to measure whales in images I collected using drones (Fig. 1) (Torres and Bierlich, 2020) (see this blog for some fieldwork memoirs!). The software is designed to be flexible, simple to use, and customizable without knowledge of scripting languages. Using MorphoMetriX, measurements are made in pixels and then multiplied by the ground sampling distance (GSD) to convert to standard units (e.g., meters) (Fig. 2A). GSD represents the distance on the ground each pixel represents (i.e., the linear size of the pixel) and therefore sets the scale of the image (i.e., cm per pixel). Figure 2A describes how GSD is dependent on the camera sensor, focal length lens, and altitude. Thus, drones equipped with different cameras and focal length lenses will have inherent differences in GSD as altitude increases (Fig. 2B). A larger GSD increases the length each pixel represents in a photo and results in a lower resolution image, potentially obscuring important features in the photo and introducing measurement error.

Figure 1. An example of a Pacific Coast Feeding Group gray whale measured in MorphoMetriX (Torres & Bierlich, 2020).
Figure 2: Overview of photogrammetry methods and calculating ground sampling distance (GSD). A) Photogrammetry methods for how each image is scaled to convert measurements in pixels to standard units (e.g., meters). Altitude is the distance between the camera lens and whale (usually at the surface of the water). Figure from Torres and Bierlich (2020). B) The exact (not accounting for distortion or altitude error) ground sampling distance (GSD) for two drone platforms commonly used to obtain morphological measurements of cetaceans. The difference in GSD between the P4Pro and Inspire 1 is due to the difference in sensor width and focal lengths of the cameras used. Figure from Bierlich et al. (2021).

Obtaining accurate altitude information is a key component in obtaining accurate measurements. All drones are equipped with a barometer, which measures altitude from changes in pressure. In general, barometers usually yield low accuracy in the altitude recorded, particularly for low-cost sensors commonly found on small, off-the-shelf drones (Wei et al., 2016). Dawson et al. (2017) added a laser altimeter (i.e., LightWare SF11/C, https://www.mouser.com//datasheet//2//321//28054-SF11-Laser-Altimeter-Manual-Rev8-1371857.pdf) to a drone, which yields higher accuracy in the altitude recorded. Since then, several studies have adopted use of a laser altimeter to study different species of baleen whales (i.e., Gough et al., 2019; Christiansen et al., 2018).

The first chapter of my dissertation, which was published last year in Marine Ecology Progress Series, compared the accuracy of several drones equipped with different camera sensors, focal length lenses, and a barometer vs. laser altimeter (Bierlich et al., 2021). We flew each drone over a known sized object floating at the surface and collected images at various altitudes (between 10 – 120 m). We used the known size of the floating object to determine the percent error of each measurement at each altitude. We found that 1) there is a lot of variation in measurement error across the different drones when using a barometer to measure altitude and 2) using a laser altimeter dramatically reduces measurement error for each drone (Fig. 3).

Figure 3. The % error for measurements from different drones. Black dashed line represents 0% error (true length = 1.48 m). The gray dashed lines represent under- and over-estimation of the true length by 5% (1.41 and 1.55 m, respectively).

These findings are important because if a study is analyzing measurements that are from more than one drone, the uncertainty associated with those measurements must be taken into account to know if measurements are reliable and comparable. For instance, let’s say we are comparing the body length of two different populations and found that population A is significantly longer than population B. From looking at Figure 3, that significant difference in length between population A and B could be unreliable as the difference may be due to the bias introduced by the type of drone, camera sensor, focal length lens, and whether a barometer or laser altimeter was used for recording altitude. In other words, without incorporating uncertainty associated with each measurement, how do you trust your measurement? 

Hence, the National Institute of Standards and Technology (NIST) states that a measurement is complete only when accompanied by a quantitative statement of its uncertainty (Taylor & Kuyatt, 1994). In our Bierlich et al. (2021) study, we develop a Bayesian statistical model where we use the measurements of the known-sized object floating at the surface (what was used for Fig. 3) as training data to predict the lengths of unknown-sized whales. This Bayesian approach views data and the underlying parameters that generated the data (such as the mean or standard deviation) as random, and thus can be described by a statistical distribution. Using Bayes’ Theorem, a model of the observed data (called the likelihood function), is combined with prior knowledge pertaining to the underlying parameters (called the prior probability distribution) to form the posterior probability distribution, which serves as updated knowledge about the underlying parameter. For example, if someone told me they saw a 75 ft blue whale, I would not be phased. But if someone told me they saw a 150 ft blue whale, I would be skeptical – I’m using prior knowledge to determine the probability of this statement being true. 

The posterior probability distribution produced by this Bayesian approach can also serve as new prior information for subsequent analyses. Following this framework, we used the known-sized objects to first estimate the posterior probability distribution for error for each drone. We then used that posterior probability distribution for error specific to each drone platform as prior information to form a posterior predictive distribution for length of unknown-sized whales. The length of an individual whale can then be described by the mean of this second posterior predictive distribution, and its uncertainty defined as the variance or an interval around the mean (Fig. 4). 

Figure 4. An example of a posterior predictive distribution for total length of an individual blue whale. The black bars represent the uncertainty around the mean value (the black dot) – the longer black bars represent the 95% highest posterior density (HPD) interval, and the shorter black bars represent the 65% HPD interval. 

For over half a decade, the GEMM Lab has been collecting drone images of Pacific Coast Feeding Group (PCFG) of gray whales off the coast of Oregon to measure their morphology and body condition (see GRANITE Project Blog). We have been using several different types of drones equipped with different cameras, focal length lenses, barometers, and laser altimeters. These measurements from different drones will inherently have different levels of error associated with them, so adapting these methods for incorporating uncertainty will be key to ensure our measurements are comparable and analyses are robust. To do this, we fly over a known-sized board (1 m) at the start of each flight to use as training data to generate a posterior predictive distribution for length of the an unknown-sized PCFG gray whale that we fly over (Fig. 5). Likewise, we are working closely with several other collaborators who are also using different drones. Incorporating measurement uncertainty from drones used across research labs and in different environments will help ensure robust analyses and provide great opportunity for some interesting comparisons – such as differences in gray whale body condition on their feeding grounds in Oregon vs. their breeding grounds in Baja, Mexico, and morphological comparisons with other baleen whale species, such as blue and humpback whales. We are currently wrapping up measurement from thousands of boards (Fig. 5) and whales (Fig. 1) from 2016 – 2021, so stay tuned for the results!

Figure 5. An example of a known-sized object (1 m long board) used as training data to assess measurement uncertainty. 

References

Bierlich, K.C., Schick, R.S., Hewitt, J., Dale, J., Goldbogen, J.A., Friedlaender, A.S., Johnston D.J. (2021). A Bayesian approach for predicting photogrammetric uncertainty in morphometric measurements derived from UAS. Marine Ecology Progress Series. DOI: https://doi.org/10.3354/meps13814

Christiansen F, Vivier F, Charlton C, Ward R, Amerson A, Burnell S, Bejder L (2018) Maternal body size and condition determine calf growth rates in southern right whales. Mar Ecol Prog Ser 592: 267−281

Dawson SM, Bowman MH, Leunissen E, Sirguey P (2017) Inexpensive aerial photogrammetry for studies of whales and large marine animals. Front Mar Sci 4: 366

Gough, W.T., Segre, P.S., Bierlich, K.C., Cade, D.E., Potvin, J., Fish, F. E., Dale, J., di Clemente, J., Friedlaender, A.S., Johnston, D.W., Kahane-Rapport, S.R., Kennedy, J., Long, J.H., Oudejans, M., Penry, G., Savoca, M.S., Simon, M., Videsen, S.K.A., Visser, F., Wiley, D.N., Goldbogen, J.A. (2019). Scaling of swimming performance in baleen whales. Journal of Experimental Biology222(20).https://doi.org/10.1242/jeb.204172  

Taylor, B. N., and Kuyatt, C. E. (1994). Guidelines for Evaluating and Expressing the Uncertainty of NIST Measurement Results. Washington, DC: National Institute of Standards and Technology. 1–25.

Torres, W.I., & Bierlich, K.C. (2020). MorphoMetriX: a photogrammetric measurement GUI for morphometric analysis of megafauna. Journal of Open Source Software5(45), 1825. https://doi.org/10.21105/joss.01825  

Wei S, Dan G, Chen H (2016) Altitude data fusion utilizing differential measurement and complementary filter. IET Sci Meas Technol (Singap) 10: 874−879

The costs and benefits of automated behavior classification

Clara Bird, PhD Student, OSU Department of Fisheries, Wildlife, and Conservation Sciences, Geospatial Ecology of Marine Megafauna Lab

“Why don’t you just automate it?” This is a question I am frequently asked when I tell someone about my work. My thesis involves watching many hours of drone footage of gray whales and meticulously coding behaviors, and there are plenty of days when I have asked myself that very same question. Streamlining my process is certainly appealing and given how wide-spread and effective machine learning methods have become, it is a tempting option to pursue. That said, machine learning is only appropriate for certain research questions and scales, and it’s important to consider these before investing in using a new tool.

The application of machine learning methods to behavioral ecology is called computational ethology (Anderson & Perona, 2014). To identify behaviors from videos, the model tracks individuals across video frames and identifies patterns of movement that form a behavior. This concept is similar to the way we identify a whale as traveling if it’s moving in a straight line and as foraging if it’s swimming in circles within a small area (Mayo & Marx, 1990, check out this blog to learn more). The level of behavioral detail that the model is able to track  depends on the chosen method (Figure 1, Pereira et al., 2020). These methods range from tracking each animal as a simple single point (called a centroid) to tracking the animal’s body positioning in 3D (this method is called pose estimation), which range from providing less detailed to more detailed behavior definitions. For example, tracking an individual as a centroid could be used to classify traveling and foraging behaviors, while pose estimation could identify specific foraging tactics. 

Figure 1. Figure from Pereira et al. (2020) illustrating the different methods of animal behavior tracking that are possible using machine learning.

Pose estimation involves training the machine learning algorithm to track individual anatomical features of an individual (e.g., the head, legs, and tail of a rat), meaning that it can define behaviors in great detail. A behavior state could be defined as a combination of the angle between the tail and the head, and the stride length. 

For example, Mearns et al. (2020) used pose estimation to study how zebrafish larvae in a lab captured their prey. They tracked the tail movements of individual larvae when presented with prey and classified these movements into separate behaviors that allowed them to associate specific behaviors with prey capture (Figure 2). The authors found that these behaviors occurred in a specific sequence, that the behaviors kept the prey within the larvae’s line of sight, and that the sequence was triggered by visual cues.  In fact, when they removed the visual cue of the prey, the larvae terminated the behavior sequence, meaning that the larvae are continually choosing to do each behavior in the sequence, rather than the sequence being one long behavior event that is triggered only by the initial visual cue. This study is a good example of the applicability of machine learning models for questions aimed at kinematics and fine-scale movements. Pose estimation has also been used to study the role of facial expression and body language in rat social communication (Ebbesen & Froemke, 2021). 

Figure 2. Excerpt from figure 1 of Mearns et al. (2020) illustrating (A) the camera set up for their experiment, (B) how the model tracked the eye angles and tail of the larvae fish, (C) the kinematics extracted from the footage. In panel (C) the top plot shows how the eyes converged on the same object (the prey) during prey capture event, the middle plot shows when the tail was curved to the left or the right, and the bottom plot shows the angle of the tail tip relative to the body.

While previous machine learning methods to track animal movements required individuals to be physically marked, the current methods can perform markerless tracking (Pereira et al., 2020). This improvement has broadened the kinds of studies that are possible. For example, Bozek et al., (2021) developed a model that tracked individuals throughout an entire honeybee colony and showed that certain individual behaviors were spatially distributed within the colony (Figure 3). Machine learning enabled the researchers to track over 1000 individual bees over several months, a task that would be infeasible for someone to do by hand. 

Figure 3. Excerpt from figure 1 of Bozek et al., (2021) showing how individual bees and their trajectories were tracked.

These studies highlight that the potential benefits of using machine learning when studying fine scale behaviors (like kinematics) or when tracking large groups of individuals. Furthermore, once it’s trained, the model can process large quantities of data in a standardized way to free up time for the scientists to focus on other tasks.

While machine learning is an exciting and enticing tool, automating behavior detection via machine learning could be its own PhD dissertation. Like most things in life, there are costs and benefits to using this technique. It is a technically difficult tool, and while applications exist to make it more accessible, knowledge of the computer science behind it is necessary to apply it effectively and correctly. Secondly, it can be tedious and time consuming to create a training dataset for the model to “learn” what each behavior looks like, as this step involves manually labeling examples for the model to use. 

As I’ve mentioned in a previous blog, I came quite close to trying to study the kinematics of gray whale foraging behaviors but ultimately decided that counting fluke beats wasn’t necessary to answer my behavioral research questions. It was important to consider the scale of my questions (as described in Allison’s blog) and I think that diving into more fine-scale kinematics questions could be a fascinating follow-up to the questions I’m asking in my PhD. 

For instance, it would be interesting to quantify how gray whales use their flukes for different behavior tactics. Do gray whales in better body condition beat their flukes more frequently while headstanding? Does the size of the fluke affect how efficiently they can perform certain tactics? While these analyses would help quantify the energetic costs of different behaviors in better detail, they aren’t necessary for my broad scale questions. Consequently, taking the time to develop and train a pose estimation machine learning model is not the best use of my time.

That being said, I am interested in applying machine learning methods to a specific subset of my dataset. In social behavior, it is not only useful to quantify the behaviors exhibited by each individual but also the distance between them. For example, the distance between a mom and her calf can be indicative of the calves’ dependence on its mom (Nielsen et al., 2019). However, continuously measuring the distance between two individuals throughout a video is tedious and time intensive, so training a machine learning model could be an effective use of time. I plan to work with an intern this summer to develop a machine learning model to track the distance between pairs of gray whales in our drone footage and then relate this distance data with the manually coded behaviors to examine patterns in social behavior (Figure 4).  Stay tuned to learn more about our progress!

Figure 4. A mom and calf pair surfacing together. Image collected under NOAA/NMFS permit #21678

Did you enjoy this blog? Want to learn more about marine life, research, and conservation? Subscribe to our blog and get a weekly alert when we make a new post! Just add your name into the subscribe box on the left panel.  

References

Anderson, D. J., & Perona, P. (2014). Toward a Science of Computational Ethology. Neuron84(1), 18–31. https://doi.org/10.1016/j.neuron.2014.09.005

Bozek, K., Hebert, L., Portugal, Y., Mikheyev, A. S., & Stephens, G. J. (2021). Markerless tracking of an entire honey bee colony. Nature Communications12(1), 1733. https://doi.org/10.1038/s41467-021-21769-1

Ebbesen, C. L., & Froemke, R. C. (2021). Body language signals for rodent social communication. Current Opinion in Neurobiology68, 91–106. https://doi.org/10.1016/j.conb.2021.01.008

Mayo, C. A., & Marx, M. K. (1990). Surface foraging behaviour of the North Atlantic right whale, Eubalaena glacialis , and associated zooplankton characteristics. Canadian Journal of Zoology68(10), 2214–2220. https://doi.org/10.1139/z90-308

Mearns, D. S., Donovan, J. C., Fernandes, A. M., Semmelhack, J. L., & Baier, H. (2020). Deconstructing Hunting Behavior Reveals a Tightly Coupled Stimulus-Response Loop. Current Biology30(1), 54-69.e9. https://doi.org/10.1016/j.cub.2019.11.022

Nielsen, M., Sprogis, K., Bejder, L., Madsen, P., & Christiansen, F. (2019). Behavioural development in southern right whale calves. Marine Ecology Progress Series629, 219–234. https://doi.org/10.3354/meps13125

Pereira, T. D., Shaevitz, J. W., & Murthy, M. (2020). Quantifying behavior to understand the brain. Nature Neuroscience23(12), 1537–1549. https://doi.org/10.1038/s41593-020-00734-z

Of snakes and whales: How food availability and body condition affect reproduction

Clara Bird, PhD Student, OSU Department of Fisheries, Wildlife, and Conservation Sciences, Geospatial Ecology of Marine Megafauna Lab

Over six field seasons the GEMM lab team has conducted nearly 500 drone flights over gray whales, equaling over 100 hours of footage. These hours of footage are the central dataset for my PhD dissertation, so it’s up to me to process them all. This process can be challenging, tedious, and daunting, but it is also quite fun and a privilege to be the one person who gets to watch all the footage. It’s fascinating to get to know the whales and their behaviors and pick up on patterns. It motivates me to get through this video processing step and start doing the data analysis. Recently, it’s been especially fun to notice patterns that I’ve seen mentioned in the literature. One example is adult social behavior. 

There are two categories of social behavior that I’m interested in studying: maternal behavior, defined as interactions between a mom and its calf, and general social behaviors, defined as social interactions between non-mom/calf pairs. In this blog I’ll focus on general social behaviors, but if you’re interested in maternal behavior check out this blog. General social behavior, which I’ll refer to as social behavior moving forward, includes tactile interactions and promiscuous behaviors (Torres et al. 2018; Clip 1). While gray whales in the PCFG range are primarily foraging, researchers have observed increases in social behavior towards the end of the foraging season (Stelle et al., 2008; Torres et al., 2018). We think that this indicates that the whales are starting to focus less on feeding and more on breeding. This tradeoff of foraging vs. socializing time is interesting because it comes at an energetic cost.

Clip 1. Example of social interaction between a male and female gray whale off the coast of Oregon, USA. Collected under NOAA/NMFS permit #21678

Broadly, animals need to balance the energetic demands of survival with those of reproduction. They need to reproduce to pass on their genes, but reproduction is energetically demanding, and animals also need to survive and grow to be able to reproduce. The decision to reproduce is costly because reproduction requires energetic investment and time investment since animals do not forage (gaining energy) when they are socializing. Consequently, only animals with sufficient energy reserves (i.e., body condition) to invest in reproduction actually engage in reproduction. Given these costs associated with reproduction, we expect to see a relationship between social behavior and body condition (Green, 2001) with mainly animals in good body condition engaging in social behavior because these animals have sufficient reserves to sustain the cost. Furthermore, since body condition is an indicator of foraging success and prey availability, environmental conditions can also affect social behavior and reproduction through this pathway. 

Rahman et al. (2014) used a lab experiment to study the relationship between nutritional stress and male guppy courtship behavior (Figure 1). In their experiment they tested for the effects of both decreased diet quantity and quality on the frequency of male courtship behaviors. Rahman et al (2014) found that individuals in the low-quantity group were significantly smaller than those in the high-quality group and that diet quantity had a significant effect on the frequency of courtship behaviors. Males fed a low-quantity diet performed fewer courtship behaviors. Interestingly, there was no significant effect of diet quality on courtships behavior, although there was some evidence of an interaction effect, which suggests that within the low-quantity group, males fed with high-quality food performed more courtship behaviors that those fed with low-quality food. This study is interesting because it shows how foraging success (diet quantity and quality) can affect courting behavior. 

Figure 1. A guppy (Rahman et al., 2013)

However, guppies are not the ideal species for comparison to gray whales because gray whales and guppies have quite different life history traits. A more fitting comparison would be with an example species with more in common with gray whales, such as viviparous capital breeders. Viviparous animals develop the embryo inside the body and give live birth. Capital breeders forage to build energy reserves and then rely on those energy reserves during reproduction. Surprisingly, I found asp vipers to be a good example species for comparison to gray whales.

Asp vipers (Figure 2) are viviparous snakes who are considered capital breeders because they forage prior to hibernation, and then begin reproduction immediately following hibernation without additional foraging. Naulleau & Bonnet (1996) conducted a field study on female asp vipers to determine if there was a difference in body condition at the start of the breeding season between females who reproduced or not during that season. To do this they marked individuals and measured their body condition at the start of the breeding season and then recaptured those individuals at the end of the breeding season and recorded whether the individual had reproduced. Interestingly, they found that there was a strongly significant difference in body condition between females that did and did not reproduce. In fact, they discovered that no female below a certain body condition value reproduced, meaning that they found a body condition threshold for reproduction. 

Figure 2. An asp viper

Additionally, a study on water pythons found that their body condition threshold for reproduction shifted over time in response to prey availability (Madsen & Shine, 1999). These authors found that females lowered their threshold after several consecutive years of poor prey availability. These studies are really exciting to me because they address questions that the GRANITE project team is interested in tackling.

Understanding the relationship between body condition and reproduction in gray whales is an important puzzle piece for our work. The aim of the GRANITE project is to understand how the effects of stressors on individual whales scales up to population level impacts (read Lisa’s blog to learn more). Reproduction rates play a big role in population dynamics, so it is important to understand what factors affect reproduction. Since we’re studying these whales on their foraging grounds, assessing body condition provides an important link between foraging behavior and reproduction. 

For example, if an individual’s response to a stressor is to forage less, that may lead to poorer body condition, meaning that they may be less likely to reproduce. While reduced reproduction in one individual may not have a big effect on the population, the same response from multiple individuals could impact the population’s dynamics (i.e., increasing or decreasing abundance). Understanding these different relationships between behavior, body condition, and reproduction rates is a big undertaking, but it’s exciting to be a member of the GRANITE team as this strong group of scientists works to bring together different data streams to work on this big picture question. We’re all deep into data processing right now so stay tuned over the next few years to learn more about gray whale social behavior and to find out if fat whales are more social than skinny whales. 

Did you enjoy this blog? Want to learn more about marine life, research, and conservation? Subscribe to our blog and get weekly updates and more! Just add your name into the subscribe box on the left panel.  

References

Green, A. J. (2001). Mass/Length Residuals: Measures of Body Condition or Generators of Spurious Results? Ecology82(5), 1473–1483. https://doi.org/10.1890/0012-9658(2001)082[1473:MLRMOB]2.0.CO;2

Madsen, T., & Shine, R. (1999). The adjustment of reproductive threshold to prey abundance in a capital breeder. Journal of Animal Ecology68(3), 571–580. https://doi.org/10.1046/j.1365-2656.1999.00306.x

Naulleau, G., & Bonnet, X. (1996). Body Condition Threshold for Breeding in a Viviparous Snake. Oecologia107(3), 301–306.

Rahman, M. M., Kelley, J. L., & Evans, J. P. (2013). Condition-dependent expression of pre- and postcopulatory sexual traits in guppies. Ecology and Evolution3(7), 2197–2213. https://doi.org/10.1002/ece3.632

Rahman, M. M., Turchini, G. M., Gasparini, C., Norambuena, F., & Evans, J. P. (2014). The Expression of Pre- and Postcopulatory Sexually Selected Traits Reflects Levels of Dietary Stress in Guppies. PLOS ONE9(8), e105856. https://doi.org/10.1371/journal.pone.0105856

Stelle, L. L., Megill, W. M., & Kinzel, M. R. (2008). Activity budget and diving behavior of gray whales (Eschrichtius robustus) in feeding grounds off coastal British Columbia. Marine Mammal Science24(3), 462–478. https://doi.org/10.1111/j.1748-7692.2008.00205.x

Torres, L. G., Nieukirk, S. L., Lemos, L., & Chandler, T. E. (2018). Drone up! Quantifying whale behavior from a new perspective improves observational capacity. Frontiers in Marine Science5(SEP). https://doi.org/10.3389/fmars.2018.00319