Population parameters and heterogeneity in survival rates of Indo-Pacific humpback dolphins in a heavily urbanized coastal region of southeast China: implications for conservation

1 Introduction

Indo-Pacific humpback dolphins (Sousa chinensis, hereafter referred to as humpback dolphins) are obligatory inshore delphinids that inhabit estuarine and coastal waters off southeast Asia and western Pacific (Jefferson and Smith, 2016). With few exceptions, most extant populations of this species number only several tens or at most a few hundred of individuals (Jefferson and Smith, 2016), likely due to their dependence on shallow-water inshore habitats which are generally patchy and limited in size (Gowans et al., 2007), as seen in other species of the genus Sousa (e.g., Karczmarski et al., 2000; Parra et al., 2006). More critically, their habitats are often heavily impacted by adverse effects of various human activities, much of which concentrate in coastal regions (Lin et al., 2016; Karczmarski et al., 2017a; Karczmarski et al., 2017b; Huang et al., 2020; Lu et al., 2023). Wherever studied across the known species range, populations have been declining in numbers in recent years and decades (Reeves et al., 2008; Jefferson and Smith, 2016; Lin et al., 2022), while conservation management measures are either absent or ineffective (e.g., Karczmarski et al., 2016; Karczmarski et al., 2017a).

In waters of the Pearl River Delta (PRD) region on the southeast coast of China, humpback dolphins are found along >1000 km of coastline from Hong Kong in the east to Xiachuan Island and Yangjiang city in the west (Chen et al., 2011; Chan, 2019; Li et al., 2019) and are thought to represent the world’s largest population of this species. While human impacts on the coastal ecology of PRD go back thousands of years (Zong et al., 2009) and have been affecting humpback dolphins for centuries (Lin et al., 2016), the severity of anthropogenic impacts increased drastically since the onset of the industrial boom in mainland China in the 1990s (Zhu et al., 2013; Wu et al., 2016; Zhang et al., 2016). A greatly accelerated process of urbanization and industrialization, along with massive coastal alterations and land reclamation have degraded much of the PRD coastal habitats, especially at the eastern reaches of the region in Lingding Bay¹ (Wu et al., 2016) which represents the main body of the Pearl River estuarine system. In the absence of effective environmental management strategy, major anthropogenic perturbations over the past decades have turned Lingding Bay (LB) into one of the most polluted and ecologically disturbed estuaries in the world (Li and Damen, 2010; Wang et al., 2016; Wu et al., 2016; Wang and Rainbow, 2020; Li et al., 2022). While much of the environmental health of LB has been severely compromised (Wang and Rainbow, 2020), any attempts at protection of its cetacean inhabitants have been further hampered by the lack of basic knowledge of their population demographic processes and structure (Karczmarski et al., 2016).

Much of the current knowledge of humpback dolphins in the estuarine system of the Pearl River comes from studies undertaken in late 1990s and early 2000s (e.g., Jefferson, 2000; Jefferson et al., 2006; Jefferson et al., 2012), which were conducted primarily in Hong Kong (HK) waters with only occasional surveys deeper into Lingding Bay. These early works and subsequent monitoring program in Hong Kong (e.g., Hung, 2015) along with preliminary line-transects across the PRD (Chen et al., 2010) have been used for over two decades as a primary reference to infer the status and trend of the entire population (e.g., Huang et al., 2012; Miller, 2016; Jefferson, 2018). This should be a cause for concern, given that much of the quantitatively replicable data were from Hong Kong (HK) waters that represent only a small fraction of Lingding Bay and even smaller fraction in the spatial scale of the PRD. Admittedly, in the absence of any other data, these were the only means available at the time that could facilitate any form of potential management and conservation effort. With time, however, it became increasingly evident that further studies of population processes of the PRD dolphins were much needed and at spatio-analytical scales that could provide a considerably greater data resolution than the early work could offer (Karczmarski et al., 2016). Some such studies have recently been completed and indicate a distinct socio-spatial and demographic structure of the dolphin population across the PRD (Chan, 2019; Guo et al., 2020; Chan et al., 2022a), and further fine-scale socio-spatial structure within LB where multiple clusters interact socially but remain considerably dissimilar spatially (Karczmarski et al. unpublished²). Various environmental factors also vary spatially in LB. Numerous man-made stressors differ in their type (e.g., habitat loss, pollutant load, maritime traffic, underwater noise) and intensity (e.g., extent of habitat destruction, intensity of fisheries pressure) across the anthropogenic seascape of LB (Chen et al., 2013; Wang et al., 2016; Wu et al., 2016; Yu et al., 2016; Gui et al., 2017; Pine et al., 2017; Gui et al., 2019; Guo et al., 2022; Li et al., 2022). Correspondingly, there is a considerable spatial heterogeneity in dolphin prey resources, an apparent result of heavy fishing pressure in the region (e.g., Duan et al., 2009; Lin et al., 2019; Zhou et al., 2019). Prey quantity/availability and size/quality appear to vary spatially across both PRD and LB, likely affecting the energetic costs of foraging (Zhang et al., 2019; Ho, 2021; Lin et al., 2021) which may have implications on individual fitness. Spatial heterogeneity in both resource availability and anthropogenic pressure have been attributed to spatial heterogeneity in demographic process in other taxa elsewhere (e.g., Reid et al., 2008; Harris et al., 2021; Ritterson et al., 2021). Consequently, management plans and conservation priorities in the PRD and/or LB have to be perceptively formulated with careful considerations of population structure and habitat utilization of the dolphins, as well as the quality of key dolphin habitats and levels of anthropogenic disturbance (Karczmarski et al., 2016; Karczmarski et al., 2017a). Otherwise, conservation resources and effort may be misplaced to areas of low conservation value (Freitag and Jaarsveld, 1997; Karczmarski et al., 2016).

Early studies of humpback dolphins in China applied primarily line-transect techniques to collect information such as relative abundance and distribution (e.g., Jefferson, 2000; Hung and Jefferson, 2004; Chen et al., 2008; Chen et al., 2010). Recent studies have increasingly used mark-recapture techniques for more advanced demographic analyses (e.g., Chan and Karczmarski, 2017; Chen et al., 2018; Zeng et al., 2020; Chan et al., 2022a; Lin et al., 2022; Chan and Karczmarski, 2024). Mark-recapture population modeling techniques not only deliver robust estimates of demographic parameters with considerable accuracy (e.g., Hastings et al., 2011; Pace et al., 2017; Kendall et al., 2019; Oosthuizen et al., 2019a; Schleimer et al., 2019; Bright Ross et al., 2022), but provide also powerful tools to identify the underlying driver(s) of ecological and demographic processes (e.g., Horswill et al., 2014; Oosthuizen et al., 2019b; Coxon et al., 2022; Gabriele et al., 2022; Marneweck et al., 2022; Jordaan et al., 2023) and offer a flexible framework for population monitoring that is both effective and reliable, and applicable across a wide range of taxa and management needs (e.g., Atkinson et al., 2021; Chan et al., 2022b; Goldenberg et al., 2022; Szott et al., 2022; Verborgh et al., 2022).

In the study presented here, we built a coordinated photo-identification (photo-ID) database by cross-matching multi-year data collected across Lingding Bay by research teams based in Hong Kong and mainland China. With the application of mark-recapture modelling techniques, we quantify vital population parameters of the LB humpback dolphins, and investigate the effects of age-class and individual spatial preferences on the probabilities of survival and recruitment. With the latest population size estimates delivered in this study, we construct a recent trajectory of the population by comparing previous abundance estimates (from late 1990s) with our current findings. Our study quantifies spatial variations and temporal changes of demographic parameters of humpback dolphins inhabiting LB waters, which improves our understanding of their population processes. More importantly, it provides valuable baseline information for delineating conservation management priorities in the large and complex estuarine system of Lingding Bay.

2 Materials and methods

2.1 Study area

The second largest river in China by freshwater discharge, the Pearl River, with an annual runoff of 320 billion m³ (Lu et al., 2009) enters the South China Sea through eight major distributaries in a region known as the Pearl River Delta (PRD). Four of the distributaries enter Lingding Bay (LB) and form the main body of the Pearl River estuarine system, in the east of the PRD (Figure 1). The metropolitan area that surrounds LB (including Shenzhen, Dongguan, Guangzhou, Zhongshan and Zhuhai in mainland China, and the Special Administrative Regions of Hong Kong and Macao) comprises 97.5% of the human population of the PRD economic zone, making LB one of the most urbanized and industrialized estuaries in the world (see Karczmarski et al., 2016 for more details). The present study focuses on LB, including the administrative sectors of Guangdong Province (mainland China) and Hong Kong (Figure 1). As LB represents the main body of the extensive estuarine system of the Pearl River, it is also referred to as the “Pearl River Estuary (PRE)”. It should not be mistaken, however, with what some earlier studies (e.g. Jefferson, 2000; Jefferson and Hung, 2004) referred to as PRE, as much of the earlier work was based primarily (though not exclusively) in Hong Kong waters, a small section of a much larger LB within still larger and physiographically complex PRD coastal system (e.g., Karczmarski et al., 2016; Lin et al., 2021; Chan et al., 2022a).

Figure 1 The study area (solid red line) in the Lingding Bay (LB), at the eastern reaches of the Pearl River Delta (PRD) region, southeast China. The thick black dashed line indicates the administrative border between Guangdong Province and Hong Kong (HK). The study area in Guangdong waters is further divided into inner Lingding Bay (iLB) and outer Lingding Bay (oLB) along the East-West shipping lane across the bay (dashed red line); see also ‘Data structure’. The Hong Kong-Zhuhai-Macao Bridge (HZMB) displayed in this figure represents the construction status at the time of this study, when the eastern section of the bridge in Hong Kong waters was not yet finished (grey dashed line).

2.2 Data collection and photographic identification

From January 2011 to December 2015, when sea conditions and weather allowed (sea state ≤ 3 in Beaufort scale and visibility ≥ 1 km), boat-based photo-identification surveys were carried out concurrently in Guangdong and Hong Kong waters, following internationally recognized protocol (e.g., Hammond et al., 1990; Karczmarski et al., 2022a; Karczmarski et al., 2022b). Once a group of dolphins was encountered, the research vessels approached the animals at speed ≤6 km/h. Photographs of dolphin dorsal fins and upper bodies were taken with at least two sets of high-speed image-stabilized digital SLR cameras (Canon EOS 1D III/IV/X) each equipped with 100-400mm f/4.5-5.6 image-stabilized zoom lens. A conscious effort was made to photo-capture both right and left sides of all group members irrespective of their distinctiveness, age-class or behavior (Chan and Karczmarski, 2017; Chan and Karczmarski, 2024). The geographic location of each encounter was recorded using global positioning system (GPS). Each group was followed for at least 10 minutes or until all group members were photographed.

To ensure equal “catchability” of all individuals, only the features on the dorsal fin and dorsal ridge that were well-visible above the water (such as the shape of dorsal fin, pattern of notches, pigmentation of the dorsal ridge and fin, and prominent injury marks) were used for individual identification (Chan and Karczmarski, 2017; Lin et al., 2018); and the image quality (Q) and individual distinctiveness (D) were rated independently (Chan et al., 2022a). The image quality of ID-photographs was rigorously assessed with a score ranging from Q1 to Q100 (adopted from Karczmarski et al., 2005) based on multiple criteria, such as how well exposed and in-focus the image was, the contrast against the background, the angle of the animal to the camera, and the proportion of dolphin body photographed above the water. The distinctiveness of individuals was rated from D0 to D5, with D0 representing individuals with no identifiable features and D5 representing highly distinctive individuals (Friday et al., 2000). A consistent application of our rigorous quality criteria along with the use of long-term reliable individually distinctive features (e.g., pigmentation of the dorsal fin and ridge, pattern of notches) and considerable survey effort (frequent surveys, whenever the weather/sea conditions allowed) facilitated an effective monitoring of any potential individual mark-change (after Karczmarski et al., 2005).

All photographs of individuals with distinctiveness D≥1 and image quality Q≥60 were matched and catalogued using the photo-ID data management program DISCOVERY (Gailey and Karczmarski, 2012; Chan et al., 2022c). To minimize any potential bias from mis-identification due to insufficient image quality and/or individual distinctiveness, only high-quality ID-photographs (Q≥70) of highly distinctive individuals (D≥3) were used for all subsequent analyses. Calves were excluded from the analytical dataset as the majority of them lack reliably recognizable features (Chan and Karczmarski, 2017; Lin et al., 2018; Guo et al., 2020; Chan and Karczmarski, 2024), but were retained in the mark ID-ratio calculation to give inference of the entire population (see section 2.6 for further details). Juveniles and adults were categorized as two age-class groups as they generally have different demographic characteristics (Currey et al., 2009; Manlik et al., 2016); see Chan and Karczmarski (2017) for detailed description/definition of calves, juveniles and adults.

2.3 Data structure

The individual sighting histories were pooled by calendar year to provide annual estimates of population parameters. While mark-recapture models assume equal demographic parameters (such as apparent survival rates and capture/recapture probabilities) for all marked individuals, this assumption is often violated in studies of marine mammals, especially in the presence of population substructure and individual heterogeneity in habitat use (Peel et al., 2016). A recent study of socio-spatial dynamics of humpback dolphins in LB (Karczmarski et al. unpublished²) identified three ‘core areas’ within the present study area (each representing a relatively un-urbanized coastal habitat), namely waters surrounding islands in the inner Lingding Bay (iLB; around Neilingding and Qi’Ao Islands) and the outer Lingding Bay (oLB; around Green and Triangle Islands), and waters off the west and southwest coast of Lantau Island in Hong Kong (HK). Although the dolphins using these three core habitats are inter-connected socially, their individual range appears restricted, with considerable long-term fidelity to their respective core areas (Karczmarski et al. unpublished²), similarly as described recently for humpback dolphins in western reaches of the PRD (Chan, 2019; Chan et al., 2022a). The variation in individual spatial preference (area/habitat use), together with the differences in the intensity of anthropogenic stressors across the region (Karczmarski et al., 2016), may lead to heterogeneity in both survival and recapture estimates, which had to be accounted for in our analytical design (see Discussion).

To examine the potential heterogeneity in population parameters, as well as the variability that might stem from the differences in survey effort, the study area was divided into three sub-regions (Figure 1): (i) Hong Kong (HK) demarcated by the HK-Guangdong administrative border, (ii) the inner LB (iLB) and (iii) outer LB (oLB), north (iLB) and south (oLB) of the East-West shipping lane between Zhuhai/Macao and Hong Kong. The East-West shipping lane demarcates also the transitional area of environmental features such as salinity (Wang and Lin, 2006), prey resource (Pine et al., 2018; Lin et al., 2019) and much of human activities (Karczmarski et al., 2016). The index of individual spatial preference of a sub-region (RP_ij) was calculated for all identified dolphins as follows:

where r_ij and r_i.total represent the number of sightings of individual i in sub-region j (either HK, iLB, or oLB) and in the whole study area, respectively; while D_j and D_total represent the survey effort (days; Table 1) in area j and in the whole study area, respectively.

Table 1 Survey effort (days) spent in the three sub-regions of Lingding Bay (LB).

Subsequently, agglomerative hierarchical clustering (AHC) was performed based on the index of individual spatial preference RP_ij using four clustering methods (‘complete’, ‘average’, ‘single’, and ‘ward.D’) in the R package “hclust” (Scrucca et al., 2016). The best clustering method resulting in the highest agglomerative coefficient value was selected, and the number of clusters (K) was then identified using the silhouette algorithm by calculating the optimum of the total sum of square of agglomerative coefficient value (the squared difference between each individual from the overall mean of clusters) with K ranging from 1 to 10. Thereafter, the mark-recapture dataset was structured with individuals grouped according to the resultant clusters of spatial (sub-region) preference, which allows further examinations of the potential heterogeneity. An equilateral ternary plot was constructed to illustrate how clustered (or scattered) the dolphins are in their regional preferences. As the RP indices were corrected for the survey effort (see Equation 1), the values have to be standardized to have them displayed in the same scale on all three axes of the equilateral triangular plot. The standardized RP_ij indices used for the ternary plot were therefore rescaled across RP_HK, RP_iLB and RP_oLB as follows:

2.4 Goodness-of-fit (GOF) tests

Prior to mark-recapture modelling, GOF tests were performed to estimate the variance inflation factors ($\widehat{c}$, the level of overdispersion of data) of various saturated Cormack-Jolly-Seber (CJS) models (Cormack, 1964; Jolly, 1965; Seber, 1965) using the bootstrap and median- $\widehat{c}$methods (Anderson et al., 1994). We started with the standard time-dependent model with the time-varying effect (‘t’) on both apparent survival rates (φ) and recapture probabilities (p), and added the age-class effect (‘Age’) and/or the individual regional preference effect (‘3site’), with the exception of age-class effect on recapture probabilities, as juveniles and adults had the same chance of being photographed under the rigorous field protocol. The value of $\widehat{c}>1$ indicates overdispersion of data while $\widehat{c}<3$ indicates acceptable fitness to the data, and the model with the lowest $\widehat{c}$ was used as the starting model.

2.5 Modeling apparent survival rates and recapture probabilities

Mark-recapture analyses were performed using program MARK (White and Burnham, 1999). The most parsimonious model, which shows best fit to the data, was selected based on the Akaike’s Information Criterion (AIC_c) or Quasi-AIC_c (QAIC_c, after the correction of $\widehat{c}$) (Burnham and Anderson, 2002). If multiple models were similarly supported (ΔQAIC_c< 2), weighted averages of parameters were calculated according to the AIC_c/QAIC_c weight (w_i) of models (Buckland et al., 1997).

To estimate the apparent survival rates (φ) and recapture probabilities (p) of humpback dolphins, CJS models were constructed with various effects applied on both parameters. Besides the time-varying effect (‘t’), constant effect (‘.’) and age-class effect (‘Age’), regional preference effects were included as individual covariates to test the assumptions that: (1) individuals grouped in the three clusters differ in their apparent survival/recapture probabilities (‘3site’); and/or (2) dolphins seen more frequently in Guangdong waters (iLB + oLB; ‘2site_a’) or (3) in outer reaches of the LB (oLB + HK; ‘2site_b’) exhibit similar apparent survival/recapture probabilities but differ from other individuals. Survey effort effect (the total number of survey days per year; ‘effort’) was incorporated as an additional covariate for recapture probabilities.

2.6 Modeling super-population size and entry probabilities

POPAN formulation of Jolly-Seber model (Schwartz and Arnason, 1996; hereafter referred to as POPAN model) was used to estimate the super-population size (N), which represents the number of animals that used the study area as part of their range at any time during the study period, and the probability of entry (PENT), which represents the probability of an individual entering the super-population during the study period as either birth, immigration, or new identification. In the POPAN analysis, we applied the same effects on apparent survival rates and capture probabilities as in the most parsimonious CJS model, and further modeled the effects of time, regional preference and survey effort on the probabilities of entry. As GOF tests are not available for POPAN models, the value of $\widehat{c}$ based on that of CJS model was adopted to account for the over-dispersion of data.

Since mark-recapture analyses can only estimate the number of marked individuals (N), the total super-population size (including both marked and unmarked individuals, ${\widehat{\text{N}}}_{\text{T}}$) was projected by the mark ID-ratio $(\widehat{\theta })$ following: ${\widehat{\text{N}}}_{\text{T}}=\widehat{\text{N}}/\widehat{\theta }$, where the mark ID-ratio (the proportion of marked individuals) was estimated from a pool of 100,000 randomly selected images. The variance of ${\widehat{\text{N}}}_{\text{T}}$ was calculated following Urian et al. (2015) as given in Equation 3:

The upper and lower bounds 95% confidence intervals of ${\widehat{\text{N}}}_{\text{T}}$ were estimated with ${\widehat{\text{N}}}_{\text{T}}^{\text{lower}}={\widehat{\text{N}}}_{\text{T}}/C$ and ${\widehat{\text{N}}}_{\text{T}}^{\text{upper}}={\widehat{\text{N}}}_{\text{T}}\times C$, where C was calculated as given in Equation 4:

following Burnham et al. (1987).

3 Results

3.1 Data summary

During the 5-year study period, photo-ID data across LB were collected in 564 survey days. A total of 1287 groups of humpback dolphins were encountered, resulting with 7269 cumulative sighting records of various individuals and over 400,000 ID-images of sufficient quality (Q ≥ 60). A total of 809 highly marked (D ≥ 3) individuals were catalogued using high-quality (Q ≥ 70) images (Figure 2A), of which 107, 199 and 187 dolphins were photographed only in HK, iLB and oLB, respectively; while 236 individuals were seen in any two of the three sub-regions and 80 individuals were found in all three sub-regions across LB (Figure 2B). The mean number of sightings was 11.1 while the most frequently seen individual was photographed 73 times; and 655 individuals (81%) were photographed more than once. The rate of identifying new individuals was similar in all three sub-regions (HK, iLB and oLB regional discovery curves) and remained so throughout the study (Figure 2A). The overall cumulative number of newly identified dolphins (cumulated discovery curve) climbed sharply in the first year and continued to increase, albeit at a considerably lower rate through the end of the study period (Figure 2A).

Figure 2 (A) The overall cumulative discovery curve (dark grey line) and three regional discovery curves (green, red, and blue lines) of well-marked (D ≥ 3) Indo-Pacific humpback dolphins photo-identified with high-quality (Q ≥ 70) images in Lingding Bay (LB) during five years from 2011 to 2015. The small number of newly identified dolphins (black bars) and the high number of recaptured individuals per month (light grey bars) indicate good coverage of the population in the present study. (B) Graphic summary of sightings of the 809 well-marked (D ≥ 3) humpback dolphins photo-captured in LB between 2011 and 2015, of which 107, 199 and 187 were photographed only in HK (green circle), iLB (red circle) and oLB (blue circle), respectively; while 80 individuals were seen in all three sub-regions across LB.

3.2 Agglomerative hierarchical clustering analysis

All four AHC methods estimated agglomerative coefficient close to one (‘complete’: 0.997; ‘average’: 0.994; ‘single’: 0.998; ‘ward.D’: 0.999), providing a strong indication of a spatial structure among dolphins inhabiting Lingding Bay. Using the best performing AHC method (‘ward.D’), three clusters were identified (Figure 3) corresponding to individual dolphin preferences (regional preference, RP) of specific sub-regions (Table 2). Although a small number of individuals displayed similar preference to more than one sub-region (at the center of Figure 3C), our clustering approach have structured the mark-recapture dataset sufficiently well to account for the over-dispersion of data (see further details in Section 3.3).

Figure 3 Agglomerative hierarchical clustering (AHC) analysis of spatial (sub-region) preferences of Indo-Pacific humpback dolphins in Lingding Bay (LB). The dendrogram (A) indicates three clusters of individuals identified using the ‘ward.D’ clustering method (the color of clusters corresponds to colors in Figures 2A, B, and panel (C)), where the optimal number of clusters (B) is determined using the silhouette algorithm. The ternary plot (C) displays the distribution of the regional preference of dolphins to the three sub-regions, inner LB (iLB), outer LB (oLB) and Hong Kong (HK), where the coloured dots indicate individuals of the three clusters. The numeric values of individual regional preference (RP) indices shown in panel (C) were standardised according to Equation 2 (see Methods section 2.3) in order to have them displayed in the same scale on all three axes of the equilateral trianglar plot.

Table 2 Summary statistics of the regional preference (RP) indices.

3.3 GOF test

When only time-varying effect (‘t’) was considered, the fit of the dataset was poor $(\widehat{c}>3)$ with the highest level of over-dispersion of data among the models tested (Table 3). A clear asymmetry in the deviance residuals (Supplementary Figure S1) suggested that the lack of fitness was largely due to unaddressed data structure. The fit of data was improved when age-class effect (‘Age’) was incorporated to apparent survival rate parameters, and even more so when individuals were clustered into three groups according to their regional preferences (‘3site’) ($\widehat{c}<3$; Table 3). Further improvement was only meagre when both effects of age-class and regional preference were incorporated. However, in order to include both effects to further examine the related effects in the model selection process, the model with all three effects (time-varying, age-class, and regional preference) was selected as the starting model (Table 3; see also Sections 3.4 and 3.5). To take a conservative approach, the higher $\widehat{c}$ value of 3.30 using the bootstrap GOF method was adopted (even though $\widehat{c}$ by the median- $\widehat{c}$ method was well below 3) to correct for QAIC_c during the model selection.

Table 3 Results of the goodness-of-fit (GOF) test of the 5-year (2011-2015) mark-recapture dataset of humpback dolphins in Lingding Bay; φ: apparent survival rate, p: recapture probability.

3.4 Apparent survival rates and recapture probabilities

The most parsimonious (best-fit) CJS candidate model presented the QAIC_c weight evidence ratio (w₁/w₂) of nearly 15 times more likely than the next model (ΔQAIC_c > 5, Table 4) and was therefore used for the estimation of parameters. The result of model selection suggested that the recapture probabilities varied over time (‘t’), but with little impacts of survey effort (‘effort’, model #3). Adult dolphins presented higher apparent survival rates than juveniles (‘Age’), while both survival and recapture estimates showed spatial differences (‘2site_a’, the same for iLB and oLB, but different from HK; Table 4). Dolphins that were seen more frequently in Hong Kong waters exhibited higher survival rates (HK, adults: φ = 0.980, SE = 0.005, 95% CI = 0.968-0.987; juveniles: φ = 0.926, SE = 0.018, 95% CI = 0.883-0.954) than those seen primarily in Guangdong waters (iLB and oLB, adults: φ = 0.919, SE = 0.012, 95% CI = 0.893-0.939; juveniles: φ = 0.748, SE = 0.039, 95% CI = 0.666-0.816) (Figure 4A). The difference in survival estimates between adults and juveniles was also more prominent for the iLB and oLB clusters (individuals with preference to Guangdong waters) than the HK cluster (Figure 4A). When averaged across sub-regions based on the covariate of individual regional preference, the apparent survival rates were estimated at 0.943 (SE = 0.008, 95% CI = 0.929-0.958) and 0.815 (SE = 0.025, 95% CI = 0.766-0.866) for adults and juveniles, respectively. While the overall recapture probabilities for all individuals remained high throughout the study period (p = 0.762-0.861 at all times), the estimates for the HK cluster (p = 0.844-0.915) were noticeably higher than for iLB and oLB clusters (p = 0.629-0.830). Moreover, the recapture rates in HK were relatively stable, while the estimates in iLB and oLB decreased across years (Figure 4B).

Table 4 Selection of 11 Cormack-Jolly-Seber (CJS) candidate models in estimating the apparent survival rates (φ) and recapture probabilities (p) of Indo-Pacific humpback dolphins in Lingding Bay during the period 2011-2015.

Figure 4 Estimates of (A) apparent survival rates and (B) recapture probabilities from the most parsimonious Cormack-Jolly-Seber (CJS) model. The apparent survival rates (A) show differences between age-classes (adults vs. juveniles) as well as spatially across the LB (iLB+oLB vs. HK); while the recapture rates (B) show both temporal (across years) and spatial (iLB+oLB vs. HK) differences. The overall apparent survival rates and recapture probabilities represent the averaged estimates across the sub-regions of LB. Error bars indicate standard errors (SE) of estimates.

3.5 Population size and entry probabilities

Adopting the effects on apparent survival rates and recapture probabilities of the most parsimonious CJS model (see Section 3.4 and model #1 in Table 4), the two best-fit POPAN candidate models accounted for 90% of model support, which suggested a mixture of age-class (‘Age’) and effort (‘effort’) effects in the estimation of entry probabilities (PENT) (Table 5). As both models (models #1 & #2) received comparable support (ΔQAIC_c< 2), model-averaging was performed for the estimates of both entry probabilities and super-population size (N). Entry probabilities for both juveniles (PENT_juvenile = 0.076-0.147) and adults (PENT_adult = 0.029-0.057) decreased across years (Table 6) and were used to derive the annual entrance of marked individuals (N_E), which ranged from 10 to 20 juveniles and 21-40 adults per year across the study period (Table 6). Super-population size estimates (N) suggested 142 (CV = 0.064, 95% CI = 124-160) and 713 (CV = 0.015, 95% CI = 692-734) highly marked (D ≥ 3) juveniles and adults, respectively. After correction for the mark-ID ratio (θ = 92.0%, SE = 0.27%), the total super-population size (N_T) was estimated at 928 (CV = 0.061, 95% CI = 823-1046) humpback dolphins, including 154 juveniles and 774 adults. The annual super-population size estimates (after mark-ID ratio correction) fluctuated between 708 (CV = 0.127, 95% CI = 553-907) and 750 (CV = 0.077, 95% CI = 645-871) dolphins (Table 7).

Table 5 Selection of 13 POPAN candidate models considered in estimating the apparent survival rates (φ) and capture probabilities (p), entry probabilities (PENT) and super-population size (N) of Indo-Pacific humpback dolphins in Lingding Bay during the period 2011-2015.

Table 6 Model-averaging estimates of entry probabilities (PENT) and numbers of annual entrance (N_E) from two best-fit POPAN candidate models.

Table 7 Model-averaging estimates of annual super-population size (N_Ti; after correction for mark-ID ratio) from two best-fit POPAN candidate models.

4 Discussion

With the application of photo-identification mark-recapture techniques, we quantified population parameters of Indo-Pacific humpback dolphins inhabiting Lingding Bay (LB) – the main body of the Pearl River Estuary (PRE) and the most urbanized and industrialized part of the greater Pearl River Delta (PRD) region. The overall super-population size was estimated at 928 individuals, with generally low survival rates. The annual estimates – likely indicative of the latest number of humpback dolphins in LB (see further discussion) – fluctuated between 708 and 750 individuals. Both the overall and annual estimates generated by our study are considerably lower than previously published (Jefferson, 2000; Chen et al., 2010; Hung, 2015) abundance estimates (see 4.5 ‘Demographic trajectory’). The dolphins exhibited heterogeneous apparent survival rates across the complex anthropogenic seascape of LB, with the rates notably lower outside Hong Kong territorial waters.

4.1 The dataset

Heterogeneity in data structure (e.g., heterogeneity in capture probabilities), if unaccounted for, may prevent mark-recapture models from accurately fitting the sighting history data and it may introduce biases to the estimates of population parameters (Sandercock, 2006). Therefore, the robustness of models and the potential violation of the underlying assumptions needs to be carefully assessed. Given the relatively long lifespan of humpback dolphins (~38 years; Jefferson et al., 2012) as compared to the duration of our study, it was to be expected that entries of marked individuals, both adults and juveniles, declined gradually over time, dropping by half in 2014 (see Section 3.5), which corresponds to the semi-asymptotic discovery curve (Figure 2). Given the considerable population size with limited albeit continuous recruitment, the discovery curve of dolphins in LB is not expected to reach a plateau even with additional years of data. As such, the annual super-population estimates (Table 7), especially in the later years of this study that covered most individuals with a cumulated survey effort, are likely to accurately reflect the dolphin numbers in waters of LB (for comparison, see Zeng et al., 2020). Even though some LB humpback dolphins were seen moving occasionally to and from neighboring Modao Estuary, further to the west of the present study area (Chan, 2019; Chan et al., 2022a), the transition rates are exceptionally low (Chan, 2019; see also further) and potential implications on recapture probabilities, if any, can be considered negligible (Chan, 2019; Chan et al., 2022a).

We found no indication that dolphin movement may cause biases in the estimation of apparent survival rates. Firstly, movement of humpback dolphins outside LB across the PRD is limited, with less than 2% of individuals seen occasionally in waters off Macao but no further west (Chan, 2019), suggesting that downward-bias associated with permanent emigration is negligible. Secondly, dolphin movement within LB is not likely to impact estimates of survival rates, as our field surveys covered almost the entire LB across all three sub-regions. For example, prior to this study, Chan and Karczmarski (2017) modeled the demographic parameters of humpback dolphins seen in Hong Kong waters and their estimate of adult apparent survival rates (φ = 0.980) is highly consistent with the apparent survival rates of adults from the HK cluster in the present study. While the work by Chan and Karczmarski (2017) was spatially restricted to Hong Kong territorial waters only, in the present study (with a substantially larger survey area) the analytical dataset included sightings of individuals from the HK cluster in other sub-regions of LB as well, further lowering the likelihood of biases from movement and/or emigration.

Finally, recapture probabilities remained high throughout the study period and for all individuals (ranging from 0.762 to 0.861; Figure 4), minimizing the potential bias in the estimation of population parameters stemming from insufficient survey effort (Lebl et al., 2011; Lin et al., 2018).

4.2 Survival rates

Survival rates of humpback dolphins in LB are generally low, especially in individuals inhabiting Guangdong waters of LB where the survival rates are substantially lower than those estimated for other conspecific populations (Wang et al., 2012, Chen et al., 2018; Zeng et al., 2020) and several other known populations of coastal dolphins (e.g., Currey et al., 2009; Daura-Jorge et al., 2013; Lane et al., 2015). The mark-recapture survival estimates generated by our study are comparable to earlier estimates of non-calf survival rate (0.930) derived from static life-table analyses (Huang et al., 2012; Karczmarski et al., 2017a), where over 90% of samples originated from stranded carcasses recovered in LB between 1995 and 2009 (Jefferson et al., 2012). Despite the methodological differences, the continuously low survival estimates indicate that humpback dolphins inhabiting LB waters must have been under considerable stress over the past decades.

The differences in the estimated apparent survival rates of dolphins from the three sub-regions of LB are well pronounced, which is rarely reported at such a fine spatial scale (few tens of kilometers) for highly mobile animals as cetaceans. Survival rates of adults inhabiting primarily the Guangdong section of LB were estimated to be 6.2% lower than those of dolphins in HK waters, and this difference more than tripled (19.2%) for juveniles. It seems therefore apparent that humpback dolphins with spatial preference for different sub-regions of LB experience different levels of stress, either natural or anthropogenic, affecting their survival.

The individual home range of humpback dolphins (averaged at ~99.5 km², Hung and Jefferson, 2004) comprises only a small fraction of LB estuarine waters (~2100 km²). While various anthropogenic stressors may vary in frequency and intensity across LB, the limited spatial range of individuals and their obligatory inshore distribution (Jefferson and Karczmarski, 2001; Jefferson and Smith, 2016) renders them susceptible to impacts from localized anthropogenic sources at extents corresponding to individual ranging behavior. In general, individuals living in suboptimal habitats are likely to display lower survival probabilities than those in less disturbed areas (Reid et al., 2008; Harris et al., 2021; Ritterson et al., 2021). Similarly, considerable environmental disruption, even if localized, will affect to a greater extend individuals that have used that area as part of their range compared to those whose home range lay outside the disturbed zone, and the likely implications on their survival probability may differ accordingly (Lin et al., 2023).

Since the early 1990s, the PRD economic zone has undergone rapid development and urbanization, followed by a substantial increase in human population which tripled within ten years (Ouyang et al., 2005) and continues to grow. Growing human population increases the need for food and space, which in the PRD has led to expansion of fisheries and subsequent increase in the fishery resource exploitation (Chen et al., 2013), and coastal land reclamation (especially in urban centers at the north and west of LB; Lei et al., 2002; Karczmarski et al., 2016; Zhang et al., 2016). The artificial coastline in these areas is thought to reduce fish spawning, while pollutants and suspended sediment particles from intensive dredging (which are notably higher in the Guangdong section of LB, compared to Hong Kong waters; Hong et al., 1999; Zhou et al., 2004; Chen et al., 2006; Ip et al., 2007), have weakened the fish stock recruitment by affecting the spawning rates and larval survival rates (Cai and Li, 2011). Correspondingly, in the past decade, there has been a major decline (over 70%) in fish biomass in the northernmost LB (Yu et al., 2016), the most degraded area of this estuarine ecosystem. By comparison, the coastal zone of west and southwest Lantau Island in Hong Kong underwent far less development and retains relatively natural inshore habitat (Karczmarski et al. unpublished²), representing apparently the least degraded habitat for humpback dolphins at the outer rim of LB (Karczmarski et al., 2016; Or, 2016).

The spatially-differing and locally depleted prey resources likely hamper the foraging efficiency of humpback dolphins, especially in the north and west of LB (Lin et al., 2021). This corresponds with the dolphins’ decreasing trophic level in the west of LB (Zhang et al., 2019; Ho, 2021), where individuals have to resort to feeding on much smaller prey items (Lin et al., 2021), likely undermining individual fitness. On the other hand, trophic level of humpback dolphins in Hong Kong waters appears relatively stable (Ho, 2021).

The exceptionally high level of pollutants imposes another threat to humpback dolphins in this region (Jefferson et al., 2006; Jefferson et al., 2012; Gui et al., 2014; Gui et al., 2019). Organochlorines are of particular concerns (Jefferson et al., 2006) as these endocrine disruptors are well-known for their adverse effects on reproductive and immune systems of mammals (Wells et al., 2005; Murphy et al., 2018). While it is difficult to assess the direct impacts on reproductive success across the spatial scale of the current study, high level of organochlorines likely contributes to immunosuppression which in turn may affect the dolphins’ survival, as indicated by the overall low survival rates across LB.

Intense maritime traffic in LB adds to the multitude of anthropogenic pressure on the dolphins. While a major international shipping lane lies in the east of LB (next to the Guangdong-HK administrative border), domestic ferry lanes connecting the main metropolitan centers (i.e., Hong Kong, Macau, Zhuhai, Shenzhen) crisscross the estuary. Moreover, small motor vessels that do not follow shipping lanes often cruse at high speeds in inshore waters frequented by the dolphins, causing considerable injuries that can be fatal at times (Chan and Karczmarski, 2019). The impacts of vessel collisions, behavioral disturbance and underwater noise pollution stemming from various sea traffic can be considerable (e.g., Lusseau, 2005; Bejder et al., 2006; Lusseau, 2006; Rako et al., 2013; Pine et al., 2021) and it is thought to be so for humpback dolphins in LB (Wilson et al., 2008; Jefferson and Smith, 2016; Chan and Karczmarski, 2019; Ho et al., 2024). A recent study, however, suggests that prey availability outweighs the risk posed by maritime traffic in determining the ranging behavior of LB dolphins (Pine et al., 2017). While the impact of boat traffic on dolphins’ fitness should not be underestimated, it is possible, perhaps even likely, that the pressure from other human activities, particularly habitat degradation and loss (Karczmarski et al., 2016), and depletion of prey resources in coastal waters (Lin et al., 2021), are more critical in determining the survival of the LB humpback dolphins, as reflected in the spatial variation of apparent survival rates across the estuary.

4.3 Recruitment

To minimize the potential bias of misidentification, only highly distinctive juveniles and adults were included in the analytical database. As such, new identifications resulted from either photographing individuals that have not yet been encountered before or those that were previously unmarked but gained distinctive marks over time (mostly likely juveniles). Therefore, the estimation of entry probabilities (PENT) in POPAN model considers a mixture of entries, including real immigration as well as previously non-distinctive individuals turning distinctive, without distinguishing one from the other (Arnason and Schwarz, 1999). True recruitment rates of young individuals (the combined outcome of successful birth and calf survival) are likely comparable to, if not lower than the estimated PENT for juveniles, which remained low throughout the study period (PENT_juvenile: 0.06-0.10). Given the very high mark-ID ratio (θ = 92.0%), individual identification based on multiple features (using individual-specific pattern of notches on the dorsal fin and dark spots on the upper body), and limited movements of individuals beyond the present study area (see earlier), adult immigration or mark change (change in appearance that is sufficiently large to confuse previously catalogued individual-ID) are unlikely sources of new entries.

Low recruitment rate may be ascribed to low birth rate or calf survival, or – in the case of LB humpback dolphins – more likely a combined effect of both (Jefferson et al., 2012). Most notably, the declining prey resources in LB increase the nutritional stress of humpback dolphins (Lin et al., 2021), which may consequently force the female dolphins to lower their reproductive effort as a trade-off for survival (Gunderson, 1997; Ford et al., 2010). As the fish resource deterioration in LB predates the present study for nearly a decade (see earlier section of the Discussion), the dolphins may have experienced low recruitment rates for many years (Jefferson et al., 2012), which corresponds with the aging population structure noted recently for the LB humpback dolphins (Guo et al., 2020).

Although we cannot quantitatively estimate the net birth and calf survival in the present study, the previously proposed fecundity of 0.286 (Huang et al., 2012), which is unusually high compared to other delphinid populations (0.09-0.23) (e.g. Herzing, 1997; Kogi et al., 2010; Fruet et al., 2015), is very likely overestimated given the current status of the LB dolphins. It is unlikely that a heavily impacted population of long-lived slow-reproducing species would have such high birth and calf survival rates, especially when adult survival rates – which are generally more resilient to environmental fluctuation than survival rates of young individuals (Rotella et al., 2012) – are concerningly low. Therefore, given the latest evidence presented here, earlier population risk assessments (Huang et al., 2012; Miller, 2016; Karczmarski et al., 2017a) based on the estimate of fecundity available at the time, albeit suggesting already alarmingly declining trend, were likely over-optimistic nonetheless. As aged odontocete females reproduce at lower rates (Marsh and Kasuya, 1986; Fruet et al., 2015), the LB humpback dolphins will likely continue to decline unless the majority of currently immature females become reproductive in an improved environment with considerably lower level of stress (Guo et al., 2020).

4.4 Demographic trajectory over past two decades

Since the economic “opening-up” in mainland China in 1980s, the PRD region has been transformed into a highly urbanized industrial hub of the country. The anthropogenic pressure on the estuarine system of LB has been growing ever since, and increasingly rapidly in recent decades (Chen et al., 2013; Lin et al., 2016). Demographic responses to such environmental change by its apex predator, the humpback dolphins, can be particularly revealing of long-term implications. As our mark-recapture study delivers the latest population estimates for the LB dolphins, it offers valuable insights into their current status and recent demographic trends. The overall super-population size estimate (N_T = 928, 95% CI = 823-1046) represents all individuals that have used LB waters as part of their range during the time-period of this study, thus including all new entries as well as individuals that have died during that period. Given the alarmingly low survival rates (see earlier), the number of dolphins that were photographed and died during the study period is likely considerable, and thus the overall super-population size estimate over-represents the latest number of extant animals. Instead, it provides the closest inference of the initial population size during the present study period (year 2011), as it includes all identified individuals that were already in the area at the onset of the study and survived at least until encountered and photo-identified. On the other hand, the last annual super-population size estimate (N_{T_2015} = 723; 95% CI = 592-884) (year 2015) reflects most accurately the latest number of humpback dolphins in waters of LB, as most of new entries have likely been captured with the increasing survey effort while the individuals that have died were excluded.

While there are no other mark-recapture estimates across the entire LB area, two relative abundance estimates, both based on line-transect techniques, were available prior to this study (Jefferson, 2002; Chen et al., 2010; Hung, 2015; see Figure 5). Caution has to be taken, however, when comparing mark-recapture and line-transect estimates given the fundamental differences in the survey and analytical techniques. The estimates produced by these two techniques are only comparable when surveys, especially line-transects surveys, cover the entire range of a population (Daura-Jorge and Simões-Lopes, 2016; Lin et al., 2018), as it was the case with Chen et al. (2008) and Chen et al. (2018) estimates of humpback dolphins in Xiamen Bay. In our case, as the previous line-transect estimates covered a comparable geographic range as our current study, a comparison of the estimated number of dolphins may provide an indication of a demographic trajectory over the past two decades.

Figure 5 The inferred trajectory of Indo-Pacific humpback dolphin abundance in Lingding Bay from 1996 to 2015. The hollow circle and square indicate the super-population size estimate and the annual population size estimate of the last survey year of the present study, respectively. Black solid circles indicate two relative abundance estimates from previous line-transect studies (Jefferson, 2000; Jefferson, 2002; Chen et al., 2010; Hung, 2015). Error bars indicate the 95% confidence interval (CI) of the abundance/(super-)population size estimates. While the previous studies suffer from low precision (very large CI), if their estimates are taken for their face value, exponential regressions with least square fitting indicate accelerated rates of decline from -0.007 during 1996-2005 to -0.059 during 2005-2015.

In the 1990s, it was first estimated that there were 1382 dolphins in waters of LB (inferred from Jefferson, 2000; Jefferson, 2002), although the precision of the estimate was very low (95% CI = 708-2864). Because the data from Guangdong waters were available only for the period from November 1997 to November 1998, this estimate reflects the abundance in 1998. The second survey, which followed the exact same survey routes as Jefferson (2000), was conducted in Guangdong waters from February 2005 to January 2006 (Chen et al., 2010). Combining the abundance estimate from Chen et al. (2010) with that obtained in the same year in Hong Kong waters (Hung, 2015) gives an abundance estimate of 1314 dolphins across the LB area in 2005, albeit again with poor confidence interval (95% CI = 718-2444).

It has to be pointed out that the wide range of confidence interval of the early abundance estimates (as compared to mark-recapture estimates presented in this study; Figure 5) makes it very challenging to assess the true demographic trend. However, imperfect as they are, in the absence of any other estimates in the region, we can only infer the apparent change in dolphin numbers based on the information available. If the lack of precision of the earlier line-transect estimates were to be neglected, exponential least square fitting indicates an instantaneous rate of decline of -0.7% from 1998 to 2005, and an accelerated declining rate of -5.9% from 2005 to 2015 (Figure 5).

The magnitude of acceleration in the population decline (-0.7% vs. -5.9%) is likely exaggerated, especially as the projected rate of decline during the earlier phase based on the previous highly imprecise abundance estimates (1998-2005) was likely underestimated (Huang and Karczmarski, 2014) and lower than that reported by Huang et al. (2012) (r = -0.025). However, the accelerated declining rate in recent years is plausible, even likely, and corresponds with the ever-growing anthropogenic pressure in the PRD region, the apparently aging dolphin population (Guo et al., 2020), and low survival rates of humpback dolphins presented in this study. This puts the effectiveness of current conservation measures in LB (and across the greater PRD) in question, as discussed at length in Karczmarski et al. (2016) and Karczmarski et al. (2017a). Furthermore, the lack of confidence in projecting demographic change using the existing abundance estimates underscores the urgent need for a region-wide long-term monitoring program that could deliver accurate and precise population estimates. Only then will we be able to better understand the dolphins’ demographic response to environmental change and how it determines their long-term viability in waters of Lingding Bay.

5 Concluding comments

The inland region surrounding Lingding Bay has long been central to agriculture and industrial development throughout Chinese history (Weng, 2007; Zong et al., 2009). Currently, the urban and industrial region of the PRD accounts for at least 36% of the total trade of China (Ma and Zhang, 2009) and has become the world’s largest metropolitan area by both land area and human population size (World Bank, 2015). When faced with priority economic goals, environmental protection that is fundamental to preservation of species and habitats (Storch, 1997; Gardner et al., 2007) becomes extremely challenging (Karczmarski et al., 2017a). Recently, an ambitious development plan locally known as the ‘Greater Bay Area’ has been formulated and aims at expanding the economic development goals across the greater PRD region at unprecedented scales, which raises concerns about the fate of the region’s coastal and estuarine ecosystem in not-so-distant future.

In such an environmentally challenging situation, it is particularly important to strategize the conservation of local biodiversity with the best available science. An ongoing study (Karczmarski et al. unpublished²) has identified three major core areas of humpback dolphin habitat within LB, all of them located in inshore waters off relatively natural, uninhabited, and anthropogenically unaltered shores. Coastal waters off Lantau Island in particular comprise the largest continuous dolphin habitat in the entire LB, which matches the spatial variation in survival rates of humpback dolphins described in the present study. It is therefore of the utmost importance that all remaining dolphin core habitats in LB, majority of which are within Hong Kong waters, be given the highest conservation priority and be protected with an effective management plan as a matter of urgency. In the long-term, it is equally critical to identify and mitigate the major threats in Guangdong waters of LB, so to improve dolphin survival rates and recover the population’s healthy age structure. At present, however, in the absence of quantitative information on anthropogenic stress, it remains impractical to assess separately various sources of impacts and the dolphins’ apparent demographic response. Instead, a region-wide coordinated conservation plan is much needed to secure the ecological viability of dolphin key habitats and maintain the minimum viable number of humpback dolphins in LB waters. A continuous effort of photo-ID surveys combined with spatial analyses of habitat and resource selection is particularly important as it offers both research and management tools; it provides means to maintain vigilance of the demographic trend and data to assess the main stressors in a quantitative management-applicable manner that can inform the decision-making process of prioritizing conservation effort based on sound scientific evidence.

Data availability statement

The data used in this article may be made available by the authors upon a reasonable request.

Ethics statement

Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because the present study was conducted through distant observation with no harm to the animals, and also in accordance with the local legislation and institutional requirements.

Author contributions

WL, LK and SC provided equal intellectual contribution. WL: conceptualization, data collection and management, methodology, formal analysis, data visualization, manuscript writing (initial draft), funding acquisition. LK: conceptualization and implementation of a multi-year mark-recapture study and analytical framework, data collection, validation of results, manuscript writing/revisions/edits, funding acquisition, project management, supervision of research. SC: conceptualization, data collection and management, methodology, data visualization, manuscript review/revisions/edits, funding acquisition. RZ: data collection and data processing, funding acquisition. Y-WH, YM: data collection and data processing. All authors contributed to the article and approved the submitted version.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Guangdong Natural Science Foundation grant number 2018A030313870; General Research Fund (GRF) of the Research Grants Council (RGC) of Hong Kong grant number HKU768110M and HKU17100015M; Marine Ecology Enhancement Fund (MEEF) grant number MEEF2017015, MEEF2017015A, MEEF2017015B and MEEF2017015C; Ocean Park Conservation Foundation Hong Kong (OPCFHK) grant number AW04.1617, MM01.1921, MM01.1415 and MM03.1213, the Alashan Society of Entrepreneurs and Ecology (SEE) and the Biodiversity Investigation, Observation and Assessment Program (2019-2023) of Ministry of Ecology and Environment of China (No. 9-38-27-1), and Hong Kong Airport Authority (AAHK) and Hong Kong International Airport Business Community funding CWD1_1415.

Acknowledgments

This project was partially performed during WL’s post-doctoral appointment at the University of Hong Kong. WL wishes to thank Philip Hammond for valuable suggestions at an early stage of this study; and Ruilian Zhou, Ye Lin and Ge Lin for their invaluable support.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2024.1252661/full#supplementary-material

Footnotes

1. ^ Lingding Bay is also referred to as the Pearl River Estuary (PRE; e.g., Jefferson 2000), or the PRE proper, as it represents the main body of the Pearl River estuarine system.
2. ^ Karczmarski, L., Chan, S.C.Y., Ho, Y-W., Lin, W. (unpublished manuscript). Socio-spatial dynamics of a coastal delphinid in an anthropogenically degraded estuarine seascape.

References