Zooming in: A closer look at bottlenose dolphin distribution patterns off of San Diego, CA

By: Alexa Kownacki, Ph.D. Student, OSU Department of Fisheries and Wildlife, Geospatial Ecology of Marine Megafauna Lab

Data analysis is often about parsing down data into manageable subsets. My project, which spans 34 years and six study sites along the California coast, requires significant data wrangling before full analysis. As part of a data analysis trial, I first refined my dataset to only the San Diego survey location. I chose this dataset for its standardization and large sample size; the bulk of my sightings, over 4,000 of the 6,136, are from the San Diego survey site where the transect methods were highly standardized. In the next step, I selected explanatory variable datasets that covered the sighting data at similar spatial and temporal resolutions. This small endeavor in analyzing my data was the first big leap into understanding what questions are feasible in terms of variable selection and analysis methods. I developed four major hypotheses for this San Diego site.

The study species: common bottlenose dolphin (Tursiops truncatus) seen along the California coastline in 2015. Image source: Alexa Kownacki.

Hypotheses:

H1: I predict that bottlenose dolphin sightings along the San Diego transect throughout the years 1981-2015 exhibit clustered distribution patterns as a result of the patchy distributions of both the species’ preferred habitats, as well as the social nature of bottlenose dolphins.

H2: I predict there would be higher densities of bottlenose dolphin at higher latitudes spanning 1981-2015 due to prey distributions shifting northward and less human activities in the northerly sections of the transect.

H3: I predict that during warm (positive) El Niño Southern Oscillation (ENSO) months, the dolphin sightings in San Diego would be distributed more northerly, predominantly with prey aggregations historically shifting northward into cooler waters, due to (secondarily) increasing sea surface temperatures.

H4: I predict that along the San Diego coastline, bottlenose dolphin sightings are clustered within two kilometers of the six major lagoons, with no specific preference for any lagoon, because the murky, nutrient-rich waters in the estuarine environments are ideal for prey protection and known for their higher densities of schooling fishes.

Data Description:

The common bottlenose dolphin (Tursiops truncatus) sighting data spans 1981-2015 with a few gap years. Sightings cover all months, but not in all years sampled. The same transect in San Diego was surveyed in a small, rigid-hulled inflatable boat with approximately a two-kilometer observation area (one kilometer surveyed 90 degrees to starboard and port of the bow).

I wanted to see if there were changes in dolphin distribution by latitude and, if so, whether those changes had a relationship to ENSO cycles and/or distances to lagoons. For ENSO data, I used the NOAA database that provides positive, neutral, and negative indices (1, 0, and -1, respectively) by each month of each year. I matched these ENSO data to my month-date information of dolphin sighting data. Distance from each lagoon was calculated for each sighting.

Figure 1. Map representing the San Diego transect, represented with a light blue line inside of a one-kilometer buffered “sighting zone” in pale yellow. The dark pink shapes are dolphin sightings from 1981-2015, although some are stacked on each other and cannot be differentiated. The lagoons, ranging in size, are color-coded. The transect line runs from the breakwaters of Mission Bay, CA to Oceanside Harbor, CA.

Results: 

H1: True, dolphins are clustered and do not have a uniform distribution across this area. Spatial analysis indicated a less than a 1% likelihood that this clustered pattern could be the result of random chance (Fig. 1, z-score = -127.16, p-value < 0.0001). It is well-known that schooling fishes have a patchy distribution, which could influence the clustered distribution of their dolphin predators. In addition, bottlenose dolphins are highly social and although pods change in composition of individuals, the dolphins do usually transit, feed, and socialize in small groups.

Figure 2. Summary from the Average Nearest Neighbor calculation in ArcMap 10.6 displaying that bottlenose dolphin sightings in San Diego are highly clustered. When the z-score, which corresponds to different colors on the graphic above, is strongly negative (< -2.58), in this case dark blue, it indicates clustering. Because the p-value is very small, in this case, much less than 0.01, these results of clustering are strongly significant.

H2: False, dolphins do not occur at higher densities in the higher latitudes of the San Diego study site. The sightings are more clumped towards the lower latitudes overall (p < 2e-16), possibly due to habitat preference. The sightings are closer to beaches with higher human densities and human-related activities near Mission Bay, CA. It should be noted, that just north of the San Diego transect is the Camp Pendleton Marine Base, which conducts frequent military exercises and could deter animals.

Figure 3. Histogram comparing the latitudes with the frequency of dolphin sightings in San Diego, CA. The x-axis represents the latitudinal difference from the most northern part of the transect to each dolphin sighting. Therefore, a small difference would translate to a sighting being in the northern transect areas whereas large differences would translate to sightings being more southerly. This could be read from left to right as most northern to most southern. The y-axis represents the frequency of which those differences are seen, that is, the number of sightings with that amount of latitudinal difference, or essentially location on the transect line. Therefore, you can see there is a peak in the number of sightings towards the southern part of the transect line.

H3: False, during warm (positive) El Niño Southern Oscillation (ENSO) months, the dolphin sightings in San Diego were more southerly. In colder (negative) ENSO months, the dolphins were more northerly. The differences between sighting latitude and ENSO index was significant (p<0.005). Post-hoc analysis indicates that the north-south distribution of dolphin sightings was different during each ENSO state.

Figure 4. Boxplot visualizing distributions of dolphin sightings latitudinal differences and ENSO index, with -1,0, and 1 representing cold, neutral, and warm years, respectively.

H4: True, dolphins are clustered around particular lagoons. Figure 5 illustrates how dolphin sightings nearest to Lagoon 6 (the San Dieguito Lagoon) are always within 0.03 decimal degrees. Because of how these data are formatted, decimal degrees is the easiest way to measure change in distance (in this case, the difference in latitude). In comparison, dolphins at Lagoon 5 (Los Penasquitos Lagoon) are distributed across distances, with the most sightings further from the lagoon.

Figure 5. Bar plot displaying the different distances from dolphin sighting location to the nearest lagoon in San Diego in decimal degrees. Note: Lagoon 4 is south of the study site and therefore was never the nearest lagoon.

I found a significant difference between distance to nearest lagoon in different ENSO index categories (p < 2.55e-9): there is a significant difference in distance to nearest lagoon between neutral and negative values and positive and neutral years. Therefore, I hypothesize that in neutral ENSO months compared to positive and negative ENSO months, prey distributions are changing. This is one possible hypothesis for the significant difference in lagoon preference based on the monthly ENSO index. Using a violin plot (Fig. 6), it appears that Lagoon 5, Los Penasquitos Lagoon, has the widest variation of sighting distances in all ENSO index conditions. In neutral years, Lagoon 0, the Buena Vista Lagoon has multiple sightings, when in positive and negative years it had either no sightings or a single sighting. The Buena Vista Lagoon is the most northerly lagoon, which may indicate that in neutral ENSO months, dolphin pods are more northerly in their distribution.

Figure 6. Violin plot illustrating the distance from lagoons of dolphin sightings under different ENSO conditions. There are three major groups based on ENSO index: “-1” representing cold years, “0” representing neutral years, and “1” representing warm years. On the x-axis are lagoon IDs and on the y-axis is the distance to the nearest lagoon in decimal degrees. The wider the shapes, the more sightings, therefore Lagoon 6 has many sightings within a very small distance compared to Lagoon 5 where sightings are widely dispersed at greater distances.

 

Bottlenose dolphins foraging in a small group along the California coast in 2015. Image source: Alexa Kownacki.

Takeaways to science and management: 

Bottlenose dolphins have a clustered distribution which seems to be related to ENSO monthly indices, and likely, their social structures. From these data, neutral ENSO months appear to have something different happening compared to positive and negative months, that is impacting the sighting distributions of bottlenose dolphins off the San Diego coastline. More research needs to be conducted to determine what is different about neutral months and how this may impact this dolphin population. On a finer scale, the six lagoons in San Diego appear to have a spatial relationship with dolphin sightings. These lagoons may provide critical habitat for bottlenose dolphins and/or for their preferred prey either by protecting the animals or by providing nutrients. Different lagoons may have different spans of impact, that is, some lagoons may have wider outflows that create larger nutrient plumes.

Other than the Marine Mammal Protection Act and small protected zones, there are no safeguards in place for these dolphins, whose population hovers around 500 individuals. Therefore, specific coastal areas surrounding lagoons that are more vulnerable to habitat loss, habitat degradation, and/or are more frequented by dolphins, may want greater protection added at a local, state, or federal level. For example, the Batiquitos and San Dieguito Lagoons already contain some Marine Conservation Areas with No-Take Zones within their reach. The city of San Diego and the state of California need better ways to assess the coastlines in their jurisdictions and how protecting the marine, estuarine, and terrestrial environments near and encompassing the coastlines impacts the greater ecosystem.

This dive into my data was an excellent lesson in spatial scaling with regards to parsing down my data to a single study site and in matching my existing data sets to other data that could help answer my hypotheses. Originally, I underestimated the robustness of my data. At first, I hesitated when considering reducing the dolphin sighting data to only include San Diego because I was concerned that I would not be able to do the statistical analyses. However, these concerns were unfounded. My results are strongly significant and provide great insight into my questions about my data. Now, I can further apply these preliminary results and explore both finer and broader scale resolutions, such as using the more precise ENSO index values and finding ways to compare offshore bottlenose dolphin sighting distributions.

Remote Sensing Applications

By Leila Lemos, PhD candidate

Fisheries and Wildlife Department, OSU

 

I am finally starting my 3rd and last year of my PhD. Just a year left and yet so many things to do. As per department requirements, I still need to take some class credits, but what classes could I take? In this short amount of time it is important to focus on my research project and on what could help me better understand the many branches of the project and what could improve my analyses. Thinking of that, both my advisor (Dr. Leigh G. Torres) and I agreed that it would be useful for me to take a class on remote sensing. So, I could learn more about this field, as well as try to include some remote sensing analyses in my project, such as sea surface temperature (SST) and chlorophyll (i.e., as a productivity indicator) conditions over the years we have collected data on gray whales off the Oregon coast.

 

Our photogrammetry data indicates that whales gradually increased their body condition over the feeding seasons of 2016 and 2018, while 2017 is different. Whales were still looking skinny in the middle of the season, and we were not collecting many fecal samples up to that point (indicating not much feeding). These findings made us wonder if this was related to delayed seasonal upwelling events and consequently low prey availability. These questions are what motivated me the most to join this class so that we might be able to link environmental correlates with our observations of gray whale body condition.

Figure 01: Skinny body condition state of the gray whale “Pancake” in August 2017.
Source: Leila S. Lemos

 

If we stop to think about what remote sensing is, we have already been implementing this method in our project since the beginning, as my favorite definition for remote sensing is “the art of collecting information of objects or phenomenon without touching it”. So, yes, the drone is a type of sensor that remotely collects information of objects (in this case, whales).

Figure 02: Drone remotely collecting information of a whale in September 2018. Drone in detail. Collected under NOAA/NMFS permit #16111.
Source: Leila Lemos

 

However, satellites, all the way up in the space, are also remotely sensing the Earth and its objects and phenomena. Even from thousands of km above Earth, these sensors are capable of generating a great amount of detailed data that is easily and freely accessible (i.e., NASA, NOAA), and can be used for multiple applications in different fields of study. Satellites are also able to collect data from remote areas like the Antarctica and the Arctic, as well as other areas that are not easily reached by humans. One important application of the use of satellite imagery is wildlife monitoring.

For example, satellite data was used to detect variation in the abundance of Weddell seals (Leptonychotes weddellii) in Erebus Bay, Antarctica (LaRue et al., 2011). Because this is a well-studied seal population, the object of this study was to test if satellite imagery could produce reliable abundance estimates. The authors used high-resolution (0.6 m) satellite imagery (from satellites Quick-Bird-2 and WorldView-1) to compare counts from the ground with counts from satellite images in the same locations at the same time. This study demonstrated a reliable methodology for further studies to replicate.

Figure 03: WorldView-1 image (0.6 m resolution) of Weddell seals hauled out east of Inaccessible Island, Erebus Bay, Antarctica.
Source: LaRue et al. (2011).

 

Satellite imagery was also applied to estimate colony sizes of Adélie penguins in Antarctica (LaRue et al., 2014). High-resolution (0.6 m) satellite imagery combined with spectral analysiswas used to estimate the sizes of the penguin breeding colonies. Ground counts were also used in order to check the reliability of the applied method. The authors then created a model to predict the abundance of breeding pairs as a function of the habitat, which was identified terrain slope as an important component of nesting density.

The identification of whales using satellite imagery is also possible. Fretwell et al. (2014)pioneered this method by successfully identifing Southern Right Whales (Eubalaena australis) in the Golfo Nuevo, Península Valdés, in Argentina in satellite images. By using very high-resolution satellite imagery (50 cm resolution) and a water penetrating coastal band that was able to see deeper into the water column, the researchers were able to successfully identify and count the whales (Fig. 04). The importance of this study was very significant, since this species was extensively hunted from the 17ththrough to the 20thcentury. Since then, the species has shown a strong recovery, but population estimates are still at <15% of historical estimates. Thus, being able to use new tools to identify, count and monitor individuals in this recovering population is a great development, especially in remote and hard to reach areas.

Figure 04: Identification of Southern Right Whales by using imagery from the WorldView2 satellite in the Golfo Nuevo Bay, Península Valdés, Argentina.
Source: Fretwell et al. (2014).

 

Polar bears (Ursus maritimus) have also been studied in the Foxe Basin, in Nunavut and Quebec, Canada (LaRue et al., 2015). Researchers used high-resolution satellite imagery in an attempt to identify and count the bears, but spectral signature differences between bears and other objects were insufficient to yield useful results. Therefore, researchers developed an automated image differencing, also known as change detection, that identifies differences between remotely sensed images collected at different times and “subtract of one image from another”. This method correctly identified nearly 90% of the bears. The technique also generated false positives, but this problem can be corrected by a manual review.

Figure 05 shows the difference in resolution of two types of satellite imagery, the panchromatic (0.6 m resolution) and the multispectral (2.4 m resolution). LaRue et al. (2015)decided not to use the multispectral imagery due to resolution constraints.

Figure 05: Polar Bears on panchromatic (0.6 m resolution) and multispectral (2.4 m resolution) imagery.
Source: LaRue et al. (2015).

 

A more recent study is being conducted by my fellow OSU Fisheries and Wildlife graduate student, Jane Dolliveron breeding colonies of three species of North Pacific albatrosses (Phoebastria immutabilis, Phoebastria nigripes, and Phoebastria albatrus)(Dolliver et al., 2017). Jane is using high-resolution multispectral satellite imagery (DigitalGlobe WorldView-2 and -3) and image processing techniques to enumerate the albatrosses. They are also using albatross species at multiple reference colonies in Hawaii and Japan (Fig. 06) to determine species identification accuracy and required correction factor(s). This will allow scientists to accurately count unknown populations on the Senkakus, which are uninhabited islands controlled by Japan in the East China Sea.

Figure 06: Satellite image of a colony of short-tailed albatrosses (Phoebastria albatrus) in Torishima, Japan, 2016.
Source: Satellite image provided by Jane Dolliver.

 

Using satellite imagery to count seals, penguins, whales, bears and albatrosses is just the start of this rapidly advancing technology. Techniques and resolutions are continuously improving. Methods can also be applied to many other endangered species, especially in remote areas, providing data on presence, abundance, annual productivity, population estimates and trends, changes in distribution, and breeding ground usage.

Other than directly monitoring wildlife, satellite images can also provide information on the environmental variables that can be related to wildlife presence, abundance, productivity and distribution.

Gentemann et al. (2017), for example, used satellite data from NASA to analyze SST variations along the west coast of the United States from 2002 to 2016. The NASA Jet Propulsion Laboratory produces global, daily, 1 km, multiscale ultra-high resolution, motion-compensated analysis of SST, and incorporates SSTs from eight different satellites. Researchers were able to identify warmer than usual SSTs (also called anomalies) along the Washington, Oregon, and California coasts from January 2014 to August 2016 (Fig.07) relative to previous years. This marine heat wave started in the Gulf of Alaska and ended in Southern California, where SST reached a maximum temperature anomaly of 6.2°C, causing major disturbances and substantial economic impacts.

Figure 07: Monthly SST anomalies in the West Coast of United States, from January 2014 to August 2016.
Source: Gentemann et al. (2017).

 

Changes in SST and winds may alter events such as the coastal upwelling that supplies nutrients to sustain a whole food chain. A marine heat-wave event as described by Gentemann et al. (2017)could have significant impacts on the health of the marine ecosystem in the subsequent season (Gentemann et al., 2017).

These findings may even relate to our questions regarding the poor gray whale body condition we noticed in 2017: this marine heat wave that lasted until August 2016 along the US west coast could have impacted the ecosystem in the subsequent season. However, I must conduct a more detailed study to determine if this heat wave was related or if another oceanographic process was involved.

So, whether remotely sensed data is generated by satellites, drones, thermal imagery, robots (as I previously wrote about), or another type of technology, it can have important  and informative applications to monitor wildlife or environmental variables associated with their ecology and biology. We can take advantage of remotely sensed technology to aid wildlife conservation efforts.

 

References

Dolliver, J., et al., Multispectral processing of high resolution satellite imagery to determine the abundance of nesting albatross. Ecological Society of America, Portland, OR, United States., 2017.

Fretwell, P. T., et al., 2014. Whales from Space: Counting Southern Right Whales by Satellite. Plos One. 9,e88655.

Gentemann, C. L., et al., 2017. Satellite sea surface temperatures along the West Coast of the United States during the 2014–2016 northeast Pacific marine heat wave. Geophysical Research Letters. 44,312-319.

LaRue, M. A., et al., 2014. A method for estimating colony sizes of Adélie penguins using remote sensing imagery. Polar Biology. 37,507-517.

LaRue, M. A., et al., 2011. Satellite imagery can be used to detect variation in abundance of Weddell seals (Leptonychotes weddellii) in Erebus Bay, Antarctica. Polar Biology. 34,1727–1737.

LaRue, M. A., et al., 2015. Testing Methods for Using High-Resolution Satellite Imagery to Monitor Polar Bear Abundance and Distribution. Wildlife Society Bulletin. 39,772-779.