Remotely Sensing the Biophysical Drivers of Sardinella aurita Variability in Ivorian Waters

The coastal regions of the Gulf of Guinea constitute one of the major marine ecosystems, producing essential living marine resources for the populations of Western Africa. In this region, the Ivorian continental shelf is under pressure from various anthropogenic sources, which have put the regional fish stocks, especially Sardinella aurita, the dominant pelagic species in Ivorian industrial fishery landings, under threat from overfishing. Here, we combine in situ observations of Sardinella aurita catch, temperature, and nutrient profiles, with remote-sensing ocean-color observations, and reanalysis data of wind and sea surface temperature, to investigate relationships between Sardinella aurita catch and oceanic primary producers (including biomass and phenology of phytoplankton), and between Sardinella aurita catch and environmental conditions (including upwelling index, and turbulent mixing). We show that variations in Sardinella aurita catch in the following year may be predicted, with a confidence of 78%, based on a bilinear model using only physical variables, and with a confidence of 40% when using only biological variables. However, the physics-based model alone is not sufficient to explain the mechanism driving the year-to-year variations in Sardinella aurita catch. Based on the analysis of the relationships between biological variables, we demonstrate that in the Ivorian continental shelf, during the study period 1998–2014, population dynamics of Sardinella aurita, and oceanic primary producers, may be controlled, mainly by top-down trophic interactions. Finally, based on the predictive models constructed here, we discuss how they can provide powerful tools to support evaluation and monitoring of fishing activity, which may help towards the development of a Fisheries Information and Management System.


Introduction
Upwelling ecosystems provide more than 40% of the world fisheries catch, although they occupy less than 3% of the ocean's surface [1].The coastal regions of the Gulf of Guinea (GG) constitute one of the major marine ecosystems, producing essential living marine resources for the populations of Western Africa [2].Among these coastal ecosystems, that of Côte d'Ivoire, and its upwelling, presents unusual characteristics compared with the upwelling systems of the bordering regions of Mauritania in the North, and Central Guinea in the South [3].The originality of the Ivorian upwelling system lies in the existence of two well-differentiated cooling seasons [4,5]: a minor upwelling during the months of January to February, followed by a major upwelling during the months of June to September.The major upwelling is characterized by lower sea surface temperature (SST), between 22 and 23 • C, higher chlorophyll (Chl) concentration, between 1.5 and 3 mg Chl m −3 , and strong winds, inducing turbulence that reaches up to 350 m 3 s −3 in the month of July [6].
Since the 1980s, the Ivorian continental shelf has been subject to strong anthropogenic pressures associated with increased population size, as well as economic development, urbanization, wastewater discharge, and other localized pollution, habitat degradation, unsustainable fishing practices, and global climate change.Cumulatively, these changes have put the regional fish stocks (especially Sardinella) under threat from overfishing [7,8].Consequently, an assessment of the regional risks and vulnerabilities of the marine environment is essential for responsible management of fishery resources.This requires that the biological and physical variables influencing fishery resources be identified.
Understanding recruitment mechanisms in marine fish populations is a major challenge in fisheries research.The importance of oceanic processes has become increasingly apparent in determining recruitment levels [9,10].However, it is the recent advancement of remote-sensing technology, providing biological and physical observations with sufficient spatial and temporal resolution, which has permitted us to verify fish recruitment mechanisms (some proposed more than a century ago), as well as to explore operational metrics for ecosystem-based assessment and management frameworks for the resource concerned [11,12].
Several hypotheses based on food availability and predation have been proposed to explain the dominant mechanisms driving recruitment variability [13].First, Hjort [14] suggested that the larval feeding phase is a critical life period, and that its coincidence with high prey abundance is a necessary condition for survival and recruitment success of the fish larvae (i.e., so-called, the "critical period" hypothesis).Lasker [15,16] provided experimental and field evidence in support of the critical period hypothesis for the northern anchovy (Engraulis mordax), and introduced the concept of stable ocean conditions, relating larval fish food aggregation to stability of the water column and its variations in relation to oceanographic and meteorological conditions (i.e., so-called the "stable ocean" hypothesis; [17]).Later, Cushing [18] reported that a match or mismatch, between the timing of peak spawning activity of fish and the timing of phytoplankton bloom, could be a major determinant of interannual variability in fish larvae survival (i.e., the "match-mismatch" hypothesis).Using remote-sensing technology, Platt et al. [11] have made operational tests of the match-mismatch hypothesis, and shown that a high proportion of the variance between years in the larval recruitment of haddock (Melanogrammus aeglefinus) off Nova Scotia could be accounted for by local anomalies in the phenology of phytoplankton.Also in the North Atlantic Ocean, Koeller et al. [19] have shown that the recruitment of a commercially important crustacean (Pandalus borealis) is tightly coupled to the timing of the phytoplankton bloom.This coupling is observed because the early stages of fish and crustacean larvae feed primarily on phytoplankton.Thus, if phytoplankton (i.e., source of food) is not available when the larvae are released, their survival rate decreases, due to starvation during their most vulnerable stage.By contrast, other hypotheses, such as that proposed by Iles and Sinclair [20], and then Sinclair [21], argued that physics predominates over food chain processes in the control of fish populations (i.e., the "member/vagrant" hypothesis).In the latter hypothesis, maintenance of the spatial integrity of any marine fish population is regarded as the most important factor, and spawning times are adapted to the physical dispersive characteristics of an area, rather than to the biological production characteristics.Retention of early life history stages in favorable oceanographic conditions is fundamental for survival.However, there is evidence in the literature of both dispersal and food chain processes influencing year-class survival [13].Lasker [17] concluded that physical processes may affect the characteristics of the food chain, such that long periods of stable oceanic conditions lead to fine scale particle aggregations, enhancing food availability and survival of larval stages.Finally, Bakun [22] reported that the spawning peak of pelagic fishes in upwelling areas coincides with a seasonal minimum in offshore Ekman transport.
Like any coastal upwelling ecosystems, the Ivorian marine environment is favorable to the development of phytoplankton populations, whose biomass can be characterized by their chlorophyll concentrations [23][24][25][26][27][28].To date, most studies on phytoplankton concentration in Ivorian waters have been based on in situ measurements [23][24][25][26][29][30][31], with sparse and discontinuous coverage, due to limited cruise sampling frequency.To our knowledge, only few studies have used satellite observations of ocean color, wind, and SST to monitor and explain phytoplankton variations in the region [28].Maximum biomass of phytoplankton is found near the coast, and decreases gradually towards the open ocean [27,28,30,[32][33][34].The high phytoplankton concentrations observed at the two capes (Cap des Palmes at the Liberian border and Cap des Trois-Pointes at the Ghanaian border) are attributed to the effect of the Guinea current, which creates a cyclonic eddy that continuously renews the waters at these two locations [35].Between the two capes, Roy et al. [36] and Koné et al. [37] have reported high abundances of pelagic fish populations, including Sardinella aurita (S. aurita).
The species S. aurita is known to live in shoals, which are highly migratory, and which represent an important component of the tertiary biomass of the Ivorian continental shelf.This species breeds throughout the year, with marked seasonal peaks [36]: South of Cape Verde, it breeds mainly during the upwelling period with a major peak in May and June [38].According to Oliver-Massuti [39], Sardinella follows areas of higher planktonic density.Komarovsky [40] has concluded similarly for S. aurita along the Israeli coasts in the Mediterranean Sea.Interannual variability in catches of Sardinella appears to be closely linked to the variations of chlorophyll concentration off the continental shelf of Côte d'Ivoire [6,41,42].Highest catches are generally reported during the cold season between June and September [43], when concentrations of chlorophyll pigments are highest.
Several studies [7,32,34,44] have provided information on the retention and recruitment areas for S. aurita in the marine environment of Côte d'Ivoire using physical variables.However, the information accumulated, so far, does not provide appropriate decision-support tools to fishermen, nor to the department for fishing regulation (DPH), especially in the estimation and forecast of the temporal variation of this species.The aim of the present work is to advance knowledge of the relationships between S. aurita and oceanic primary producers (including phytoplankton phenology, and phytoplankton biomass) as the main source of food for the fish larvae, and between S. aurita and environmental indices (i.e., upwelling index and wind-induced turbulent mixing).These relationships are then used to construct diagnostic models and examine the predictability of year-to-year variations in S. aurita catch along the Ivorian continental shelf.Finally, we explore how these models may be used to inform fisheries management, and discuss their potential use in an operational context.

Ocean Color Remote Sensing Data
The level-3 data of chlorophyll-a concentration in seawater at 4 × 4 km and 8 day resolution were obtained from the European Space Agency Ocean Color Climate Change Initiative ESA OC-CCI project at https://rsg.pml.ac.uk/thredds/catalog-cci.html for the period September 1997-December 2014 [45].Version 3 of the OC-CCI chlorophyll data product comprising merged and bias-corrected MERIS, MODIS-Aqua, VIIRS, and SeaWiFS sensor observations was downloaded.The chlorophyll product had been generated using a blended combination of OCI, OC3, and OC5 retrieval algorithms, depending on water class memberships [46].The remote-sensing reflectance data for use in the chlorophyll algorithms, are generated using SeaDAS atmospheric-correction algorithm [47] for SeaWiFS and VIIRS, and POLYMER atmospheric-correction algorithm [48] for MERIS and MODIS-Aqua.The merging of multiple sensor observations and the improvement in the atmospheric-correction algorithm, POLYMER, have permitted us, for the first time, to enhance significantly the data coverage in Ivorian waters, from ~35% data coverage with single sensor SeaWiFS or MODIS to ~95% data coverage with OC-CCI merged product ( [28], their Figure 4).Furthermore, reduction in the number of missing data has been shown to be important to improve accuracy and precision in the estimation of phenological indices [49].

Sardinella Fish Catch Dataset
Sardinella aurita data used in this study have been obtained for the period from 1998 to 2014 from the statistical and documentation department of the Abidjan Fisheries Direction.The data are available in the yearbook of fishery statistics produced for the Ministry of Production of Animals (MPA).The landings of S. aurita are available for the fishing zone of Abidjan.The S. aurita landings data from the Abidjan Fisheries Direction have been compared with the global capture of round sardines (S. aurita) from the Food and Agriculture Organization (FAO).The latter data were downloaded from http://www.fao.org/fishery/statistics/global-capture-production/en for the waters of Côte d'Ivoire over the period 1998-2014.The two datasets are plotted in Figure 1.The coefficient of regression is estimated to r = 0.92, p < 0.0001, which suggests that the interannual variability in the landings of S. aurita observed in the fishing zone of Abidjan (see box region location in Figure 2) is representative of the landings for the whole of the Ivorian waters.FAO datasets have been examined before in studies focusing on the links between fisheries trends and biophysical/climate factors (e.g., [50][51][52][53]).Trends in stock abundances are usually detected by analyzing catch or landings per unit of effort (CPUE or LPUE), rather than pure catch or landings data [54].The main reason is that un-normalized fishery landings do not necessarily represent stock abundance, and that fluctuations in these datasets could be also influenced by other factors.However, it has been argued that trends in catch data are consistent with biomass trends [55].In our study, the use of two different forms of landing data, which are consistent with each other, and with one of them representing a broader region, guarantees that the patterns detected are not artefacts produced by the sampling scheme or changes in manner of recording by a specific country.~95% data coverage with OC-CCI merged product ( [28], their Figure 4).Furthermore, reduction in the number of missing data has been shown to be important to improve accuracy and precision in the estimation of phenological indices [49].

Sardinella Fish Catch Dataset
Sardinella aurita data used in this study have been obtained for the period from 1998 to 2014 from the statistical and documentation department of the Abidjan Fisheries Direction.The data are available in the yearbook of fishery statistics produced for the Ministry of Production of Animals (MPA).The landings of S. aurita are available for the fishing zone of Abidjan.The S. aurita landings data from the Abidjan Fisheries Direction have been compared with the global capture of round sardines (S. aurita) from the Food and Agriculture Organization (FAO).The latter data were downloaded from http://www.fao.org/fishery/statistics/global-capture-production/en for the waters of Côte d'Ivoire over the period 1998-2014.The two datasets are plotted in Figure 1.The coefficient of regression is estimated to r = 0.92, p < 0.0001, which suggests that the interannual variability in the landings of S. aurita observed in the fishing zone of Abidjan (see box region location in Figure 2) is representative of the landings for the whole of the Ivorian waters.FAO datasets have been examined before in studies focusing on the links between fisheries trends and biophysical/climate factors (e.g., [50][51][52][53]).Trends in stock abundances are usually detected by analyzing catch or landings per unit of effort (CPUE or LPUE), rather than pure catch or landings data [54].The main reason is that unnormalized fishery landings do not necessarily represent stock abundance, and that fluctuations in these datasets could be also influenced by other factors.However, it has been argued that trends in catch data are consistent with biomass trends [55].In our study, the use of two different forms of landing data, which are consistent with each other, and with one of them representing a broader region, guarantees that the patterns detected are not artefacts produced by the sampling scheme or changes in manner of recording by a specific country.

Wind Data
The ERA Interim reanalysis products of zonal velocity U (i.e., the component of the horizontal wind positive for a west to east flow), and meridional velocity V (i.e., the component of the horizontal wind positive for a south to north flow) were obtained from the European Centre for Medium-Range Weather Forecast (ECMWF) at http://apps.ecmwf.int/datasets/data/interim-full-moda/levtype=sfc/at 12.5 × 12.5 km and monthly resolution over the period 1998-2014.

Sea-Surface Temperature Data
The Sea Surface Temperature (SST) data of the ERA Interim reanalysis were obtained from the European Centre for Medium-Range Weather Forecast (ECMWF) at http://apps.ecmwf.int/datasets/data/interim-full-moda/levtype=sfc/at 12.5 × 12.5 km and monthly resolution over the period 1998-2014.

In-Situ Data of Temperature and Nutrient Concentrations
Monthly climatologies of temperature, concentrations of nitrate and phosphate were obtained at 1° × 1° resolution from the global oceans from the World Ocean Atlas Climatology [56] at https://www.nodc.noaa.gov/OC5/woa13/.

Wind Data
The ERA Interim reanalysis products of zonal velocity U (i.e., the component of the horizontal wind positive for a west to east flow), and meridional velocity V (i.e., the component of the horizontal wind positive for a south to north flow) were obtained from the European Centre for Medium-Range Weather Forecast (ECMWF) at http://apps.ecmwf.int/datasets/data/interim-full-moda/levtype=sfc/at 12.5 × 12.5 km and monthly resolution over the period 1998-2014.

Sea-Surface Temperature Data
The Sea Surface Temperature (SST) data of the ERA Interim reanalysis were obtained from the European Centre for Medium-Range Weather Forecast (ECMWF) at http://apps.ecmwf.int/datasets/data/interim-full-moda/levtype=sfc/ at 12.5 × 12.5 km and monthly resolution over the period 1998-2014.

In-Situ Data of Temperature and Nutrient Concentrations
Monthly climatologies of temperature, concentrations of nitrate and phosphate were obtained at 1 • × 1 • resolution from the global oceans from the World Ocean Atlas Climatology [56] at https: //www.nodc.noaa.gov/OC5/woa13/.

Estimation of Phenological Indices
The phenological indices for timings of initiation, peak, termination, and duration of the main phytoplankton-growing period have been estimated each year at each 4 km × 4 km pixel using the method described in Racault et al. [57].The amplitude and timing of the chlorophyll peak concentration were estimated between the months of April to October of year t.These two months usually correspond to the time of lowest chlorophyll concentration in the seasonal cycle of chlorophyll in the region (Figure 2).The threshold criterion of long-term (1997-2014) median chlorophyll concentration plus 20% was calculated.It was tuned to account for the seasonal pattern of variations in chlorophyll concentration (which may vary by up to ~one order of magnitude in some years) in Ivorian waters.We carried out sensitivity analyses and observed that decreasing (increasing) the threshold from 5% to 30% tends to advance (delay) the timing of initiation.Similar results were reported in sensitivity analyses of threshold criteria carried out in the North Atlantic [58] and global oceans [57].The threshold of median plus 20% was shown to be most representative of the interannual variations in timing of initiation.Starting from the timing of peak chlorophyll, the timings of initiation and termination were estimated as the times when the chlorophyll concentration fell below the threshold, on either side of the peak.The duration was calculated as the number of days between the timings of initiation and termination.

Estimation of Upwelling Index (UI)
The upwelling index was calculated following Bakun [59].The displacement of surface waters forced by winds, and deflected by the Coriolis force, is effective only to a certain depth.The integrated flow between the ocean surface and this depth of "limited" motion is known as Ekman transport, and it moves perpendicularly to the wind direction [60].The study of Bakun [59], based on the previous study of Ekman [60], defined the "upwelling index".This index estimates the volume of upwelled flow, and its transport depends on latitude and wind stress.Thus, the upwelling index summarizes the oceanwards flow of surface waters per km of coastline (in m 3 s −1 km −1 ), in the west-east (Q x ) and north-south (Q y ) directions.The Ivorian coast presents an orientation west-east, hence the equation in the direction west-east (Q x ) is used in this study, such that where τ y is the wind stress component estimated from the wind speed, as follows: where ρ a = 1.2 kg m −3 (air density), C d = 1.4 × 10 −3 (empirical resistance coefficient), and u and v are wind components in m•s −1 .

Estimation of Wind-Induced Turbulent Mixing
A proxy for turbulence was calculated as the cube of wind stress [9].This provides information on wind-induced turbulent mixing, and gives an idea of vertical stability of the sea.

Seasonal Variability of Biophysical Drivers in Ivorian Waters
The seasonal variations of chlorophyll pigment concentration and SST in the fishing zone of Abidjan are shown in Figure 2a,b and Figure 3a.High phytoplankton blooms with chlorophyll concentrations of ~1.7 mg m −3 are observed in coastal waters during the period from July to September, which corresponds to the period of intense upwelling and low SST.Phytoplankton variability follows seasonal fluctuations in SST.
The SST value increases from 27.5 to 29 • C from January to April, followed by a gradual decrease from April to August, to a minimum value of 24 • C. Finally, the SST increases until December to reach the maximum value of 28.5 • C. We observe that cold water is present over the entire Ivorian coastal zone during the major cold season, with an average temperature of 24 • C.During this season, the coldest surface temperatures are observed in the central and eastern parts of the littoral between 5 and 2 • W during the month of August.Warmer temperatures are found in the western part between longitudes 7 • W and 8 • W during the month of July.
The seasonal variation in coastal upwelling in the fishing zone of Abidjan is shown in Figure 2c.The upwelling index increases from 500 m 3 s −1 km −1 in January to 1100 m 3 s −1 km −1 in March, then decreases slightly in May to 750 m 3 s −1 km −1 .From June, there is a rapid and maximal increase in upwelling index, reaching values of 2000 m 3 s −1 km −1 in July-August.Finally, values decrease gradually to reach 500 m 3 s −1 km −1 in December.In the fishing zone of Abidjan, wind-induced turbulent mixing (Figure 2d) follows the same seasonal variations as the upwelling index.Turbulent mixing increases from 10 m 3 s −3 in January, reaching 35 m 3 s −3 in March-April-May.Then, turbulent mixing values increase again, reaching a maximum of 80 m 3 s −3 in July-August, before gradually decreasing to values of 10 m 3 s −3 in December.The variations in coastal upwelling index and turbulent mixing are consistent with in situ measurements of thermocline depth, showing shallow depths between 0 and 50 m (upper and lower bounds, respectively) from July to September during the major cold season, and between 10 and 60 m from January to March during the minor cold season.The thermocline deepens to between 15 and 65 m (upper and lower bounds, respectively), on average, between April and June, and between October and December.
The SST value increases from 27.5 to 29 °C from January to April, followed by a gradual decrease from April to August, to a minimum value of 24 °C.Finally, the SST increases until December to reach the maximum value of 28.5 °C.We observe that cold water is present over the entire Ivorian coastal zone during the major cold season, with an average temperature of 24 °C.During this season, the coldest surface temperatures are observed in the central and eastern parts of the littoral between 5 and 2°W during the month of August.Warmer temperatures are found in the western part between longitudes 7°W and 8°W during the month of July.
The seasonal variation in coastal upwelling in the fishing zone of Abidjan is shown in Figure 2c.The upwelling index increases from 500 m 3 s −1 km −1 in January to 1100 m 3 s −1 km −1 in March, then decreases slightly in May to 750 m 3 s −1 km −1 .From June, there is a rapid and maximal increase in upwelling index, reaching values of 2000 m 3 s −1 km −1 in July-August.Finally, values decrease gradually to reach 500 m 3 s −1 km −1 in December.In the fishing zone of Abidjan, wind-induced turbulent mixing (Figure 2d) follows the same seasonal variations as the upwelling index.Turbulent mixing increases from 10 m 3 s −3 in January, reaching 35 m 3 s −3 in March-April-May.Then, turbulent mixing values increase again, reaching a maximum of 80 m 3 s −3 in July-August, before gradually decreasing to values of 10 m 3 s −3 in December.The variations in coastal upwelling index and turbulent mixing are consistent with in situ measurements of thermocline depth, showing shallow depths between 0 and 50 m (upper and lower bounds, respectively) from July to September during the major cold season, and between 10 and 60 m from January to March during the minor cold season.The thermocline deepens to between 15 and 65 m (upper and lower bounds, respectively), on average, between April and June, and between October and December.
The seasonal distribution of depth-integrated nutrient concentration showed a rapid increase in nitrate from January to March, reaching maximal concentration in April-May, with values of 0.55 mol m −2 on average.Nitrate concentration decreases to a minimum of 0.32 mol m −2 in October, which is then followed by increasing concentration from November onwards.Seasonal distribution in phosphate showed lower concentrations from March to May, and then from July to August, with mean values of 0.033 and 0.025 mol m −2 , respectively.Phosphate concentration increases, again, from December to February, reaching values of 0.041 mol m −2 , on average.Water column oceanographic profiles are from the World Ocean Atlas 2013 [56].
The seasonal distribution of depth-integrated nutrient concentration showed a rapid increase in nitrate from January to March, reaching maximal concentration in April-May, with values of 0.55 mol m −2 on average.Nitrate concentration decreases to a minimum of 0.32 mol m −2 in October, which is then followed by increasing concentration from November onwards.Seasonal distribution in phosphate showed lower concentrations from March to May, and then from July to August, with mean values of 0.033 and 0.025 mol m −2 , respectively.Phosphate concentration increases, again, from December to February, reaching values of 0.041 mol m −2 , on average.


Temporal variations in S. aurita catch, chlorophyll, and timing of initiation In the fishing zone of Abidjan, timing of initiation of phytoplankton growth shows a synchronous relationship with S. aurita catch, i.e., when S. aurita catch is high, the timing of initiation is delayed, during the period from 1998 to 2014 (Figure 6a).Whereas analysis of S. aurita catch and concentration of chlorophyll shows an inverse relationship between the two variables (Figure 6b), i.e., when S. aurita catch is high, there is a large decrease in the concentration of chlorophyll and viceversa.The latter relationship is observed in most of the years during the period from 1998 to 2014, except for the years 1999-2000 and from 2012 to 2014, when both variables change synchronously.• Temporal variations in S. aurita catch, chlorophyll, and timing of initiation In the fishing zone of Abidjan, timing of initiation of phytoplankton growth shows a synchronous relationship with S. aurita catch, i.e., when S. aurita catch is high, the timing of initiation is delayed, during the period from 1998 to 2014 (Figure 6a).Whereas analysis of S. aurita catch and concentration of chlorophyll shows an inverse relationship between the two variables (Figure 6b), i.e., when S. aurita catch is high, there is a large decrease in the concentration of chlorophyll and vice-versa.The latter relationship is observed in most of the years during the period from 1998 to 2014, except for the years 1999-2000 and from 2012 to 2014, when both variables change synchronously.


Temporal variations in S. aurita catch, chlorophyll, and timing of initiation In the fishing zone of Abidjan, timing of initiation of phytoplankton growth shows a synchronous relationship with S. aurita catch, i.e., when S. aurita catch is high, the timing of initiation is delayed, during the period from 1998 to 2014 (Figure 6a).Whereas analysis of S. aurita catch and concentration of chlorophyll shows an inverse relationship between the two variables (Figure 6b), i.e., when S. aurita catch is high, there is a large decrease in the concentration of chlorophyll and viceversa.The latter relationship is observed in most of the years during the period from 1998 to 2014, except for the years 1999-2000 and from 2012 to 2014, when both variables change synchronously. ) in year t.Note that in (a), the relationship is shown for phytoplankton bloom timing in year t and S. aurita catch in year t + 1, as we are investigating the influence of the bloom timing (food availability) on the S. aurita recruitment success, which is shown in the catch success in following year.On the other hand, in (b), the relationship is shown for chlorophyll concentration and S. aurita catch both in the same year, t, as we are investigating grazing ) in year t.Note that in (a), the relationship is shown for phytoplankton bloom timing in year t and S. aurita catch in year t + 1, as we are investigating the influence of the bloom timing (food availability) on the S. aurita recruitment success, which is shown in the catch success in following year.On the other hand, in (b), the relationship is shown for chlorophyll concentration and S. aurita catch both in the same year, t, as we are investigating grazing pressure of S. aurita on chlorophyll concentration in the same year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure 2).• Temporal variations in S. aurita catch, turbulent mixing, and upwelling index In the fishing zone of Abidjan, S. aurita catch is positively related to wind-induced turbulent mixing and upwelling index (Figure 9a,b).Decreasing trends are shown in these physical variables from 2000 to 2011, associated with a decreasing trend in the quantity of S. aurita caught.Furthermore, in the first years of the time series 1999-2000, the upwelling and turbulent mixing are at their highest level, and sardine catch is high, whereas, in the years 2010 and 2011, upwelling and turbulent mixing are at their lowest, and low catches of sardines are observed.

Empirical Relationships between Sardinella aurita and Biophysical Variables
• S. aurita catch and first trophic level Our analysis in the fishing zone of Abidjan shows significant linear relationships between S. aurita catch and variables of chlorophyll concentration and timing of initiation (Figure 10a,b).Catches of S. aurita are inversely related to the concentrations of chlorophyll, such that low concentration of chlorophyll in year t is associated with high S. aurita catch in same year.Interannual variations in the concentration of chlorophyll can explain 48% of the variations in S. aurita catch.This negative relationship may indicate an important grazing pressure from Sardinella larvae on phytoplankton.A positive linear relationship is observed between the timing of initiation of phytoplankton growth in year t and S. aurita catch in year t + 1.Early (late) timing of initiation of phytoplankton growth is associated with low (high) S. aurita catch in the following year.If the timing of phytoplankton bloom occurs earlier (negative anomaly) than average conditions, a mismatch may be observed between food availability (phytoplankton) and the arrival of Sardinella larvae.Under such conditions, survival of Sardinella larvae may be low, which may then translate into low adult fish recruitment and consequent low catch in year t + 1.By contrast, delay in the timing of phytoplankton blooms (positive anomalies) may bring more favorable conditions for the survival of Sardinella larvae (good timing of food availability), which may then allow for high adult fish recruitment and their catch in year t + 1.


Temporal variations in S. aurita catch, turbulent mixing, and upwelling index In the fishing zone of Abidjan, S. aurita catch is positively related to wind-induced turbulent mixing and upwelling index (Figure 9a,b).Decreasing trends are shown in these physical variables from 2000 to 2011, associated with a decreasing trend in the quantity of S. aurita caught.Furthermore, in the first years of the time series 1999-2000, the upwelling and turbulent mixing are at their highest level, and sardine catch is high, whereas, in the years 2010 and 2011, upwelling and turbulent mixing are at their lowest, and low catches of sardines are observed.Note that in (a), the relationship is shown for upwelling index and wind-induced turbulent mixing in year t, and S. aurita catch in year t + 1, as we are investigating the influence of the ocean physics on the S. aurita recruitment success, which is shown in the catch success in following year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure 2).Note that in (a), the relationship is shown for upwelling index and wind-induced turbulent mixing in year t, and S. aurita catch in year t + 1, as we are investigating the influence of the ocean physics on the S. aurita recruitment success, which is shown in the catch success in following year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure 2).

• S. aurita catch and ocean physics
Our analysis in the fishing zone of Abidjan shows significant positive linear relationships between S. aurita catch, and the dynamics of the physical variables of upwelling index and wind-induced turbulent mixing (Figure 10c,d).For moderate turbulent mixing (from June to December) and favorable upwelling (annual mean), there is generally good recruitment of S. aurita on the Ivorian continental shelf.High upwelling index may increase nutrients and food availability to the surface waters, bringing favorable feeding conditions for the Sardinella, and leading to higher catches.In addition, the turbulent mixing of the months from June to December presents favorable conditions for a strong increase in nutrients in the coastal regions.In this context, high wind-induced turbulent mixing is shown to be positively related to S. aurita catch, with turbulent mixing explaining 69% of the variations in S. aurita catch from 1998 to 2014.Note that in (a), the relationship is shown for chlorophyll concentration and S. aurita catch both in the same year t, as we are investigating grazing pressure on chlorophyll concentration in that year.Similarly, in (c-f) the relationships are shown for variables in the year t, as we are investigating the influence of ocean physics on the concentration and timing of the first trophic level in the same year, and in (i,j), as we are investigating autocorrelation between physical variables and between biological variables.However, in (b,g,h), the relationship is shown for biophysical variables in year t and S. aurita catch in year t + 1, as we are investigating the influence of the biophysical variables on the S. aurita catch success in following year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure 2).Note that in (a), the relationship is shown for chlorophyll concentration and S. aurita catch both in the same year t, as we are investigating grazing pressure on chlorophyll concentration in that year.Similarly, in (c-f) the relationships are shown for variables in the year t, as we are investigating the influence of ocean physics on the concentration and timing of the first trophic level in the same year, and in (i,j), as we are investigating autocorrelation between physical variables and between biological variables.However, in (b,g,h), the relationship is shown for biophysical variables in year t and S. aurita catch in year t + 1, as we are investigating the influence of the biophysical variables on the S. aurita catch success in following year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure 2).

• First trophic level and ocean physics
Our analysis in the fishing zone of Abidjan shows significant positive linear relationships between phytoplankton and ocean physics (Figure 10e-h).On the one hand, the high upwelling index and wind-induced turbulent mixing observed in the present analysis allow sufficient nutrient supply for phytoplankton to grow and develop in a bloom, such that bottom-up trophic interactions may not be a limiting factor.On the other hand, favorable environmental conditions may lead to enhanced grazing pressure by Sardinella larvae, exerting important top-down control on the phytoplankton (low chlorophyll concentration), such that increased turbulent mixing and upwelling index may be associated with good recruitment of S. aurita larvae, increasing grazing pressure on phytoplankton, leading to lower chlorophyll concentration, and a delayed bloom initiation, i.e., time when chlorophyll reaches the threshold of median plus 20% (Figure 10b).

Linear Regression Analysis and Model Equation
Based on the empirical relationships demonstrated between S. aurita catch and physical-biological variables (Section 3.3.1),diagnostic models for prediction of S. aurita catch in the fishing zone of Abidjan can now be explored.To identify the best combination of predictors, we performed an exhaustive search over all possible combinations of input variables derived from two biological and two physical variables (i.e., there are 15 possible combinations out of four variables).The bilinear regression model using only the biological variables (timing of initiation of phytoplankton growth and chlorophyll concentration) can explain 40% of the variations in S. aurita catch, whereas the model using only the physical variables (upwelling index and turbulent mixing) shows that in year t, these two variables may allow us, with a confidence of 78%, to predict S. aurita catch in year t + 1 (Figure 11e,i).The latter score is also observed in the multilinear regression model using chlorophyll concentration, upwelling index and turbulent mixing as predictors (Figure 11m), as well as in the multilinear regression model using timing of initiation of phytoplankton growth, chlorophyll concentration, upwelling index, and turbulent mixing as predictors (Figure 11o).However, statistical analysis of the coefficient estimates indicates that only the variables turbulent mixing and upwelling index significantly improve the fit of these models.
Based on these results, we provide the predictive equation using only the physical variables, as follows: However, although the physical variables may be statistically better predictors of the variations in S. aurita catch in year t + 1, they have no direct influence on the trophic interactions of this species.To assess how much the variations in S. aurita catch in year t + 1 can be explained by changes at the base of the marine food chain (i.e., phytoplankton food availability), we need to consider predictive models using biological variables as predictors (Figure 11a,b,e).We found that phytoplankton chlorophyll concentration and bloom timing initiation alone can explain 44% and 31% of S. aurita catch in year t + 1, respectively (Equations ( 2) and ( 3)), and that the combination of these two variables can explain 40% of the variations in S. aurita catch (Equation ( 4)).The equations describing these linear models are given as follows:

Seasonal and Interannual Variability in Biophysical Conditions
The SST along the Ivorian coastline is characterized by large seasonal variations [61,62].Based on these thermal variations, four marine seasons have been defined: the minor cold season from January to February; the major warm season from March to June, with a maximum temperature centered on the month of May; the major cold season from July to September; and finally, the minor warm season from November to December.Our results showed that the peaks of chlorophyll concentration, upwelling index, and wind-induced turbulent mixing generally coincide with low SST.In support of our results, similar observations have been previously reported in several studies based on in situ measurements [61][62][63][64][65] and a hydrodynamic regional model [37], describing the seasonal development of the upwelling (the resurgence of cold water) off the Côte d'Ivoire.Highest intensity of upwelling is observed from July to September.During this period, upwelling enriches the ocean with nutrients, with immediate influence on primary production [25,31,66].Nutrient enrichment has been shown to control the development of phytoplankton, such that chlorophyll and nitrate concentrations may be linearly related during the months of July and August [67], and the phenomenon of upwelling can be localized by the presence cold seawater that is rich in phytoplankton.In situ measurements of nitracline depth-integrated nitrate concentration and phosphatocline depth-integrated phosphate concentration (Figure 3c,e respectively) indicate that phytoplankton consume nutrients (decreasing concentration in July to August).However, nutrient depletion is not observed, suggesting that nutrient availability may not be limiting phytoplankton growth during the major cold season, and other factors may play an important role, including grazing and light availability.Influence of grazing is discussed in detail in Section 4.2.The role played by the light availability in controlling phytoplankton production in Ivorian waters has been previously assessed by Herbland and Voituriez [67], who demonstrated that the depth of the thermocline may be inversely related to the slope of the chlorophyll-a concentration and in vivo fluorescence.In the region, previous studies [68][69][70] have also reported that when the wind is at the origin of the upwelling, it causes mixing that can be associated with reduced light availability within the mixed layer, which can become one of the factors limiting production.Further increase in wind speed may also increase the speed of surface currents, which can drive phytoplankton-rich waters towards the open ocean.Such changes in wind conditions may change availability of food (phytoplankton) for fish larvae and affect fish recruitment (see Section 4.2).
Although upwelling is the main source of enrichment in the euphotic zone of the littoral, the contribution of the discharges of the main rivers should not be ignored [30].According to Binet [23], and Herbland and Le Loeuff [25], during the great warm season, terrigenous contributions may be provided to the Ivorian coastal zone.These nutrient-rich inputs come from runoff from the land surrounding the coastal lagoons, and peak in October [71].They are especially important in the eastern part of the coast, and may increase the concentration of phytoplankton in the area [29].

Influence of Biophysical Conditions on the Recruitment of Sardinella Aurita
Coastal upwelling ecosystems are known to be highly productive [18,72].Generally, fish avoid breeding in areas dominated by strong currents and water column mixing.The spawning sites are usually located downstream from the main upwelling centers, and breeding is mainly tuned to avoid the peak of the upwelling season [73].On the Ivorian continental shelf, we show that SST is low at the coast, and increases towards the open sea, whereas chlorophyll concentration is high at the coast and decreases towards the open sea.In this area, phytoplankton and other biological components retained at the coast provide favorable conditions for the reproduction of S. aurita.Our results are in good agreement with those of Roy [73], who reported, during a study of upwelling-induced retention area off Senegal, that the coastal zone represents a very favorable environment for fish reproduction.The author indicated that eggs and larvae might develop in productive and relatively stable coastal waters.In the Ivorian continental shelf, upwelling and coastal retention of larvae occur simultaneously and, together, provide favorable conditions for fish reproduction.
In the present study, we have shown a significant linear relationship between upwelling index and wind-induced turbulent mixing strength, suggesting the presence of an Ekman-type upwelling system, according to the definition by Cury and Roy [74].The latter authors report the existence of an "optimal environmental window" for the recruitment of sardines, which is characterized by moderate values of turbulent mixing and upwelling index.Along the Ivorian continental shelf, during the period 1998-2014, turbulent mixing and upwelling index remained within the range of moderate values, 800 to 1400 m 3 s −3 , and 30 to 60 m 3 s −1 km −1 , respectively (Figure 10), which are commensurate with values identified by Cury and Roy [74], suggesting that throughout the period of study, the ocean physical conditions were "optimal" for fish recruitment.In addition, during the major upwelling, nutrient concentration is shown to be sufficient (Figure 3d,f; [71]), supporting the growth of phytoplankton.In these conditions, bottom-up control may not be the limiting factor for the recruitment of Sardinella larvae, and the observed year-to-year variations in S. aurita catch may be driven by top-down control (Figure 12a,b).
Remote Sens. 2017, 9, x FOR PEER REVIEW 17 of 27 waters.In the Ivorian continental shelf, upwelling and coastal retention of larvae occur simultaneously and, together, provide favorable conditions for fish reproduction.
In the present study, we have shown a significant linear relationship between upwelling index and wind-induced turbulent mixing strength, suggesting the presence of an Ekman-type upwelling system, according to the definition by Cury and Roy [74].The latter authors report the existence of an "optimal environmental window" for the recruitment of sardines, which is characterized by moderate values of turbulent mixing and upwelling index.Along the Ivorian continental shelf, during the period 1998-2014, turbulent mixing and upwelling index remained within the range of moderate values, 800 to 1400 m 3 s −3 , and 30 to 60 m 3 s −1 km −1 , respectively (Figure 10), which are commensurate with values identified by Cury and Roy [74], suggesting that throughout the period of study, the ocean physical conditions were "optimal" for fish recruitment.In addition, during the major upwelling, nutrient concentration is shown to be sufficient (Figure 3d,f; [71]), supporting the growth of phytoplankton.In these conditions, bottom-up control may not be the limiting factor for the recruitment of Sardinella larvae, and the observed year-to-year variations in S. aurita catch may be driven by top-down control (Figure 12a,b).In the Ivorian shelf waters, S. aurita spawns during or after the cold season (June to September), and the highest abundance of larvae is observed in August [75].The diet of Sardinella larvae may include small crustaceans (copepods), but, along the coasts of Senegal, Congo [76,77], Côte d'Ivoire [78], and Ghana [79], their main source of food has been shown to be phytoplankton.In the Congolese fishery of S. aurita, Ghéno and Fontana [80] have observed that the variability in spawning in the previous seasons is reflected in a strong variability of the catch rates of the fishery.Furthermore, enhanced S. aurita catch in the Ivorian fishery has been observed when the onset of the phytoplankton bloom occurs early in April-May, just preceding the main spawning and recruitment period of S. aurita between June and October (with a maximum in July-August; [81,82]).During this period, the vulnerability of fish stocks to fishing activities is particularly high [36,[83][84][85].These reports are consistent with our findings of a positive relationship between the phytoplankton bloom timing and the S. aurita catch in the following year, which may be further explained as follows (Figure 10): in the years when the phytoplankton flowering begins early (e.g., negative anomalies of ~2-4 weeks; Figure 10a), seasonal chlorophyll concentration (i.e., food availability) is high when Sardinella larvae are released, and the larvae have sufficient food in their critical period to eat on and grow.Such conditions are favorable for successful recruitment of Sardinella larvae, and lead to a higher Sardinella catch in the year t + 1.In the years when phytoplankton flowering is delayed (e.g., positive anomalies of ~2-4 weeks; Figure 10a), the amount of food available (i.e., chlorophyll concentration) to the larvae is low, and larvae may die by starvation.In such conditions, Sardinella larval recruitment is low, and the Sardinella catch in the following year may remain low.In the Ivorian shelf waters, S. aurita spawns during or after the cold season (June to September), and the highest abundance of larvae is observed in August [75].The diet of Sardinella larvae may include small crustaceans (copepods), but, along the coasts of Senegal, Congo [76,77], Côte d'Ivoire [78], and Ghana [79], their main source of food has been shown to be phytoplankton.In the Congolese fishery of S. aurita, Ghéno and Fontana [80] have observed that the variability in spawning in the previous seasons is reflected in a strong variability of the catch rates of the fishery.Furthermore, enhanced S. aurita catch in the Ivorian fishery has been observed when the onset of the phytoplankton bloom occurs early in April-May, just preceding the main spawning and recruitment period of S. aurita between June and October (with a maximum in July-August; [81,82]).During this period, the vulnerability of fish stocks to fishing activities is particularly high [36,[83][84][85].These reports are consistent with our findings of a positive relationship between the phytoplankton bloom timing and the S. aurita catch in the following year, which may be further explained as follows (Figure 10): in the years when the phytoplankton flowering begins early (e.g., negative anomalies of ~2-4 weeks; Figure 10a), seasonal chlorophyll concentration (i.e., food availability) is high when Sardinella larvae are released, and the larvae have sufficient food in their critical period to eat on and grow.Such conditions are favorable for successful recruitment of Sardinella larvae, and lead to a higher Sardinella catch in the year t + 1.In the years when phytoplankton flowering is delayed (e.g., positive anomalies of ~2-4 weeks; Figure 10a), the amount of food available (i.e., chlorophyll concentration) to the larvae is low, and larvae may die by starvation.In such conditions, Sardinella larval recruitment is low, and the Sardinella catch in the following year may remain low.
During the period of study 1998-2014, year-to-year variations in Sardinella larvae are not available to allow direct validation of the match-mismatch hypothesis [86] between Sardinella larvae and phytoplankton bloom timing.However, we also observed negative relationships between phytoplankton chlorophyll concentration in year t and Sardinella catch in the same year t, suggesting that Sardinella larvae exert a significant grazing pressure on phytoplankton.Hence, to provide an indirect validation of the hypothesis about the relationship between Sardinella recruitment and phytoplankton availability, we next estimate, empirically, the carbon consumption rate by Sardinella larvae in the area based on the relationship demonstrated in Section 3.3.1;and we then compare this estimate with in situ consumption rates of Sardinella larvae reported in the literature.
• Estimation of carbon consumption rate in the study area Figure 10a shows, for different years, the total sardine catch (a proxy for abundance) plotted against the annual average of chlorophyll concentration in the same years.In the years when the abundance of sardines is high, the chlorophyll is low (and vice versa), implying a role for grazing on chlorophyll by sardines.From the slope of the regression line (i.e., 22 × 10 3 tons sardines per mg Chl m −3 ), we could infer that an increase in sardine population of 22 × 10 3 tons would be associated with a decrease in chlorophyll of 1 mg m −3 , averaged over the study area.We interpret this decrease as the average food supply required to sustain the increase in sardine abundance, and make the assumption that the 1 mg m −3 of chlorophyll represents the average daily food consumption of 22 × 10 3 tons.Finally, we proceed to estimate the carbon consumption rate in the study area of the Abidjan fishing zone, which is characterized by a surface of 2 × 10 12 m 2 , and mixed layer depth of 60 m (the latter value is based on vertical profiles of chlorophyll concentration measured in front of Abidjan during the sampling campaign reported in Herbland and Voituriez [67].By multiplying the consumption rate of 1 mg Chl m −3 day −1 by the depth of the mixed layer and surface of the study area, we calculate the total amount of chlorophyll consumed to be 4.4 × 10 16 mg Chl year −1 .Assuming a carbon (C) to Chl ratio of 50 [87], we obtain a monthly consumption rate of carbon of 1.9 × 10 17 mg C month −1 or 2.5 × 10 17 µg C h −1 .

• Estimation of larval abundance and consumption rate
Assuming that adult Sardinella weighs, on average, 229 g [88], the number of Sardinella in the water volume of the study area would be 5 × 10 12 individuals or 0.04 individuals per m 3 .We hypothesize that half of the individual Sardinella are females, and each female produces an average of 80,000 larvae [89].Then, we obtain the total number of larvae, 2.1 × 10 17 , in the water volume of the study area, or about 1700 larvae per m 3 .Finally, the carbon consumption rate may be divided by the total number of larvae, which gives 1.2 µg C larva −1 h −1 .The latter result is within the same order of magnitude compared to previously published figures by Caldeira et al. [90], who reported carbon intake rates for sardine species Sardina pilchardus in European waters between 4.5 and 7 µg C larva −1 h −1 .The relatively lower carbon consumption rate observed here for Sardinella aurita may be due to the fact that (1) the European Sardina pilchardus is likely to present anatomical and physiological differences, as well as living in different environmental conditions, compared with the species Sardinella aurita studied here in Ivorian waters; and (2) our calculation is based on total Sardinella catch, and this value is likely to underestimate the total number of Sardinella individuals present in the study area (i.e., not all Sardinella may have been caught).

Relevance of Remote-Sensing Observations and Diagnostic Models for Sardinella Fisheries Management and Operational Applications
Since the intensification of fishing activity and the increase in awareness of the need to monitor fish stocks closely, evaluation methods and reporting of fish catches have been progressively standardized at national and international levels.To date, 441 of the 584 commercial stocks of marine animals throughout the world are routinely evaluated [91].
With a continental shelf of 12,000 km 2 , the Côte d'Ivoire is among the poorest countries in fishery resources in the area of the Gulf of Guinea.Small pelagics, whose abundance depends on environmental conditions, such as the seasonal appearance of coastal cold water, are estimated to be 81,000 tons between Togo, Ghana, and Côte d'Ivoire [92].In recent years, several analyses have reported the collapse and scarcity of a large number of commercially important fish stocks across these countries [93][94][95].The continuous increase in fishing pressure beyond self-renewal capability of already limited fish stocks, is identified as the main cause of declining abundance of the fishery resource, and represents the greatest threat to the sustainable development of fishing activity in the region [96,97].
Sardinella aurita plays a significant role in the profitability of the fishing industry, and is mainly responsible for supplying local consumers and neighboring regions with protein-rich, high-quality food at a low price.However, the structures responsible for the management of the Ivorian fishery face several difficulties.There is a lack of information in the sector as a whole, with weak decision-making tools, weak policy and regulatory frameworks, and poverty in fishing communities.To tackle these issues, we propose that the biophysical metrics of phytoplankton phenology, phytoplankton biomass, turbulence, and upwelling index, identified in the present study, be incorporated into stock assessments within an ecosystem-based management framework for Côte d'Ivoire.The metrics may be readily, objectively, and freely estimated from remote-sensing and reanalysis products, and they meet operational requirements in terms of spatiotemporal resolution, repeat frequency, and cost-effectiveness [12].Furthermore, we have demonstrated that, applied in a multilinear diagnostic model, these biophysical metrics may allow us, with a confidence of 78%, to predict the S. aurita catch in the following year.
In this context, we suggest the following approach towards an implementation of a Fisheries Information and Management System for Ivorian waters: (1) On a monthly basis, retrieval of data metrics, mapping of metric anomalies, and plotting of time-series metrics; (2) On a yearly basis, estimation of annual mean S. aurita catch, annual mean upwelling index, June to December mean wind-induced turbulent mixing, annual mean chlorophyll concentration, calculation of long-term chlorophyll threshold criterion and timing of initiation of the phytoplankton growing period; (3) Analysis of empirical relationships between S. aurita catch and the four biophysical variables, using data for all available years from 1998 to present; (4) If the empirical relationships are significant, estimate the new parameters of the multilinear diagnostic model and forecast, with confidence levels, the S. aurita catch in the following year; If the empirical relationships are not significant, further analyses and search of optimal diagnostic model need to be performed in consultation with experts from CURAT or other institutes, as relevant; (5) Produce an information bulletin providing an annual summary of biophysical metric status, including anomaly maps (such as presented for instance in Figure 4, Figure 5, Figure 7, and Figure 8), time-series plots (for instance, as in Figures 6 and 9), and the estimate, with a confidence level, of S. aurita catch forecast for the following year; (6) Review the information bulletin with the representatives in charge of the S. aurita stock assessment and issue advice on protection measures, as required; (7) Provide evidenced-based management advice to government agencies, enabling them to make informed decisions.

Conclusions
In this study, a regression model based on physical variables has been shown to explain 78% of the catch of Sardinella aurita in the following year (year t + 1) when the annual mean upwelling index and mean wind-induced turbulent mixing in the months of Jun to Dec of the current year (year t) are known.On the other hand, the regression model based on biological variables of the timing of phytoplankton bloom initiation and chlorophyll concentration has been shown to explain 40% of the catch of Sardinella aurita in the following year.Although the physical variables may, statistically, be better predictors of the variations in Sardinella catch in the following year, we have highlighted that, during the period 1998-2014, the upwelling and turbulence strengths have remained of moderate intensity, suggesting non-limiting nutrient supply to primary producers and trophic interactions, primarily driven by top-down control.In this context, analysis and incorporation of biological variables (into predictive models) is essential to identify alterations in trophic interactions, and to understand the fluctuations in commercially-important fisheries, such as Sardinella.Moreover, even though the predictive models developed here are based on nearly two decades of data, their parameters may be influenced by interdecadal variability and long-term regime shifts [98][99][100][101].Examining the relationships between interdecadal climate oscillations and biological variables, such as phytoplankton phenology, as well as between other physical variables, such as precipitation and river discharge, and biological variables, may provide further insights toward improving projection of Sardinella catch in the Ivorian waters.This would bring further strength and robustness to the predictive tools, and to the proposed approach towards implementation of a Fisheries Information and Management System for Ivorian waters, which would help to lay the foundation for long-term economic growth and sustainability for the benefit of this species, and for society itself.

Figure 1 .
Figure 1.Comparison of Sardinella aurita landings (tons) in the fishing zone of Abidjan produced by the Ministry of Production of Animals (MPA) at the Abidjan Fisheries Direction with the corresponding landings for all Ivorian waters produced by the Food and Agriculture Organization (FAO) for the period 1997-2014.

Figure 1 .
Figure 1.Comparison of Sardinella aurita landings (tons) in the fishing zone of Abidjan produced by the Ministry of Production of Animals (MPA) at the Abidjan Fisheries Direction with the corresponding landings for all Ivorian waters produced by the Food and Agriculture Organization (FAO) for the period 1997-2014.

Figure 3 .
Figure 3. Seasonal variations of (a) thermocline depth (m), (b) surface temperature (°C), (c) nitracline depth (m), (d) nitracline depth-integrated nitrate concentration (mol m −2 ), (e) phosphatocline depth (m), and (f) phosphatocline depth-integrated phosphate concentration (mol m −2 ).Left panel: blue line is upper bound of the layer and turquoise line is lower bound of the layer.All observations are based on in situ measurements collected in Abidjan fishing zone at the point of coordinate 5.5°W, 4.5°N.Water column oceanographic profiles are from the World Ocean Atlas 2013[56].

Figure 3 .
Figure 3. Seasonal variations of (a) thermocline depth (m), (b) surface temperature ( • C), (c) nitracline depth (m), (d) nitracline depth-integrated nitrate concentration (mol m −2 ), (e) phosphatocline depth (m), and (f) phosphatocline depth-integrated phosphate concentration (mol m −2 ).Left panel: blue line is upper bound of the layer and turquoise line is lower bound of the layer.All observations are based on in situ measurements collected in Abidjan fishing zone at the point of coordinate 5.5 • W, 4.5 • N.Water column oceanographic profiles are from the World Ocean Atlas 2013[56].

3. 2 . 1 .
Phytoplankton Indices • Spatial distribution in chlorophyll concentration anomalies The year-to-year variations of chlorophyll concentration over the period 1998 to 2014 are represented in Appendix A Figure A1 as annual means, expressed in mg Chl m −3 , and in Figure 4 as the relative anomalies, expressed as percent change from the 1998-2014 mean.In the fishing zone of Abidjan (box area in Figure 4), negative anomalies are prominent from 1998 to 2002, then in 2007, and finally, from 2012 to 2013.The lowest relative chlorophyll anomalies (−40%) are observed in 1998, 1999 and 2001.In 2002, 2007, 2012, and 2013, the negative anomalies are in the range of −20% to −5%.In 2006, 2010, 2011, and 2014, the relative chlorophyll concentration anomalies are close to zero on average, with values lying between −5 and +5%.The positive relative chlorophyll anomalies are dominant in 2003-2005, 2008, and 2009, with a maximum percentage value of +35% in 2008.Remote Sens. 2017, 9, x FOR PEER REVIEW 8 of 27 3.2.Interannual Variability of Biophysical Drivers in Ivorian Waters 3.2.1.Phytoplankton Indices  Spatial distribution in chlorophyll concentration anomalies The year-to-year variations of chlorophyll concentration over the period 1998 to 2014 are represented in Appendix Figure A1 as annual means, expressed in mg Chl m −3 , and in Figure 4 as the relative anomalies, expressed as percent change from the 1998-2014 mean.In the fishing zone of Abidjan (box area in Figure 4), negative anomalies are prominent from 1998 to 2002, then in 2007, and finally, from 2012 to 2013.The lowest relative chlorophyll anomalies (−40%) are observed in 1998, 1999 and 2001.In 2002, 2007, 2012, and 2013, the negative anomalies are in the range of −20% to −5%.In 2006, 2010, 2011, and 2014, the relative chlorophyll concentration anomalies are close to zero on average, with values lying between −5 and +5%.The positive relative chlorophyll anomalies are dominant in 2003-2005, 2008, and 2009, with a maximum percentage value of +35% in 2008.

Figure 4 .
Figure 4. Maps of relative anomalies in annual mean chlorophyll concentration expressed in percentage in Ivorian waters.Chlorophyll concentration anomalies are calculated using the Ocean Color Climate Change Initiative (OC-CCI) product for the period 1998 to 2014.Blue (red) color indicates lower (higher) chlorophyll concentration relative to the 17 year mean.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan. Spatial distribution in anomalies of phytoplankton phenology Year-to-year variations in timing of initiation of phytoplankton growth over the period 1998 to 2014 are shown in Appendix Figure A2, and the anomalies in Figure 5. Three types of phytoplankton growth initiation times were observed from 1998 to 2014, namely early initiation (blue color in Figure 5), late initiation (indicated in red), and normal initiation (white).In the fishing zone of Abidjan, earlier than average timing of initiation between −4 and −2 weeks is observed in the years 2003-2005, 2010, 2013, and 2014, whereas in 2008, the phytoplankton growth started particularly early (−6 weeks).Moderate delay in the initiation of phytoplankton growth between +1 and +4 weeks is observed in 2007, 2009, and 2012.Longer delay, with greater than +4 weeks, are estimated for the first

Figure 4 . 27 Figure 5 .
Figure 4. Maps of relative anomalies in annual mean chlorophyll concentration expressed in percentage in Ivorian waters.Chlorophyll concentration anomalies are calculated using the Ocean Color Climate Change Initiative (OC-CCI) product for the period 1998 to 2014.Blue (red) color indicates lower (higher) chlorophyll concentration relative to the 17 year mean.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure 6 .
Figure 6.Interannual variability during the period 1998 to 2014 in (a) Sardinella aurita catch (tons) in year t + 1, and timing of initiation (weeks) in year t; and (b) Sardinella aurita catch (tons) in year t and chlorophyll concentration (mg m −3 ) in year t.Note that in (a), the relationship is shown for

Figure 5 .
Figure 5. Maps of anomalies of timing of initiation (8-Day period) of the phytoplankton growth in Ivorian waters.The timing of initiation is calculated using the OC-CCI chlorophyll product and applying the threshold method presented in Racault et al. [57].Anomalies are estimated for the period 1998 to 2014.Blue (red) color indicates earlier (later) phytoplankton growth initiation compared with the 17 year mean.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

3. 2
.2. Physical Indices • Spatial distribution in anomalies of upwelling index Upwelling index is an important variable in the determination of fishing areas on the coasts.The distribution of the upwelling index for the period from 1998 to 2014 in Ivorian waters is presented in Appendix A Figure A3, and the anomalies in Figure 7.In 2002, 2003, 2005, and 2007, upwelling index is close to average condition, and anomalies oscillate between −100 and 100 m 3 s −1 km −1 .In the area of the fishing zone of Abidjan, negative anomalies are observed in 1999, 2004, 2006, and from 2008 to 2012, with lowest values of −500 m 3 s −1 km −1 , while positive upwelling index anomalies are observed in the years 1998, 2000, 2001, 2013, and 2014 with maximum values of +500 m 3 s −1 km −1 .

Figure 6 .
Figure 6.Interannual variability during the period 1998 to 2014 in (a) Sardinella aurita catch (tons) in year t + 1, and timing of initiation (weeks) in year t; and (b) Sardinella aurita catch (tons) in year t and chlorophyll concentration (mg m −3) in year t.Note that in (a), the relationship is shown for phytoplankton bloom timing in year t and S. aurita catch in year t + 1, as we are investigating the influence of the bloom timing (food availability) on the S. aurita recruitment success, which is shown in the catch success in following year.On the other hand, in (b), the relationship is shown for chlorophyll concentration and S. aurita catch both in the same year, t, as we are investigating grazing

Figure 6 .
Figure 6.Interannual variability during the period 1998 to 2014 in (a) Sardinella aurita catch (tons) in year t + 1, and timing of initiation (weeks) in year t; and (b) Sardinella aurita catch (tons) in year t and chlorophyll concentration (mg m −3) in year t.Note that in (a), the relationship is shown for phytoplankton bloom timing in year t and S. aurita catch in year t + 1, as we are investigating the influence of the bloom timing (food availability) on the S. aurita recruitment success, which is shown in the catch success in following year.On the other hand, in (b), the relationship is shown for chlorophyll concentration and S. aurita catch both in the same year, t, as we are investigating grazing pressure of S. aurita on chlorophyll concentration in the same year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure2).
Remote Sens. 2017, 9, x FOR PEER REVIEW 10 of 27 pressure of S. aurita on chlorophyll concentration in the same year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure2).

3. 2 . 2 .
Physical Indices  Spatial distribution in anomalies of upwelling index Upwelling index is an important variable in the determination of fishing areas on the coasts.The distribution of the upwelling index for the period from 1998 to 2014 in Ivorian waters is presented in Appendix Figure A3, and the anomalies in Figure 7.In 2002, 2003, 2005, and 2007, upwelling index is close to average condition, and anomalies oscillate between −100 and 100 m 3 s −1 km −1 .In the area of the fishing zone of Abidjan, negative anomalies are observed in 1999, 2004, 2006, and from 2008 to 2012, with lowest values of −500 m 3 s −1 km −1 , while positive upwelling index anomalies are observed in the years 1998, 2000, 2001, 2013, and 2014 with maximum values of +500 m 3 s −1 km −1 .

Figure 7 .
Figure 7. Maps of upwelling index anomalies (m 3 s −1 km −1 ).Anomalies are estimated for the period 1998 to 2014.Blue (red) color indicates stronger (weaker) upwelling conditions compared to the 17 year mean.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.Spatialdistribution in wind-induced turbulent mixing anomaliesThe year-to-year variations in annual mean wind-induced turbulent mixing for the period 1998 to 2014 are shown in Appendix FigureA4.Wind-induced turbulent mixing anomalies are presented in Figure8, and show similar patterns of variations to those of upwelling index anomalies (Figure7).The years of weak or strong upwelling index anomalies correspond with the years of weak or strong turbulent mixing anomalies.In the fishing zone of Abidjan, waters appear more turbulent in 2000 and 2001, with anomalies exceeding 18 m 3 s −3 , and more stable in 2008, with values of anomalies up to −30 m 3 s −3 .

Figure 7 . 27 Figure 8 .
Figure 7. Maps of upwelling index anomalies (m 3 s −1 km −1 ).Anomalies are estimated for the period 1998 to 2014.Blue (red) color indicates stronger (weaker) upwelling conditions compared to the 17 year mean.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure 9 .
Figure 9. Interannual variability during the period 1998 to 2014 in (a) Sardinella aurita catch (tons) in year t + 1 and upwelling index (annual mean in m 3 s −1 km −1 ) in year t; and (b) Sardinella aurita catch (in tons) in year t + 1 and wind-induced turbulent mixing (mean June-December in m 3 s −3 ) in year t.

Figure 8 .
Figure 8. Maps of wind-induced turbulent mixing anomalies (m 3 s −3 ).Anomalies are estimated for the period 1998 to 2014.Blue (red) color indicates stronger (weaker) turbulent conditions compared with the 17 year mean.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure 8 .
Figure 8. Maps of wind-induced turbulent mixing anomalies (m 3 s −3 ).Anomalies are estimated for the period 1998 to 2014.Blue (red) color indicates stronger (weaker) turbulent conditions compared with the 17 year mean.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure 9 .
Figure 9. Interannual variability during the period 1998 to 2014 in (a) Sardinella aurita catch (tons) in year t + 1 and upwelling index (annual mean in m 3 s −1 km −1 ) in year t; and (b) Sardinella aurita catch (in tons) in year t + 1 and wind-induced turbulent mixing (mean June-December in m 3 s −3 ) in year t.Note that in (a), the relationship is shown for upwelling index and wind-induced turbulent mixing in year t, and S. aurita catch in year t + 1, as we are investigating the influence of the ocean physics on the S. aurita recruitment success, which is shown in the catch success in following year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure2).

Figure 9 .
Figure 9. Interannual variability during the period 1998 to 2014 in (a) Sardinella aurita catch (tons) in year t + 1 and upwelling index (annual mean in m 3 s −1 km −1 ) in year t; and (b) Sardinella aurita catch (in tons) in year t + 1 and wind-induced turbulent mixing (mean June-December in m 3 s −3 ) in year t.Note that in (a), the relationship is shown for upwelling index and wind-induced turbulent mixing in year t, and S. aurita catch in year t + 1, as we are investigating the influence of the ocean physics on the S. aurita recruitment success, which is shown in the catch success in following year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure2).

Figure 10 .
Figure10.Plots of empirical relationships between interannual variations in Sardinella aurita catch, phytoplankton chlorophyll, phytoplankton timing of initiation, turbulent mixing, and upwelling index.Annual mean chlorophyll concentration in mg m −3 ; timing of initiation of phytoplankton growth in days; June to December mean wind-induced turbulent mixing in m 3 s −3 ; and annual mean upwelling index in m 3 s −1 km −1 .Note that in (a), the relationship is shown for chlorophyll concentration and S. aurita catch both in the same year t, as we are investigating grazing pressure on chlorophyll concentration in that year.Similarly, in (c-f) the relationships are shown for variables in the year t, as we are investigating the influence of ocean physics on the concentration and timing of the first trophic level in the same year, and in (i,j), as we are investigating autocorrelation between physical variables and between biological variables.However, in (b,g,h), the relationship is shown for biophysical variables in year t and S. aurita catch in year t + 1, as we are investigating the influence of the biophysical variables on the S. aurita catch success in following year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure2).

Figure 10 .
Figure 10.Plots of empirical relationships between interannual variations in Sardinella aurita catch, phytoplankton chlorophyll, phytoplankton timing of initiation, turbulent mixing, and upwelling index.Annual mean chlorophyll concentration in mg m −3 ; timing of initiation of phytoplankton growth in days; June to December mean wind-induced turbulent mixing in m 3 s −3 ; and annual mean upwelling index in m 3 s −1 km −1 .Note that in (a), the relationship is shown for chlorophyll concentration and S. aurita catch both in the same year t, as we are investigating grazing pressure on chlorophyll concentration in that year.Similarly, in (c-f) the relationships are shown for variables in the year t, as we are investigating the influence of ocean physics on the concentration and timing of the first trophic level in the same year, and in (i,j), as we are investigating autocorrelation between physical variables and between biological variables.However, in (b,g,h), the relationship is shown for biophysical variables in year t and S. aurita catch in year t + 1, as we are investigating the influence of the biophysical variables on the S. aurita catch success in following year.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure2).

Figure 11 .
Figure 11.Comprehensive assessment of predictive performance of diagnostic models using all possible combinations of input variables.(a-o) Regression analyses between observed Sardinella aurita catch and predicted Sardinella aurita catch based on combinations of four predictive variables, i.e., anomalies of timing of initiation, annual mean chlorophyll concentration, annual mean upwelling index, and June to December mean wind-induced turbulent mixing in year t.The regression statistics are reported using adjusted R 2 , which increases only if added variables are relevant ones.Significances are calculated using F-test.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure2).

Figure 11 .
Figure 11.Comprehensive assessment of predictive performance of diagnostic models using all possible combinations of input variables.(a-o) Regression analyses between observed Sardinella aurita catch and predicted Sardinella aurita catch based on combinations of four predictive variables, i.e., anomalies of timing of initiation, annual mean chlorophyll concentration, annual mean upwelling index, and June to December mean wind-induced turbulent mixing in year t.The regression statistics are reported using adjusted R 2 , which increases only if added variables are relevant ones.Significances are calculated using F-test.All data were averaged over the area of the fishing zone of Abidjan (please see box-area location in Figure2).

Figure 12 .
Figure 12.Conceptual diagram for the year-to-year variations in Sardinella aurita catch in relation to first trophic level and ocean physics along the Ivorian continental shelf.(a) Year when Sardinella catch is high; (b) Year when Sardinella catch is low.Lowercase letters a-j refer to the empirical relationships presented in Figure 10.

Figure 12 .
Figure 12.Conceptual diagram for the year-to-year variations in Sardinella aurita catch in relation to first trophic level and ocean physics along the Ivorian continental shelf.(a) Year when Sardinella catch is high; (b) Year when Sardinella catch is low.Lowercase letters a-j refer to the empirical relationships presented in Figure 10.

Figure A1 .
Figure A1.Maps of annual mean chlorophyll concentration in mg Chl m −3 for the years 1998 to 2014.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure A2 .
Figure A2.Maps of timing of initiation of phytoplankton growth in days for the years 1998 to 2014.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure A2 .
Figure A2.Maps of timing of initiation of phytoplankton growth in days for the years 1998 to 2014.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure A3 .
Figure A3.Maps of annual mean upwelling index in m 3 s −1 m −1 for the years 1998 to 2014.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure A4 .
Figure A4.Maps of annual mean wind-induced turbulent mixing in m 3 s −3 for the years 1998 to 2014.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure A3 . 27 Figure A3 .
Figure A3.Maps of annual mean upwelling index in m 3 s −1 m −1 for the years 1998 to 2014.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure A4 .
Figure A4.Maps of annual mean wind-induced turbulent mixing in m 3 s −3 for the years 1998 to 2014.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.

Figure A4 .
Figure A4.Maps of annual mean wind-induced turbulent mixing in m 3 s −3 for the years 1998 to 2014.The boxed area against the coast indicates the geographical location of the fishing zone of Abidjan.