Journal of the Marine Biological Association of the United Kingdom, 2012, 92(5), 1117 – 1126. doi:10.1017/S0025315411001718 # Marine Biological Association of the United Kingdom, 2012 Copepoda Harpacticoida community of a rocky shore under the influence of upwelling (Arraial do Cabo, southeastern Brazil) v.c. sarmento1, l.m. lage2 and p.j.p. santos1 1 Universidade Federal de Pernambuco, Centro de Ciências Biológicas, Departamento de Zoologia, Avenida Prof. Morais Rêgo s/n, 50670-420, Recife, PE, Brazil, 2Instituto de Estudos do Mar Almirante Paulo Moreira, Rua Kioto, 253, Praia dos Anjos, 28930-000, Arraial do Cabo, RJ, Brazil Upwelling can determine important changes in community structure on rocky shores. However, studies on the meiofauna in areas of upwelling are scarce and do not include analysis on the species level for Copepoda Harpacticoida. The present study aimed to describe the Harpacticoida community in a region under the influence of coastal upwelling (Arraial do Cabo, Rio de Janeiro, Brazil). The hypothesis tested was that temporal differences in the fauna would be greater at the rocky shore more affected by upwelling. Samples were collected from the sublittoral of two rocky shores: Sonar, which is more exposed to upwelling, in January and June 2004; and Pedra Vermelha, which is a sheltered shore, in March and September 2002. Each sample consisted of four replicates collected by scraping epilithic algae along with the associated sediment. Weekly data on surface water temperature and nitrate content denoted stronger upwelling in 2002 compared to 2004. Nine genera of Copepoda Harpacticoida are reported for the first time for the Brazilian coast. Univariate indices identified no differences between rocky shores or seasons, whereas multivariate analysis indicated significant differences in assemblages between shores and between seasons. The strong variation in physicochemical conditions associated with upwelling favoured the dominance of opportunistic species, such as Parastenhelia spinosa. At the Sonar rocky shore, temporal differences were significantly stronger than at Pedra Vermelha, thereby confirming the initial hypothesis. These results indicated that the occurrence of upwelling had a very important role structuring the Harpacticoida assemblage at the species level. Keywords: algae, meiobenthos, biodiversity, seasonality, isocommunities, South America, upwelling, Arraial do Cabo Submitted 26 May 2011; accepted 5 September 2011; first published online 30 January 2012 INTRODUCTION Rocky shores are ecosystems that support a wide variety of marine animals, such as crustaceans and juvenile fish of great ecological and economic importance (Henriques & Almada, 1998; Thompson et al., 2002). These shores can be widely covered by algae forming a phytal environment that has a characteristic composition and distribution of flora and fauna (Giere, 2009). On rocky shores, environmental gradients of nutrients and productivity have major effects on the type and strength of interactions in communities by establishing strong correlations between species distribution and physical factors, such as temperature, wave exposure, degree of sedimentation and upwelling (Schiel, 2004). Most rocky shores are strongly influenced by adjacent coastal basins hydrography, especially in areas where upwelling occurs. The occurrence of upwelling, defined as the rise of the richest deep, cold water to the surface, affects the community structure on rocky shores, as interactions between organisms on all trophic levels are strongly influenced by the primary production of the ecosystem (Menge, 2000). Corresponding author: V.C. Sarmento Email: visnu.ubi@gmail.com Biodiversity and density distribution patterns for marine benthic communities in areas under the influence of upwelling have been reported in various locations around the world, mainly for macrobenthos. For the macrofauna of the Atlantic coast of tropical Africa, Le Loeuff & Cosel (1998) report that species richness is greater in regions of upwelling when compared to typical tropical regions. On the north-western coast of the Iberian Peninsula, Lastra et al. (2006) also report that the number of macrofauna species increases significantly on beaches near the coastal upwelling zones. On the western coast of southern Africa, Bosman et al. (1987) report that intertidal communities cannot be separated clearly into two groups based on presence or absence of upwelling. Considering temporal variations Witman & Smith (2003) report that the turnover of the community of epifaunal invertebrates at Gordon Rocks in the Galapagos Marine Reserve is extremely fast for a tropical sublittoral community of a natural hard substrate. A twofold increase in the number of species and density in the space of a year in this location of upwelling in the Galapagos is among the most rapid increases in diversity among any sublittoral epifaunal community. For a community of marine algae, Ormond & Banaimoon (1994) confirm that there is a very marked seasonal pattern of seaweed growth, with periods of maximal growth occurring in summer/autumn, when there are high amounts of nutrients due to the seasonal upwelling that occurs along the southern Arabian coast. Temporary episodes of downwelling/upwelling 1117 1118 v.c. sarmento et al. are also known to significantly affect the supply and settlement of invertebrate larval stages (Almeida & Queiroga, 2003; Ma, 2005). The meiofauna community is represented by benthic metazoans that are trapped between meshes of 0.044 mm (or 0.062 mm) and 0.5 mm (or 1 mm) (Giere, 2009). This community plays an important role in the transfer of energy in benthic systems, serving as food for macrobenthos, fish and other organisms (Coull, 1988). Its distribution pattern associated with sediment or algal substrate is very complex due to the influence of geological, chemical, hydrodynamic and various biological processes (Giere, 2009). Among the groups of meiofauna, Copepoda Harpacticoida is regularly the most abundant taxon in the phytal environment (Hicks, 1977a; Coull et al., 1983; Hall & Bell, 1993), with high diversity values (Hicks, 1985). According to Wells (2007), there are about 4300 species in 589 genera and 56 families described for this important taxon. Studies on the meiofauna community in areas that undergo periodic upwelling are rare and include only one study on major taxonomic groups (Soltwedel, 1997) and a few for selected taxa. In one of the few studies on meiofauna (using meshes of 0.250 and 0.177 mm and thus addressing partially only the meiofauna organisms) in upwelling areas at Cabo Frio, Brazil, Machado et al. (2005) suggest that both Ostracoda species composition and community structure are probably controlled by the presence of cold water associated with intense upwelling. The authors further state that the Ostracoda fauna may be characterized as a cold water assemblage. Nearly in the same area (Arraial do Cabo, Brazil), Fonsêca-Genevois et al. (2004) studied Nematoda, concluding that the considerable genus richness is related to upwelling, which determines significant temporal variations in the assemblage structure of this group. To our knowledge, there are no studies evaluating the composition or structure of the Copepoda Harpacticoida assemblage on the species level in phytal communities in areas of upwelling. Arraial do Cabo in southeastern Brazil is the extreme point of inflection of the Brazilian coast and experiences seasonal coastal upwelling, which is a well-known phenomenon associated with east-northeast winds. This phenomenon is more frequent during the spring and summer and its high primary production supports important local fishing activities (Valentin, 1984). The aim of the present study was to describe, for the first time, the specific composition of the community of meiofaunal Harpacticoida copepods in a phytal environment situated in a region under the influence of seasonal coastal upwelling. Spatial variations between two shores with different levels of exposure to the upwelling and temporal variations between two main seasons of the year are also described. The data were used to test the hypothesis that temporal differences in the community of the Harpacticoida fauna would be greater on the shore more affected by upwelling. Arraial do Cabo depends directly on meteorological conditions. The wind pattern is responsible for the distribution of water masses, which consist of the Brazil Current, the Coastal Water and the deep cold South Atlantic Central Water (SACW). Under the influence of prevailing eastnortheast winds, the surface water flows to the open sea, allowing the deep water (approximately 300 m) to ascend. Upon entering the euphotic layer, this cold, deep water (SACW), with temperatures below 188C and nitrate concentrations about 10 mg-at l21, influences the structure of the tropical coastal ecosystem. This common pattern is inverted during the passage of cold fronts from the south. There are two main austral seasons: spring –summer, with prevailing winds favourable to upwelling, and autumn– winter, with the frequent passage of cold polar fronts and rapid changes in wind patterns unfavourable to upwelling. Despite its primarily seasonal nature, with higher frequency of occurrence in spring –summer, upwelling events are not restricted to this period and occur sporadically when the winds are favourable (Valentin, 1984; Valentin et al., 1987). In this region, upwelling is characterized by the occurrence of a local phytoplankton bloom near the coast (above 6 mg Chl-a l21) with a duration of less than 24 hours. Although the amount of chlorophyll produced is small in comparison to other regions, it is considered high for a tropical region and supports the local food chain as well as important fishery activities (Valentin & Coutinho, 1990). For the present study, two rocky shore sampling stations in Arraial do Cabo were chosen: Sonar and Pedra Vermelha (Figure 1). Sonar is located in the main area of upwelling occurrence in the outer portion of the Pontal do Atalaia. It has a high degree of declivity, with about 10 m in length from surface to bottom (maximal depth of 8 m), and high exposure to waves and winds. Pedra Vermelha is located in the inner portion of Cabo Frio Island and is a sheltered location, with a lower impact from the upwelling. This shore has a rocky slope about 25 m long from top to bottom and MATERIALS AND METHODS Study area and coastal upwelling The study area is located in Arraial do Cabo (Rio de Janeiro, southeastern Brazil) between 22858′ and 23800′ 50′′ S and 41857′ 40′′ and 42803′ W (Figure 1). The local hydrology at Fig. 1. Map of study area showing sample station ‘S’ at Sonar (located in Pontal do Atalaia) and ‘PV’ at Pedra Vermelha (located on Cabo Frio Island) in Arraial do Cabo, southeastern Brazil. harpacticoida under upwelling influence its maximal depth is also 8 m. At both shores the composition of epilithic algae turfs in depths extending to 2 m are dominated by Jania spp. and Amphiroa sp. with other species such as Centroceras clavulatum, Ceramium spp., Polysiphonia spp. and Cladophora spp. occurring with lower frequency. or autumn – winter (sporadic upwelling events). The Mann – Whitney test was used to determine differences in nitrate values between seasons. Considering that the upwelling effect over surface temperature variable affected mainly its variance, an F test was used to determine differences between seasons. The level of significance was set at P , 0.05 for all analyses. Sampling and sample processing Weekly data of surface water temperature and nitrate collected at Pedra Vermelha rocky shore over a two-month period prior to the meiofauna sampling (data provided by the Instituto de Estudos do Mar Almirante Paulo Moreira) were used to characterize upwelling seasonality in the area. At Pedra Vermelha, samples were taken in March (spring – summer, season with greater frequency of upwelling events) and September 2002 (autumn –winter season, when upwelling rarely occurs). At Sonar, samples were taken in January (spring –summer season) and June 2004 (autumn – winter season). At both rocky shores, samples were collected at high tide in depths extending to 2 m. Four replicates were performed for each rocky shore and season by SCUBA diving. The collection method consisted of scraping the substrate with the aid of acrylic bottles, gathering the epilithic algae and trapped sediment. Meiofauna samples were washed with filtered tap water through geological sieves with mesh sizes of 0.5 mm and 0.044 mm. The meiofauna retained between sieves was fixed in 4% formalin. Under stereomicroscope, fifty individuals of Copepoda Harpacticoida were selected from each replicate, when available, and placed in Eppendorf tubes with 70% alcohol. The identification of Copepoda Harpacticoida was performed under optical microscope following the taxonomic keys of Lang (1948, 1965), Huys et al. (1996) and Wells (2007) as well as publications with specific descriptions. RESULTS Coastal upwelling characterization Upwelling events were characterized by an increased concentration of nitrate in the surface water two months prior to sampling in March 2002 and January 2004 (spring and summer, respectively). There were significantly higher mean nitrate values in the periods preceding the samplings in March 2002 than those obtained in the periods preceding the samplings in September 2002 (U ¼ 2.75, P , 0.01). The same pattern was observed in 2004, with significantly higher mean nitrate values in the periods preceding the samplings in January 2004 than those obtained in the periods preceding the samplings in June 2004 (U ¼ 2.59, P , 0.01) (Figure 2A). Mean surface water temperature obtained from weekly measurements two months prior to the samplings demonstrated no apparent difference between spring – summer and autumn – winter averages (Figure 2B). However, the variance in temperature (s2 ¼ 5.6) in the period with a higher incidence of upwelling was higher than the variance (s2 ¼ 0.6) in the period with sporadic upwelling in 2002 (F7,8 ¼ 9.33, P , 0.01); no significant difference was found among variances in 2004. These data clearly indicate that the upwelling events were stronger during 2002 compared to 2004. Statistical analysis Statistical analysis was based on the methods proposed by Clarke & Warwick (1994) using the Primerw 5 program (Plymouth Routines in Multivariate Ecological Research). The Shannon – Wiener (H′ , using log2), Pielou’s evenness (J′ ), rarefaction ES(50) and species richness indices were calculated. Two-way analysis of variance (ANOVA) was used to determine significant differences for calculated ES(50) values both between shores and seasons. Bartlett’s test was used to determine the homogeneity of variances. The Bray – Curtis measure was applied to standardized data to assess the similarity between replicates, since the samples were semi-quantitative. Multi-dimensional scaling (MDS) was used to represent the similarity matrix. Two-way analysis of similarity (ANOSIM) was used to determine significant differences in the structure of the Copepoda Harpacticoida assemblage between the two rocky shores and sampling seasons. Similarity percentage (SIMPER) analysis was applied to determine which species were responsible for the similarities within rocky shores/seasons and dissimilarities between rocky shores/seasons. The t-test was used to determine whether dissimilarities between seasons are more important at Pedra Vermelha or Sonar. This test was also used to determine whether the heterogeneity between replicates of the same rocky shore is more pronounced in spring – summer (upwelling season) Fig. 2. Weekly nitrate (A) and surface water temperature (B) measurements (mean + variance) of two months prior to meiofauna samplings in 2002 and 2004 in Arraial do Cabo, southeastern Brazil. 1119 1120 v.c. sarmento et al. Copepoda Harpacticoida composition A total of 700 individuals of Harpacticoida copepods were analysed, 80% of which were identified at species level. The remaining 20% were copepodites, malformed or broken animals, for which the determination of species was not possible. Fourteen families, thirty-six genera and fifty-six different species were recorded for Arraial do Cabo (Table 1). A review of the literature (Reid, 1998; Kihara, 2003; Vasconcelos, 2008; Wandeness, 2009) indicates that nine of the 36 genera found in this study and the Amphiascus species group pacificus are reported for the first time for Brazil: Ameiropsyllus Bodin, 1979, Amphiascopsis Gurney, 1927, Psyllocamptus T. Scott, 1899, Ameiropsis Sars, 1907, Dactylopodamphiascopsis Lang, 1944, Idomene Philippi, 1843, Harpacticella Sars, 1908, Eupelte Claus, 1860 and Perissocope Brady, 1910. The dominant families in the study area were Miraciidae (21%), Parastenheliidae (15%), Ameiridae (13%) and Laophontidae (13%). Laophontidae had the highest number of species. The species of highest relative abundance were Parastenhelia spinosa (15.2%), Orthopsyllus linearis (10.3%), Ameira sp. 2 (9.8%), Amonardia sp. (5.8%), Paralaophonte congenera (5.6%) and Ectinosoma sp.1 (5.5%). Considering the rocky shores separately, the species with the highest relative abundance at Pedra Vermelha were Amonardia sp., Orthopsyllus linearis, Paralaophonte congenera, Parastenhelia spinosa, Ectinosoma sp. 1, Esola vervoorti and Harpacticus sp. 1. At Sonar, three species dominated the assemblage: Parastenhelia spinosa, Ameira sp. 2 and Orthopsyllus linearis (Figure 3). Harpacticoida community structure variation The Copepoda Harpacticoida assemblage for both rocky shores of Arraial do Cabo had high values of species richness (32 species at Pedra Vermelha, 39 at Sonar; total of 56 for the region), evenness (0.78 for Sonar, 0.81 for Pedra Vermelha; overall equitability ¼ 0.79) and diversity (4.04 for Pedra Vermelha, 4.15 for Sonar; overall diversity ¼ 4.56). ANOVA found no significant differences in ES(50) values between rocky shores (mean values being 11.88 for both shores), seasons or the interaction between these factors (P . 0.05). The mean ES(50) value was 11.9 considering all replicates. The MDS ordination analysis (Figure 4) of the Copepoda Harpacticoida fauna revealed a clear separation between the rocky shores and sampling seasons. Table 1. List of Copepoda Harpacticoida species identified in the phytal environment of Arraial do Cabo, southeastern Brazil. Order Harpacticoida Sars, 1903 Suborder Oligoarthra Lang, 1944 Family Laophontidae T. Scott, 1905 Echinolaophonte armiger (Gurney, 1927) Esola vervoorti Huys & Lee, 2000 Heterolaophonte aff. ströemi (Baird, 1834) Laophontinae sp. 1 Laophontinae sp. 2 Laophontinae sp. 3 Laophontinae sp. 4 Laophontinae sp. 5 Laophontinae sp. 6 Laophontinae sp. 7 Laophontinae sp. 8 Laophontinae sp. 9 Laophonte cornuta Philippi, 1840 Loureirophonte catharinensis Jakobi, 1953 Paralaophonte brevirostris (sensu Yeatman, 1970) Paralaophonte congenera (Sars, 1908) Family Miraciidae Dana, 1846 Amphiascus (minutus) sp. Amphiascus (pacificus) sp. 1 Amphiascus (pacificus) sp. 2 Amonardia sp. Amphiascus (varians) sp. Amphiascoides sp. Amphiascopsis cinctus (Claus, 1866) Dactylopodamphiascopsis latifolius (Sars, 1909) Diosaccus sp. Diosaccinae sp. 1 Diosaccinae sp. 2 Paramphiascella sp. Robertgurneya sp. Robertsonia sp. Family Tegastidae Sars, 1904 Parategastes sphaericus (Claus, 1863) Family Louriniidae Monard, 1927 Lourinia armata Claus, 1866 Family Dactylopusiidae Lang, 1936 Dactylopusia pectenis Pallares, 1975 Dactylopusia tisboides Claus, 1863 Dactylopusia sp. Paradactylopodia sp. Family Harpacticidae Dana, 1846 Harpacticus sp. 1 Harpacticus sp. 2 Harpacticella sp. Perissocope sp. Family Ectinosomatidae Sars, 1903 Ectinosoma sp. 1 Ectinosoma sp. 2 Sigmatidium sp. Family Peltidiidae Claus, 1860 Eupelte sp. 1 Eupelte sp. 2 Family Tisbidae Stebbing, 1910 Scutellidium sp. Tisbe sp. Family Canthocamptidae Brady, 1880 Mesochra sp. Family Orthopsyllidae Huys, 1990 Orthopsyllus linearis Claus, 1866 Family Parastenheliidae Lang, 1936 Parastenhelia spinosa Fischer, 1860 Family Pseudotachidiidae Lang, 1936 Idomene sp. Family Ameiridae Boeck, 1865 Ameira sp. 1 Ameira sp. 2 Ameiropsyllus sp. Ameiropsis sp. Psyllocamptus sp. harpacticoida under upwelling influence Fig. 3. Relative abundance (%) of dominant species of Harpacticoida at (A) Sonar (S) and (B) Pedra Vermelha (PV) according to sampling season (spring – summer and autumn – winter) in Arraial do Cabo, southeastern Brazil. The pattern observed in the MDS ordination was confirmed by ANOSIM. Significant differences in the Copepoda Harpacticoida community structure were detected between rocky shores (Rglobal ¼ 0.818, P , 0.001, number of permutations ¼ 1225) and between seasons (Rglobal ¼ 0.734, P , 0.001, number of permutations ¼ 1225). The separation between spring– summer and autumn– winter was more evident at Sonar, with an average similarity of 18.5% between replicates from the two seasons, compared to an average of 35% for Pedra Vermelha (t ¼ 5.06, P , 0.01). Although the separation between rocky shores in the MDS is clear for each season, the Bray – Curtis similarity values between replicates from the same rocky shore in the same sampling season indicate lower heterogeneity in spring– summer (average similarity ¼ 51%) when compared to autumn –winter (36%) (t ¼ 3.58, P ¼ 0.002). Table 2 displays the results of the SIMPER analysis, revealing Esola vervoorti (23.09%), Ectinosoma sp. 1 (18.92%) and Orthopsyllus linearis (12.11%) to be the species that most contributed to the similarity among replicates at Pedra Vermelha in spring –summer. On this same rocky shore during the season when upwelling is rare (autumn –winter), the species with the greatest contribution were Paralaophonte congenera (24.87%), Parastenhelia spinosa (19.22%) and Orthopsyllus linearis (10.66%). At Sonar, P. spinosa stood alone, contributing 53.94% of the average similarity among the replicates in spring –summer. In autumn –winter, Ameira sp. 2 (36.87%), Mesochra sp. (9.07%) and Amphiascoides sp. 1 (9.07%) were the most important species contributing to the similarity among replicates at Sonar. At Pedra Vermelha, Esola vervoorti (12.8%), Paralaophonte congenera (10.57%), Amonardia sp. (9.11%), Harpacticus sp. 1 (9.06%), Diosaccus sp. (7.49%) and Orthopsyllus linearis (7.38%) were the species that most contributed to the dissimilarity in the spring –summer versus autumn –winter replicates. At Sonar, the species that contributed most to the dissimilarity among replicates between seasons were Parastenhelia spinosa (20.26%), Ameira sp. 2 (9.32%), Orthopsyllus linearis (9.31%), Amphiascus (pacificus) sp. 2 (6.77%), Sigmatidium sp. (6.37%) and Mesochra sp. (4.22%) (Table 2). DISCUSSION Fig. 4. Multidimensional scaling ordination between Pedra Vermelha (diamonds) and Sonar (squares) samples in spring – summer (SS) and autumn – winter (AW) in Arraial do Cabo, southeastern Brazil. Although Copepoda Harpacticoida is a very diverse group, it has not been the subject for detailed studies on the Brazilian coast until recently. Only 68 genera had been recorded by 1998 (Reid, 1998). While the number of monographic studies, mostly not published, on the taxonomy of this group has increased in recent years (Kihara, 2003; Vasconcelos, 2008; Wandeness, 2009), 25% of the total number of genera identified in the present study are new records for the Brazilian coast. Although this paper is the first to address the Copepoda Harpacticoida assemblage on the species level in a phytal environment in Brazil, previous authors have cited the occurrence of this taxon on the genus level in a few phytal environments. Studying the meiofauna associated with the algae Sargassum cymosum on the Lázaro shore of the State of São Paulo, Curvêlo (1998) reports the genera Harpacticus, Scutellidium, Amphiascus and Ectinosoma as the main representatives of the order Harpacticoida. Jakobi (1953) found numerous individuals of the genera Harpacticus, Parastenhelia and Paralaophonte in seaweed on the Itapocorói and Porto Belo beaches in the State of Santa Catarina. The Harpacticoida fauna found in the phytal 1121 1122 v.c. sarmento et al. Table 2. Cumulative percentage (Cum.%) of species contribution to average similarity for Pedra Vermelha and Sonar rocky shores in spring–summer (SS) and autumn–winter (AW) in Arraial do Cabo, southeastern Brazil. Pedra Vermelha—SS (Average similarity 5 50.13) Esola vervoorti Ectinosoma sp. 1 Orthopsyllus linearis Amonardia sp. Cum.% 23.09 42.01 54.13 65.95 Sonar—SS (Average similarity 5 51.97) Parastenhelia spinosa Orthopsyllus linearis Amphiascus (pacificus) sp. 2 Ameira sp. 2 Cum.% 53.94 73.42 87.25 93.28 Pedra Vermelha—AW (Average similarity 5 41.67) Paralaophonte congenera Parastenhelia spinosa Orthopsyllus linearis Ectinosoma sp. 1 Cum.% 24.87 44.09 54.75 65.41 Sonar—AW (Average similarity 5 30.23) Ameira sp. 2 Mesochra sp. Amphiascoides sp. 1 Perissocope sp. Cum.% 36.87 45.94 55.01 63.16 Pedra Vermelha—SS 3 AW (Average similarity dissimilarity 5 65.03) Esola vervoorti Paralaophonte congenera Amonardia sp. Harpacticus sp. 1 Diosaccus sp. Orthopsyllus linearis Cum.% 12.80 23.38 32.49 41.55 49.04 56.41 Sonar—SS 3 AW (Average dissimilarity 5 81.52) Parastenhelia spinosa Ameira sp. 2 Orthopsyllus linearis Amphiascus (pacificus) sp. 2 Sigmatidium sp. Mesochra sp. Cum.% 20.26 29.57 38.88 45.65 52.02 56.24 environment of Arraial do Cabo was also similar in terms of genera to that reported by other authors in different parts of the world, such as Ria Deseado, Argentina (Pallares, 1968), Ashtamudi in India (Arunachalam & Nair, 1988), Cook Strait (Hicks, 1977b), Wellington (Hicks, 1986) and Island Bay (Coull & Wells, 1983) in New Zealand, Robin Hood’s Bay and St Abbs in England (Hicks, 1980), South Carolina (Coull et al., 1983) and Florida (Walters & Bell, 1994) in the USA, British Columbia in Canada (Webb & Parsons, 1992) and Port Phillip Bay in Australia (Jenkins et al., 2002). The term ‘isocommunity’ was first proposed by Thorson (1957 cited in Por, 1964) for macrofauna. The concept of isocommunity or parallel ecological communities suggests that similar substrates, although geographically separated, are inhabited or colonized by the same set of dominant genera, although the species composition may vary between sites. Por (1964) concluded that this hypothesis was even more pronounced on the meiofauna level although Gheskiere et al. (2005) studying the Nematoda community of sandy beaches suggested that the concept of parallel ecological communities is only supported for the upper beach zones. This concept was suggested by Hicks (1977b, 1980) as being applicable to Copepoda Harpacticoida communities in the phytal environment. Thus, the Harpacticoida community in Arraial do Cabo, compared to that of other locations in Brazil and regions of the world, reinforces the hypothesis of geographical parallelism or isocommunities for the phytal environment. The presence of a large number of cosmopolitan species with high relative abundance (38% of total individuals) in this study should be stressed, with Orthopsyllus linearis, Dactylopusia tisboides, Laophonte cornuta, Parategastes sphaericus, Paralaophonte brevirostris and Paralaophonte congenera standing out. These species have been cited as having considerable morphological variation or likely to belong to a complex of species with as yet unresolved taxonomy (see Wells (2007) for comments on the taxonomy of these species). The high degree of species richness and evenness for Copepoda Harpacticoida in Arraial do Cabo is reflected in the diversity value (H′ ¼ 4.56 bits), which is higher than most diversity values reported for the phytal environment (Hicks, 1980 and values calculated from Hicks, 1977b) as well as most of those found in sediments in other coastal environments (Hicks, 1980; Kihara, 2003; Gheerardyn et al., 2008). Table 3 displays the maximal diversity values (H′ log2) calculated for Copepoda Harpacticoida species in different phytal and sedimentary substrates available in the literature. The results displayed in Table 3 reveal that phytal environments have diversity indices similar to those of other coastal environments. As suggested by Huys et al. (1996), some of the highest diversity values were obtained from phytal samples including the associated sediment (Laminaria holdfast and epilithic algae). However, the other two highest values were obtained in infralittoral sediment (sandy sediments and coral fragments). The value from Kihara (2003) was obtained at a number of stations with a large variation of sediment texture (stations with the dominant sediment class ranging from silt to coarse sand), which may explain the high diversity. Gheerardyn et al. (2008) classified the Harpacticoida fauna associated with coral fragments as mainly composed of genera typically found in phytal assemblages. Similarly to what occurs with the fauna on holdfasts of macroalgae, ‘the sediment trapped by coral fragments may provide a microhabitat for sediment-dwellers, while the complex microtopography of the coral branches may be a suitable substratum for true epibenthic or even “phytal” harpacticoids’. The composition of the epilithic algae tufts at the two rocky shores of Arraial do Cabo has a great abundance of Jania sp. and Amphiroa sp., both of which are calcareous articulated algae, with erect, well-branched stems that can hold the sediment and provide a structurally complex habitat for Harpacticoida fauna. According to Hicks (1985), areas in which algal cover is performed by different species and large amounts of filamentous epiphytes have much greater spatial heterogeneity when compared with most homogeneous and monospecific areas. In these complex environments, the harpacticoida under upwelling influence Table 3. Maximal diversity values (H′ log2) calculated for Copepoda Harpacticoida assemblages on phytal and sedimentary substrates available in the literature (Sed., sediment; calc., calculated). H′ (log2) Source Phytal substrate Ceramium Cladophora Corallina fronds Corallina fronds Ecklonia fronds Enteromorpha fronds Epilithic algae Fucus Gigartina Halodule wrightii Halophila ovalis Hp. ovalis Hp. stipulacea Heterozostera tasmanica Laminaria Laminaria holdfast L. ochroleuca Palmaria Pterocladia fronds Syringodium isoetifolium Thalassia hemprichii Ulva Xiphophora fronds Zonaria fronds Zostera 2.42 2.75 2.89 3.55 4.07 1.50 4.56 2.15 2.63 3.54 2.18 4.07 3.98 3.35 2.49 5.17 3.33 1.49 3.63 4.72 3.63 2.73 2.97 1.72 1.00 Hicks (1980) Hicks (1980) Hicks (1980) calc. from Hicks (1977b) calc. from Hicks (1977b) calc. from Hicks (1977b) Present study Hicks (1980) Hicks (1980) calc. from De Troch et al. (2001) calc. from De Troch et al. (2001) calc. from Arunachalam & Nair (1988) calc. from De Troch et al. (2001) calc. from Jenkins et al. (2002) Hicks (1980) Moore (1973b) in Hicks (1980) calc. from Arroyo et al. (2006)1 Hicks (1980) calc. from Hicks (1977b) calc. from De Troch et al. (2001) calc. from De Troch et al. (2001) Hicks (1980) calc. from Hicks (1977b) calc. from Hicks (1977b) calc. from Hicks (1986) Sediment substrate Coarse sand Coral fragments Coral gravel Fine sand Mean to coarse sand Mud Sand Sand–silt Sed. adjacent to Hp. minor Sed. adjacent to Hd. uninervis Sed. adjacent to S. isoetifolium Sed. adjacent to T. hemprichii Silt to coarse sand 3.3 4.39 3.75 3.38 3.71 2.5 2.2 3.9 3.15 3.49 2.7 3.31 4.88 Hartzband & Hummon (1974) in Hicks (1980) Gheerardyn et al. (2008) Gheerardyn et al. (2008) Moore (1979) in Hicks (1980) Gheerardyn et al. (2008) Marcotte & Coull (1974) in Hicks (1980) Coull & Fleeger (1977) in Hicks (1980) Coull (1970) in Hicks (1980) calc. from De Troch et al. (2008)2 calc. from De Troch et al. (2008)2 calc. from De Troch et al. (2008)2 calc. from De Troch et al. (2008)2 Kihara (2003) 1 , values calculated for species of Thalestridae, 2, value calculated for identification to the genus level. different microhabitats favour the establishment of animals with different habits and/or morphological adaptations, which can thus coexist, foraging different niches in the seaweed. This, in turn, allows a greater distribution of food sources and a loss of competitive interactions, with the consequent establishment of rich, diverse Copepoda Harpacticoida fauna (Hicks, 1977b; Arroyo et al., 2006). In Arraial do Cabo, the high diversity values were due to the strong spatial heterogeneity and structural complexity of epilithic algae that cover the Sonar and Pedra Vermelha rocky shores. Although the rocky shores sampled in this study differ in wave exposure and slope, the indices of diversity, richness and evenness calculated separately for each shore were similar and the ANOVA was unable to detect significant differences between rocky shores or seasons for ES(50) values. These results indicate that univariate indices were not sensitive enough to detect differences between the shores or seasons. This result contrasts with that reported by Toefy et al. (2003) studying Foraminifera associated with Gelidium pristoides in False Bay, South Africa, who found greater species abundance and diversity in the exposed area, where the algae had a larger size and more trapped sediment in comparison to the fauna at the sheltered beach. However, it should be stressed that the more exposed shore in the present study had a higher total number of species than the sheltered shore. Unlike the univariate indices, the ANOSIM revealed significant differences in the Harpacticoida assemblage between the two rocky shores of Arraial do Cabo. This result contrasts with the finding described by Somerfield & Jeal (1995), who studied Halacaridae on rocky shores in Ireland and noted that, although the alga and barnacle community available as substrate was different between exposed and sheltered beaches, the association of Halacaridae was similar. The difference between the rocky shores in Arraial do Cabo may be related to physical differences, such as slope and wave exposure, or to the fact that the rocky shores are located in areas with different degrees of exposure to upwelling events. 1123 1124 v.c. sarmento et al. While the entire region of Arraial do Cabo is influenced by coastal upwelling, the effect of this phenomenon is significantly stronger at Sonar than Pedra Vermelha (Guimaraens & Coutinho, 1996). This may be explained by the fact that Sonar is located in the outer portion of the Pontal do Atalaia, which is the main area of upwelling, whereas Pedra Vermelha is located on Cabo Frio Island and is sheltered, thereby it experiences less intense effects of upwelling. In Arraial do Cabo, Guimaraens & Coutinho (1996) found distinct phytal communities during upwelling events, with a correlation between the gradient of upwelling intensity (areas with greater or lesser exposure to this phenomenon) and flora distribution. The benthic flora survey recorded elements with warm temperate affinities predominant in the summer at Sonar and algae with tropical affinities in more sheltered areas (such as Pedra Vermelha). Changes of the benthic flora affect the meiobenthic harpacticoids composition since there is a close relationship among both associations (Hicks, 1977b). Moreover, exposure differences between both rocky shores may also affect both the phytal substrate and the meiofaunal association as was earlier observed by Dommasnes (1968) and Gibbons (1988). Among the species with the highest relative abundance in the present Harpacticoida assemblage, Esola vervoorti and Diosaccus sp. occur only at Pedra Vermelha, while Amphiascus (pacificus) sp. 2, Sigmatidium sp. and Lourinia armata were unique to Sonar. Among these species, Esola vervoorti and Amphiascus (pacificus) sp. 2 occurred only in spring –summer and Sigmatidium sp. and Amphiascoides sp. 1 occurred only in samples collected in autumn –winter. On both rocky shores, there was a significant change in species composition from autumn– winter (period of sporadic upwelling) to spring– summer (period of intense upwelling). This pattern was especially pronounced on Sonar, the rocky shore more exposed to upwelling although samplings on this shore occurred during 2004, a year with weaker upwelling events according to both the surface water temperature and nitrate data that denote stronger upwelling in 2002 (Pedra Vermelha samplings). At Sonar in autumn– winter, the species Ameira sp. 2, Mesochra sp. and Amphiascoides sp. 1 together accounted for 55% of the average similarity between replicates. This percentage changed strongly in spring –summer, with Parastenhelia spinosa alone accounting for approximately 54% of average similarity. The association between upwelling enrichment with benthic communities dominated by opportunistic species has been suggested. Studying the macrofauna community on the Galician coast of northwestern Spain, Lopez-Jamar et al. (1992) found that enrichment from upwelling production reaches the seabed in pulses and, thus, benthic organisms (mainly small surface detritus-feeding polychaetes) are specialized in exploiting such irregular events. Studying the benthos in southeastern Brazil, De Léo & Pires-Vanin (2006) suggest that megabenthic community distribution in this region is closely linked to the seasonality of the SACW thermal front, with abrupt changes in species composition and the dominance of key species. This was found to be especially true for Cabo Frio, where a high dominance of top carnivore species, such as the crustacean Portunus spinicarpus and the sea-star Astropecten brasiliensis occurs during upwelling events (spring – summer 2002). In a study carried out on the sediment of the inner part of Cabo Frio Island, Fonsêca-Genevois et al. (2004) found differences in the dominance of Nematoda genera among samples from winter (September –October), dominated by Pseudosteineria, and summer (March), dominated by Viscosia. Unlike the present study and the above examples indicating a strong effect from periodicity or upwelling intensity on faunal changes, Kelaher & Castilla (2005) studied macrofaunal assemblages in coralline algal turf on seven rocky shores in northern Chile and found that there were no patterns to suggest that mesoscale variation in upwelling intensity affects either faunal assemblages or local habitat characteristics. Nevertheless, the assemblages varied significantly between shores due to differences in the complexity and heterogeneity of the habitat structure, which have greater influence on faunal assemblages in mat-like habitats on rocky shores than environmental variables associated with mesoscale coastal upwelling. The considerable variation in physicochemical conditions associated with upwelling (temperature reduction and increase in the concentration of nutrients) and the resulting unidirectional modification in phytal communities (Guimaraens & Coutinho, 1996) may favour the dominance of species adapted to such conditions, as seems be the case with Parastenhelia spinosa. Upwelling is suggested as the main factor contributing to the stronger temporal dissimilarity at Sonar, thereby confirming the initial hypothesis that the Harpacticoida assemblage at the rocky shore more affected by upwelling exhibits greater variation between spring – summer and autumn –winter samples. The peculiar hydrological characteristics of upwelling sites determine important ecological changes in benthic communities (Guimaraens & Coutinho, 1996; De Léo & Pires-Vanin, 2006; Lastra et al., 2006) and allow the occurrence of a more diverse fauna in comparison to that of areas in which this phenomenon does not occur (Le Loeuff & Cosel, 1998; Machado et al., 2005). As suggested by Witman & Smith (2003), due to these characteristics, upwelling areas warrant special attention in the planning of marine reserves to ensure the conservation of biodiversity. ACKNOWLEDGEMENTS The authors thank IEAPM for supplying the abiotic data. P.J.P. Santos acknowledges a research fellowship (305609/ 2004-1) and V.C. Sarmento an MSc scholarship from CNPq. L.M. Lage acknowledges an MSc scholarship from CAPES. Thanks are also due to Mr Richard Boike for the English revision. Thanks are also due to anonymous referees for providing comments that improved the manuscript. REFERENCES Almeida M.J. and Queiroga H. (2003) Physical forcing of onshore transport of crab megalopae in the northern Portuguese upwelling system. Estuarine, Coastal and Shelf Science 57, 1091–1102. Arroyo L.N., Maldonado M. and Walters K. (2006) Within- and between-plant distribution of harpacticoid copepods in a North Atlantic bed of Laminaria ochroleuca. Journal of the Marine Biological Association of the United Kingdom 86, 309–316. Arunachalam M. and Nair N.B. (1988) Harpacticoid copepods associated with the seagrass Halophila ovalis in the Ashtamudi Estuary, southwest coast of India. Hydrobiologia 167/168, 515–522. harpacticoida under upwelling influence Bosman A.L., Hockey P.A.R. and Siegfried W.R. (1987) The influence of coastal upwelling on the functional structure of rocky intertidal communities. Oecologia 72, 226 –232. Henriques M. and Almada V.C. (1998) Juveniles of non-resident fish found in sheltered rocky subtidal areas. Journal of Fish Biology 52, 1301–1304. Clarke K.R. and Warwick R.M. (1994) Change in marine communities: an approach to statistical analysis and interpretation. Plymouth: Plymouth Marine Laboratory. Hicks G.R.F. (1977a) Species composition and zoogeography of marine phytal harpacticoid copepods from Cook Strait, and their contribution to total phytal meiofauna. New Zealand Journal of Marine and Freshwater Research 11, 441–69. Coull B.C. (1988) Ecology of the marine meiofauna. In Higgins R.P. and Thiel H. (eds) Introduction to the study of meiofauna. Washington, DC: Smithsonian Institution Press, pp. 18–38. Coull B.C. and Wells J.B.J. (1983) Refuges from fish predation: experiments with phytal meiofauna from the New Zealand rocky intertidal. Ecology 64, 1599–1609. Coull B.C., Creed E.L., Eskin R.A., Montagna P.A., Palmer M.A. and Wells J.B.J. (1983) Phytal meiofauna from the rocky intertidal at Murrels Inlet, South Carolina. Transactions of the American Microscopical Society 102, 380 –389. Curvêlo R.R. (1998) A meiofauna vágil associada à Sargassum cymosum C. Agardh, na praia do Lázaro, Ubatuba, SP. MSc dissertation. São Paulo University, São Paulo, Brazil. De Léo F.C. and Pires-Vanin A.M.S. (2006) Benthic megafauna communities under the influence of the South Atlantic Central Water intrusion onto the Brazilian SE shelf: a comparison between an upwelling and a non-upwelling ecosystem. Journal of Marine Systems 60, 268 – 284. De Troch M., Fiers F. and Vincx M. (2001) Alpha and beta diversity of harpacticoid copepods in a tropical seagrass bed: the relation between diversity and species’ range size distribution. Marine Ecology Progress Series 215, 225–236. Hicks G.R.F. (1977b) Species associations and seasonal population densities of marine phytal harpacticoid copepods from Cook Strait. New Zealand Journal of Marine and Freshwater Research 11, 621–643. Hicks G.R.F. (1980) Structure of phytal harpacticoid copepod assemblages and the influence of habitat complexity and turbidity. Journal of Experimental Marine Biology and Ecology 44, 157 –192. Hicks G.R.F. (1985) Meiofauna associated with rock shore algae. In Moore P.G. and Seed R. (eds) The ecology of rocky coasts. New York: Columbia University Press, pp. 36–56. Hicks G.R.F. (1986) Distribution and behaviour of meiofaunal copepods inside and outside seagrass beds. Marine Ecology Progress Series 31, 159–170. Huys R., Gee J.M., More C.G. and Hammond R. (1996) Marine and brackish water Harpacticoid Copepods. Part 1: keys and notes for identification of the species. In Barnes R.S.K. and Crothers J.H. (eds) Synopses of the British fauna (New Series) No. 51. Shrewsbury: Field Studies Council, pp.1–352. Jakobi H. (1953) Novos Laophontidae (Copepoda–Crustacea) da costa brasileira. Dusenia 4, 47–60. Jenkins G.P., Walker-Smith G.K. and Hamer P.A. (2002) Elements of habitat complexity that influence harpacticoid copepods associated with seagrass beds in a temperate bay. Oecologia 131, 598–605. De Troch M., Melgo-Ebarle J.L., Angsinco-Jimenez L., Gheerardyn H. and Vincx M. (2008) Diversity and habitat selectivity of harpacticoid copepods from seagrass beds in Pujada Bay, the Philippines. Journal of the Marine Biological Association of the United Kingdom 88, 515–526. Kelaher B.P. and Castilla J.C. (2005) Habitat characteristics influence macrofaunal communities in coralline turf more than mesoscale coastal upwelling on the coast of Northern Chile. Estuarine, Coastal and Shelf Science 63, 155–165. Dommasnes A. (1968) Variations in the meiofauna of Corallina officinalis L. with wave exposure. Sarsia 34, 117–124. Kihara T.C. (2003) Diversidade dos copépodes harpacticóides da meiofauna marinha do litoral norte do estado de São Paulo. PhD thesis. São Paulo University, São Paulo, Brazil. Fonsêca-Genevois V., Santos G.A.P., Castro F.J.V., Botelho A.P., Almeida T.C.M. and Coutinho R. (2004). Biodiversity of marine nematodes from an atypical tropical coast area affected by upwelling (Rio de Janeiro, Brazil). Meiofauna Marina 13, 37–44. Lang K. (1948) Monographie der Harpacticiden: I: 1– 896; II: 897–1682. Lund, Stockholm, Nodiska Bökhandeln Håkan Ohlsson Booksellers, 1682 pp. Gibbons M.J. (1988) The impact of sediment accumulations, relative habitat complexity and elevation on rocky shore meiofauna. Journal of Experimental Marine Biology and Ecology 122, 225–241. Lang K. (1965) Copepoda Harpacticoidea from the Californian Pacific Coast. Kungliga Svenska Vetenskapsademiens. Handlinger, Fjrde Serien 10, 1 –560. Giere O. (2009) Meiobenthology: the microscopic motile fauna of aquatic sediments. 2nd edition. Berlin: Springer-Verlag. Lastra M., La Huz R., Sanchez-Mata A.G., Rodil I.F., Aerts K., Beloso S. and Lopez J. (2006) Ecology of exposed sandy beaches in northern Spain: environmental factors controlling macrofauna communities. Journal of Sea Research 55, 128–140. Gheerardyn H., De Troch M., Ndaro S.G.M., Raes M., Vincx M. and Vanreusel A. (2008) Community structure and microhabitat preferences of harpacticoid copepods in a tropical reef lagoon (Zanzibar Island, Tanzania). Journal of the Marine Biological Association of the United Kingdom 88, 747 –758. Le Loeuff P. and Cosel R. (1998) Biodiversity patterns of the marine benthic fauna on the Atlantic coast of tropical Africa in relation to hydroclimatic conditions and paleogeographic events. Acta Oecologica 19, 309–321. Gheskiere T., Vincx M., Urban-Malinga B., Rossano C., Scapini F. and Degraer S. (2005) Nematodes from wave-dominated sandy beaches: diversity, zonation patterns and testing of the isocommunities concept. Estuarine, Coastal and Shelf Science 62, 365–375. López-Jamar E., Cal R.M., González G., Hanson R.B., Rey J., Santiago G. and Tenore K.R. (1992) Upwelling and outwelling effects on the benthic regime of the continental shelf off Galicia, NW Spain. Journal of Marine Research 50, 465–488. Guimaraens M.A. and Coutinho R. (1996) Spatial and temporal variation of benthic marine algae at Cabo Frio upwelling region, Rio de Janeiro, Brasil. Aquatic Botany 52, 283 –299. Ma H. (2005) Spatial and temporal variation in surfclam (Spisula solidissima) larval supply and settlement on the New Jersey inner shelf during summer upwelling and downwelling. Estuarine, Coastal and Shelf Science 62, 41–53. Hall M.O. and Bell S.S. (1993) Meiofauna on the seagrass Thalassia testudinum: population characteristics of harpacticoid copepods and associations with algal epiphytes. Marine Biology 116, 137–146. Machado C.P., Coimbra J.C. and Carreño A.L. (2005) The ecological and zoogeographical significance of the sub-recent Ostracoda off 1125 1126 v.c. sarmento et al. Cabo Frio, Rio de Janeiro State, Brazil. Marine Micropaleontology 55, 235–253. Menge B.A. (2000) Top-down and bottom-up community regulation in marine rocky intertidal habitats. Journal of Experimental Marine Biology and Ecology 250, 257–289. Ormond R.F.G. and Banaimoon S.A. (1994) Ecology of intertidal macroalgal assemblages on the Hadramout coast of Southern Yemen, an area of seasonal upwelling. Marine Ecology Progress Series 105, 105–120. Pallares R.E. (1968) Copépodos marinos de la Rı́a Deseado (Santa Cruz, Argentina). I. Physis 27, 1–125. Valentin J.L., André D.L. and Jacob S.A. (1987) Hydrobiology in the Cabo Frio (Brazil) upwelling: two dimensional structure and variability during a wind cycle. Continental Shelf Research 7, 77–88. Valentin J.L. and Coutinho R. (1990) Modelling maximum chlorophyll in the Cabo Frio (Brazil) upwelling: a preliminary approach. Ecological Modelling 52, 103 –113. Vasconcelos D.M. (2008) Distribuição dos Copepoda Harpacticoida da meiofauna em área de talude no litoral de Sergipe, Brasil. PhD thesis. Federal University of Pernambuco, Recife, Brazil. Por F.D. (1964) A study of the Levantine and Pontic Harpacticoida (Crustacea, Copepoda). Zoologische Verhandelingen 64, 1 –128. Walters K. and Bell S.S. (1994) Significance of copepod emergence to benthic, pelagic, and phytal linkages in a subtidal seagrass bed. Marine Ecology Progress Series 108, 237–249. Reid J.W. (1998) Maxillopoda –Copepoda Harpacticoida. In Young P.S. (ed.) Catalogue of Crustacea of Brazil. Rio de Janeiro: Série Livros 6. Museu Nacional, pp. 75–127. Wandeness A.P. (2009) Ecologia e taxonomia da associação de Copepoda Harpacticoida no talude da bacia de Campos, RJ, Brasil. PhD thesis. Federal University of Pernambuco, Recife, Brazil. Schiel D.R. (2004) The structure and replenishment of rocky shore intertidal communities and biogeographic comparisons. Journal of Experimental Marine Biology and Ecology 300, 309–342. Webb D.G. and Parsons T.R. (1992) Winter–spring recruitment patterns of epiphytic harpacticoid copepods in a temperate-zone seagrass bed. Marine Ecology Progress Series 82, 151–162. Soltwedel T. (1997) Meiobenthos distribution pattern in the tropical East Atlantic: indication for fractionated sedimentation of organic matter to the sea floor? Marine Biology 129, 747–756. Wells J.B.J. (2007) An annotated checklist and keys to the species of Copepoda Harpacticoida (Crustacea). Zootaxa 1568, 1 –872. Somerfield P.J. and Jeal F. (1995) Vertical-distribution and substratum association of Halacaridae (Acari, Prostigmata) on sheltered and exposed Irish shores. Journal of Natural History 29, 909 –917. Thompson R.C., Crowe T.P. and Hawkins S.J. (2002) Rocky intertidal communities: past environmental changes, present status and predictions for the next 25 years. Environmental Conservation 29, 168–191. Toefy R., McMillan I.K. and Gibbons M.J. (2003) The effect of wave exposure on the foraminifera of Gelidium pristoides. Journal of the Marine Biological Association of the United Kingdom 83, 705–710. Valentin J.L. (1984) Analyses des paramètres hidrobiologiques dans la remontée de Cabo Frio (Brésil). Marine Biology 82, 259–276. and Witman J.D. and Smith F. (2003) Rapid community change at a tropical upwelling site in the Galápagos Marine Reserve. Biodiversity and Conservation 12, 25–45. Correspondence should be addressed to: V.C. Sarmento Universidade Federal de Pernambuco Centro de Ciências Biológicas Departamento de Zoologia Avenida Prof. Morais Rêgo s/n, 50670-420, Recife, PE, Brazil email: visnu.ubi@gmail.com