Chunling Zhang, Danyang Wang, Zhenfeng Wang. Fishery analysis using gradient-dependent optimal interpolation[J]. Acta Oceanologica Sinica, 2022, 41(2): 116-126. doi: 10.1007/s13131-021-1895-y
Citation: Fengxia Zhou, Shuangling Wang, Han Fang, Jiani He, Liang Ye, Zhaohai Ding, Cuiting Li, Fajin Chen, Xuan Lu, Chunqing Chen, Yafei Meng. Distribution and sources of sedimentary organic matter in different aquaculture areas of northeastern Zhanjiang Bay using stable carbon and nitrogen isotopes[J]. Acta Oceanologica Sinica, 2024, 43(6): 38-48. doi: 10.1007/s13131-024-2377-9

Distribution and sources of sedimentary organic matter in different aquaculture areas of northeastern Zhanjiang Bay using stable carbon and nitrogen isotopes

doi: 10.1007/s13131-024-2377-9
Funds:  The National Natural Science Foundation of China under contract No. 42276047; the Guangdong Basic and Applied Basic Research Foundation under contract Nos 2023A1515110473 and 2021A1515110172; the Program for Scientific Research Start-up Funds of Guangdong Ocean University under contract No. R17058; the National College Student Innovation and Entrepreneurship Training Program Project under contract No. 202310566007.
More Information
  • Corresponding author: fjchen@gdou.edu.cn
  • Received Date: 2023-11-15
  • Accepted Date: 2024-05-07
  • Available Online: 2024-07-22
  • Publish Date: 2024-06-30
  • Zhanjiang Bay is a major aquaculture area in China with many types of mariculture products (such as oysters, fish, and shrimp). The culture area and shrimp output in Zhanjiang Bay are ranked first in China. We investigated the total organic carbon (TOC), total nitrogen (TN), TOC/TN ratio, and stable isotopes (δ13C and δ15N) of the fish and shrimp feed, fish and shrimp feces, and sedimentary organic matter (SOM) in and around different aquaculture areas of northeastern Zhanjiang Bay to study the impact of aquaculture activities on SOM. The average TOC contents of fish and shrimp feed were 39.20% ± 0.91% and 39.29% ± 0.21%, respectively. The average TOC content in the surface sediments of the oyster culture area, the mixed (fish and shrimp) culture area, and the cage fish farm area were 0.66%, 0.88% ± 0.10%, and 0.58% ± 0.19%, respectively, which may indicate that mixed culture had a greater impact on SOM. The relatively high TOC and TN contents and relatively low TOC/TN ratios, and δ15N values in the upper layer of the core sediment in the mixed culture area could also support the significant influence of mixed culture. The average δ13C and δ15N values of fish and shrimp feed were −20.6‰ ± 2.2‰ and 1.8‰ ± 1.2‰, respectively, which were different from the isotopic values of SOM in the study area. δ13C and δ15N values for SOM in different aquaculture areas were different from those of nearby reference stations, probably reflecting the influence of aquaculture. The δ13C and δ15N values in the oyster culture area (−25.9‰ and 6.0‰, respectively) seemed to have reduced δ13C and enriched δ15N relative to those of the reference station (−24.6‰ and 5.8‰, respectively). This may reflect the influence of organic matter on oyster culture. The δ15N value of the station in the mixed culture area (7.1‰ ± 0.4‰) seemed to be relatively enriched in δ15N relative to that of the reference station (6.6‰). Sedimentation and the subsequent degradation of organic matter from mixed cultures may have contributed to this phenomenon. The surface sediment at the cage fish farm area seemed to be affected by fish feces and primary production based on the indication of δ13C and δ15N values. The sediment core at the mixed culture region (NS6) had lower TOC/TN ratios and more positive δ13C and δ15N values than the sediment core at the oyster culture area, suggesting a higher proportionate contribution of marine organic matter in the mixed culture area. In summary, oyster culture, mixed culture, and cage fish culture in northeastern Zhanjiang Bay had a certain degree of impact on SOM, and mixed culture had more significant influences on SOM based on the high TOC contents and the significant vertical variations of TOC/TN ratio and δ15N value in the sediment of this area. This study provides new insights into the impact of aquaculture activities on SOM content.
  • Analyzing the environmental characteristics of fish activity is important in fishery oceanography. The basic properties of seawater, such as temperature and salinity, have important effects on fish behavior: different types of fish have different temperature and salinity requirements (Song and Wu, 2011; Lan et al., 2017). In the past decades, sea surface observation data were obtained by satellite sensors at various spatial and temporal scales, which is helpful for fishery research in the open ocean (Laurs and Fiedler, 1985; Fiúza, 1990; Kluger et al., 2016; Lan et al., 2017). However, many important processes and structures associated with fishery resource distribution occur below the surface or at greater depths. As a major component of the ocean observation system, the increasing Argo array has allowed, for the first time, a near-global sampling of the upper 2 000 m of ice-free oceans with relatively small bias since 2005 (Roemmich and Gilson, 2009; Johnson et al., 2015; Liu et al., 2020). The obtained subsurface environmental data containing multiple variables and can be used to better understand the relationships between oceanic environments and fishing conditions. However, the spatial–temporal mismatch between fishery and environment data means that the latter cannot be used effectively in fishery research. Presently, the integration of available environmental data sources into fishery data is a challenging but urgent task in fisheries oceanography.

    Tuna species are critical commercial open-ocean fish and are of great ecological importance. For example, the yellowfin tuna (YFT) can swim at high speeds and exhibit deep-diving behavior (diving to over 500 m) (Schaefer et al., 2007, 2009). They dive into deep depths for food and spend about 30% of the daytime below the thermocline. It is believed that the long-distance vertical movement of the YFT may weaken the correlation between the YFT catch rate and sea-surface temperature (Zagaglia et al., 2004; Schaefer et al., 2011, 2014). Many studies have indicated that the longline catch per unit effort of the YFT had a closer relationship with the subsurface water temperature rather than the surface water temperature. Therefore, the vertical temperature structure, especially the thermocline distribution, is a vital factor in tuna fisheries (Langley et al., 2009; Lan et al., 2011; Yang et al., 2015). Nevertheless, research on the relationship between oceanic interior characteristics and fishery activities is rare because of the lack of synchronous environmental observations with catch data. Some studies rely on inversion models, climate-gridded products, or finding the nearest synchronous observation (Song and Wu, 2011; Klemas and Yan, 2014; Yang et al., 2015). Subsurface data retrieved from satellites or climate products have the problem of spatial–temporal mismatch. When the distance is significant, the nearest synchronous observation may describe the environment imprecisely or unilaterally. It is necessary to develop an approach to solve this problem.

    Optimal interpolation (OI) is an objective analysis method proposed by Gandin (1963). OI equations can be easily derived as a variation problem to find a vector of estimates that minimizes the total error variance of the field being estimated (Wikle and Berliner, 2007). The estimator is unbiased and linear between the observations and the first estimate (a priori estimation) of the field. It has been widely used to calculate the gridded analysis data for ocean numerical models (Bonanno et al., 2019; Burgess and Webster, 1980; De Feis et al., 2020). Zhang et al. (2013) proposed the gradient-dependent correlation scale method, which uses the horizontal gradient change of oceanographic elements to estimate the correlation scales in the OI. This scheme, named gradient-dependent OI in the following, effectively extracted shorter wavelength information in zones with a large gradient (Zhang et al., 2015). It provides a new perspective for constructing real-time environmental information of fishing sites using a point-to-point method. The optimal weighted average is expected to provide results close to real observations. Then, we provide a data basis for better study of the fish living environment.

    In this respect, this study applies the gradient-dependent OI to obtain high-resolution subsurface data from the Argo array according to YFT fishery dates and locations in the western and central Pacific Ocean (WCPO). The results were used to determine the most suitable environmental conditions for YFT fishing. The fishery and Argo observations and the principle of the method are described in Section 2. In Section 3, we verify the results and apply them to validate the relationship between fish activities and the subsurface seawater environment. We discussed the suitable temperature and salinity vertical distribution for YFT in Section 4. Finally, in Section 5, we present our conclusions.

    Daily fishery data from the WCPO purse seiners were obtained between January 1, 2008 and August 31, 2017 from 14 ocean-going fishing vessels (Jinhui 1, 2, 3, 6, 7, 8, 9, 18, and 58, Pohnpei1, LOJET, LOMETO, LOMALO, and MAJURO1) of the Shanghai Kaichuang Fishery Co., Ltd. (SKFC). These data included the date, longitude, latitude, operation type (free-swimming and fish-aggregating device used), ship name, daily catch per ship, and other production information. The amount of free-swimming was small, and we selected the fishery data of the fish aggregating device (FAD) used only (Fig.1a). There were 6354 fishery points used in this study. These fishing vessels have similar fishing capacities. Therefore, this study ignored the influence of ship properties on fishery catch data. We adopted the daily output per ship as the catch as the measure of the resource abundance.

    Figure  1.  Location of YFT fishery points and the available Argo observations. a. The YFT fishing locations in the WCPO and the locations of fishing locations and Argo profiles during August 2017; b. the frequency statistic of observation number surrounding each fishing point.

    Two categories of environmental data were collected in this study: Argo profiles and WOA18 gridded data. The former were derived from the China Argo Real-Time-Data Center (ftp://ftp.argo.org.cn/pub/ARGO/global/). These include temperature, salinity, and pressure, which are processed in real-time and delayed-mode quality control (Tong et al., 2003; Liu et al., 2006). We collected Argo observations for the same period as the YFT catch data in the WCPO. A radial distance of 500 km was set to ensure that more than five Argo temperature and salinity profiles were obtained at every fishery point each day. Only profiles located within the specified range were used to establish subsurface environmental information for fishery points at depths between 5 m and 1500 m. A total of 62726 Argo profiles surround the 6354 fishery points. Figure 1a displays the locations of fishery points in August 2017 and their corresponding Argo profiles within the effective radius. The frequency statistics of the available Argo observations surrounding every fishery point are shown in Fig. 1b. All fishery points have more than 10 and sometimes as many as 35 environmental profiles within the specified range. Most fishery points have 15–23 available profiles. The World Ocean Atlas 2018 (WOA18) gridded data were used to calculate the horizontal gradients and provide the first guess field via bilinear interpolation at each fishery point. Two gridded datasets, longitude, latitude, pressure, temperature, and salinity were downloaded from the Ocean Climate Laboratory of NODC (https://www.nodc.noaa.gov/OC5/woa18/). These data have a pixel resolution of (1/4)° and span the area from 30°S to 30°N and 120°E to 140°W.

    The purpose of the scheme was to obtain an optimal analysis value for the goat point based on least squares theory. The analysis results at the fishery point was used as the background value plus the innovation value weighted by the optimal weights. As for the fishery point influenced by M observations, the standard analysis equation is given by Eq. (1). A critical part of the scheme is to estimate the optimal weights. According to the minimum variance theory, the optimal weight can be determined by solving Eq. (2) (Kalnay, 2003):

    $$ v_i^a = v_i^b + \sum\limits_j^M {{w_{ij}}\delta y_j^o} {\rm{,}} $$ (1)
    $$ \sum\limits_{j = 1}^M {{w_{ij}}} {\mu _{jk}} + {\eta _k}{w_{ik}} = {\mu _{ik}},\mathop {}\nolimits_{}^{} k = 1, \cdot \cdot \cdot ,M {\rm{,}} $$ (2)

    where $ v_i^a $ is the analysis value and $ v $ can be any environmental variable, such as temperature or salinity. The first estimate $ v_i^b $ was given by the WOA18 monthly data via bilinear interpolation. The subscript i denotes the number of fishery sampling points, and j and k denote the number of available Argo profile sites. For observational increments, $ \delta y_j^o{{ = }}y_j^o{{ - }}H(v_j^b) $, the observational operator H converted the background field into the first guesses of the observation $ y_j^o $. In this study, a radial distance was set to ensure that only the Argo profiles located within a specified range surrounding each fishery point were counted. Each increment has an optimal weight $ {w_{ik}} $ associated with the background error correlation, $ {\mu _{jk}} $ and $ {\mu _{ik}} $. where $ {\mu _{jk}} $, $ {\mu _{ik}} $ are the correlations between the background errors at the two observational points, j and k, and at the fishery and observation points i and k, respectively. $ \eta $ is the mean square of the relative observation errors compared with the background errors, which was frequently “tuned” to vary the weights of the observation. We used the value of $ \eta {{ = 0}}{{.25}} $ for the Pacific Ocean in this study (Zhang et al., 2015). The correlations are usually assumed to follow a Gaussian exponential function (Zhang et al., 2013):

    $$ {\mu _{ik}}\sim\exp \left[ { - \frac{{{{({x_i} - {x_k})}^2}}}{{{{({L_\phi }/{G_x})}^2}}} - \frac{{{{({y_i} - {y_k})}^2}}}{{{{({L_\phi }/{G_y})}^2}}}} \right] , $$ (3)
    $$ {G}_{x}=1+\frac{|\partial v/\partial x|}{E|\partial v/\partial x|}, $$ (4)
    $$ {G}_{y}=1+\frac{|\partial v/\partial y|}{E|\partial v/\partial y|} , $$ (5)

    where x and y are the longitude and latitude, respectively, and $ {L_\phi }/G $ depend on the Rossby radius of deformation and the horizontal gradient changes, which define the background error correlation scale. $ {L_\phi } $ is the correlation scale that can be obtained from the scale parameter multiplied by the cosine function of the latitude ($ \phi $) at fishery sampling point i. The radial distance was set to 500 km when constructing high-resolution profiles to ensure that a sufficient number of observations were considered. The parameter G, calculated using the WOA18 climatological data, is associated with the horizontal gradient (s) at location i. It contains a zonal component ($ G_x $) and a meridional component ($ G_y $). The Gaussian functions presented in Eqs (3–5) provide anisotropic background error correlations at each fishery point.

    According to the algorithm described in Eqs (1–5), we obtained 6354 temperature and salinity profiles based on Argo observations. Each profile was obtained on the same time (day) and location as the YFT fishery data from January 1, 2008 to August 31, 2017, in the WCPO. The vertical depths (5–1500 m) were divided into unequally spaced 24 layers.

    We randomly selected 30 Argo T/S profiles as the “true” observations and constructed the analysis profiles according to the Argo observation locations using the gradient-dependent OI shown in Section 2.2. The true observations at each location were not included when we constructed the analysis profiles using the Argo observation profiles located within the specified range. These analysis profiles were then used to verify the algorithm by comparing them with the corresponding “true” observations. The verification was carried by calculating the root mean square errors (RMSEs) and differences between analysis and “true” observations. The RMSEs of analysis results and those of nearest profiles at the depth range of 5 m to 1 500 m were shown in Fig. 2. Figure 3 displays the temperature and salinity differences between analysis results and “true” observations at depths of 5 m to 500 m.

    Figure  2.  Root mean square errors (RMSEs) of 30 sample points from the analysis results (indicated by dots) and nearest profiles (indicated by circles) at a depth range of 5 m to 1500 m.
    Figure  3.  Temperature and salinity differences between 5 m to 500 m at each fishery point. a and c show the vertical difference distribution, whereas b and d are the difference frequency at each depth.

    It can be seen that both temperature and salinity have larger RMSEs and differences at the upper depth than at the deep depth. The temperature differences were no more than 0.50°C and salinity differences were less than 0.01 above 500 m. Figures 3a and c indicate that the differences were near zero in the middle and deep layers. The maximum temperature RMSE is approximately 0.4°C. If the nearest temperature profiles are adopted, the corresponding RMSEs exceed 0.4°C at depths of 50 m to 300 m, which are the principal activity layers of YFT (Fig. 2a). As for the salinity RMSE (Fig. 2b), the difference was more significant. The salinity RMSEs of the analysis results present a value of less than 0.02 at a depth of 5–1500 m. At the nearest profiles, the salinity RMSEs are all larger than 0.04 above 300 m.

    In other words, the analysis results were very close to the observations. Near the surface, the temperature and salinity differences were relatively significant. There were one-third profiles whose temperature differences are in 0.4–0.5°C at a depth of 5 m (Fig. 3b). At a depth of 10 m, 11 profiles in 30 had salinity differences of 0.006–0.01(Fig. 3d). The lack of surface observations in most Argo profiles made the true values scarce at 5 m, and the accuracy was lower. The salinity observations with more significant variances because of the salinity sensor drift at the near-surface would influence the accuracy of the analysis. However, most temperature and salinity differences at depths greater than 150 m were less than 0.2°C and 0.002, respectively. At depths deeper than 300 m, temperature RMSEs are less than 0.1°C, and salinity RMSEs are close to 0 for the analysis results. The nearest RMSEs are large, with a maximum RMSE of approximately 0.9°C for temperature and 0.17 for salinity, as shown in Fig. 2.

    The temperature and salinity profiles of 6354 constructed using Eqs (1–5) were applied to discuss the relationship between the YFT catch data and several oceanic environmental factors, including temperature, thermocline, and salinity, using frequency statistics and contrastive analysis. Most of the YFT catch values corresponding to these profiles were concentrated in 1–50 t. To give the general vertical characteristics clearly, we interpolated the temperature and salinity of each water layer to YFT catch values of 1−50 t. The temperature and salinity sections with the change in catch are displayed in Figs 4 and 5. The YFT suitable temperature and salinity were statistically significant at 5 m, 150 m, and 300 m, as shown in Fig. 6. Finally, the relationship between the thermocline parameters and YFT catch is shown in Fig. 7 and Table 1.

    Figure  4.  Temperature section at the depth range from 5 m to 100 m (a), 5 m to 500 m (b), and temperature curves at 5 m and 100 m (c) for different catch. The temperatures at every water layer were interpolated to the catch sequence. C. I indicates the contour intervals.
    Figure  5.  Salinity section at the depth range from 5 m to 100 m (a), 5 m to 500 m (b), and salinity curves at 5 m and 150 m (c) for different catch. The salinities at every water layer were interpolated to the catch sequence. C. I indicates the contour intervals.
    Figure  6.  The statistic of suitable temperature and salinity for YFT corresponding different catch group. Temperature at 5 m, 150 m, and 300 m are indicated by a, c, and e, respectively. Salinity at 5 m, 150 m and 300 m are indicated by b, d, and f, respectively.
    Figure  7.  Statistics for the thermocline variable matching each fishery point. a. The upper bound depth (indicated by h1) and the deeper bound depth (indicated by h2); b. the thermocline thickness (represented by thk) and thermocline strength represented by temperature gradient (indicated by gth).
    Table  1.  Thermocline parameters for different catch (t) groups
    <10 t 10−20 t 20−30 t >30 t
    h1/m 80−90 80−90 90−100 90−100
    h2/m 350−400 350−400 350−400 300−350
    thk/m 350−360 260−270 280−290 240−250
    gth/(°C·m−1) 0.14−0.15 0.14−0.15 0.12−0.13 0.12−0.13
    Note: The upper and deeper bound depth of the thermocline are represented by h1 and h2, respectively. Thk and gth indicate the thickness and gradient representing strength of the thermocline, respectively.
     | Show Table
    DownLoad: CSV

    Figure 4b shows the general characteristics of the temperature distribution corresponding to different catches. At the fishery ground, the near-surface temperature was approximately 30°C. The temperature remained high (>27°C) up to a depth of 100 m. At a depth of 300 m, the temperature decreased significantly to 12°C. There was an evident thermocline corresponding to each catch value. When the catches were larger than 25 t, the temperature isopleths tended to be smooth. In contrast, the temperature fluctuated greatly with the catch change when the catch was small. The vertical temperature characteristics at the upper mixed layer (<100 m) are clearly displayed in Fig. 4a. It can be seen that many fishery points with lower catches had surface temperature values of 30°C at depths lower than 50 m and had temperature values of 28°C at depths of 50–100 m. At the fishery points with higher catches, their temperature was concentrated in the range of 28.5–29.5°C at depths of 5–100 m. The temperature curves shown in Fig. 4c clearly demonstrate the similar trends: (1) when the catches were higher (>25 t), the temperatures tended to be stable. Otherwise, the temperatures had large fluctuations at the YFT fishery points with lower catches. (2) At the near-surface, the temperatures of the fisher points with higher catches were slightly colder than those with a lower catch. (3) At all the fishery points, the vertical temperature change had a value of approximately 17°C at a depth of 5–300 m.

    Figure 5 shows that the relationships between the YFT catch weights and seawater salinities had similar characteristics with temperature. Overall, the salinities presented high values of 35–36 at depths of 50–300 m (Fig. 5b). The near-surface salinity was approximately 34.5–35.5. Most salinities derived from higher catches (>25 t) were 34.5–35 at depths lower than 50 m. At the fishery points with lower catches, the near-surface salinities were concentrated in 34–35 (Fig. 5a) with salinity values near 35.5 at a depth of 150 m (Fig. 5c). When the catches were higher than 25 t, the salinity differences between the depths of 5 m and 150 m were less than 0.5. The vertical salinity change was smaller than that in the lower catch fishery points.

    According to their weights, the YFT catches (daily output per ship) were divided into two groups: higher catch (>25 t) and lower catch (<25 t). Three representative layers (5 m, 150 m, and 300 m) were selected, and five parts were divided at each depth according to the temperature and salinity scale. Figure 6 displays the statistical results for suitable temperature and salinity for the YFT.

    At the near-surface (<5 m), the suitable temperatures and salinity were 28–29°C and 34.5–35.0, respectively (see Figs 6a and b). The higher catches had smaller variations, and the lower catches had a large number of higher temperatures (> 30°C) and lower temperatures (<27°C). When the salinities were in the range of 35.0–35.5, the catches increased. However, the catches were lower with a salinity larger than 35.5 or smaller than 34.0. Suitable temperatures were 22–23°C at a depth of 150 m. This can be represented by both catch groups (Fig. 6c). Figure 6d shows that the suable salinities had a larger value of exceeding 35.6 at a depth of 150 m. When the salinities were between 35.4 and 35.6, the YFT were frequently active. At a depth of 300 m, the most suitable temperatures and salinities were 11–12°C and 34.5–35.0, respectively, just as shown in Figs 6e and f.

    Thermoclines can be clearly observed in the upper ocean where the YFT migrates mainly. The thermocline, or the vertical temperature gradient, is an important factor influencing the YFT activities. Figure 7 shows the characteristics of the thermocline parameters corresponding to different catchments. The depth and strength of the thermocline were calculated using the maximum angle method (Chu and Fan, 2011). Here, we considered only the seasonal thermocline that existed in the top layer. The thermocline parameters were obtained first according to the Argo temperature profiles surrounding each fishery point. These parameters were then merged with the fishery points using the method described in Section 2.2. The thermocline parameters mainly included the upper bound depth (indicated by h1 in Fig. 7), the deeper bound depth (indicated by h2 in Fig. 7), and the temperature gradient that indicates the thermocline strength (indicated in the gth panel in Fig. 7). The thermocline thickness is represented by symbol thk.

    Combined with Fig. 1, these fishery points were primarily located in the tropical western Pacific. The location of the thermocline was evident. Figure 7a shows that most of the upper bounds of the thermocline (h1) are 50–150 m. At points whose catches exceeded 25 t, the upper bound of the thermocline was approximately 100 m, which was deeper than that of the lower catch points. The deeper bound depth of the thermocline (h2) was between 250 m and 450 m. Most of the large catch points had a value of 300–400 m. The mean thermocline thickness was approximately 250 m at most fishery points. Fishery catches were lower when the thermocline thickness exceeded 300 m. As can be seen in Fig. 7b, the stronger the thermocline, the lower the catch. For the higher catch fishery points, most of them had temperature gradients (gth) of 0.12–0.15°C/m. For the fishery points with catches lower than 25 t, their corresponding temperature gradients (gth) were very large (>0.16°C/m) or very small (< 0.10°C/m). The temperature gradients of the thermocline were larger than 0.08°C/m at all of the fishery points.

    Based on their weight, YFT catches were divided into four groups: < 10 t, 10–20 t, 20–30 t, and > 30 t. In Table 1, matching the four catch clusters (< 10 t, 10–20 t, 20–30 t, and > 30 t), the thermocline parameter ranges with the largest proportion are given. The h1, h2, thk, and gth results were given by statistics with intervals of 10 m, 50 m, 10 m, and 0.01°C/m, respectively. The table shows that the suitable upper (h1) and deeper bound depths of the thermocline (h2) were 90–100 m and 300–350 m, respectively. The corresponding suitable thickness was 240–250 m. The larger the thickness and gradient, the smaller the catches became. The most suitable temperature gradients for YFT were 0.12–0.13°C/m.

    The YFT is a warm-water fish that needs a specific water temperature to inhabit and spawn. Studies have shown that the highest yield of YFT comes from areas where the surface temperature is approximately 29°C (Wang et al., 2014). The statistical results indicated that the most suitable temperatures and salinities for the YFT were 28–29°C at 5 m. Most fishery points with high catch values had temperature values of 29–29.5°C. When the near-surface temperatures were above 30°C or below 28°C, there were few catches. The vertical temperature decreased with increasing depth. The temperature was reduced to 22–23°C at depths of 150 m, where the YFT spent 5% of their time staying in the water layer. The temperature difference between the surface and 150 m was approximately 8°C. This is consistent with previous research results (Schaefer et al., 2007; Weng et al., 2009).

    Salinity is not as critical as temperature. The suitable near-surface salinities were found to be 34.5–35.0. The higher the catch, the smaller the change in the vertical salinity. These conclusions are similar to those obtained in previous studies (Song and Wu, 2011; Yang et al., 2015). Most of the time, YFT lived in the warm water layer below 100 m and migrated to food in the upper ocean. The tolerance time of YFT to environmental changes is limited. They were more suitable for less vertical temperature and salinity changes, as shown in this study.

    The YFT can break through the temperature thermocline and enter the deep water layers, even deeper than 1000 m (Schaefer et al., 2011). However, it was also pointed out that the purpose of YFT frequently entering deep water was foraging. Thermocline may influence the diurnal vertical distribution of the YFT by affecting the diurnal vertical distribution of food organisms. Most of the time, the YFT lived in water layers lower than 100 m but rarely below 10 m in the western and central Pacific Ocean (Leroy et al., 2007). Compared with sea surface environmental information, the thermocline had a more significant effect on YFT catch. At the same time, the YFT catches were small when the thickness and gradient of the thermocline were very large. At the moment of thermocline formation, the temperature changed dramatically in the vertical direction. This is not conducive to forming a good fishing ground. As shown in Table 1, the catches were lower when the temperature gradient exceeded 0.13°C/m. The central fishing ground was concentrated in an area with an upper bound depth of ~90 m. Moreover, tuna were visual and opportunistic predators, and light was weak in deep water, which reduced the ability of YFT to hunt. It is difficult to form a fishing ground where the deep boundary of the thermocline is deep.

    A new data assimilation algorithm called gradient-dependent OI was proposed in this study. The purpose of this study was to verify the effectiveness of gradient-dependent OI in constructing near-real-time (daily) subsurface (5–1500 m) environmental information for fishery data based on Argo profiles. We constructed temperature and salinity profiles corresponding to YFT dates and locations. These fishery points spanned 10 years, with more than 6000 locations in the WOPC. Most of the relative analysis errors were near the observational errors. Compared with the Argo observations, the derived RMSEs were less than 0.4°C for temperature and 0.02 for salinity. In places where the temperature and salinity changed rapidly, the analysis errors were relatively large. The most significant differences between analysis and observation were less than ±0.5°C for temperature and ±0.01 for salinity.

    The most suitable temperatures and salinities for YFT at the near-sea surface were 28–29°C and 34.5–35.0 at the near-sea surface. As the depth increased, the temperatures decreased to 22–23°C and 11–12°C at depths of 150 m and 300 m, respectively, and suitable salinities changed to higher than 35.6°C at a depth of 150 m. The upper and lower bound depths of the thermocline were 90–100 m and 300–350 m, respectively. The thermocline characteristics were prominent, with a mean temperature gradient exceeding 0.08°C/m. It is not surprising to see complex environmental effects on biological populations. Nevertheless, the simple application demonstrated that the method presented in this study is effective. The novel approach proposed in this paper can provide many types of near-real-time environmental information based on Argo observations. Other factors, including dissolved oxygen and chlorophyll, will be analyzed in future studies.

    The research is supported by the China National Scientific Seafloor Observatory. We are grateful to the SKFC of China for providing us with the data from the WCPO purse seine fishery, the China Argo Real-time Data Center for the near-real-time observations, and the National Oceanographic Data Center for providing us with the climate data. We thank Wiley Editing Services (www.wileyauthors.com/eeo/preparation) for editing this manuscript.

  • Anderson B, Zhang Li, Wang Huining, et al. 2017. Sedimentary carbon and nitrogen dynamics reveal impact of human land-use change on Kawainui Marsh, O‘ahu, Hawai‘i. Pacific Science, 71(1): 17–27, doi: 10.2984/71.1.2
    Bouchet V M P, Sauriau P G. 2008. Influence of oyster culture practices and environmental conditions on the ecological status of intertidal mudflats in the Pertuis Charentais (SW France): A multi-index approach. Marine Pollution Bulletin, 56(11): 1898–1912, doi: 10.1016/j.marpolbul.2008.07.010
    Bouillon S, Connolly R M, Lee S Y. 2008. Organic matter exchange and cycling in mangrove ecosystems: recent insights from stable isotope studies. Journal of Sea Research, 59(1-2): 44–58, doi: 10.1016/j.seares.2007.05.001
    Brandes J A, Devol A H. 2002. A global marine-fixed nitrogen isotopic budget: implications for Holocene nitrogen cycling. Global Biogeochemical Cycles, 16(4): 1120
    Chen Fajin, Huang Chao, Lao Qibin, et al. 2021a. Typhoon control of precipitation dual isotopes in southern China and its palaeoenvironmental implications. Journal of Geophysical Research: Atmospheres, 126(14): e2020JD034336, doi: 10.1029/2020JD034336
    Chen Fajin, Lu Xuan, Song Zhiguang, et al. 2021b. Coastal currents regulate the distribution of the particulate organic matter in western Guangdong offshore waters as evidenced by carbon and nitrogen isotopes. Marine Pollution Bulletin, 172: 112856, doi: 10.1016/j.marpolbul.2021.112856
    Cifuentes L A, Sharp J H, Fogel M L. 1988. Stable carbon and nitrogen isotope biogeochemistry in the Delaware Estuary. Limnology and Oceanography, 33(5): 1102–1115, doi: 10.4319/lo.1988.33.5.1102
    Cole M L, Valiela I, Kroeger K D, et al. 2004. Assessment of a δ15N isotopic method to indicate anthropogenic eutrophication in aquatic ecosystems. Journal of Environmental Quality, 33(1): 124–132, doi: 10.2134/jeq2004.1240
    Das S K, Routh J, Roychoudhury A N, et al. 2008. Elemental (C, N, H and P) and stable isotope (δ15N and δ13C) signatures in sediments from Zeekoevlei, South Africa: A record of human intervention in the lake. Journal of Paleolimnology, 39(3): 349–360, doi: 10.1007/s10933-007-9110-5
    Franco-Nava M A, Blancheton J P, Deviller G, et al. 2004. Particulate matter dynamics and transformations in a recirculating aquaculture system: Application of stable isotope tracers in seabass rearing. Aquacultural Engineering, 31(3-4): 135–155, doi: 10.1016/j.aquaeng.2004.01.003
    He Biyan, Dai Minhan, Huang W, et al. 2010. Sources and accumulation of organic carbon in the Pearl River Estuary surface sediment as indicated by elemental, stable carbon isotopic, and carbohydrate compositions. Biogeosciences, 7(10): 3343–3362, doi: 10.5194/bg-7-3343-2010
    Jiang Zengjie, Fang Jianguang, Mao Yuze, et al. 2009. Assessment on the sediment quality condition of cage culture area in Nansha Bay. Environmental Science and Management (in Chinese), 34(6): 159–163
    Jiang Zengjie, Fang Jianguang, Mao Yuze, et al. 2012. Identification of aquaculture-derived organic matter in the sediment associated with coastal fish farming. Journal of Fishery Sciences of China (in Chinese), 19(2): 348–354
    Kaiser D, Unger D, Qiu Guanglong. 2014. Particulate organic matter dynamics in coastal systems of the northern Beibu Gulf. Continental Shelf Research, 82: 99–118, doi: 10.1016/j.csr.2014.04.006
    Ke Zhixin, Chen Danting, Liu Jiaxing, et al. 2020. The effects of anthropogenic nutrient inputs on stable carbon and nitrogen isotopes in suspended particulate organic matter in Jiaozhou Bay, China. Continental Shelf Research, 208: 104244, doi: 10.1016/j.csr.2020.104244
    Kendall C. 1998. Tracing nitrogen sources and cycling in catchments. In: Kendall C, McDonnell J J, eds. Isotope Tracers in Catchment Hydrology. Amsterdam: Elsevier Science, 519–576
    Krishnamurthy R V, Bhattacharya S K, Kusumgar S. 1986. Palaeoclimatic changes deduced from 13C/12C and C/N ratios of Karewa Lake sediments, India. Nature, 323(6084): 150–152, doi: 10.1038/323150a0
    La Rosa T, Mirto S, Favaloro E, et al. 2002. Impact on the water column biogeochemistry of a Mediterranean mussel and fish farm. Water Research, 36(3): 713–721, doi: 10.1016/S0043-1354(01)00274-3
    Lamb A L, Wilson G P, Leng M J. 2006. A review of coastal palaeoclimate and relative sea-level reconstructions using δ13C and C/N ratios in organic material. Earth-Science Reviews, 75(1–4): 29–57, doi: 10.1016/j.earscirev.2005.10.003
    Lao Qibin, Chen Fajin, Jin Guangzhe, et al. 2023a. Characteristics and mechanisms of typhoon-induced decomposition of organic matter and its implication for climate change. Journal of Geophysical Research: Biogeosciences, 128(6): e2023JG007518, doi: 10.1029/2023JG007518
    Lao Qibin, Liu Sihai, Ling Zheng, et al. 2023b. External dynamic mechanisms controlling the periodic offshore blooms in Beibu Gulf. Journal of Geophysical Research: Oceans, 128(6): e2023JC019689, doi: 10.1029/2023JC019689
    Lao Qibin, Lu Xuan, Chen Fajin, et al. 2023c. A comparative study on source of water masses and nutrient supply in Zhanjiang Bay during the normal summer, rainstorm, and typhoon periods: Insights from dual water isotopes. Science of the Total Environment, 903: 166853, doi: 10.1016/j.scitotenv.2023.166853
    Lao Qibin, Wu Junhui, Chen Fajin, et al. 2022a. Increasing intrusion of high salinity water alters the mariculture activities in Zhanjiang Bay during the past two decades identified by dual water isotopes. Journal of Environmental Management, 320: 115815, doi: 10.1016/j.jenvman.2022.115815
    Lao Qibin, Zhang Shuwen, Li Zhiyang, et al. 2022b. Quantification of the seasonal intrusion of water masses and their impact on nutrients in the Beibu Gulf using dual water isotopes. Journal of Geophysical Research: Oceans, 127(7): e2021JC018065, doi: 10.1029/2021JC018065
    LaZerte B D. 1983. Stable carbon isotope ratios: implications for the source of sediment carbon and for phytoplankton carbon assimilation in Lake Memphremagog Quebec. Canadian Journal of Fisheries and Aquatic Sciences, 40(10): 1658–1666, doi: 10.1139/f83-192
    Li Jiacheng, Cao Ruixue, Lao Qibin, et al. 2020. Assessing seasonal nitrate contamination by nitrate dual isotopes in a monsoon-controlled bay with intensive human activities in South China. International Journal of Environmental Research and Public Health, 17(6): 1921, doi: 10.3390/ijerph17061921
    Liao Weisen, Hu Jianfang, Zhou Haoda, et al. 2018. Sources and distribution of sedimentary organic matter in the Beibu Gulf, China: Application of multiple proxies. Marine Chemistry, 206: 74–83, doi: 10.1016/j.marchem.2018.09.006
    Liu Jingjing. 2019. Sources of sedimentary organic matter in and around cage fish farming of Poyang Lake using stable carbon and nitrogen isotopes (in Chinese) [dissertation]. Nanchang: Nanchang University
    Liu Guofeng, Xu Pao, Wu Ting, et al. 2018. Present condition of aquaculture nitrogen and phosphorus environmental pollution and future development strategy. Jiangsu Journal of Agricultural Sciences (in Chinese), 34(1): 225–233
    Liu Sai, Yang Qian, Yang Shu, et al. 2014. The long-term records of carbon burial fluxes in sediment cores of culture zones from Sanggou Bay. Haiyang Xuebao (in Chinese), 36(8): 30–38
    Lu Xuan, Zhou Fengxia, Chen Fajin, et al. 2020. Spatial and seasonal variations of sedimentary organic matter in a subtropical bay: Implication for human interventions. International Journal of Environmental Research and Public Health, 17(4): 1362, doi: 10.3390/ijerph17041362
    Lu Xuan, Zhou Xin, Jin Guangzhe, et al. 2022. Biological impact of Typhoon Wipha in the coastal area of western Guangdong: A comparative field observation perspective. Journal of Geophysical Research: Biogeosciences, 127(2): e2021JG006589, doi: 10.1029/2021JG006589
    McGhie T K, Crawford C M, Mitchell I M, et al. 2000. The degradation of fish-cage waste in sediments during fallowing. Aquaculture, 187(3-4): 351–366, doi: 10.1016/S0044-8486(00)00317-3
    Meyers P A. 1997. Organic geochemical proxies of paleoceanographic, paleolimnologic, and paleoclimatic processes. Organic Geochemistry, 27(5-6): 213–250, doi: 10.1016/S0146-6380(97)00049-1
    Miyake Y, Wada E. 1971. The isotope effect on the nitrogen in biological, oxidation-reduction-reactions. Records of Oceanographic Works in Japan, 11: 1–6
    Pan Zhe, Gao Qinfeng, Dong Shuanglin, et al. 2019. Effects of abalone (Haliotis discus hannai Ino) and kelp (Saccharina japonica) mariculture on sources, distribution, and preservation of sedimentary organic carbon in Ailian Bay, China: Identified by coupling stable isotopes (δ13C and δ15N) with C/N ratio analyses. Marine Pollution Bulletin, 141: 387–397, doi: 10.1016/j.marpolbul.2019.02.053
    Pancost R D, Boot C S. 2004. The palaeoclimatic utility of terrestrial biomarkers in marine sediments. Marine Chemistry, 92(1-4): 239–261, doi: 10.1016/j.marchem.2004.06.029
    Redfield A C, Ketchum B H, Rechards F A. 1963. The influence of organisms on the composition of seawater. In: Hill M N, ed. The Composition of Sea-Water, Comparative and Descriptive Oceanography. New York: Wiley Interscience, 26–27
    Ren Lihua, Zhang Jihong, Niu Yali, et al. 2015. Stable isotope evidence for the sediment impacts on biodeposits from long-line cultured Crassostrea gigas in Sungo Bay. Marine Sciences (in Chinese), 39(11): 79–85
    Rubio-Portillo E, Villamor A, Fernandez-Gonzalez V, et al. 2019. Exploring changes in bacterial communities to assess the influence of fish farming on marine sediments. Aquaculture, 506: 459–464, doi: 10.1016/j.aquaculture.2019.03.051
    Sarkar A, Chakraborty P, Nath B N. 2016. Distribution and nature of sedimentary organic matter in a tropical estuary: An indicator of human intervention on environment. Marine Pollution Bulletin, 102(1): 176–186, doi: 10.1016/j.marpolbul.2015.11.013
    Schubert C J, Calvert S E. 2001. Nitrogen and carbon isotopic composition of marine and terrestrial organic matter in Arctic Ocean sediments: implications for nutrient utilization and organic matter composition. Deep-Sea Research Part I: Oceanographic Research Papers, 48(3): 789–810, doi: 10.1016/S0967-0637(00)00069-8
    Srithongouthai S, Tada K. 2017. Impacts of organic waste from a yellowtail cage farm on surface sediment and bottom water in Shido Bay (the Seto Inland Sea, Japan). Aquaculture, 471: 140–145, doi: 10.1016/j.aquaculture.2017.01.021
    Talbot M R, Livingstone D A. 1989. Hydrogen index and carbon isotopes of lacustrine organic matter as lake level indicators. Palaeogeography, Palaeoclimatology, Palaeoecology, 70(1–3): 121–137
    Wada E, Hattori A. 1991. Nitrogen in the Sea: Forms, Abundances, and Rate Processes. Boca Raton, FL: CRC Press
    Wang Maolan, Lai Jianping, Hu Ketu, et al. 2014. Compositions and sources of stable organic carbon and nitrogen isotopes in surface sediments of Poyang Lake. China Environmental Science (in Chinese), 34(4): 1019–1025
    Wang Maolan, Zhao Liyue, Wan Yangjie, et al. 2022. Tracing the organic matter source of cage culture sediments based on stable carbon and nitrogen isotopes in Poyang Lake, China. Marine Pollution Bulletin, 182: 113943, doi: 10.1016/j.marpolbul.2022.113943
    Wu Xiaoyi, Yang Yufeng. 2011. Heavy metal (Pb, Co, Cd, Cr, Cu, Fe, Mn and Zn) concentrations in harvest-size white shrimp litopenaeus vannamei tissues from aquaculture and wild source. Journal of Food Composition and Analysis, 24(1): 62–65, doi: 10.1016/j.jfca.2010.03.030
    Yamada Y, Yokoyama H, Ishihi Y, et al. 2003. Historical feeding analysis in fish farming based on carbon and nitrogen stable isotope ratio in sediment. Fisheries Science, 69(1): 213–215, doi: 10.1046/j.1444-2906.2003.00609.x
    Ye Lixun, Ritz D A, Fenton G E, et al. 1991. Tracing the influence on sediments of organic waste from a salmonid farm using stable isotope analysis. Journal of Experimental Marine Biology & Ecology, 145(2): 161–174
    Yokoyama H, Abo K, Ishihi Y. 2006. Quantifying aquaculture-derived organic matter in the sediment in and around a coastal fish farm using stable carbon and nitrogen isotope ratios. Aquaculture, 254(1-4): 411–425, doi: 10.1016/j.aquaculture.2005.10.024
    Zhang Caixue, Lin Hongsheng, Sun Xingli. 2012. Source and distribution of organic matter in surface sediments of a typical bay in Guangdong Province, China. Journal of Tropical Oceanography (in Chinese), 31(6): 62–68
    Zhang Ling, Yin Kedong, Wang Lu, et al. 2009. The sources and accumulation rate of sedimentary organic matter in the Pearl River Estuary and adjacent coastal area, Southern China. Estuarine, Coastal and Shelf Science, 85(2): 190–196
    Zhou Fengxia, Gao Xuelu, Yuan Huamao, et al. 2018. The distribution and seasonal variations of sedimentary organic matter in the East China Sea shelf. Marine Pollution Bulletin, 129(1): 163–171, doi: 10.1016/j.marpolbul.2018.02.009
    Zhou Xin, Jin Guangzhe, Li Jiacheng, et al. 2021. Effects of typhoon Mujigae on the biogeochemistry and ecology of a Semi-Enclosed Bay in the Northern South China Sea. Journal of Geophysical Research: Biogeosciences, 126(7): e2020JG006031, doi: 10.1029/2020JG006031
    Zhuang Yanpei, Li Yangjie, Chen Ling, et al. 2023. Biogeochemical and physical controls on ammonium accumulation on the Chukchi shelf, western Arctic Ocean. Marine Environmental Research, 190: 106084, doi: 10.1016/j.marenvres.2023.106084
  • Relative Articles

  • Cited by

    Periodical cited type(4)

    1. Mengli Zhang, Chunling Zhang, Kefeng Mao, et al. Experimental construction of eddy real-time structure based on gradient-dependent OI in the Kuroshio-Oyashio confluence region. Progress in Oceanography, 2024, 224: 103262. doi:10.1016/j.pocean.2024.103262
    2. Zenghong Liu, Xiaogang Xing, Zhaohui Chen, et al. Twenty years of ocean observations with China Argo. Acta Oceanologica Sinica, 2023, 42(2): 1. doi:10.1007/s13131-022-2076-3
    3. Chunling Zhang, Danyang Wang, Zenghong Liu, et al. Global Gridded Argo Dataset Based on Gradient-Dependent Optimal Interpolation. Journal of Marine Science and Engineering, 2022, 10(5): 650. doi:10.3390/jmse10050650
    4. Huizan Wang, Yan Chen, Weimin Zhang. Some progress on ocean data assimilation in China: Introduction of the special section “Ocean Data Assimilation”. Acta Oceanologica Sinica, 2022, 41(2): 1. doi:10.1007/s13131-021-1982-0

    Other cited types(0)

  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(7)  / Tables(3)

    Article Metrics

    Article views (199) PDF downloads(9) Cited by(4)
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return