Introduction

The brown smoothhound (Mustelus henlei) is an epibenthic elasmobranch ranging in shallow waters of the east Pacific Ocean from the Gulf of California to Coos Bay, Oregon, in the northern hemisphere, and from Ecuador to Peru in the southern hemisphere (Love 1996). Brown smoothhounds are found on a diversity of substrates along the continental slope, but show a preference for shallow muddy or sandy bays, such as Tomales Bay in northern California, utilizing them for feeding and refuge from larger predators (Love 1996). Throughout their range, Tomales Bay being no exception, brown smoothhounds are known to feed heavily on benthic organisms including crabs, shrimp, and polychaete worms (Haeseker and Cech 1993; Love 1996). Brown smoothhounds remain in bays throughout most of the year, but may leave them during the winter as the waters cool and rains reduce the salinity (Love 1996; Hopkins and Cech 2003).

Prior to this study, two other elasmobranch species had been tracked in Tomales Bay. Coarse-scale, population-based tracking studies indicated that time of day was the most significant influence on the movements of bat ray (Myliobatus californica) [Matern et al. 2000], and that movement direction in leopard sharks (Triakis semifasciata) was related to tide (Ackerman et al. 2000). Fine-scale, individual-based tracking (Klimley et al. 2005) demonstrated that bat rays exhibit two modes of travel: rapid and directional, and slow and non-directional. While some of the bat rays in Klimley et al. (2005) behaved consistently with the behavioral thermoregulation hypothesis presented by Matern et al. (2000), some did not, emphasizing the importance of being cognizant of individual behavior when interpreting the results of population-based telemetry studies. There have been no past studies of the movement patterns of the brown smoothhound.

Our objective was to relate the directionality and rate of movement of brown smoothhounds to environmental variables using both coarse-scale and fine-scale tracking. Coarse-scale tracking consisted of locating multiple individuals sequentially during an environmental-based sampling period (e.g. day or night, high or low tide), while fine-scale tracking consisted of constantly following a single individual for an extended period while recording the same environmental variables. We reasoned that the movement patterns of brown smoothhounds would resemble those of sympatric leopard sharks because of the overlap in their morphology (Setka and Cech 1994), habitat selection (Love 1996), and food habits (Haeseker and Cech 1993; Webber and Cech 1998). As observed in leopard sharks (Ackerman et al. 2000), we hypothesized that the brown smoothhounds would: (1) move during incoming and high tides towards the inner bay where food would be more available in the submersed littoral zone of the inner bay, and (2) that their rate of movement would be greater during nighttime than during daytime.

Materials and methods

Study site

We tracked smoothhound sharks in Tomales Bay, located north of San Francisco Bay in Marin County, California. It is a 20-km-long straight bay, averaging 1.4 km in width and 3 m in depth (Fig. 1). The cumulative area of the bay is ca. 24 km2. The bay has been divided into an inner and outer bay based on physical and chemical characteristics (Smith et al. 1991; Hopkins and Cech 2003). Longitudinal gradients of the physical and chemical properties of the water in Tomales Bay are large relative to lateral and vertical gradients, and vertical stratification is weak (Smith et al. 1991). The water in the shallow inner bay (12 km) has large seasonal and daily shifts in temperature, whereas the deeper outer bay (8 km) is more influenced by oceanic conditions and therefore the water temperature remains more constant (Matern et al. 2000). Exchange of water between the inner and outer bay during the summer is minimal, and the inner bay tends to be 2–6°C warmer and 1–3 ppt more saline than the outer bay at this time of year (Smith et al. 1991). A deep channel exists on the west side of the outer bay and is flanked by a sandbar, whereas the inner bay consists of extensive expanses of intertidal mud flats (Smith et al. 1991).

Fig. 1
figure 1

Bathymetric map of Tomales Bay and its position in California

Fish collection

We collected six adult brown smoothhounds, one male and five females, on 24–30 June 2004 (Table 1). The sharks were captured using gillnets (100 m x 3 m, 6 cm stretch mesh) set for 20–60 min perpendicular to the western shoreline of the inner and outer regions of Tomales Bay. We immediately placed captured brown smothhounds into a cooler on the boat with oxygenated sea water until transmitter implantation (5 min-1.5 h).

Table 1 Capture information for brown smoothhounds tracked in Tomales Bay. (TL = total length, M = male, F = female)

Transmitter implantation

Each shark was anesthetized by immersion in water saturated with CO2, made by mixing 40 ml concentrated acetic acid and 170 g sodium bicarbonate into 20 L of sea water (Peake 1988). Immediately after gill ventilation movements ceased (3–5 min), the shark was removed from the cooler, placed in a supine position on a portable operating table, and gills were irrigated with sea water (alternating, 1-L beakers of water with the anesthetic and sea water). We made a 2.5-cm, axial incision adjacent to the ventral midline and inserted an ultrasonic transmitter (Lotek, CAFT16-1), which was 54-mm long and 16-mm in diameter and weighed 18.9 g in water, in the peritoneal cavity and closed the incision with silk sutures and cyanoacrylate glue. After the surgery (duration <10 min), we recorded the shark’s mass, length, and sex (Table 1). Once the shark recovered from anesthesia, it was immediately released at the capture site. Three days were allowed for the individual to recover from the stress of capture and anesthesia before we began tracking.

Telemetry: coarse-scale

Coarse-scale tracking was conducted over 2.5 weeks from 29 June to 15 July 2004 during day and night and all tidal stages (incoming, high, outgoing, and low). Each shark’s transmitter emitted a unique signal, which was detected by a submerged directional hydrophone (Lotek, LHP-1) attached to the boat gunwale. The signal was converted to a radio frequency by an upconverter (Lotek, UUCN-150) and decoded by a radio receiver (Lotek, SRX-400). The receiver recorded which shark was being tracked and digitally stored its identity and its geographic coordinates obtained from a global positioning system (GPS; Garmin, GPSMAP 76S; 5 m accuracy). The range of transmitter detection was ca. 500 m, depending on ambient noise and the bay’s surface conditions. When searching for sharks, transects were made 500 m offshore along one side of the bay, submerging the hydrophone to search for the transmitters’ signal at 500-m distance intervals. Upon detecting a transmitter’s signal, we proceeded toward its source until directly overhead, indicated by the strong and non-directional nature of the transmitter’s signal (Winter 1996). We then recorded the transmitter’s identity, its position, time of day, and measured five hydrographic variables: depth (with a sounding line), substrate (using indicator clay on a sounding-line weight), and salinity, temperature, and dissolved oxygen at the surface and at ca. 3.3 m depth with a meter and probe (YSI, Model 85). We then left the shark and resumed searching for other individuals. We attempted relocation of each shark multiple times during each sampling session, allowing frequent calculations of movement directions and rates.

Data analysis: coarse-scale

Brown smoothhound positions obtained during coarse-scale tracking were projected onto a bathymetric map of Tomales Bay, using a 5 m grid digital elevation model. The digital chart was obtained from the GIS Service Center of the Information Technology Division of the California Department of Fish and Game, located in Sacramento, California. The map was displayed on a PC using a geographic information system (GIS) (ESRI, ArcView 3.2), following Klimley et al. (2005).

Tidal stages were determined using tide tables and were split into four stages: incoming, high, outgoing and low. High and low tides were defined as 1 h before and after the extreme high or low. The diel cycle was split into day and night based on the times of sunrise and sunset according to the U.S. Naval Observatory for Point Reyes Station, California. The distance to the ocean (the mouth of Tomales Bay) was calculated for every shark location using the Geofunc Microsoft Excel add-in. Using SigmaStat 3.5 statistical software (Systat Software, Inc.), two-way ANOVAs were employed to analyze the effects of the diel and tidal cycles on shark distance-to-ocean (ie. movement direction), followed by a Tukey pairwise multiple comparison test when significant differences were detected. To boost our sample size per treatment for the distance-to-ocean analyses we grouped tidal-stage shark positions into high/incoming and low/outgoing stages.

Shark movement rates were calculated by dividing distance traveled between successive positions by the time interval between position readings. The movement rate data was analyzed in cm sec−1 after a natural-log transformation, but reported in the text as m sec−1 to aid comprehension. Two-way ANOVAs were then used to analyze differences in movement rates between tidal stages and between diel stages, followed by a Tukey pairwise multiple comparison test when significant differences were detected.

Telemetry: fine-scale

In fine-scale tracking, we continuously tracked a single shark (BS2 or BS3 depending on the tracking session) for two 20–21 h intervals and one 4-h interval over a period from 1–13 July 2004. Tracking intervals were chosen to include as many combinations of tidal and diel stages as possible. After locating our target shark for the tracking session, we positioned our boat directly above it and followed it continuously, recording the time of day and the geographical coordinates every 6–15 min with a few exceptions (i.e. the shark could not be detected or we returned to the dock for refueling/crew exchange). During fine-scale tracking sessions, the same five hydrographic variables were measured every hour.

Data analysis: fine-scale

Brown smoothhound positions recorded during fine-scale tracking were imported into a GIS (ESRI, ArcView 3.2) in the same manner described for coarse-scale analysis. In order to describe the general patterns of fine-scale brown smoothhound movements, consecutive points were selected from each track and classified by polygraphic shape: points in no apparent order, but separated, were classified as continuous, non-uniform swimming (CS) movement; points that were linear were classified as straight-line swimming (SS) movement; points that were very close to one another (<10 m h−1) were classified as non-movement (NM). We used software (Oriana 3) to calculate circular statistics that quantified the directionality of the polygraphic shapes in each tracking session, indicated by the coefficient of angular concentration (Rayleigh’s ‘r’). An r value of 1 indicates maximum directionality, whereas 0 indicates omni-directionality. Furthermore, Rayleigh’s uniformity test was used to determine whether the distributions of directional movements were significantly non-uniform. Because of our low, fine-scale tracking sample size (N =2), we did not analyze the effects of the tidal and diel cycles on the distribution of BS2 and BS3 as we did with our coarse-scale tracking data.

Our methods for computing fine-scale rates of movement and the statistical analyses of these data were the same as described above for the coarse-scale analysis, except that the fine-scale data could not be transformed to fit a normal distribution. Thus, a Kruskal-Wallis ANOVA on ranks was used, followed by Dunn’s pairwise multiple comparison test when significant differences were detected. The data from tracks BS3-A and -B were combined before analysis.

Results

Coarse-scale tracking

We gathered a total of 121 locations for all six smoothhounds during coarse-scale tracking. Whereas there was no apparent pattern to the sharks’ movement directions based on the diel cycle (Fig. 2a), a pattern did exist for the tidal cycle (Fig. 2b). Mean shark positions were located significantly further from the mouth of the bay during high/incoming tides than during low/outgoing tides (P = 0.014) (Table 3). Though, BS6 was the only individual to exhibit this directional movement at a significant level in Tukey multiple pairwise comparisons (P = 0.01; Table 4).

Fig. 2
figure 2

Mean locations for individual sharks through the diel (on the left; part a) and tidal (on the right; part b) cycles overlain on a bathymetric map of a section of Tomales Bay. The sample size for each symbol is indicated on the map

Table 2 Two-way ANOVA results from analyses of coarse-scale tracking of brown smoothhound movements in Tomales Bay. Rate of movement (ROM) data were converted to cm s−1 then natural-log transformed for the analyses. Means reported are least square means. Bold text indicates a statistically significant result.
Table 3 Two-way ANOVA results from analyses of coarse-scale tracking of brown smoothhound movements in Tomales Bay. Rate of movement (ROM) data were converted to cm s−1 then natural-log transformed for the analyses. Means reported are least square means. Bold text indicates a statistically significant result.
Table 4 Tukey’s test results for the analysis of brown smoothhound positions in the Tomales Bay relative to the tidal cycle. Means reported are least square means. Bold text indicates a statistically significant result. Overall P-value = 0.014. (H/I = high/incoming tides, L/O = low/outgoing tides)

The sharks’ movement rates ranged over 0.01–0.34 m sec−1, with an overall mean of 0.09 m sec−1. The sharks exhibited significantly different rates of movement in night (mean = 0.10 m sec−1) vs. day (mean = 0.07 m) (P = 0.024) and across the tidal cycle (P = 0.003) (Table 2). However, the only significant difference in movement rates between tidal stages occurred between outgoing and low tides (P = 0.008) and outgoing and incoming tides (P = 0.029) (Table 5). Movement rates between all other tidal combinations were statistically indistinguishable (Table 5). Due to the time gap between positions during coarse-scale tracking, these data represent minimal movement rates.

Table 5 Tukey’s test results for the coarse-scale rate of movement analysis through the tidal cycle. Overall P-value = 0.003. Bold text indicates a statistically significant result. (O = outgoing tide, L = low tide, H = high tide, I = incoming tide)

Hydrographic gradients we measured were small throughout the study site for the duration of our study. Dissolved oxygen levels (mean: 8.5 mg O2 L−1) were quite homogeneous throughout the bay with the lowest reading (7.5 mg O2 L−1) being well above stressful levels (Davis 1975; Lai et al. 1990). Similarly, temperature (mean = 19.2°C, range = 17.7–19.9°C), salinity (mean = 34.2 ppt, range = 33.8–34.4 ppt), and substrate composition (95% of measurements were mud) all showed minimal variation within the bay. However, because we only measured these variables in waters used by the sharks and not what was available but not used, we cannot draw any correlations between the movements of the sharks and these hydrographic variables.

Fine-scale tracking

Two brown smoothhounds were tracked on a fine scale over three tracking periods. BS3 was tracked from 10:15 to 14:11 on 1 July 2004 and from 14:54 on 7 July to 11:46 on 8 July 2004. BS2 was tracked from 13:24 on 12 July to 09:13 on 13 July 2004. Both sharks displayed the three patterns previously defined: CS, SS, NM. Figure 3 shows the fine-scale track of BS3 in its entirety, exhibiting examples of all three patterns of movement. We determined the degree of directionality to the segments of the tracks that we segregated into the three categories (Table 6). As expected, no directionality was indicated for the segments of the track exhibiting little movement while more directionality was exhibited while actively swimming. Shark BS2 exhibited no significant directionality to NM1, CS1, and CS2 (Table 6). The same shark showed a significant degree of directionality during SS1 and near-significant directionality during SS2. Similarly, the r-values for SS1 and SS2 were considerably higher than those of NM1, CS1, CS2. BS3 showed no significance in directionality for NM2, CS3, CS4, CS5, or SS3. However, BS3 did show significant directionality during SS4 and SS5. Furthermore, SS4 and SS5 had greater r-values than CS4 or CS5.

Fig. 3
figure 3

A fine-scale track of brown smoothhound BS3 in Tomales Bay. Dark circles designate nighttime positions; white circles daytime positions. This figure illustrates all three swimming classifications exhibited during our fine-scale tracking: 1) continuous, non-uniform swimming (CS), 2) straight-line swimming (SS), and 3) no movement (NM). The rectangular inset shows an expanded view of NM2 when the shark appeared to be motionless. The shark first exhibited Straight-line Swimming (SS5) at sunset during the change of the tide to outgoing. The inlaid circular histogram indicates strong directionality during SS5. The shark performed continuous, non-uniform swimming (CS5) from night to day

Table 6 Directionality of movements for individual brown smoothhounds tracked on a fine scale based on classification of movement polygons. Bold text indicates a statistically significant result. (A = Day 1, B = Day 2, CS = continuous, non-uniform swimming, SS = straight-line swimming, NM = non-movement, r = Rayleigh’s coefficient of angular concentration, P = probability of uniform distribution)

Both BS2 and BS3 displayed movement rates drastically higher for nighttime than during daytime, but only BS2 exhibited differences between tidal stages (Table 7). For BS2, movement rates averaged 0.18 and 0.12 m sec−1 for nighttime and daytime, respectively. For BS3, movement rates averaged 0.20 m sec−1 and 0.10 m sec−1 for nighttime and daytime, respectively. BS2’s mean movement rates during outgoing and low tides were near equal at 0.19 m sec−1 and 0.19 m sec−1, and both were significantly higher than those during incoming tide (O vs. I: Diff of ranks = 28.97, Q = 4.834, P < 0.05; L vs I: Diff of ranks = 27.81, Q = 3.805, P < 0.05; O vs. L: Diff. of ranks = 1.16, Q = 0.18, P > 0.05). Only one location of BS2 was obtained during high tide, thus, high tide was excluded in the analysis of its movement rates. Unlike BS2, BS3’s rates of movement did not differ among tidal stages. Only two locations of BS3 were obtained during low tide, thus, low tide was excluded in the analysis of its rates of movement.

Table 7 Kruskal-Wallis ANOVA on ranks results for the fine-scale tracking of sharks BS2 and BS3 in Tomales Bay. Due to low sample sizes, high tide was excluded in the analysis of movement rates of BS2 and low tide was excluded in the analysis of movement rates of BS3. Bold text indicates a statistically significant result. (A = Day 1, B = Day 2, ROM = rate of movement)

BS2 was located over substrates consisting of mud (85% of samples), gravel/algae (2%), mud/algae (5%), gravel (5%), algae (2%) and sand (1%) in water column depths of 2.6–4.7 m. At the surface, the salinity recorded over the shark was relatively constant at 34 ppt, while dissolved oxygen ranged over 89.2–99.5% of air saturation, and temperature ranged over 19.3–20.1°C. At 3.3 m depth, salinity stayed constant at 34 ppt, while dissolved oxygen saturation fluctuated over 88.9–99.8%, and temperature ranged over19.4–20.1°C.

BS3 spent 95% of its time over a muddy substrate and 5% of its time over a sand/clay substrate in depths ranging from 2.0 to 9.6 m. At the surface, salinity was 34 ppt, dissolved oxygen fluctuated over 95.4–101.7% of air saturation, and temperature ranged over 16.5–19.6°C. At 3.3 m below the surface, salinity was 34 ppt, dissolved oxygen fluctuated over 95.4–107.0%, and temperature ranged over 16.5–19.4°C.

Discussion

Our results from these two, distinct methods of tracking elucidated the complex movement patterns of brown smoothhounds in Tomales Bay, California. For coarse-scale tracking data, the tidal cycle appeared to be the dominant factor affecting their distribution, whereas both the diel and tidal cycles affected their movement rates. At the individual level, (i.e. fine-scale tracking) the diel cycle affected the movement rates of both BS2 and BS3, while the tidal cycle influenced the movement rate of only BS2. Their movements were contained within waters with weak temperature, salinity, and dissolved oxygen gradients typical of the inner bay in summer (Smith et al. 1991). Tracking studies of other elasmobranchs in Tomales Bay also endeavored to relate movements to diel and tidal cycles. Matern et al. (2000) found a distinct periodicity to bat rays’ movements between the inner bay and the outer bay, hypothesized to be related to the large temperature gradient between the two areas, but their movement rates did not correlate with either diel or tidal cycles. In contrast, Ackerman et al. (2000) found that leopard sharks showed strong movements with tides and increased their movement rates at night compared to day. As hypothesized, the brown smoothhounds displayed similar movement patterns to those of leopard sharks. Interestingly, brown smoothhound and leopard sharks often associate with the same substrates (Love 1996; Ackerman et al. 2000), have similar diet compositions (Haeseker and Cech 1993; Webber and Cech 1998), and share almost identical blood and muscle characteristics (Setka and Cech 1994).

Sharks were located farthest toward the inner bay during high and incoming tides. These movements into freshly flooded mud flats on incoming and high tides are consistent with what is known about their foraging ecology. Brown smoothhounds are known to feed heavily on benthic organisms including crabs, shrimp, and polychaete worms (Haeseker and Cech 1993; Love 1996), many of which are widely available in the tidal mud flats when they are flooded. As the sharks moved inward and outward they generally swam in shallow water, stayed near shore, and may have been foraging. Alternatively their behavior also may reflect predator avoidance (e.g., minimizing contact with larger sharks). Sevengill sharks are present in the bay and they grow to be twice as long as brown smoothhounds, which have been recorded as their prey (Love 1996).

The similarities in the body design and feeding habits of smoothhounds to other sharks that move with tides supports our findings that brown smoothhound sharks also may move in relation to tidal stage. In temperate, estuarine habitats on the eastern U.S. coast, the dusky shark and the sandbar shark moved predominately in the direction of tidal currents (Huish and Benedict 1978; Medved and Marshall 1983). Ackerman et al. (2000) explained the tidally influenced movement on leopard sharks as a means of reducing energy expenditure. Our findings suggest that brown smoothhounds’ movements in Tomales Bay are associated with tide to some degree, but not to the degree of the leopard shark. The movement rates through the tidal cycle observed in this study are more difficult to interpret.

From the coarse-scale data, smoothhounds exhibited a maximum rate of movement on an outgoing tide, mid-range rates of movement during high tide, and minimum rates of movement during incoming and low tide. The observed behavior may perhaps be explained by their foraging ecology. The mid-range rates of movement exhibited during high tides may reflect active foraging behavior over mud flats. During an outgoing tide it may be beneficial to a shark to rapidly move off these flats to avoid being stranded. At low tide brown smoothhounds have generally retreated to deeper water where they may await their next foraging opportunity while digesting their last meal, thus exhibiting slower rates of movement. Looking at the individual level, both sharks we tracked exhibited some patterns not clear at the population level. Although fine-scale tracking likely yielded our most accurate movement rates, our sample sizes across the tidal spectrum were inadequate to fully interpret those results. However, across the diel cycle, our sample size doubled per treatment as treatment variables halved versus the tidal cycle (i.e. day and night), permitting comparison to the coarse-scale results.

As we expected, at a coarse-scale, our brown smoothhounds showed a significant increase in their movement rates at night versus day. At a fine scale, the difference in movement rate across the diel cycle was even more distinct. The sharks showed a marked increase in their movement rates at night, with BS2 moving 50% faster and BS3 moving 100% faster than during the day. These rates were nearly equivalent to those found in leopard sharks by Ackerman et al (2000). Based on our field observations during both tracking modes, the higher resolution of the individual tracking method seems to be the more accurate and precise measure of the smoothhounds’ actual movement rates over the diel cycle.

Movement polygons

Among the movement patterns obtained from fine-scale tracking, CS was the most common. Our primary hypothesis for the motive behind this movement pattern was foraging. In our field observations, and when plotted on the bathymetric map, continuous non-uniform swimming was generally displayed just off the shoreline and on mudflats where smoothhounds are known to feed (Haeseker and Cech 1993).

Straight-line (SS) swimming was defined as a period of movement in which individuals showed directionality by a high ‘r’ coefficient of concentration. Using fine-scale tracking, Klimley (1993) described highly directional swimming in scalloped hammerhead sharks (Sphyrna lewini) in the Sea of Cortez and in bat rays (Myliobatis californica) in Tomales Bay (Klimley et al. 2005). Regarding the cues that brown smoothhounds use for directing its movement, Tomales Bay’s high turbidity precludes effective use of vision for directionality, suggesting cues other than vision are at play. Brown smoothhound SS occurred more often during a moving tide, as opposed to a slack tide, in Tomales Bay (SS1, SS2, SS5), perhaps suggesting the use of rheotactic orientation and passive transport. However, Klimley (1993) showed that the heading of a flow-oriented shark would change unless maintained by forward momentum. Klimley et al. (2005) suggested that bat rays swimming in straight-line manner may have been guided by the strong linear magnetic field that runs parallel with the length of Tomales Bay, which was created by a large geological (San Andreas) fault. It is plausible that brown smoothhounds use the same magnetic cues.

The combination of late maturity, low fecundity, and slow growth in most shark species has made them especially vulnerable to over-exploitation, and shark populations around the world are declining at a rapid rate due to overfishing (Musick 1999). In Tomales Bay, female brown smoothhounds are known to birth young throughout the summer (Hopkins and Cech 2003). A better understanding of the brown smoothhounds’ movements provides natural resource managers with more information to aid in managing and protecting California’s small commercial fishery for brown smoothhounds. While the small scale of this study precludes a detailed assessment of the environmental variables influencing the movements of brown smoothhounds, it is apparent that variables such as tide and diel cycle influence the distribution of the species within inshore habitats. Further large-scale studies are necessary to tease apart the interactions of brown smoothhounds with their environment and to assess the possible impact of anthropogenic activities.