Dynamics of fish larvae recruitment in the hydrographic basin of the Paraguay River in western Brazil

Abstract The conservation of fish assemblages depends on the longitudinal and lateral connectivity between riverine habitats, in particular during the breeding season and the initial development. This study investigated the composition and spatio-temporal structure in the ichthyoplankton of the hydrographic basin of the Paraguay River in western Brazil to identify the local spawning grounds and the importance of the longitudinal connectivity of economically valuable migratory species. Data were collected at 10 sites between two breeding seasons (2017/2018 and 2018/2019). Were collected 8,635 larvae, representing 55 taxa in 25 families, including the migratory species Brycon hilarii, Hemisorubim platyrhynchos, Piaractus mesopotamicus, Prochilodus lineatus, Pseudoplatystoma spp., Salminus brasiliensis, Sorubim lima, and Zungaro jahu, which are important fishery resources, with the highest larval densities being recorded between November and January. The Sepotuba, Paraguay, Jauru, and Cabaçal rivers were the areas of greatest connectance, and should thus be considered critical for the conservation of the longitudinal connectivity of this fluvial system, indicating that the migratory species spawn upriver. More efficient fisheries management mechanisms are needed, respecting the spawning period of migratory species, maintaining quality and longitudinal connectivity between habitats, and characteristics necessary for successful larval recruitment.


INTRODUCTION
Rivers are complex aquatic ecosystems that typically encompass a considerable heterogeneity of lotic and lentic habitats, as well as backwaters and floodplains, which may vary in their distribution and connectivity (Nannini et al., 2012).The fish exploit these fluvial biotopes during different stages of their life cycle, in particular during spawning and initial development (Kallasvuo et al., 2017), which guarantee larval recruitment.However, the rivers of the Neotropical region are coming under increasing pressure from hydroelectric dams, which have impacted many hydrographic basins on a large scale (Suzuki et al., 2016).Hydroelectric power plants impose different barriers on the dispersal of many fish, in particular, long-distance migrants and rheophilic species.The impoundment of the river impedes movements upriver from the dam, while the reservoir affects primarily the movements downriver from the dam, through the reduction of lotic environments and river flow changes (Pelicice et al., 2015).
During the reproduction season, long-distance migrants search for spawning habitats, which are typically found in upriver locations (Makrakis et al., 2012).This results in longitudinal movements along with the river system toward its headwaters, where the fish spawn in the main channel of the river and the eggs drift passively downriver, with the larvae hatching downstream (Lopes et al., 2014), and subsequently disperse further until they reach an area appropriate for their early development (Baumgartner et al., 2004).The success of the larval recruitment will depend on the drift of the lotic waters as far as swampy riparian habitats, where the larvae encounter feeding resources and refuges from predators (Suzuki et al., 2009).The spawning period of the longdistance migrants coincides with the flood pulse and relatively high temperatures (Zaniboni-Filho et al., 2017), when new habitats become available and floodplain lakes are connected to the main channel of the river through the inundation of marginal environments (Junk et al., 1989).In this scenario, the recruitment of the juveniles into the adult population will depend on the preservation of both the spawning grounds and the nurseries in which the larvae grow and develop (Bower et al., 2015;Teixeira et al., 2019), as well as the connectivity between these different habitats within the fluvial system (Branco et al., 2016).
The Brazilian upper Paraguay River basin in encompasses two distinct sectorsthe plateau and the Pantanal floodplain.The plateau rivers fed the Pantanal wetlands' flood pulse, forming a complex aquatic ecosystem.Modern-day complication is the increasing impoundment of the tributaries of the Paraguay River for the implantation of hydroelectric power plants, which not only affects the hydrological regime of the Pantanal but also interrupts the connectivity of the fluvial systems, which reduces the access of fish to the nursery habitats and feeding grounds on the floodplains (Winemiller et al., 2016;Alho, Reis, 2017).This reinforces the need to identify the key areas of larval recruitment of migratory fish, given that these species require long stretches of wellpreserved riverine habitat in which to complete their life cycle.It is also important to identify the importance of the connectivity of the habitats used by the fish species during the initial phases of their life cycle, in particular in areas subject to impoundment for hydroelectric projects (Agostinho et al., 2002).The temporal variation in using of the different habitats is an additional facet of the life cycle of the species, which helps to define the spawning period and provides important guidelines for the establishment of protective measures for the breeding stocks.
In this context, the spatio-temporal analysis of ichthyoplankton distribution provides important insights into the location of the spawning grounds and nursery areas of fish species (Baumgartner et al., 2004), and in particular, the dispersal process, which is fundamental to the reproductive success of these species (Barthem et al., 2014).Furthermore, this study uses the abundance of fish larvae to evaluate the reproductive patterns of different species in the Brazilian upper Paraguay River basin.Specifically, the study (i) analyzes the temporal and spatial distribution of the fish larvae, (ii) determines the taxonomic composition of the fish larvae, with emphasis on long-distance migrants and species with a high commercial value, (iii) detects the spawning grounds of the species, in particular the most valuable fishery resources, and the dispersal of their offspring, and (iv) demonstrates the importance of longitudinal connectivity for the reproduction of commercially valuable migratory species.We predicted that the fish larvae disperse passively within the fluvial system according to the migratory patterns of the spawning stock, and tested two hypotheses: (i) that there will be a higher density of fish larvae at the sampling sites further downstream within the Paraguay River basin (i.e., on the Cabaçal [site CAB2], Jauru, Sepotuba [site SEP3], and Paraguay rivers), and (ii) that the migratory species will depend on the longitudinal connectivity of the fluvial system to ensure their spawning migrations, with the larvae subsequently dispersing to the marginal areas of the floodplain.

MATERIAL AND METHODS
Study area and sampling.In Brazil, there are two different sectors of the upper Paraguay River basin, the plateau and the Pantanal floodplain, with distinct sets of natural and hydrological resources (Alho, 2008).The plateau is located to the north and east of the upper Paraguay River basin, and covers an area of 214,802 km 2 , about 59% of the basin in Brazil (Galdino et al., 2006).This area encompasses the headwaters of the principal rivers that flow toward the Pantanal floodplain.The main tributaries of the upper Paraguay River basin include the Paraguay, Jauru, Sepotuba, Cuiabá, and São Lourenço rivers.This region has a seasonal humid tropical climate, with two welldefined seasons, rainy and dry, and mean annual precipitation of 1,400 mm, most that falls between November and March (Fantin-Cruz et al., 2011).
The study area is located on the Paraguay River and its right-margin tributaries on the plateau, a region with several future hydroelectric projects still at the planning stage (Silva et al., 2007;Souza et al., 2017;Calheiros et al., 2018).A total of 10 sampling sites were established within this area (Fig. 1) on the Cabaçal (sites CAB1 and CAB2), Formoso (FOR), Jauru (JAU), Juba (JUB), Sepotuba (sites SEP1, SEP2, and SEP3), Vermelho (VERN), and Paraguay (PAR) rivers, all in the state of Mato Grosso, in western Brazil.Collections considered each water body distinct from each other due to the motionless characteristics of the ichthyoplankton, with one site in the tributary and another in the main river, upstream of the tributary mouth.However, the overlap of the FOR/SEP1, JUB/SEP2, and SEP3/PAR sites corresponds to the proximity of the collected water bodies (see Fig. 1).On the two rivers with impoundments -the Jauru and Juba -the sampling sites were located in non-flooded stretches downriver from the dams because of the loss of connectivity due to the projects.Sites CAB2, JAU, SEP3, and PAR were considered representing downriver areas in the river system studied here (Tab.S1).
To investigate the spatial-temporal variation, monthly samplings were carried out in each of the 10 sites, in two consecutive breeding seasons, November-March 2017/2018 (RP1) and October-March 2018/2019 (RP2), respecting the reproductive cycle of the regional ichthyofauna (Ziober et al., 2012).Ichthyoplankton was collected using conical-cylindrical plankton nets with a mesh size of 0.5 mm and 0.38 m 3 diameter, with a flowmeter attached to measure the volume of filtered water.To encompass the spatial variation inside each site, the surface layer of the water column (right and left margins, and center channel) and at the bottom were sampled using nets submerged for 10 and 3 min, respectively.For bottom samples, we use adapted conical-cylindrical plankton nets coupled to a sled-shaped iron structure for sites over 1.5 m deep.The sampling was performed at 6 h intervals (00, 06, 12 and 18h) for to include the nictemeral variation.These samples were grouped to represent the spatio-temporal variation inside the sites.Thus, for each month sampled, 10 replicates were obtained.Ichthyoplankton was fixed in 4% formalin, and buffered with calcium carbonate (CaCO 3 ) to conservation of the biological material.In the laboratory, the larvae were separated and quantified using a stereomicroscope.The abundance of the larvae was standardized to a volume of 10 m 3 of filtered water (Tab.S2), following Tanaka (1973), as modified by Nakatani et al. (2001).The larvae were identified using the regressive developmental sequence technique, which involves the morphological comparison of the smaller individuals  River system connectivity and recruitment fish with more developed forms, using traits such as body shape, the presence of barbels, the sequence of the formation of the fins, the relative position of the anal opening, and the number of vertebrae/myomeres or fin rays (Ahlstrom, Moser, 1976;Nakatani et al., 2001), which was supported by literature (Oldani 1979 a,b;Nascimento, Araújo-Lima, 2000;Nakatani et al., 2001;Oliveira et al., 2012;Andrade et al., 2014Andrade et al., , 2016;;Taguti et al., 2015;Garcia et al., 2016).Larvae were identified to the lowest possible taxonomic level (larvae identified to at least genus included some long-distance migrant species) (Suzuki et al., 2004).Once identified, the larvae were assigned to one of the four development stages (yolk-sac, preflexion, flexion, and postflexion) following Ahlstrom et al. (1976), as modified by Nakatani et al. (2001), considering the flexion of the notochord and its supporting elements.
Data analysis.Spatio-temporal distribution.We used a longitudinal design, where the same sampling units (10 sites) were sampled over time (months).The analyzes were conducted considering this design to adjust the degrees of freedom since samples in the same sample unit over time are not independent.After testing the normality and homoscedasticity of the data, a nonparametric test with inter-sample dependence (nested month by sampling site) was used.The differences in larval density (i.e., individuals/10 m 3 ) among sites and between periods were verified using the Friedman (repeated-measures analysis of variance -ANOVA) nonparametric test (MacFarland, Yates, 2016), to avoid temporal pseudoreplication, with Dunn's post hoc test for pairwise comparisons (Dunn, 1964).
Taxonomic composition.We evaluated the differences in the taxonomic composition of the samples using a PERMANOVA based on a Bray-Curtis matrix, with 9999 permutations, considering all the larvae together, independent of the level of their taxonomic identification.A permutation test of the homogeneity of multivariate dispersal was applied to test the significance of the differences in the taxonomic composition of the samples at a spatio-temporal scale (Anderson et al., 2006).The dispersal of the data on the composition of the larval assemblage was used to calculate the homogeneity of the multivariate dispersal (PERMDISP), considering the dissimilarity of a sample concerning its distance from the centroid of the entire group of samples.The values were then configured spatially using a Principal Coordinates Analysis (PCoA).The spatial autocorrelation in the taxonomic composition of the larvae with the distance between the sampling sites was verified using a partial Mantel test.Euclidean distance, with the geographical coordinates transformed into UTM, was used as a matrix of geographic distance, while the Bray-Curtis distance was applied for the taxonomic composition.
Spawning.An indicator species analysis (IndVal; Dufrene, Legendre, 1997) was used to identify associations at spatial (Cáceres, 2022) and temporal scales with 999 permutations.The larval development stages of the commercially valuable migratory species were related to the spatial scale, which permitted the inference of the possible location of the spawning grounds and the dispersal of the fry, while the temporal scale was evaluated using only the larvae that were identified to species.
Longitudinal connectivity.We used a bipartite network approach to evaluate the importance of longitudinal connectivity for economically-valuable migratory species in the upper Paraguay basin (Mateus et al., 2004;Alho, Reis, 2017).This type of approach 7/21 ni.bio.br| scielo.br/niTatiane P. Sousa, Andréa Bialetzki and Lucia Aparecida F. Mateus is typically used for the analysis of pollinator or predator-prey networks (Dormann et al., 2008(Dormann et al., , 2009)), but may also be applied to other types of interaction.Here, we used metrics for quantitative matrices, considering the density of the migratory species identified at each sampling site.We then compiled a bipartite plot, which shows the relationship between the density of species and the sampling sites, to visualize the interaction network.
To evaluate whether the interruption of longitudinal connectivity, by impoundment, for example, impedes the access of fish to spawning grounds upriver from a sampling site, we simulated the elimination of connections beginning with the sites with a greater frequency of connections and estimated the coefficient of extinction, 'a' (Memmott et al., 2004).We assumed that the proportion of larvae (y) reflects an exponential function of the proportion of extinction (disconnected) sites, whose slope is estimated by coefficient 'a'.The discontinuity of the sites affects the greater the value of 'a' (Dormann et al., 2009).We repeated the procedure, this time beginning with the elimination of the sites with the lowest density of connections, as well as the random discontinuity of the sites.
We also evaluated the robustness of the system to the loss of connectivity.Burgos et al. (2007) developed a quantitative measure of robustness, based on a parameter, R, which is defined as the area under the extinction curve.Intuitively, R = 1 corresponds to a curve that declines gently to the scenario in which almost all the sites are eliminated (unconnected).This is consistent with a very robust system in which, for example, most of the species persist, even after a large proportion of the sites have been disconnected.By contrast, R = 0 represents a curve that declines abruptly as soon as any one connection is lost.This is consistent with a fragile system in which, for example, most of the species lose all their interactions (that is, the potential to reach their spawning grounds) and become extinct, even when only a small proportion of the sites have become disconnected (Burgos et al., 2007;Dormann et al., 2009).
The estimates of 'a' and 'R' were compared with the null models, which were generated based on 1000 simulations, using a randomization algorithm that maintains constant the sum of the lines and columns of the matrix (Oksanen et al., 2019).All the analyses were run in the R environment (R Development Core Team, 2020), using the indicspecies (Cáceres, 2022), FD, lmPerm (Wheeler, Torchiano, 2016), vegan (Oksanen et al., 2019), and bipartite (Dormann et al., 2020) packages, considering a significance level of p < 0.05.

RESULTS
Spatio-temporal distribution.These samples resulted in 8,635 larvae.The Friedman test indicates significant effects of spatial scale (X 2 (9) = 27.644;p = 0.001) and the sample month (X 2 (5) = 27.248;p < 001) on the larval density.Dunn's post hoc test (Tab.S1) showed higher significant larval densities (p < 0.05) in downstream areas, between CAB2, JAU, SEP3, PAR, and VERN and all the other sampling sites.Larval densities recorded in November, December, January, and February were all significantly different (p < 0.05) from October and March.The reproductive peak occurred between November and February in both reproductive periods with the highest larval densities, indicating the onset of spawning between October and late March (months with the lowest larval density, Fig. S3).The results of the PERMANOVA revealed significant differences in the total taxonomic composition among the sample sites (F 9;26 = 1.62; p < 0.001) and at a temporal scale (F 5;70 = 1.40; p = 0.005).However, no significant variation was found in the heterogeneity of the species composition at either spatial or temporal scales (Figs.2-3).The partial Mantel test detected a significant correlation between the matrix of the taxonomic composition of the larvae and the geographical distance between sampling sites (r = 0.23; p < 0.001), as well as for the composition of the commercially valuable migratory species (r = 0.19; p = 0.003), which indicates that sites closer together have a more similar composition of larval taxa.River system connectivity and recruitment fish Spawning.The commercially valuable migratory species presented significant patterns of the association at both spatial and temporal scales (Tab.2; p < 0.05).The flexion larval stage of B. hilarii was associated with the Paraguay River.Brycon hilarii larvae in yolk-sac stage and Pseudoplatystoma spp. in flexion stage were both associated with site SEP3.The preflexion stage of P. were associated with the sites JAU/ PAR and the yolk-sac stage of P. mesopotamicus, with sites JAU/SEP3, while the yolksac stage of Z. jahu were associated with sites SEP1/SEP2.The preflexion larvae of H. platyrhynchos were associated with sites SEP3/VERN.On a temporal scale, P. mesopotamicus was associated with November and December.Considering only the variation in the larval development stages (p < 0.05), the pre and postflexion larvae were associated with sites CAB2/JAU/PAR/SEP3 (COMB1), whereas the larvae in flexion were associated with sites CAB2/JAU/PAR/SEP3/VERN (COMB2) in the Tab. 2.
Longitudinal connectivity.The migratory-local larval density network (Fig. 4) is affected by the loss in connectivity among the sites, with low values of the extinction coefficient 'a', not only for the scenario in which the sites with the highest frequency of interactions are disconnected (a = 0.72; a̅ null = 3.35; sd = 0.24; p < 0.0001) but also for the scenarios in which the disconnection begins with the sites of lower density (a = 0.718; a̅ null = 3.36; sd = 0.23; p < 0.001) or is completely random (a = 0.55; a̅ null = 3.44; sd = 0.21; p < 0.001).In all cases, the estimate of the value of 'a' is much lower than that expected    by chance (p < 0.0001), which shows that losing of connectivity is likely to have a major impact on the density of migratory species larvae.Besides to the sensitivity of these species to extinction, the network indicated low levels of robustness when the loss of connectivity begins with the sites with the highest frequency of interaction (R = 0.49), which reflects the vulnerability of the system to the loss of connectivity in even a small proportion of the sites.By contrast, when the disconnection begins with the sites with a low frequency of interaction, the system is more robust (R = 0.78).This measure of robustness (R) was not different from that expected by chance, however, whether losing connectivity started with the sites with greater densities of connections (R null = 0.46; sd = 0.027; p = 0.85) or with those of lower density (R null = 0.74; sd = 0.010; p = 0.99).When the simulation of the loss of connectivity was random, the robustness was intermediate and, once again, not different from that expected by chance (R = 0.67; R null = 0.63; sd = 0.012; p = 0.99).River site 1; 2 = Cabaçal River site 2; FOR = Formoso River; 4 = Jauru River; 6 = Paraguay River; 7 = Sepotuba River site 1; 8 = Sepotuba River site 2; 9 = Sepotuba River site 3; 10 = Vermelho River.

DISCUSSION
Our findings show the existence of significant spatio-temporal variation in the density of larvae over the course of the breeding season, with a fundamental association with geographic patterns, both for the assemblage and for the migratory species alone.The high larval densities at the downstream sites sampled within the study area (CAB2, JAU, PAR, and SEP3) highlight the importance of these rivers for larval recruitment and, ultimately, the preservation of the local fish fauna.Furthermore, the high density of the larvae collected in the study suggests that many taxa use spawning sites within this area between November and February, with the larvae subsequently dispersing to the habitats further downriver.Sampling took place during the period when most species of regional ichthyofauna are breeding (Resende, 2003).Despite considering the theoretical references, the non-inclusion of annual collection can be interpreted as a limitation of this work.
The presence of the initial embryonic development stages (yolk-sac and preflexion) indicates that the species are spawning near the sampling sites and confirms the drift of the larvae to nursery areas on the floodplain.It is also important to note that the species composition varied among the sites, with migratory species, including important fishery resources, being indicators of specific sampling sites, in particular, Paraguay and Sepotuba rivers (SEP3).It is worth mentioning that the sampling took place only in the rainy season and it is not possible to estimate the species that reproduce in the dry season, which does not guarantee a complete understanding of the reproduction of the ichthyofauna.Then, the interaction network that encompasses both species and 15/21 ni.bio.br| scielo.br/niTatiane P. Sousa, Andréa Bialetzki and Lucia Aparecida F. Mateus sampling sites highlighted the fragility of the system and the importance of preserving the longitudinal connectivity of the fluvial system to guarantee the reproduction of commercially valuable migratory species.
The seasonal variation in the density of larvae in the study area reveals that the breeding season begins in October (the austral spring), and peaks between November and February, with a minor peak in March, representing the end of the breeding season.The groups of fish that make up the ichthyofauna of the Pantanal have several reproductive strategies, with the commercially valuable migratory species of the order Characiformes spawning at the onset of the breeding season, and those of the order Siluriformes spawning later in the season (Vazzoler, 1996).Even considering the absence of annual sampling or at other times of the year, the sampling in this work reflected the reproductive period for most species in the Paraguay River basin, as in other works in the literature.Clearly, each group of fish presents distinct patterns of breeding and larval development, which allow the whole assemblage to reproduce successfully at the end of the breeding season (Silva et al., 2019).The decline in larval density observed in March may be related to the last phase of the breeding season in the region (Tondato et al., 2010;Reynalte-Tataje et al., 2011;Ziober et al., 2012).These temporal data are thus essential for establishing protective measures that ensure the survival of spawning stocks, particularly through the implementation of adequate closed seasons by the regional fishery management authorities to protect these stocks.
The abundance of larvae of migratory species in downstream sites and taxonomic similarities with geographic space are evidence of longitudinal connectivity.The rivers constitute a dendritic network of channels connected by the downstream flow and distinct ecological processes, which may influence the dynamics of both the populations and community through the movement of the different species through the network (Grant et al., 2007).The reproductive success of the fish depends on the preservation of the connectivity of the different bodies of water to ensure the dispersal of individuals within the fluvial system (Branco et al., 2016;Humphries et al., 2020).The migratory fish species require adequate conservation measures that include the maintenance of uninterrupted stretches of river (Pelicice, Agostinho, 2008) and the longitudinal connectivity of spawning grounds and the lateral connectivity of areas of larval recruitment, as well as other types of habitats (Pelicice et al., 2015).However, our findings showed clearly that losing longitudinal connectivity induces the loss of species, reflecting the vulnerability of the network to extinction (low 'a' values).In this context, the sites with the highest level of interaction among the commercially valuable migratory species (SEP3, PAR, and CAB2) and larval abundance (JAU) should act as focal sites for the conservation of longitudinal connectivity in the upper Paraguay basin.
Despite the sampling points being downstream of the dams, these projects become physical barriers to reproductive migration as well as ecological barriers to ichthyoplankton.These impoundments interrupt the movements of migratory species during the spawning period (Fernandes et al., 2009;Barzotto et al., 2015) and retain the ichthyoplankton in the reservoir, isolating them from the habitats necessary to complete their life cycle (Pelicice, Agostinho, 2008;Pelicice et al., 2015;Winemiller et al., 2016;Brambilla et al., 2020).Additionally, the fluvial system was not robust ('R') to the loss of connectivity, even in widely distinct scenarios, with a lesser or greater density of connections or even a random arrangement.The reduced resilience of the network 16/21 ni.bio.br| scielo.br/ni River system connectivity and recruitment fish should be considered a key point for the future implantation of new hydroelectric projects because so many species spawn upriver from the sampling sites.Given this, the alteration and fragmentation of habitats may cause the loss of local species, which cannot be substituted by the process of dispersal, resulting in a simultaneous loss of species richness and an increase in the dissimilarity of the communities (Edge et al., 2017).
The associations observed at the spatio-temporal scale highlight the need to preserve the quality, complexity, and longitudinal connectivity of the region's rivers to ensure adequate larval development and dispersal to the floodplain downstream.Migratory species can be used as indicators to identify the areas with characteristics intrinsically important for the adequate development and recruitment of the larvae, and the different stages of larval development.Given this, the administrators of freshwater systems must consider the migratory behavior of the fish into account, with their use of habitats within areas of overlap, to ensure the conservation of the largest possible number of species (Bower et al., 2015) through the preservation of the longitudinal and lateral connectivity among the different habitats (Branco et al., 2012;Pelicice et al., 2015;Cruz et al., 2020).It will also be essential to integrate scientific and traditional knowledge in the development of public policies that ensure the most effective fishery management strategies, to guarantee the reduction of conflicts of interest and socioeconomic inequalities (Schulz et al., 2019), as well as the reproductive success and the maintenance of the stocks of the region's ichthyofauna.
The larval densities observed in the sampling sites, therefore, corroborate the hypothesis that the highest densities are observed in the areas further downstream (CAB2, JAU, SEP3, PAR, and VERN) due to passive dispersion of ichthyoplankton.This suggests that larvae can find refuge and growth areas toward the floodplain, favoring the recruitment of ichthyofauna.Although the Jauru River presents hydropower reservoirs upstream of the sampled area, the dam-free stretch offers suitable conditions for fish reproduction.However, losing connectivity in the dendritic network will have a great impact on the density of migratory larvae, corroborating the second hypothesis of this study.However, the simulation of losing connectivity of the fluvial system reflects in reduced population recruitment, interfering with the population dynamics and in the fishery production in the region.

FIGURE 3 |
FIGURE 3 | Results of the Principal Coordinates Analysis (PCoA) based on the Bray-Curtis dissimilarity matrix of the taxonomic composition of the fish larvae in relation to the months sampled in the upper Paraguay River basin in Mato Grosso State, western Brazil.OCT = October (black circle), NOV = November (red circle), DEC = December (green circle), JAN = January (black triangle), FEB = February (red triangle), MAR = March (green triangle).