Remotesensing 10 01604
Remotesensing 10 01604
Article
Mapping Fishing Activities and Suitable Fishing
Grounds Using Nighttime Satellite Images
and Maximum Entropy Modelling
Rollan C. Geronimo 1,2, * , Erik C. Franklin 1,3 , Russell E. Brainard 2 , Christopher D. Elvidge 4 ,
Mudjekeewis D. Santos 5 , Roberto Venegas 2 and Camilo Mora 1
1 Department of Geography, University of Hawai’i at Mānoa, Honolulu, HI 96822, USA;
[email protected] (C.M.); [email protected] (E.C.F.)
2 Ecosystem Sciences Division, Pacific Islands Fisheries Science Center, National Marine Fisheries Service,
U.S. National Oceanic and Atmospheric Administration, Honolulu, HI 96818, USA;
[email protected] (R.E.B.); [email protected] (R.V.)
3 Hawaii Institute of Marine Biology, School of Ocean and Earth Science and Technology,
University of Hawai’i at Mānoa, Kaneohe, HI 96744, USA; [email protected]
4 Earth Observation Group, National Oceanic and Atmospheric Administration’s National Centers
for Environmental Information, Boulder, CO 80305, USA; [email protected]
5 National Fisheries Research and Development Institute, Quezon City 1103, Metro Manila, Philippines;
[email protected]
* Correspondence: [email protected]; Tel.: +1-808-956-8465
Received: 13 August 2018; Accepted: 28 September 2018; Published: 9 October 2018
Abstract: Fisheries surveys over broad spatial areas are crucial in defining and delineating appropriate
fisheries management areas. Yet accurate mapping and tracking of fishing activities remain largely
restricted to developed countries with sufficient resources to use automated identification systems
and vessel monitoring systems. For many countries, the spatial extent and boundaries of fishing
grounds are not completely known. We used satellite images at night to detect fishing grounds in
the Philippines for fishing gears that use powerful lights to attract coastal pelagic fishes. We used
nightly boat detection data, extracted by U.S. NOAA from the Visible Infrared Imaging Radiometer
Suite (VIIRS), for the Philippines from 2012 to 2016, covering 1713 nights, to examine spatio-temporal
patterns of fishing activities in the country. Using density-based clustering, we identified 134 core
fishing areas (CFAs) ranging in size from 6 to 23,215 km2 within the Philippines’ contiguous maritime
zone. The CFAs had different seasonal patterns and range of intensities in total light output, possibly
reflecting differences in multi-gear and multi-species signatures of fishing activities in each fishing
ground. Using maximum entropy modeling, we identified bathymetry and chlorophyll as the main
environmental predictors of spatial occurrence of these CFAs when analyzed together, highlighting
the multi-gear nature of the CFAs. Applications of the model to specific CFAs identified different
environmental drivers of fishing distribution, coinciding with known oceanographic associations for
a CFA’s dominant target species. This case study highlights nighttime satellite images as a useful
source of spatial fishing effort information for fisheries, especially in Southeast Asia.
1. Introduction
Monitoring and mapping of fishing activities are critical components of planning and management
for marine fisheries [1]. Fishing location data have been used to identify and delineate fishing
grounds [2], improve stock assessments [1,3], estimate fishing effort [4], and evaluate the impact
Table 1. Major fishing gears utilizing lights to attract fishes in the Philippines [37]
Figure 1. Flow of application of VIIRS boat detection data to map out potential fishing grounds for
small pelagic fishes in the Philippines. VBD: VIIRS boat detection; CFA: core fishing area; HDBSCAN:
hierarchical density-based spatial clustering of applications with noise.
clustering of applications with noise (HDBSCAN) [43] using different parameter combinations of
‘minimum cluster size’ (i.e., 3 to 30) and ‘minimum samples’ (i.e., 3 to 35), resulting in 1200 different
clustering combinations (Supplementary Figure S1). We selected the clustering parameter combination
that yielded the least number of points classified as ‘noise’ while maintaining reasonable clustering of
known small pelagic fish fishing areas (e.g., Northeastern Palawan round scad fishing grounds; Tayabas
Bay fisheries; and Basilan strait sardine fisheries). Concave hull polygons were estimated for each
cluster of points using a Delaunay triangulation algorithm to visualize the clusters. Finally, using high
resolution images from Google Earth and shipping traffic data from https://www.marinetraffic.com/
as reference (Supplementary Figure S2), we developed a raster mask to remove polygons in areas with
possible high occurrence of non-fishing related shipping activities (e.g., around busy harbors, shipping
lanes, and areas with mining activities near the coast) since VBD points in these high-traffic areas
could include a large number of non-fishing vessels that turn on additional deck and navigational
lights. This mask was also used for the MaxEnt modeling of suitable fishing grounds based on annual
VBD occurrence. Annual raster overlays were done using R (v. 3.5.1, R Development Core Team) [44],
while hierarchical density-based spatial clustering of applications with noise [45] and concave hull
estimations were ran using Python 2.7 (www.python.org).
2.3. Clustering CFAs Based on Monthly Patterns and Radiance Values of VBDs
Each CFA can be characterized by the composition of fishing gears operating in the area and the
seasonal patterns in fishing. Multiple light-assisted fishing gears (Table 1) most likely operate in each
fishing ground. Radiance values can be used as a relative measure of light source intensity with higher
radiance corresponding to brighter lights [46]. A narrow radiance range could imply dominance of
a single type of fishing gear while a wide range could mean that multiple gears are operating in the
area. Seasonal patterns in fishing activities, on the other hand, could be defined by the mean nightly
VBD by month, averaged across years within each CFA.
For radiance values, we applied histogram-valued data analysis based on the L2 Wasserstein
metric between distributions on log-transformed radiance per CFA [47]. We limited this analysis to VBD
points within CFAs occurring during nights with less than 50% moon illumination and less than 50%
cloud cover since both factors are known to affect recorded radiance beyond this percentage [21,46]
(Supplementary Figure S3). Clustering was performed using the R Package “HistDAWass” [48].
For seasonal patterns in fishing activities within each CFA, we used the mean nightly VBD by month,
averaged across years, and standardized the values between 0 and 1, corresponding to the minimum
and maximum mean nightly VBD count by month. CFAs were grouped based on hierarchical clustering
using Ward’s linkage method [49].
fishes are also caught using other gears that do not employ powerful lights, we conservatively treated
the VBD locations as presence-only data and randomly derived pseudo-absences or background
locations outside the CFAs. MaxEnt is one of the few models available that can handle presence-only
data and has been used extensively for species distribution modeling [53,54]. It estimates the
relative suitability of one place versus another by comparing the conditional density of covariates at
presence sites with the unconditional density of covariates across the study area [55]. The probability
of occurrences are returned which is relative to the ‘typical’ conditions generated from available
presence points.
We applied MaxEnt at an aggregated 4 km resolution to match the resolution of remotely-sensed
environmental datasets used as predictors (Table 2). The model domain is the Philippines’ 24 nautical
mile contiguous zone bounded from 4◦ to 22◦ North and from 116◦ to 128◦ East. ‘Presence’ locations
were defined as the centroid coordinates of 4-km grid cells that had VBD data and were inside any
of the CFAs for the periods 2013 to 2016. VBD data for 2012 were incomplete (i.e., starting only in
April 2012) so they were not used in the MaxEnt models. Background points (~10,000) were randomly
selected within the model domain. MaxEnt’s regularization parameter, controlling for over-fitting,
and the selection of feature classes (e.g., linear, quadratic, product, threshold, or hinge) [55] were
tuned for each model using the ENMevaluate function in the ENMeval package in R [56]. For each
environmental predictor, the ‘auto feature’ selection method was used to allow MaxEnt to determine
the appropriate feature class to use for each environmental predictor in each model with ‘maximum
iterations’ set at 5000.
We collated environmental covariates commonly associated with the distribution of small pelagic
fishes [57–60] (Table 2; Supplementary Figure S4). Gridded bathymetry was downloaded from the
General Bathymetric Chart of the Oceans (GEBCO) website. Monthly netCDF files for the other
environmental predictors were downloaded from the NOAA Coast Watch’s Environmental Research
Division’s Data Access Program archives for the years 2013 to 2016. Climatological and annual
means were calculated for each predictor across the four years. All raster layers were projected to
a Behrmann’s equal area projection at 4km resolution using bilinear interpolation under the ‘raster’
package in R [61].
Table 2. Data used for the MaxEnt modeling of core fishing areas
Original
Predictor Units Source
Resolution
NOAA Earth Observation Group
VIIRS boat detection 742m –
(https://ngdc.noaa.gov/eog/viirs/download_boat.html)
GEBCO (https://www.gebco.net/data_and_products/
Bathymetry 30 arc-sec meters
gridded_bathymetry_data/gebco_30_second_grid/)
Sea surface temperature 4 km ◦C MODIS-AQUA from NOAA CoastWatch
Surface chlorophyll a 4 km mg/cm3 (https://coastwatch.pfeg.noaa.gov/erddap/index.html)
1/12◦ (9
Sea surface salinity psu
km) HYCOM + NCODA Global 1/12◦ Analysis
Mixed layer thickness (at 1/12◦ (9 (https://coastwatch.pfeg.noaa.gov/erddap/index.html)
meters
density change of 0.03 kg/m3 ) km)
Sea surface height 1/12◦ (9km) meters
WaveWatch III Global Wave Model
Significant wave height 0.5◦ meters (https://coastwatch.pfeg.noaa.gov/erddap/griddap/
NWW3_Global_Best.html)
The exact location of fishing activities, even within the same CFA, often changes across years.
To account for annual variability of dynamic predictors on annual aggregated locations of VBD inside
CFAs, we developed annual MaxEnt models from 2013 to 2016, using respective annual VBD data
and annual averages for each predictor. Since the CFAs most likely differ in terms of the fishing fleet
composition and targeted fish species, we selected two CFAs for which the dominant target species
group is known and where most fishing gears that use lights are almost exclusively targeting these
species. These CFAs are (a) the round scad fishery in Northeastern Palawan island (CFA # 106 in Figure
3) and (b) the sardine fisheries in the Sulu Archipelago (CFA # 70).
To identify other areas within the Philippines’ 24 nautical mile contiguous zone with
environmental characteristics similar to the delineated CFAs, we generated predicted CFA suitability
maps for all 10 MaxEnt models (Table 3). Mean predicted suitability values per model was generated
from suitability maps produced in each of the cross-validation runs (n = 10 per model). The resulting
suitability distributions were compared with the location of existing CFAs.
each individual variable and measuring changes in training or testing data gain. Response curves
were generated and examined to derive the preferred range of environmental variables of potential
fishing areas.
3. Results
VIIRS boat detection (VBD) locations within the Philippine’s 24-nautical mile contiguous zone
are highly clustered with specific areas being frequently visited (Figure 2a). Spatial patterns are also
apparent in terms of mean radiance values with most VBD locations having radiance values greater
than 2.0 nanoW/cm2 /sr (Figure 2b). VBD with less than 2.0 radiance values are clearly clustered
together in some areas. Across years, the spatial patterns of VBD aggregations within the Philippines
has not change substantially from 2012 to 2016.
Figure 2. VIIRS boat detection (VBD) data from April 1, 2012 to December 31, 2016 aggregated on
a 700 × 700m resolution grid showing (a) number of VBD detection per pixel over 1,713 nights and
(b) mean radiance value. The 24 nautical mile contiguous zone shapefile for the Philippines came from
Flanders Marine Institute’s Marine Regions database [40].
The 134 CFAs ranged in size from 6 km2 to 23,215 km2 with a mean size of 700 km2 . The maximum
number of VBD for a given night was detected in March 3, 2014 with 1,124 lit-up pixels within the entire
Philippines Contiguous Zone, 87% or 974 of these pixels were within CFAs. Per CFA, the maximum
number of lights detected for a single night ranged from 2 to 227.
Figure 3. Core fishing areas (n = 134) identified from HDBSCAN and after removing clusters in areas
that overlap with significant vessel traffic and close to major ports and piers, including mining jetties.
Details for each CFAs and corresponding names can be found in Supplementary Table S1.
Seasonal patterns in fishing activity data varied per CFA (Figure 4 and Supplementary Figure S5-a).
Cluster analysis applied to averaged monthly VBDs showed two different groups based on the months
when peak fishing activity occur, generally corresponding with the monsoon months: (a) northeast monsoon
(November to April) and (b) southwest monsoon (May to October).
Remote Sens. 2018, 10, 1604 10 of 23
Figure 4. Hierarchical clustering of core fishing areas based on scaled mean nightly VBD counts per
month using Ward’s linkage method showing two groups based on months of peak fishing activity:
(a) southwest monsoon (May to October) and (b) northeast monsoon (November to April).
Mean radiance values per CFA ranged from 1 to 166 nanoWatts/cm2 /sr (Figure 5; Supplementary
Table S1). Hierarchical clustering on radiance distribution values per CFA identified two distinct
groups (h cutoff = 5): (a) low radiance values (group mean = 3.4 nanoWatts/cm2 /sr) and (b) high
radiance values (group mean = 61.2 nanoWatts/cm2 /sr; Figure 5 and Supplementary Figure S5-b).
Remote Sens. 2018, 10, 1604 11 of 23
Figure 5. Hierarchical clustering of log-10 histograms of radiance values by CFA for nights with less
than 50% moon illumination and cloud cover. Colors represent probability density of radiance bin per
CFA. High-level groupings (around blue box) are based on spread of values with group (a) indicating
high kurtosis, narrow range, and low radiance values while group (b) are high radiance values with
large spread and even some bimodality.
Remote Sens. 2018, 10, 1604 12 of 23
Figure 6. Density distribution of environmental predictors inside CFAs (red) versus outside (blue)
within the Philippines’ 24 nautical mile contiguous zone. Bathymetry and chlorophyll a values log-10
transformed for clarity. Non-transformed values were used in the MaxEnt models.
All MaxEnt models performed better than random (i.e., AUC > 0.5 and TSS > 0.4; Table 4). Among
all the models, the one without bathymetry had the lowest AUC and TSS scores. The full climatology,
monsoonal, and annual models all had closely similar accuracy scores. The CFA-specific models
performed much better than the other models, with high AUC and TSS scores (greater than 0.9 for
both), highlighting the unique environmental attributes of these CFAs.
Remote Sens. 2018, 10, 1604 13 of 23
Table 4. Model tuned parameters and evaluation scores for the different MaxEnt models
Regularization Feature
Model AUC TSS
Parameter Classes *
Full climatology 25 LQHP 0.88 0.59
Full–no bathymetry 1 H 0.84 0.52
Northeast monsoon 2 LQHP 0.87 0.58
Southwest monsoon 17 LQH 0.87 0.58
2013 13 LQHP 0.88 0.61
2014 1 LQHP 0.89 0.62
2015 1 LQHP 0.88 0.60
2016 3 LQHP 0.87 0.58
Northeast Palawan 23 LQHP 0.99 0.90
Sulu 18 LQHP 0.98 0.94
* L: linear; Q: quadratic; H: hinge; P: product.
Figure 7. The relative contribution of different environmental predictors in determining the overall fit of the
MaxEnt models to the VBD presence data used (normalized to percentages). Values averaged across two
replicates of five-fold cross-validation runs of each MaxEnt model. Error bars are ±1 standard deviation.
Chl-a: chlorophyll a; MLD: mixed layer thickness (by density); SSH: sea surface height; SSS: sea surface
salinity; SST: sea surface temperature; SWH: significant wave height.
Remote Sens. 2018, 10, 1604 14 of 23
Based on the response curves from the MaxEnt full climatology model, the logistic output or
probability of occurrence of a site becoming a core fishing area increased with depth from the shallowest
depths up to around 100 m then decreased rapidly to depths of 500 m and gradually declined further
with depths deeper than 4000 m. The probability also peaked with chlorophyll a values between 0.3 and
2.0 mg/cm3 . Frequently visited fishing sites also occurred more in areas with sea surface temperatures
between 28.6 ◦ C and 30.6 ◦ C. The effect of sea surface salinity, sea surface height, and mixed layer
thickness in discriminating CFA versus non-CFA areas was smaller compared to the first three variables
mentioned as indicated by the almost flat response across their entire range of values.
The predicted suitability of marine areas for light-assisted fishing based on the climatology and
annual MaxEnt models are shown in Figure 8. Predicted suitable areas with environmental conditions
typical of the CFAs were found in shallow, nearshore areas around existing CFAs. Predicted suitable
areas did not vary greatly across monsoons and years. For the CFA-specific MaxEnt models, we found
highly constrained areas of high suitability. However, the predicted areas corresponded with other
CFAs that were not included in the model (e.g., CFA # 6, 18, 80-94, and 102 for the Northeast Palawan
round scad CFA and CFA # 39, 53, 59, 67, 69, 71, and 72 for the Sulu sardine CFA).
Figure 8. Predicted suitability maps for core fishing areas for the selected MaxEnt models. Color scale
represents degree of suitability relative to ‘typical’ conditions found inside CFAs with 1 being highly
suitable areas and 0 are unsuitable.
Remote Sens. 2018, 10, 1604 15 of 23
4. Discussion
The VIIRS nighttime imaging capabilities significantly enhanced those of the DMSP satellites,
particularly in terms of radiometric and spatial resolution, allowing for the nightly extraction of
locations and radiance information of maritime lights [21,32]. We applied this freely accessible and
novel dataset to identify, delineate, and characterize core fishing areas (CFAs) and determine the
suitability of other parts of Philippine coastal waters to light-assisted fisheries and their target species.
Many of the large CFAs identified are known primary fishing grounds for various coastal pelagic
fishes. The Sulu CFA (#70) in Southern Mindanao is dominated by sardine purse seine fishing. It is one of
the largest sources of the country’s sardine (Sardinella lemuru) supply and has been subjected to seasonal
fishery closures between November and March since 2011 due to concerns about overfishing and stock
declines [69]. The Northeast Palawan CFA (#106) and those west of Palawan (e.g., CFAs # 102 to 117) are
known round scad (Decapterus macrosoma and Decapterus russelli) fishing areas, also subjected to annual
seasonal closure since 2015 [70]. Visayan Sea 1 (#37) and 2 (#40) CFAs are also part of a small pelagic fish
seasonal closure established back in the 1970s primarily for sardines (Sardinella gibbosa).
A similar study was done globally, identifying clusters or ‘agglomeration areas’ of marine lights
from cloud-free monthly composites of the VIIRS day/night band images [1]. The authors limited
the presentation of their results to only the top 100 agglomeration areas by size. Within this list, they
identified six (6) agglomeration areas for the Philippines. All these six agglomeration areas were also
identified as CFAs in our study. Twenty-one of the 134 CFAs are larger than 1000 km2 . Unfortunately,
information on the Philippine agglomeration areas were not available in Zhao et al. [71] since they only
presented statistics for the top 50 of their top 100 agglomeration areas. The difference in the number of
agglomeration areas and CFAs identified for the Philippines could be due to Zhao et al.’s [71] use of
only cloud-free monthly composites and a fixed distance threshold (i.e., 10 km) to group pixels.
Bathymetry, chlorophyll a, and sea surface temperature have all been cited as predictors of small
pelagic fish spatial distribution in previous studies. Catches and fishing grounds for the Indian
mackerel (Rastrelliger kanagurta) in Indonesia correlated significantly with chlorophyll a and sea surface
temperature [72]. Bottom depth and sea surface temperature explained variability in European sardine
(Sardinella pilchardus) presence in the Mediterranean Sea while bottom depth and surface chlorophyll
concentration predicted European anchovy (Engraulis encrasicolus) distribution [73–75]. These three
variables, along with salinity, are also used in mapping global distributions of fish species through
Aquamaps (www.aquamaps.org) [76,77]. These variables influence suitable habitats for small pelagic
fishes through zooplankton production and availability [78].
Bathymetry consistently ranked highest in discriminating between CFA and non-CFA sites,
accounting for the most variability not captured in the other variables. Chlorophyll accounted for
most of the variability in the MaxEnt model that excluded bathymetry indicating some correlation
between chlorophyll a and bathymetry (Pearson’s correlation r = 0.25; p < 0.01). Indeed, chlorophyll
a tended to decline with increasing depth. However, zooming in on some of the CFAs shows how
VBD data track, with surprising precision, bathymetry, and depth contours using a 1-km resolution
bathymetry data (Figure 9). Removing bathymetry as a predictor resulted in a lower accuracy model
(AUC from 0.88 to 0.84 and TSS from 0.59 to 0.52; Table 4) and broader predicted areas of suitability
(Figure 8). Aggregating data at a 4-km resolution removes these fine-scale spatial patterns in nighttime
lights and fishing sites. Finer resolution environmental data could explain more of the variability in
the distribution of nighttime lights in CFAs. Unfortunately, current available monthly remote sensing
data for chlorophyll and sea surface temperatures for this region are limited to a 4-km resolution
(i.e., MODIS images) while finer oceanographic models are not available for the entire study area.
Remote Sens. 2018, 10, 1604 16 of 23
Figure 9. CFAs overlay on a bathymetry map of the Philippines showing most of the CFAs are located
in less than 200 meter depths. Inset maps show how VBD points align with depth contours for
(a) Northeast Palawan and (b) Visayan Sea.
The dominance of bathymetry in distinguishing CFA from non-CFA sites could be attributed to
(1) the depth preference of target species and (2) the multi-species nature of fisheries within CFAs.
A search in Fishbase (www.fishbase.org) on depth information for the dominant small pelagic fishes
caught in Philippines waters shows that commonly targeted small pelagic fishes by the Philippine
fishing fleet are often found in depths of less than 200 m [63]. In addition, CFAs most likely have
a large mix of different fishing gears that target various species of small pelagic fishes. The presence of
lights in deeper areas, although not dense enough or frequently visited to classify as a core fishing
area, indicate that fishing vessels using lights are able to fish further offshore and not constrained
by vessel capability. Since different species or species group tend to respond differently to water
column characteristics (e.g., SST and chlorophyll i concentrations), the common denominator in their
environmental preference becomes bottom depth, a static variable which then gets captured as the
main predictor of CFA versus non-CFA sites. If the CFAs could be differentiated according to their
Remote Sens. 2018, 10, 1604 17 of 23
dominant target species, the effect of dynamic predictors in differentiating CFAs from non-CFAs for
that particular species could be bigger than bathymetry.
One of the main uses of the identified core fishing areas and habitat suitability map in this study
is to delineate fisheries management areas as part of implementing ecosystem approaches to fisheries
management (EAFM) [8,9]. While we have provided here a way to identify frequently visited fishing
grounds and predict habitat suitability of non-fishing grounds based on nighttime lights activity and
environmental attributes, delineation of meaningful fisheries management unit as part EAFM should
be done in consultation with relevant stakeholders [2]. Some of the CFAs could be grouped into larger
fisheries management areas based on similarities in their target species and fishing gears used. In the
absence of actual fishing data available for each CFA, the seasonality of fishing activities and the
distribution of VBD radiance values could be used as a proxy for these characteristics. For example,
the northern Palawan CFAs (e.g., #102 to 117) can be grouped into one FMA since they share the
same seasonality in fishing activities (during northeast monsoon; Figure 4) and radiance signature
(Figure 5; see also Supplementary Figure S5 for maps of the CFAs color-coded according to cluster
membership). The CFAs west of Mindoro island (#s 122 to 125) could also be considered as a separate
fisheries management area given their similarities in both seasonality of fishing activities (during
northeast monsoon) and radiance signature (mostly low radiance values). Of course, other fishing
activities and practices should also be considered in defining appropriate fisheries management areas
since the VBD only covers fishing vessels that use lights to attract fish. As Jennings and Lee [2] stated,
the ultimate delineation and definition of ‘fishing grounds’ would be a societal decision of choosing
cut-offs and satisfying potential trade-offs in the use of limited marine space.
The VBD data could also provide important information on the levels of fishing effort over time.
The estimated number of fishing vessels using lights in the Philippines from the VBD data (~1,000) is
most likely an underestimate since the VIIRS instrument takes only one image per night and for most
nights, there is always a part of the country that is obscured by thick cloud cover. Summing up the
maximum nightly VBD per CFA for selected larger, well-separated CFAs (i.e., #s 19, 66, 70, 75, 106,
122) already results in 925 estimated light boats or fishing vessels with on-board lights. In comparison,
there are an estimated 6,371 commercial fishing vessels (i.e., 3 gross tons and heavier) in the country as
of 2007 [79]. Despite this potential underestimation of absolute fishing effort, the quarterly pattern of
nighttime fishing lights still tracks closely the national statistics of landed catches of small pelagic fish
(Supplementary Figure S6). Monthly VBD data show a secondary peak between the third and fourth
quarter of each year which is not captured in the quarterly landed catch data. This secondary peak in
the VBD data is most likely associated with the varying seasonal patterns of fishing activities among
CFAs with some CFAs’ peak VBD counts occurring around this period (Figure 4). Subject to further
ground validation, the temporal trends in VBD data could be used as a proxy for relative fishing effort
at the national level for small pelagic fish catch, a key input for normalizing catch data.
in typical ocean conditions to most of the CFAs where lights are found (Figure 3; e.g., northeast tip
of Mindanao).
We ran the maximum entropy model under the assumption that areas where fishers often fish
(e.g., CFAs) are also areas of high relative abundance for the target species group, in this case, small pelagic
fishes, and that environmental factors can be used to predict habitat suitability [80]. At a resolution of 742 m,
it is surprising that many of the cells where nighttime lights are found in the last five years are revisited
over the years. We hypothesized that environmental factors played a large role in the spatial structure of
the nighttime lights, particularly for these frequently-visited and high-density sites. However, other factors
also influence fishers’ decisions in selecting fishing grounds such as projected fishing costs, estimated risk,
and adherence to traditional fishing grounds [50]. These are not captured in our model.
The utility of the VBD data in informing ecosystem approaches to fisheries management
could be enhanced by studies that can identify the relationships between radiance values and
fishing characteristics (e.g., fishing gear or strategy). This would help translate the VBD products
from remote sensing outputs to the more common fisheries metrics (e.g., number of vessels,
fishing gears, and/or species). In addition, the effect of cloud cover on detection probability and
limits should be investigated further especially when utilizing the VBD data for near real-time
applications. By aggregating the boat detection data annually or over several years, we minimized
the potential spatial bias brought about by spatial differences in the frequency of cloudy nights.
Supplementary Figure S7 in the supplementary file shows that lights are still detected even in areas
with high frequencies of cloudy nights. Despite these limitations, having nightly information
on distributions of light-assisted fishing activities freely accessible and available to fisheries
managers partially fills the crucial spatial data gap for implementing ecosystem approaches to
fisheries management.
5. Conclusions
Any marine resource management requires an understanding of the spatial and temporal patterns
of both resource use and target resources and the underlying processes. Yet collecting information
on spatial and temporal patterns of resource use at relevant scales for locally-adapted action remains
sparse, especially in marine fisheries in developing countries. We used the VIIRS boat detection
data service from NOAA to identify high density fishing grounds and predict areas within the
Philippines’s contiguous zone with similar environmental conditions in aid of implementing ecosystem
approaches to fisheries management. Nighttime lights are extensively detected in the Philippine coastal
waters. These lights spatially cluster and we have identified 134 distinct core fishing areas. Aggregating
all boat detections inside CFAs, these areas are strongly distinguished from non-CFA areas mainly
by bathymetry, owing most probably to the multi-gear nature of fishing and the varying fishing
operation strategies employed within each CFA. In two of the largest CFAs which we know are
each targeting round scads and sardines, variable importance in the MaxEnt models corresponded
with environmental parameters previously known to be associated with these species’ occurrence.
This highlights the potential model improvement if the catch composition and fishing fleets within
each CFA can be characterized. In addition, using finer spatial resolution environmental predictors
could enhance model precision and differentiate suitable fishing areas even within the same CFA.
The VIIRS boat detection data is provided freely and for public use. Its applications can be further
expanded by conducting ground-truthing and calibration studies to relate VBD radiance values to
common fishery statistics (e.g., number of vessels, types of vessels, comparative fishing effort based on
boat size, etc.). However, even with the current dataset, having nightly coverage of fishing activities,
even if just a subset of the entire fishing sector, provides an invaluable data stream of fishing effort
for data-deficient fisheries, especially in Southeast Asia. Our results could also be used as a basis for
establishing fisheries management areas in the Philippines and our approach could be applied to other
countries that have high numbers of light-employing fishing vessels.
Remote Sens. 2018, 10, 1604 19 of 23
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1, Figure S1: HDBSCAN
parameter combination; Figure S2: Marine traffic density map; Figure S3: Mean radiance versus percentage cloud
cover and moon illumination; Figure S4: Mean climatology maps of environmental predictors; Figure S5: Core
fishing areas mapped according to seasonal and radiance clusters; Figure S6: Comparison of VBD trend and
reported quarterly production of small pelagic fish in the Philippines; Figure S7. Spatial frequency of cloud cover
in relation to VIIRS boat detection data; Table S1: Characteristics of core fishing areas.
Author Contributions: Conceptualization, R.C.G., R.E.B., and C.M.; Data curation, C.D.E. and M.D.S.; Formal
analysis, R.C.G. and E.C.F.; Funding acquisition, R.E.B.; Methodology, R.C.G., E.C.F., C.D.E., M.D.S., R.V., and C.M.;
Supervision, C.M.; Writing—original draft, R.C.G.; Writing—review and editing, E.C.F., R.E.B., C.D.E., M.D.S.,
R.V., and C.M.
Acknowledgments: We thank Raffy Ramiscal and regional offices of the Philippine’s Bureau of Fisheries and
Aquatic Resources, the National Fisheries Research and Development Institute, and the National Stock Assessment
Program for comments on the manuscript and sharing their knowledge about nighttime fishing practices during
the various workshops conducted under the BFAR-USAID-NOAA Partnership Program. RCG was also supported
by the East-West Center. We thank Melanie Abecassis and three anonymous reviewers for their comments and
suggestions that improved our analyses and manuscript.
Conflicts of Interest: The authors declare no conflict of interest. The funders had no role in the design of the
study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to
publish the results.
References
1. Booth, A.J. Incorporating the spatial component of fisheries data into stock assessment models. ICES J. Mar. Sci.
2000, 57, 858–865. [CrossRef]
2. Jennings, S.; Lee, J. Defining fishing grounds with vessel monitoring system data. ICES J. Mar. Sci. 2012, 69,
51–63. [CrossRef]
3. Tidd, A.; Brouwer, S.; Pilling, G. Shooting fish in a barrel? Assessing fisher-driven changes in catchability
within tropical tuna purse seine fleets. Fish Fish. 2017, 1–13. [CrossRef]
4. Parnell, E.P.; Dayton, P.K.; Fisher, R.A.; Loarie, C.C.; Darrow, R.D. Spatial patterns of fishing effort off
San Diego: implications for zonal management and ecosystem function. Ecol. Appl. 2010, 20, 2203–2222.
[CrossRef] [PubMed]
5. Cabral, R.B.; Gaines, S.D.; Johnson, B.A.; Bell, T.W.; White, C. Drivers of redistribution of fishing and
non-fishing effort after the implementation of a marine protected area network. Ecol. Appl. 2017, 27, 416–428.
[CrossRef] [PubMed]
6. Walters, C. Folly and fantasy in the analysis of spatial catch rate data. Can. J. Fish. Aquat. Sci. 2003, 60,
1433–1436. [CrossRef]
7. Lorenzen, K.; Steneck, R.S.; Warner, R.R.; Parma, A.M.; Coleman, F.C.; Leber, K.M. The spatial dimension of
fisheries_putting it all in place.pdf. Bull. Mar. Sci. 2010, 86, 169–177.
8. Heenan, A.; Pomeroy, R.; Bell, J.; Munday, P.L.; Cheung, W.; Logan, C.; Brainard, R.; Yang Amri, A.;
Aliño, P.M.; Armada, N.B.; David, L.T.; Rivera-Guieb, R.; Green, S.; Jompa, J.; Leonardo, T.; Mamauag, S.S.;
Parker, B.; Shackeroff, J.; Yasin, Z. A climate-informed, ecosystem approach to fisheries management.
Mar. Policy 2015, 57, 182–192. [CrossRef]
9. Gorospe, K.D.; Michaels, W.; Pomeroy, R.; Elvidge, C.; Lynch, P.; Wongbusarakum, S.; Brainard, R.E.
The mobilization of science and technology fisheries innovations towards an ecosystem approach to fisheries
management in the Coral Triangle and Southeast Asia. Mar. Policy 2016, 74, 143–152. [CrossRef]
10. Stewart, K.R.; Lewison, R.L.; Dunn, D.C.; Bjorkland, R.H.; Kelez, S.; Halpin, P.N.; Crowder, L.B.
Characterizing fishing effort and spatial extent of coastal fisheries. PLoS ONE 2010, 5, e14451. [CrossRef]
[PubMed]
11. Kroodsma, D.A.; Mayorga, J.; Hochberg, T.; Miller, N.A.; Boerder, K.; Ferretti, F.; Wilson, A.; Bergman, B.;
White, T.D.; Block, B.A.; Woods, P.; Sullivan, B.; Costello, C.; Worm, B. Tracking the global footprint of
fisheries. Science 2018, 359, 904–908. [CrossRef] [PubMed]
12. Mccauley, B.D.J.; Woods, P.; Sullivan, B.; Bergman, B.; Jablonicky, C.; Roan, A.; Hirshfi, M.; Boerder, K.
Ending hide and seek at sea. Science 2016, 351, 1148–1150. [CrossRef] [PubMed]
13. Fett, R.W. Fishing Fleet Activities Revealed by Night-Time Data from the Defense Meteorological Satellite
Program. Mar. Fish. Rev. 1975, 47, 1972–1975.
Remote Sens. 2018, 10, 1604 20 of 23
14. Croft, T. Nighttime images of the earth from space. Sci. Am. 1978, 239, 86–98. [CrossRef]
15. Huang, Q.; Yang, X.; Gao, B.; Yang, Y.; Zhao, Y. Application of DMSP/OLS nighttime light images:
A meta-analysis and a systematic literature review. Remote Sens. 2014, 6, 6844–6866. [CrossRef]
16. Davies, T.W.; Duffy, J.P.; Bennie, J.; Gaston, K.J. The nature, extent, and ecological implications of marine
light pollution. Front. Ecol. Environ. 2014, 12, 347–355. [CrossRef]
17. Davies, T.W.; Duffy, J.P.; Bennie, J.; Gaston, K.J. Stemming the Tide of Light Pollution Encroaching into
Marine Protected Areas. Conserv. Lett. 2016, 9, 164–171. [CrossRef]
18. Elvidge, C.D.; Ziskin, D.; Baugh, K.E.; Tuttle, B.T.; Ghosh, T.; Pack, D.W.; Erwin, E.H.; Zhizhin, M. A fifteen
year record of global natural gas flaring derived from satellite data. Energies 2009, 2, 595–622. [CrossRef]
19. Elvidge, C.D.; Zhizhin, M.; Baugh, K.; Hsu, F.C.; Ghosh, T. Methods for global survey of natural gas flaring
from visible infrared imaging radiometer suite data. Energies 2016, 9, 14. [CrossRef]
20. Miller, S.D.; Straka, W.; Mills, S.P.; Elvidge, C.D.; Lee, T.F.; Solbrig, J.; Walther, A.; Heidinger, A.K.; Weiss, S.C.
Illuminating the capabilities of the suomi national Polar-orbiting partnership (NPP) visible infrared imaging
radiometer suite (VIIRS) day/night Band. Remote Sens. 2013, 5, 6717–6766. [CrossRef]
21. Elvidge, C.D.; Zhizhin, M.; Baugh, K.; Hsu, F.C. Automatic boat identification system for VIIRS low light
imaging data. Remote Sens. 2015, 7, 3020–3036. [CrossRef]
22. Straka, W.C.; Seaman, C.J.; Baugh, K.; Cole, K.; Stevens, E.; Miller, S.D. Utilization of the suomi national
polar-orbiting partnership (npp) visible infrared imaging radiometer suite (viirs) day/night band for arctic
ship tracking and fisheries management. Remote Sens. 2015, 7, 971–989. [CrossRef]
23. Waluda, C.M.; Trathan, P.N.; Elvidge, C.D.; Hobson, V.R.; Rodhouse, P.G. Throwing light on straddling
stocks of Illex argentinus: assessing fishing intensity with satellite imagery. Can. J. Fish. Aquat. Sci. 2002, 59,
592–596. [CrossRef]
24. Kiyofuji, H.; Saitoh, S.I. Use of nighttime visible images to detect Japanese common squid Todarodes
pacificus fishing areas and potential migration routes in the Sea of Japan. Mar. Ecol. Prog. Ser. 2004, 276,
173–186. [CrossRef]
25. Waluda, C.M.; Yamashiro, C.; Elvidge, C.D.; Hobson, V.R.; Rodhouse, P.G. Quantifying light-fishing for
Dosidicus gigas in the eastern Pacific using satellite remote sensing. Remote Sens. Environ. 2004, 91, 129–133.
[CrossRef]
26. Saitoh, S.; Fukaya, A.; Saitoh, K.; Semedi, B. Estimation of number of Pacific saury fishing vessels using
night-time visible images. Int. Arch. Photogramm. Remote Sens. Spatial Inf. Sci. 2010, 1013–1016.
27. Liu, Y.; Saitoh, S.I.; Hirawake, T.; Igarashi, H.; Ishikawa, Y. Detection of Squid and a Pacific Saury fishing vessels
around Japan using VIIRS Day/Night Band data. Proc. Asia-Pacific Adv. Netw. 2015, 39, 28–39. [CrossRef]
28. Alabia, I.D.; Saitoh, S.; Mugo, R.; Igarashi, H.; Ishikawa, Y.; Usui, N.; Kamachi, M.; Awaji, T. Seasonal
potential fishing ground prediction of neon flying squid (Ommastrephes bartramii) in the western and central
North Pacific. Fish. Oceanogr. 2015, 190–203. [CrossRef]
29. Syah, A.F.; Saitoh, S.I.; Alabia, I.D.; Hirawake, T. Predicting potential fishing zones for pacific saury
(Cololabis saira) with maximum entropy models and remotely sensed data. Fish. Bull. 2016, 114, 330–342.
[CrossRef]
30. Elvidge, C.D.; Ghosh, T.; Baugh, K.; Zhizhin, M. Rating the Effectiveness of Fishery Closures With Visible
Infrared Imaging Radiometer Suite Boat Detection Data. Front. Mar. Sci. 2018, 5, 1–15. [CrossRef]
31. Cabral, R.B.; Mayorga, J.; Clemence, M.; Lynham, J.; Koeshendrajana, S.; Muawanah, U.; Nugroho, D.;
Anna, Z.; Mira; Ghofar, A.; et al. Rapid and lasting gains from solving illegal fishing. Nat. Ecol. Evol. 2018, 2,
650–658. [CrossRef] [PubMed]
32. Elvidge, C.D.; Baugh, K.; Zhizhin, M.; Hsu, F.C. Why VIIRS data are superior to DMSP for mapping nighttime
lights. Proc. Asia-Pacific Adv. Netw. 2013, 35, 62–69. [CrossRef]
33. Macusi, E.D.; Babaran, R.P.; van Zwieten, P.A.M. Strategies and tactics of tuna fishers in the payao (anchored
FAD) fishery from general Santos city, Philippines. Mar. Policy 2015, 62, 63–73. [CrossRef]
34. Geronimo, R.C. Small Fish, Big Impact: Dulong fisheries of San Juan, Batangas, Philippines: A Synthesis Report;
Conservation International: Quezon City, Philippines, 2013.
35. Badan Pusat Statistik. Statistics Indonesia Statistics of Marine and Coastal Resources; BPS-Statistics Indonesia:
Jakarta, Indonesia, 2015; p. xviii + 217.
Remote Sens. 2018, 10, 1604 21 of 23
36. BFAR. Philippine Fisheries Profile 2016; Department of Agriculture – Bureau of Fisheries and Aquatic Resources:
Quezon City, Philippines, 2017; p. 69.
37. Dickson, J.O.; Alba, E.B.; Munprasit, A.; Chokesanguan, B.; Siriraksophon, S. Fishing Gear and Methods in
Southeast Asia: III. Philippines Part 2; Ruangsivakul, N., Prajakjitt, P., Dickson, J.O., Siriraksophon, S., Eds.;
Southeast Asian Fisheries Development, Training Department: Samutprakan, Thailand, 2003.
38. Hernando, A.M.J.; Flores, E.E.C. The Philippines squid fishery: A review. Mar. Fish. Rev. 1981, 43, 13–20.
39. Zaragoza, E.; Pagdilao, C.; Moreno, E. Overview of small pelagic fisheries. In In Turbulent Seas: The status of
Philippine marine fisheries; DA-BFAR (Department of Agriculture-Bureau of Fisheries and Aquatic Resources),
Ed.; Coastal Resource Management Project: Cebu City, Philippines, 2004; p. 378.
40. Claus, S.; De Hauwere, N.; Vanhoorne, B.; Souza Dias, F.; Oset García, P.; Schepers, L.; Hernandez, F.; Mees, J.;
(Flanders Marine Institute). MarineRegions.org. Available online: http://www.marineregions.orgon2016-08-05.
2016 (accessed on 20 January 2018).
41. Jennings, S.; Dulvy, N.K. Reference points and reference directions for size-based indicators of community
structure. ICES J. Mar. Sci. 2005, 62, 397–404. [CrossRef]
42. Saul, S.E.; Walter, J.F.; Die, D.J.; Naar, D.F.; Donahue, B.T. Modeling the spatial distribution of commercially
important reef fishes on the West Florida Shelf. Fish. Res. 2013, 143, 12–20. [CrossRef]
43. Campello, R.J.G.B.; Moulavi, D.; Sander, J. Density-Based Clustering Based on Hierarchical Density Estimates.
In Advances in Knowledge Discovery and Data Mining. PAKDD 2013. Lecture Notes in Computer Science;
Pei, J., Tseng, V.S., Cao, L., Motoda, H., Xu, G., Eds.; Springer: Berlin, Germany, 2013; Volume 7819,
ISBN 978-3-642-37455-5.
44. R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical
Computing, Vienna, Austria. Available online: https://www.r-project.org/ (accessed on 20 January 2018).
45. McInnes, L.; Healy, J.; Astels, S. hdbscan: Hierarchical density based clustering. J. Open Source Softw. 2017,
2, 205. [CrossRef]
46. Cao, C.; Bai, Y. Quantitative analysis of VIIRS DNB nightlight point source for light power estimation and
stability monitoring. Remote Sens. 2014, 6, 11915–11935. [CrossRef]
47. Irpino, A.; Verde, R. A New Wasserstein Based Distance for the Hierarchical Clustering of Histogram Symbolic
Data. In Data Science and Classification. Studies in Classification, Data Analysis, and Knowledge Organization; Batagelj, V.,
Bock, H., Ferligoj, A., Žiberna, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2006.
48. Irpino, A. HistDAWass: Histogram-Valued Data Analysis. R Package Version 1.0.1. 2018. Available online:
https://CRAN.R-project.org/package=HistDAWass (accessed on 20 January 2018).
49. Murtagh, F.; Legendre, P. Ward’s Hierarchical Agglomerative Clustering Method: Which Algorithms
Implement Ward’s Criterion? J. Classif. 2014, 31, 274–295. [CrossRef]
50. Girardin, R.; Hamon, K.G.; Pinnegar, J.; Poos, J.J.; Thébaud, O.; Tidd, A.; Vermard, Y.; Marchal, P. Thirty
years of fleet dynamics modelling using discrete-choice models: What have we learned? Fish Fish. 2017, 18,
638–655. [CrossRef]
51. Phillips, S.J.; Dudík, M.; Schapire, R.E. [Internet] Maxent Software for Modeling Species Niches and Distributions
(Version 3.4.1). Available online: http://biodiversityinformatics.amnh.org/open_source/maxent/ (accessed on
4 August 2018).
52. Phillips, S.J.; Dudík, M. Modeling of species distribution with Maxent: New extensions and a comprehensive
evalutation. Ecograpy 2008, 31, 161–175. [CrossRef]
53. Morales, N.S.; Fernández, I.C.; Baca-González, V. MaxEnt’s parameter configuration and small samples:
Are we paying attention to recommendations? A systematic review. PeerJ 2017, 5, e3093. [CrossRef]
[PubMed]
54. Robinson, N.M.; Nelson, W.A.; Costello, M.J.; Sutherland, J.E.; Lundquist, C.J. A Systematic Review of
Marine-Based Species Distribution Models (SDMs) with Recommendations for Best Practice. Front. Mar. Sci.
2017, 4, 1–11. [CrossRef]
55. Elith, J.; Phillips, S.J.; Hastie, T.; Dudík, M.; Chee, Y.E.; Yates, C.J. A statistical explanation of MaxEnt for
ecologists. Divers. Distrib. 2011, 17, 43–57. [CrossRef]
56. Muscarella, R.; Galante, P.J.; Soley-Guardia, M.; Boria, R.A.; Kass, J.M.; Uriarte, M.; Anderson, R.P. ENMeval:
An R package for conducting spatially independent evaluations and estimating optimal model complexity
for Maxent ecological niche models. Methods Ecol. Evol. 2014, 5, 1198–1205. [CrossRef]
Remote Sens. 2018, 10, 1604 22 of 23
57. Merino, G.; Barange, M.; Mullon, C. Climate variability and change scenarios for a marine commodity:
Modelling small pelagic fish, fisheries and fishmeal in a globalized market. J. Mar. Syst. 2010, 81, 196–205.
[CrossRef]
58. Zwolinski, J.P.; Emmett, R.L.; Demer, D.A. Predicting habitat to optimize sampling of Pacific sardine
(Sardinops sagax). ICES J. Mar. Sci. 2011, 68, 867–879. [CrossRef]
59. Peck, M.A.; Reglero, P.; Takahashi, M.; Catalán, I.A. Life cycle ecophysiology of small pelagic fish and
climate-driven changes in populations. Prog. Oceanogr. 2013, 116, 220–245. [CrossRef]
60. Kaplan, I.C.; Williams, G.D.; Bond, N.A.; Hermann, A.J.; Siedlecki, S.A. Cloudy with a chance of sardines:
Forecasting sardine distributions using regional climate models. Fish. Oceanogr. 2016, 25, 15–27. [CrossRef]
61. Hijmans, R.J. raster: Geographic Data Analysis and Modeling. R package version 2.6-7. 2017. Available
online: https://CRAN.R-project.org/package=raster (accessed on 20 January 2018).
62. Thuiller, W.; Georges, D.; Engler, R.; Breiner, F. biomod2: Ensemble Platform for Species Distribution Modeling. R
package version 3.3-15/r728. Available online: https://R-Forge.R-project.org/projects/biomod/2017 (accessed
on 20 January 2018).
63. Trinidad, A.C.; Pomeroy, R.S.; Corpuz, P.V.; Aguero, M. Bioeconomics of the Philippine small pelagics fishery;
ICLARM: Los Banos, Philippines, 1993; ISBN 971-8709-38-X.
64. Elith, J.; Leathwick, J.R. Species Distribution Models: Ecological Explanation and Prediction Across Space
and Time. Annu. Rev. Ecol. Evol. Syst. 2009, 40, 677–697. [CrossRef]
65. Elith, J.; Graham, C.H.; Anderson, R.P.; Dudík, M.; Ferrier, S.; Guisan, A.; Hijmans, R.J.; Huettmann, F.;
Leathwick, J.R.; Lehmann, A.; Li, J.; Lohmann, L.G.; Loiselle, B.A.; Manion, G.; Moritz, C.; Nakamura, M.;
Nakazawa, Y.; Overton, J.M.; Townsend Peterson, A.; Phillips, S.J.; Richardson, K.; Scachetti-Pereira, R.;
Schapire, R.E.; Soberón, J.; Williams, S.; Wisz, M.S.; Zimmermann, N.E. Novel methods improve prediction
of species’ distributions from occurrence data. Ecography. 2006, 29, 129–151. [CrossRef]
66. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions.
Ecol. Modell. 2006, 190, 231–259. [CrossRef]
67. Thuiller, W.; Lafourcade, B.; Engler, R.; Araújo, M.B. BIOMOD—A platform for ensemble forecasting of
species distributions. Ecography. 2009, 32, 369–373. [CrossRef]
68. Landis, J.R.; Koch, G.G. The Measurement of Observer Agreement for Categorical Data. Biometrics 1977, 33,
159. [CrossRef] [PubMed]
69. Rola, A.C.; Narvaez, T.A.; Naguit, M.R.A.; Elazegui, D.D.; Brillo, B.B.C.; Paunlagui, M.M.; Jalotjot, H.C.;
Cervantes, C.P. Impact of the closed fishing season policy for sardines in Zamboanga Peninsula, Philippines.
Mar. Policy 2018, 87, 40–50. [CrossRef]
70. Diego, T.J.A.S.; Fisher, W.L. Trends in the capture fisheries in Cuyo East Pass, Philippines. Int. J. Fish. Aquat. Stud.
2014, 1, 57–72.
71. Zhao, X.; Li, D.; Li, X.; Zhao, L.; Wu, C. Spatial and seasonal patterns of night-time lights in global ocean
derived from VIIRS DNB images. Int. J. Remote Sens. 2018, 1161, 1–31. [CrossRef]
72. Nurdin, S.; Mustapha, M.A.; Lihan, T. The Relationship between Sea Surface Temperature and Chlorophyll-
a Concentration in Fisheries Aggregation Area in the Archipelagic Waters of Spermonde Using Satellite
Images. In AIP Conference Proceedings; AIP Publishing: Melville, NY, USA, 2014; Volume 1571, pp. 466–472.
73. Bellido, J.M.; Brown, A.M.; Valavanis, V.D.; Giráldez, A.; Pierce, G.J.; Iglesias, M.; Palialexis, A. Identifying
essential fish habitat for small pelagic species in Spanish Mediterranean waters. Hydrobiologia 2008, 612,
171–184. [CrossRef]
74. Schismenou, E.; Tsoukali, S.; Giannoulaki, M.; Somarakis, S. Modelling small pelagic fish potential spawning
habitats: Eggs vs spawners and in situ vs satellite data. Hydrobiologia 2017, 788, 17–32. [CrossRef]
75. Tugores, P.; Giannoulaki, M.; Iglesias, M.; Bonanno, A.; Tičina, V.; Leonori, I.; MacHias, A.; Tsagarakis, K.;
Díaz, N.; Giráldez, A.; Patti, B.; De Felice, A.; Basilone, G.; Valavanis, V. Habitat suitability modelling for
sardine Sardina pilchardus in a highly diverse ecosystem: The Mediterranean Sea. Mar. Ecol. Prog. Ser. 2011,
443, 181–205. [CrossRef]
76. Kaschner, K.; Kesner-Reyes, K.; Garilao, C.; Rius-Barile, J.; Rees, T.; Froese, R. AquaMaps: Predicted range
maps for aquatic species. Available online: https://www.aquamaps.org/ (accessed on 20 January 2017).
77. Kesner-Reyes, K.; Kaschner, K.; Kullander, S.; Garilao, C.; Baril, J.; Froese, R. AquaMaps: Algorithm and
data sources for aquatic organisms. Available online: https://www.aquamaps.org/am.pdf (accessed on
20 January 2017).
Remote Sens. 2018, 10, 1604 23 of 23
78. Saraux, C.; Fromentin, J.M.; Bigot, J.L.; Bourdeix, J.H.; Morfin, M.; Roos, D.; Van Beveren, E.; Bez, N. Spatial
Structure and Distribution of Small Pelagic Fish in the Northwestern Mediterranean Sea. PLoS ONE 2014,
9, e111211. [CrossRef] [PubMed]
79. BFAR. Philippine Fisheries Profile 2012; Department of Agriculture – Bureau of Fisheries and Aquatic Resources:
Quezon City, Philippines, 2013; p. 69. Available online: https://www.bfar.da.gov.ph/files/img/photos/
2012FisheriesProfile(Finalcopy)(1).pdf (accessed on 20 January 2017).
80. Alabia, I.D.; Saitoh, S.I.; Hirawake, T.; Igarashi, H.; Ishikawa, Y.; Usui, N.; Kamachi, M.; Awaji, T.; Seito, M.
Elucidating the potential squid habitat responses in the central North Pacific to the recent ENSO flavors.
Hydrobiologia 2016, 772, 215–227. [CrossRef]
© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access
article distributed under the terms and conditions of the Creative Commons Attribution
(CC BY) license (http://creativecommons.org/licenses/by/4.0/).