- 1Department of Environmental Biology and Fisheries Science, National Taiwan Ocean University, Keelung, Taiwan
- 2Center of Excellence for the Oceans, National Taiwan Ocean University, Keelung, Taiwan
The swimming crabs is a crucial predator species in benthic habitats and a high value in commercial fishery industries in subtropical and tropical Asia. The climate variability caused by El Niño–Southern Oscillation (ENSO) events has substantial impacts on the catch and habitat of this species. In this study, a weighted habitat suitability index (HSI) model was constructed using logbooks and voyage data records from Taiwanese crab vessels (2013–2019) with the addition of environmental variables to examine the influence of ENSO events on catch rates (CRs) and habitat suitability for Charybdis feriatus, Portunus pelagicus, and Portunus sanguinolentus in the Taiwan Strait (TS). The autumn (September–October) is the major fishing season for catching these three swimming crab species in the TS. A high CR of P. sanguinolentus was observed across the TS, whereas high CRs of P. pelagicus and C. feriatus were recorded in areas in the southern and northern TS, respectively, during autumn. Moreover, the CRs for C. feriatus and P. pelagicus were higher (>7.0 and >8.0 kg/h) during La Niña events, with the increase being more than 40.0% compared with the CRs under normal and El Niño events in autumn. For P. sanguinolentus, the CRs were higher during both La Niña and El Niño events (>8.0 kg/h) compared with normal years. The high CRs for C. feriatus and P. sanguinolentus during autumn in La Niña years co-occurred with high sea temperature and low salinity, whereas the high CR of P. pelagicus co-occurred with high sea temperature and high salinity. Furthermore, the high CRs for C. feriatus and P. pelagicus were observed in areas with high HSI in the La Niña years but were distributed more widely with a lower HSI during normal and El Niño years. The low CRs for C. feriatus and P. pelagicus during normal and El Niño years and the low CR for P. sanguinolentus in normal years during autumn were highly consistent with substantial shrinkage of suitable habitats. Our findings suggest that ENSO events strongly affected the catch and habitat suitability of C. feriatus, P. pelagicus, and P. sanguinolentus during autumn in the TS.
Introduction
The Taiwan Strait (TS) is located in the tropical to subtropical western Pacific region and is an important channel in the west Pacific Ocean for transporting water and chemical constituents from the East China Sea to the South China Sea. Crab species in the TS are diverse, with more than 250 species found in East Asia (Ng et al., 2001). The crucifix crab Charybdis feriatus, blue swimming crab Portunus pelagicus, and three-spot swimming crab Portunus sanguinolentus are three of the most crucial commercial crustacean species in the TS (Hsueh and Hung, 2009; Naimullah et al., 2020). C. feriatus and P. pelagicus are commonly distributed in the north and south of the TS, respectively, whereas P. sanguinolentus is widely distributed across the TS (Naimullah et al., 2020).
Charybdis feriatus, widely distributed in the Indo-Pacific region, is usually found in sublittoral zones on muddy and sandy bottoms and rocky and stony coasts, including coral reef flats, at depths of approximately 10–60 m (Abelló and Hispano, 2006). P. pelagicus thrives in sandy to sandy–muddy substrates in shallow waters (depth, 50 m), including areas near reefs and mangroves as well as in seagrass and algal beds; it is commonly found in tropical and subtropical estuaries and nearshore habitats (Hosseini et al., 2012). P. sanguinolentus is widely distributed in the Indo-Pacific region from the eastern coast of South Africa to waters around Hawaii and typically inhabits sandy oceanic habitats up to a depth of 30 m (Rasheed and Mustaquim, 2010).
For efficient resource management, understanding the spatial and temporal variability of key environmental variables in commercial fishing grounds is crucial to help identify species distributions and their habitat preferences. This information can be used for implementing sustainable fisheries through harvest strategy planning (Kruse et al., 2010), ecosystem management (Szuwalski and Punt, 2013; Kunsook et al., 2014; Szuwalski et al., 2021), spatial management (Frid et al., 2016), and bycatch reduction (Morris et al., 2011). Environmental variables – such as water temperature, salinity, water depth, and increased productivity – may be critical factors driving changes in population patterns (de Lestang et al., 2010; Spencer et al., 2019; Naimullah et al., 2020), reproductive and recruitment strategies (Johnston et al., 2011; de Andrade et al., 2015), and spatiotemporal distributions (De Anda-Montañez et al., 2016).
Climate change is having a considerable impact across marine ecosystems, latitudes, and trophic levels (Scheffers et al., 2016) and is predicted to redistribute the global fishery catch (Cheung et al., 2009). The effects of climate change on global fishery can be grouped into two broad categories: changes in stock productivity, which affect yields and potential profits, and changes in stock distribution, which influence the regions in which fish can be caught and who might catch them (Erauskin−Extramiana et al., 2019). The most extensively investigated climatic events having an impact on fishery are those with an interannual scale, such as the El Niño–Southern Oscillation (ENSO; Lehodey et al., 2006). The ENSO occurs in the tropical Pacific Ocean and is caused by an interaction between surface wind pressure and sea surface temperature (SST) changes (McPhaden et al., 2006); this intermediary leads to climate fluctuations between years. In La Niña events, the SST is lower than normal in the eastern and central tropical Pacific, whereas during El Niño events, it is higher (Lau and Nath, 2006). ENSO events have been reported to have various positive and negative effects on the marine population in the TS (Lan et al., 2014; Ho et al., 2016).
Crustacean fisheries are considered crucial in many countries. Climate change poses severe environmental and economic threats (Azra et al., 2020). Moreover, the impact of climate changes on the fluctuations of fish populations, including crabs (Meynecke et al., 2012a,b; Chandrapavan et al., 2019; Dvoretsky and Dvoretsky, 2020), and on fisheries has been increasingly reported and continues to be critical (Lehodey et al., 2006; Lehodey et al., 2020). However, the effects of climate change on crab communities in the waters of Taiwan have not been sufficiently documented. Therefore, this study investigated the effects of ENSO events on the catch rates (CRs) and habitat preferences of C. feriatus, P. pelagicus, and P. sanguinolentus in the TS. The results can help implement sustainable crab fisheries in the future through improved fishery management and conservation planning in the TS.
Materials and Methods
Conceptual Model
The weighted habitat suitability index (HSI) was calculated using environmental factors – the SST, bottom temperature (BT), sea surface height (SSH), chlorophyll-a (Chl-a), and salinity. The HSI model coupled with a climatic and environmental variability analysis can provide insights into climate-driven habitat suitability variation under different anomalous events for swimming crab species. To elucidate the mechanisms through which climate change affects crab habitat suitability, oceanographic characteristics and suitable habitat areas for the three swimming crabs in the TS during various ENSO events were compared (Figure 1).
Swimming Crab Fishery Data
Fishery data on traps from 2013 to 2019 were collected from the voyage data records and logbooks of 218 Taiwanese crab vessels (ranges from 61 to 118 vessels across years with 5–200 metric tons) operating in the TS (Supplementary Table 1 and Supplementary Figure 1). The fishing data comprised daily fishing positions for 0.1° spatial grids, including latitude and longitude, fishing date, soaking time, total catch (kg), and species. C. feriatus, P. pelagicus, and P. sanguinolentus were selected for analysis and other species recorded were excluded. Swimming crabs were collected using circular crab traps with frozen mackerel as bait placed in the center of the trap.
Monthly CRs were calculated as the cumulative weight of each crab species caught from all fishing vessels for a month divided by the cumulative soaking time of all set lines in that month to adjust for fishing effort. Monthly CRs were calculated for individual 0.1° grid squares across the study region by using the following expression:
where i represents the latitude and j represents the longitude of each 0.1° spatial grid square.
We constructed the distribution maps of CRs for the three swimming crab species using Interface Descriptive Language version 7.0.
Oceanic Niño Index
To obtain tropical climate signals, we considered the Oceanic Niño Index (ONI), an ENSO-related index and examined the relationship between the mean CRs of C. feriatus, P. pelagicus, and P. sanguinolentus. The ONI is currently used by National Oceanic and Atmospheric Administration to provide forecasts of El Niño and La Niña activity and is defined as the 3-month running mean of SST anomalies in the Niño 3.4 region of the equatorial Pacific (5°N–5°S and 120°–170°W). An El Niño and La Niña event is defined to occur when the ONI index higher than +0.5°C and lower than -0.5°C for 5 consecutive months, respectively. An ONI in the range ±0.5°C is defined as a normal year. The definition and the ONI data were sourced from the NOAA climate prediction center1.
Marine Environmental Data
Marine environmental data – daily SST, BT, SSH, and sea surface salinity (SSS) – were extracted from the Copernicus Marine Environment Monitoring Service2 from 2013 to 2019 at a horizontal resolution of 1/12°. The SST, SSH, and SSS were obtained across the upper 0.5 m of surface water and potential temperature at the seafloor of each grid for the BT.
Monthly Chl-a data were downloaded from the NASA Aqua satellite3, which has a sensor that detects the concentration of Chl-a in the world’s oceans. The level 3 data map image of Chl-a data has a monthly temporal resolution and a spatial resolution of 4.6 km (at the equator).
The daily environmental data were then calculated mean monthly on a spatial grid of resolution 0.1° to fit the fishery data. Intercorrelations among environmental data were tested before application of statistical models to determine the correlation coefficients among the variables. The test for intercorrelations among all environmental variables obtained low correlation coefficients (r2 < 0.25).
The relationship of each swimming crab CRs with the environmental variables in ENSO events on the basis of crab trap fishery data locations from 2013 to 2019 in the TS were analyzed using a one-way ANOVA. Pairwise post hoc comparisons were conducted using Tukey’s tests. For the one-way ANOVA, we tested and found that required assumptions were met with regard to the homogeneity of variance, normal distribution of errors, and independence of errors.
Weighted Habitat Suitability Index Model Development
The habitat suitability modeling approach helps efficiently obtain the distribution pattern of potential habitats for marine species, including crustaceans (Chen et al., 2019). The HSI is a numerical index between 0 and 1, where 0 indicates an unsuitable habitat and 1 represents an optimal habitat. Suitable habitat is closely related to the area that had the highest catches (Chang et al., 2019). By considering CRs > 60% as the high catches, the relationships between the CRs for three swimming crabs and each environmental variable were presented as curves of the suitability index (SI, 0 to 1). The formula for the SI was as follows:
where catch ratesij denotes the cumulative fishing efforts in the ith interval of the range of environmental variable j, and max(catch ratesij) indicates the maximum cumulative fishing efforts. The relationships between the SI and environmental variables for each swimming crab species were determined as follows (Yu et al., 2018):
where a and b are the model parameters to be estimated, and i denotes the environmental variables considered in this study.
The fitted SI models were combined into an integrated HSI model by using the arithmetic mean model; this approach is the most commonly used and yields optimal performance (Silva et al., 2016; Yu et al., 2018). Studies have reported that setting weights for different environmental variables can significantly influence HSI modeling.
Partial least squares regression (PLSR) is a bilinear calibration method that compresses data by reducing many measured collinear spectral variables to a few non-correlated principal components. PLSR is a valuable tool for exploring datasets with numerous predictor variables (Carrascal et al., 2009). PLSR considers the latent structure of the predictor variables (CRs > 60%) and response variables (five environmental parameters) and selects the model with predictor variables that best describe the variation in the response variables (Wold et al., 2001). PLSR variable influence on projection (VIP) scores were generated to assess the relative importance of the predictors in the overall model, and the relative weighting was based proportionally on the VIP score.
Therefore, in this study, five weights based on VIP scores were used for each swimming crab species in the HSI model development. The weighted HSI arithmetic mean model can be described as follows:
where SIi is the suitability index or value associated with environmental variable i, n is the number of environmental variables considered in this study, and wi is the weight based on the VIP scores of each environmental variable for each swimming crab species. The spatial distribution of the weighted HSI obtained from the models in the major fishing seasons was mapped, and the influence of ENSO events to habitat suitability of C. feriatus, P. pelagicus, and P. sanguinolentus was investigated.
Results
Time-Series Trends and Spatial Distribution of Catch Rates
The monthly mean catch and CRs in the TS for each species exhibited similar trends (Figures 2A–C). The highest monthly mean catch and CRs for the three swimming crabs were observed from September to November (autumn) and then gradually decreased until June to July. Because the major fishing season with the highest catch and CRs for C. feriatus, P. pelagicus, and P. sanguinolentus was the autumn, this season was selected as the study period to determine the effects of climate events.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/763543/fmars-08-763543-HTML/image_m/fmars-08-763543-g002.jpg)
Figure 2. Monthly mean catch (grey bar) and CRs (black line) of Charybdis feriatus (A), Portunus pelagicus (B), and Portunus sanguinolentus (C).
The CR distribution maps for C. feriatus revealed that higher CRs were distributed in the northern part of the TS from 25 to 26°N and 120 and 122°E during autumn (Figure 3A). For P. pelagicus, higher CRs were distributed further to the south of the TS between 22 and 24°N and extended to the northeast at 25–26°N (Figure 3B). P. sanguinolentus was widely distributed across the TS, with higher CRs extending to northeastern areas during autumn (Figure 3C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/763543/fmars-08-763543-HTML/image_m/fmars-08-763543-g003.jpg)
Figure 3. CRs (kg/h) of (A) Charybdis feriatus, (B) Portunus pelagicus, and (C) Portunus sanguinolentus from 2013 to 2019 during autumn.
Effects of El Niño–Southern Oscillation Events on Catch Rates of Swimming Crabs
Based on the definition of ENSO events, La Niña events occurred during the autumns of 2016 and 2017, and El Niño events occurred during the autumns of 2014, 2015, and 2018 (Figure 4). Compared with the mean CRs of ENSO events, the CRs for C. feriatus and P. pelagicus were higher (>7.0 and >8.0 kg/h) in La Niña events and were more than 40.0% higher compared with the CRs in normal and El Niño events during autumn (Figure 4). For P. sanguinolentus, the effect of ENSO events was not clear; the CRs were higher in both La Niña and El Niño events (>8.0 kg/h) compared with those in normal years.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/763543/fmars-08-763543-HTML/image_m/fmars-08-763543-g004.jpg)
Figure 4. ENSO events and mean CRs for Charybdis feriatus, Portunus pelagicus, and Portunus sanguinolentus during autumn in the TS.
The mean CRs were further calculated for the northern TS, middle TS, and southern TS (Figure 3) to understand the effect of ENSO events on the distribution of three swimming crabs during autumn (Figure 5). The CRs for C. feriatus were higher in the northern TS than in the middle or southern TS, especially in 2017, a La Niña year (Figure 5A). The lowest CRs for C. feriatus in the northern TS occurred in the El Niño events of 2015, and the CRs were lower than 1.7 kg/h in the southern TS throughout the study period. For P. pelagicus, the mean CRs were higher in the southern TS and lower in the middle and northern TS (Figure 5B). In the southern TS, the highest CR for P. pelagicus was recorded in the La Niña year 2016, whereas the lowest CRs occurred in the normal year 2019. CRs < 1.0 kg/h were recorded for P. pelagicus in middle and northern TS for each climate event from 2013 to 2019. The CRs for P. sanguinolentus were higher in the northern and southern TS throughout the study period than in the middle TS, especially in 2017, a La Niña year (Figure 5C). The lowest CRs for P. sanguinolentus (<3 kg/h) occurred in the middle TS during the normal year 2019.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/763543/fmars-08-763543-HTML/image_m/fmars-08-763543-g005.jpg)
Figure 5. Mean CRs for Charybdis feriatus (A), Portunus pelagicus (B), and Portunus sanguinolentus (C) in the ENSO events during autumn in three TS areas: the northern TS (25–27°N and 118–122°E), middle TS (24–25°N and 118–121°E), and southern TS (22–24°N and 118–121°E).
Relationship Between Swimming Crab Catch and Environmental Variables
Obtained through PLSR, the relative weighting of important environmental variables of the CRs for C. feriatus, P. pelagicus, and P. sanguinolentus based on VIP scores are presented in Table 1. The BT explained the highest VIP score for C. feriatus, followed by the SSH, SSS, and Chl-a, whereas the SST explained the lowest (Table 1). Similarly, for P. sanguinolentus, the BT explained the highest VIP score, followed by the Chl-a, SSH, SST, and SSS. By contrast, for P. pelagicus, SSS explained the highest VIP score, followed by the SST, BT, and SSH, with Chl-a explaining the lowest score.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/763543/fmars-08-763543-HTML/image_m/fmars-08-763543-t001.jpg)
Table 1. Relative weighting of the sea surface temperature (SST), bottom temperature (BT), sea surface height (SSH), sea surface salinity (SSS), and chlorophyll-a (Chl-a) for Charybdis feriatus, Portunus pelagicus, and Portunus sanguinolentus based on the VIP scores.
The box-and-whisker plots illustrate the relationship of each swimming crab species with the environmental variables during autumn in ENSO events on the basis of crab trap fishery data locations from 2013 to 2019 in the TS (Figure 6). Statistically, ANOVA showed significant effects of environmental variables during autumn in ENSO events on CRs of all swimming crab species (Figure 6). The results revealed that high CRs for C. feriatus in La Niña co-occurred with higher SST, BT, and SSH and lower Chl-a and SSS (Figures 6A–F). However, lower CRs for C. feriatus in normal and El Niño years revealed no significant differences and co-occurred with lower SST, BT, and SSH and higher Chl-a.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/763543/fmars-08-763543-HTML/image_m/fmars-08-763543-g006.jpg)
Figure 6. Box-and-whisker plots showing the effects of mean environmental variables on CRs for three commercial swimming crabs during autumn in different ENSO events on the basis of crab trap fishery data location: (A–F) Charybdis feriatus, (G–L) Portunus pelagicus, and (M–R) Portunus sanguinolentus. The horizontal black line in the box denotes the mean value, and the bottom and top of the box indicate the standard error. The ends of the whiskers represent the 95% confidence interval. The horizontal red line is the mean value during autumn from 2013 to 2019. One-way ANOVA indicate significant differences at right top of each box and different subscripts above the box-and-whisker structure indicate significant differences detected in Tukey’s tests (p < 0.05).
For P. pelagicus, the CRs were higher in La Niña years and co-occurred with higher SST, BT, SSH, SSS, and Chl-a (Figures 6G–L). Similar to C. feriatus, no significant differences were noted in lower CRs for P. pelagicus in the normal and El Niño years, and they co-occurred with lower SST, SSH, SSS, and Chl-a. Moreover, for P. sanguinolentus, higher CRs were recorded in La Niña and El Niño years and co-occurred with higher BT and SSH, whereas the SST and Chl-a were higher in La Niña years and the SSS was not significantly different between various ENSO events (Figures 6M–R).
Suitability Indices of Environmental Variables
The fitted SI model was established for environmental variables in autumn, and all the models were significant (p < 0.01) with low root-mean-square errors (RMSEs) and high correlation coefficients (Table 2). Different optimal ranges for each environmental variable were identified for each swimming crab species when the SI curve was used (Figure 7). The optimal range for SSTs with higher SI (>0.6) for C. feriatus and P. sanguinolentus varied from approximately 25.0 to 29.0°C, whereas that for P. pelagicus varied from 27.0 to 30.0°C (Figure 7A). For BTs, the optimal ranges were approximately from 19.0 to 26.0°C for C. feriatus and P. sanguinolentus and from 24.0 to 27.0°C for P. pelagicus (Figure 7B). The optimal range for SSHs for C. feriatus and P. sanguinolentus was approximately from 0.52 to 0.64 m, whereas that for P. pelagicus was from 0.64 to 0.76 m (Figure 7C). Moreover, the optimal range for SSSs was 32.8–33.7 PSU for C. feriatus and P. sanguinolentus and 33.1–34.0 PSU for P. pelagicus (Figure 7D). However, the suitable Chl-a range was higher for C. feriatus than for P. pelagicus and P. sanguinolentus (0.7–1.4 mg/m3; Figure 7E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/763543/fmars-08-763543-HTML/image_m/fmars-08-763543-t002.jpg)
Table 2. Autumn fitted SI model for each environmental variable for Charybdis feriatus, Portunus pelagicus, and Portunus sanguinolentus in the TS.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/763543/fmars-08-763543-HTML/image_m/fmars-08-763543-g007.jpg)
Figure 7. Fitted SI curves based on the relationship between the fishing effort for Charybdis feriatus, Portunus pelagicus, and Portunus sanguinolentus and each environmental variable: (A) sea surface temperature (SST), (B) bottom temperature (BT), (C) sea surface height (SSH), (D) sea surface salinity (SSS), and (E) chlorophyll-a (Chl-a). The intersections of the SI curve and the dashed line (SI = 0.6) denote the optimal range for each environmental variable.
Habitat Suitability Variation With El Niño–Southern Oscillation Events
The study of HSI maps revealed that suitable habitat areas (HSI > 0.6) of C. feriatus in normal (2013) and El Niño (2014) years were smaller in the northern TS (Figures 8A,B). By contrast, the suitable habitat areas of C. feriatus were higher in the northern TS and extended to northern areas (north to 26°N) in La Niña year 2017 (Figure 8C), with higher CRs and a higher HSI, indicating that fishing locations moved further north in the TS.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/763543/fmars-08-763543-HTML/image_m/fmars-08-763543-g008.jpg)
Figure 8. Mean CRs of Charybdis feriatus (A–C), Portunus pelagicus (D–F), and Portunus sanguinolentus (G–I) overlaid with selected suitable habitats (in high-CR years for each species) during autumn in the different ENSO events.
For P. pelagicus, the selected habitat maps show suitable areas in the southern TS (south to 24°N) in each ENSO event (Figures 8D–F). However, the suitable habitat area was smaller during autumn in the normal (2013) and El Niño (2014) years than in the La Niña years. Moreover, a smaller HSI indicated that fishing locations moved to the northern TS (25–24°N) during the normal and El Niño years (Figures 8D,E).
The HSI maps reveal that the suitable habitat areas of P. sanguinolentus during autumn in a normal year (2013) were lower in the northern TS and further to the south of the TS (Figure 8G), whereas the suitable habitat area was higher in the northern TS during autumn in the La Niña years (2017) (Figure 8I). By contrast, the suitable habitat of P. sanguinolentus was widely distributed across the TS during autumn in El Niño years (Figure 8H).
A further comparison of the interannual variability of CRs with HSI in ENSO events (Figure 9) revealed that higher CRs for C. feriatus and P. pelagicus were distributed in the higher HSI areas between 25 and 26°N and between 22.5 and 23.5°N, respectively, in the La Niña years (Figures 9A,B). However, the higher CRs areas were distributed more widely further north of 26.5°N for C. feriatus and 23.5°N for P. pelagicus during normal and El Niño years. The high CRs for P. sanguinolentus were widely distributed in the TS in each climate event (Figure 9C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/763543/fmars-08-763543-HTML/image_m/fmars-08-763543-g009.jpg)
Figure 9. Mean trends in CRs and SI per latitude of Charybdis feriatus (A), Portunus pelagicus (B), and Portunus sanguinolentus (C) during different climate events in the TS.
Discussion
High Catch Seasons and Distribution Patterns
This study demonstrated that autumn is a crucial season for obtaining large catches and high CRs for C. feriatus, P. pelagicus, and P. sanguinolentus in the TS. Furthermore, high CRs for P. sanguinolentus were widely distributed in the TS, whereas those for P. pelagicus and C. feriatus were found in two main areas in the southern and northern TS during autumn. The optimal season for crab fishing is dependent on various factors such as population size, food availability, relevant environmental factors, and the type of gear used (Olsen et al., 2019; Spencer et al., 2019; Naimullah et al., 2020). Variation in CRs is most likely due to environmental factors affecting recruitment success over time and the geographic distribution of swimming crabs in the TS.
The decrease in catches and CRs for C. feriatus, P. pelagicus, and P. sanguinolentus from March to July may have been related to the spawning seasons and recruitment process. Brachyuran crabs inhabiting tropical waters usually breed throughout the year, whereas those found in temperate waters breed only in certain months (Warner, 1977). The spawning season of the swimming crab in the East China Sea was recorded from March to July (Song et al., 1989; Wu et al., 2007). Moreover, because of low water temperatures during winter and spring in the East China Sea, berried females do not hatch until late April (Song et al., 1989; Xie et al., 2002). The berried swimming crab mostly remains buried in the sand and carries its eggs until the larvae hatch as zoeae, which drift about in the water current (Rasheed and Mustaquim, 2010). During spawning seasons, female swimming crabs have a period of low catchability when they cease feeding and become inactive in preparation for the molt and are not easily catchable with baited traps or drop nets (Soundarapandian et al., 2013). Therefore, we suggest that autumn is the preferred season for fishing of these three swimming crab species in the TS.
The higher CRs for P. pelagicus and C. feriatus in the southern and northern TS during autumn may have been due to the habitat and environment being conducive to spawning, recruitment, and growth. After the hatched zoea transform into juvenile crabs, most of those in coastal habitats settle there or migrate to estuaries as conditions in the estuaries become favorable. Moreover, some studies have shown that swimming crabs that inhabit marine embayments do not leave these marine environments to spawn (Potter and De Lestang, 2000; Polity et al., 2011). The wide distribution of P. sanguinolentus in the TS may be primarily influenced by sea currents. Most P. sanguinolentus larvae hatch as zoea, which drift in the water current before inhabiting benthic habitats and growing into juvenile crabs (Rasheed and Mustaquim, 2010). Throughout the year, the strong Kuroshio Branch Current and South China Sea Current in the TS during autumn may affect the movement of P. sanguinolentus in the northern area. P. sanguinolentus typically inhabit the sandy and muddy bottom in shallow areas (Sumpton et al., 1989; Carpenter et al., 1997), whereas adult and berried females often migrate to deeper waters for spawning (Yang et al., 2014), which results in a widespread distribution of this species.
Influence of Environmental Variables on Catch Rates of Swimming Crabs
Understanding the relationships between environmental variables and fundamental processes is crucial to explaining variation in CRs. Although many environmental variables affected the abundance of the three investigated swimming crabs in the coastal waters of Taiwan (Naimullah et al., 2020), each factor contributed differently to the stock dynamics. BT had a positive impact (>21.7°C) on the habitat suitability for C. feriatus and P. sanguinolentus in autumn, implying that rising water temperatures are favorable to the formation of potentially high-quality habitats. Temperature is a crucial factor influencing the physiological and ecological properties of brachyuran crabs (Azra et al., 2020). The variation in water temperature influences molt timing and frequency, growth rates, size at maturity, spawning period, larval development, geographic range, and survival (Defeo and Cardoso, 2002; Ikhwanuddin et al., 2012; Green et al., 2014; Ryer et al., 2016; Waiho et al., 2016; Chandrapavan et al., 2019).
Water temperature is known to modify reproduction and maturation in crustaceans, particularly crabs (Olson et al., 2018). Cooler temperatures generally retard growth and delay maturity, causing animals, including crabs, to begin maturation when they are larger (Leffler, 1972; Azra et al., 2020). By contrast, rearing crabs at higher water temperatures may benefit their growth, molting, and time to maturity (Kuhn and Darnell, 2019). Reducing the size of mature female crabs may negatively affect reproduction because body size is a factor limiting fertility (Green et al., 2014). Environmental variations may retard brachyuran crab growth (Azra et al., 2020). High CRs for C. feriatus and P. sanguinolentus in autumn were also found when the BT was higher during La Niña years, which affected the catch and life cycle processes of these species. Moreover, the La Niña event in Northern Territory Australia showed that intense rainfall and warm waters may explain between 30 and 40% of the Scylla serrata catch variability (Meynecke et al., 2012a).
The SSS, in the range 33.1–34.0 PSU, had a positive impact on habitat quality for P. pelagicus and was higher in La Niña years in the fishing positions. Batoy et al. (1987) also discovered that salinity was one of the most critical environmental factors affecting the reproductive cycle of P. pelagicus species, with higher salinity more favorable for breeding. Furthermore, the optimal range of the SSS was 32.8–33.7 PSU for C. feriatus and P. sanguinolentus in the TS. In Haizhou Bay, salinity significantly influenced the distribution of C. bimaculata, indicating that the preferred salinity range was 29–31 PSU (Luan et al., 2018). By contrast, P. sanguinolentus females mainly inhabit deeper waters and prefer higher salinity than male crabs (Campbell, 1986). Many crustacean fisheries are located in coastal and estuarine areas where the salinity fluctuates due to rainfall and freshwater in-flow (Meynecke et al., 2012a,b). Crabs are generally restricted by a lower salinity threshold; therefore, both adults and larvae avoid low salinity by modifying their behavior such as migration (Gibson and Najjar, 2000; Charmantier et al., 2001).
The Chl-a concentration explained one of the vital environmental factors for P. sanguinolentus. Signa et al. (2008) determined that the density of the swimming crab Polybius henslowii is strongly related to the Chl-a concentration, suggesting higher accumulation in locations with higher production. Moreover, the highest CRs for P. sanguinolentus in the TS were mostly in waters with Chl-a concentration >0.5 mg/m3 (Naimullah et al., 2020). The essential fishing ground was mainly related to the upwelling area consisting of a high Chl-a concentration and low SST that carried large amounts of nutrients from the bottom layer to the ocean’s surface layer (Hu and Wang, 2016) and is a crucial factor affecting the diet of P. sanguinolentus in TS. The SSH, a proxy for the current and temperature balance, was one of the critical factors affecting the swimming crab density. Liu et al. (2001) reported that changes in thermocline depth were related to the SSH. In addition, an increase in the SSH would increase the temperature, thus affecting the crab habitat. A greater SSH (>0.60 m) in La Niña years yielded a positive impact on the CRs for the three swimming crabs in the TS. Srisurichan et al. (2005) reported a similar outcome for spiny lobsters; they reported a positive response of CRs for legal-sized spiny lobsters to the SSH. They noted that swell might produce conditions favoring better protection from predators by increasing bottom turbidity and food availability.
El Niño–Southern Oscillation Events and Suitable Habitats
A modeling approach for evaluating the quality of crustacean species’ habitats is attracting considerable attention and has been increasingly applied in fisheries research (Luan et al., 2018; Chen et al., 2019; Naimullah et al., 2020). The El Niño and La Niña events in the Pacific Ocean represent considerably complex climate variability (Wang et al., 2017). The present study found that ENSO events have substantial impacts on the environmental conditions of fishing grounds for C. feriatus, P. pelagicus, and P. sanguinolentus and change the suitability of habitats in the TS. Higher CRs were observed in La Niña years than in normal and El Niño years for C. feriatus and P. pelagicus, whereas no significant effect of ENSO events was noted on the CRs and distributions of P. sanguinolentus. Furthermore, higher CRs for C. feriatus and P. pelagicus were concentrated in the higher-HSI areas in the La Niña years but distributed more widely in areas with lower HSI during normal and El Niño years. These findings indicate that fishing vessels require more time and oil to search the fishing grounds for C. feriatus and P. pelagicus during normal and La Niña years.
Climate-induced habitat suitability changes may strongly influence fisheries (Kuczynski et al., 2018). For example, in the southern TS, pelagic species resources are susceptible to climate change, whereas benthic species are mainly insensitive to climatic factors (Hsiao et al., 2021). This variability may explain changes in CRs in suitable and optimal habitat areas caused by ENSO events in the TS. In each season, the Kuroshio Branch proportions are the highest during El Niño periods, and this leads to higher salinity than that during La Niña events in the TS (Huang et al., 2015). Concentrations of most nutrients, hydrogen ions, and Chl-a are higher in the autumn of a La Niña year than in the autumn of a normal or El Niño year in the southern TS (Huang et al., 2015). However, the SST was found to be lower in an El Niño summer than in a La Niña summer along the central TS northeastward and up to 24°N (Kuo and Ho, 2004). Thus, we suggest that higher nutrient levels and sea temperature lead to greater SSH in autumn and more extensive optimal habitat areas, higher CRs, and greater growth of C. feriatus and P. pelagicus in the La Niña year in the TS. However, the CRs for P. sanguinolentus are not influenced by ENSO events due to its wide distribution and optimal habitat areas and these species being able to strongly adapt to changing environments.
Numerous environmental factors affect crabs throughout their complex life cycle (Green et al., 2014), and lagged effects on swimming crabs have been reported (Meynecke et al., 2012a; Chandrapavan et al., 2019). Meynecke et al. (2012a) indicated that the best environmental predictor of S. serrata catch variability was the Southern Oscillation Index with a 6–9-month lag; negative impacts on catches were observed after heavy and prolonged flooding with a 1–2-year lag effect. However, in Western Australia, trap CRs and the monthly commercial catch from both trap and trawl fleets declined to historically low levels almost 6 months after extreme marine heatwave phenomena in the 2011 La Niña (Chandrapavan et al., 2019). In addition, swimming crabs mostly have a maximum lifespan of approximately 3 years, and the abundance of crabs is highly dependent on the success of recruitment and the survival of a small number of year classes (Sukumaran and Neelakantan, 1997). High fishing pressure on crab stocks, unsuitable habitats, and climate change year can reduce a year’s swimming crab catch (Meynecke et al., 2012a,b; Chandrapavan et al., 2019). Although our study only focused on the season in which swimming crab production is highest, we expect that negative or positive climate change events in previous seasons would influence the swimming crab CRs in the following season, which should be further studied.
Conclusion
Catch trends for the three crucial commercial swimming crab species available in different climate events from 2013 to 2019 during the major fishing season in the TS revealed different patterns and variations with drastic decreases or increases, indicating a notable impact of climate change and regime shifts. The different patterns among the three swimming crab species strongly suggest that each species responds differently to climate change. We discovered that La Niña events were positively correlated with higher CRs during autumn for C. feriatus and P. pelagicus, whereas higher CRs were observed for P. sanguinolentus in both La Niña and El Niño years than in normal years. Moreover, different ENSO events influenced the habitat suitability of C. feriatus, P. pelagicus, and P. sanguinolentus in the TS. The lower CRs for C. feriatus and P. pelagicus in normal and El Niño years and for P. sanguinolentus in normal years during autumn were highly consistent with substantial shrinkage of suitable and optimal habitats.
These findings will enable the implementation of sustainable crab fisheries in the future through harvest strategy planning, ecosystem management, and spatiotemporal management. We suggest that other environmental factors – such as bottom salinity, sediment type, and organic matter content – should be added in future modeling to improve predictions because the natural habitat of crabs is typically the bottom of the sea, and a species’ habitat considerably affects its CRs. Moreover, in future studies, scholars should analyze the variables in other seasons and record the carapace size and sex of swimming crabs during the study period to more comprehensively investigate the influence of environmental factors on crab habitats and distributions.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
MN and K-WL designed the study. M-AL provided the datasets. MN and Y-LW processed and analyzed the datasets and worked on interpretations of the results with K-WL. MN wrote the manuscript. K-WL and M-AL reviewed the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was financially supported by the Council of Agriculture (107AS-11.2.1-FA-F1 and 108AS-9.2.1-FA-F1) and the National Science Council (MOST 108-2611-M-019-007 and MOST 109-2611-M-019-017).
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.
Acknowledgments
We are grateful to the Taiwan Ocean Conservation and Fisheries Sustainability Foundation for providing data from Taiwanese crab vessels. Wallace Academic Editing edited this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2021.763543/full#supplementary-material
Footnotes
- ^ https://origin.cpc.ncep.noaa.gov
- ^ https://marine.copernicus.eu
- ^ https://oceancolor.gsfc.nasa.gov
References
Abelló, P., and Hispano, C. (2006). The capture of the indo-Pacific crab Charybdis feriata (Linnaeus, 1758) (Brachyura: Portunidae) in the Mediterranean Sea. Aquat. Invasions 1, 13–16. doi: 10.3391/ai.2006.1.1.4
Azra, M. N., Aaqillah−Amr, M. A., Ikhwanuddin, M., Ma, H., Waiho, K., Ostrensky, A., et al. (2020). Effects of climate−induced water temperature changes on the life history of brachyuran crabs. Rev. Aquac. 12, 1211–1216. doi: 10.1111/raq.12380
Batoy, C. B., Sarmago, J. F., and Pilapil, B. C. (1987). Breeding season, sexual maturity and fecundity of the blue crab, Portunus pelagicus (L.) in selected coastal waters in Leyte and vicinity, Philippines. Ann. Trop. Res. 9, 157–177.
Campbell, G. R. (1986). Size at sexual maturity and occurrence of ovigerous females in three species of commercially exploited portunid crabs in SE Queensland. Proc. R. Soc. Qld. 97, 79–87.
Carpenter, K. E., Krupp, F., and Jones, D. A. (1997). Living Marine Resources of Kuwait, Eastern Saudi Arabia, Bahrain, Qatar, and the United Arab Emirates. Rome: Food & Agriculture Org.
Carrascal, L. M., Galván, I., and Gordo, O. (2009). Partial least squares regression as an alternative to current regression methods used in ecology. Oikos 118, 681–690. doi: 10.1111/j.1600-0706.2008.16881.x
Chandrapavan, A., Caputi, N., and Kangas, M. I. (2019). The decline and recovery of a crab population from an extreme marine heatwave and a changing climate. Front. Mar. Sci. 6:510. doi: 10.3389/fmars.2019.00510
Chang, Y. J., Lan, K. W., Walsh, W. A., Hsu, J., and Hsieh, C. H. (2019). Modelling the impacts of environmental variation on habitat suitability for Pacific saury in the Northwestern Pacific Ocean. Fish. Oceanogr. 28, 291–304. doi: 10.1111/fog.12408
Charmantier, G., Haond, C., Lignot, J., and Charmantier-Daures, M. (2001). Ecophysiological adaptation to salinity throughout a life cycle: a review in homarid lobsters. J. Exp. Biol. 204, 967–977. doi: 10.1242/jeb.204.5.967
Chen, T. Y., Hwang, G. W., Mayfield, A. B., Chen, C. P., and Lin, H. J. (2019). The development of habitat suitability models for fiddler crabs residing in subtropical tidal flats. Ocean Coast. Manag. 182:104931. doi: 10.1016/j.ocecoaman.2019.104931
Cheung, W. W. L., Lam, V. W. Y., Sarmiento, J. L., Kearney, K., Watson, R. E. G., Zeller, D., et al. (2009). Large−scale redistribution of maximum fisheries catch potential in the global ocean under climate change. Glob. Chang. Biol. 16, 24–35. doi: 10.1111/j.1365-2486.2009.01995.x
De Anda-Montañez, J. A., Martínez-Aguilar, S., Balart, E. F., Zenteno-Savín, T., Méndez-Rodríguez, L., Amador-Silva, E., et al. (2016). Spatio-temporal distribution and abundance patterns of red crab Pleuroncodes planipes related to ocean temperature from the Pacific coast of the Baja California Peninsula. Fish. Sci. 82, 1–15. doi: 10.1007/s12562-015-0938-8
de Andrade, L. S., Frameschi, I. F., Castilho, A. L., da Costa, R. C., and Fransozo, A. (2015). Can the pattern of juvenile recruitment and population structure of the speckled swimming crab Arenaeus cribrarius (Decapoda: Brachyura) be determined by geographical variations? Mar. Ecol. 36, 950–958. doi: 10.1111/maec.12188
de Lestang, S., Bellchambers, L. M., Caputi, N., Thomson, A. W., Pember, M. B., Johnston, D. J., et al. (2010). “Stock-recruitment-environment relationship in a Portunus pelagicus fishery in Western Australia,” in Biology and Management of Exploited Crab Populations under Climate Change, eds G. H. Kruse, G. L. Eckert, R. J. Foy, R. N. Lipcius, B. SainteMarie, D. L. Stram, et al. (Fairbanks, AK: Alaska Sea Grant, University of Alaska Fairbanks), 317–334.
Defeo, O., and Cardoso, R. S. (2002). Macroecology of population dynamics and life history traits of the mole crab Emerita brasiliensis in Atlantic sandy beaches of South America. Mar. Ecol. Prog. Ser. 239, 169–179. doi: 10.3354/meps239169
Dvoretsky, A. G., and Dvoretsky, V. G. (2020). Effects of environmental factors on the abundance, biomass, and individual weight of juvenile red king crabs in the Barents Sea. Front. Mar. Sci. 7:726. doi: 10.3389/fmars.2020.00726
Erauskin−Extramiana, M., Arrizabalaga, H., Hobday, A. J., Cabré, A., Ibaibarriaga, L., Arregui, I., et al. (2019). Large−scale distribution of tuna species in a warming ocean. Glob. Chang. Biol. 25, 2043–2060. doi: 10.1111/gcb.14630
Frid, A., McGreer, M., and Stevenson, A. (2016). Rapid recovery of Dungeness crab within spatial fishery closures declared under indigenous law in British Columbia. Glob. Ecol. Conserv. 6, 48–57. doi: 10.1016/j.gecco.2016.01.002
Gibson, J. R., and Najjar, R. G. (2000). The response of Chesapeake Bay salinity to climate−induced changes in streamflow. Limnol. Oceanogr. 45, 1764–1772. doi: 10.4319/lo.2000.45.8.1764
Green, B. S., Gardner, C., Hochmuth, J. D., and Linnane, A. (2014). Environmental effects on fished lobsters and crabs. Rev. Fish Biol. Fish. 24, 613–638. doi: 10.1007/s11160-014-9350-1
Ho, C. H., Chen, J. L., Nobuyuki, Y., Lur, H. S., and Lu, H. J. (2016). Mitigating uncertainty and enhancing resilience to climate change in the fisheries sector in Taiwan: Policy implications for food security. Ocean Coast. Manag. 130, 355–372. doi: 10.1016/j.ocecoaman.2016.06.020
Hosseini, M., Vazirizade, A., Parsa, Y., and Mansori, A. (2012). Sex ratio, size distribution and seasonal abundance of blue swimming crab, Portunus pelagicus (Linnaeus, 1758) in Persian Gulf Coasts, Iran. World Appl. Sci. J. 17, 919–925.
Hsiao, P. Y., Shimada, T., Lan, K. W., Lee, M. A., and Liao, C. H. (2021). Assessing summertime primary production required in changed marine environments in upwelling ecosystems around the Taiwan Bank. Remote Sens. 13:765. doi: 10.3390/rs13040765
Hsueh, P. W., and Hung, H. T. (2009). Temporal and spatial reproductive patterns of subtidal brachyuran crabs in coastal waters of Taiwan. Crustaceana 82, 449–465. doi: 10.1163/156854009x404743
Hu, J., and Wang, X. H. (2016). Progress on upwelling studies in the China seas. Rev. Geophys. 54, 653–673. doi: 10.1002/2015rg000505
Huang, T. H., Chen, C. T. A., Zhang, W. Z., and Zhuang, X. F. (2015). Varying intensity of Kuroshio intrusion into Southeast Taiwan Strait during ENSO events. Cont. Shelf Res. 103, 79–87. doi: 10.1016/j.csr.2015.04.021
Ikhwanuddin, M., Azra, M. N., Talpur, M. A., Abol-Munafi, A. B., and Shabdin, M. L. (2012). Optimal water temperature and salinity for production of blue swimming crab, Portunus pelagicus 1st day juvenile crab. Aquac. Aquar. Conserv. Legis. 5, 4–8.
Johnston, D., Harris, D., Caputi, N., and Thomson, A. (2011). Decline of a blue swimmer crab (Portunus pelagicus) fishery in Western Australia—History, contributing factors and future management strategy. Fish. Res. 109, 119–130. doi: 10.1016/j.fishres.2011.01.027
Kruse, G. H., Zheng, J., and Stram, D. L. (2010). Recovery of the Bristol Bay stock of red king crabs under a rebuilding plan. ICES J. Mar. Sci. 67, 1866–1874. doi: 10.1093/icesjms/fsq136
Kuczynski, L., Legendre, P., and Grenouillet, G. (2018). Concomitant impacts of climate change, fragmentation and non−native species have led to reorganization of fish communities since the 1980s. Glob. Ecol. Biogeogr. 27, 213–222. doi: 10.1111/geb.12690
Kuhn, A. A., and Darnell, M. Z. (2019). Elevated temperature induces a decrease in intermolt period and growth per molt in the lesser blue crab Callinectes similis Williams, 1966 (Decapoda: Brachyura: Portunidae). J. Crustac. Biol. 39, 22–27. doi: 10.1093/jcbiol/ruy089
Kunsook, C., Gajaseni, N., and Paphavasit, N. (2014). A stock assessment of the blue swimming crab Portunus pelagicus (Linnaeus, 1758) for sustainable management in Kung Krabaen Bay, Gulf of Thailand. Trop. Life Sci. Res. 25, 41–59.
Kuo, N. J., and Ho, C. R. (2004). ENSO effect on the sea surface wind and sea surface temperature in the Taiwan Strait. Geophys. Res. Lett. 31:L13309. doi: 10.1029/2004GL020303
Lan, K. W., Lee, M. A., Zhang, C. I., Wang, P. Y., Wu, L. J., and Lee, K. T. (2014). Effects of climate variability and climate change on the fishing conditions for grey mullet (Mugil cephalus L.) in the Taiwan Strait. Clim. Change 126, 189–202. doi: 10.1007/s10584-014-1208-y
Lau, N. C., and Nath, M. J. (2006). ENSO modulation of the interannual and intraseasonal variability of the East Asian monsoon—A model study. J. Clim. 19, 4508–4530. doi: 10.1175/JCLI3878.1
Leffler, C. W. (1972). Some effects of temperature on the growth and metabolic rate of juvenile blue crabs, Callinectes sapidus, in the laboratory. Mar. Biol. 14, 104–110. doi: 10.1007/bf00373209
Lehodey, P., Alheit, J., Barange, M., Baumgartner, T., Beaugrand, G., Drinkwater, K., et al. (2006). Climate variability, fish, and fisheries. J. Clim. 19, 5009–5030. doi: 10.1175/JCLI3898.1
Lehodey, P., Bertrand, A., Hobday, A. J., Kiyofuji, H., McClatchie, S., Menkès, C. E., et al. (2020). “ENSO impact on marine fisheries and ecosystems,” in El Niño Southern Oscillation in a Changing Climate, eds M. J. McPhaden, A. Santoso, and W. Cai (Hoboken, NJ: John Wiley & Sons, Inc), 429–451. doi: 10.1002/9781119548164.ch19
Liu, Q., Jia, Y., Liu, P., Wang, Q., and Chu, P. C. (2001). Seasonal and intraseasonal thermocline variability in the central South China Sea. Geophys. Res. Lett. 28, 4467–4470. doi: 10.1029/2001GL013185
Luan, J., Zhang, C., Xu, B., Xue, Y., and Ren, Y. (2018). Modelling the spatial distribution of three Portunidae crabs in Haizhou Bay, China. PLoS One 13:e0207457. doi: 10.1371/journal.pone.0207457
McPhaden, M. J., Zebiak, S. E., and Glantz, M. H. (2006). ENSO as an integrating concept in earth science. Science 314, 1740–1745. doi: 10.1126/science.1132588
Meynecke, J. O., Grubert, M., Arthur, J. M., Boston, R., and Lee, S. Y. (2012a). The influence of the La Niña-El Niño cycle on giant mud crab (Scylla serrata) catches in Northern Australia. Estuar. Coast. Shelf Sci. 100, 93–101. doi: 10.1016/j.ecss.2012.01.001
Meynecke, J. O., Grubert, M., and Gillson, J. (2012b). Giant mud crab (Scylla serrata) catches and climate drivers in Australia–a large scale comparison. Mar. Freshw. Res. 63, 84–94. doi: 10.1071/MF11149
Morris, A. S., Wilson, S. M., Dever, E. F., and Chambers, R. M. (2011). A test of bycatch reduction devices on commercial crab pots in a tidal marsh creek in Virginia. Estuaries Coast. 34, 386–390. doi: 10.1007/s12237-010-9330-1
Naimullah, M., Lan, K. W., Liao, C. H., Hsiao, P. Y., Liang, Y. R., and Chiu, T. C. (2020). Association of environmental factors in the Taiwan Strait with distributions and habitat characteristics of three swimming crabs. Remote Sens. 12:2231. doi: 10.3390/rs12142231
Ng, P. K., Wang, C. H., Ho, P. H., and Shih, H. T. (2001). An annotated checklist of brachyuran crabs from Taiwan (Crustacea: Decapoda). Natl. Taiwan Mus. 11, 1–86. doi: 10.11646/zootaxa.4502.1.1
Olsen, L., Herrmann, B., Sistiaga, M., and Grimaldo, E. (2019). Effect of gear soak time on size selection in the snow crab pot fishery. Fish. Res. 214, 157–165. doi: 10.1016/j.fishres.2019.02.005
Olson, A. P., Siddon, C. E., and Eckert, G. L. (2018). Spatial variability in size at maturity of golden king crab (Lithodes aequispinus) and implications for fisheries management. R. Soc. Open Sci. 5:171802. doi: 10.1098/rsos.171802
Polity, I. A., Muhamad, J. H., Long, S. M., and Bolong, A. M. A. (2011). Fecundity of blue swimming crab, Portunus pelagicus Linnaeus, 1758 from Sematan fishing district, Sarawak coastal water of South China Sea. Borneo J. Resour. Sci. Technol. 1, 46–51. doi: 10.33736/bjrst.262.2011
Potter, I. C., and De Lestang, S. (2000). Biology of the blue swimmer crab Portunus pelagicus in Leschenault Estuary and Koombana Bay, south-western Australia. J. R. Soc. West. Aust. 83, 443–458.
Rasheed, S., and Mustaquim, J. (2010). Size at sexual maturity, breeding season and fecundity of three-spot swimming crab Portunus sanguinolentus (Herbst, 1783)(Decapoda, Brachyura, Portunidae) occurring in the coastal waters of Karachi, Pakistan. Fish. Res. 103, 56–62. doi: 10.1016/j.fishres.2010.02.002
Ryer, C. H., Ottmar, M., Spencer, M., Anderson, J. D., and Cooper, D. (2016). Temperature-dependent growth of early juvenile southern Tanner crab Chionoecetes bairdi: implications for cold pool effects and climate change in the southeastern Bering Sea. J. Shellfish Res. 35, 259–267. doi: 10.2983/035.035.0128
Scheffers, B. R., De Meester, L., Bridge, T. C., Hoffmann, A. A., Pandolfi, J. M., Corlett, R. T., et al. (2016). The broad footprint of climate change from genes to biomes to people. Science 354:aaf7671. doi: 10.1126/science.aaf7671
Signa, G., Cartes, J. E., Solé, M., Serrano, A., and Sánchez, F. (2008). Trophic ecology of the swimming crab Polybius henslowii Leach, 1820 in Galician and Cantabrian Seas: Influences of natural variability and the Prestige oil spill. Cont. Shelf Res. 28, 2659–2667. doi: 10.1016/j.csr.2008.08.008
Silva, C., Andrade, I., Yáñez, E., Hormazabal, S., Barbieri, M. Á, Aranis, A., et al. (2016). Predicting habitat suitability and geographic distribution of anchovy (Engraulis ringens) due to climate change in the coastal areas off Chile. Prog. Oceanogr. 146, 159–174. doi: 10.1016/j.pocean.2016.06.006
Song, H. T., Ding, Y. P., and Xu, Y. J. (1989). Population component characteristics and migration distribution of the Portunus trituberculatus in the coast water of Zhejiang Province. Mar. Sci. Bull. 8, 66–74.
Soundarapandian, P., Varadharajan, D., Ilavarasan, N., Kumar, J., and Kumar, A. (2013). Mating behaviour of flower crab, Charybdis feriata (Linnaeus). J. Mar. Sci. Res. Dev. 3:127. doi: 10.4172/2155-9910.1000127
Spencer, D. M., Doubell, M. J., Brown, I. W., Rodriguez, A. R., Lee, S. Y., and Lemckert, C. J. (2019). Environmental indices for spanner crab (Ranina ranina) catch rates depend on regional oceanographic features. Estuar. Coast. Shelf Sci. 228:106361. doi: 10.1016/j.ecss.2019.106361
Srisurichan, S., Caputi, N., and Cross, J. (2005). Impact of lunar cycle and swell on the daily catch rate of western rock lobster (Panulirus cygnus) using time series modelling. N. Z. J. Mar. Freshwater Res. 39, 749–764. doi: 10.1080/00288330.2005.9517350
Sukumaran, K. K., and Neelakantan, B. (1997). Age and growth in two marine portunid crabs, Portunus (Portunus) sanguinolentus (Herbst) and Portunus (Portunus) pelagicus (Linnaeus) along the southwest coast of India. Indian J. Fish. 44, 111–131.
Sumpton, W. D., Smith, G. S., and Potter, M. A. (1989). Notes on the biology of the portunid crab, Portunus sanguinolentus (Herbst), in subtropical Queensland waters. Mar. Freshw. Res. 40, 711–717. doi: 10.1071/MF9890711
Szuwalski, C., Cheng, W., Foy, R., Hermann, A. J., Hollowed, A., Holsman, K., et al. (2021). Climate change and the future productivity and distribution of crab in the Bering Sea. ICES J. Mar. Sci. 78, 502–515. doi: 10.1093/icesjms/fsaa140
Szuwalski, C. S., and Punt, A. E. (2013). Fisheries management for regime-based ecosystems: a management strategy evaluation for the snow crab fishery in the eastern Bering Sea. ICES J. Mar. Sci. 70, 955–967. doi: 10.1093/icesjms/fss182
Waiho, K., Fazhan, H., and Ikhwanuddin, M. (2016). Size distribution, length–weight relationship and size at the onset of sexual maturity of the orange mud crab, Scylla olivacea, in Malaysian waters. Mar. Biol. Res. 12, 726–738. doi: 10.1080/17451000.2016.1200726
Wang, C., Deser, C., Yu, J. Y., DiNezio, P., and Clement, A. (2017). “El Niño and southern oscillation (ENSO): a review,” in Coral Reefs of the Eastern Tropical Pacific, eds P. W. Glynn, D. P. Manzello, and I. C. Enochs (Dordrecht: Springer Netherlands), 85–106.
Wold, S., Sjöström, M., and Eriksson, L. (2001). PLS-regression: a basic tool of chemometrics. Chemometr. Intell. Lab. Syst. 58, 109–130. doi: 10.1016/s0169-7439(01)00155-1
Wu, X. G., Yu, Z. Y., Cheng, Y. X., He, S. S., Yang, X. Z., Lu, J. F., et al. (2007). Effect of four groups of live feeds on larval development, growth (from Z_4 to Megalopa) and fatty acid composition of Eriocheir sinensis. J. Fish. Sci. China 14, 911–918.
Xie, Z., Liu, H., and Feng, L. (2002). “Culture technology of swimming crab, Portunus trituberculatus,” in The Breeding and Culture of Economic Marine Crabs, ed. Z. Y. Lin (Beijing: Chinese Agriculture Press). 1–98.
Yang, C. P., Li, H. X., Li, L., Xu, J., and Yan, Y. (2014). Population Structure, morphometric analysis and reproductive biology of Portunus sanguinolentus (Decapoda: Brachyura: Portunidae) in Honghai Bay, South China Sea. J. Crustac. Biol. 34, 722–730. doi: 10.1163/1937240X-00002273
Keywords: crustacean, Portunidae, climate change, habitat modeling, crab fishery
Citation: Naimullah M, Wu Y-L, Lee M-A and Lan K-W (2021) Effect of the El Niño–Southern Oscillation (ENSO) Cycle on the Catches and Habitat Patterns of Three Swimming Crabs in the Taiwan Strait. Front. Mar. Sci. 8:763543. doi: 10.3389/fmars.2021.763543
Received: 24 August 2021; Accepted: 27 September 2021;
Published: 12 October 2021.
Edited by:
Jose Luis Iriarte, Austral University of Chile, ChileReviewed by:
Alexander G. Dvoretsky, Murmansk Marine Biological Institute, RussiaKelsey Shea Swieca, Oregon State University, United States
Copyright © 2021 Naimullah, Wu, Lee and Lan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kuo-Wei Lan, [email protected]