Occurrence and abundance of young mullet Mugil liza (Teleostei: Mugilidae) in the surf zone along the southern coast of Brazil

INTRODUCTION

Estuaries and coastal lagoons are considered highly productive areas and are important refuge and nursery zones for several fish species with high economic and/or environmental value (Boehlert and Mundy 1988Boehlert G.W., Mundy B.C. 1988. Roles of behavioral and physical factors in larval and juvenile fish recruitment to estuarine nursery areas. In: American Fisheries Society Symposium 3: 1-67.). Migratory marine species spawn in offshore waters, regularly enter estuaries in substantial numbers, particularly as juveniles, and use the nearshore marine waters as an alternative habitat (Boehlert and Mundy 1988Boehlert G.W., Mundy B.C. 1988. Roles of behavioral and physical factors in larval and juvenile fish recruitment to estuarine nursery areas. In: American Fisheries Society Symposium 3: 1-67.). In southern Brazil, juveniles of many commercially and recreationally important fish species (Micropogonias furnieri [Desmarest, 1823], Brevoortia pectinate [Jenyns, 1842], Pogonias courbina [Linnaeus, 1766] and Mugil liza Valenciennes, 1836) use the estuarine areas as nurseries (Seeliger et al. 1998Seeliger U., Odebrecht C., Castello J.P. 1998. Os Ecossistemas Costeiro e Marinho do Extremo Sul do Brasil. Ed. Ecoscientia, Rio Grande, 337 pp., Mai et al. 2019Mai A.C.G., Albuquerque C.Q., Lemos V.M., et al. 2019. Coastal zone use and migratory behavior of the southern population of Mugil liza in Brazil. J. Fish. Biol. https://doi.org/10.1111/jfb.14120 ).

The Lebranche mullet, Mugil liza (Teleostei, Mugilidae), occurs from the Gulf of Mexico to the coast of Argentina. Using the mitochondrial control region, Heras et al. (2016)Heras S., Maltagliati F., Fernández M.V., Roldán M.I. 2016. Shaken not stirred: A molecular contribution to the systematics of genus Mugil (Teleostei, Mugilidae). Integr. Zool. 11: 263-281. https://doi.org/10.1111/1749-4877.12173 described two highly divergent M. liza populations: a northern population, which is distributed from Cuba to the state of Rio de Janeiro, Brazil, and a southern population, which is distributed from Rio Grande do Sul, Brazil, to Argentina. Mai et al. (2014)Mai A.C.G., Minõ C.I., Marins L.F.F., et al. 2014. Microsatellite variation and genetic structuring in Mugil liza (Teleostei: Mugilidae) populations from Argentina and Brazil. Est. Coast. Shelf. Sci. 149: 80-86. https://doi.org/10.1016/j.ecss.2014.07.013 and Lemos et al. (2017)Lemos V.M., Monteiro-Neto C., Cabral H., Vieira J.P. 2017. Stock identification of tainha (Mugil liza) by analyzing stable carbon and oxygen isotopes in otoliths. Fish. Bull. 115: 201-206. https://doi.org/10.7755/FB.115.2.7 showed that the southern population encompasses the Rio Grande do Sul state and actually ranges from the state of São Paulo in southeastern Brazil (23°S) to Argentina.

Mugil liza is commercially exploited in Argentina by a small-scale fishery that operates mainly in Samborombón Bay. Commercial catches were between 5.4 and 78.8 t from 2000 to 2010, with a maximum capture of 194.0 t in 2004 (Castellini et al. 2017Castellini D.L., Brown D., Lajud N.A., et al. 2017. Juveniles recruitment and daily growth of the southern stock of Mugil liza (Actinopterygii; Fam. Mugilidae): new evidence for the current life-history model. J. Mar. Biol. Assoc. U.K. 99: 204-215. https://doi.org/10.1017/S0025315417001904 ). In the southern states of Brazil (Rio Grande do Sul, Santa Catarina and Paraná) and São Paulo, the M. liza southern population fishery is a very important resource (Lemos et al. 2014Lemos V.M., Varela Jr. A.S., Schwingel P.R., et al. 2014. Migration and reproductive biology of Mugil liza (Teleostei : Mugilidae) in south Brazil. J. Fish. Biol. https://doi.org/10.1111/jfb.12452 , 2016Lemos V.M., Troca D.F.Á., Castello J.P., Vieira J.P. 2016. Tracking the southern Brazilian schools of Mugil liza during reproductive migration using VMS of purse seiners. Lat. Am. J. Aquat. Res. 44: 238-246. https://doi.org/10.3856/vol44-issue2-fulltext-5 , de Abreu-Mota et al. 2018de Abreu-Mota M.A., Medeiros R.P., Noernberg M.A. 2018. Resilience thinking applied to fisheries management: perspectives for the mullet fishery in Southern-Southeastern Brazil. Reg. Environ. Change 18: 2047-2058. https://doi.org/10.1007/s10113-018-1323-9 ), with a maximum reported capture of 13600 t in 2007 (de Abreu-Mota et al. 2018de Abreu-Mota M.A., Medeiros R.P., Noernberg M.A. 2018. Resilience thinking applied to fisheries management: perspectives for the mullet fishery in Southern-Southeastern Brazil. Reg. Environ. Change 18: 2047-2058. https://doi.org/10.1007/s10113-018-1323-9 ). The Patos Lagoon estuary, located in the state of Rio Grande do Sul, is considered the main M. liza nursery area in Brazil (Herbst and Hanazaki 2014Herbst D.F., Hanazaki N. 2014. Local ecological knowledge of fishers about the life cycle and temporal patterns in the migration of mullet (Mugil liza) in Southern Brazil. Neotrop. Ichthyol. 12: 879-890. https://doi.org/10.1590/1982-0224-20130156 , de Abreu-Mota et al. 2018de Abreu-Mota M.A., Medeiros R.P., Noernberg M.A. 2018. Resilience thinking applied to fisheries management: perspectives for the mullet fishery in Southern-Southeastern Brazil. Reg. Environ. Change 18: 2047-2058. https://doi.org/10.1007/s10113-018-1323-9 ).

The southern population of M. liza performs an annual reproductive migration from Argentina, Uruguay and estuaries of southern Brazil, especially the Patos Lagoon estuary, to offshore spawning areas located on the Brazilian coast between northern Santa Catarina and Paraná (28°-26°S) (Lemos et al. 2014Lemos V.M., Varela Jr. A.S., Schwingel P.R., et al. 2014. Migration and reproductive biology of Mugil liza (Teleostei : Mugilidae) in south Brazil. J. Fish. Biol. https://doi.org/10.1111/jfb.12452 , 2016Lemos V.M., Troca D.F.Á., Castello J.P., Vieira J.P. 2016. Tracking the southern Brazilian schools of Mugil liza during reproductive migration using VMS of purse seiners. Lat. Am. J. Aquat. Res. 44: 238-246. https://doi.org/10.3856/vol44-issue2-fulltext-5 ). Mugil liza is a single spawner and spawning occurs between May and July at temperatures of 19°C to 21°C in offshore waters (Lemos et al. 2014Lemos V.M., Varela Jr. A.S., Schwingel P.R., et al. 2014. Migration and reproductive biology of Mugil liza (Teleostei : Mugilidae) in south Brazil. J. Fish. Biol. https://doi.org/10.1111/jfb.12452 ). According to the life cycle model proposed by Vieira (1991)Vieira J.P. 1991. Juvenile mullets (Pisces: Mugilidae) in the Estuary of Lagoa dos Patos, RS, Brazil. Copeia 409-418. https://doi.org/10.2307/1446590 , M. liza larvae and early juveniles are neustonic and are carried, mainly by wind-driven currents, from offshore waters to the coast, where they concentrate in surf zones prior to recruitment to estuarine environments. This use of the coastal zone and migratory behaviour characterize M. liza as a marine estuarine-dependent species (Mai et al. 2019Mai A.C.G., Albuquerque C.Q., Lemos V.M., et al. 2019. Coastal zone use and migratory behavior of the southern population of Mugil liza in Brazil. J. Fish. Biol. https://doi.org/10.1111/jfb.14120 ).

Martin and Drewry (1978)Martin F.D., Drewry G.E. 1978. Development of fishes of the Mid-Atlantic Bight. U.S. Dept. Inter., Fish and Wildl. Serv. recognized a pre-juvenile stage for mullets starting at 11 mm total length (TL), and the transformation from the larval to the pre-juvenile stage ends at approximately 13 mm TL and involves a change from predaceous to iliophagous (detritivorous) feeding habits (Martin and Drewry 1978Martin F.D., Drewry G.E. 1978. Development of fishes of the Mid-Atlantic Bight. U.S. Dept. Inter., Fish and Wildl. Serv., Vieira 1991Vieira J.P. 1991. Juvenile mullets (Pisces: Mugilidae) in the Estuary of Lagoa dos Patos, RS, Brazil. Copeia 409-418. https://doi.org/10.2307/1446590 ). Young mullets, which first appear in small schools along the coasts and in the estuaries, measure 18 to 28 mm in length (Vieira 1991Vieira J.P. 1991. Juvenile mullets (Pisces: Mugilidae) in the Estuary of Lagoa dos Patos, RS, Brazil. Copeia 409-418. https://doi.org/10.2307/1446590 , González-Castro and Minos 2016González-Castro M., Minos G. 2016. Sexuality and reproduction of Mugilidae. In: Crosetti D., Blaber S. (eds) Biology, ecology and culture of mullets (Mugilidae). CRC Press. Boca Raton. pp. 227-263.). In the coastal area of the state of Rio Grande do Sul, young mullets are frequently present and abundant in the surf zone throughout the year (Ramos and Vieira 2001Ramos L.A., Vieira J.P. 2001. Composição específica e abundância de peixes de zonas rasas dos cinco estuários do Rio Grande do Sul, Brasil. Bol. Inst. Pesca. 27: 109-121., Rodrigues et al. 2015Rodrigues F.L., Cabral H., Vieira J.P. 2015. Assessing surf zone fish assemblage variability in southern Brazil. Mar. Freshwater. Res. 66: 106-119. https://doi.org/10.1071/MF13210 ), which is also the case on the coast of Argentina (Acha 1990Acha E.M. 1990. Estudio anatômico-ecologico de la lisa (Mugil liza) durante su primer año de vida. Frente. Marit. 7: 37-43., Castellini et al. 2017Castellini D.L., Brown D., Lajud N.A., et al. 2017. Juveniles recruitment and daily growth of the southern stock of Mugil liza (Actinopterygii; Fam. Mugilidae): new evidence for the current life-history model. J. Mar. Biol. Assoc. U.K. 99: 204-215. https://doi.org/10.1017/S0025315417001904 ). However, in the northern area of the southern population (state of São Paulo) (Mai et al. 2014Mai A.C.G., Minõ C.I., Marins L.F.F., et al. 2014. Microsatellite variation and genetic structuring in Mugil liza (Teleostei: Mugilidae) populations from Argentina and Brazil. Est. Coast. Shelf. Sci. 149: 80-86. https://doi.org/10.1016/j.ecss.2014.07.013 ), the presence of young mullets in the surf zone occurs shortly after the reproductive period, and individuals are less abundant or absent in the surf zone during the remaining months of the year (Favero and Dias 2015Favero J.M. del, Dias J.F. 2015. Juvenile fish use of the shallow zone of beaches of the Cananéia - Iguape coastal system, southeastern Brazil. Braz. J. Oceanogr. 63: 103-114. https://doi.org/10.1590/S1679-87592015081806302 ). It seems that there is a clinal decrease in surf zone occurrence from the northern to the southern limits of the southern population.

Surf zones of sandy beaches are important recruitment and nursery zones and migratory paths to other nearshore habitats (Strydom 2003Strydom N.A. 2003. Occurrence of larval and early juvenile fishes in the surf zone adjacent to two intermittently open estuaries, South Africa. Environ. Biol. Fish. 66: 349-359. https://doi.org/10.1023/A:1023949607821 , McLachlan and Brown 2006McLachlan A., Brown A. 2006. The Ecology of Sandy Shores. Academic Press, 392 pp.). Larvae and juvenile of estuary-dependent marine fishes move from offshore spawning grounds (initially passively but later by active migration) to surf zones (Boehlert and Mundy 1988Boehlert G.W., Mundy B.C. 1988. Roles of behavioral and physical factors in larval and juvenile fish recruitment to estuarine nursery areas. In: American Fisheries Society Symposium 3: 1-67.). Despite the available information on recruitment of young fishes into estuaries worldwide, there is a paucity of information describing the factors driving the recruitment process (Strydom 2003Strydom N.A. 2003. Occurrence of larval and early juvenile fishes in the surf zone adjacent to two intermittently open estuaries, South Africa. Environ. Biol. Fish. 66: 349-359. https://doi.org/10.1023/A:1023949607821 ). The aim of the present study was to compare the abundance of young M. liza (≤30 mm TL) in the adjacent surf zones of six estuaries between 28.5°S and 33.7°S in southern Brazil and to determine which environmental variables are related to their occurrence and abundance.

MATERIALS AND METHODS

Study area

 ⌅

The study area encompassed the coastal surf zone of southern Brazil from Laguna, Santa Catarina state (29°S, 48°W), to Chuí, Rio Grande do Sul state (34°S, 53°W) (Fig. 1). These coastal areas are characterized by sandy beaches exposed to wave action with a wide dissipative to intermediate surf zone (Calliari and Klein 1995Calliari L.J., Klein A.H.F. 1995. Variabilidade das praias oceânicas entre Cassino e Chuí, RS. IO/USP. 11: 55-63.). Tidal oscillations along the southern Brazilian coast are defined by their low amplitude, with an average tidal amplitude of 0.47 m and minimal influence on estuarine circulation (Seeliger et al. 1998Seeliger U., Odebrecht C., Castello J.P. 1998. Os Ecossistemas Costeiro e Marinho do Extremo Sul do Brasil. Ed. Ecoscientia, Rio Grande, 337 pp.).

Fig. 1.  - Map of the six sampling sites: Site 1, Chuí; Site 2, Rio Grande; Site 3, Mostardas; Site 4, Tramandaí; Site 5, Passo de Torres; and Site 6, Laguna. Numbers (1 to 12) represent sampling points where beach seine samples were obtained.

The regional climate is temperate with well-defined seasons and is under the control of the South Atlantic anticyclone high-pressure centre. The latitudinal movement of the anticyclone centre and the passage of polar frontal systems influence and modify the seasonal cycle of the climate (Seeliger et al. 1998Seeliger U., Odebrecht C., Castello J.P. 1998. Os Ecossistemas Costeiro e Marinho do Extremo Sul do Brasil. Ed. Ecoscientia, Rio Grande, 337 pp.). These systems are responsible for dominant winds from the northeast quadrant (NE, mean speed 5 m s^-1) throughout the year and winds from the southwest quadrant (SW, mean speed 8 m s^-1) during the passage of polar frontal systems, which are common during the winter. The mean annual rainfall (1200-1500 mm) is related to the pattern and frequency of these systems and to the effects of the El Niño-Southern Oscillation cycle (Seeliger et al. 1998Seeliger U., Odebrecht C., Castello J.P. 1998. Os Ecossistemas Costeiro e Marinho do Extremo Sul do Brasil. Ed. Ecoscientia, Rio Grande, 337 pp.).

The complex hydrographic system of the coastal plain in southern Brazil includes several lagoons and estuaries. Our study area encompassed six estuarine systems: 1) Chuí Creek (33.73°S), 2) the Patos Lagoon estuary (32.20°S), 3) the Peixe Lagoon (31.15°S), 4) the Tramandaí-Armazém Lagoon complex (29.98°S), 5) the Mampituba River (29.32°S), and 6) the Imaruí-Mirim-Santo Antônio Lagoon complex (28.49°S) (Fig. 1, Table 1).

Table 1.  - Summary of environmental variables used in the Mugil liza occurrence and abundance models.

Site	Lat.-Long.	Estuarine area (km2 )	Sampling point number	Tide	Distance from estuary mouth (km)	Salinity		Temperature (°C)		Transparency (cm)
						min - max	mean ± sd	min - max	mean ± sd	min -max	mean ± sd
1. Chuí	33.73°S, 53.36°W	0.3a	1	0.6e	0.3	24.0-30.3	28.1±1.7	13.2-28.0	19.8±4.6	10-50	33.6±13
			2		1.4	26.4-30.1	28.5±1.2	14.9-26.0	20.0±4.0	10-100	39.0±23
2. Rio Grande	32.20°S, 52.17°W	971.0a	3	0.47f	8.8	15.8-32.6	28.5±4.5	12.0-28.0	20.3±4.4	10-60	30.0±16
			4		0.7	17.3-33.0	28.1±4.6	11.4-28.5	20.1±4.5	10-100	40.3±23
3. Mostardas	31.15°S, 50.80°W	45.0b	5	0.45g	13	29.0-36.4	32.0±2.6	14.9-26.0	20.2±3.1	15-60	31.0±17
			6		13.9	29.1-37.5	32.4±2.8	15.2-25.7	20.3±3.4	10-70	34.0±18
4. Tramandaí	29.98°S, 50.12°W	30.0a	7	0.3h	3.6	25.2-36.7	30.9±3.0	15.4-24.8	20.6±2.7	0-100	44.6±30
			8		0.4	22.6-36.9	30.8±3.8	15.4-25.0	21.0±2.6	15-70	45.8±14
5. P. Torres	29.32°S, 49.71°W	0.5c	9	0.5i	0.2	6.4-37.8	23.4±9.3	16.6-26.0	20.7±2.9	20-60	40.3±13
			10		0.8	15.5-38.2	29.8±5.0	15.0-23.5	20.1±3.0	0-45	29.1±14
6. Laguna	28.49°S, 48.76°W	184.0d	11	0.6j	0.3	28.0-36.8	32.5±2.3	16.1-23.0	19.1±2.3	35-100	71.9±24
			12		2.2	25.2-37.4	31.9±3.2	16.3-23.0	19.1±2.3	30-100	64.6±25

Estuarine area and mouth width obtained from the literature were a (Ramos and Vieira 2001Ramos L.A., Vieira J.P. 2001. Composição específica e abundância de peixes de zonas rasas dos cinco estuários do Rio Grande do Sul, Brasil. Bol. Inst. Pesca. 27: 109-121.), b (Loebmann and Vieira 2005Loebmann D., Vieira J.P. 2005. Distribuição espacial e abundância das assembléias de peixes no Parque Nacional da Lagoa do Peixe, Rio Grande do Sul, Brasil. Ver. Bras. Zool. 22: 667-675. https://doi.org/10.1590/S0101-81752005000300023 ) and d (Fonseca and Netto 2006); the rest were calculated with the ArcMap tool (c). Tides were obtained from the literature: e (Hartmann and Pereira 2001Hartmann C., Pereira L.E. 2001. Arroio Chuí: influencia antrópica na manutenção de um ambiente estuarino. Uberlândia, Sociedade e Natureza 13: 27-36.), f (Seeliger et al. 1998Seeliger U., Odebrecht C., Castello J.P. 1998. Os Ecossistemas Costeiro e Marinho do Extremo Sul do Brasil. Ed. Ecoscientia, Rio Grande, 337 pp.), g (Schwarzbold and Schäfer 1984Schwarzbold A., Schäfer A. 1984. Gênese e Morfologia das Lagoas Costeiras do Rio Grande do Sul-Brasil. Amazoniana 9: 87-104.), h (Lira et al. 1976Lira L., Martins I., Araujo T.H., Dehnardt E.A. 1976. Nota prévia sobre o comportamento da cunha salina no Estuário de Tramandaí RS. Ciências Biológicas 3: 115-126.), i (D’Aquino et al. 2011D’Aquino C.D.A., Andrade Neto J.S.D., Barreto G.A.M., Schettini C.A. 2011. Caracterização oceanográfica e do transporte de sedimentos em suspensão no estuário do Rio Mampituba, SC. Rev. Bras. Geofis. 29: 217-230. https://doi.org/10.1590/S0102-261X2011000200001 ) and j (Pimenta 1958Pimenta J. 1958. A faixa costeira meridional de Santa Catarina. Boletim no. 176, Rio de Janeiro: Divisão de Geologia e Mineralogia do DNPM.).

1. Chuí Creek, Rio Grande do Sul (RS), Brazil, is located on the border between Brazil and Uruguay. The Chuí Creek drainage basin has an area of approximately 340 km² and an estuarine area of 0.5 km² (Ramos and Vieira 2001Ramos L.A., Vieira J.P. 2001. Composição específica e abundância de peixes de zonas rasas dos cinco estuários do Rio Grande do Sul, Brasil. Bol. Inst. Pesca. 27: 109-121.). Chuí Creek arises in an area of wetlands and runs almost parallel to the coast in a southerly direction for approximately 66 km until it reaches the Atlantic Ocean. Chuí Creek has suffered from anthropic modifications over time, including the change in its original course, the use of its fresh water and the construction of a pair of jetties (1500 m long) in the late 1970s to demarcate the border between Brazil and Uruguay. At low flow, the stream has a maximum width of 100 m and a depth ranging between 1.5 and 3 m. The amplitude of the astronomical tide is low (0.5 m) and the wind regime determines the saltwater flow into the system (Hartmann and Pereira 2001Hartmann C., Pereira L.E. 2001. Arroio Chuí: influencia antrópica na manutenção de um ambiente estuarino. Uberlândia, Sociedade e Natureza 13: 27-36.).

2. The Patos Lagoon estuary is located in Rio Grande, RS. The Patos Lagoon, with a total area of ~10000 km² (draining a large hydrological basin of 200000 km²), is one of the largest choked coastal lagoons in the world. It extends for 270 km along the coastline of southern Brazil. It was formed by multiple sand barrier complexes and can be divided into biological units. The estuarine region, with a total area of 900 km² (approximately 10% of the lagoon) is characterized by large shallow areas (<2 m deep). The estuary is connected to the South Atlantic via a 20 km channel (0.5-3000 m in width), and the mouth is stabilized by a pair of jetties. Owing to its morphological characteristics and the proximity of an amphidromic point, the tides in the region have low amplitude (mean 0.4 m). The wind patterns and precipitation in the watershed determine the salinity regime in the estuary (Seeliger et al. 1998Seeliger U., Odebrecht C., Castello J.P. 1998. Os Ecossistemas Costeiro e Marinho do Extremo Sul do Brasil. Ed. Ecoscientia, Rio Grande, 337 pp.).

3. The Peixe Lagoon is located in Mostardas, RS. It is situated in the narrow sandy strip between the Patos Lagoon and the Atlantic Ocean and was recognized as a National Park in 1986. It has an area of 34000 ha, including not only the Peixe Lagoon but also surrounding ecosystems such as beaches and small freshwater lakes. The Peixe Lagoon covers approximately 43 km² with a depth ranging from 0.3 to 2 m, and the lagoon has an intermittent opening to the ocean (Schwarzbold and Schäfer 1984Schwarzbold A., Schäfer A. 1984. Gênese e Morfologia das Lagoas Costeiras do Rio Grande do Sul-Brasil. Amazoniana 9: 87-104.).

4. The Tramandaí-Armazém Lagoon complex is located in Tramandaí, RS. The Tramandaí River Basin has an area of 3145 km² and the Tramandaí River is the primary effluent. From the estuary area (29°55’S to 30°00’S and 50°06’W to 50°11’W), the connection with the Atlantic Ocean is via a narrow, short channel (55 m) that is stabilized by a single jetty along the left margin. The amplitude of the astronomical tide in this region is low (0.3 m) (Lira et al. 1976Lira L., Martins I., Araujo T.H., Dehnardt E.A. 1976. Nota prévia sobre o comportamento da cunha salina no Estuário de Tramandaí RS. Ciências Biológicas 3: 115-126.).

5. The Mampituba River is located in Torres, RS, and Passo de Torres, Santa Catarina. The hydrographic basin of the Mampituba River (29°11’ to 29°26’S and 49°42’ to 50°12’W) has an area of 1224 km² that empties into the Atlantic Ocean. The Mampituba River is small (62 km long) and has an average flow of 19 m³ s^-1. The average depth is 4 m and the mouth is stabilized by a pair of jetties (north jetty, 100 m and south jetty, 260 m). The tidal pattern is microtidal and is primarily governed by the wind regime (D’Aquino et al. 2011D’Aquino C.D.A., Andrade Neto J.S.D., Barreto G.A.M., Schettini C.A. 2011. Caracterização oceanográfica e do transporte de sedimentos em suspensão no estuário do Rio Mampituba, SC. Rev. Bras. Geofis. 29: 217-230. https://doi.org/10.1590/S0102-261X2011000200001 ).

6. The Imaruí-Mirim-Santo Antônio Lagoon complex is located in Laguna, Santa Catarina. The lagoon system has an area of approximately 220 km² and has the Tubarão River as its main tributary. The communication between the lagoon system and the ocean occurs through a fixed bar in the Santo Antônio Lagoon that is approximately 200 m wide and 1200 m long. The tide is of low amplitude, and hydrodynamics are primarily governed by the wind (Pimenta 1958Pimenta J. 1958. A faixa costeira meridional de Santa Catarina. Boletim no. 176, Rio de Janeiro: Divisão de Geologia e Mineralogia do DNPM.).

RESULTS

Significant differences in salinity (ANOVA, F=22.69; p<0.05), temperature (ANOVA, F=10.74; p<0.05) and water transparency (ANOVA, F=19.06 p<0.05) were observed between sites and months over the study period (Fig. 2). The mean values of seawater surface temperature differed among the six sampling sites but showed seasonal fluctuations with a clear annual cycle: they declined steadily from April to August (monthly average 21°C-13.5°C), and thereafter the monthly average temperature increased to 18.5°C by September (Fig. 2A). Site 2 (Rio Grande) had higher temperatures in the warmer months and lower temperatures in the colder months. Site 6 (Laguna) had lower temperatures in the warmer months (Fig. 2A). Salinity showed no seasonal pattern (Fig. 2B). The lowest salinity values were recorded at Sites 5 (Passo de Torres) and 6 (Laguna) (the lowest salinity value was recorded in August at Site 5). The highest salinities were found from May to July (Fig. 2B). Significant differences were observed among sampling sites for water transparency, and sampling Site 6 (Laguna) showed the highest water transparency values in January, February, June, September and November (Fig. 2C). The nMDS plot was used to compare the similarity in the environmental factors of sampling sites. It is possible to observe three groups in the nMDS plot: Site 1 (Chuí), separated from the other sites; Sites 2 (Rio Grande) and 3 (Mostardas); and Sites 4 (Tramandaí), 5 (Passo de Torres) and 6 (Laguna) (Fig. 3).

A total of 6493 M. liza individuals were captured (Table 2). Significant differences in the number of individuals captured per beach seine sample (CPUE) were observed between the collection sites. Rio Grande (Site 2) had the highest abundance of M. liza (CPUE=26.4 individuals/haul; ANOVA, F=4.98, p<0.05), representing more than 52% of the total of individuals caught at the six collection sites. No significant differences in abundance were observed among the remaining sites; 1) Chuí, 3) Mostardas, 4) Tramandaí, 5) Passo de Torres, and 6) Laguna (ANOVA, F=10.73, p>0.05). Among the 12 sampling points, only points 3 and 4 (Site 2, Rio Grande) showed significant differences in mullet abundance (ANOVA, F=14.15; p<0.05).

The size distribution of the M. liza individuals varied between 12 mm and 126 mm TL. Individuals 30 mm TL were more abundant than larger individuals (TL>30 mm) at all sites (ANOVA, F=13.050, p<0.05) (Fig. 4). At all sites, the total number of individuals from classes 15 to 30 (≤30 mm TL) represented more than 96% of all mullets collected, with the exception of Site 3 (Mostardas), where they represented 77% (Table 2). Individuals in size class 15 mm TL (10 to 15 mm TL) occurred occasionally at Chuí (Site 1), Tramandaí (Site 4) and Laguna (Site 6), but individuals from size class 25 mm TL were the modal length at all sampling sites (Table 2). Mullets ≤30 mm TL occurred in all 12 months sampled at Sites 2 and 4 and were present among 9 to 11 of all months sampled at the other sampling sites (Table 2).

The best model for presence/absence of young mullets (i.e. lower AIC) included three significant variables (site, month and water temperature). Those variables explained 33.5% of the total deviance of the model. The water temperature alone explained more than half of the total deviance (23%) (Table 3). Young mullets were present in all seasons, but the probability of occurrence was higher in winter (May to July) and was associated with low temperatures (Fig. 5A, B). The highest probability of occurrence of young mullets was observed at Site 2 (Rio Grande) (Fig. 5C). The best GLM for the abundance data of M. liza included three explanatory variables, sampling site, month and water temperature (Table 3), and explained 39% of the total deviance of the model. Similar to the GLM model for presence/absence, the abundance GLM model showed that water temperature was responsible for more than half of the explained deviance (21%). Young mullets were observed in the surf zone in all months, and abundance was a positive relationship of temperature (Fig. 5D, E). Mullet abundance was highest at Site 2 (Rio Grande) (Fig. 5F).

DISCUSSION

Surf zone habitats form a productive ecosystem, serving important functions as nursery areas, major fishery zones and migratory pathways. The importance of this area for species dependent on the estuary has long been known (Whitfield 1989Whitfield A.K. 1989. Ichthyoplankton in a southern African surf zone: nursery area for the postlarvae of estuarine associated fish species? Est. Coast. Shelf. Sci. 29: 533-547. https://doi.org/10.1016/0272-7714(89)90009-7 ). Vieira (1991)Vieira J.P. 1991. Juvenile mullets (Pisces: Mugilidae) in the Estuary of Lagoa dos Patos, RS, Brazil. Copeia 409-418. https://doi.org/10.2307/1446590 proposed a model for the early life history of M. liza in southern Brazil, suggesting that young individuals use the surf zone temporarily after the larval neustonic life stages and before recruitment to the estuarine environment. Young mullets remain present throughout the year in the surf zone over the entire sampling area. The surf zone is therefore an important environment for young mullets, and long periods in these areas could provide better opportunities for young mullets to find and penetrate estuaries.

The recruitment to estuaries of species that spawn at sea has two major phases (Boehlert and Mundy 1988Boehlert G.W., Mundy B.C. 1988. Roles of behavioral and physical factors in larval and juvenile fish recruitment to estuarine nursery areas. In: American Fisheries Society Symposium 3: 1-67.). The first is the arrival at and movement through the coastal zone. This phase is associated with physical factors responsible for the passive transport of organisms, such as tides, oceanographic fronts and coastal currents (Boehlert and Mundy 1988Boehlert G.W., Mundy B.C. 1988. Roles of behavioral and physical factors in larval and juvenile fish recruitment to estuarine nursery areas. In: American Fisheries Society Symposium 3: 1-67.). The longshore currents in the study region follow the seasonal wind pattern. Strong winds from the south and southwest predominate between April and August (especially during the passage of cold fronts), causing strong currents towards the northeast (Seeliger et al. 1998Seeliger U., Odebrecht C., Castello J.P. 1998. Os Ecossistemas Costeiro e Marinho do Extremo Sul do Brasil. Ed. Ecoscientia, Rio Grande, 337 pp.). The importance of this wind pattern for the reproductive migration of adult individuals is already known (Lemos et al. 2014Lemos V.M., Varela Jr. A.S., Schwingel P.R., et al. 2014. Migration and reproductive biology of Mugil liza (Teleostei : Mugilidae) in south Brazil. J. Fish. Biol. https://doi.org/10.1111/jfb.12452 ). However, during most of the year (September to March), there is a predominance of northeast winds and coastal currents with low intensity but high frequency towards the southwest (Tozzi and Calliari 2000Tozzi H.D.M., Calliari L.J. 2000. Morfodinâmica da praia do Cassino, RS. Pesquisas em Geociências 27: 29-42. https://doi.org/10.22456/1807-9806.20176 ). These coastal currents can be important for the transport of young mullets along the coast towards the south.

The second phase of recruitment to estuaries is the aggregation of individuals near estuary mouths (Boehlert and Mundy 1988Boehlert G.W., Mundy B.C. 1988. Roles of behavioral and physical factors in larval and juvenile fish recruitment to estuarine nursery areas. In: American Fisheries Society Symposium 3: 1-67.). Many of important commercial fish species recruit in estuaries, presumably in response to certain olfactory cues transported into the marine environment from either estuarine or freshwater sources (Wasserman and Strydom 2011Wasserman R.J., Strydom N.A. 2011. The importance of estuary head waters as nursery areas for young estuary-and marine-spawned fishes in temperate South Africa. Est. Coast. Shelf. Sci. 94: 56-67. https://doi.org/10.1016/j.ecss.2011.05.023 ). A strong aggregation of juvenile mugilids was observed in response to the presence of estuarine discharges in surf zones (Strydom 2003Strydom N.A. 2003. Occurrence of larval and early juvenile fishes in the surf zone adjacent to two intermittently open estuaries, South Africa. Environ. Biol. Fish. 66: 349-359. https://doi.org/10.1023/A:1023949607821 ). In southern Brazil, Rodrigues et al. (2015)Rodrigues F.L., Cabral H., Vieira J.P. 2015. Assessing surf zone fish assemblage variability in southern Brazil. Mar. Freshwater. Res. 66: 106-119. https://doi.org/10.1071/MF13210 observed a positive effect of the weaker morphodynamics and the occurrence of an estuarine plume at the beach on the abundance of M. liza juveniles in the surf zone. Several environmental and ecological attributes were assessed in the present study to explain the occurrence and abundance of young mullets in the surf zone, and the abundance near the Patos Lagoon estuary mouth (Site 2) was significantly higher than that at the other sites sampled. The Patos Lagoon is the largest estuary in the study area, and the plume is related to the size of the drainage basin and the high discharge of the lagoon, which are important factors for the high sediment input into the adjacent coastal zone (Figueiredo and Calliari 2006Figueiredo S.A.D., Calliari L.J. 2006. Sedimentologia e suas implicações na morfodinâmica das praias adjacentes às desembocaduras da linha de costa do Rio Grande do Sul. Gravel. 4: 73-87.). In contrast, the sediment input via estuarine discharge was not significant at the beaches adjacent to the Mampituba River, the Tramandaí-Armazém Lagoon Complex and Chuí Creek (Figueiredo and Calliari 2006Figueiredo S.A.D., Calliari L.J. 2006. Sedimentologia e suas implicações na morfodinâmica das praias adjacentes às desembocaduras da linha de costa do Rio Grande do Sul. Gravel. 4: 73-87.). Although water salinity and transparency were not important factors in explaining the occurrence and availability of young M. liza throughout the study area, the sediment input via estuarine discharge can be an important factor for explaining the greater abundance and occurrence of young mullet near the mouth of Patos Lagoon.

The best models for presence/absence and for abundance of young mullets demonstrate that temperature was the factor that most contributed to the explanation of variance. However, temperature showed a negative relationship with the occurrence and a positive relationship with the abundance of young mullet in the study area. Young mullets were present in all seasons, but the probability of presence was higher in winter (May to July) and was associated with low temperatures. This can be explained by the reproductive migration of M. liza, which occurs at the end of the austral autumn and the beginning of the austral winter (from April to July), and the occurrence of one peak spawning in June between northern Santa Catarina and Paraná (26°S) (Lemos et al. 2014Lemos V.M., Varela Jr. A.S., Schwingel P.R., et al. 2014. Migration and reproductive biology of Mugil liza (Teleostei : Mugilidae) in south Brazil. J. Fish. Biol. https://doi.org/10.1111/jfb.12452 ). Similar to the GLM model for presence/absence, the GLM model for abundance showed that the water temperature was responsible for more than half of the explained deviance, but mullet abundance was positively associated with temperature. The abundance of young mullets remained at the same level over a wide temperature range, between 16°C and 24°C.

Juveniles of M. liza have been observed to remain in the surf zone for long periods in southern Brazil and Argentina (Acha 1990Acha E.M. 1990. Estudio anatômico-ecologico de la lisa (Mugil liza) durante su primer año de vida. Frente. Marit. 7: 37-43., Vieira 1991Vieira J.P. 1991. Juvenile mullets (Pisces: Mugilidae) in the Estuary of Lagoa dos Patos, RS, Brazil. Copeia 409-418. https://doi.org/10.2307/1446590 , Castellini et al. 2017Castellini D.L., Brown D., Lajud N.A., et al. 2017. Juveniles recruitment and daily growth of the southern stock of Mugil liza (Actinopterygii; Fam. Mugilidae): new evidence for the current life-history model. J. Mar. Biol. Assoc. U.K. 99: 204-215. https://doi.org/10.1017/S0025315417001904 ). The specific size at which the mullets start to appear in surf zones and enter estuaries varies between authors (González-Castro and Minos 2016González-Castro M., Minos G. 2016. Sexuality and reproduction of Mugilidae. In: Crosetti D., Blaber S. (eds) Biology, ecology and culture of mullets (Mugilidae). CRC Press. Boca Raton. pp. 227-263.). In the present study, we observed that individuals smaller than 15 mm TL occur occasionally throughout the study area. However, individuals ≤30 mm TL represented more than 96% of all mullets collected (with the exception of Site 3, where they represented 77%) and were present in the surf zone at least 9 of the 12 months of the year. This result is in accordance with Vieira (1991), who demonstrated that after reaching approximately 20 mm TL, the young M. liza gradually migrate to the bottom in the surf zone and begin to feed on benthic organisms and mineral particles. At this stage, individuals occupy the surf zone before recruitment to the estuarine environment (Vieira 1991Vieira J.P. 1991. Juvenile mullets (Pisces: Mugilidae) in the Estuary of Lagoa dos Patos, RS, Brazil. Copeia 409-418. https://doi.org/10.2307/1446590 ).

From the larvae to juvenile stage, the feeding habits change from zooplanktophagous to iliophagous (Acha 1990Acha E.M. 1990. Estudio anatômico-ecologico de la lisa (Mugil liza) durante su primer año de vida. Frente. Marit. 7: 37-43., Vieira 1991Vieira J.P. 1991. Juvenile mullets (Pisces: Mugilidae) in the Estuary of Lagoa dos Patos, RS, Brazil. Copeia 409-418. https://doi.org/10.2307/1446590 ). The surf zone occurrence of young mullets appears to be less important at beaches with strong morphodynamics, and in this situation the change in feeding habits seems to occur inside the estuaries (Blaber and Whitfield 1977Blaber S.J.M., Whitfield A.K. 1977. The feeding ecology of juvenile mullet (Mugilidae) in southeast African estuaries. Biol. J. Linn. Soc. 9: 277-284. https://doi.org/10.1111/j.1095-8312.1977.tb00270.x ). The coastal zone in our study area is characterized by an extensive dissipative beach (Calliari and Klein 1995Calliari L.J., Klein A.H.F. 1995. Variabilidade das praias oceânicas entre Cassino e Chuí, RS. IO/USP. 11: 55-63., Ramos and Vieira 2001Ramos L.A., Vieira J.P. 2001. Composição específica e abundância de peixes de zonas rasas dos cinco estuários do Rio Grande do Sul, Brasil. Bol. Inst. Pesca. 27: 109-121.), and at dissipative beaches the transition between feeding habits may occur slowly and in the surf zone (Blaber and Blaber 1980Blaber S.J.M., Blaber T.G. 1980. Factors affecting the distribution of juvenile estuarine and in shore fish. J. Fish. Biol. 17: 143-162. https://doi.org/10.1111/j.1095-8649.1980.tb02749.x ). Considering that M. liza in southern Brazil show only one peak of spawning in June and young mullets occur year-round, it is possible that the young mullets that remained in the surf zone exhibited a lower growth rate than those that were already recruited to the estuary. Thus, long periods in the surf zone could provide better opportunities for young mullets to find and penetrate estuaries, but this would imply lower growth rates. However, comparative studies with different salinities (surf zone and estuarine areas) must be performed to analyse the available information on the biology and physiology of the young stages of M. liza.

Acha (1990)Acha E.M. 1990. Estudio anatômico-ecologico de la lisa (Mugil liza) durante su primer año de vida. Frente. Marit. 7: 37-43. and Castellini et al. (2017)Castellini D.L., Brown D., Lajud N.A., et al. 2017. Juveniles recruitment and daily growth of the southern stock of Mugil liza (Actinopterygii; Fam. Mugilidae): new evidence for the current life-history model. J. Mar. Biol. Assoc. U.K. 99: 204-215. https://doi.org/10.1017/S0025315417001904 also reported young mullets along the Argentinian coast throughout the year. González-Castro et al. (2011)González-Castro M., Macchi G.J., Cousseau M.B. 2011. Studies on reproduction of the mullet Mugil platanus Günther, 1880 (Actinopterygii, Mugilidae) from the Mar Chiquita coastal lagoon, Argentina: similarities and differences with related species. Ital. J. Zool. 78: 343-353. https://doi.org/10.1080/11250003.2010.549154 described a strong reproductive period in the austral winter (April-May), similar to that reported for southern Brazil, and a second period in the austral summer (November-December), which was not observed by Lemos et al. (2014)Lemos V.M., Varela Jr. A.S., Schwingel P.R., et al. 2014. Migration and reproductive biology of Mugil liza (Teleostei : Mugilidae) in south Brazil. J. Fish. Biol. https://doi.org/10.1111/jfb.12452 in southern Brazil. Those two periods of reproduction were used by Castellini et al. (2017)Castellini D.L., Brown D., Lajud N.A., et al. 2017. Juveniles recruitment and daily growth of the southern stock of Mugil liza (Actinopterygii; Fam. Mugilidae): new evidence for the current life-history model. J. Mar. Biol. Assoc. U.K. 99: 204-215. https://doi.org/10.1017/S0025315417001904 to explain the occurrence of young mullets with a minimum size of between 19 and 33 mm SL along the Argentinian coast throughout the year. However, Lemos et al. (2014)Lemos V.M., Varela Jr. A.S., Schwingel P.R., et al. 2014. Migration and reproductive biology of Mugil liza (Teleostei : Mugilidae) in south Brazil. J. Fish. Biol. https://doi.org/10.1111/jfb.12452 reviewed and described the reproductive biology of M. liza in southern Brazil and showed that the peak in the gonadosomatic index (IG) values for females caught in marine waters (by industrial purse seine fisheries) occurred in June (mean±sd IG♀=11.5±3.1) and was closely related to the presence of hyaline oocytes. The gonadosomatic index values for fish caught in Argentina are lower (4 to 6) according to González-Castro et al. (2011)González-Castro M., Macchi G.J., Cousseau M.B. 2011. Studies on reproduction of the mullet Mugil platanus Günther, 1880 (Actinopterygii, Mugilidae) from the Mar Chiquita coastal lagoon, Argentina: similarities and differences with related species. Ital. J. Zool. 78: 343-353. https://doi.org/10.1080/11250003.2010.549154 , who did not identify the presence of hyaline oocytes. The hydration phase during final maturation appears to be important in the production of pelagic (buoyant) mullet eggs, and the presence of hyaline oocytes has been associated with high salinity and sea surface temperatures of 19 to 21°C (Lemos et al. 2014Lemos V.M., Varela Jr. A.S., Schwingel P.R., et al. 2014. Migration and reproductive biology of Mugil liza (Teleostei : Mugilidae) in south Brazil. J. Fish. Biol. https://doi.org/10.1111/jfb.12452 ). Recently, Morado et al. (2021)Morado C.N., de Andrade-Tubino M.F., Araújo F.G. 2021. Local ecological knowledge indicates: There is another breeding period in the summer for the mullet Mugil liza in a Brazilian tropical bay. Ocean. Coastal. Manage. 205: 105569. https://doi.org/10.1016/j.ocecoaman.2021.105569 describes a second breeding period in the summer for the mullet M. liza in a Brazilian tropical bay (Sepetiba Bay coast of the State of Rio de Janeiro). The summer spawning was not foreseen in the biological literature for this place, but the local ecological knowledge of artisanal fishermen associated this spawning season with a possible resident population (Morado et al. 2021Morado C.N., de Andrade-Tubino M.F., Araújo F.G. 2021. Local ecological knowledge indicates: There is another breeding period in the summer for the mullet Mugil liza in a Brazilian tropical bay. Ocean. Coastal. Manage. 205: 105569. https://doi.org/10.1016/j.ocecoaman.2021.105569 ). More studies must be carried out to support the hypothesis of a local summer spawning event in southern Brazil.

In summary, the results of our study demonstrate that the probability of occurrence of young M. liza throughout the sampling area is higher in winter (May to July) and associated with low temperatures; however, young mullets remain present throughout the year in the surf zone over the entire sampling area. Site 2, adjacent to the Patos Lagoon estuary, had the highest probability of occurrence of all collection sites and the greatest abundance (n=3402) of young M. liza, which represented 52% of the total of individuals (n=6493) caught among the six sites sampled. Although other estuaries along the southern Brazilian coast are also used by mullets as nursery grounds (Ramos and Vieira 2001Ramos L.A., Vieira J.P. 2001. Composição específica e abundância de peixes de zonas rasas dos cinco estuários do Rio Grande do Sul, Brasil. Bol. Inst. Pesca. 27: 109-121., Loebmann et al. 2005Loebmann D., Vieira J.P. 2005. Distribuição espacial e abundância das assembléias de peixes no Parque Nacional da Lagoa do Peixe, Rio Grande do Sul, Brasil. Ver. Bras. Zool. 22: 667-675. https://doi.org/10.1590/S0101-81752005000300023 , Favero and Dias 2015Favero J.M. del, Dias J.F. 2015. Juvenile fish use of the shallow zone of beaches of the Cananéia - Iguape coastal system, southeastern Brazil. Braz. J. Oceanogr. 63: 103-114. https://doi.org/10.1590/S1679-87592015081806302 ), our results indicate that the Patos Lagoon estuary can be considered the largest nursery area for the southern population of M. liza, and corroborate the perception of local fishermen that Patos Lagoon is the largest nursery area for M. liza (Herbst and Hanazaki 2014Herbst D.F., Hanazaki N. 2014. Local ecological knowledge of fishers about the life cycle and temporal patterns in the migration of mullet (Mugil liza) in Southern Brazil. Neotrop. Ichthyol. 12: 879-890. https://doi.org/10.1590/1982-0224-20130156 ). Comparative studies between surf zone and estuary growth rates, based on daily increment growth in otoliths, morphological development and changes in feeding habits in the M. liza ontogeny, should be performed to support some of the hypotheses raised in this work.