Red Sea MODIS Estimates of Chlorophyll a and Phytoplankton Biomass Risks to Saudi Arabian Coastal Desalination Plants

: Harmful algal blooms (HABs) and the high biomass associated with them have afﬂicted marine desalination plants along coastal regions around the world. Few studies of HABs have been conducted in the Red Sea, where desalination plants along the Saudi Arabian Red Sea coast provide drinking water for millions of people. This study was conducted along the Saudi Arabian Red Sea coast from 2014 to 2015 to assess the potential for using Moderate Resolution Imaging Spectroradiometer (MODIS) remote sensing of chlorophyll a (Chl a) or ﬂuorescence line height (FLH) to identify risks for biofouling at these desalination plants. Ship-based surveys of phytoplankton were conducted along the Saudi Arabian coastline offshore of desalination plants at Jeddah, Al Shoaibah and Al Qunfudhuh to assess the density of phytoplankton populations and identify any potential HAB species. Ship-based surveys showed low to moderate concentrations of phytoplankton, averaging from 1800–10,000 cells L − 1 at Jeddah, 2000–11,000 cells L − 1 at Al Shoaibah and 1000–20,500 cells L − 1 at Al Qunfudhuh. Sixteen different species of potentially toxigenic HABs were identiﬁed through these surveys. There was a good relationship between ship-based total phytoplankton counts and monthly averaged coastal MODIS Chl a (R 2 = 0.49, root mean square error (RMSE) = 0.27 mg m − 3 ) or FLH (R 2 = 0.47, RMSE = 0.04 mW m − 2 µ m − 1 sr − 1 ) values. Monthly average near shore Chl a concentrations obtained using MODIS satellite imagery were much higher in the Red Sea coastal areas at Al Qunfudhuh (maximum of about 1.3 mg m − 3 ) than at Jeddah or Al Shoaibah (maximum of about 0.4 and 0.5 mg m − 3 , respectively). Chlorophyll a concentrations were generally highest from the months of December to March, producing higher risks of biofouling desalination plants than in other months. Concentrations decreased signiﬁcantly, on average, from April to September. Long-term (2005–2016) monthly averaged MODIS Chl a values were used to delineate four statistically distinct zones of differing HAB biomass across the entire Red Sea. Sinusoidal functions representing monthly variability were ﬁt to satellite Chl a values in each zone (RMSE values from 0.691 to 0.07 mg m − 3 , from Zone 1 to 4). December to January mean values and annual amplitudes for Chl a in these four sinusoidal functions decreased from Zones 1–4. In general, the greatest risk of HABs to desalination occurs during winter months in Zone 1 (Southern Red Sea), while HAB risks to desalination plants in winter months are low to moderate in Zone 2 (South Central Red Sea), and negligible in Zones 3 (North Central) and 4 (Northern). Chattonella marina, Amphi-dinium sp., Alexandrium minitum, Dinophysis caudata, Dinophysis miles, sp., Peridinium sp., Goniodoma sp., Protoperidinium depressum, Protoperidinium divergens, Protoperidinium curtipes, Protoperidinium bipes, Protoperidinium ovum , and Protoperidinium ovatum


Introduction
Algal blooms frequently occur in coastal sea waters, examples include the Gulf of Mexico [1], south Florida coastlines [2], and Arabian Gulf [3]. These are primarily caused by nutrient enrichment, often associated with anthropogenic sources or natural upwelling. Algal blooms are also affected by temperature, wind, and salinity [4][5][6][7][8][9]. Sea color changes can be in response to chlorophyll a (Chl a) content associated with blooms. Harmful algal blooms (HABs) are characterized by high cell biomass that can cause biofouling of desalination plant filters that must be removed to protect the quality of drinking water. HABs are produced by several classes of microalgae, which are characterized by different colors ranging from blue-green to golden-brown to red. These classes may include diatoms, dinoflagellates, raphidophytes, cyanobacteria, prymnesiophytes, pelagophytes, and silicoflagellates [10].
Few studies have been conducted to measure the extent and severity of HABs in the Red Sea. An exception is the study of Genin et al. [11], which documented a serious incident involving algal and phytoplankton blooms in the Gulf of Aqaba as a result of the Mt. Pinatubo volcanic eruption. There is a serious need to better understand the potential frequency and severity of HABs in the Red Sea, especially along the eastern coast of the Red Sea bordering Saudi Arabia, which has three major desalination plants located at Jeddah, Al Shoaibah (100 km south of Jeddah), and Al Qunfudhuh (420 km south of Jeddah). These desalination plants provide drinking water for nearly 5 million residents. High biomass from coastal Red Sea HABs at these desalination plants could cause biofouling of the filters that protect drinking water supply [12]. Therefore, it is critically important to develop techniques to assess the frequency and severity of HABs in the eastern Red Sea along the Saudi Arabian coast.
Three contrasting methods have been developed to assess HABs in coastal regions of the world, including: (a) satellite optical remote sensing based on indirect measures of HAB biomass [2], (b) remote sensing of ecological characteristics such as sea surface temperature that affect HAB growth [13], or (c) association of HABs with other species, or with specific patterns in upwelling, wind, and ocean currents [14]. The focus of the present study is on optical remote sensing of Red Sea phytoplankton based on MODIS Aqua satellite imagery.
MODIS collects imagery at a revisit frequency of 16 days along the same ground track with a minimum pixel resolution of 250 m. However, the level three MODIS product for chlorophyll a is provided by NASA at a pixel size of 4.6 km. Since algal blooms vary significantly in density over distances on the order of 10 m, MODIS estimates of Chl a are coarser in scale than that measured with ship-based sampling of marine phytoplankton. Satellite imagery averages variability across a pixel and offers data that are more representative of larger spatial areas than ship-based sampling. This can be beneficial for coarser scale analyses, to better portray spatial and temporal evolution of HAB blooms and potential risks to desalination plants along the coast.
MODIS collects remote sensing bands in the ultraviolet (UV), purple, blue, green, and red wavelengths. These wavelengths have proven to be useful for detection of algal blooms in marine environments [2,10]. Two major optical remote sensing techniques for detecting algal blooms have emerged using satellite-based MODIS Aqua imagery. These include measurements of either: (a) Chl a, or (b) Fluorescence line height (FLH). Satellite-based Chl a algorithms have been developed for the detection of blooms of the genus Karenia, Gymnodinium catenatum, Lingulodinium polyedra, Chattonella sp., and noxious diatoms [15,16]. More generally, satellite-based Chl a estimates have been used to map algal blooms in the East China or Yellow Seas [15], the Gulf of Arabia [17], and the Red Sea [18]. FLH has proven successful for mapping harmful algal blooms associated with high ship-based Chl a concentrations and a red tide in previous studies from the southwestern Florida coast [2]. The feasibility of using MODIS FLH to identify HABs in the Red Sea has not previously been assessed [19]. The objective of this study is to evaluate the feasibility of using MODIS Chl a and FLH remote sensing products to assess phytoplankton biomass densities, and identify the locations and seasonal timing of the risks they pose to desalination plants along eastern coastal regions in the Red Sea.

Ship-Based Coastal Red Sea Surveys
Ship-based surveys were conducted offshore of Jeddah, Al Shoaibah and Al Qunfudhuh, Saudi Arabia (Figure 1 A YSI EXO2 sonde was used to measure sea surface temperature, Chl a, nitrate-N and turbidity. Results from the shipbased Chl a measurements were, however, unreliable due to a malfunction of the probe. Three replicate samples were collected from the sea surface during ship-based surveys at each of three distances (150, 750, and 1250 m) offshore of the three desalination plant locations at Jeddah, Al Shoaibah and Al Qunfudhuh ( Figure 1). Bathymetry of the Red Sea at these offshore distances ranged from 50 m to over 100 m depth. Results from these three offshore distances were averaged at each of the three locations for better comparison with 4.6 km scale MODIS Aqua satellite imagery pixels. A Niskin water sampler was used to collect three replicates of 1.7 L water samples at the surface from the three offshore distances at each coastal location [20]. Seven hundred mL was retained for chemical analysis of nitrogen (ammonium-, nitrite-and nitrate-N) in the lab. For determination of phytoplankton taxa [21] and quantitative determination of phytoplankton densities (counts), 1 L aliquots were preserved in amber-color bottles using a 4% iodine solution and stored in the dark at 4 • C pending analysis. The 1 L samples were concentrated via sedimentation in 1000 mL graduated cylinders according to the method detailed in Dolan and Marasse [22]. After settling for 4 to 8 days, the top 950 mL of the sample was slowly siphoned off with a small bore (0.5 cm diameter) tubing. The 50 mL remaining concentrate was used to count phytoplankton densities under an inverted microscope after settling for 24 h using Utermohl's method [23]. The entire surface of the settling chamber was examined at 400× magnification.

Chlorophyll a
The MODIS Aqua satellite platform operated by NASA [24] compiles and distributes monthly satellite composite images, which were acquired for the time period from Decem- Chlorophyll a data were acquired at a 4.6 km resolution as standard mapped images from the MODIS Aqua platform's level 3 dataset. Near-surface concentrations of Chl a in the MODIS aqua algorithm are calculated by NASA using an empirical relationship derived from ratios of blue-to-green band reflectances. MODIS radiance values were corrected by NASA using a calibration process described by [25] that improves accuracy along scan edges. Satellite data were typically not available on the exact dates of ship-based sample collection. Monthly composite satellite data acquired from NASA were compared with ship-based phytoplankton counts using linear regression. Although monthly composite satellite data do not correspond to an exact date for ship-based sample collections, they nevertheless represent spatial and temporal differences that occur from month to month or across coastal locations (Jeddah, Al Shoaibah and Al Qunfudhuh) with Saudi Arabian desalination plants.

Fluorescence Line Height
Fluorescence, the main deactivation pathway for photosystem II, depends on chlorophyll concentration as well as quantum yield and water clarity and has been used as a proxy for phytoplankton concentration [2,13]. MODIS fluorescence bands focus on the red part of the spectrum where chlorophyll fluorescence dominates the total signal, and therefore suffers little from other interferences (CDOM, shallow bottom, and the atmosphere). Monthly composite images for FLH were collected from NASA's MODIS Aqua satellite platform as 4.6 km standard mapped images [26] for the time period from December 2014 to June 2015. MODIS bands 13, 14, and 15 (centered at 667, 678, and 748 nm, respectively, with 10 nm bandwidth) were used by NASA to estimate the FLH [27]. A baseline is first formed between radiances for bands 13 and 15, and this baseline is then subtracted from band 14 radiance to obtain the FLH (FLH = Band 14 − (70/81) × Band 13 − (11/81) × Band 15). FLH satellite data from these NASA composites were compared with ship-based phytoplankton counts using linear regression.

Long-Term Red Sea Phytoplankton Biomass Mapping
Red Sea Phytoplankton Biomass Provinces have been delineated along lines of equal latitude by Raitsos et al. [18] based on Chl a concentrations measured with MODIS satellite imagery. This study refines these Provinces for the entire Red Sea based on long-term (Oct-May 2005-2016) spatial data trends in MODIS Chl a values. Data from June-September were not included, because phytoplankton densities are typically low during those months.
Our goal in this refinement was to better reflect spatial and temporal patterns in phytoplankton densities that are not aligned with the generalized latitudinal boundaries proposed by Raitsos et al. [18]. In particular, this refinement is intended to better represent spatial and temporal patterns in phytoplankton biomass near Saudi Arabian coastal areas where major desalination plants are located. Monthly mean MODIS values of Chl a from 2005-2016 were calculated, and refined curvilinear boundaries based on these data were hand-digitized to delineate the four Red Sea Phytoplankton Biomass Zones. Chlorophyll a data within these boundaries were fit with a sinusoidal function [28] Asin (ωt + ϕ) + m, where A is the fitted amplitude of Chl a (mg m −3 ), ω is frequency (2π/24), t is month (1 = January 12 = December), ϕ is phase shift (π) and m is the average fitted maximum value for Chl a (mg m −3 ). Accuracy of this model was determined by computing a root mean square error (RMSE). A one-way analysis of variance (ANOVA) test for mean differences was calculated on the data within each digitized boundary after checking the three assumptions for use of ANOVA.

Ship-Based Phytoplankton Densities
Ship-based surveys generally showed low to moderate concentrations of phytoplankton, ranging from 1800-10,000 cells L −1 at Jeddah, 2000-11,000 cells L −1 at Al Shoaibah and 1000-20,500 cells L −1 at Al Qunfudhuh. Over 100 phytoplankton species were identified, belonging to the Cyanobacteria, Dictyocophyceae, Dinophyceae and Bacillariophyceae groups.

Accuracy of MODIS Chl a and FLH Values for Estimating Total Phytoplankton
Ship-based surveys of phytoplankton in the Red Sea in this study were compared with MODIS Aqua satellite-based estimates of Chl a and FLH for the Jeddah, Al Shoaibah and Al Qunfudhuh sites (Figure 2). A good relationship was observed between ship-based total phytoplankton counts and monthly averaged MODIS Chl a or FLH estimates (Figure 2a  Monthly average nearshore FLH measurements obtained using MODIS satellite imagery in the Red Sea for this study were generally higher at Al Qunfudhuh than at Jeddah or Al Shoaibah (Figure 3, orange lines in right panels). At Al Qunfudhuh, the highest FLH values (0.17 mW m −2 µm −1 sr −1 ) occurred during February, and were more than double FLH values at Al Shoaibah and Jeddah during the same month. FLH values were generally negligible or zero in the months of April and May at all locations (Figure 3), except for Al Qunfudhuh in May. During these months densities of phytoplankton were always less than 3156 cells L −1 , and toxigenic HAB densities were always less than 756 cells L −1 .

Phytoplankton Biomass Zone Delineation
Phytoplankton biomass zone boundaries delineated by [18] were refined in the present study by incorporating information about spatial and temporal patterns in phytoplankton densities across the entire Red Sea based on long-term (Oct-May 2005-2016) trends in MODIS Chl a values (Figure 4). Temporal patterns in Chl a ( Figure 4) were accurately represented with sinusoidal functions for each Zone (RMSE ranging from 0.691 to 0.07 mg m −3 , in decreasing order from Zone 1 to 4) that peaked in the months of December to January and were lowest in the months of June to July. December to January mean values and annual amplitudes for Chl a in these four sinusoidal functions decreased from Zone 1 to  The spatial representation of Phytoplankton Biomass Zones ( Figure 5) indicates that Zones 1, 3 and 4 cover a majority of the Red Sea, while Zone 2 is essentially a transitional area between Zones 1 and 3. In contrast with a previous study on phytoplankton zone delineation [18], the refined zones in the present study show a marked northern extension of zone boundaries along coastal regions along both the eastern and western coastlines of the Red Sea. This northward extension has important implications for protection of coastal desalination plants, as these coastal areas have higher phytoplankton biomass than would be inferred using previously defined zones [18] that are delineated using horizontal latitude-based boundaries. MODIS satellite imagery in the Red Sea (exemplified by trends during 2015 as shown in Figure 6) shows a pattern for high HAB biomass that develops offshore of Yemen in Phytoplankton Zone 1 near the Gulf of Aden during October. This high phytoplankton biomass subsequently migrates in the northeasterly direction towards Al Qunfudhuh during November, December, and January before decreasing in March. Based on data collected during our ship excursions near the western Saudi Arabian coastline (Section 3.1), nearly half the phytoplankton biomass in Zone 1 during the months of January, February and March may be comprised of toxigenic species. Phytoplankton Zone 2, which includes much of the Red Sea shoreline, reaches a peak HAB biomass in December and January. Phytoplankton Zone 3 reaches a maximum in January (Figure 4), while Zone 4 peaks in February (Figure 4). However, HAB biomass is minimal in Zones 3 and 4.

Discussion
In comparison with results from this study, ship-based surveys conducted along the southwestern Florida coastline during an HAB event [2] showed moderate (10 4 to 10 5 cells L −1 ) to high (>10 5 cells L −1 ) concentrations of Karenia brevis. The highest concentrations of phytoplankton in the present Red Sea study were in the low part of the high range in the Florida study [2]. During a severe 2008 HAB event in the Arabian Gulf [3,[30][31][32], very high (>10 7 cells L −1 ) concentrations of Cochlodinium polykrikoides caused the closure of five desalination plants and resulted in severe fish mortality, similar to effects of severe HAB events in other parts of the world [33][34][35][36][37][38][39][40][41]. Concentrations of HABs in the 2008 Arabian Gulf event [3] were much higher than any phytoplankton concentrations measured in the Red Sea during the present study. Thus, the risk from HABs to coastal desalination plants is smaller along the Red Sea Saudi Arabian coastline than along the Arabian Gulf coastline.
Based on results presented in Figure 2 above, MODIS satellite Chl a values are a good predictor of total phytoplankton in the eastern coastal region of the Red Sea. Accuracy of the Chl a model presented in Figure 2a (R 2 = 0.49, RMSE = 0.27 mg m −3 ) was much better than the Chl a model (R 2 = 0.27) developed for coastal Florida in a 30 km study area during December of 2004 [2]. Accuracy of the FLH model presented in Figure 2b (R 2 = 0.47, RMSE = 0.04 mW m −2 µm −1 sr −1 ) was somewhat weaker than the FLH model (R 2 = 0.58) developed in the Arabian Gulf and Gulf of Oman from 2006 to 2011 [19]. These comparisons provide confidence that the Chl a and FLH models in the present study are an accurate measure of phytoplankton abundance in the Red Sea. However, monthly composite images from MODIS were not always temporally well aligned with ship-based measurements of phytoplankton (e.g., during May and June of 2015, when the ship excursions were a week later or a week earlier, respectively, than the MODIS imagery used to construct monthly composites). This mismatch in temporal alignment could have contributed to a slight weakening of the Chl a and FLH models developed in this study. During severe algal blooms observed in the Arabian Gulf starting in August 2008 [30][31][32], Chl a concentrations often exceeded 2 mg m −3 , much higher than the values observed in the present study. During the severe algal blooms that were observed in a previous study located along the southwestern Florida coastal region starting in October 2004 [2], FLH values often exceeded 0.2 mW m −2 µm −1 sr −1 . FLH values at Jeddah or Al Shoaibah were generally lower than these levels, and rarely exceeded 0.1 mW m −2 µm −1 sr −1 .
In cases where a single toxigenic HAB species was responsible for red tides (e.g., in Florida [2] and Arabian Gulf [19]), satellite estimates of FLH were superior to Chl a at detecting toxigenic HAB species. The relationship between FLH and total phytoplankton (or indirectly the relatively constant proportion of toxigenic HABs) was very good in our study, but not as strong as the relationship between Chl a and total phytoplankton (Figure 2a vs. Figure 2b). However, there were more than 15 toxigenic HAB species in our study, and this may be why FLH did not perform as well as in the Florida or Arabian Gulf studies at identifying toxigenic HAB species.
Ship-based surveys of phytoplankton in the offshore areas of Jeddah, Al Shoaibah and Al Qunfudhuh were limited in spatial extent and were time consuming to conduct. Satellite remote sensing, in contrast, covers the entire Red Sea at a pixel resolution of 4.6 km on a side (area of 21.16 km 2 ), and images are collected at regular time intervals using a consistent methodology. This suggests the possibility of monitoring phytoplankton biomass in the Red Sea using satellite remote sensing, in order to have early warning of phytoplankton densities that could pose risks to Saudi Arabian coastal desalination plants.
High levels of turbidity can interfere with MODIS Chl a estimates [42]. However, measured turbidities in the coastal Red Sea areas from the present study were relatively low. This observation is consistent with previous findings showing that the entire Red Sea basin (excepting the Farasan Islands in the Southern Red Sea) is classified as optically clear Case-I waters [43]. There are no major rivers discharging sediment or colored dissolved organic matter into the eastern side of the Red Sea. Average annual ship-based turbidity levels in the present study at Jeddah and Al Shoaibah were 0.64 and 0.61 Formazin Nephelometric Units (FNU), respectively, with values that increased to only 1.0 FNU during the months of January through March. These highest FNU values are at the low end of FNU values measured in the Gulf of Finland during August of 1997 [44], which ranged from 1 to 8 FNU. FNU values at Jeddah and Al Shoaibah are low, implying that turbidity is unlikely to interfere with MODIS Chl a estimates in most months. Average annual ship-based turbidity at Al Qunfudhuh was 1.1 FNU, and highest turbidity readings (>2 FNU) at this location occurred during April to August, when ship-based total phytoplankton densities were generally lowest. Turbidity at Al Qunfudhuh is likely to interfere only slightly with MODIS estimates of Chl a during the months from April to August, when lowest total phytoplankton densities were observed. In future studies, it is recommended to use a MODIS two-band (red and near infrared) model [45] for Chl a estimates in the Southern Red Sea, to further reduce concerns over potentially small interferences from turbidity [42], particularly at Al Qunfudhuh.
The Red Sea is recognized as an oligotrophic marine water body [18], with nutrient limitations from nitrate-N becoming less severe from the northern regions to the southern regions. Seasonal trends in total phytoplankton counts from this study are consistent with strong winds from the southeast in October, bringing an influx of cold, low salinity, nutrient rich water to the southern Red Sea from the Gulf of Aden [18,46] In the present study, ship-based nitrate-N concentrations at Al Qunfudhuh increased from 0.15 mg L −1 in August to 1.3 mg L −1 in November [21]. This influx of nutrients initiates a marked increase in phytoplankton populations from summer to late fall in the portion of the southern Red Sea that includes Al Qunfudhuh. Seasonal trends in phytoplankton offshore of Jeddah and Al Shoaibah generally occur for the same reasons as trends at Al Qunfudhuh, but the influx of south easterly winds and cold, nutrient rich water from the Gulf of Aden are greatly diminished at Jeddah and Al Shoaibah. This trend, combined with lower sea surface temperatures [46] at Jeddah and Al Shoaibah than Al Qunfudhuh (ship-based SST values were about 2 • C smaller at Jeddah and Al Shoaibah than at Al Qunfudhuh (29.7 • C) during November), causes productivity at Jeddah and Al Shoaibah to be smaller than at Al Qunfudhuh. Wind reverses direction in summer, causing nutrient deficiencies that limit phytoplankton growth in the southern Red Sea. SST increases on average by 6 • C from winter to summer in any given portion of the Red Sea [18,46], but more importantly, there is a marked increase in SST from southeast to northwest during any given season, driving increased phytoplankton biomass from southeast to northwest.
The lowest phytoplankton biomass at Al Qunfudhuh (as well as Jeddah and Al Shoaibah) typically occurs during the late summer months. An exception to this pattern of phytoplankton biomass in relation to seasonal patterns in wind and SST, were the uncharacteristically high phytoplankton counts during June. This pattern is explained by Raitsos et al. [18] as being caused by anti-cyclonic eddies that temporarily bring colder, nutrient rich waters [46] or nutrient-rich wind-blown dust to the Al Shoaibah and Al Qunfudhuh regions. Remote sensing offers the potential for early detection of these upwelling anomalies in order to better protect Saudi Arabian coastal desalination plants.
Sea water desalination provides the ultimate source for fresh water in Arabian region [47], which has domestic, industrial and agricultural uses. However, HABs impose great challenges for production of desalinized water [47][48][49]. Monitoring and early warning of HAB events is key to managing desalination processes and water production [50]. Remote sensing analysis has been extensively explored to establish reliable early warning tools to predict HAB events in this study and other parts of the world [30][31][32][51][52][53]. In future, it would be useful to combine remote sensing with more general models for risk of HAB biomass to coastal desalination plants in the Red Sea that simultaneously account for spatial and temporal patterns in satellite-based SST and Chl a concentrations, along with measurements of wind speed and direction.

Conclusions
Satellite-based phytoplankton risk zones for the onset of HAB biomass risks in marine coastal regions are greatly needed to enhance the protection of desalination plants, as well as for the protection of marine habitats, fisheries, and recreation areas. HAB presence and severity has been rarely studied in the Red Sea, and has never been studied using satellite based FLH products. MODIS satellite imagery for Chl a values were used to assess HAB risks and delineate four Phytoplankton Biomass Zones with distinctly different spatial and temporal patterns of Chl a values along Saudi Arabian coastal regions in the Red Sea for the first time. These coastal regions include Jeddah, Al Shoaibah and Al Qunfudhuh, which each rely heavily on drinking water supplies produced from desalination plants. HAB risks were significant offshore of Al Qunfudhuh and in the Red Sea south of Al Qunfudhuh (Phytoplankton Biomass Zone 1) during fall and winter months. Therefore, this region exhibits serious HAB risks in winter months that could produce biofouling threats to nearby coastal desalination plants. In contrast, Phytoplankton Biomass Zone 2, which is near the populated coastal Saudi Arabian areas north of Al Qunfudhuh (e.g., Jeddah and Al Shoaibah), showed low to moderate levels of HAB biomass during fall and winter months. Remote sensing of Chl a values shows promise for use in estimating the timing and magnitude serious HAB onset in the Red Sea, just as it has shown promise for the same purpose in the Arabian Gulf, southwestern coastal Florida region, and East China and Yellow Seas.

Data Availability Statement:
The data presented in this study are available on request from the corresponding authors. The data are not publicly available due to privacy concerns.