REVISTA BIO CIENCIAS

Materiales y Métodos

Methods and Materials

Se realizaron cuatro arrastres superficiales zooplanctónicos diurnos, estacionalmente de   primavera a invierno del 2002, se utilizó un máximo de cinco sitios de muestreo que estuvieron situados a lo largo del canal principal (11.5 km) mediante la metodología estándar descrita por Smith y Richardson (1977);  utilizando una lancha con motor fuera de borda y  una red tipo Zeppelín (cónica), de luz de malla 505 µm por 1.50 m de longitud y 0.60 m de diámetro de la boca, equipada con flujómetro digital para medir el volumen de agua filtrada. Los arrastres fueron realizados con una duración de 10 minutos y bajo condiciones de marea alta, apoyándose en las tablas de marea emitidas por el Departamento de Oceanografía Física del CICESE (Centro de Investigación Científica y Educación Superior de Ensenada) y a 10 cm por debajo de la superficie del agua para evitar tomar materia orgánica suspendida, de forma simultánea, en cada sitio de muestreo se obtuvieron los datos de temperatura y salinidad superficial por medio de un termómetro de inmersión graduado (0.1 °C de precisión) y un refractómetro de campo (ATAGO s/mill-E, de 0-100 ‰  al 1 ‰  ups de precisión), respectivamente. El material colectado fue preservado con formol al 4% y 20 mL de borato de sodio. Las larvas de peces fueron separadas del resto de los organismos y colocadas en frascos (400 y 500 mL) debidamente etiquetados, las larvas fueron identificadas hasta el nivel de especie cuando era posible. Las principales fuentes bibliográficas utilizadas fueron Sumida et al., (1985), Moser et al., (1993) Moser et al., (1994), Moser (1996), Watson (1996), Leis y Carsons-Ewart (2000). Finalmente, las larvas identificadas se contabilizaron  y estandarizaron a un volumen de 1000m³ mediante la siguiente fórmula:

Donde E=abundancia estandarizada, n=número de organismos y  v=volumen filtrado del agua. Se calculó el índice de diversidad y equidad de Shannon H’ (Magurran, 1988) para cada uno de los sitios de muestreo y las estaciones del año, mediante la fórmula H’=-S Pi ln Pi, donde Pi  la abundancia por especie proporcional del total encontrado. Asimismo fue aplicado el análisis de cluster entre las variables: especies-sitios de muestreo y especies-estaciones del año, con la finalidad de conocer el grado de afinidad o similitud entre éstas (Ludwing y Reynolds 1988), el software utilizado fue el Peimer-ELtd.

Four daytime superficial zooplanktonic trawls were performed seasonally from spring to winter 2002; a maximum of five sampling sites were used located along the main channel (11.5 km) by using the standardized methodology described by Smith y Richardson (1977), utilizing a boat with outboard motor and a Zeppelin type net (conic), a mesh size of 505 µm with a length of 1.50 m and mouth width of 0.60 m which was equipped with a digital flow meter to calculate the volume of filtered water. Trawls were made with a duration of 10 minutes and under high tide conditions, based on the tide tables issued by the Physics Oceanography Department of CICESE (Scientific Investigation Center and Superior Education of Ensenada) and at 10 cm under water surface to avoid taking organic suspended matter, simultaneously, in each sampling site, temperature and superficial salinity data were obtained through a graduated immersion thermometer (precision of 0.1 °C) and a field refractometer (ATAGO s/mill-E, from 0-100 ‰  to 1 ‰  psu accuracy) respectively. Collected material was preserved with formaldehyde at 4 % and 20 mL of sodium borate. Fish larvae were separated from the rest of the organisms and placed in jars (400 and 500 mL) duly labeled; larvae were identified with the species level when possible. Main bibliographic references used were Sumida et al., (1985), Moser et al., (1993) Moser et al., (1994), Moser (1996), Watson (1996), Leis and Carsons-Ewart (2000). Finally, identified larvae were counted and standardized to a volume of 1000 m³ through the next formula:

Where E=standardized abundance, n=number of organisms and  v=filtered volume of water. The Shannon H’ diversity and equity index (Magurran, 1988) was calculated for each of the sampling sites and the year seasons, by means of the formula H’= -S Pi ln Pi, where Pi abundance per proportional species from the total found. In addition, a cluster analysis was applied between the variables: species – sampling sites and species – seasons of the year, in order to know the degree of affinity or similarity between them (Ludwing and Reynolds 1988), the software used was Peimer- ELtd.

Resultados y Discusión

Results and Discussion

La variación de la temperatura y salinidad promedio por estación del año registró en primavera el valor más alto de la temperatura (34.6 °C) y el más bajo fue para el verano (23.4 °C), en cambio para la salinidad el valor más alto se presentó en verano (34.6 ups) y el más bajo en invierno (24.2 ups), cabe mencionar que durante el otoño no se obtuvieron registros (Figura 2a). La variación de dichos parámetros  por sitio de muestreo, presentaron homogeneidad entre sí, aún así el valor mayor (30 °C) de temperatura promedio fue registrado en los sitios  4 y 5, y el menor en el sitio 2 (28.3 °C), sin embargo, las fluctuaciones de la temperatura entre sitios fue menor a un grado, misma situación fue observada en la variación de la salinidad entre sitios (variación de 0.5  ups), aún cuando el valor mayor (30.8 ups)  se presentó en el sitio 2, en tanto que el valor menor (29.5 ups) fue registrado en el sitio 3 (Figura 2b).  El elenco ictioplanctónico correspondiente al período de estudio (2002) estuvo representado por 11 familias, con 13 géneros y nueve especies. Cuyas densidades estuvieron representadas principalmente por  Engraulis mordax (215.28 org/1000 m3), con una abundancia relativa de 40.2 %, seguida de Eucinostomus sp. (133.78 org/1000m³),con el 22.68 %, en tanto que Dormitator  latifrons (76.88 org/1000m³), presentó el 13.14 % y finalmente Ophistonema libertate y  Achirus mazatlanus con bajas densidades (49.59 org/1000 m³  respectivamente), presentaron abundancias relativas de 7.47 y 7.21 % respectivamente. Las densidades del resto de las especies fluctuaron de 33.83 a 1.53 org/1000 m³  cuyas abundancias relativas registradas fueron de 5.66 a  0.25%  (Figura 3, Tabla 1).

Referente a la variación en tiempo, otoño fue la estación del año en la que se registró la mayor densidad larval (269.1 org./1000 m³), seguido de primavera y verano (175.3 y 127.63 org./1000 m³ respectivamente), en tanto que invierno apenas si alcanzó una mínima densidad (29.21 org./1000m³) (Figura 4).

En lo que respecta a la variación espacial el sitio 3 fue el punto con mayor densidad larval (176.83 org./1000m³), seguido de los sitios 5 y 2 (155.31 y 141.47 org./1000m³ respectivamente), en tanto que los sitios 1 y 4 reportaron las densidades más bajas de todo el periodo de estudio (72.27 y 55.35 org./1000 m³ respectivamente) (Figura 5). Cuya mayor densidad larval fue registrada para  E. mordax y A mazatlanus en el sitio 3,  con un dominio absoluto de E. mordax en los sitios 5 y 2, no así para los sitios 1 y 4 cuyos mayores registros fueron para las especies D. latifrons y Eucinostomus  sp. respectivamente. Las variaciones espacio-temporales en la temperatura, salinidad y abundancia del ictioplancton indican que la abundancia larval se encuentra en estrecha relación con el patrón estacional influyendo  a su vez en el ciclo reproductivo de estas especies, ya que las mayores abundancias se presentaron en otoño y primavera (269.1 y 175.3 org./1000m³ respectivamente) con un amplio intervalo de temperatura y salinidad (29.1 °C a 34.6 °C y hasta 32.8 ups) con valores mayores de ambas variables en primavera, no así para el verano, lo que sugiere la influencia de la época de lluvias, abundantes en el otoño y nulas en primavera. Donde las mayores densidades larvales se presentan con amplios intervalos de temperatura (29 °C a 34.6 °C) y salinidad (34.6 ups), Margalef (1969), Livingston (1984) y Contreras Espinoza (1993)  señalan que en las lagunas costeras el zooplancton en general se asocia a la variabilidad del sistema en términos de salinidad y temperatura; Por su parte Navarro-Rodríguez et al., (2006), señalan que en este mismo sistema (2001) las variaciones de la abundancia espacio temporal estuvieron influenciadas por las fluctuaciones de temperatura (29.8-31.6 °C) y la salinidad (11-34.8 ups), presentándose altas concentraciones en otoño y menores en primavera asociadas a las épocas frías y cálidas respectivamente. Indicando además que los picos de la época reproductiva de la mayoría de las especies (nueve en primavera y siete en otoño) principalmente la de E. mordax (especie más abundante) se presentan principalmente en una época cálida con abundantes lluvias. Y otro más en la época de estiaje, asociándose posiblemente a las altas concentraciones de salinidades (32.8 ups) puesto que la gran mayoría de las especies son de hábitos marinos, lo que además puede asociarse también al transporte por las corrientes, patrones de circulación locales así como a los efectos de mareas. Lo que sugiere que posiblemente  se presenta el mismo patrón en 2001 y 2002 con ciertas variaciones en las densidades de las especies, de la temperatura y salinidad, lo que permitirá a la posteridad la caracterización de patrones bien establecidos de dicha área. Evans y Granger (1990) y  Contreras Espinoza  (1993) señalan que la temperatura está relacionada íntimamente con la reproducción ya que en aguas cálidas, se llegan a presentar varias generaciones al año, en tanto que en las zonas templadas o frías las condiciones limitan a una o dos generaciones. Por otro lado, Navarro-Rodríguez et al., (2004) indican que en El Salado, sistema cercano al área de estudio (15 km), presenta un patrón de distribución diferente puesto que se registraron las mayores densidades (12,000 org./1000 m-3) con el incremento de la temperatura, pero con una disminución notoria en la salinidad (20-30 °C y 6,5 ups); señalan que las variaciones de la abundancia espacio-temporal estuvieron influenciadas tanto por las condiciones de marea como por la época del año, debido a que las mayores concentraciones de larvas fueron asociadas a la presencia de lluvias durante en el verano, en tanto que, en primavera se presentaron los registros más bajos durante todo el periodo de estudio asociado con la época de sequía. Es posible que estas diferencias se deban a las características naturales propias de cada sistema y las continuas modificaciones a las que se encuentran sujetos.  Algunos autores mencionan, que los sistemas estuarinos se encuentran fuertemente influenciados por las mareas y las fluctuaciones en el aporte de agua dulce de los sistemas hídricos adyacentes, además de los vientos y la geomorfología, ocasionando de alguna manera la sustitución de unas especies por otras, incluso dentro de una misma familia, hasta terminar con el remplazo de familias completas, o bien, incrementando la distribución de algunas otras especies (Torres-Orozco,1994; Sánchez Ramírez 1997; Contreras Espinoza, 2003; Kenneth y Fahay 2010). Aunado a esto, cabe hacer mención que la especie E. mordax fue muy abundante en otoño y primavera lo que se asocia a su ciclo reproductivo puesto que desova entre diciembre y junio por lo que sus densidades larvales se ven reflejadas  en primavera y principalmente en otoño, aun cuando es una especie que es capturada generalmente muy al norte del Pacífico mexicano incluyendo el Golfo de California (Whitehead et al., 1995),  no es común encontrarla en estas latitudes; sin embargo, no significa que no pueda presentarse en estas áreas, puesto que la Corriente Costera de Costa Rica durante junio y julio sigue la costa de Centro América y México (donde se le conoce como Corriente Mexicana), hasta alcanzar la zona de transición del Pacífico Norte, frente al extremo de la península de Baja California y, que conjuntamente giran para realimentar el sistema ecuatorial (Instituto Oceanográfico del Pacífico, 1995; Badan, 1997;  Navarro-Rodríguez et al.,  2001), lo que favorece de alguna manera la distribución de dicha especies hasta esta área de estudio. De acuerdo a los intervalos de confianza del análisis de diversidad temporal se observó que primavera es la estación del año con el mayor nivel de diversidad (1.3 bits/ind.) (Figura 6), con un valor de la riqueza de 4.8 (no. de especies) (Figura 7) y una equidad de 0.8 (Figura 8). En tanto que en invierno se refleja la diversidad más baja del período de estudio (0.3 bits/ind.) (Figura 6), con una riqueza específica menor (1 no. de especies). Sin embargo, Navarro-Rodríguez et al., (2004) mencionan que en el estero El Salado lado, la diversidad fue mínima ya que se obtuvieron registros del 97.8 % para Dormitator latifrons, en tanto que el 2.2 % restante correspondió al resto de al especies. Es posible que estas diferencias se deban a toda una serie de factores ambientales propios de los sistemas, y principalmente, a las actividades antropogénicas sobre todo a las de infraestructura turística, lo que desfavorece en gran medida a las diferentes especies que utilizan estos cuerpos de agua ya sea como área de refugio, de alimentación, reproducción y crianza, lo que concuerda con diversos autores (Saldierna y Vera, 1993; Núñez- Moreno, 1996; Funes et al., 1998; Miller y Kendal, 2009; Cowan et al., 2013),   puesto que las  especies marinas fueron las que dominaron en el sistema El Quelele (menos impactado) en tanto que en el estero El Salado (altamente impactado), la especie dominante es principalmente de hábitos estuarinos, a pesar de que ambos cuerpos de agua se localizan en la misma zona climática (semicálida subhúmeda y se encuentran bajo el mismo régimen de mareas. Respecto a los índices de diversidad temporal Navarro-Rodríguez et al., (2006), señalan que en otoño se registró el valor máximo (1.21 bits/ind.) mientras que, en invierno se reportaron los valores mínimos (1.10 bits/ind.), en tanto, en el presente trabajo, los intervalos de confianza del análisis de diversidad temporal mostraron que en primavera  se presenta  el mayor nivel de diversidad (1.3 bits/ind.), en tanto que en invierno se refleja la diversidad más baja del período de estudio (0.3 bits/ind.). Lo que sugiere que esta variación en la diversidad pueda estar asociada a los hábitos y ciclos reproductivos de las especies, así como a las altas concentraciones de salinidad puesto que son especies principalmente marinas.

Referente al análisis de cluster aplicado a las variables especies-sitios de muestreo mostró una mayor similitud o afinidad entre los sitios 3 y 5 (72.2 %), mientras que en los sitios 1 y 2 la similitud fue menor (67.4 %), el sitio 4 es el que menor afinidad mostró con el resto de los sitios (Figure  9a). Referente al análisis entre las  variables especies-estación del año  se observa un grupo con una afinidad del 60.52 % entre primavera y otoño en tanto verano e invierno no presentaron  ningún tipo de afinidad con el resto de las estaciones (Figura 9b).  Cabe mencionar que aún cuando las afinidades entre los sitios de muestreo 3-5 y 1-2 compartieron generalmente las mismas especies, el sitio 4 aún cuando comparten algunas especies,  fue totalmente diferente al resto de los sitios este presenta una menor densidad, diversidad y equidad, posiblemente atribuido a su  ubicación y a que se encuentra bajo un mayor efecto de impacto antropogénico que el resto de los sitios. Sin embargo, Navarro-Rodríguez et al., (2006) señalan que en este mismo sistema (2001) los resultados del análisis de similitud aplicado en las variables especies-sitios de muestreo mostró una mayor afinidad entre los sitios 2 y 3 mientras que en los sitios 1 y 5, la afinidad fue menor; el sitio 4 es el que menos afinidad presentó con el resto de los sitios, situación completamente diferente en el 2002 ya que los sitios afines cambian completamente a excepción del sitio 4, que en ambos periodos presenta afinidad diferente con los otros cuatro puntos. Estas diferencias se atribuyen principalmente a las modificaciones que ha venido presentando este cuerpo de agua en los diferentes periodos de estudio, así como a los diferentes tiempos de la toma de las muestras, aún cuando se establecen en las mismas épocas estacionales.   

Asimismo fue observado entre estaciones del año afinidades entre primavera y otoño cuya afinidad pudiese deberse a la presencia de especies marinas principalmente de E. mordax, Eucinostomus sp., A. mazatlanus y T. rhodopus.  y que utilizan este cuerpo de agua en alguna etapa de su desarrollo o bien debido a la tolerancia de las variaciones en la salinidad e incluso a su ciclo reproductivo (Fischer et al., 1995). Por su parte Navarro-Rodríguez et al., (2006) señalan que este mismo sistema (2001)se observaron dos grupos con afinidades entre las estaciones frías (otoño-invierno) y las estaciones cálidas (primavera-verano) indicando que las especies fueron predominante mente marinas con una mayor densidad de E. mordax y  T. rhodopus, lo que puede indicar que la presencia de especies marinas es asociada al transporte por las corrientes, patrones de circulación y efectos de mareas. Esto mismo ha sido indicado por Álvarez-Cadena et al., (1984), Day y Yáñez Arancibia (1985) y Contreras Espinoza (1993).

Las variaciones espacio-temporales en la temperatura, salinidad y abundancia del ictioplancton indican que la abundancia larval se encuentra en estrecha relación con el patrón estacional influyendo a su vez en el ciclo reproductivo de estas especies, en otoño y primavera con registros de temperatura promedio moderadas y bajas salinidades asociada a las épocas frías y cálidas. Por otro lado los picos de la época reproductiva de la mayoría de las especies pero principalmente el de E. mordax  se presentan particularmente en una época cálida con abundantes lluvias lo que favorece al aumento de la productividad del sistema incrementando los nutrientes y la materia orgánica en suspensión, por lo que hay una mayor disponibilidad de alimento para los estadios  larvales y juveniles de peces. Y otro más en la época de estiaje, asociado a las altas concentraciones de salinidades  puesto que la mayoría de las especies son de hábitos marinos.

Average temperature and salinity variation per season of the year registered the highest value of temperature (34.6 °C) in spring, and the lowest in summer (23.4 °C); in contrast, the highest value for salinity was present in summer (34.6 psu) and the lowest in winter (24.2 ups); it should be noted that during the fall no registers were obtained (Figure 2a). Variation of such parameters per sampling site presented homogeneity among themselves, even if the highest value (30 °C) of average temperature was registered in the sites 4 and 5, and the lowest in the site 2 (28.3 °C); however, temperature fluctuations among sites were lower than a degree, same situation was observed in the variation of salinity among sites (variation of 0.5 psu), even when the highest value (30.8 psu) was presented in site 2, while the lowest value (29.5 psu) was registered in site 3 (Figure 2b). The ichthyoplanktonic group corresponding to the period of study (2002) was represented by 11 families, with 13 genres and new species. Their densities were represented mainly by Engraulis mordax (215.28 org/1000m³), with relative abundance of 40.2 %, followed by Eucinostomus sp. (133.78 org/1000 m³),with 22.68 %, while Dormitator latifrons 76.88 org/1000 m3), presented 13.14 % and finally, Ophistonema libertate and Achirus mazatlanus with low densities (49.59 org/1000m³  respectively), presented relative abundance of 7.47 and 7.21 % respectively. Densities of the rest of the species fluctuated from 33.83 to 1.53 org/1000m³ and their relative abundances were from 5.66 to 0.25 % (Figure 3, Table 1).

Regarding time variation, fall was the season of the year where highest larvae density was registered (269.1 org./1000m³), followed by spring and summer (175.3 and 127.63 org./1000m³ respectively), while a minimum density was reached in winter (29.21 org./1000m³) (Figure 4).

In respect of spatial variation, site 3 was the highest larvae density point ((176.83 org./1000 m³), followed by sites 5 and 2 (155.31 and 141.47 org./1000 m³ respectively), while sites 1 and 4 reported the lowest densities from the total of the study period (72.27 and 55.35 org./1000 m³ respectively) (Figure 5). The highest larvae density was registered by E. mordax and A mazatlanus in site 3, with absolut dominium of E. mordax in sites 5 and 2, not so for the sites 1 and 4, where the highest records were D. latifrons and Eucinostomus sp. Respectively. Spatial and temporal variations in temperature, salinity and abundance of ichthyoplankton indicate that larvae abundance is deeply related with the seasonal pattern, influencing the reproductive cycle of these species, since the highest abundance was presented in fall and spring (269.1 and 175.3 org./1000 m³ respectively), with a wide interval of temperature and salinity (29.1 °C to 34.6 °C and up to 32.8 psu) with higher values for both variables in spring, not so for summer, which suggests the influence of the rainy season, abundant in the fall and null in spring. When highest larvae densities are present with wide intervals of temperature (29 °C to 34.6 °C) and salinity (34.6 psu), Margalef (1969), Livingston (1984) and Contreras Espinoza (1993) state that in coast lagoons, zooplankton is associated with the variability of the system in terms of salinity and temperature. On the other hand, Navarro-Rodríguez et al., (2006), state that in this same system (2001), spatial – temporal variations were influenced by temperature fluctuations (29.8-31.6 °C) and salinity (11-34.8 psu), presenting high concentrations in fall and minor in spring, associated to cold and warm times, respectively. Furthermore, the indicate that the peaks of the reproductive season in most of the species (nine in spring and seven in the fall), mainly the E. mordax (the most abundant species), present mainly in a warm season with heavy rain. And one more in the dry season, associated probably with the high salinity concentrations (32.8 psu) since the great majority of species are of marine habits, which can also be associated to the transport of ocean currents, local circulation patterns as well as tide effects. It suggests that the same pattern is presented in 2001 and 2002 with certain variations in species densities, temperature and salinity, which will allow in a short future the characterization of well stablished patterns in the area. Evans y Granger (1990) and Contreras Espinoza (1993) state that temperature is deeply related with reproduction, since in warm waters, several generations are presented per year, while in warm or cold zones, conditions limit it to one or two generations. On the other hand, Navarro-Rodríguez et al., (2004) indicate that in El Salado, close system to the area of study (15 km), there is a different distribution pattern, since the highest densities (12,000 org./1000 m-3) were registered, with an increase in temperature but a decrease in salinity (20-30 °C and 6,5 psu); they state that variations in spatial – temporal abundance were influenced for both tide conditions and season of the year, since major concentration of larvae were associated to rain presence during the summer, while in spring, the lowest records during the study period associated with dry season were presented. It is possible that these differences are due to natural characteristics of each system and continuous modifications to which they are subject. Some authors mention that estuary systems are deeply influenced by tides and fluctuation in the freshwater input from the adjacent hydraulic systems, in addition to the winds and the geomorphology, resulting somehow in the substitution of some species for others, even within the same family, until they reach a replacement of complete families, or, increasing the distribution of some other species especies (Torres-Orozco,1994; Sánchez Ramírez 1997; Contreras Espinoza, 2003; Kenneth and Fahay, 2010). In addition, the species E. mordax was very abundant in fall and spring, which is associated to its reproductive cycle, because it spawns between December and June, therefore its larvae densities are reflected in spring and mainly in the fall, even when it is a species captured generally in the north of the Mexican Pacific, including the Gulf of California (Whitehead et al., 1995). It is not common to find it in these latitudes; however, it does not mean that in can be present in these areas, because the Costa Rica Coastal Current during June and July follows the coast of Central America and Mexico (where it is known as the Mexican Current), until it reaches the transition zone of the North Pacific, in front of the tip of the Baja California Peninsula and which together rotate in order to replenish the Equatorial System (Instituto Oceanográfico del Pacífico, 1995; Badan, 1997;  Navarro-Rodríguez et al., 2001), which favors the distribution of such species to this area of study. According to confidence intervals of the temporal diversity analysisis, it was observed that spring is the first season of the year with the highest diversity level (1.3 bits/ind.) (Figure 6), with richness value of 4.8 (no. of species) and equity of 0.8 (Figure 8). While the lowest diversity of the period of study is reflected in winter (0.3 bits/ind.) (Figure  6), with less specific richness (1 no. of species). However, Navarro-Rodríguez et al., (2004) mention that in the stuary El salado, diversity was mínimum, since records of 97.8 % were obtained for Dormitator latifrons, while the rest 2.2 % corresponded to the remaining species. It is possible that these differences are due to a series of environmental factors from the own systems, and mainly, to anthropogenic activities on the touristic infrastructure, which works against the different species that use these bodies of water, either as shelter, feeding, reproduction or breeding area, which agrees with several authors (Saldierna and Vera, 1993; Núñez-Moreno, 1996; Funes et al., 1998; Miller and Kendal, 2009; Cowan et al., 2013), since marine species dominated the El Quelele system (with less impact), while the estuary El Salado (with higher impact), the dominant species is mainly from of estuary habitats, despite the fact that both bodies of water are located in the same climate zone (semi-warm  semi-humid and under the same regime of tides). Regarding the temporal diversity indexes, Navarro-Rodríguez et al., (2006), state that during the fall the maximum value was registered (1.21 bits/ind.) while in winter minimum values were reported (1.10 bits/ind.). In this study, confidence intervals of the temporal diversity analysis showed that in spring the highest level of diversity is present (1.3 bits/ind.), while in winter the lowest diversity of the study period is reflected (0.3 bits/ind.). This suggests that the variation in diversity can be associated to the habits and reproductive cycles of species, as well as high salinity concentrations, being species primarily marine.

In respect of the cluster analysis applied to the variables species – sampling sites, a higher similarity or affinity was shown between sites 3 and 5 (72.2 %), while in sites 1 and 2 similarity was less (67.4 %), site 4 is the one that less affinity showed with the rest of the sites (Figure 9a). Regarding the analysis between the species – season of the year variables, a group with an affinity of 60.52 % is observed between spring and fall, while they did not present any type of affinity in summer and winter with the rest of the seasons (Figure 9b). It should be noted that even when affinities between sampling sites 3-5 and 1-2 generally shared the same species, site 4 was completely different from the rest of the sites, it present lower density, diversity and equity, probably due to its location and the fact that it is found under a higher effect of anthropogenic impact than the rest of the sites. However, Navarro-Rodríguez et al., (2006) state that in this same system (2001), the results of the similarity analysis applied in the variables species – sampling sites, showed higher affinity between sites 2 and 3, while in sites 1 and 5 affinity was lower; site 4 is the one that less affinity presented with the rest of the sites, situation completely different from 2002, since related sites change completely excepting site 4, which in both periods presented different affinity with the other four points. These differences are mainly attributed to the modifications that this body of water has presented in the different periods of study, as well as the different times in the taking of samples, even when they are established in the same seasonal times.

In addition, it was observed that between seasons of the year, affinities between spring and fall, which can be due to the presence of marine species, mainly E. mordax, Eucinostomus sp., A. mazatlanus and T. rhodopus, that use this body of water at some stage of their development or due to the tolerance of variations in salinity, or even due to their reproductive cycle (Fischer et al., 1995). Navarro-Rodríguez et al., (2006) state that in this system (2001) two groups with affinities between cold seasons (fall – winter) and warm seasons (spring–summer) were observed, indicating that the species were mostly marine, with a higher density of E. mordax and T. rhodopus, indicating that the presence of marine species is associated to transportation of currents, circulation patterns and tide effects. The same has been indicated by Alvarez-Cadena et al., (1984), Day and Yáñez Arancibia (1985) and Contreras Espinoza (1993).

Spatial–temporal variations in temperature, salinity and abundance of ichthyoplankton indicate that larvae abundance is strongly related with the seasonal pattern, influencing in the reproductive cycle of these species, in fall and spring, with moderated and low salinity records of average temperature, associated to cold and warm timings. On the other hand, peaks of reproductive season of most species, but mainly of E. mordax, are present in a particularly warm season with heavy rain, which favors the increase of productivity of the system, increasing nutrients and organic matter in suspension, so there is a higher availability of food for fish larvae and juvenile stays. And one more in dry weather timing, associated to high saline concentrations, since most of the species belong to marine habits.