The anchoveta (Engraulis ringens) is widely distributed throughout the Humboldt Current (4°30′-44°S). In recent years, its eggs and larvae have also been found inside fjords and channels of northern Patagonia, close to the southern limit of the central-south Chilean fishery zone. Currently, it is unclear whether these southern individuals constitute an independent subpopulation. This study analysed the elemental composition of otoliths from 102 specimens from central Chile and northern Patagonia using laser ablation inductively coupled plasma mass spectrometry. The results indicated that the elemental composition of the otolith cores (Mg, Pb, Zn, Ba) differed significantly between sites, revealing the existence of two discrete spawning zones, one in central Chile and one in northern Patagonia. However, the low significant differences of elemental signatures of otolith edges suggest that either individuals from both areas move between spawning areas and mix at certain periods of the year, or they represent pocket units that form part of a larger stock that moves along the coast.
RESUMEN
La anchoveta (Engraulis ringens), está ampliamente distribuida a lo largo de la Corriente de Humboldt (4°30′-44°S). En años recientes, sus huevos y larvas han sido encontradas dentro de fiordos y canales del norte de la Patagonia, cerca del límite sur de la zona de pesca centro-sur chilena. Actualmente, no está claro si estos individuos del sur constituyen una subpoblación independiente. Este estudio analiza la composición elemental de otolitos de 102 especímenes de Chile central y Patagonia norte utilizando Ablación Láser con Espectrometría de Masas con Plasma Acoplado Inductivamente. Los resultados indican que la composición elemental (Mg, Pb, Zn, Ba) de los centros de los otolitos difieren significativamente revelando la existencia de dos zonas discretas de desove, una en Chile central y la otra en la Patagonia norte. Sin embargo, las escasas diferencias en las señales de los elementos en los bordes de los otolitos sugieren que los individuos de ambas áreas pueden moverse entre zonas de desove y mezclarse en ciertos periodos del año, o bien, que ellos representan unidades pequeñas que forman parte de un stock mayor que se mueve a lo largo de la costa.
KEYWORDSLA-ICP-MSanchovyHumboldt Currentsmall pelagic fishesPatagoniaPALABRAS CLAVELA-ICP-MSanchovetacorriente de Humboldtpequeños peces pelágicosPatagoniaINTRODUCTION
The anchoveta, Engraulis ringens (Jenyns, 1842), is a coastal pelagic anchovy that is widely distributed throughout the eastern Pacific Ocean, from northern Peru (4°30′S) to Chiloé Island in southern Chile (44°14′S) (Serra 1983, Pauly and Tsukayama 1987), and is the target of one of the largest pelagic fisheries in the world. This species forms very dense shoals and is strongly influenced by biotic and abiotic environmental factors at all stages of its life cycle (Alheit and Ñiquen 2004). A number of early life history characteristics and reproductive traits, such as egg size, larval size and yolk sac size of recently hatched larvae of the anchoveta, have been described to vary along most of its range of distribution (Llanos-Rivera and Castro 2004, 2006, Castro et al. 2009). At its southern limit of distribution (41°-44°S), however, few studies of biology and ecology of adult anchoveta are available. The southern zone of the species distribution, named as the central-southern Chilean fishery zone for administrative purposes, covers a wide latitudinal range (33°S and 44°S; Fig. 1) in which several historic spawning locations have been identified along the coast (Cubillos et al. 2005, 2006). In recent years, early developmental stages of anchoveta were found in the interior sea of Chiloé Island (northern Chilean Patagonia; 41°-43°S) (Landaeta and Castro 2006, Landaeta et al. 2009, Castro et al. 2015). Asynchrony and significant differences in spawning season have been documented within the central-southern Chilean fishery zone of anchoveta, where spawning occurs from July to December, with higher gonadosomatic indices from August to November below latitude 38°S and from September to October above latitude 38°S (Cubillos et al. 2009, Claramunt et al. 2014). Furthermore, significantly larger females have been found above latitude 38°S (Cubillos et al. 2009).
Evidence suggesting some spawning stock and population connectivity in this central-southern fishery has been proposed. Soto-Mendoza et al. (2012) used individually based models and concluded that a significant northward transport of eggs and larvae occurs from Corral (39°52’S) towards central Chile (36°-37°S). Other methodological approaches, including analyses of morphometrics, electrophoretic protein genetics, DNA and parasites (Ferrada et al. 2002, Valdivia et al. 2007, George-Nascimento et al. 2015), have been used to assess spawning stock structure of anchoveta from northern to central Chile but have not been conclusive (and have sometimes been contradictory) in terms of the number of stock units reported. The stock structure of this fishery remains poorly defined. Furthermore, it is unknown whether the individuals found in the interior/exterior sea of Chiloé Island constitute a separate unit (Zúñiga and Canales 2014) or are part of the central-southern fishery unit that migrates during the annual cycle.
Understanding of fish population structures and determination of discrete populations or spawning stock units is essential for sustainable management of fisheries resources (Sinclair 1988, Sinclair and Smith 2002, Cadrin et al. 2005). Spawning stock identification is a complex task in the pelagic environment. Most common methods for distinguishing between marine populations include analyses of morphometric variation, parasitic fauna, biological markers (gonadosomatic index and sexual maturity curve) and genetic markers (Cadrin et al. 2005). In recent years, trace element composition of otoliths has become a common and effective tool for identifying the spatial heterogeneity of populations and has helped identify migratory patterns and differentiate fish stocks (Hicks et al. 2010, Morais et al. 2010, Avigliano and Volpedo 2016). This technique is based on the fact that otoliths originate during the embryonic stage and grow continuously throughout the life of the fish. As they grow, otoliths show characteristic morphologic and chemical composition patterns. Otolith growth that occurred during the pelagic larval stage can be distinguished from growth during the juvenile or later adult period. The ability to correlate a particular part of an otolith with a discrete point in time is one of the most valuable aspects of this structure (Campana 1999, Chang and Geffen 2013, Górski et al. 2015). Furthermore, otoliths do not seem to experience chemical reabsorption for the majority of elements and therefore reflect the physico-chemical characteristics of water, forming a permanent record of the past environmental conditions experienced by fish (Bath et al. 2000, Rooker et al. 2003, Thorrold et al. 2007).
The objective of this study was to evaluate whether otolith elemental composition could indicate changes or persistence in habitats occupied by anchovetas and thereafter help to determine whether the central-south fishery of Chile can be divided into sub-units corresponding to different areas of spawning and growth or, alternatively, the level of mixing of specimens throughout the area shows that such sub-units did not exist. To address this question, the elemental composition of otoliths from individuals caught off central-southern Chile in 2014 and 2015 was analysed.
METHODSSample locations and fish collection
The study was conducted in the central-southern Chile fishery zone (33°S-44°S) at two locations: Talcahuano (Thno., the central zone) and the Inner Sea of the Chiloé (ISCh, northern Patagonia). Mean satellite seasonal sea surface temperature was obtained along the central-southern Chilean coast, covering both sampling locations for the year of fish collections. Also, for longer term comparisons, an 11-year series (2005-2015) of mean sea surface temperature during the main spawning months (Thno, August-November; ISCh, September-October; Claramunt et al. 2014, Cubillos et al. 2009) in the Talcahuano zone (36.48°S 72.98°W) and at the Inner Sea of Chiloé zone (42.10°S 72.94°W) was obtained from satellite data. The satellite data were obtained from MODIS Aqua Globales (4 km pixel resolution) and downloaded from OceanColor Web (http://oceancolor.gsfc.nasa.gov/), and the images were processed in SeaDas (versión 7.3.2).
At both locations, specimens of Engraulis ringens were obtained from commercial catches in autumn and spring 2014 and summer 2015; only the general catch area was reported, as the exact catch locations are commercially sensitive (Fig. 1). Fish samples were selected at random. Subsequently, sex, weight and total length (TL) were recorded for each fish that was selected for otolith analysis. The fish were sorted into two age groups following the age-length key developed by Araya et al. (2008) and otolith ring counts. Group 1 included adult fish (11.5 and 13.4 cm TL; one year old) and Group 2 included adults larger than 13.4 cm TL (>two years old) (Table 1). Selected specimens were stored in a freezer (–20°C) until analysis at the Fisheries Oceanography and Larval Ecology Laboratory (LOPEL) of the University of Concepción, Chile. Individuals obtained from the central zone and northern Patagonia were similar in size, ranging from 13.1 to 16.1 cm in the central zone and 13.1 cm to 15.6 cm in northern Patagonia. However, during April and January only individuals two years old or older were present, with the largest specimens being observed in April in the central zone (Table 1).
Engraulis ringens fishing zones (central and northern Patagonia) where the samples were obtained.
Month of collection, age group (1, adults 1 year; 2, adults >2 years), catch zone, mean (±S.D.), total length (TL) and number of Engraulis ringens individuals used for otolith microchemical analyses.
Month
Group
central zone
northern Patagonia zone
No. of fish
TL (cm)
No. of fish
TL (cm)
APR 2014
2
8
16.1±1.1
10
15.1±1.6
OCT 2014
1
7
14.7±1.7
10
13.1±0.6
2
10
16.1±1.1
5
15.7±1.4
NOV 2015
1
10
13.1±0.6
2
12.7±0.7
2
11
15.7±1.4
10
15.6±1.4
JAN 2015
2
10
15.7±1.4
9
15.1±1.6
Otolith preparation
After the sagittal otoliths had been extracted from each fish and any adhering tissue had been cleaned away by rinsing with Milli-Q water, the dry samples were stored in 0.5 mL polypropylene micro tubes. The left otolith of each individual was then mounted on a glass slide with epoxy resin and sonicated. The slides were stored in polyethylene bags for later microchemical analysis.
Microchemical analyses
The elemental composition of each otolith was analysed using laser ablation inductively coupled mass spectrometry (LA-ICP-MS) at the Centre for Trace Element Analysis at the University of Otago in Dunedin, New Zealand (the centre uses an ultraclean laboratory for analyses of this type). Depth profiles were obtained along ablation transects from the distal otolith surface (convex side), through the core, to the proximal surface (concave side) (Macdonald et al. 2008, Warburton et al. 2016). This was done using an ICP-MS Agilent 7500cs coupled to an ASI Resolution M-50 laser ablation system powered by a Coherent 193 nm ArF excimer laser.
The otoliths were mounted on a glass slide and the core was located visually through a video imaging system at 400× magnification. Each sample was ablated continuously in a cross section of the unpolished otolith using a 75 micron laser beam. The laser was operated at 5 Hz with an on-sample fluence of 2.2 J cm–2. Ablation occurred in an atmosphere of pure helium to minimize the re-condensation of the materials and potential elemental fractioning (Eggins et al. 1998). For each otolith ablation, data were collected for up to 450 seconds to ensure that the transect completely traversed the core. To reduce the delay time associated with cleaning, no signal smoothing collector was used (Woodhead et al. 2008). Count rate data were collected for the ions of the following isotopes: 24Mg, 31P, 27Al, 43Ca, 55Mn, 60Ni, 66Zn, 88Sr, 63Cu, 85Rb, 138Ba and 208Pb. Concentrations of these elements are usually sufficiently high in the environment (Campana 1999). NIST 612 calibration glass and MACS-3 otolith reference material were run, bracketing groups of four otoliths.
Data analyses
Raw count rate data were reduced to molar ratios to calcium using the IOLITE software package (Iolitesoftware.com version 2.5). The means of the MACS-3 analyses were typically within 5% of the recommended values, with a precision generally better than 3% (relative standard deviation) for most elements. The relationship between Ca normalized concentrations of the elements Mg, Al, P, Mn, Ni, Zn, Cu, Rb, Sr, Ba and Pb was recorded. Mg:Ca and P:Ca concentration are known to decrease throughout the ontogeny of sardine, because of physiological factors influencing the deposit rate of these elements on the otolith (Javor and Dorval 2016). This happened to also be the case in our data from anchoveta otoliths, as Mg:Ca and P:Ca were always strongly correlated over time (Pearson product-moment correlation p<00.1 for each transect). Hence, Mg:Ca and P:Ca often showed a “butterfly pattern” on the otolith transects, with a high concentration at emergence, slowly decreasing over time and mirroring on each side of the nucleus. Therefore, we used the concentration of these elements to locate the nucleus on the transects (Fig. 2).
Mg:Ca and P:Ca ratios in relation to the ablation depth of the otolith of Engraulis ringens. Core and edge regions that were averaged for natal origin analyses are indicated.
The element:calcium ratios were averaged for core and edge regions of each of the studied otoliths. For both edges and cores, 20 s of data were averaged, corresponding to 15 micrometres. Fifteen micrometre averages were used to ensure that only the early larval period was covered in core samples (Plaza and Cerna 2015). Mg:Ca and P:Ca were always strongly correlated and therefore only Mg:Ca was used in the statistical analyses. A permutational multivariate analysis of variance (PERMANOVA, Anderson 2001) was used to determine differences in the element concentrations (Mg, Al, Mn, Ni, Zn, Cu, Rb, Sr, Ba, Pb) of the otoliths between catch zones and period of capture. This provided an analysis of variance for a set of explanatory factors based on dissimilarity measures, thereby allowing differences to be tested at a multivariate level and across a wide range of distributions of empirical data. The analyses were based on Euclidean distance matrices. The tests were performed separately for microchemical concentrations of the otolith edge and core.
Principal components analysis was used to illustrate differences in the multi-elemental composition of cores and edges of otoliths among samples (Mg, Al, Mn, Ni, Zn, Cu, Rb, Sr, Ba and Pb). Subsequently, data of each element were tested for normality and homogeneity of variance using the Shapiro-Wilk and Levene tests (original and log(x+1)-transformed data). Only ratios of Sr:Ca, Mg:Ca, Ba:Ca, Zn:Ca, Rb:Ca and Mn:Ca were normally distributed (Shapiro-Wilk test, p>0.05) for both original and log-transformed data. However, only the Sr:Ca ratio met the assumption of homogeneity of variance (Levene test, p>0.05). Therefore, differences between zones were tested using a non-parametric test (Kruskal-Wallis) followed by a Wilcoxon test. The differences were tested separately for cores and edges.
RESULTS
Differences between locations and among seasons were observed in the mean seasonal sea surface temperature along the central-southern coast of Chile. In the Talcahuano area, temperatures were about 1° to 2°C higher than in the Inner Sea of Chiloé area in all seasons (Fig. 3A). The 11-year time series of mean sea surface temperatures in the peak spawning months of anchoveta at each location also showed that mean temperature values fluctuated inter-annually, but mean temperatures in those months were consistently higher in the Talcahuano area than in the Inner Sea of Chiloé (Thno, 11.2°C-12.5°C vs. ISCh, 10.4°C-11.1°C; χ2 0.05,1, 22.80, p<0.05) (Fig. 3B).
Seasonal satellite-obtained mean sea surface temperature along the central-southern Chilean coast showing the two zones (squares) of fish collection during the study (A). Time series (2005-2015) of mean surface temperature during the peak months of spawning in the two zones (Talcahuano, 36.48°S 72.98°W; Inner Sea of Chiloé, 42.10°S 72.94°W). Shaded area represents standard deviations (B).
This study used a total of 102 adult Engraulis ringens specimens. Otoliths were extracted for laser ablation from 29 specimens that were one-year-old and from 73 that were two years old or older. In total, 56 specimens came from central Chile and 46 from northern Patagonia (Table 1). The concentrations of different trace elements at the otolith edges showed no significant differences in elemental concentrations for any of the factors (zone, month, group and sex; Table 2, Fig. 4). Furthermore, univariate tests for each of the elements separately revealed significant differences only for Mg (p<0.05), which showed higher concentrations in the central zone than in northern Patagonia (Table 3, Fig. 5).
Results of the permutational multivariate analysis of variance based on a Euclidean distance matrix of the standardized concentrations of element/calcium ratios for Mg, Al, Mn, Ni, Zn, Sr, Cu, Rb, Pb and Ba in edges and cores of otoliths from Engraulis ringens caught off central Chile and northern Patagonia (obtained using 9999 permutations of residuals under a reduced model).
Source
Edges
Cores
SS
Pseudo-F
p (perm)
No. perm
SS
Pseudo-F
p (perm)
No. perm
Zone
2E-07
1.3
NS
9912
5E-05
36.7
<0.001
9911
Month
2E-06
3.3
NS
9956
3E-05
6,7
<0.001
9957
Group
2E-07
1.4
NS
9935
5E-07
0.4
NS
9906
Sex
4E-09
0.03
NS
9929
2E-07
0.1
NS
9916
Zone × Month
9E-07
1.8
NS
9952
1E-05
3.4
NS
9952
Zone × Group
2E-07
1.4
NS
9921
2E-07
0.1
NS
9895
Zone × Sex
2E-07
1.1
NS
9951
9E-07
0.6
NS
9905
Month × Group
1E-07
0.7
NS
9927
2E-07
0.2
NS
9902
Month × Sex
5E-07
1.0
NS
9948
4E-06
1.0
NS
9940
Group × Sex
5E-08
0.3
NS
9936
2E-06
1.1
NS
9917
Zone × Month × Group
3E-08
0.2
NS
9925
3E-08
0.02
NS
9912
Zone × Month × Sex
3E-07
0.5
NS
9955
2E-06
0.4
NS
9945
Zone × Group × Sex
6E-08
0.3
NS
9932
3E-07
0.2
NS
9902
Month × Group × Sex
2E-07
0.9
NS
9924
3E-07
0.2
NS
9906
Residual
1E-05
1E-04
Principal component analysis ordination model of the otolith cores and edges with superimposed vectors of the element/calcium ratios. Element codes: Mg, magnesium; Zn, zinc; Pb, lead; Ba, barium; Sr, strontium; Al, aluminium; Mn, manganese; Ni, nickel; Cu, copper; Rb, rubidium. Symbols identify the place of capture.
Post-hoc univariate pairwise Wilcoxon test comparisons of element/calcium ratios in anchoveta otolith cores and edges between central Chile and northern Patagonia. Statistically significant values are indicated in bold.
Element ratio
Edge
Core
Adjusted p-value
Adjusted p-value
Al:Ca
0.191
0.691
Ba:Ca
0.681
1.0e-10
Cu:Ca
0.321
0.510
Mg:Ca
1.4e-6
2.4e-10
Mn:Ca
0.079
0.236
Ni:Ca
0.078
0.226
Pb:Ca
0.121
0.027
Rb:Ca
0.220
0.087
Sr:Ca
0.651
0.432
Zn:Ca
0.816
0.063
Box plot for Mg:Ca, Ba:Ca, Pb:Ca, Zn:Ca and Rb:Ca ratios for the edges and cores of the otoliths in the central zone and northern Patagonia.
The PERMANOVA revealed significant differences in core elemental concentrations among zones and months (PERMANOVA; pseudo-F=36.7, p<0.001, pseudo-F=6.7, p<0.001) but not among age groups and sexes (Table 2, Fig. 4). April was the only month that did not differ between zones in pairwise comparisons (p>0.05; Table 4).
PERMANOVA pairwise tests of differences in multi-elemental concentrations in otolith cores between zones and months of sampling; t and p values were obtained using 9999 permutations of residuals under a reduced model.
Level
Zones
t
p (perm)
Unique perm
APR
central Chile, northern Patagonia
0.1
NS
9953
OCT
central Chile, northern Patagonia
6.5
<0.001
9911
NOV
central Chile, northern Patagonia
4.4
<0.001
9882
JAN
central Chile, northern Patagonia
3.1
0.001
9940
Univariate analyses revealed differences between zones for some of the elements. Specifically, post-hoc comparisons indicated significant differences between zones for Ba, Mg and Pb (p<0.05) and strong trends for Rb and Zn (p<0.09) (Table 3, Fig. 5). Cores of otoliths of specimens from the central zone had higher concentrations of Ba, Mg, Pb and Zn and lower concentrations of Rb than those of specimens from northern Patagonia.
DISCUSSION
The objective of this study was to evaluate whether the elemental composition of otoliths could serve as indicators of changes or of persistence in habitats occupied by the anchoveta Engraulis ringens and, thereafter, to determine whether the central-south fishery of Chile could be divided into sub-units corresponding to different areas of spawning and growth. The results showed significant differences in the elemental composition of otolith cores between the two fishing zones that form the central-south fishery, suggesting that there are two distinct spawning areas and that most individuals spent their early life stages in the area where they hatched. However, low differences in the elemental composition of otolith edges in some short periods indicated the possibility that adults from different origins mixed during the year.
In the otolith cores, the most important differences in terms of elemental composition were higher Ba, Mg, Pb and Zn concentrations in the central Chile zone than in northern Patagonia. In some fish species, magnesium has been positively associated with seawater temperature (Schuchert et al. 2010, Javor and Dorval 2017), which is consistent with the sea surface latitudinal temperature gradient along the Chilean coast, where warmer waters were found off central Chile (>12°C) and colder waters occurred farther south (<11.5°C). This pattern is consistent with our results of differences in sea surface temperatures among seasons and during the 11-year time series during the months of peak spawning at both locations, and has also been reported previously in other studies (Cubillos et al. 2015, Lillo et al. 2015). Similar to patterns observed in juvenile sardines Sardinops sagax (Javor and Dorval 2017), here we report a strong correlation of Mg:Ca and P:Ca ratios for anchoveta Engraulis ringens. These concentrations were strong discriminants between sites with significant water temperature differences, and decreased throughout the ontogeny. Therefore, ontogenetic effects should be always considered when spatial patterns and stock structure of anchoveta are evaluated on the basis of microchemical signatures of otoliths.
We found higher concentrations of Ba in otolith cores from fish collected in central Chile (upwelling zone) than in those collected from northern Patagonia (a non-upwelling zone). Ba has a nutrient-like vertical distribution with much higher concentrations in deeper waters (Lea et al. 1989, Nozaki 2001), and its concentrations in surface waters significantly increase during upwelling events (Hsieh and Henderson 2017). Therefore, the higher Ba:Ca ratios we found in fish in the central zone are most probably an effect of upwelling. These findings are concordant with previous observations of Ba:Ca in seawater and fishes in other upwelling areas where increased Ba:Ca has been detected under upwelling conditions (Wheeler et al. 2016, Hampton et al. 2018). A similar nutrient-like distribution might be expected for Sr (De Villiers 1999), though we found no significant differences in Sr:Ca rations in otoliths between central Chile and northern Patagonia. Wheeler et al. (2016) concluded that individual ratios of Sr:Ca in fish otoliths showed little variation in the California upwelling system and did not serve as valuable proxies of oceanic conditions.
Regarding other elements, off central Chile, Ahumada (1995) found higher metal concentrations (Pb, Zn, Cu, Cr, Ni and Cd) from anthropogenic origins in San Vicente Bay (central Chile, 36°44′S, 73°09′W). Reports from other bays in the same area were consistent with these results. For example, Pb was the predominant trace metal found in the seawater of Concepción Bay (central Chile; 36°40′S, 73°01′W) (Carrera et al. 1993). The greater concentrations of Pb and Zn recorded in the anchoveta otolith cores from central Chile in this study might be a result of high heavy metals recorded for this area traced to nearby iron, steel and petrochemical industries. Zinc concentrations in otoliths have previously been used to indicate particular habitats because they are influenced by fish diet (Ranaldi and Gagnon 2008) and tends to be higher in urban coastal zones due to anthropogenic sources (Sturrock et al. 2012).
Unlike the elemental concentrations of the otolith cores of anchoveta, elemental concentrations of the otolith edges did not differ between zones. There are at least three potential explanations for these core differences and edge similarities in elemental composition. Differences among cores and similar edge concentrations might occur if the adult habitat is environmentally heterogeneous, wide in extension and covering inshore and offshore zones, but the spawning habitat where small larvae occur is restricted to the coastal zone given their limited swimming capacities. These across-shelf ontogenetic differences in elemental composition have been observed in other anchovy species in other systems (e.g. E. encrasicolus in the Gulf of Cadiz; Catalán et al. 2014). Under this scenario, if environmental conditions inshore in central Chile and northern Patagonia differ significantly but are similar offshore, the otolith edge elemental concentrations might not differ between individuals collected in the two zones, while their otolith core concentrations might. Heavy metal concentrations are usually reported to be several times higher in coastal than in oceanic waters due to greater continental contributions from urban areas or sites of high industrial concentration (Luoma 1990, Giordano et al. 1992, French 1993). In our study, besides Ba, these higher concentrations allowed us to infer the presence of two discrete spawning grounds. Prior research cruises (September-October 2009, 2010 and 2012; Cubillos et al. 2010, 2011, 2013) verified that, spatially, spawning grounds occur in coastal areas not exceeding 15 nautical miles offshore. However, larger individuals could be found as far as 30 to 40 nautical miles from the coastal spawning grounds, as depicted from the current industrial fishing areas. Consequently, across-shelf movements would homogenize elemental composition differences in the otolith edges of older fishes.
A second explanation might be that the greater swimming capacity of the larger individuals might allow them to migrate north and south, could they could sometime access other spawning grounds where they would mix with specimens from local spawning and nursery areas. If this movement were frequent and involved large fractions of the populations, differences in the core elemental composition would tend to decrease significantly, and any difference in edge elemental composition (such as Mg, Table 4) indicative of different environmental temperature regimes probably associated with latitude should vanish. In this study during autumn (April), only adult fish older than one year of age were analysed (Table 1), and no differences were observed in core elemental composition between individuals from the two analysed zones. This suggests that adult individuals from both zones might mix before returning to the spawning habitats. This pattern of segregation of population units and a short period of mixing in some seasons is common in a number of fish species displaying homing behaviour associated with estuarine conditions, e.g. Engraulis encrasicolus L. in the Mediterranean Sea (Aldanondo et al. 2010, Guidetti et al. 2013).
A third potential explanation is that these two spawning units represent pocket populations of a broader stock that migrates and that might shift from year to year depending on annual variations in temperature, prey abundance or other undetermined factors. Unfortunately, because the length of our data series is relatively short (one year), we cannot make strong inferences on the restricted or periodical movements across several years. In other systems, elemental differentiation has been observed in cohorts of small pelagic fishes in some years but not in others (Guidetti et al. 2013). While changes in the main spawning locations along the central-southern fishery zone have been reported at 5- to 10-year periods or longer (Bernal et al. 1983, Serra 1983, Castro et al. 1997), shorter-term (inter-annual) changes in migration patterns have not been reported for anchoveta units along their southern distribution range.
In summary, this study of the elemental composition of otoliths provides evidence of two discrete spawning areas of Engraulis ringens along the central-south fishery zone: one in central Chile and one in northern Patagonia, and mixing of adults during some non-spawning periods of the year. This finding agrees with the spatial and structural asynchrony in the reproductive processes in the two areas. Future studies should assess whether this pattern is consistent year to year or varies over longer time periods.
AKNOWLEDGEMENTS
This study was financed by the CONICYT Basal Programme COPAS Sur-Austral PIA PFB31. During the study, the first author was supported by graduate fellowships from CONICYT and COPAS Sur-Austral CONICYT PIA PFR31. LC was also supported by COPAS Sur-Austral CONICYT PIA APOYO CCTE AFB170006 and FONDECYT 1161131. The authors thank J. Betzhold, G. Vásquez, and R. Riffo for collections at the artisanal fishery port and vessels and colleagues at the Laboratorio de Oceanografía Pesquera y Ecología Larval (LOPEL) at the Universidad de Concepción for technical support.
REFERENCESAhumadar.Bahías: Áreas de uso múltiple un enfoque holístico del problema de contaminación1995Número Especial5968AlheitJ.ÑiquenM.Regime shift in the Humboldt Current ecosystem200460201222https://doi.org/10.1016/j.pocean.2004.02.006AldanondoN.CotanoU.TiepoloM.Growth and movement patterns of early juvenile European anchovy (Engraulis encrasicolus L.) in the Bay of Biscay based on otolith microstructure and chemistry201019196208https://doi.org/10.1111/j.1365-2419.2010.00537.xAndersonM.A new method for non-parametric multivariate analysis of variance2001263246https://doi.org/10.1111/j.1442-9993.2001.01070.pp.xArayaM.CubillosL.PeñaililloJ.2008ChileFondo de Investigación Pesquera138Informe final Proyecto FIP/2004-38AviglianoE.VolpedoA.A Review of the Application of Otolith Microchemistry Toward the Study of Latin American Fishes201624369384https://doi.org/10.1080/23308249.2016.1202189BathG.ThorroldS.JonesC.Strontium and barium uptake in aragonitic otoliths of marine fish20006417051714https://doi.org/10.1016/S0016-7037(99)00419-6BernalP.A.RoblesF.L.RojasO.Variabilidad física y biológica en la región meridional del sistema de corrientes Chile-Peru1983291683711CadrinS.FriedlandK.WaldmanJ.2005LondonElsevier Academic Press719CampanaS.Chemistry and composition of fish otoliths: pathways, mechanisms and applications1999188263297https://doi.org/10.3354/meps188263CastroL.QuiñonesR.ArancibiaH.1997Fondo de Investigación Pesquera. Subsecretaría de Pesca235Informe Técnico FIP-IT/96-11CastroL.R.ClaramuntG.KrautzM.C.Egg trait variations in anchoveta Engraulis ringens: A maternal effect to changing environmental conditions in contrasting spawning habitats2009381237248https://doi.org/10.3354/meps07922CastroL.SotoS.LlanosA.2015ChileFondo de Investigación Pesquera414FIP 2013-17 Final Research Report (in Spanish)CarreraM.RodríguezV.AhumadaR.Metales trazas en la columna de agua y sedimentos blandos en Bahía de Concepción, Chile. Determinación mediante voltametría de redisolución199328151163CatalánI.A.Pérez-MayolS.AlvarezI.Daily otolith growth and ontogenetic geochemical signatures of age-0 anchovy (Engraulis encrasicolus) in the Gulf of Cadiz (SW Spain)201415781789https://doi.org/10.12681/mms.819ChangM.GeffenA.Taxonomic and geographic influences on fish otolith microchemistry201314458492https://doi.org/10.1111/j.1467-2979.2012.00482.xClaramuntG.CubillosL.A.CastroL.Variation in the spawning periods of Engraulis ringens and Strangomera bentincki off the coasts of Chile: A quantitative analysis201416096102https://doi.org/10.1016/j.fishres.2013.09.010CubillosL.CastroL.OyarzúnC.2005ChileInstituto de Fomento Pesquero130Informe Técnico Proyecto FIP-IT/2004-03CubillosL.CastroL.ClaramuntG.2006ChileInstituto de Fomento Pesquero147Informe Técnico Proyecto FIP-IT/2005-02CubillosL.PedrazaM.CanalesM.2009ChileInstituto de Fomento Pesquero141Informe final Proyecto FIP/2006-13CubillosL.CastroL.ClaramuntG.2010ChileInstituto de Fomento Pesquero129Informe Final Proyecto FIP/2009-08CubillosL.CastroL.ClaramuntG.2011ChileInstituto de Fomento Pesquero125Informe Final Proyecto FIP/2010-02CubillosL.CastroL.ClaramuntG.2013ChileInstituto de Fomento Pesquero151Informe Final Proyecto FIP/2012-09CubillosL.CastroL.ClaramuntG.2015ChileInstituto de Fomento Pesquero212De VilliersS.Seawater strontium and Sr/Ca variability in the Atlantic and Pacific oceans. Earth Planet1999171623634https://doi.org/10.1016/S0012-821X(99)00174-0EgginsS.KinsleyL.ShelleyJ.Deposition and fractionation processes during atmospheric pressure laser sampling for analysis by ICP-MS1998127-129278286https://doi.org/10.1016/S0169-4332(97)00643-0FerradaS.HernándezK.MontoyaR.Estudio poblacional del recurso anchoveta (Engraulis ringens Jenyns 1842) (Clupeiformes, Engraulidae), mediante análisis de ADN200266243248https://doi.org/10.4067/S0717-65382002000200022FrenchP.Post-industrial pollutant levels in contemporary Svern estuary intertidal sediments compared to pre-industrial levels1993263035https://doi.org/10.1016/0025-326X(93)90594-AGeorge-NascimentoM.OlivaM.Fish population studies using parasites from the Southeastern Pacific Ocean: considering host population changes and species body size as sources of variability of parasite communities20151422535https://doi.org/10.1017/S0031182014001127GuidettiP.PetrilloM.De BenedettoG.The use of otolith microchemistry to investigate spawning patterns in European anchovy: a case study in the Eastern Ligurian Sea (NW Mediterranean)201313914https://doi.org/10.1016/j.fishres.2012.10.015GiordanoP.MusmeciL.CiaralliP.Total content and sequential extractions of Hg, Cd and Pb in coastal sediments19924010421050GórskiK.De GruijterC.TanaR.Variation in habitat use along the freshwater-marine continuum by grey mullet Mugil cephalus at the southern limits of its distribution20158710591071https://doi.org/10.1111/jfb.12777HamptonS.L.MaloneyC.L.van der LingenC.D.Spatial and temporal variability in otolith elemental signatures of juvenile sardine off South Africa2018188109116https://doi.org/10.1016/j.jmarsys.2018.02.001HicksA.ClossG.SwearerS.Otolith microchemistry of two amphidromous galaxiids across an experimental salinity gradient: a multi-element approach for tracking diadromous migrations20103948697https://doi.org/10.1016/j.jembe.2010.07.018HsiehY.T.HendersonG.M.Barium stable isotopes in the global ocean: Tracer of Ba inputs and utilization2017473269278https://doi.org/10.1016/j.epsl.2017.06.024JavorB.J.DorvalE.Stability of trace elements in otoliths of juvenile Pacific sardine Sardinops sagax201657109123JavorB.J.DorvalE.Composition and inter-annual variability in trace element profiles of Pacific sardine otoliths20175895104LandaetaM.F.CastroL.R.Seasonal variability in the spatial patterns of ichthyoplanktonic assemblages from the fjord zone of austral Chile200629107127LandaetaM.MuñozM.I.CastroL.R.Seasonal and short-term variability in the vertical distribution of ichthyoplankton in a stratified fjord of southern Chile2009322742LeaD.W.ShenG.T.BoyleE.A.Coralline barium records temporal variability in Equatorial Pacific upwelling1989340373376https://doi.org/10.1038/340373a0LilloS.LangC.MolinaE.2015ChileInstituto de Fomento Pesquero161Informe Final Proyecto Convenio desempeño 2014LuomaS.Processes affecting metal concentrations in estuarine and coastal sedimentsFurnesR.RainbowP.1990New YorkCRC Press5166Llanos-RiveraA.CastroL.R.Latitudinal and seasonal egg size variations of the anchoveta Engraulis ringens off the Chilean Coast2004102207212Llanos-RiveraA.CastroL.R.Inter-population differences in temperature effects on Engraulis ringens yolk-sac larvae2006312245253https://doi.org/10.3354/meps312245MacdonaldJ.ShelleyM.CrookD.A Method for Improving the Estimation of Natal Chemical Signatures in Otoliths200813716741682https://doi.org/10.1577/T07-249.1MoraisP.BabalukJ.CorreiaA.Diversity of anchovy migration patterns in a European temperate estuary and in its adjacent coastal are: Implications for fishery management201064295303https://doi.org/10.1016/j.seares.2010.04.001NozakiY.Elemental distributionSteeleJ.H.ThorpeS.A.TurekianK.K.2001San Diego, CAAcademic Press840845https://doi.org/10.1006/rwos.2001.0402PaulyD.TsukayamaI.1987ICLARM Studies and Reviews351PlazaG.CernaF.Validation of daily microincrement deposition in otoliths of juvenile and adult Peruvian anchovy Engraulis ringens201586203216https://doi.org/10.1111/jfb.12561RanaldiM.GagnonM.Zinc incorporation in the otoliths of juvenile pink snapper (Pagrus auratus Forster): The influence of dietary versus waterborne sources20083605662https://doi.org/10.1016/j.jembe.2008.03.013RookerJ.SecorD.ZdanowiczV.Identification of Atlantic bluefin tuna (Thunnus thynnus) stocks from putative nurseries using otolith chemistry2003127584https://doi.org/10.1046/j.1365-2419.2003.00223.xSerraJ.R.Proceedings of the Expert Consultation to Examine Changes in Abundance and Species Composition of Neritic Fish ResourcesSharpG.CsirkeJ.Changes in the abundance of pelagic resources along the Chilean coast1983291255284SinclairM.1988Seattle, WashingtonWashington Sea Grant Program252SinclairM.M.SmithT.D.The notion that fish species form stocks2002215297304Soto-MendozaS.ParadaC.CastroL.Modelling transport and survival of anchoveta eggs and yolk-sac larvae in the coastal zone off central-southern Chile: Assessing spatial and temporal spawning parameters201292178191https://doi.org/10.1016/j.pocean.2011.07.001SchuchertP.AlexanderI.ArkhipkinE.Traveling around Cape Horn: Otolith chemistry reveals a mixed stock of Patagonian hoki with separate Atlantic and Pacific spawning grounds20101028086https://doi.org/10.1016/j.fishres.2009.10.012SturrockA.TruemanC.DarnaudeA.Can otolith elemental chemistry retrospectively track migrations in fully marine fishes201281766795https://doi.org/10.1111/j.1095-8649.2012.03372.xThorroldS.ZacherlC.LevinL.Population Connectivity and Larval Dispersal, Using Geochemical Signatures in Calcified Structures2007208089https://doi.org/10.5670/oceanog.2007.31ValdiviaI.ChávezA.OlivaE.Metazoan parasites of Engraulis ringens as tools for stock discrimination along the Chilean coast20077015041511https://doi.org/10.1111/j.1095-8649.2007.01429.xWarburtonM.ReidM.StirlingC.Validation of depth-profiling LA-ICP-MS in otolith applications201674572581https://doi.org/10.1139/cjfas-2016-0063WheelerS.G.RusselA.D.FehrerenbacherJ.S.Evaluating chemical signatures in a costal upwelling region to reconstruct water mass associations of settlement-stage rockfishes2016550191206https://doi.org/10.3354/meps11704WoodheadJ.HellstromJ.PatonC.A guide to depth profiling and imaging applications of LA-ICP-MSSylvesterP.Short Course402008Quebec City, QC, CanadaMineralogical Association of Canada135145ZúñigaM.J.CanalesC.2014ChileInstituto de Fomento Pesquero81Cuarto Informe Proyecto 2.5