Ethmalosa fimbriata (Bowdich 1825), a Clupeid Fish That Exhibits Elevated Batch Fecundity in Hypersaline Waters

Little is known about the concerted influence of temperature and salinity on the fecundity of clupeid fishes. Due to a globally changing climate, both physical parameters might act as stressors, severely affecting the reproductive potential of clupeid fish populations inhabiting tropical estuaries. Differences in relative batch fecundities, the gonado-somatic index, and the condition index of bonga shad (Ethmalosa fimbriata) were analysed in individual females sampled at the Senegalese coast and inside the inverse Sine Saloum estuary, where salinity increases upstream in all seasons. Multiple linear regression models on fecundity and gonadal energy storage show that clupeids can adapt towards increasing their reproductive investment at temperatures (26–30 ◦C) and salinities (42–51), which by far exceed marine conditions, in an effort to maximize recruitment success. This reproductive strategy, however, is accompanied by a trade-off between reproductive effort and somatic growth, which ultimately limits the species’ reproductive potential inside the estuary. The observed high variability in batch fecundities might be a viable mechanism to adjust to fluctuating and rather extreme environmental conditions. Understanding the spawning biology of exploited clupeid fishes in drastically changing environments is crucial for evaluating the reproductive potential of stocks at the outer reach of their physiological performance curve.


Introduction
While spawning traits of marine fish species in marine, brackish, or fresh waters are extensively studied [1,2], little is known about reproductive characteristics of fishes under hypersaline conditions exposed to salinities >40. Such extreme conditions require elevated maintenance costs for osmoregulation and are especially challenging for the mature female [3,4], as well as for its eggs and larvae [5,6]. Spawning as such is metabolically one of the most demanding activities in the lifetime of female fishes [7,8], making reproduction under hypersaline conditions an energetically expensive endeavour [9].
Adaptive phenotypic plasticity is a valuable mechanism for populations to deal with rapid environmental fluctuations such as changes in salinity [10]. Clupeids are generally able to cope with extreme environmental conditions by varying their life history traits via either genetic changes or phenotypic plasticity [11,12]. For example, Baltic sprat (Sprattus sprattus balticus, Schneider 1908) exhibit a decrease in batch fecundities along a declining salinity gradient within the brackish Baltic Sea [13][14][15][16][17].
Temperature not only has significant effects on metabolic rates and fecundities in clupeid fishes, but also on egg and larval development [18][19][20]. Under conditions of extreme temperature and salinity, adult survival rates of marine fishes, as well as egg and larval development, are impaired [21,22]. At the boundary of its physiological performance curve an adult fish is thus forced to budget between survival and reproductive capacity [23]. Therefore, in a heterogeneous environment a species' reproductive potential may exhibit strong variations spatially as well as temporally [19,24].
Little is known about the concerted influence of temperature and salinity on an individual's condition and growth, and consequently on the amount of eggs produced [16,25,26]. Water temperature has been identified as a factor influencing fecundity in marine fishes. The influence can be either passive through its regulatory function of phytoplankton abundances, and thus food availability, and/or active through altering the individual's metabolic rate [8,27,28]. In mature fish, feeding conditions directly affect the individual's condition and the amount of energy that can be put into reproduction [8,29,30]. Furthermore, osmoregulation processes are energetically expensive and may therefore interfere with species' reproductive traits at elevated salinities [4,31]. All of these variables potentially affect subsequent recruitment success and have to be considered when evaluating a species' reproductive potential at the boundary of its physiological optimum range [8,32].
A variety of studies lead to the conclusion that the bonga shad (Ethmalosa fimbriata, Bowdich 1825) is locally adapted to cope with a wide range of environmental conditions in terms of its morphology, growth, and reproductive traits [9,25,[33][34][35][36]. In different areas along the West African coast, E. fimbriata spawns either in the sea [34], or in estuaries and lagoons [36]. A physiology adapted to euryhaline conditions allows the species to cope with varying salinities in these diverse habitats, allowing it to tolerate salinities from around 0 to 97 [9] and even to reproduce in waters with salinities up to 66 [37]. The species is distributed in West African Atlantic waters from Mauritania to Angola [38,39], and despite its marine origin, is highly abundant in estuaries, deltas, and lagoons [37]. It is one of the most important fishes targeted by artisanal fisheries in Cameroon, Ivory Coast, Nigeria, and Senegal [40][41][42].
Especially in Senegalese waters, E. fimbriata has to endure a wide range of varying environmental settings. Former studies demonstrated that the species spawns all year round in the lower and middle reaches of the Saloum River [43,44], an inverse hypersaline environment with monthly fluctuations in temperature and salinity [45]. The rather extreme environmental conditions in the Sine Saloum estuary were used to investigate the reproductive response of a clupeid fish at the outer reach of its physiological performance curve. Due to a globally changing climate the inversion of estuaries in other parts of the dry tropics can be expected in the near future [45][46][47]. Thus, understanding the spawning biology and adaptation potential of commercially important fish species is the key to predict inevitable impacts on the ecosystems and local fisheries.
The current study focused on how this hypersaline environment influences the spawning biology of female E. fimbriata and how the species is adapted to cope with the aforementioned extreme salinity conditions. More specifically, it was attempted to reveal possible energetic trade-offs accompanied by inhabiting and spawning in hypersaline environments. Consequently, variations in somatic condition, gonado-somatic index, and in the batch fecundity of female E. fimbriata were investigated, along a spatial and temporal salinity and temperature gradient at the Senegalese coast and within the hypersaline Saloum estuary.

Environmental Parameters
Night-time surface water temperatures at the three sampling sites steadily increased throughout the sampling period, peaking in October. Surface water temperatures were generally higher at Foundiougne than at Joal, with intermediate values recorded at Djifer. At the beginning of the sampling period (February), salinity was determined to be 39 at Foundiougne as well as at Djifer. In the following months, salinity increased steeply at Foundiougne, eventually peaking in July at 54 until it dropped to 43 in September and October. At Djifer, salinity peaked at 43 in July and August before dropping down to 32 in September. The tropical open ocean's salinity at Joal negligibly fluctuated around 37 throughout the entire sampling period (Figure 1).

Length-Weight Relationship
A total of 1108 female E. fimbriata were sampled for this study; 919 specimens were identified as stage I-IV and the ovaries of 189 individuals stored hydrated oocytes (Table 1). Different power functions were calculated to describe the length-weight relationships of females sampled at the three sampling sites ( Figure 2). Comparison of linear growth models revealed that females sampled inside the Saloum River (Foundiougne) weighed significantly less at a certain length than their counterparts sampled at the river's mouth (Djifer) and at the Atlantic coast (Joal) (Analysis of Covariance (ANCOVA), F(3, 1107) = 3837.104, p < 0.0001; Tukey Honest Significant Difference (HSD), p < 0.001). The length-weight relationship of all sampled females across all sampling sites could be described by the following function (Student's t-test, p < 0.05, r 2 = 0.90, n = 1108):

Length-Weight Relationship
A total of 1108 female E. fimbriata were sampled for this study; 919 specimens were identified as stage I-IV and the ovaries of 189 individuals stored hydrated oocytes (Table 1). Different power functions were calculated to describe the length-weight relationships of females sampled at the three sampling sites ( Figure 2). Comparison of linear growth models revealed that females sampled inside the Saloum River (Foundiougne) weighed significantly less at a certain length than their counterparts sampled at the river's mouth (Djifer) and at the Atlantic coast (Joal) (Analysis of Covariance (ANCOVA), F (3,1107) = 3837.104, p < 0.0001; Tukey Honest Significant Difference (HSD), p < 0.001). The length-weight relationship of all sampled females across all sampling sites could be described by the following function (Student's t-test, p < 0.05, r 2 = 0.90, n = 1108): (1)

The Condition Index and the Gonado-Somatic Index
A statistical analysis revealed not only significant monthly differences in condition index (K) for each station but also significant differences between stations for each month apart from February. At Djifer, the condition of female fish (stage I-IV) peaked in February as well as in July (Analysis of variance (ANOVA), F (8, 286) = 7.50, p < 0.0001; Tukey HSD, p < 0.05). Comparably, the mean K of females sampled at Foundiougne was significantly higher in February and in July than in most other months, except in March and June (ANOVA, F (8, 367) = 10.92, p < 0.0001; Tukey HSD, p < 0.05). Females were in a significantly better condition at Djifer than at Foundiougne in all months except for October (Tukey HSD, p < 0.05). The mean K of females sampled at Joal was significantly higher than the one of individuals sampled at Foundiougne in all months except for March, April, and July (Tukey HSD, p < 0.05). Still, the mean condition of immature females was always lowest in specimens sampled inside the Saloum River's middle reaches at Foundiougne (Figure 3).

The Condition Index and the Gonado-Somatic Index
A statistical analysis revealed not only significant monthly differences in condition index (K) for each station but also significant differences between stations for each month apart from February. At Djifer, the condition of female fish (stage I-IV) peaked in February as well as in July (Analysis of variance (ANOVA), F(8, 286) = 7.50, p < 0.0001; Tukey HSD, p < 0.05). Comparably, the mean K of females sampled at Foundiougne was significantly higher in February and in July than in most other months, except in March and June (ANOVA, F(8, 367) = 10.92, p < 0.0001; Tukey HSD, p < 0.05). Females were in a significantly better condition at Djifer than at Foundiougne in all months except for October (Tukey HSD, p < 0.05). The mean K of females sampled at Joal was significantly higher than the one of individuals sampled at Foundiougne in all months except for March, April, and July (Tukey HSD, p < 0.05). Still, the mean condition of immature females was always lowest in specimens sampled inside the Saloum River's middle reaches at Foundiougne (Figure 3). The condition index (Khyd) of stage V showed a decreasing trend at all sampling sites during the progression of the sampling period (February to October). The mean condition of female spawners was significantly higher in February than in September (ANOVA, F(8, 111) = 3.39, p < 0.01) (see Appendix Table A1 for Tukey HSD results). Comparably, Khyd at Joal (ANOVA, F(2, 30) = 0.46, p = 0.64) and Djifer (ANOVA, F(4, 45) = 1.82, p = 0.14) was higher in March than in all other sampling months ( Figure 4). The condition index (K hyd ) of stage V showed a decreasing trend at all sampling sites during the progression of the sampling period (February to October). The mean condition of female spawners was significantly higher in February than in September (ANOVA, F (8, 111) = 3.39, p < 0.01) (see Appendix A Table A1 for Tukey HSD results). Comparably, K hyd at Joal (ANOVA, F (2, 30) = 0.46, p = 0.64) and Djifer (ANOVA, F (4, 45) = 1.82, p = 0.14) was higher in March than in all other sampling months ( Figure 4).  At Djifer (ANOVA, F (4, 45) = 7.33, p < 0.001) and Foundiougne (ANOVA, F (8, 111) = 2.05, p < 0.05), the mean gonado-somatic index (I G ) of stage V females reached a significant minimum in April and subsequently increased again (Appendix A Table A1). At Joal (ANOVA, F (2, 30) = 0.025, p = 0.98), I G remained constant throughout the three months in which spawning females could be sampled.

Batch Fecundities
Fecund females were sampled in the size ranges (L T ) of 23.1 cm-27.5 cm, 20.8 cm-26.9 cm, and 19.9 cm-28.3 cm at Joal, Djifer, and Foundiougne, respectively. The variance in the number of produced oocytes per spawner was generally high throughout the three sampling sites: absolute batch fecundity (ABF) at Joal ranged between 3316 hydrated oocytes in a 25.4 cm individual and 66,375 ripe oocytes in a 26.1 cm long female, both sampled in May. At Djifer, a 23.3 cm female sampled in April exhibited the lowest ABF with 2200 oocytes, while in February a similarly sized individual (23.5 cm) showed an ABF of 39,490 hydrated oocytes. At Foundiougne ABF ranged between 3160 oocytes in a 25.5 cm female and as much as 85,410 oocytes in a 24.8 cm specimen, both sampled in July. An effect of the female's size on the number of produced oocytes was controlled for by applying ANCOVA on the linear regressions between total length (L T ) and ABF, as well as ovary free body weight (OFBW) and ABF. The assumption of parallel lines was met. The intercept of L T and ABF was significantly higher in females sampled at Foundiougne and Djifer than at Joal (ANCOVA, F (3, 188) = 18.4923, p < 0.0001, Tukey HSD, p < 0.0001). Also, the intercept of OFBW and ABF was significantly higher at Foundiougne than at Djifer and Joal (ANCOVA, F (3, 188) = 16.8477, p < 0.0001, Tukey HSD, p < 0.01). Therefore, ABF values were always highest at Foundiougne for the same female size.
In spite of considerable inter-individual variation, significant monthly variations in relative batch fecundity (RBF) could be observed. Whereas no monthly differences in mean RBF were detected in females sampled at Joal (ANOVA, F (2, 30) = 0.37, p = 0.69), fecundity was significantly lower in February than in April at Djifer. During the consecutive months mean RBF increased again at this sampling site: RBF was determined to be significantly higher in July than in all other sampling months (ANOVA, F (4, 45) = 8.45, p < 0.0001). At Foundiougne a similar pattern emerged; RBF was significantly higher in October than in April (ANOVA, F (8, 111) = 2.04, p < 0.05) (Figure 4; Appendix A Table A1).

Correlation of Reproductive Parameters with Temperature and Salinity
The applied multiple linear regression models revealed concerted effects of temperature and salinity explaining 90%, 25%, and 97% of the observed variances in RBF, K hyd , and I G , respectively. Temperature and salinity had opposing effects on RBF. While temperature had a positive and quadratic effect on the square root transformed data on RBF, the relationship between salinity was negative and quadratic ( Table 2). The apex point of the function was determined for waters warmer than 30 • C and at salinities of around 46. The model results on I G also revealed mirroring effects of temperature and salinity on the ratio of gonad weight to female body weight. In both models the intercept did not pose as a significant term. Conversely, temperature had a linear and negative effect on K hyd in E. fimbriata, while salinity exhibited a negative and quadratic effect on this somatic energy storage parameter. A minimum in K hyd became apparent in females spawning in waters with surface temperatures around 30 • C at a salinity of 46 (Appendix B Figure A4). Table 2. Results of the multiple linear regression models of relative batch fecundity (RBF), condition index (K hyd ), and gonado-somatic index (I G ) against linear and quadratic terms of surface water temperature (T, T 2 ) as well as surface water salinity (S, S 2 ). Estimates of the parameters for the least adequate models, including only significant terms, are shown, together with their respective standard errors and p-values. RSE: Residual Standard Error.

Discussion
Understanding the spawning biology of fish species in extreme environments, such as inverse hypersaline estuaries, is crucial for evaluating the reproductive potential of stocks at the outer ranges of their physiological capacity [8,32]. For the first time in marine fishes, the optimal salinity for egg production was observed to be in hypersaline waters (salinity 42-51). Obtained results on fecundity variations in E. fimbriata suggest that the species is well adapted to extreme fluctuations in environmental conditions through highly variable reproductive output. However, the observed increase in reproductive potential under hypersaline conditions is accompanied by a trade-off between batch fecundity, gonadal energy storage, and somatic growth.
The length spectrum of females sampled for the present study is comparable to the ones obtained in former studies conducted in the same geographical area [25,48]. The observed differences in length/weight relationships revealed that females sampled in the Saloum River's middle reaches (Foundiougne) generally had a lower weight at a certain length when compared to their counterparts sampled at the Saloum River's mouth (Djifer) and at the Senegalese Coast (Joal). Females (stage I-IV) caught at Foundiougne were in a significantly worse nutritional condition (K) in all months when compared to specimens caught at Joal and Djifer. In specimens sampled inside the estuary seasonal differences in K became apparent with notable peaks in February (cool and dry season) and July (warm and dry season). This may hint towards a shift from a marine to a hypersaline phytoplankton community, and thus to temporal differences in food item quality/abundance [49,50]. Also, females spawning in the estuary's middle reaches were exhibiting lower condition indices (K hyd ) (1.37 ± 0.08) than the ones reproducing in the open ocean (1.56 ± 0.12). It has been proposed before that mature fish might travel upstream to spawn, whereas juveniles and immature females may instead tend to move downstream to forage under marine conditions [51,52]. An energetically expensive migratory behaviour towards spawning grounds has been described for the anadromous American shad (Alosa sapidissima, Wilson 1811) [7,53,54] and for the Pontic shad (Alosa immaculata, Bennett 1835) [55]. The expected travelling distances of female E. fimbriata in West African waters, however, are much smaller than those reported for the Alosa species in North America or Eastern Europe and are therefore not necessarily suitable to explain this discrepancy in somatic energy budget. High energetic costs for osmoregulation processes under hypersaline conditions within the middle reaches of the Saloum River and/or differences in stock densities between the sea and the estuary might rather serve as explanations for the observed differences in the condition index [25,56]. Still, elevated stock sizes appear to have negative effects on fecundity in clupeids [16].
Spawning females were observed at Foundiougne during the entire study period, at Djifer from February until July, and at Joal from March until May. These results are in good agreement with the species' temporally limited spawning periods defined in past studies [25,43]. It has been suggested before that seasonally varying hydrographic conditions (changes in water temperature and salinity) play a critical role in the timing and intensity of spawning in clupeid fishes [57][58][59]. Locally varying spawning times therefore hint towards favourable hydrographic conditions for offspring survival at the coast after the onset of upwelling at the beginning of the year. Given the species' prolonged spawning season, however, conditions are most likely favourable for spawning and early life stages survival in the Saloum River's middle reaches all year around. Thus, an examination of E. fimbriata's batch fecundity in relation to respective physical parameters could shed some light on observed energetic disparities and the species' spawning strategy.
High variances in ABF were previously reported for E. fimbriata throughout the entire West African coast [36,60,61]. Our results further underline the species' high phenotypic plasticity in this particular reproductive trait [62,63]. Further, females sampled at Foundiougne produced significantly more hydrated oocytes at a certain length/weight than their counterparts sampled at Djifer and Joal. Thus, females spawning inside the Saloum River's middle reaches were more fecund for the same size. A closer look at the observed spatial and monthly variability in RBF might explain the observed differences in ABF. Comparable to past studies, calculated values for mean RBF also varied widely between months [25,44,52]. Females sampled at elevated salinities were producing four times more oocytes (294 ± 137 oocytes g −1 ) than their oceanic counterparts (74 ± 65 oocytes g −1 ). Clupeid fishes are known to adjust reproductive strategies to local conditions [7,64,65] and thus produce egg quantities and egg sizes that favour offspring survival under the environmental conditions in which they are to develop [66,67].
Length-specific fecundities of experimentally starved Pacific herring (Clupea pallasii pallasii, Valenciennes 1847) were similar to those of fed fish, but the unfed fish had higher weight-specific fecundities, corresponding to a greater loss of somatic tissue during captivity [68]. Consequently, calculated RBF was highest among the unfed/starving fish. In the case of E. fimbriata, elevated RBF as an adaptive response to starvation is unlikely, since no contradictory effect of female size on ABF could be observed. Additionally, primary production within tropical mangrove ecosystems has been amply described as comparably high and thus a shortage in food items is unlikely [69,70].
For estimating the correlation of spawning parameters and water temperature we used satellite-derived sea surface temperatures, as well as in situ surface temperatures, for inland waters. Due to the range in observed surface water temperatures from 18 to 31 • C we accepted the satellite-derived values with an accuracy of ±0.5 • C for our purposes [71]. The multiple linear regression models for RBF and the I G indicate that E. fimbriata is producing increasing numbers of oocytes, and in consequence the ratio of gonad weight to ovary-free body weight is increasing along a rising temperature and salinity gradient. The negative, quadratic effects of temperature and salinity on RBF and I G , on the other hand, illustrate that elevated metabolic rates and osmoregulation processes play a crucial role in limiting the species' reproductive potential under hypersaline conditions. Because K hyd is negatively correlated with temperature and salinity, a trade-off effect between batch fecundity/gonadal energy storage and the female's somatic condition can be concluded. Further evidence is provided by the fact that the mean monthly condition of stage V females decreased significantly while mean RBF increased in concert with rising temperature and salinity gradients throughout most of the sampling period. At Foundiougne, spawner condition only improved again in October when salinity decreased due to the prevailing rainy season. These observations lead to the conclusion that the reproductive investment of well-adapted clupeids is not necessarily dependant on fish size, but rather on the respective spawning environment. High variability in reproductive output through phenotypic plasticity seems to allow fast adjustment to drastically varying environments, and may therefore be beneficial for recruitment in clupeids [7,16].
Embryonic developmental times and hatching rates are not only affected by water temperatures, but also by salinities [73]. While former studies showed that the survival of early life history stages is possible in either high temperature [74] or high salinity regimes [75], egg development, hatching, and larval survival are prone to impairment under combined high temperature and high salinity conditions [21,22]. Certain temperature and salinity combinations might trigger reproduction in clupeids inhabiting variable hypersaline environments, comparable to the isolated effect of temperature on temperate [58], or even on other, tropical species [19]. The observed reproductive strategy may therefore aim towards spawning in a favourable environmental window [76,77], ensuring the survival of well-adapted offspring through enhanced quality in food items and/or predator avoidance during the consecutive weeks [78]. On the contrary, it may also be hypothesised that this strategy might overcompensate for high offspring mortality under these rather extreme conditions [79].
In consequence, the spawning effort combined with elevated metabolic rates under high temperature and energy allocation to the osmoregulation processes under hypersaline conditions are most likely too high to be compensated for by nutritional intake and are therefore at the expense of somatic growth in clupeids [4,20]. This trade-off mechanism will ultimately limit the reproductive potential of clupeid fishes inhabiting estuarine environments impacted by global change. The observed high variability in the reproductive potential of E. fimbriata seems to be a viable strategy in order to adjust to the drastically changing and rather extreme environments such as hypersaline estuaries. Future examinations of food webs in hypersaline environments, in combination with laboratory experiments on egg survival and developmental rates under different temperature and salinity regimes, are necessary to understand the underlying mechanisms. Even further, future application of a Dynamic Energy Budget model may foster our understanding of the observed discrepancies. It would allow for a closer investigation of energy allocations to determine to which degree they are associated with water temperature and salinity in this species [80].

Sampling Sites
Located 100 km south of Dakar (between 13 • 55 and 14 • 10 N, and 16 • 03 and 16 • 50 W), the Sine Saloum delta is comprised of three main branches, the Saloum, the Diomboss, and the Bandiala ( Figure 5). An extended dry season is characteristic for this region, cool from November to March, and warm from April to June. The wet and warm season is rather short, lasting from July to October.
Total rainfall has been steadily decreasing since the 1920's, with a particularly severe decline since 1961. Combined effects of tides, reduced river run-off, intense evaporation, and shallow waters in the estuary's mouth have led to high salinity and an overall inversion of the salinity gradient [45]. Consequently, salinity increases upstream in all seasons in the Saloum River. In the river's upper reaches, surface water temperatures thus exceed 30 • C and salinity levels are two to three times higher than those of seawater by the end of the dry season [81].
Sampling took place at the Senegalese coast and inside the Saloum River from February to October 2014, during the peak of E. fimbriata's spawning season [25,52]. Three different sites were sampled: Joal (Senegalese coast, 14

Environmental Parameters
To get an overview of the seasonal changes of environmental conditions, water temperatures and salinity were measured. Because the upwelling area in front of Senegal (mean water depth on the shelf <20 m) and the Saloum River (mean depth <10 m) are well mixed water bodies, surface water temperatures and salinities were employed [82][83][84]. Moderate-resolution Imaging Spectroradiometer (MODIS), satellite-derived (Aqua, level 2, 0.1 degrees) sea surface temperatures (https://oceancolor.gsfc.nasa.gov/) were assessed for Joal (20 km radius) and Djifer once per sampling week to enhance spatial coverage by avoiding cloud cover. As no remote sensing data for inland waters were available, the Saloum River's surface water temperatures were recorded in situ once per sampling week with a digital thermometer (ama-digit ad 15th; precision 0.4%; accuracy 0.4%). For all stations, salinity was measured with a handheld refractometer (Aqua Medic; precision 0.7%, accuracy 0.2%) using the Practical Salinity Scale (PSS-78). Monthly means in MODIS satellite-derived (Aquarius, level 3, 0.5 degrees) sea surface salinities (https://oceancolor.gsfc.nasa.gov/) were used to double-check the coastal ocean's salinity in a radius of 20 km around Joal.

Fish Sampling
Fish were caught with gill nets (32-36 mm mesh size) by local fishermen, and immediately stored on crushed ice after landing. To ascertain the species' batch fecundity, ca. 1000 fish per sampling site and month were examined in order to find stage V females. Staging of fish was conducted macroscopically [60]. In order to discriminate between mature females with ovaries containing fully hydrated oocytes (stage V) and females that had recently spawned (spent, stage VI), further descriptions on maturity stages in clupeid fishes were consulted [85,86]. Running ripe females were not sampled, since they might have released part of the egg batch. Overall, the ovaries of 189 specimens were dissected and transferred to a 4% borax buffered formaldehyde and freshwater liquid.
Fish were measured by fork length to avoid falsifications due to fin damage. Accordingly, 100 fish with intact caudal fins were measured twice, by fork length as well as total length, to obtain a length coefficient for subsequent length data transformation. Conversions of fish fork length (LF) to

Environmental Parameters
To get an overview of the seasonal changes of environmental conditions, water temperatures and salinity were measured. Because the upwelling area in front of Senegal (mean water depth on the shelf <20 m) and the Saloum River (mean depth <10 m) are well mixed water bodies, surface water temperatures and salinities were employed [82][83][84]. Moderate-resolution Imaging Spectroradiometer (MODIS), satellite-derived (Aqua, level 2, 0.1 degrees) sea surface temperatures (https://oceancolor. gsfc.nasa.gov/) were assessed for Joal (20 km radius) and Djifer once per sampling week to enhance spatial coverage by avoiding cloud cover. As no remote sensing data for inland waters were available, the Saloum River's surface water temperatures were recorded in situ once per sampling week with a digital thermometer (ama-digit ad 15th; precision 0.4%; accuracy 0.4%). For all stations, salinity was measured with a handheld refractometer (Aqua Medic; precision 0.7%, accuracy 0.2%) using the Practical Salinity Scale (PSS-78). Monthly means in MODIS satellite-derived (Aquarius, level 3, 0.5 degrees) sea surface salinities (https://oceancolor.gsfc.nasa.gov/) were used to double-check the coastal ocean's salinity in a radius of 20 km around Joal.

Fish Sampling
Fish were caught with gill nets (32-36 mm mesh size) by local fishermen, and immediately stored on crushed ice after landing. To ascertain the species' batch fecundity, ca. 1000 fish per sampling site and month were examined in order to find stage V females. Staging of fish was conducted macroscopically [60]. In order to discriminate between mature females with ovaries containing fully hydrated oocytes (stage V) and females that had recently spawned (spent, stage VI), further descriptions on maturity stages in clupeid fishes were consulted [85,86]. Running ripe females were not sampled, since they might have released part of the egg batch. Overall, the ovaries of 189 specimens were dissected and transferred to a 4% borax buffered formaldehyde and freshwater liquid.
Fish were measured by fork length to avoid falsifications due to fin damage. Accordingly, 100 fish with intact caudal fins were measured twice, by fork length as well as total length, to obtain a length coefficient for subsequent length data transformation. Conversions of fish fork length (L F ) to total length (L T ) were conducted using the equation L T = 1.0007 × L F + 3.9 (Student's t-test, p < 0.05, r 2 = 0.83).

The Condition Index and the Gonado-Somatic Index
In order to access the nutritional status of immature female fish (Stage I-IV), a condition index (K) Equation (2) was calculated for each individual using wet weight (W W , ±0.1 g), total length (L T , nearest mm) [87], and b of the length-weight relationship. The K hyd (Equation (3)) and the I G (Equation (4)) of hydrated female spawners (Stage V) were based on the ovary-free body weight (W OF , ±0.1 g) and the ovary weight (W Ov , ±0.0001 g) [88]. These indices were used as proxies for the somatic and gonadal energy storage of female spawners.

Absolute Batch Fecundity and Relative Batch Fecundity
E. fimbriata is an indeterminate batch spawner [44]. Therefore, measuring the number of hydrated oocytes produced in a single spawning batch is the only useful method to determine fecundity [89]. ABF was consequently assessed gravimetrically using the hydrated oocyte method for indeterminate spawners [16,90,91]. ABF was divided by W OF to calculate the RBF [13]. Due to high variability in E. fimbriata's batch fecundities, we re-analysed the 10% of the sampled ovaries containing the least amount of oocytes. We measured the diameters of seemingly hydrated oocytes within the questionable samples and compared these values with the diameter values of hydrated oocytes stored in ovaries with high absolute batch fecundities, which were sampled in the same months using ANOVA techniques. Ovaries containing oocytes with significantly lower diameters were rejected. Only hydrated oocytes with a diameter greater than 0.80 mm were taken into account for fecundity analysis. Additionally, ovaries containing oocytes without significant differences in diameter were re-examined macroscopically.

Statistical Analysis
The relationships between L T and W W were described using power functions. Data on L T and W W were log-transformed in order to test for spatial differences in the resulting linear growth models. The assumption of parallel lines was tested and intercepts were subsequently compared using ANCOVA [92]. The parameters of the linear regressions for L T and ABF, as well as for OFBW and ABF, were also compared using ANCOVA. Regression coefficients for the relationships between FL and L T , as well as L T and W W were tested for significance using Student's t-tests. Differences in K were tested using ANOVA, with both sampling month as well as sampling location as independent variables. Additionally, for each sampling location, differences in K hyd , I G , and RBF were tested applying ANOVA with the sampling month as an independent variable. Given the fact that ANOVA techniques are quite robust against small sampling sizes if the assumption of homogeneity of variances is met [93], data on I G and RBF were square root transformed [16,94]. In cases where significant differences were detected, post hoc multiple comparisons were performed. The employed Tukey HSD test calculates p-values that have been corrected for the number of independent pair-wise comparisons that are possible given the number of factor levels. Statistical analyses were carried out using JMP 11.0.1 (SAS Institute Inc., Cary, NC, USA, www.jmp.com).
Since tolerance curves describing physiological processes are a common tool in ecology, predictive models were used, which allowed a comparison between E. fimbriata's spawning preferences and condition optimum [77,95,96]. After interactions between independent variables were tested, we employed multiple linear regressions to reveal the respective responses of RBF, I G , and K hyd to temperature (T, T 2 ) and salinity (S, S 2 ). The resulting minimal adequate models were tested for the constancy of variances (residuals against fitted values; left hand side plots in Appendix B Figures A1-A3) and the normal distribution of their residuals (normal quantile-quantile plot; right hand side plots in Appendix B Figures A1-A3; and Anderson-Darling tests and Cramer-von-Mises tests with a significance threshold of p = 0.1). Where any of these two assumptions were violated, we transformed the response variable and re-ran the multiple regression. This was the case for all three dependent variables, and thus square-root transformations were applied. All final minimal adequate models show constancy of variance and normal distribution of their residuals (Appendix B Figures A1-A3). Models were built in RStudio: Integrated Development for R 1.0.44 (RStudio, Inc., Boston, MA, USA, www.rstudio.com).

Conclusions
Albeit hydrographic conditions in the middle reaches of the Saloum River are favourable for spawning in E. fimbriata all year around, temporal and spatial variations in spawner condition, gonado-somatic index, and batch fecundity could be observed and linked to heterogeneous environmental conditions. Our results show that E. fimbriata has to budget a significant amount of energy towards elevated metabolic rates and osmoregulation processes in an effort to spawn under high temperature/hypersaline conditions. The observed trade-off mechanism between the somatic energy budget and reproductive investment ultimately limits the species' reproductive potential inside the estuary. Due to its high adaptive potential and euryhaline physiology, E. fimbriata is, so far, likely to benefit from the severe impacts of global change on its spawning habitat by potentially outcompeting other pelagic fish species. However, a further decrease in precipitation and a concerted elevation of ambient salinities will have negative impacts on the species' habitat size, reproductive potential, and eventually on its stock size.
Appendix A Table A1. Range in total length (L T ), mean values for condition factor (K hyd ), gonado-somatic index (I G ), relative batch fecundity (RBF), and the range in absolute batch fecundity (ABF) of female E. fimbriata by sampling site and month. Letters indicate Tukey Honest Significant Difference post hoc test results, values not sharing the same letters are significantly different from each other (p < 0.05).

Site
LT (