Larval Fish Community in the Northwestern Iberian Upwelling System during the Summer Period

: The Galician shelf (northwestern Iberian Peninsula) is a highly dynamic area with an important multi-species ﬁsheries industry that exploits resources from several habitats, characterized by being not only highly diverse, rich, and productive but also seasonally and interannually variable. Early life stages of different species are distributed throughout the year, with ﬂuctuating abundances and community composition. Likewise, the inﬂuence of environmental factors and processes on larval production and survival remains unknown. Sampling was carried out in July 2012, and all the larvae obtained were identiﬁed to establish the speciﬁc composition of the community in a summer upwelling scenario. The results show no zonation in the species distribution, a consequence of the mixing effects of the upwelling and eddies, with high diversity but low abundance, which render in a slight predominance of a few species. Due to the dependence of planktonic populations on upwelling events, which was not highly pronounced in 2012, we cannot conclude that this was a typical conformation of the Galician summer larval ﬁsh community, but it is a ﬁrst approach to comprehend the community composition.


Introduction
The fluid environment in which marine populations live offers a wide variety of ways for individuals to disperse within and among populations.The extent of successful dispersal is one of the major determinants of population dynamics but is poorly understood for most marine species.Understanding the drivers of fish larval dispersal is a bio-physical problem that comprises processes that influence offspring production, growth, development, and survival, as well as advection, diffusion, and other physical properties of water circulation and their interactions with larval traits (e.g., vertical migrations), and that operate at various scales [1].
Understanding the influence of environmental variables on fish larval ecology is even more relevant in highly dynamic areas, such as those affected by upwelling events.The Galician shelf, located in the northwestern corner of the Iberian Peninsula, marks the northern boundary of the Iberian-Canary current upwelling system.During summer, wind typically blows southward along the coast, inducing upwelling events and associated southward currents [2].In addition, a subsurface front occurs off Cape Finisterre between two modes (subtropical and subpolar) of the eastern North Atlantic Central Water (ENACW) [3].These modes can mix in an area of convergence, a situation that is intensified by northerly winds [4][5][6][7].
Shoreline also plays a role in upwelling conditions, and changes in orientation between the western and northern coasts modulate wind direction and intensity and are mainly Oceans 2021, 2 701 responsible for the differences observed along the two coasts; upwellings are generally more prevalent in the west and more discontinuous and distant from the coast (occurring near the edge of the continental shelf) in the north [8][9][10].The Galician margin is characterized by a system of terraces up to 30 km wide that form a stepped slope followed by an abrupt lower slope affected by large-scale rotational failures [11], which differs from the surrounding Cantabrian and Portuguese margins.
The predictive models for environmental changes in the area hypothesize a future scenario in which the frequency and intensity of upwelling events will increase due to the increasing northerly winds, especially in the region of the northwestern coast [9].Changes in upwelling frequency or intensity have consequences for ecosystem productivity and composition.For example, the reduction in the intensity and length of upwelling events in the area during the last 40 years has had a significant impact on the abundance, distribution, and species composition of zooplankton, although this influence has been delayed for several years [12][13][14].Linking the environmental conditions to the characteristics of the larval fish community in Galicia could be important for improving the understanding of variations in fish recruitment in present and future environmental scenarios.The cold nutrient-rich deeper water pumped by the upwelling from the ENACW generates a large amount of primary productivity [15], which supports the high amount of fishery and aquaculture activity in this region.The Galician fisheries industry is multi-specific and exploits a resource comprising several habitats (from coastal pelagic to demersal, benthic, or oceanic species), with high diversity, richness, and biomass, and seasonal and interannual variability in abundance and spawning seasons [16,17].Hence, comprehensive larval fish community (LFC) knowledge related to the spatial and temporal structures of the community and the associated environmental factors is required.
Several studies were accomplished in the past, most of them focused on more restricted areas of the same region [18], adjacent regions [19,20], different seasons [21][22][23][24][25], or even with different methodologies [26].This study contributes to the knowledge of the region adding seasonal information of the summer conditions and with a more detailed set of data, given the high level of segmentation of the sample grid in relation with other studies that may comprehend the same area.The ecosystem approach to fisheries has been advocated as one way forward because it considers the holistic complexity of linkages across ecosystems, identifies conflicts between ecosystem services, and considers the impacts of fisheries on marine ecosystems.[27].However, it needs a comprehensive understanding of ecosystem structure and function-the ecosystem approach to fisheries advances as fast as the empirical support of science allows.
In fisheries ecology, the unit of study is the population; communities are composed of fish populations and operate within ecosystems that comprise all other levels of ecological organization [28].Given the lack of successful attempts to manage fisheries based on single populations (maximum sustainable yield), fisheries management has turned its attention to the highest level of ecological organization [29,30].Therefore, there is an interest in matching the practical management scales with those of ecosystem dynamics, while maintaining awareness that there is enormous variability between all the key components of fisheries ecosystems [27].
The causes of variation in recruitment have not led to annual predictive power, however, predictions at a community (or larval fish assemblage) and ecosystem-level dynamics are more powerful [30].Thus, it is critical to understand the composition and distribution of fisheries and other ecological setups related to ecosystem functioning (e.g., phenology, population dynamics, trophic relationships) because it provides relevant information on fisheries management in the context of climate change.
The goal of this study is to describe the composition and structure of the larval fish community in Galicia during a summer upwelling period which will serve as a reference for current and future fishery management efforts.We hypothesize that upwelling has the dominant influence on the horizontal and vertical distributions of the LFC.

Materials and Methods
The CRAMER1207 survey was conducted aboard the Spanish research vessel Cornide de Saavedra along the Galician coast from 17 to 31 July 2012.Ninety-two stations were established in a sampling grid with the stations being distributed along fifteen transects perpendicular to the coastline and separated by eight nautical miles (nm), with a distance between stations of four nm (Figure 1).The transects extended from the 50 m isobath to the 500 m isobath, although due to weather conditions several stations were not sampled.Sampling elapsed for 7 days around the clock, from south to north.Temperature, salinity, and fluorescence were measured at every station with a conductivity-temperature-depth sensor (CTD, SBE25, Seabird Electronics, Inc., Bellevue, WA, USA), coupled to a Turner fluorometer.CTD's casts and net hauls were performed at depths of 200 m at deeper stations or 5-10 m above the bottom at shallower stations.
Oceans 2021, 2, FOR PEER REVIEW 3 for current and future fishery management efforts.We hypothesize that upwelling has the dominant influence on the horizontal and vertical distributions of the LFC.

Materials and Methods
The CRAMER1207 survey was conducted aboard the Spanish research vessel Cornide de Saavedra along the Galician coast from 17 to 31 July 2012.Ninety-two stations were established in a sampling grid with the stations being distributed along fifteen transects perpendicular to the coastline and separated by eight nautical miles (nm), with a distance between stations of four nm (Figure 1).The transects extended from the 50 m isobath to the 500 m isobath, although due to weather conditions several stations were not sampled.Sampling elapsed for 7 days around the clock, from south to north.Temperature, salinity, and fluorescence were measured at every station with a conductivity-temperature-depth sensor (CTD, SBE25, Seabird Electronics, Inc., Bellevue, WA, USA), coupled to a Turner fluorometer.CTD's casts and net hauls were performed at depths of 200 m at deeper stations or 5-10 m above the bottom at shallower stations.In the laboratory, all fish larvae were sorted and identified to the lowest possible taxonomic level.When identification at the species level was impossible, individuals were identified to the family level.Selected taxa (those with ≥20 larvae caught, both day and night) were photographed and, when possible, measured for standard length (±0.01 mm) using ImageJ (v.1.53a).In the laboratory, all fish larvae were sorted and identified to the lowest possible taxonomic level.When identification at the species level was impossible, individuals were identified to the family level.Selected taxa (those with ≥20 larvae caught, both day and night) were photographed and, when possible, measured for standard length (±0.01 mm) using ImageJ (v.1.53a).
Zooplankton was also sampled at each station from 200 m depth through vertical hauls using a CalVET net (25 cm diameter, 53 µm mesh size) equipped with an electronic flowmeter at its mouth.At shallower stations, hauling started 10 m above the bottom.Samples were filtered upon collection through 55 and 200 µ mesh sieves to separate micro (55-200 µm) and mesozooplankton (>200 µm).Both fractions were frozen in liquid nitrogen.In the lab, both fractions were dried in an oven at 60 • C to obtain an estimate of their weight [31].Then, micro-and mesozooplankton weight were standardized with the total volume of filtered water to obtain biomass m −2 .Mesozooplankton abundance (zooAbd) was calculated in the laboratory from the multinet samples for every depth stratum in each station by using a semi-automatic image analysis technique [32].For that purpose, a subsample of 5 mL from each sample was stained with 0.1% eosin for 24 h and scanned; the resulting images were processed using Zoolmage and ImageJ software.

Environmental Data
Chlorophyll fluorescence, temperature, and salinity at 10 m depth (referred to as chlorophyll; Chlor, mg•m −3 ), sea surface temperature (SST, • C), and sea surface salinity (SSS, practical salinity scale) were extracted from the CTD data and used for most statistical analyses.Thermocline was obtained from the CTD data with the R packages oce and rLake-Analyzer.The micro-and mesozooplankton biomass (MiB and MeB) values were referred as mg•m −2 , mesozooplankton counts (Abdzoo) were standardized to ind•m −3 , and depth stratum were also integrated for the sampled water column to obtain a value per station.Dynamic height (DH) was integrated over the water column and calculated from vertical profiles of temperature, salinity, and pressure using the 400 m depth as the reference level of no motion.At shallower stations, or when data were recorded only to 200 m, the density anomaly at the closest 400 m station was assigned to the deepest level sampled by the CTD.As DH and integrated water column density are inversely proportional [33], areas of high DH correspond to low salinity, warm seawater, and anticyclonic regions [34], whereas locations with low DH correspond to salty, cool seawater, and cyclonic eddy regions, showing gradients in frontal regions [35,36].Geostrophic velocities (GVs) were obtained by the first derivative of the DH profiles analyzed on a regular grid of 3 × 3 nautical miles and extracting manually the closest value for each sampling point.Geostrophic velocities were used as an indicator of eddy boundaries and frontal regions because GVs should be higher in these regions [37,38].Finally, spiciness was estimated, which is defined as a state variable and constructed to characterize water masses and indicate double-diffusive stability [39], being higher in warm and salty (spicy) waters [40].Spiciness was calculated using R software v.3.5.1 [41] and the package oce [42].
We used published data on upwelling and wind regimes to understand the environmental scenario prevailing during the survey.The upwelling index and information about eddies was obtained from Instituto Español de Oceanografía [43].Winds regime were obtained from the 'Puertos del Estado' database [44], specifically from the buoy located at Cabo Vilano (43.29 • N and 9.12 • W).
The maps of the physical and biological variables were constructed from the fitted variogram and posterior kriging of the values per station with the R packages gstat [45] and automap [46].

Community
Fish larval abundance was standardized to the number of larvae found beneath a 10 m 2 area of sea surface [ind•10 m −2 ].Similarly, larval densities were standardized to individuals 1000 m −3 and calculated using flowmeter measurements [47].Larval fish diversity (Shannon-Wiener index) and species richness were calculated for each station.Densities were also calculated for each depth stratum.

Horizontal Distribution
Several regression models were tested to assess the influence of the biological and environmental variables (MiB, MeB, Chlor, zooAbd, depth, SSS, SST, DH, and GV) on the parameters of the LFC: larval abundance, diversity (Shannon index), and richness of species at every sampled station.These parameters were obtained integrating the number of larvae at every depth stratum in each sampling station.
Data exploration was carried out after testing for collinearity [48,49].When the correlation between pairs of variables was >0.6, one of the variables was removed from analyses.Finally, the covariables included in the models were depth, SST, SSS, GV, Chlor, MeB, and MiB.Chlor and depth were log transformed to reduce the influence of extremely high values.The fish larval abundance was modeled using the general additive model (GAM), given the nonlinear behavior of some covariates in relation to the response variables.Richness and abundance were adjusted to a negative binomial distribution with a logistic link function, and diversity was adjusted to a normal distribution.To model larval abundance, we added the volume of water filtered at every station as an offset in the equation.The final model was selected following a forward stepwise procedure based on the AIC minimization.The model assumptions were verified by plotting the residuals against the fitted values, each covariate in the model, and each covariate not included in the model.The R package mgcv [50] was used to fit the model.The same procedure was used for the other LFC parameters.
Cluster and ordination methods were used to analyze the structure of the LFC using the matrix of the larval fish abundances, but selected taxa had to have abundances greater than 0.2 larvae•10 m −2 and be present at more than 5% of the stations, resulting in a matrix of twenty-one species.The abundance data were log (x + 1) transformed to dampen the influence of the most abundant species prior to obtaining the dissimilarity (Bray-Curtis) matrix [51,52].R package NbClust [53] was used to find the more appropriated number of clusters, while ANOSIM (analysis of similarities) was performed to test the significance of the clusters.Hierarchical agglomerative clustering with average linking in conjunction with non-metric multidimensional ordination (nMDS) was used to identify assemblages.
The relationship between environmental factors and community structure was assessed with canonical correspondence analysis (CCA).The environmental data matrix included the values of the biological and oceanographic variables at a depth of 10 m.The selection of variables included in the CCA followed a forward stepwise procedure.Significance (p < 0.05) was tested with an ANOVA-like permutation test.Only those variables that significantly explained the species distribution pattern were included in the model.The canonical axes were also tested for significance with the same permutation test.All the ordinations were performed using the vegan package in R [54,55].
The relative length distribution of larvae was explored to assess the direction that the most influential currents were having in the LFC.Lengths were standardized to each species maximum length registered during the entire survey, as a proxy of age to infer the spawning area since smaller specimens will likely be nearer to the spawning area.

Vertical Distribution
A matrix with the densities of species with more than 20 individuals (adding the day and night collections) was created for the analysis of vertical and daily (day/night) distributions.This resulted in a matrix of fifteen species, five depth strata, and two light regimes.The depth strata were 0-20 m, 21-40 m, 41-60 m, 61-100 m, and >100 m.For the light regime analysis, the day period was defined as 7:30 am till 9:30 pm GMT, and the night period was defined as 10:30 pm till 6:30 am GMT, considering the sunrise and sunset time in the study area during sampling.Samples between 6:30-7:30 am and 9:30-10:30 pm were considered as transitional and discarded for this analysis.The weighted mean depths (WMDs) of the larvae in each (MultiNet) haul were calculated as the center of masses of the larval distribution: where n i is the density of fish larvae [ind•1000 m 3 ] in the ith stratum, and d i is the mid-depth of the ith stratum [56].
The amplitude of the diel vertical migrations (DVM) was calculated as the difference between the average WMDs for the day and night periods (DVM= WMD day − WMD night ).Positive values of DVM (type I) correspond to species that move towards the surface during the night, whereas negative values indicate downward movement at night (type II) [57].A T-test was used to test for the significance of the DVM for each species.
The larval densities were used to calculate the Bray-Curtis matrix of distances from the log x + 1-transformed data, these distances were used for the clustering and nMDS.ANOSIM was performed to test if the vertical structure of larval assemblages existed depending on depth and time of day.
Differences between depth strata and day/night abundances were tested using the permutational multivariate analysis of variance using distance matrices (PERMANOVA) with the R package vegan.
Because DVM is often a size-related (ontogenetic) phenomenon, with larvae starting to migrate after yolk sac absorption [57], differences in larval length between day and night and depth strata were also assessed for the selected taxa by using a two-factor ANOVA.When significant differences were found, the Tukey post hoc test was conducted to identify which groups were significantly different from the other.The ANOSIM, NMDS, and distance matrix calculations were done using the R package vegan.
The standard deviation and standard error were used to describe the dispersion of environmental and biological data.

Environmental Variables
There was a north-prevailing wind regime during the sampling period, increasing its intensity towards the end of the cruise, although in the preceding days of the survey, northerly and southerly pulses were alternating.The change of wind direction resulted in an average upwelling index (UI) for the month of July below the historical average value for this period (21 m 3 •s −1 •km −1 and 370 m 3 •s −1 •km −1 , respectively).However, the dominant upwelling conditions prevailed during the cruise and during the four previous days-the mean UI was 191 m 3 .s −1 .km−1 .
The hydrographic structure along the sampled area changed due to the wind regime and coastal orientation.Upwelling events occurred south and north of the grid area during the survey, although they were slightly weaker in the north.In the Artabro Gulf area, the average temperatures were higher than those in the surroundings and linked to the higher GV, which suggests the presence of an anticyclonic eddy.Nevertheless, the upwelling event was not strong enough to completely break the stratification (except partially in the inner shelf), as the water column was stratified in the upper ~60 m, according to CTD's profiles, except when the wind-driven mixed layer reached the top 30 m.
The SST showed an along-shelf temperature front with colder waters occupying the inner shelf; it varied between 13.2 • C next to Cape Finisterre and 19.5 • C west of Cape Ortegal (Figure 2a).The lowest SSS (35.6) were associated with freshwater masses from river runoff near Cape Ortegal, although other coastal low values were found south of Artabro Gulf (35.7).The maximum SSS (35.8) was recorded near the Cape Finisterre over the 200 m isobath (Figure 2d).Temperature and salinity at different depths show cross-shelf gradients.Temperature gradient was from minimum to maximum, from coast, outwards.It ranged from 12.5 • C to 17.8 • C (between bathymetric 10 m and 150 m) and reached its minimum north of Finisterre cape and maximum in the north of the area, over the 200-500 bottom isobaths (Figure 2a,b).Salinity had minute variations with depth, except off Finisterre cape where salinity reached 36.9 at 50 m depth, over the 200-500 isobath (Figure 2e).except off Finisterre cape where salinity reached 36.9 at 50 m depth, over the 200-500 isobath (Figure 2e).A front separates the colder and fresher waters of the inner shelf from the warmer and saltier offshore waters.This front is well reflected by the spiciness, which ranged from 2.3 to a maximum of 5.5 kg.m −3 offshore of Cape Finisterre, that reveals the presence of a subsurface thermohaline front (Figure 2f).
The gradient between the cold freshwater and warm salty water resulted in an area of low dynamic heights along the inner shelf (over grounds <100 m) and a geostrophic current flowing southwestward along the 100 m isobath.The dynamic height ranged from −3.1 to 3.4 cm in the whole area, and geostrophic velocities were weak (<10 cm.s −1 in the upper levels and <5 cm.s −1 at depths under 100 m) but exceeded 30 cm.s −1 at some points near Cape Ortegal (Figure 1).The general currents in the area informed of an anticyclonic mesoscale feature associated with high temperature and salinity, stacking in the gulf during the entire period.A front separates the colder and fresher waters of the inner shelf from the warmer and saltier offshore waters.This front is well reflected by the spiciness, which ranged from 2.3 to a maximum of 5.5 kg.m −3 offshore of Cape Finisterre, that reveals the presence of a subsurface thermohaline front (Figure 2f).
The gradient between the cold freshwater and warm salty water resulted in an area of low dynamic heights along the inner shelf (over grounds <100 m) and a geostrophic current flowing southwestward along the 100 m isobath.The dynamic height ranged from −3.1 to 3.4 cm in the whole area, and geostrophic velocities were weak (<10 cm•s −1 in the upper levels and <5 cm•s −1 at depths under 100 m) but exceeded 30 cm•s −1 at some points near Cape Ortegal (Figure 1).The general currents in the area informed of an anticyclonic mesoscale feature associated with high temperature and salinity, stacking in the gulf during the entire period.

Biological Variables
Regarding the biological variables, the average superficial chlorophyll fluorescence was 8.4 ± 7.7 mg•m −3 , microzooplankton biomass was 1033.1 ± 220.4 mg•m −2 , meso-zooplankton biomass was 2661.5 ± 916.2 mg•m −2 , and mesozooplankton abundance was 3320.4 ± 1472 ind•m −3 .Primary production (Chlor) at 10 m depth was higher near the coast than in the shelf break.The mean Chlor at the stations below the 100 m isobath was 12.3 ± 9.5 mg•m −3 and 5.8.1 ± 4.8 mg•m −3 in stations between 100 and 200 m isobaths (Figure 3a).The MiB maximum (1836.3mg•m −2 ) was located below the 100 m isobaths off the Finisterre coast and in a small area of the Artabro coast (Figure 3b), while the MeB maximum (4990.5 mg•m −2 ) was observed between the 100-200 m isobaths in the Finisterre area and over the 200 m isobath in front of the Artabro coast (Figure 3c).Zooplankton abundance showed its maximum (7405.2ind•m −3 ) between 200 and 500 isobaths north of Finisterre cape and in some stations in the Artabro Gulf (6945.1 ind•m −3 ).Concerning the vertical distribution of chorophyll, the maximum was concentrated in the 20-60 m layer (74.1-0.6 mg•m −3 , max-min), while the maximum zooplankton abundance was mainly in the 0-60 m layer (3635.9-101.2ind•m −3 , max-min).

Biological Variables
Regarding the biological variables, the average superficial chlorophyll fluorescence was 8.4 ± 7.7 mg•m −3 , microzooplankton biomass was 1033.1 ± 220.4 mg•m −2 , mesozooplankton biomass was 2661.5 ± 916.2 mg•m −2 , and mesozooplankton abundance was 3320.4 ± 1472 ind•m −3 .Primary production (Chlor) at 10 m depth was higher near the coast than in the shelf break.The mean Chlor at the stations below the 100 m isobath was 12.3 ± 9.5 mg•m −3 and 5.8.1 ± 4.8 mg•m −3 in stations between 100 and 200 m isobaths (Figure 3a).The MiB maximum (1836.3mg•m −2 ) was located below the 100 m isobaths off the Finisterre coast and in a small area of the Artabro coast (Figure 3b), while the MeB maximum (4990.5 mg•m −2 ) was observed between the 100-200 m isobaths in the Finisterre area and over the 200 m isobath in front of the Artabro coast (Figure 3c).Zooplankton abundance showed its maximum (7405.2ind.m −3 ) between 200 and 500 isobaths north of Finisterre cape and in some stations in the Artabro Gulf (6945.1 ind.m −3 ).Concerning the vertical distribution of chorophyll, the maximum was concentrated in the 20-60 m layer (74.1-0.6 mg•m −3 , maxmin), while the maximum zooplankton abundance was mainly in the 0-60 m layer (3635.9-101.2ind•m −3 , max-min).

Descriptors and Structure of the Larval Fish Community
A total of 2189 larvae belonging to 64 taxa and grouped into 31 families were collected and identified (Table 1).The percentage of unidentified larvae (2.4%) was due to a lack of early life history descriptions of the regional species, damaged specimens, and to a major extent due to the early development stage of the larvae and lack of characteristic features.The abundance per station ranged from 5.7 to 749.6 larvae•10 m −2 , and averaged 217.2 ± 189.1 larvae•10 m −2 , with a maximum of 139 larvae collected in one sampling station (Figure 4a).The most diverse family represented was Sparidae, with six species (one identified to the genus level), and Gobiidae was the most abundant (specimens were grouped into the family level due to the high difficulty in their classification).Two coastal species, Trachurus trachurus (pelagic) and Serranus cabrilla (demersal), were the most abundant and

Descriptors and Structure of the Larval Fish Community
A total of 2189 larvae belonging to 64 taxa and grouped into 31 families were collected and identified (Table 1).The percentage of unidentified larvae (2.4%) was due to a lack of early life history descriptions of the regional species, damaged specimens, and to a major extent due to the early development stage of the larvae and lack of characteristic features.The abundance per station ranged from 5.7 to 749.6 larvae•10 m −2 , and averaged 217.2 ± 189.1 larvae•10 m −2 , with a maximum of 139 larvae collected in one sampling station (Figure 4a).The most diverse family represented was Sparidae, with six species (one identified to the genus level), and Gobiidae was the most abundant (specimens were grouped into the family level due to the high difficulty in their classification).Two coastal species, Trachurus trachurus (pelagic) and Serranus cabrilla (demersal), were the most abundant and ubiquitous species, and the unique mesopelagic (oceanic slope) species that was among those with relative abundance higher than 1% was Maurolicus muelleri (Table 1).

708
Oceans 2021, 2, FOR PEER REVIEW ubiquitous species, and the unique mesopelagic (oceanic slope) species that was am those with relative abundance higher than 1% was Maurolicus muelleri (Table 1).4a).Regarding diversit the community, the highest value was recorded between Cape Prior and Cape Orte and for richness the maximum was registered near the Artabro coast (Figure 4b,c).species richness values ranged from 0 to 20 (11.9 ± 4.7), and the Shannon index ran from 0 to 2.7 (1.9 ± 0.5).
The mean standard length of T. trachurus was 3 ± 1.3 mm, ranging from 1.1 to 7.0 m for S. cabrilla it was 3.7 ± 1.2 mm, ranging from 1.4 to 6.8 mm; for M. muelleri it was 7 ± mm, ranging from 1.6 to 11.2 mm; for C. julis it was 3.2 ± 1.3 mm, ranging from 1.2 to mm.  4a).Regarding diversity of the community, the highest value was recorded between Cape Prior and Cape Ortegal and for richness the maximum was registered near the Artabro coast (Figure 4b,c).The species richness values ranged from 0 to 20 (11.9 ± 4.7), and the Shannon index ranged from 0 to 2.7 (1.9 ± 0.5).
The mean standard length of T. trachurus was 3 ± 1.3 mm, ranging from 1.1 to 7.0 mm; for S. cabrilla it was 3.7 ± 1.2 mm, ranging from 1.4 to 6.8 mm; for M. muelleri it was 7 ± 2.1 mm, ranging from 1.6 to 11.2 mm; for C. julis it was 3.2 ± 1.3 mm, ranging from 1.2 to 10.7 mm.The species that define the summer LFC did not show any horizontal structure (Figure 5); they were not differentially grouped, even among coastal-, shelf-, or slopespawning species as indicated by cluster analysis.The method used to find the best grouping for the dissimilarity matrix identified two clusters as the best result for the matrix of dissimilarities (Figure 5a,c), which was corroborated by the nMDS ordination (Figure 5b).Groups separated at a similarity level of 85%.ANOSIM was performed as well and confirmed this grouping (R statistic = 0.8, p < 0.005).ANOSIM was also used to test the three clusters option, but results were not supported (R statistic = 0.3, p < 0.005).

Family Engraulidae Engraulis encrasicolus
Ee 26.9 2.7 The species that define the summer LFC did not show any horizontal structure (Figure 5); they were not differentially grouped, even among coastal-, shelf-, or slope-spawning species as indicated by cluster analysis.The method used to find the best grouping for the dissimilarity matrix identified two clusters as the best result for the matrix of dissimilarities (Figure 5a,c), which was corroborated by the nMDS ordination (Figure 5b).Groups separated at a similarity level of 85%.ANOSIM was performed as well and confirmed this grouping (R statistic = 0.8, p < 0.005).ANOSIM was also used to test the three clusters option, but results were not supported (R statistic = 0.3, p < 0.005).In general, no differential pattern was evident when mapping the most abundant taxa, except for S. pilchardus (Figure 6b) and E. encrasicolus (Figure 6e), which were concentrated in shallow waters and almost entirely near to Cape Ortegal, respectively.Serranus cabrilla was widespread over the study area but dominating north of Artabro Gulf (Figure 6c), while M. muelleri had no clear distribution pattern (Figure 6f).Trachurus trachurus and C. julis were more abundant in the Artabro Gulf, although C. julis showed a more northward distribution (Figure 6a).
In general, no differential pattern was evident when mapping the most abundant taxa, except for S. pilchardus (Figure 6b) and E. encrasicolus (Figure 6e), which were concentrated in shallow waters and almost entirely near to Cape Ortegal, respectively.Serranus cabrilla was widespread over the study area but dominating north of Artabro Gulf (Figure 6c), while M. muelleri had no clear distribution pattern (Figure 6f).Trachurus trachurus and C. julis were more abundant in the Artabro Gulf, although C. julis showed a more northward distribution (Figure 6a).From the ontogenic point of view, the horizontal distribution of the larvae of the most abundant species did not show any trend, except for S. pilchardus which appears to be concentrated by the coast (Figure 6b).Overall, the horizontal distribution of larvae according to relative length showed that smaller larvae accumulated in the northern area, from the Artabro Coast to Cape Ortegal, while bigger larvae were distributed between Cape Vilán and Cape Ortegal and towards the slope (Figure 7).From the ontogenic point of view, the horizontal distribution of the larvae of the most abundant species did not show any trend, except for S. pilchardus which appears to be concentrated by the coast (Figure 6b).Overall, the horizontal distribution of larvae according to relative length showed that smaller larvae accumulated in the northern area, from the Artabro Coast to Cape Ortegal, while bigger larvae were distributed between Cape Vilán and Cape Ortegal and towards the slope (Figure 7).In contrast, according to the multinet samples results, the vertical distribution of larvae appeared well structured, showing higher abundances in the first 40 m.The PER-MANOVA analysis (R 2 = 0.27, p ≤ 0.001) and Tukey post hoc confirm these differences (Table 2).Mesozooplankton abundance was higher in the first 60 m (min-max, mean: 101.2-3635.9,988.7 ind.m −3 ), both day (101.2-3635.9,913.3 ind.m −3 ) and night (155.0-3566.1,1056.8 ind.m −3 ) (Table 3).Most of the species were almost absent from the 60-200 m depth stratum (except some myctophids) during the day and night.The vertical distribution of the most abundant species was not significantly correlated with mesozooplankton, except for A. thori and S. pilchardus (Table 4).For their distribution in the water column, cluster and ordination analyses (Figure 5d,e) showed four groups (ANOSIM R statistic = 0.95, p < 0.005) of larvae at a similarity level of 40% (Figure 5d,e pilchardus, C. aper, and T. draco.However, the ANOSIM detected only a slight segregation along the water column (ANOSIM R statistic = 0.3, p > 0.005), and no differences in community composition between day and night (R statistic = 0.06 p < 0.005).No significant correlation between the vertical migrations of the pool of fish larvae and the mesozooplankton was observed, except for M. muelleri, C. julis, and S. cabrilla that performed significant type I DVMs (i.e., upward movement at night) (Table 5).The mean vertical displacement was 7.8 ± 6.5 m for the selected species, with M. muelleri being the species that had a wider migration (27 m).In contrast, according to the multinet samples results, the vertical distribution of larvae appeared well structured, showing higher abundances in the first 40 m.The PERMANOVA analysis (R 2 = 0.27, p ≤ 0.001) and Tukey post hoc confirm these differences (Table 2).Mesozooplankton abundance was higher in the first 60 m (min-max, mean: 101.2-3635.9,988.7 ind•m −3 ), both day (101.2-3635.9,913.3 ind•m −3 ) and night (155.0-3566.1,1056.8 ind•m −3 ) (Table 3).Most of the species were almost absent from the 60-200 m depth stratum (except some myctophids) during the day and night.The vertical distribution of the most abundant species was not significantly correlated with mesozooplankton, except for A. thori and S. pilchardus (Table 4).For their distribution in the water column, cluster and ordination analyses (Figure 5d,e) showed four groups (ANOSIM R statistic = 0.95, p < 0.005) of larvae at a similarity level of 40% (Figure 5d,e No significant correlation between the vertical migrations of the pool of fish larvae and the mesozooplankton was observed, except for M. muelleri, C. julis, and S. cabrilla that performed significant type I DVMs (i.e., upward movement at night) (Table 5).The mean vertical displacement was 7.8 ± 6.5 m for the selected species, with M. muelleri being the species that had a wider migration (27 m).Regarding the analysis of length distributions by depth (Figure 8) and day/night period, the results of the two-factor ANOVA showed significant differences in day/night larval length for only two species.In the case of M. Muelleri, larvae caught at night were significantly larger than larvae caught during the daytime, while the inverse pattern was observed for S. cabrilla.Differences in larval length in relation to depth were significant for five taxa; C. julis, M. muelleri, S. pilchardus, and T. trachurus were larger at deeper strata, while the opposite trend was observed for L. crocodilus.The interaction between daytime and depth was significant only for T. trachurus (Table 6).

Relationship between the LFC and Environmental Variables
The regression models constructed to explain the variation in the abundance, diversity, and richness of the fish larvae community in relation to the biotic and abiotic variables explained 66.1%, 59%, and 51.5% of the variability, respectively.In the case of fish larvae abundance, depth, abundance of zooplankton, and GV were the most important drivers followed by SST.For diversity, SSS and GV were the most relevant drivers followed by SST.For richness, only SSS and GV had a significantly influence in setting its variability (Table 7).In general, fish larvae abundance, diversity, and richness increase with increasing temperature, salinity, and geostrophic velocities, while zooplankton abundance has a positive relationship with fish larvae abundance, despite decreasing with depth. --- 6.0 ± 0.6 ns ns ns P. acarne From a multivariate perspective, the CCA resulted in a simplified model (significant p < 0.01) that included three environmental variables: GV, SST, zooplankton abundance (Figure 9).However, the proportion of total variance explained by the environmental constraints was only 18% from which the first two constrained axis explained 11%, and 7%, respectively.Among the species more influenced by temperature were some oceanic myctophids (B.glaciale, M. punctatum, C. maderensis), while the abundance of coastal species, such as S. pilchardus (Figure 6b), was more influenced by higher values of chlorophyll.However, in general, there was no clear pattern or evident association between species and environmental conditions, and none of the variables showed great influence on the species distribution, as reflected by the small percentage of variation explained by the constrains.From a multivariate perspective, the CCA resulted in a simplified model (sig p < 0.01) that included three environmental variables: GV, SST, zooplankton abu (Figure 9).However, the proportion of total variance explained by the environmen straints was only 18% from which the first two constrained axis explained 11%, a respectively.Among the species more influenced by temperature were some myctophids (B.glaciale, M. punctatum, C. maderensis), while the abundance of coas cies, such as S. pilchardus (Figure 6b), was more influenced by higher values of chlor However, in general, there was no clear pattern or evident association between and environmental conditions, and none of the variables showed great influence species distribution, as reflected by the small percentage of variation explained constrains.

Environmental Conditions and Their Influence on the Summer LFC
The abiotic conditions in the Artabro Gulf pointed that this zone was a transition area between two water masses during the summer of 2012, as indicated by the sea surface temperature and salinity and of spiciness which showed a subsurface front off Cape Finisterre.In front of Finisterre cape there was a point which can be reflecting the emergence of water from the subtropical ENACW (Figure 2f), defined as a spicy water [4,15,53].
The upwelling conditions were fully established during the sampling period (predominant northern winds), despite that northerly and southerly wind pulses were alternating in the weeks preceding the survey.Usually, the response time between wind conditions and oceanographic changes in the Galician coast is of 3 days [58].The water column, which had conditions consistent with summer stratification in temperate regions [59], had a mixed layer of variable depth that was narrower in the Artabro Gulf area.Additionally, an anticyclonic eddy in front of the Gulf was constant during the cruise as indicated by the higher geostrophic velocities along the northwestern coast and in a southward surface flow that moved parallel to the coast and changing direction with coastal orientation (Figure 1).
Both CCA and GAMs elucidated the influence of the environment in the structure and descriptors of the LFC, pointing to the upwelling and currents as the main forces conforming the summer community.In both analyses, zooplankton abundance showed a positive relationship with the larvae abundance and the LFC species distribution.
Primary production was higher along the coast than offshore, especially along the Artabro coast next to Ortegal, where the highest superficial chlorophyll concentrations reached 40 mg•m −3 .Regarding secondary production, Fletcher [60] analyzed the egg distribution in the same area and period, and marked the cape next to Ortegal (Estaca de Bares) as the point with the highest egg abundance and Finisterre cape as the point with the lowest egg abundance.These results coincided with the larval distribution seen in this work.An analysis of the LFC in a different environmental regime (spring) also found the same areas to have the highest abundances of eggs and larvae [22], coincident with a cross-shelf frontal region.
This studies about egg distribution plus the distribution of larvae by relative length, point to the Artabro coastal area as a spawning ground for many of the species detected in the present study.This distribution highlights the existence of a current that transports fish eggs southwards from the northern spawning area, and following the clockwise direction of the eddy and the offshore superficial current of the upwelling, which was stronger in the south.Young larvae were gathered by the currents and transported southwestwards with the surface current as they developed.The Artabro coast coincides with the change in orientation of the Galician coast (north to northwest), which makes it more sensitive to northern winds and hence more prone to upwelling events in summer.The presence of a single larval fish assemblage was already seen in spring and then was attributed to offshore Ekman transport associated with a coastal upwelling event, which homogenized the LFC in the cross-shelf direction [22].
For the upwelling areas, the vertical position of fish larvae throughout the day and night also determines if they are retained in shallow and productive waters or advected offshore [61].Larvae with near-surface distributions are more susceptible to offshore transport associated with coastal upwelling than larvae with deeper distributions, which are moved shoreward [61].Thus, larval transport across the shelf depends on the horizontal location of the spawning ground and vertical position at which eggs were spawned.The vertical summer distribution of the LFC in the surface layer suggests that most taxa with neritic spawning (e.g., S. pilchardus, E. encrasicolus, T. trachurus) should spread over the shelf, and oceanic larvae (Sternoptychidae: M. muelleri, Mychtophidae: M. punctatum), which spawn at greater depths, should be driven shoreward by the cold bottom flow [24].It was been noted that the highest densities of eggs were above the mixed layer [60].Eggs from neritic species (S. pilchardus, E. encrasicolus, T. trachurus) were found in higher abundancies at depths lower than 50 m, while eggs from mesopelagic species, such as M. muelleri, were found at depths between 50 and 150 m on the Artabro coast [60].The horizontal distribution of the most abundant species seems to indicate that larvae were transported passively southwards from a spawning area along the northwestern coast.
Regarding the cross-shelf distribution, the abundance of neritic species (Gobiidae, Blenidae, Cupleidae) decreased from the coast to the 200 m isobaths, while oceanic species had the opposite distribution.It had been reported that the spatial distribution of larvae could be quite heterogeneous in coastal areas with a highly variable shelf structure, and that the existence of a shelf slope front associated with a current flow parallel to the coast would contribute to larval concentrations of mesopelagic and shelf species over the shelf break [62], as we observed in this study.Thus, the front detected in the Galician shelf in the summer of 2012 could act as a barrier preventing offshore dispersion.The importance of frontal regions for fish spawning and the concentration of fish larvae has already been documented [30,63].Settlement and recruitment of the coastal species will be favored by retention structures; here lies the importance of the knowledge of seasonal environmental conditions to understand the mechanisms that affect fish larval stages.

Seasonal LFC in Galician Waters
The LFC was highly diverse and rich, with lower abundances than in winter or spring and without dominant species [23,24].The LFC was made up of a mixture of coastal, neritic, and oceanic mesopelagic species, which coincides to a great extent with the composition and diversity of coastal and transitional assemblages found in the central Cantabrian Sea in summer [52].The Galician summer assemblage seems to differ from the spring assemblage.Diversity was slightly higher in the summer of 2012 compared with spring (summer: 0-2.7, mean 1.9, spring: 0-3.1, mean 1.8) but richness lightly lower (summer: 0-20, mean 11.9, spring: 1-32, mean 13.8), and abundance considerably lower (summer: 5.7-749.6,mean 217.2 larvae•10 m −2 , spring: 4.6-26147.3,mean of 939.8 larvae 10m −2 ) [25].However, due to the markedly different seasonal environmental conditions and the impact that it has on communities and poor knowledge about its dynamics, it is impossible to evaluate if such differences or similarities are permanent or due to specific seasonal fluctuations or to species protracted reproductive period or adaptions.
In the literature, there are some descriptions about the spawning season for some of the dominant species.The spawning period of T. trachurus begins in March-April and ends November [64], while for S. cabrilla it begins in February and ends in July in southern latitudes [65].Regarding E. encrasicolus, its spawning season increases with decreasing latitude from Bay of Biscay to Gulf of Cádiz in April-August to April-November [66].Sardina pilchardus has the main spawning period during October-June, peaking in December-March, but with longer duration and earlier peak at lower latitudes [67].Unfortunately, there are no reproductive studies for most species-including C. julis, P. acarne, P. pilicornis, and T. draco in the Galician waters.Regarding the alongshore composition, there are different structures of the LFC between the northern and northwestern coasts of Galicia in the spring of 2012 [23], which were not evident in the summer.Yet, more data should be collected to obtain more robust conclusions.

Conclusions
The larval fish community of the Galician coast seems to be shaped by the interaction between fish spawning location, concurrence of two water masses, coastal upwelling, and an anticyclonic eddy.This interaction resulted in a lack of structure in the across and along the shelf.The eddy circulation may have transported fish eggs and larvae from a spawning area in the north, which with the help of winds (and consequent upwelling) pushed larvae southwards and offshore while avoiding oceanic dispersion.This scenario kept the larvae retained in the same area, either if they were spawned by shelf-, coastal-, or slope-spawning species.Similarly, changes in the seasonal environmental conditions could result in major changes affecting interspecific competition and reproductive success, which finally would be reflected in recruitment and stock abundance.

Figure 1 .
Figure 1.Distribution of the sampled stations (black circles).The pink-shaded area is the Artabro Gulf, an important spawning area for many commercial species.The main capes in the area are shown in magenta.Blue arrows are the geostrophic velocity of currents at each sampling station.Ichthyoplankton sampling was performed with a MultiNet (Midi, 0.50 m 2 mouth opening) consisting of 5 nets of 200 μm mesh.The MultiNet was equipped with a Scanmar depth sensor as well as an electronic flowmeter located at the mouth.The depth strata, defined by the different depths at which each net was open, were 0-20 m, 21-40 m, 41-60 m, 61-100 m, and 100-200 m.At the stations located in depths shallower than 200 m the programmed terminal depth was adjusted and/or the number of strata reduced.The net was towed obliquely at 2 knots.The duration of the hauls was approximately one hour, and the mean volume filtered for each depth stratum was 25 L for the first three strata, 57 L for the fourth stratum, and 87 L for those at depths deeper than 100 m.Samples were preserved in a 4% seawater/formalin buffered solution with borax.In the laboratory, all fish larvae were sorted and identified to the lowest possible taxonomic level.When identification at the species level was impossible, individuals were identified to the family level.Selected taxa (those with ≥20 larvae caught, both day and night) were photographed and, when possible, measured for standard length (±0.01 mm) using ImageJ (v.1.53a).

Figure 1 .
Figure 1.Distribution of the sampled stations (black circles).The pink-shaded area is the Artabro Gulf, an important spawning area for many commercial species.The main capes in the area are shown in magenta.Blue arrows are the geostrophic velocity of currents at each sampling station.Ichthyoplankton sampling was performed with a MultiNet (Midi, 0.50 m 2 mouth opening) consisting of 5 nets of 200 µm mesh.The MultiNet was equipped with a Scanmar depth sensor as well as an electronic flowmeter located at the mouth.The depth strata, defined by the different depths at which each net was open, were 0-20 m, 21-40 m, 41-60 m, 61-100 m, and 100-200 m.At the stations located in depths shallower than 200 m the programmed terminal depth was adjusted and/or the number of strata reduced.The net was towed obliquely at 2 knots.The duration of the hauls was approximately one hour, and the mean volume filtered for each depth stratum was 25 L for the first three strata, 57 L for the fourth stratum, and 87 L for those at depths deeper than 100 m.Samples were preserved in a 4% seawater/formalin buffered solution with borax.In the laboratory, all fish larvae were sorted and identified to the lowest possible taxonomic level.When identification at the species level was impossible, individuals were identified to the family level.Selected taxa (those with ≥20 larvae caught, both day and night) were photographed and, when possible, measured for standard length (±0.01 mm) using ImageJ (v.1.53a).Zooplankton was also sampled at each station from 200 m depth through vertical hauls using a CalVET net (25 cm diameter, 53 µm mesh size) equipped with an electronic flowmeter at its mouth.At shallower stations, hauling started 10 m above the bottom.Samples were filtered upon collection through 55 and 200 µ mesh sieves to separate

Figure 2 .
Figure 2. Temperature at depths of (a) 10 m, (b) 50 m, and (c) 100 m, salinity at depths of (d) 10 m and (e) 50 m, and (f) spiciness.The plot of salinity at 100 m is not shown due to its extremely low standard deviation (2.9).

Figure 2 .
Figure 2. Temperature at depths of (a) 10 m, (b) 50 m, and (c) 100 m, salinity at depths of (d) 10 m and (e) 50 m, and (f) spiciness.The plot of salinity at 100 m is not shown due to its extremely low standard deviation (2.9).

Figure 4 .
Figure 4. Distribution of the LFC parameters of (a) abundance, (b) diversity, and (c) richness.The maximum abundance (749.6 larvae•10 m −2 ) of larvae was recorded on the Artabro coast and it was mainly due to the contribution of the family Gobiidae (161.8 larvae•10 m −2 ) and non-identified larvae (275 larvae•10 m −2 ).The second highest value of abundance (668.3 larvae•10 m −2 ) was found south of Artabro Gulf, where T. trachurus contributed the most to the total larvae abundance (244 larvae•10 m −2 ) (Figure4a).Regarding diversity of the community, the highest value was recorded between Cape Prior and Cape Ortegal and for richness the maximum was registered near the Artabro coast (Figure4b,c).The species richness values ranged from 0 to 20 (11.9 ± 4.7), and the Shannon index ranged from 0 to 2.7 (1.9 ± 0.5).The mean standard length of T. trachurus was 3 ± 1.3 mm, ranging from 1.1 to 7.0 mm; for S. cabrilla it was 3.7 ± 1.2 mm, ranging from 1.4 to 6.8 mm; for M. muelleri it was 7 ± 2.1 mm, ranging from 1.6 to 11.2 mm; for C. julis it was 3.2 ± 1.3 mm, ranging from 1.2 to 10.7 mm.

Figure 5 .
Figure 5. (a) Dendrogram, (b) NMDS ordination and (c) map of the clustering of the horizontal species composition, (d) dendrogram, and (e) nMDS ordination of the vertical species composition.

Figure 5 .
Figure 5. (a) Dendrogram, (b) NMDS ordination and (c) map of the clustering of the horizontal species composition, (d) dendrogram, and (e) nMDS ordination of the vertical species composition.

Figure 8 .
Figure 8. Length distribution of fish larvae at varying depths.Line: median; box: 1st and 3rd quartiles; whiskers: maximum and minimum values; black dots: outliers.The results of the Tukey post hoc test (α= 0.01) for differences in larval length of the most abundant taxa by depth are indicated by the letters a and b to define the groups that are significantly different from each other; groups sharing a letter are not significantly different.

Figure 8 .
Figure 8. Length distribution of fish larvae at varying depths.Line: median; box: 1st and 3rd quartiles; whiskers: maximum and minimum values; black dots: outliers.The results of the Tukey post hoc test (α= 0.01) for differences in larval length of the most abundant taxa by depth are indicated by the letters a and b to define the groups that are significantly different from each other; groups sharing a letter are not significantly different.

Figure 9 .
Figure9.Canonical correspondence analysis (CCA) biplot for the environmental variables ( larval fish taxa (red labels), and sampling stations (grey dots).Arrows indicate the rela portance (length) and correlation (angle with axis) between each variable retained in the mo the canonical axes.The significant environmental variables were temperature (Temp), zoop abundance (Abdzoo), and geostrophic velocity (GeosV).

Figure 9 .
Figure9.Canonical correspondence analysis (CCA) biplot for the environmental variables (arrows), larval fish taxa (red labels), and sampling stations (grey dots).Arrows indicate the relative importance (length) and correlation (angle with axis) between each variable retained in the model and the canonical axes.The significant environmental variables were temperature (Temp), zooplankton abundance (Abdzoo), and geostrophic velocity (GeosV).

Table 1 .
List of larval fish species grouped by families collected off the Galician coast during the Cramer1207 research survey.The families are ordered by decreasing abundance.Code: taxon code; Oc: percentage of occurrence (%); RA: relative abundance (%).The taxa written in blue are those with relative abundances >0.2 larvae 10•m −2 and that were present in >5% of the stations, therefore were included in subsequent analyses.

Table 2 .
Tuckey post hoc test for the differences between larval abundances, per strata (p < 0.05).

Table 2 .
Tuckey post hoc test for the differences between larval abundances, per strata (p < 0.05).

Table 4 .
Pearson's correlation coefficients between the vertical distributions (weighted mean depth, WMDs) of fish larvae and mesozooplankton, where the * stands for p < 0.01.

Table 5 .
Results of the day and night weighted mean depth analysis and respective t-test analysis.The data show the diel vertical migration (DVM) of the dominant species.Positive DVM indicates species ascending at night and descending during the day (DVM type I).Negative DVM indicates species descending at night and ascending during the day (DVM type II).The ** stands for p < 0.01 and * < 0.1.

Table 6 .
Mean (± standard error) of larval length (mm) and number of individuals per depth stratum and time (D: day, N: night).The results of the two-factor ANOVA for differences in the larval length of the most abundant taxa between time and depth are also shown.Legend: ns not significant.

Table 7 .
Parameters of the best regression model for explaining the relative influence of the environmental variables on the larval fish community descriptors.Abiotic (logarithm of depth (logDepth), sea surface temperature (SST), sea surface salinity (SSS), geostrophic velocity (GV), and biotic variable (logarithm of chlorophyll (LogChlor) and zooplankton abundance (Abdzoo)) had an influence on diversity (Shannon-Wiener index), species richness, and abundance of fish larvae.Legend: Coeff = regression coefficents, SE = Standard error, p = p-value, Edf = Estimated degrees of freedom for the predictor with smoothers.