Printer Friendly

Population estimation, distribution, and habitat preference of Irrawaddy dolphins Orcaella brevirostris (Owen in Gray, 1866) in the Brunei Bay, Malaysian waters.

Delphinids (Delphinidae) represent one of the most social and complex groups of mammals, which have cognitive competencies present just in some mammal species (Whitehead et al. 2000). As top-level predators, they play a unique role in the structure and function of marine communities (Whitehead et al. 2000). They also have long life spans, late maturity, low reproduction rates and extended parental care (Taylor 2002). Therefore, they are incapable of enduring the increasing rates of anthropogenic mortality (Beasley 2007). These characteristics lead to slow rates of population growth and vulnerability to rapid population declines (Taylor 2002).

Population size is a vital aspect of the ecology of any species, and its estimation represents an indispensable component in the management of wildlife (Williams et al. 2002). For instance, the estimation of sustainable reduction of animal population levels necessitates the knowledge of the abundance and their variances (Wade 1998). In general, precise estimation of the cetacean population is difficult, costly and time-consuming (Taylor and Gerrodette 1993). The sampling design and environmental variability are crucial factors because they can affect our ability to estimate cetacean population sizes and trends of variation (Thompson et al. 2000). Careful survey plans and investigation of the assumptions are fundamental for reasonably accurate population estimates (Read et al. 2003).

Cetacean distribution is influenced by a list of environmental factors i.e. physico-chemical and climatological variables, biotic factors (competition and predation), and anthropogenic causes (fishing activities and boat traffic) (MacLeod et al. 2004). Depending on the geographical area involved, interactions between these different aspects may vary (Jefferson et al. 1993). Fluctuations in temperature and salinity of different water bodies have an especially great impact on the distribution of cetaceans (Baumgartner and Mate 2005). In addition, cetacean distribution patterns are influenced by the seasons and time of the day due to seasonal abundance of resources, presence of other species and habitat structure (Daura-Jorge et al. 2005). Substantial social knowledge is essential for group members to recognize each other, distinguish their rank in the social hierarchy, form and uphold alliances, and engage in and interpret an extensive range of social behaviors (Herman 1991). Group size and composition can fluctuate depending on the species, location, prey availability, predation risk, age, gender and reproductive status of an individual in a group (Michaud 2005).

Among the diversity of habitats that are inhabited by dolphins in the marine environment, coastal areas are considered most vulnerable to anthropogenic activities (McIntyre 1999). Thus, some species of coastal dolphins are among the most threatened species of cetaceans (Parra 2005), and the Irrawaddy dolphin, Orcaella brevirostris (Owen in Gray, 1866) is one such example. Sir Richard Owen first described this species in 1866, based on a specimen found in 1852 in the harbor of Visakhapatnam on the east coast of India (Sinha 2004). Irrawaddy dolphins have been delineated as 'facultative' river cetaceans considering their adaptability to inhabit both marine and freshwater environments (Leatherwood and Reeves 1994). According to the International Union for Conservation of Nature (IUCN) Red List of Threatened species, they are categorized as Vulnerable (Reeves et al. 2008). In general, their population is decreasing (Reeves et al. 2008). Five subpopulations of Irrawaddy dolphin are already categorized as Critically Endangered in their habitats located in the Ayeyarwady River (Smith 2004), Mahakam River (Jefferson et al. 2008), Malampaya Sound (Smith and Beasley 2004a), Mekong River (Smith and Beasley 2004b) and Songkhla Lake (Smith and Beasley 2004c). Irrawaddy dolphins inhabit the shallow, near shore tropical and subtropical rivers, and in marine waters of the Indo-West Pacific (Dolar et al. 1997). These dolphins have been observed from Visakhapatnam to the deltas of Brahmaputra and Ganges Rivers in India (James et al. 1989). They have also been observed in coastal water bodies of Bangladesh, Brunei Darussalam, Cambodia, India, Indonesia, Lao People's Democratic Republic, Malaysia, Myanmar, Philippines, Singapore, Thailand and Vietnam (Pilleri and Gihr 1974, Perrin et al. 1995, Smith et al. 1997, Stacey and Leatherwood 1997, Parra 2005, Kreb et al. 2007, Reeves et al. 2008, Sutaria and Marsh 2011, Hines et al. 2015). Several studies on Irrawaddy dolphins were conducted already in Malaysia such as distribution and sighting rate estimated by Bali et al. (2008), Jaaman (2010) and Kuit et al. (2014); abundance estimated by Minton et al. (2013), Woan et al. (2013); acoustic studied by Hoffman et al. (2016); habitat characteristics and critical areas investigated by Peter et al. (2016); and behavioral observations by Kamaruzzan et al. (2011). However, the Irrawaddy dolphin population has not been studied yet in Brunei Bay and the dolphin population in the bay is under threat by several anthropogenic activities, i.e. accidental bycatch in fisheries (Jaaman et al. 2009), water quality and habitat degradation due to land clearing for coastal development (Long 2014), untreated domestic sewage from the catchment areas (Yau 1988), intense shipping and maritime activities (Eng 1992), and overfishing (Silvestre and Garces 2004). In order to contribute to the long-term conservation and management of this species, it was necessary to carry out an ecological assessment of the population in the Bay. To this aim, the present study was conducted to achieve the following specific objectives: 1) to estimate the population size of Irrawaddy dolphin in the Brunei Bay, Malaysia; 2) to quantify group size and composition (i.e. the presence and number of calves and adults); 3) to assess the spatial occurrence of Irrawaddy dolphins; 4) to assess the daytime (morning and afternoon) and seasonal occurrences of Irrawaddy dolphins; 5) to assess how habitat abiotic factors (sea depth, surface water temperature, salinity, turbidity and nearby coast distance) may impact the occurrence of Irrawaddy dolphins.

Material and methods

Study area

The study was conducted at the Malaysian side of Brunei Bay, which is important as a nursery, foraging and transient ground for sea turtles, dugongs and coastal cetaceans (Rajmani and Marsh 2010). About 30% of Brunei Bay is in Brunei territory, and 70% is in Malaysia, the joint domain of Sarawak, Sabah and Federal Territory of Labuan. The coastal and marine ecosystems in this bay comprise mangrove forests, seagrass beds, estuarine systems, mudflats and coral reefs (Bali 2005, Bujang et al. 2006, Jaaman et al. 2010). Also, this bay has a continental slope (Ahmad-Kamil et al. 2013), where the continental shelf is 50-70 km wide and underlained by 8-10 km of siliciclastic sediments (Straub et al. 2012). At a sea depth of ~200 m the continental shelf-slope break occurs, and the seabed descends steeply to reach the floor of the Borneo Trough at a sea depth of 2800 m (Straub and Mohrig 2009). According to Mohamed and Landner (1993), in Brunei Bay the water residence time (WRT) is two weeks, stratification occurs at approximately 6 m depth; due to a 2-3 ppt salinity difference, an upwelling is observed in the inshore waters. The sea surface temperature (SST) ranges from 24.60 to 32.02oC (Hee and Suratman 2016), and air average humidity throughout the year is 82% (Hogan 2011). Brunei Bay has a high quantity of fish resources (Joseph et al. 2016) and 54 species of fish mentioned by Mohamed and Landner (1993). Fishing is the second most important economic activity in the Bay (Department of Fisheries Sabah 2008). A total of 78 species of phytoplankton and 80 different soft bottom fauna species are reported in this Bay (Mohamed and Landner 1993). For its ecological uniqueness and economic significance, Brunei Bay is a high-priority area for research and conservation of marine organisms (Joseph et al. 2016).


Line-transect boat surveys (Hiby and Hammond 1989, Beasley et al. 2013) were conducted in the Bay from April 2013 to October 2015 (Fig. 1). A 10 m long speedboat with 40 and 35 hp engines was used to survey the inner part of the Bay. The speed of the boat was approximately 12 km [h.sup.-1]. The transect lines for each year were different to cover the whole Bay (Fig. 1). Transect lines were designed within the Malaysian territory. In the surveys, limitations were imposed by local marine security during crossings of borders in state waters and the international harbor in the bay. Surveys were conducted within 5 km from the coastline during four periods of the year (January, April, July and end of September or beginning of October). In Malaysia, the four periods are corresponding to the Northeast Monsoon (November-March) and Southwest Monsoon (May--September) with two transitional periods between the two monsoons, generally known as the Inter-Monsoon seasons (IMS), occurring in April and October (Malaysian Meteorological Department 2008). For each season, 10-12 days were dedicated to the survey (depending on sea conditions) with each covering the whole bay (except the Brunei part). Surveys were generally carried out between the 7:00 to 15:00 h range due to increasing winds in the afternoon. While the boat was moving along the transect line, there was a team of five observers; one observer located at the bow of the boat searched for dolphin groups ahead; two observers on each side of the boat made observations with the aid of 7 X 50 binoculars; and two additional observers maintained a constant search of the near area around the boat with the naked eye. Parra (2005) defined a school of dolphins as the case where a group of dolphins have relatively close spatial cohesion (i.e. the distance between individuals did not exceed 100 m). Once a group was encountered, the search effort stopped and the position of the boat was marked for site information and recorded with a marine handheld GPS receiver. An electronic compass barometer was used to collect information on weather and the Beaufort Sea state from the start of the transect line. This information was recorded again whenever conditions changed during the rest of the transect line. A depth meter was used to measure sea depth (m); a portable turbidity meter to measure turbidity (Nephelometric Turbidity Unit, NTU); YSI Multiprobe meter to measure sea surface temperature ([degrees]C) and surface water salinity (ppt). The survey boat would follow a dolphin or group by keeping a safe distance of 20 to 50 m from them, to identify species, assess group size and take photographs of the dorsal fins of the adults through two cameras. A calf was defined as a particularly smaller individual which is closely associated with another larger (adult) dolphin. The number of calves was counted in a group by the naked eyes, and they were not considered for this population estimation using mark--recapture photo-ID because calves are not born with marks in their dorsal fins (Ryan et al. 2011). The calves were observed to have grown in length by 59% (96 cm to 1.53 m) and 266% (12.3 to 45 kg) in weight during the first seven months (Tas'an et al. 1980), and they achieve their adult size in 4 to 6 years (Stacey and Arnold 1999).

Data analysis

Photo identification

Right side dorsal fin photographs of adult Irrawaddy dolphins were imported in Photos 1.5 for iOS 2015, a photograph-editing program. Images were organized by date, then cropped and sorted by putting rank (Q1-Q6) based on their quality (Gowans and Whitehead 2001). The highest quality considered as Q6 is when notches and niches of the dorsal fins were clearly visible and at a good focus (Towner et al. 2013). When the quality of photo-ID decreased the rank also decreased (Towner et al. 2013). Based on dorsal fin photographic assessment criteria of Gowans and Whitehead (2001), a total of 1078 photographs (right side of dorsal fin) of Q5, Q4 and Q3 were imported into DARWIN 2.22 for IOS 2011 dorsal fin ID software (Wilkin et al. 1998, Barata and Brooks 2005, Beirao et al. 2014) (Table 1). Well-marked (having nicks in dorsal fins) photographs were considered for analysis. Any unmarked photograph of dorsal fins was not considered for mark--recapture analysis. Each dorsal fin was traced with a fixed spacing of the points along the leading and trailing edges. After that, the dorsal fin was compared to the entire catalogues and obtained a rank based on the probability of the match. When fins matched, it meant that those fins were already in the catalogue with ID codes. Before considering a new fin the matching was confirmed by looking at fin outline and marks. Once a new fin was found, it was assigned to a unique ID code (first sighting date and time) and then added to the catalogue. From the surveys 106 photographs of left side dorsal fins were captured and not considered for the analysis.

Capture--recapture analysis

To assess the abundance, data from DARWIN was transferred to Microsoft Excel. Then Program MARK 8.1 was used to perform mark--recapture analyses of the individual sighting histories of identity (White and Burnham 1999, Kreb 2004, Smith et al. 2006, Sutaria and Marsh 2011, Minton et al. 2013). We had 11 encounter occasions from 2013-2015, and it was impossible to define secondary occasions. A close population model was not a plausible assumption here for several reasons i.e. mortality could occur over the survey period, the bay is open, and Irrawaddy dolphins were travelling to the bay area of Brunei part. In those cases, the open-population POPAN parameterization with Cormack--Jolly-Seber (CJS) model was used to estimate population parameters (Schwarz and Arnason 1996, 2009). The parameter <f> refers to apparent survival rate, p is the probability of the observation, b is the probability that an individual from the super-population would enter the survey area between capture occasions, and N is the estimated population size (Schwarz and Arnason 1996). Super-population refers the sum of the observable (inside of study area) and unobservable (outside of study area) individuals (Kendall 2008). The subscripts . and t in model notation exemplify constant and time-dependent parameters respectively (Lebreton et al. 1992). Maximum likelihood models were used to estimate population parameters (Towner et al. 2013). For goodness-of-fit (GOF) tests of the CJS model, Program RELEASE was used for validating model assumptions (Burnham et al. 1987). Models were assembled for consolidations of time-dependence and consistency for each parameter, and the most suitable model was chosen using the small sample corrected Akaike information criterion ([AIC.sub.C]) (Burnham and Anderson 1998). Based on the GOF results of TEST 2+ TEST 3 in RELEASE, a post hoc variance inflation factor (C) could be calculated to conform for extra-binomial variation in the data resulting in a quasi-Akaike information criterion ([QAIC.sub.C]).

Open capture-recapture models

To avoid initiating bias in estimates of abundance, survival, recapture and entrance probability using capture--recapture methods, it is necessary that model assumptions are met. In this study the open model (POPAN parameterization with CJS) model assumptions (Lebreton et al. 1992) were the following: 1) natural marks carried by dolphins during this study should not be lost or missed, 2) natural marks carried by all individuals should be accurately identified during recaptures, 3) individuals should be released quickly after being captured, 4) sampling sessions should be of shorter duration compared to total duration of the sampling period, 5) all live marked dolphins available on each sampling occasion should have equal capture probabilities, 6) survival probabilities, and 7) individual dolphins from super-population should have higher chance to enter the survey area between capture occasions. To validate the first assumption, strict quality controls were adopted during data analysis and the only photo-ID images used to identify dolphins were high quality (Q3, Q4, Q5) with distinctive fins. Our current estimates were not produced using images of individuals with no marks so including them in future analysis might help to contextualize our current results. To validate the second assumption and minimize human errors in matching dorsal fins over the years, regular checks were conducted by only the most experienced researchers (SAJ or XZ). For the third assumption, dolphins were not physically removed during each photographic session but simply captured in a photo-identification image (Silva et al. 2009). The time spent photographing dolphins was much shorter in duration compared to the time spent searching for dolphins in between photographic sessions, so this assumption was respected. The fourth assumption was also met having sampling occasions of 10-12 days, which is much shorter of duration than the interval between occasions (three months). To verify the fifth, sixth and seventh assumptions were met using POPAN parameterization to estimate these parameters. Equal recapture probabilities were partially met because mixing between sampling periods occurred, with individuals observed leaving our study areas returning sometime later. However, movements outside of the bay (in the Brunei part) showed that some individuals have home ranges that extend beyond areas of our survey effort, making them unavailable for capture during this study. The assumption of equal survival probability was likely achieved in this study because no commercial or illegal hunting (of dolphins) occurred in Malaysian coastal shelf waters. The seventh assumption was partially met, where dolphins travelled to the Malaysian part from the bay of Brunei side, which could be attributed to foraging and socializing purposes.

Group composition

To investigate the group composition of Irrawaddy dolphins in relation to group size, the composition of groups was first calculated as percentages of adults and calves. The percentages of adults were then grouped into group sizes, i.e. small (1-5 individuals), medium (6-10 individuals) and large (11-20 individuals). For this purpose, the Kruskal-Wallis test was performed and a post hoc comparison Mann--Whitney pairwise test applied.


The dolphins' spatial occurrence map and sighting density map were generated to locate hotspots for Irrawaddy dolphins using ArGIS 9 (ArcMap ver. 9.3) by importing all sighting GPS coordinates, to test if the Irrawaddy dolphins homogenously occurred in the Bay. To compare the mean number of sightings per hours of survey effort between morning and afternoon, the Mann--Whitney test was applied. Furthermore, to investigate if all groups of Irrawaddy dolphins had the same occurrence, the mean numbers of individuals were based on the seasons [IMS (April), SWM, IMS (October) and NEM] by performing the Kruskal-Wallis test and the post hoc comparison (Mann-Whitney pairwise test). This statistical test was performed by running PAST 3.04 for iOS 2014.

Habitat preference

To investigate if this dolphin population had preferences in terms of sea depth, surface water temperature, salinity, turbidity and distance from the coast, a chi-square test of conformity was performed using R ver. 3.1.3 for iOS 2014 to compare the number of individual occurrences in different ranges of abiotic factors. For this, all the abiotic factor parameters were stratified into bins of five units (e.g. 0.1-4.99 m, 5-9.99 m for depth; 27-27.99[degrees]C, 28-28.99[degrees]C for sea surface temperature; 0.1-4.99 ppt, 5-9.99 ppt for surface water salinity; 0NTU, 0.1-9.99NTU for turbidity; 0.1-1.49 km, 1.5-2.99 km for nearby coast distance from occurrence). The Euclidean distance function (ESRI 1996) was used to compute the shortest distance to the nearest shoreline from the sighting position. However, the survey effort was randomly distributed with respect to abiotic factors, the dolphins' sighting, and occurrences. For all statistic tests, a = 5% were chosen as significance criterion.


Survey effort and distance covered by the observers

A total of 3574.92 km of surveys were undertaken in 129 days from 2013 to 2015 in the Bay with a total survey effort of 297.91 hours. Survey effort was varied during the years and seasons due to different transect lines, which were followed for different years and changing sea state (Beaufort scale) (Table 2).

Population parameters estimation (capture--recapture)

A total of 36 individual adult Irrawaddy dolphins were photo-identified from 2013 to 2015. We tested eight CJS models with POPAN parameterization, where model {<f> [p.sub.t] b} did not converge. For the remaining seven models, model choice criteria and parameter estimates are shown in Table 3. On the basis of [AIC.sub.C] scores, model {[phi] [p.sub.t] [b.sub.t]} (constant survival, time dependent probability of capture and probability of entry) was the most parsimonious, and the estimated population size (N) was 33 (95% confidence interval = 28-39, SE = 3) (Table 3). In addition, the estimated constant survival ([phi]-hat) was 0.98 (0.89-0.99, SE = 0.01), and the time dependent observation rate (p-hat) was 0.27 (0.14-0.45, SE = 0.07) (Table 4). The estimated time dependent individual entry rate from super-population (b-hat) was 0.15 (0.10-0.22, SE = 0.03) (Table 4). Based on the result of TEST 2 + TEST 3 in Program RELEASE (Table 5), a variance inflation factor of C = 0.83 was estimated and applied, reflecting under-dispersion in the data (Cooch and White 2009). There was no violation of underlying open-population mark--recapture assumptions, and all TEST results indicated equal survival probabilities (SR) and seasonal migrations (SM) among photographically captured dolphins where no significant p-value was found in the TEST results (Table 5).

Entire population estimation

According to Jaaman (2010) following the Irrawaddy dolphins was sometimes very challenging because they are shy and evasive, often changing direction or swimming away when the survey boat approached them. In Brunei Bay, 11.10% of the observed animals were not amenable to photo-ID due to several reasons: (a) sometimes dolphins were disappeared completely from water surface before capturing their dorsal fin photos due to the boat approaching, (b) poor quality dorsal fin photos (Q2, Q1) were taken due to bad weather conditions (Beaufort sea state 4). Also, 12.02% of the sampled animals' photographs (Q3-Q5) were not amenable to mark--recapture analysis due to lack of the presence of a mark in their dorsal fins. Hence the population was estimated for the 0.88 proportion that was well marked. Dividing the estimate (N) by this proportion yielded a result of 38 as an estimate of the adult population. Applying the same on 95% CI (confidence interval) gave a result of 38 (95% CI = 32-44) for the adult population. Also, once we included calves using the same approach we got an estimate of 41 (95% CI = 36-49) for the entire population including calves.

Group size and composition

A total of 47 dolphin groups were encountered during the survey period. Groups of Irrawaddy dolphins varied in size from 1 to 18 animals, with a mean of 6 (SE = 0.66). In the group composition we had observed zero to one calf in small and medium-sized groups, and one to two calves in large size groups. The mean proportion of adults and calves in the group were 91.25% [+ or -] 1.84% and 8.75% [+ or -] 1.84% respectively. No statistically significant differences were observed comparing the group composition under different group sizes. However, the observations of Irrawaddy dolphin calves in the groups are a positive indicator that the dolphins are breeding successfully in the bay area and provide hope that the population will remain stable or increase in numbers.


The survey team encountered 47 dolphin groups at 36 GPS locations (sightings). Dolphins were encountered more than once near the GPS locations of Lawas and Labuan Island. Dolphins did not occur homogeneously in the Brunei Bay, Malaysian waters. Two hotspots were identified there, one near Lawas and the other near Labuan Island. The number of sightings (sighting rate 2.8-12.3 [km.sup.-2]) (Fig. 2a) and the number of individuals' occurrence was higher in both areas (Fig. 2b). However, there were no significant differences between seasons by comparing the mean number of sightings per hours of survey effort. Also, no significant differences were found between seasons by comparing the mean number of individuals' occurrence. Regarding daytime occurrence, there were 30 sightings in the morning and six in the afternoon. No statistically significant differences were found by comparing the mean number of sightings per hours of survey effort between morning and afternoon. Survey efforts were higher in the morning than the afternoon, with 217.93 hours in the morning and 79.98 in the afternoon. Comparing the mean number of individuals per group, adults per group and calves per group, there were no significant differences observed between seasons.

Habitat preference

In the current study, the range of sea depth was 2 to 30.4m during the dolphins' encounters. Statistically significant differences ([chi square] = 40.26, df = 5, p-value < 0.001) were observed for dolphin occurrences under different depth categories. Numbers of encounters and occurrences were higher between 2 to 9.99 m depths (Table 6). The range of surface water temperatures during the survey periods was 28 to 32[degrees]C and significant differences were observed for individual occurrences under different water temperature ranges ([chi square] = 12.67, df = 5, p-value = 0.03). The highest number of encounters occurred between 29 to 31.99[degrees]C (Table 6). The range of surface water salinity was 0.27 to 27.42 ppt during the encounter period. There were no statistically significant differences observed for individual occurrences under different categories of salinity (Table 6). No significant differences were found statistically under different turbidity ranges (Table 6). Also, no encounters happened when nearby coastline distance was < 1 km and > 7 km, and maximum occurrences were from 1.5 to 4.49 km. A statistically significant difference of occurrence was found at different coastline distances from dolphins' occurrences ([chi square] = 27.23, df=4, p-value < 0.001) (Table 6).


Photo-identification was validated as an indispensable tool in the evaluation of population size, occurrence, social organization, distribution and migration patterns of many species of cetaceans since 1970 (Hammond et al. 1990, Whitehead et al. 2000). Irrawaddy dolphins are elusive and display unobtrusive behavior at the water surface; during a slow rolling dive only the upper-most dorsal surface of the animal becomes visible (Smith 2009). The absence of distinctive marks, cryptic surfacing and research vessel avoidance by the dolphins are the main problems for the photo identification techniques (Smith et al. 2004), so higher effort is needed for their photo-identification (Lloze 1973, Dhandapani 1992, Kreb 1999). Similar problems had occurred for our survey team in the Brunei Bay. In this case, the computer program DARWIN was used and found suitable for dorsal fin matching (based on fin shape and outline) and cataloging, but this program also had infrequent considerable errors to rank fins. This flaw was considered as minor when compared to the naked eye matching (Towner et al. 2013).

No historical data for the abundance of Irrawaddy dolphin exist in the Brunei Bay population. In the current study, the estimated entire population size was 41 (95%; CI = 36-49), which is comparatively higher than Songkhla Lake population but lower than the population of Sundarban Mangrove Forest, Kuching Bay, Bangpakong Estuary, Balikpapan Bay, Coastal water bodies of Bangladesh, Chilika Lake, Banten Bay, Cowie Bay and Malampaya sound (Table 7). The estimated mean group size and group size range of Irrawaddy dolphins are similar to the findings of Parra (2005) and Dolar et al. (1997), but higher than Smith et al. (2006) and Kreb and Budiono (2005a) (Table 8).

Out of 36 sightings of the dolphins, most of the encounters happened near the Lawas area and Labuan Island. One of the possible reasons could be fish availability in these areas. In fact, during the survey period the survey team observed higher numbers of fishing trawlers near these areas compare to the other sides of the bay. Dolphins have a greater risk of injury and mortality through the entanglement in the trawlers' fishing nets. According to Read (2008), fishing-related mortality of small cetaceans is considered the most severe and immediate threat worldwide. Based on DOF Sabah (2008), Brunei Bay is one of the major prawn trawling grounds of Sabah. Fishing time in the bay varies and is often determined by weather (Matsumoto 2007). In the Brunei Bay, Matsumoto (2007) identified 106 species of fishes, crustaceans, bivalves, gastropods and cephalopods, including 12 families of fishes and three families of cephalopods with a total of 29 species that are preferred by Irrawaddy dolphins as prey (Heinsohn 1979, Marsh et al. 1989). There were no dolphin occurrences in the east of the Brunei Bay, near Weston and Menumbok. The probable cause for this may be the presence of a developed industrial (offshore oil platform), mega coastal development projects and an international harbor (busy water traffics) in those areas. These infrastructures may have disturbed and chased out Irrawaddy dolphins from the areas as in Banten Bay, Indonesia (Khalifa et al. 2014), but more study is needed to determine that. Furthermore, this dolphin species exhibits high sensitivity to noise pollution, i.e. roaring sound of the speedboat and ship engines, excessive ship and boat traffic and coastal construction works (Stacey and Hvenegaard 2002, Hashim and Jaaman 2011). Noise likely distracts them because most of the dolphin activities are depending on sound perception and biosonar (Van Parijs et al. 2000, Hoffman et al. 2016). However, the dolphins show an adaptive behavior to the negative stimuli by doing deep diving, short time surfacing for breathing (Stacey and Hvenegaard 2002) and changing the travelling direction (Kreb et al. 2012).

Seasonally, the mean number of sightings did not vary significantly between seasons nor did the mean number of individuals per group vary between seasons. These results may indicate that Irrawaddy dolphins inhabit in Brunei Bay the whole year, with no seasonal differences. Similar results have been observed in Cleveland Bay, Australia (Parra 2005), as well as in Banten Bay, Indonesia for Irrawaddy dolphins (Khalifa et al. 2014). Considering daytime, the number of sightings was higher in the morning compared to the afternoon but not significantly. The factor behind this was the bad weather conditions (wind, swell and higher Beaufort level) in the afternoon, which resulted in less survey effort compared to the morning. In general, the swell of Beaufort Sea state 4 decreases the ability to sight dolphins at a distance which creates a bias. This dolphin occurs most often from medium to low water level (Smith et al. 2004). In this study, surveys have been done in the high tide time. Beasley (2007) also noticed a higher occurrence of dolphins in the slow current.

In Brunei Bay, Malaysia, Irrawaddy dolphins showed habitat preferences (e.g. 2 to 19.99 m sea depths, 29 to 31.99[degrees]C surface water temperature, and 1.5 to 4.49 km coastline distance), and the numbers of encounters were significantly different at different ranges of some abiotic factors. Dolar et al. (2002) mentioned that the distribution of Irrawaddy dolphins was restricted to sea depths of < 15 m in Malampaya Sound, Indonesia. In the East Kalimantan, this dolphin occurred at sea depths of 5-14 m (Kreb and Budiono 2005b). In Banten Bay, Indonesia, Khalifa et al. (2014) observed this dolphin at sea depths of 1-9 m. In Cleveland Bay, Australia they occurred at water depths of < 15 m depth (Parra 2005). Geographically Brunei Bay is in the tropical region and air temperature does not fluctuate much. Seasons do not change the surface water temperature much in the Bay. Peter et al. (2016) observed Irrawaddy dolphins with a mean sea surface temperature of 30.42 [+ or -] 0.61[degrees]C in Kuching Bay, Sarawak, and our result is similar to this. Smith et al. (2004) found that the dolphins prefer low salinity conditions, and in the waterways of Sundarbans Mangrove Forest of Bangladesh the occurrence of dolphins were dramatically decreased in high salinity areas (Smith et al. 2006). In Malaysian Kuching Bay, Peter et al. (2016) also found habitat preference of this dolphin for lower water salinity. In Brunei Bay water salinity is variable because there are eight river mouths entering the bay, which also bring sediments from upstream, often increasing turbidity (Howes and NWPO 1986). Also in Brunei Bay, the occurrence of dolphins was lower in high turbid waters, which is an opposite finding compared to those of other researchers (Stacey and Arnold 1999, Smith et al. 2004). Other studies linked Irrawaddy dolphins to shallow areas close to river mouths and changing tidal states, features which are likely to be associated with higher levels of turbidity (Dolar et al. 2002, 2006, Peter et al. 2016). In Cleveland Bay, Australia, 50% of the sightings of Irrawaddy dolphins were within 15 km of the coast (Parra 2005). Also, Parra et al. (2002) observed this dolphin on the east coast of Queensland within 10 km of the nearest land. These findings for nearby coastline distance from dolphins' occurrence are similar to our result in the Brunei Bay.

The present study provides the first detailed information about population abundance, spatial occurrence and habitat preferences of Irrawaddy dolphins in the Brunei Bay. This study may serve as an important baseline for future comparisons and will help researchers, conservationists, local marine park managers and policy makers in developing effective conservation and management plans for the area. Additionally, the current estimated data should also be added to inform the IUCN Red List of Threatened Species Criteria (IUCN 2012a, b), so that status assessment of Irrawaddy dolphin population in Brunei Bay could be performed. In the future, population genetics can be studied, which will help to obtain more information about the Irrawaddy dolphin population in the bay. Furthermore, a long-term research program is needed where the research team can continue (at least seasonally) to monitor the dolphins and develop a complete knowledge of the population status and trends of Irrawaddy dolphins in the bay. This bay is not only a good habitat for Irrawaddy dolphins, but also for other marine mammals (Table 9), so a long-term conservation plan and management are necessary to reduce the anthropogenic activities in the bay. For this, there are several strategies the state governments can consider:

* local governments can announce the Lawas area and Labuan Island as a Marine Protected Areas (MPA). However, conflict between bordering state and country waters may challenge this.

* Reduction of fishery activities near Lawas area and Labuan Island could be enforced.

* Speed limits for boats and ships, and specified water routes could be implemented.

* Focus group discussions with local fisherman communities could be developed to create awareness about the importance of marine mammals.

Acknowledgements--The authors would like to give special thanks to the team of the Marine Endangered Species (MES) Program, INOS-UMT for the database. The corresponding author also would like to acknowledge the European Commission for providing a scholarship (2014-2016) under the Erasmus Mundus Masters Course in Tropical Biodiversity and Ecosystems (TROPIMUNDO). Thanks to Eleonore R. A. Viez (BEVT, ULB) for helping in data analysis; Mostofa Kamrul Ahsan and Tim Marshal Hunter (Dept. of Computer Science, NDSU) for language revision. Also, thanks to Mr Ismail (Boatman) and his family members for cooking for us during the survey periods.

Permits--The authors extend their thanks to the Sarawak Forestry Department for giving permission to conduct research on biological resources in Sarawak waters [permit no. NCCD.907.4.4(JLD.11)-35].

Funding--The study was funded by Ministry of Education (MOE) Malaysia, Higher Institution Centre of Excellence (HICOE) Grant Scheme (2013-2015) of INOS-UMT; FIO-UMT surveys for marine mammals and sea turtles in the Bay of Brunei, 2015-2018 project and the China-ASEAN maritime cooperation fund.

Conflict of interest--The authors declare that they have no conflict of interest.


Ahmad-Kamil, E. I. et al. 2013. The effects of water parameters on monthly seagrass percentage cover in Lawas, east Malaysia. --Sci. World J. 2013: 1-8.

Baird, I. G. and Beasley, I. L. 2005. Irrawaddy dolphin Orcaella brevirostris in the Cambodian Mekong River: an initial survey. --Oryx 39: 301-310.

Bali, J. 2005. Marine resources of Kuala Lawas and their utilization. --Proc. 7th Ann. Workshop of the National Parks and Wildlife Division, Hornbill Vol. 7.--Sarawak For. Dept, Kuching, Sarawak, Malaysia.

Bali, J. et al. 2008. Aerial sighting rate of marine life in Sarawak waters.--7th Int. Sci. Symp. of the Intergovernmental Oceanographic Commission-Western Pacific (IOC/WESPACT), Sabah, pp. 1-13.

Barata, T. and Brooks, S. P. 2005. Dolphins who's who: a statistical perspective.--In: Kalviainen, H. et al. (eds), Image analysis. Springer, pp. 429-438.

Baumgartner, M. F. and Mate, B. R. 2005. Summer and fall habitat of North Atlantic right whales (Eubalaena glacialis) inferred from satellite telemetry.--Can. J. Fish. Aquat. Sci. 62: 527-543.

Beasley, I. L. 2007. Conservation of the Irrawaddy dolphin Orcaella brevirostris (Owen in Gray, 1866) in the Mekong River: biological and social considerations influencing management. --PhD thesis, James Cook Univ.

Beasley, I. L. et al. 2007. Review of the status and conservation of Irrawaddy dolphins Orcaella brevirostris in the Mekong Rver of Cambodia, Lao PDR and Vietnam.--In: Smith, B. D. et al. (eds), Status and conservation of freshwater populations of Irrawaddy dolphins. Wildl. Conserv. Soc, pp. 67-82.

Beasley, I. L. et al. 2013. Likely future extirpation of another Asian river dolphin: the critically endangered population of the Irrawaddy dolphin in the Mekong Rver is small and declining. --Mar. Mamm. Sci. 29: 226-252.

Beirao, L. et al. 2014. Performance of computer-assisted photographic matching of Guiana dolphins (Sotalia guianensis). --Aquat. Mamm. 40: 313-316.

Bujang, J. S. et al. 2006. Distribution and significance of seagrass ecosystems in Malaysia.--Aquat. Ecosyst. Health Manage. 9: 203-214.

Burnham, K. P. and Anderson, D. R. 1998. Model selection and inference: a practical information-theoretic approach. --Springer.

Burnham, K. P. et al. 1987. Design and analysis methods for fish survival experiments based on release--recapture.--Am. Fish Soc. Monogr. 5: 1-437.

Cooch, E. G. and White, G. C. 2009. Program MARK: a gentle introduction.--Colorado State Univ.

Daura-Jorge, F. G. et al. 2005. Seasonal and daily patterns of group size, cohesion and activity of the estuarine dolphin, Sotalia guianensis (PJ van Beneden)(Cetacea, Delphinidae), in southern Brazil.--Rev. Bras. Zool. 22: 1014-1021.

DOF, Dept of Fisheries Sabah 2008. Summary of annual fisheries statistics Sabah 2008.--<> accessed 10 February 2016.

Dhandapani, P. 1992. Status of the Irrawaddy river dolphin Orcaella brevirostris in Chilka Lake.--J. Mar. Biol. Assoc. India 34: 90-93.

Dolar, M. L. L. et al. 1997. Preliminary investigation of marine mammal distribution, abundance, and interactions with humans in the southern Sulu Sea.--In: Morton, B. and Perrin, W. F. (eds), Asian marine biology. Hong Kong Univ. Press, pp. 61-81.

Dolar, M. L. L. et al. 2002. Preliminary report on a small estuarine population of Irrawaddy dolphins Orcaella brevirostris in the Philippines.--Raffles Bull. Zool. 50: 155-160.

Dolar, M. L. L. et al. 2006. Distributional ecology of cetaceans in the central Philippines.--J. Cetacean Res. Manage. 8: 93-111.

Eng, L. P. 1992. Water quality in the coastal areas of Brunei Darussalam: status, management issues and recommendations. --In: Silvestre, G. et al. (eds), The coastal resources of Brunei Darussalam: status, utilization and management. ICLARM Conf. Proc, pp. 91-108.

Gowans, S. and Whitehead, H. 2001. Photographic identification of northern bottlenose whales (Hyperoodon ampullatus): sources of heterogeneity from natural marks.--Mar. Mamm. Sci. 17: 76-93.

Hammond, P. S. et al. 1990. Individual recognition of cetaceans: use of photo-identification and other techniques to estimate population parameters.--Int. Whaling Commission.

Hashim, N. A. N. and Jaaman, S. A. 2011. Boat effects on the behavior of Indo-Pacific humpback (Sousa chinensis) and Irrawaddy dolphin (Orcaella brevirostris) in Cowie Bay, Sabah, Malaysia.--Sains Malays. 40: 1383-1392.

Hee, Y. Y. and Suratman, S. 2016. Physico-chemical parameters profile during dry and wet seasons in southern South China Sea: Brunei Bay.--Asian J. Chem. 28: 2146-2152.

Heinsohn, G. E. 1979. Biology of small cetaceans in North Queensland waters.--The Great Barrier Reef Marine Park Authority.

Herman, L. M. 1991. What the dolphin knows, or might know, in its natural world.--In: Pryor, K. and Norris, K. S. (eds), Dolphin societies: discoveries and puzzles. Univ. of California Press, pp. 349-368.

Hiby, A. R. and Hammond, P. S. 1989. Survey techniques for estimating the abundance of cetaceans.--Rep. Int. Whaling Commission 11: 47-80.

Hines, E. M. et al. 2015. Line transect estimates of Irrawaddy dolphin abundance along the eastern Gulf Coast of Thailand. --Front. Mar. Sci. 2: 1-10.

Hoffman, J. M. et al. 2016. Description of whistles of Irrawaddy dolphins (Orcaella brevirostris) from the waters of Matang, Peninsular Malaysia.--Bioacoustics 25: 1-10.

Hogan, M. C. 2011. The encyclopedia of Earth.--<>

Howes, J. R. and N. W P. O. 1986. Evaluation of Sarawak wetlands and their importance to water birds: joint project between National Parks and Wildlife Office of the Sarawak Forest Office (NPWO) and Interwarder.--East Asia/ Pacific Shorebird Study Programme: report 4, Limbang/Lawas Districts of Brunei Bay, WWF Malaysia.

IUCN 2012a. Guidelines for application of IUCN red list criteria at regional and national levels: ver. 4.0.--IUCN Species Survival Commission, Gland.

IUCN 2012b. IUCN red list categories and criteria: ver. 3.1, 2nd edn.--IUCN Species Survival Commission, Gland.

IUCN 2015. The IUCN Red List of Threatened Species, ver. 2015-4. IUCN Species Survival Commission, Gland.

Jaaman, S. A. 2010. Marine mammals in East Malaysia: dstribution and interactions with fisheries.--VDM Verlag Dr. Muller Aktiengesellschaft & Co. KG.

Jaaman, S. A. et al. 2009. The magnitude and sustainability of marine mammal by-catch in fisheries in East Malaysia.--J. Mar. Biol. Assoc. UK 89: 907-920.

Jaaman, S. A. et al. 2010. A preliminary sighting record of shorebirds in Sundar Estuary, Lawas, Sarawak.--Proc. 10th Annu. Workshop of the National Parks and Wildlife Division, Hornbill Vol. 10, Sarawak Forestry Corporation (SFC).

James, P. S. B. R. et al. 1989. On the mortality and stranding of marine mammals and turtles at Gahirmatha, Orissa from 1983 to 1987.--J. Mar. Biol. Assoc. India 31: 28-35.

Jefferson, T. A. et al. 1993. Marine mammals of the world.--FAO species identification guide, pp. 320.

Jefferson, T. A. et al. 2008. Orcaella brevirostris (Mahakam River subpopulation). (errata version published in 2016) The IUCN Red List of Threatened Species 2008: e.T39428A988 42174. - <> accessed 11 November 2017.

Joseph, J. et al. 2016. Genetic stock compositions and natal origin of green turtle (Chelonia mydas) foraging at Brunei Bay. --Global Ecol. Conserv. 6: 16-24.

Junchompoo, C. et al. 2014. Population and conservation status of Irrawaddy dolphins (Orcaella brevirostris) in Trat Bay, Trat Province, Thailand.--Proc. Des. Symp. Conserv. Ecosyst. 2: 32-38.

Kamaruzzan, A. S. et al. 2011. Effect of water parameters on the behaviour of Indo-Pacific humpback and Irrawaddy dolphins in Cowie Bay, Sabah, Malaysia.--Borneo Sci. 28: 1-7.

Kendall, W L. 2008. The robust design.--In: Cooch, E. and White G. C. (eds), Program MARK: a gentle introduction. Colorado State Univ., pp. 15-01-15-40.

Khalifa, M. A. et al. 2014. Preliminary study on the distribution of Irrawaddy dolphin, Orcaella brevirostris, in Banten Bay. --Open J. Mar. Sci. 4: 338-343.

Kittiwattanawong, K. et al. 2007. Review of the status and conservation of Irrawaddy dolphins Orcaella brevirostris in Songkhla Lake of Thailand.--In: Smith, B. D. et al. (eds), Status and conservation of freshwater populations of Irrawaddy dolphins. Wildl. Conserv. Soc, pp. 83-89.

Kreb, D. 1999. Observations on the occurrence of Irrawaddy dolphin, Orcaella brevirostris, in the Mahakam River, East Kalimantan, Indonesia.--Z. Saugetierkunde 64: 54-58.

Kreb, D. 2004. Abundance of freshwater Irrawaddy dolphins in the Mahakam River in East Kalimantan, Indonesia, based on mark--recapture analysis of photo-identified individuals.--J. Cetacean Res. Manage. 6: 269-277.

Kreb, D. 2008. Conservation and diversity of cetaceans in and near Balikpapan Bay, East Kalimantan, Indonesia.--Technical Final Report.

Kreb, D. and Budiono 2005a. Conservation management of small core areas: key to survival of a Critically Endangered population of Irrawaddy river dolphins Orcaella brevirostris in Indonesia.--Oryx 39: 178-188.

Kreb, D. and Budiono 2005b. Cetacean diversity and habitat preferences in tropical waters of East Kalimantan, Indonesia. --Raffles Bull. Zool. 53: 149-155.

Kreb, D. et al. 2007. Status and conservation of Irrawaddy dolphins Orcaella brevirostris in the Mahakam River of Indonesia.--In: Smith, B. D. et al. (eds), Status and conservation of freshwater populations of Irrawaddy dolphins. Wildl. Conserv. Soc, pp. 53-64.

Kreb, D. et al. 2012. Abundance and threats monitoring surveys during low water levels July and September 2012.--Yayasan Konservasi RASI, Samarinda, Indonesia.

Kuit, S. H. et al. 2014. Cetacean research and a precautionary approach in developing dolphin-watching tourism in the coastal waters of the Matang mangroves.--Proc. Matang Mangrove For. Manage. Conf., Perak, Malaysia, pp. 27-39.

Leatherwood, S. and Reeves, R. R. 1994. River dolphins: a review of activities and plans of the Cetacean Specialist Group. --Aquat. Mamm. 20: 137-154.

Lebreton, J. D. et al. 1992. Modeling survival and testing biological hypotheses using marked animals: a unified approach with case studies.--Ecol. Monogr. 62: 67-118.

Lloze, R. 1973. Contribution a l'etude anatomique histologique et biologique de l' Orcaella brevirostris (Gray, 1866) (Cetacea, Delphinidae) du Mekong.--PhD thesis, LUniv. Paul Sabatier de Toulouse III.

Long, S. M. 2014. Sarawak coastal biodiversity: a current status. --Kuroshio Sci. 8: 71-84.

MacLeod, K. et al. 2004. Seasonal distribution of minke whales Balaenoptera acutorostrata in relation to physiography and prey off the Isle of Mull, Scotland.--Mar. Ecol. Prog. Ser. 277: 263-274.

Mahmud, A. I. et al. 2018. Factors influencing the behaviour of Irrawaddy dolphins Orcaella brevirostris (Owen in Gray, 1866) in Brunei Bay, Malaysia.--J. Ethol. 36: 169-180.

Malaysian Meteorological Department 2008. General climate of Malaysia.--<>.

Marsh, H. et al. 1989. Irrawaddy dolphin Orcaella brevirostris (Gray, 1866).--In: Harrison, R. J. and Ridgway, S. H. (eds), Handbook of marine mammals. Academic Press, pp. 101-118.

Matsumoto, B. M. M. 2007. Fish and fisheries resources. --In: Saleem, M. and Ejria, S. (eds), Coastal environmental profile of Brunei Bay, Sabah. Univ. Malaysia Sabah, pp. 95-133.

McIntyre, A. D. 1999. Conservation in the sea--looking ahead. --Aquat. Conserv. Mar. Freshwater Ecosyst. 9: 633-637.

Michaud, R. 2005. Sociality and ecology of the odontocetes.--In: Ruckstuhl, K. E. and Neuhaus, P. (eds), Sexual segregation in vertebrates: ecology of the two sexes. Cambridge Univ. Press, pp. 303-326.

Minton, G. et al. 2013. Population estimates and distribution patterns of Irrawaddy dolphins (Orcaella brevirostris) and Indo-Pacific finless porpoises (Neophocaena phocaenoides) in the Kuching Bay, Sarawak.--Raffles Bull. Zool. 61: 877-888.

Mohamed, M. and Landner, L. 1993. Environmental impact assessment of a mixed tropical hardwood integrated pulp and paper mill--a case study.--Environ. Impact. Assess. Rev. 13: 353-374.

Parra, G. J. 2005. Behavioural ecology of Irrawaddy, Orcaella brevirostris (Owen in Gray, 1866), and Indo-Pacific humpback dolphins, Sousa chinensis (Osbeck, 1765), in northeast Queensland, Australia: a comparative study.--PhD thesis, James Cook Univ.

Parra, G. J. and Corkeron, P. J. 2001. Feasibility of using photo-identification techniques to study the Irrawaddy dolphin, Orcaella brevirostris (Owen in Gray 1866).--Aquat. Mamm. 27: 45-49.

Parra, G. J. et al. 2002. Distribution of Irrawaddy dolphins, Orcaella brevirostris, in Australian waters.--Raffles Bull. Zool. 10: 141-154.

Perrin, W. F. et al. 1995. Osteological comparison of Bryde's whales from the Philippines with specimens from other regions.--Rep. Int. Whaling Comm. 46: 409-413.

Peter, C. et al. 2016. Identifying habitat characteristics and critical areas for Irrawaddy dolphin, Orcaella brevirostris: implications for conservation.--In: Das, I. and Tuen, A. A. (eds), Naturalists, explorers and field scientists in southeast Asia and Australasia. Springer, pp. 225-238.

Pilleri, G. and Gihr, M. 1974. Contribution to the knowledge of the cetaceans of southwest and monsoon Asia (Persian Gulf, Indus Delta, Malabar, Andaman Sea and Gulf of Siam). --Investigations of Cetacea 5: 95-153.

Rajmani, L. and Marsh, H. 2010. Using parallel regional- and local-scale initiatives to inform conservation management of rare wildlife: a case study of the dugong Dugong dugon in Sabah, Malaysia.--Endanger. Species Res. 13: 17-23.

Read, A. J. 2008. The looming crisis: interactions between marine mammals and fisheries.--J. Mamm. 89: 541-548.

Read, A. J. et al. 2003. Abundance of bottlenose dolphins in the bays, sounds, and estuaries of North Carolina.--Mar. Mamm. Sci. 19: 59-73.

Reeves, R. R. et al. 2008. Orcaella brevirostris. The IUCN Red List of Threatened Species. - < UK.2008.RLTS.T15419A4579987.en> accessed 5 January 2016.

Ryan, G. E. et al. 2011. Irrawaddy dolphin demography in the Mekong River: an application of mark--resight models. --Ecosphere 2: 58.

Schwarz, C. J. and Arnason, A. N. 1996. A general methodology for the analysis of capture--recapture experiments in open populations.--Biometrics 52: 860-873.

Schwarz, C. J. and Arnason, A. N. 2009. Jolly--Seber models in MARK.--In: Cooch, E. G. and White, G. C. (eds), MARK: a gentle introduction Program. Colorado State Univ.

Silva, M. A. et al. 2009. Estimating survival and abundance in a bottlenose dolphin population taking into account transience and temporary emigration.--Mar. Ecol. Progr. Ser. 392: 263-276.

Silvestre, G. T. and Garces, L. R. 2004. Population parameters and exploitation rate of demersal fishes in Brunei Darussalam (1989-1990).--Fish. Res. 69: 73-90.

Sinha, R. K. 2004. The Irrawaddy dolphins Orcaella brevirostris of Chilika Lagoon, India.--J. Bombay Nat. Hist. Soc. 101: 244-251.

Smith, B. D. 2004. Orcaella brevirostris (Ayeyarwady River subpopulation). The IUCN Red List of Threatened Species 2004: e.T44556A10919593. - <> accessed 11 November 2017.

Smith, B. D. 2009. Irrawaddy dolphin (Orcaella brevirostris).--In: Perrin, W F. et al. (eds), Encyclopedia of marine mammals. Academic Press, pp. 638-642.

Smith, B. D. and Beasley, I. 2004a. Orcaella brevirostris (Malampaya Sound subpopulation). The IUCN Red List of Threatened Species 2004: e.T44187A10858619. - <> accessed 11 November 2017.

Smith, B. D. and Beasley, I. 2004b. Orcaella brevirostris (Mekong River subpopulation). The IUCN Red List of Threatened Species 2004: e.T44555A10919444. - <> accessed 11 November 2017.

Smith, B. D. and Beasley, I. 2004c. Orcaella brevirostris (Songkhla Lake subpopulation). The IUCN Red List of Threatened Species 2004: e.T44557A10919695. - <> accessed 11 November 2017.

Smith, B. D. et al. 1997. Investigation of cetaceans in the Ayeyarwady River and northern coastal waters of Myanmar.--In: Morton, B. and Perrin, W F. (eds), Asian Marine Biology. Hong Kong Univ. Press, pp. 173-194.

Smith, B. D. et al. 2004. Status, ecology and conservation of Irrawaddy dolphins (Orcaella brevirostris) in Malampaya Sound, Palawan, Philippines.--J. Cetacean Res. Manage. 6: 41-52.

Smith, B. D. et al. 2005. New information on the status of finless porpoises Neophocaena phocaenoides and Irrawaddy dolphins Orcaella brevirostris in Bangladesh and Myanmar.--Scientific Committee Document SC/57/SM4. International Whaling Commission.

Smith, B. D. et al. 2006. Abundance of Irrawaddy dolphins (Orcaella brevirostris) and Ganges River dolphins (Platanista gangetica gangetica) estimated using concurrent counts made by independent teams in waterways of the Sundarbans mangrove forest in Bangladesh.--Mar. Mam. Sci. 22: 527-547.

Smith, B. D. et al. 2007. Review of the status and conservation of Irrawaddy dolphins Orcaella brevirostris in the Ayeyarwady River of Myanmar.--In: Smith, B. D. et al. (eds), Status and conservation of freshwater populations of Irrawaddy dolphins. Wildl. Conserv. Soc, pp. 21-39.

Stacey, P. J. and Arnold, P. W 1999. Orcaella brevirostris.--Mamm. Species 616: 1-8.

Stacey, P. J. and Hvenegaard, G. T. 2002. Habitat use and behavior of Irrawaddy dolphin (Orcaella brevirostris) in the Mekong River of Laos.--Aquat. Mamm. 28: 1-13.

Stacey, P. J. and Leatherwood, S. 1997. The Irrawaddy dolphin, Orcaella brevirostris: a summary of current knowledge and recommendations for conservation action.--In: Morton, B. and Perrin, W F. (eds), Asian Marine Biology. Hong Kong Univ. Press, pp. 195-214.

Straub, K. M. and Mohrig, D. 2009. Constructional canyons built by sheet-like turbidity currents: observations from offshore Brunei Darussalam.--J. Sediment Res. 79: 24-39.

Straub, K. M. et al. 2012. Architecture of an aggradational tributary submarine channel network on the continental slope offshore Brunei Darussalam.--In: Prather, B. E. et al. (eds), Application of the principles of seismic geomorphology to continental-slope and base-of-slope systems: case studies from seafloor and near-seafloor analogues. Soc. Sedimentary Geol. Spec. Publ. 99: 13-30.

Sutaria, D. and Marsh, H. 2011. Abundance estimates of Irrawaddy dolphins in Chilika Lagoon, India, using photo-identification based mark--recapture methods.--Mar. Mamm. Sci. 27: 338-348.

Taylor, B. L. 2002. Conservation biology.--In: Perrin, W. F. et al. (eds), Encyclopedia of marine mammals. Academic Press, pp. 273-276.

Taylor, B. L. and Gerrodette, T. 1993. The uses of statistical power in conservation biology: the vaquita and northern spotted owl. --Conserv. Biol. 7: 489-500.

Tas'an et al. 1980. Orcaella brevirostris (Gray, 1866) from Mahakam River.--Jaya Ancol Oceanarium, pp. 60.

Thompson, P. M. et al. 2000. Combining power analysis and population viability analysis to compare traditional and precautionary approaches to conservation of coastal cetaceans. --Conserv. Biol. 14: 1253-1263.

Tongnunui, S. et al. 2011. Preliminary investigation of Irrawaddy dolphin (Orcaella brevirostris) in the Bangpakong Estuary, Inner Gulf of Thailand.--Environ. Nat. Resour. J. 9: 48-57.

Towner, A. V. et al. 2013. Gauging the threat: the first population estimate for white sharks in South Africa using photo identification and automated software.--PloS One 8(6): e66035.

Van Parijs, S. M. et al. 2000. Sounds produced by Australian Irrawaddy dolphins, Orcaella brevirostris.--J. Acoust. Soc. Am. 108: 1938-1940.

Wade, P. R. 1998. Calculating limits to the allowable human-caused mortality of cetaceans and pinnipeds.--Mar. Mamm. Sci. 14: 1-37.

White, G. C. and Burnham, K. P. 1999. Program MARK: survival estimation from populations of marked animals.--Bird Study 4: 120-138.

Whitehead, H. et al. 2000. Science and the conservation, protection, and management of wild cetaceans.--In: Mann, J. et al. (eds), Cetacean societies: field studies of dolphins and whales. Univ. of Chicago Press, pp. 308-332.

Wilkin, D. J. et al. 1998. Query by sketch in DARWIN: digital analysis to recognize whale images on a network.--In: Yeung, M. M. et al. (eds), Storage and retrieval for image and video databases VII. SPIE, pp. 41-48.

Williams, B. K. et al. 2002. Analysis and management of animal populations: modeling estimation and decision making. --Academic Press.

Woan, T. S. et al. 2013. A preliminary study of the population size of Irrawaddy dolphins (Orcaella brevirostris) in Cowie Bay, Sabah, Malaysia.--J. Trop. Biol. Conserv. 10: 23-26.

Yau, K. H. 1988. Water quality of Brunei River and estuary. --ASEAN/US Technical Workshop on Integrated Tropical Coastal Zone Management, Singapore.

Anisul Islam Mahmud, Saifullah Arifin Jaaman, Azmi Marzuki Muda, Hairul Masrini Muhamad, Xuelei Zhang and Felicita Scapini

A. I. Mahmud ( and F. Scapini, Dipto di Biologia, Univ. degli Studi di Firenze (UNIFI), Via Madonna del Piano 6, I-50019 Sesto Fiorentino (FI), Florence, Italy.--AIM, S. A. Jaaman ( and A. M. Muda, Inst. of Oceanography and Environment (INOS), Univ. Malaysia Terengganu (UMT), Kuala Terengganu, Malaysia.--H. M. Muhamad, Key Laboratory of Underwater Acoustics & Marine Information Technology, Dept of Applied Marine Physics & Engineering College of Ocean and Earth, Xiamen Univ., Xiamen, PR China.--X. Zhang First Inst. of Oceanography, State Oceanic Administration, Qingdao, PR China.

This work is licensed under the terms of a Creative Commons Attribution 4.0 International License (CC-BY) < >. The license permits use, distribution and reproduction in any medium, provided the original work is properly cited.

doi: 10.2981/wlb.00383
Table 1. The scale applied during photographic assessment to determine 
the quality of Irrawaddy dolphins dorsal fin photos in Brunei Bay, 
based on the criteria by Gowans and Whitehead (2001).

Quality rating  Criteria

Q1              very distant, poor focus and very little flank showing,
                fin not square on
Q2              very distant photograph with little flank showing
Q3              distant photograph with little flank showing
Q4              distant photograph with most of the flank showing
Q5              close with good representation of the flank
Q6              close photograph with most of the flank showing, well
                focused and exposed image, fin square on to camera

                No. of well marked dorsal fin photos
Quality rating  imported to DARWIN 2.22

Q1              NC
Q2              NC
Q3              79
Q4              946
Q5              53
Q6              0

                No. of unmarked dorsal
Quality rating  fin photos

Q1              NC
Q2              NC
Q3              24
Q4              147
Q5              9
Q6              0

NC = not considered for the analysis.

Table 2. Conducted boat surveys for photographically identified 
Irrawaddy dolphins (adult and calf) with survey efforts from 2013-2015 
in Brunei Bay. 11.1 0% of the observed animals were not amenable to 
photo-ID. Also, 12.02% of the sampled animals' photographs (Q3-Q5) were 
not amenable to mark-recapture analysis due to lack of the presence of 
a mark in their dorsal fins.

      Boat surveys
Year  Seasons                  Observed  Captured photographically

2013  inter-monsoon (April)    10         8
      south-west monsoon       56        52
      inter-monsoon (October)   1         1
2014  north-east monsoon       47        42
      inter-monsoon (April)    16        14
      south-west monsoon        4         4
      inter-monsoon (October)  43        39
2015  north-east monsoon       21        19
      inter-monsoon (April)    39        35
      south-west monsoon       17        14
      inter-monsoon (October)  33        30

      Boat surveys             No. of dolphins
                               Captured photographically with
Year  Seasons                  well-marked dorsal fins (nicks)

2013  inter-monsoon (April)     0
      south-west monsoon       48
      inter-monsoon (October)   1
2014  north-east monsoon       38
      inter-monsoon (April)    12
      south-west monsoon        2
      inter-monsoon (October)  36
2015  north-east monsoon       17
      inter-monsoon (April)    33
      south-west monsoon       12
      inter-monsoon (October)  28

      Boat surveys             No. of dolphins
                               Newly identified
Year  Seasons                  (adult)

2013  inter-monsoon (April)     0
      south-west monsoon       12
      inter-monsoon (October)   0
2014  north-east monsoon       11
      inter-monsoon (April)     2
      south-west monsoon        1
      inter-monsoon (October)   8
2015  north-east monsoon        2
      inter-monsoon (April)     0
      south-west monsoon        0
      inter-monsoon (October)   0

                               Survey efforts (h) under different
                               Beaufort sea
      Boat surveys             states from 201 3-201 5

Year  Seasons                  B0     B1     B2     B3     Total effort

2013  inter-monsoon (April)
      south-west monsoon       23.37  44.43  8.05    8.98  84.83
      inter-monsoon (October)
2014  north-east monsoon
      inter-monsoon (April)    48.88  72.03  8.0    10.57  139.48
      south-west monsoon
      inter-monsoon (October)
2015  north-east monsoon
      inter-monsoon (April)    20.4   32.87  11.75  8.58   73.60
      south-west monsoon
      inter-monsoon (October)

Table 3. Model choice criteria and abundance estimates (N) tested in 
the mark-recapture analysis of Irrawaddy dolphins in the Brunei Bay, 
Malaysia (2013-2015), using the open-population POPAN parameterization 
in Program MARK.

Model                              [AIC.sub.C]  [DELTA][AIC.sub.C]

[phi]. [p.sub.t] [b.sub.t]           351.31            0
[[phi].sub.t] [p.sub.t] [b.sub.t]    372.34        21.02
[phi]. p. [b.sub.t].                 379.48        28.16
[[phi].sub.t] p. [b.sub.t]           394.64        44.33
[[phi].sub.t] [p.sub.t] b.           536.09       184.77
[phi]. p. b.                       17284.08     16932.77
[[phi].sub.t] p. b.                17322.49     16971.19

Model                              [AIC.sub.C] weight  Model likelihood

[phi]. [p.sub.t] [b.sub.t]         0.99997             1.00
[[phi].sub.t] [p.sub.t] [b.sub.t]  0.00003             0
[phi]. p. [b.sub.t].               0                   0
[[phi].sub.t] p. [b.sub.t]         0                   0
[[phi].sub.t] [p.sub.t] b.         0                   0
[phi]. p. b.                       0                   0
[[phi].sub.t] p. b.                0                   0

Model                              Parameters  Deviance  n

[phi]. [p.sub.t] [b.sub.t]         13             69.86  33
[[phi].sub.t] [p.sub.t] [b.sub.t]  22             66.91  32
[phi]. p. [b.sub.t].                5            116.55  33
[[phi].sub.t] p. [b.sub.t]         13            114.18  32
[[phi].sub.t] [p.sub.t] b.         19            239.07  46
[phi]. p. b.                        3          17025.43  93
[[phi].sub.t] p. b.                10          17048.26  80

                                   n with 95% confidence interval
Model                              Lower  Upper  SE for n

[phi]. [p.sub.t] [b.sub.t]         28      39     3
[[phi].sub.t] [p.sub.t] [b.sub.t]  27      37     3
[phi]. p. [b.sub.t].               28      38     3
[[phi].sub.t] p. [b.sub.t]         27      38     3
[[phi].sub.t] [p.sub.t] b.         39      55     4
[phi]. p. b.                       48     181    33
[[phi].sub.t] p. b.                79      81     1

[phi] = apparent survival; p = probability of the observation; b = 
probability that an individual from the super-population enters the 
survey area; . = constant parameter; t = time-dependent parameter; 
[AIC.sub.C] = Akaike information criterion for small sample bias; n = 
estimated population size; SE = standard error.
Model {[phi] [p.sub.t] b} did not converge.

Table 4. Model with estimated survival rate, observation rate, and 
individual entry rate from super-population in the mark-recapture 
analysis for right side fins of Irrawaddy dolphins in the Brunei Bay, 
Malaysia (2013-2015), using the open-population POPAN parameterization 
in Program MARK.

                                       [phi]-hat with 95%
                                       confidence interval
Model                       [phi]-hat  Lower  Upper  SE for [phi]-hat

[phi]. [p.sub.t] [b.sub.t]  0.98       0.89   0.99    0.01
[[phi].sub.t] [p.sub.t]
[b.sub.t]                   0.94       0.63   0.98    0.07
[phi]. p. [b.sub.t]         0.98       0.90   0.99    0.01
[[phi].sub.t] p. [b.sub.t]  0.90       0.55   0.98    0.09
[[phi].sub.t] [p.sub.t] b.  0.99       0.98   1.00   <0.001
[phi].p.b.                  0.96       0.92   0.98    0.01
[[phi].sub.t] p. b.         0.90       0.90   0.91   <0.001

                                   p-hat with 95%
                                   confidence interval
Model                       p-hat  Lower  Upper  SE for p-hat  b-hat

[phi]. [p.sub.t] [b.sub.t]  0.27   0.14   0.45    0.07          0.15
[[phi].sub.t] [p.sub.t]
[b.sub.t]                   0.28   0.15   0.46    0.08         <0.001
[phi]. p. [b.sub.t]         0.48   0.41   0.56    0.04         <0.001
[[phi].sub.t] p. [b.sub.t]  0.49   0.42   0.56    0.03         <0.001
[[phi].sub.t] [p.sub.t] b.  0.26   0.15   0.41    0.07         <0.001
[phi].p.b.                  0.61   0.58   0.64    0.02          1.00
[[phi].sub.t] p. b.         0.61   0.60   0.61   <0.001         0.99

                            b-hat with 95%
                            confidence interval
Model                       Lower   Upper   SE for b-hat

[phi]. [p.sub.t] [b.sub.t]   0.10    0.22    0.03
[[phi].sub.t] [p.sub.t]
[b.sub.t]                   <0.001  <0.001  <0.001
[phi]. p. [b.sub.t]         <0.001  <0.001  <0.001
[[phi].sub.t] p. [b.sub.t]  <0.001  <0.001  <0.001
[[phi].sub.t] [p.sub.t] b.  <0.001  <0.001  <0.001
[phi].p.b.                   0       1.00    1.00
[[phi].sub.t] p. b.         <0.001   1.00   <0.001

[phi]-hat = survival rate; p-hat = observation rate; b-hat= individual 
entry rate from super-population; SE = standard error.

Table 5. Goodness-of-fit results (program RELEASE) for the fully 
time-dependent/Cormack-Jolly-Seber model tested in mark-recapture 
analysis of individual sighting histories of Irrawaddy dolphins in the 
Brunei Bay, Malaysia (2013-2015), using open-population POPAN 
parameterization (program MARK).

Test  [chi square]  df  p-level  c

2     10.5886        8  0.2261
3      0.9804        6  0.9863
3.SR   0.6816        3  0.8775
3.SM   0.2988        3  0.9603
2+3   11.5690       14  0.6409   0.83

C = variance inflation factor; SR = survival rate; SM = seasonal 
migration; df = degrees of freedom.

Table 6. Comparison of the number of individual occurrences and 
sightings under different abiotic factors range in the Brunei Bay, 
Malaysia ([chi square]-test).

Abiotic factors       Parameters range  No. of individual occurrence

Sea depth (m)         0.1-4.99           64
                      5-9.99             74
                      10-14.99           70
                      15-19.99           53
                      25-29.99           11
                      30-34.99           10
Surface water
temperature (SWT)
([degrees]C)          27-27.99            2
                      28-28.99           12
                      29-29.99          102
                      30-30.99          102
                      31-31.99           62
                      32-32.99            2
Surface water
salinity (SWS) (ppt)  0.1-4.99           26
                      10-14.99           10
                      15-19.99           15
                      20-24.99           79
                      25-29.99          152
Turbidity (NTU)       0                  57
                      0.1-9.99          152
                      10-19.99           51
                      30-39.99           22
Nearby coastline
distance from         0.1-1.49           40
(NCDFO) (km)          1.5-2.99           98
                      3-4.49            110
                      4.5-5.99           33
                      6-7.49              1

Abiotic factors       No. of sightings

Sea depth (m)         11
Surface water
temperature (SWT)
([degrees]C)           1
Surface water
salinity (SWS) (ppt)   3
Turbidity (NTU)        8
Nearby coastline
distance from          6
(NCDFO) (km)           9

                      [chi square]-test for given probabilities
Abiotic factors       [chi square] value  df  p-value

Sea depth (m)
                      40.26               5   < 0.001
Surface water
temperature (SWT)
                      12.67               5   0.03
Surface water
salinity (SWS) (ppt)
                       8.80               4   0.06
Turbidity (NTU)
                       0.66               3   0.88
Nearby coastline
distance from
(NCDFO) (km)
                      27.23               4   < 0.001

For 20-24.99 m depth, 5-9.99 ppt salinity and 20-29.99 NTU turbidity 
ranges, we did not encounter any dolphins.

Table 7. Estimated population sizes of Irrawaddy dolphins in different 
areas. CI = confidence interval; CV = coefficient of variance; SE = 
standard error; NA= not available.

Countries                   Places

Australia                   Bowling Green Bay
Australia                   Cleveland Bay, northeast
Bangladesh                  Coastal water bodies of
Bangladesh                  Sundarbans mangrove forest
Cambodia                    Mekong River
Cambodia, LAO PDR, Vietnam  Mekong River
India                       Chilika Lake
Indonesia                   Mahakam River
Indonesia                   Banten Bay
Indonesia                   Balikpapan Bay
Malaysia                    Kuching Bay, Sarawak
Malaysia                    Cowie Bay, Sabah
Malaysia                    Brunei Bay, Malaysia
Myanmar                     Ayeyarwady River
Philippines                 Malampaya Sound
Thailand                    Songkhla Lake
Thailand                    Eastern Gulf Coast of Thailand
Thailand                    Trait Bay, Trait province
Thailand                    Bangpakong Estuary

Countries                   Estimated population size

Australia                      3
Australia                     67 (95 % CI = 51-88)
Bangladesh                  5383 (CV = 40 %)
Bangladesh                   451 (CV = 9.6 %)
Cambodia                      40
Cambodia, LAO PDR, Vietnam   125 (95% CI = 114-152)
India                        109 (CV = 7 %)
Indonesia                     70 (CV = 10 %)
Indonesia                     10-15
Indonesia                     67 (59-74)
Malaysia                     149 (CV = 28 %)
Malaysia                      31 (SE = 1.8508)
Malaysia                      41 (95 % CI = 36-49)
Myanmar                       58-72
Philippines                   77 (CV = 27.4 %)
Thailand                       1-15
Thailand                     423
Thailand                     171 (SE = 73.18)
Thailand                      20-5

Countries                   Size of the study area ([km.sup.2])

Australia                        NA
Australia                        NA
Bangladesh                    47150
Bangladesh                     1756
Cambodia                         NA
Cambodia, LAO PDR, Vietnam       NA
India                          1165
Indonesia                        NA
Indonesia                       150
Indonesia                       120
Malaysia                    ca. 520
Malaysia                        500
Malaysia                    ca.1690
Myanmar                          NA
Philippines                    2001.15
Thailand                       1040
Thailand                         NA
Thailand                         NA
Thailand                        135

Countries                   Reference(s)

Australia                   Parra and Corkeron 2001
Australia                   Parra 2005
Bangladesh                  Smith et al. 2005
Bangladesh                  Smith et al. 2006
Cambodia                    Baird and Beasley 2005
Cambodia, LAO PDR, Vietnam  Beasley et al. 2007
India                       Sutaria and Marsh 2011
Indonesia                   Kreb et al. 2007
Indonesia                   Khalifa et al. 2014
Indonesia                   Kreb 2008
Malaysia                    Minton et al. 2013
Malaysia                    Woan et al. 2013
Malaysia                    this study
Myanmar                     Smith et al. 2007
Philippines                 Smith et al. 2004
Thailand                    Kittiwattanawong et al. 2007
Thailand                    Hines et al. 2015
Thailand                    Junchompoo et al. 2014
Thailand                    Tongnunui et al. 2011

Table 8. Estimated group sizes of Irrawaddy dolphins in different 
areas. SE = standard error; SD = standard deviation.

Countries   Places                      Mean group size

Australia   Cleveland Bay               5.35 (SE= 0.35)
Bangladesh  Sundarbans Mangrove forest  2.30 (SD = 1.36)
Indonesia   Mahakam River               4.4 (SD = 2.2)
Indonesia   Malampaya Sound             5.26 (SE = 1.06)
Malaysia    Brunei Bay                  6 (SE = 0.66)

Countries   Group size range  Reference(s)

Australia   1-15              Parra 2005
Bangladesh  1-6               Smith et al. 2006
Indonesia   1-10              Kreb and Budiono 2005a
Indonesia   1-15              Dolar et al. 1997
Malaysia    1-18              this study

Table 9. Marine mammals observed in the Brunei Bay, Malaysia, during 
the survey period (2013-2015) with IUCN Red List status. No. = number; 
IUCN = International Union for Conservation of Nature. Table modified 
from Mahmud et al. 2018.

Serial no.  English name

01          Indo-Pacific bottlenose dolphin
02          Indo-Pacific finless porpoise
03          Indo-Pacific humpbacked dolphin
04          Irrawaddy dolphin
05          Killer whale
06          Dugong

Serial no.  Scientific name

01          Tursiops aduncus (Ehrenberg, 1833)
02          Neophocaena phocaenoides (G. Cuvier, 1829)
03          Sousa chinensis (Osbeck, 1765)
04          Orcaella brevirostris (Owen in Gray, 1866)
05          Orcinus orca (Linnaeus, 1758)
06          Dugong dugon (Muller, 1776)

Serial no.  IUCN (2015) Red List Status

01          data deficient
02          Vulnerable
03          Near Threatened
04          Vulnerable
05          data deficient
06          Vulnerable
COPYRIGHT 2018 Nordic Board for Wildlife Research
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2019 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Mahmud, Anisul Islam; Jaaman, Saifullah Arifin; Muda, Azmi Marzuki; Muhamad, Hairul Masrini; Zhang,
Publication:Wildlife Biology
Article Type:Report
Geographic Code:9MALA
Date:Jan 1, 2018
Previous Article:Seasonal dynamics of forage for red deer in temperate forests: importance of the habitat properties, stand development stage and overstorey dynamics.
Next Article:Population monitoring and modelling of yellow-shouldered parrot on Bonaire, Caribbean Netherlands.

Terms of use | Privacy policy | Copyright © 2022 Farlex, Inc. | Feedback | For webmasters |