Abstract

For the first time, and in a large spatial scale, the influence of ecological properties on the aquatic bird community of black water lakes in Brazilian Amazonia is evaluated. Bird surveys were conducted in 45 lakes. A total of 3626 individuals in 48 bird species were recorded; of these, 31 are aquatic, and 18 of these are primarily piscivorous. Bird richness and abundance were not significantly related to lake shape and productivity but were influenced by hydrological period (low versus high), water depth, transparency, lake isolation, and habitat richness. Matrices of bird species by lake were subjected to multivariate analyses (NMDS) to evaluate how these parameters influence bird community. The variation in bird species composition was positively correlated to lake depth and isolation and negatively correlated to water transparency and habitat richness. The results indicate that period, lake physical characteristics (depth, water transparency), isolation, and habitat richness are determinants of aquatic bird community composition in the black water lake systems of Amazonia.

1. Introduction

Ecological studies on the composition and structure of biological communities are critical to understanding the interactions between the species and populations that comprise them, to explain local and regional biodiversity [1, 2], and are essential to guiding management actions and conservation policies.

Tropical wetlands provide habitats for a wide variety of plants, fishes, birds, and mammals. Aquatic environments (rivers, lakes, and wetlands) represent 6% of Brazilian Amazonia, and over 100,000 km2 of lakes and swamps [3].

Aquatic birds are ubiquitous components of Amazonian freshwater systems. However, their role in the ecological dynamics of these systems has often been overlooked, as in other rainforest areas in the Neotropics (but see [46]). In order to evaluate the impacts of human activities on freshwater systems and improve their conservation value, the status and trends of freshwater biodiversity need to be monitored [7].

Because many aquatic birds are top predators in the aquatic environments, they can affect the distribution and abundance of fish [8], and presumably other biota. Both in tropical and temperate wetlands, aquatic bird communities have been demonstrated to be influenced by their environment (lake features such as their degree of isolation, depth, physical-chemical properties, and habitat availability) or/and by their own interaction and population dynamics (see below). For instance, in southern Brazil in fragmented wetlands, it has been recently found that some aquatic birds are more abundant in small and isolated lakes [6]. The authors defend the idea that the physicochemical conditions, type of landscape, together with food resources, and habitat might influence ecological gradients in the bird communities and therefore deserve further analysis.

Bird assemblage structure and composition can also be influenced by local ecological interactions such as predation and competition, and by large-scale environmental factors, such as the degree of habitat isolation and structure [1, 5, 6]. Considering only internal mechanisms, population structure and consequently community composition can result from a link in processes operating at the local and regional scales, which suggests that the locally most abundant species tend to have larger or extended distributions, and the replacement of individuals can be a result of immigration from the regional species pool [9]. For example, a positive relationship between the locally most abundant species and the range of their distributions across a given wetland area can indicate population dispersal as predicted by the neutral theory of biodiversity [9]. Nevertheless, weak correlations between spatial variation in bird community composition and geographical distances (distances between lakes) can be taken as evidence against the distance decay predicted by Hubbell’s theory [9]. Although the present study is within this theoretical framework and presents results on these relationships, the main purpose is not to test the theory.

Similar to those of other neotropical wetlands (the Everglades, Llanos and Pantanal), the aquatic bird communities in the Amazonian lakes are composed by grebes, cormorants, anhingas, herons, egrets, ibises, storks, waterfowl such as ducks, kingfishers, skimmers, terns, hawks and kites, rails, and jacanas [5, 6]. Some species seem to privilege lakes with shallow water, with shrubs and grasses along the shore, to forage and nest, while others use deep open waters with abundant trees in the margins. Because lake features (isolation from other lakes, water transparency, productivity) and habitat richness also vary spatially, bird species composition may also differ through many lakes over a large area, probably following gradients of spatial variation in ecological and environmental factors [5].

Only two studies of birds associated with Amazonian aquatic environments analyzed how bird community composition varies with water body characteristics and ecological gradients [4, 5]. These and other good descriptive studies [913] have dealt with bird communities of white water systems (várzea), which supposedly would have more aquatic bird species than black water (igapó) wetland systems, because they are more productive. However, investigation of other local ecological factors such as lake features (i.e., water transparency, depth, productivity, habitat richness at lake margins) and geographical position (lake isolation in relation to other lakes and to river margins) has never been undertaken for the entire Amazon region.

In a previous study, it was demonstrated that variations in white water lake shape produced ecological gradients, resulting in spatial variation in aquatic bird communities [5]. The study was made during the period of highest waters in Amazonian lakes (July). The study indicates that bird richness and abundance tend to decrease with the increase in lake shape complexity, certainly influencing the variation that bird community composition also undergoes with lake shape. The authors were surprised with these results, because they believed that the more complex the shape of a lake, the higher the diversity of microhabitats, and in consequence, the higher the species richness and abundance. However, aquatic birds are very territorial during the high water season, when most food resources are dispersed throughout the landscape. Their territories tend to be linear, located in the forest-water ecotone. In elongated lakes, opposite margins are usually closer than in roundish ones, and aquatic birds tend to defend both margins instead of only one. This inter- and intraspecific competition for space (to gain access to food resources) may make habitats more difficult to colonize by dispersers coming from other lakes. This may help explain the differences in waterbird species richness and composition between long and round lakes. Although they have not analyzed the spatial distribution of each individual species, during their surveys they noticed an “even distribution” of several egrets, herons, kingfishers, and allies in the margins of long water bodies.

My expectation is that lake shape in Anavilhanas will not influence bird communities, since most of them tend to be elliptical and not so variable (rounded to amoeboid or even dendritic) such as those that have influenced and form white water systems [5]. However, because the variation on habitat richness at lake shore may be independent of lake shape, I expect it will influence bird community.

Due to high seasonal precipitation over the Amazon basin, associated with very low inclination of the terrain, the wetlands on the margins of the large rivers are annually flooded up to 12 m depth, and this scenario may last for months [14]. Therefore, the várzeas flooded by white waters, receiving annually high amounts of sediments rich in nutrients and supporting high productivity, are historically the most heavily occupied and with the highest human densities [15]. On the other hand, the black water areas of igapó are inundated by rivers that originate in geologically older areas and therefore are acidic, with low sediment loads, and poor in nutrients (see below) resulting in low human occupation.

Together, várzeas and igapós cover about 400,000 km2 (up to 300,000 for várzea and 100,000 for igapó) among the 1.5 million km2 of wetland area estimated for the Amazonian region [16]. Both várzea and igapó have a high diversity of herbaceous plants [17] and are the richest flooded forests in tree species in the planet [14].

Because black water systems have lower productivity, it is expected that bird richness and assemblages differ from those that have been described previously for white waters, showing a gradient in community structure and composition from black to neighboring white water systems due to lake shape, and probably due to lake to lake spatial differences in habitat richness and water productivity [5]. What is not so clear yet is how bird communities respond to variation in lakes within black water systems themselves, principally in those areas with similar lake shape and no influence of white water wetlands. Although black waters are the most acidic in Amazonia and are much less productive than white waters, they still seem to have similar and high biodiversity, as indicated by their high fish and aquatic and forest bird richness [4, 18, 19]. However, the investigation on factors influencing aquatic bird communities in black water lakes needs to be made to confirm these expectations.

The black water lakes of Anavilhanas are located in islands in the second largest fresh water archipelago in the planet; the first is located upwards in the same Rio Negro, a typical representative of Amazonian black waters, and floods huge areas of forest, providing the environment for endemic and rare taxa, including recently described fish families and lizard genera.

Although more recently the knowledge has substantially increased on biology and ecology of many aquatic bird species [5, 6, 1013], studies are lacking on the function of species in black water systems and how their communities are influenced by local variation in the physical and biological features observed in black water lakes.

This is the first ecological study in the Amazon on the aquatic bird community of a system composed exclusively of black water lakes. This study is an overview of the system and investigates how richness, abundance, and species composition of aquatic bird communities differ over a large spatial scale among lakes and is locally influenced by lake features and seasonality.

The black water wetlands of the Anavilhanas Archipelago, central Brazilian Amazonia, are composed of about 70 lakes differing in size and several features, distributed over the landscape (about 140 km long and up to 20 km wide in the middle of the archipelago). Lake systems are useful to evaluate whether ecological factors (as possible underlying mechanisms) and differences in spatial scale influence aquatic bird community. Relating the environmental continuum and lake characteristics of the Anavilhanas to attributes of bird species, such as variation in presence/absence and abundance, and/or species composition across the archipelago, may help understand how lake ecological components influence their community.

In this study, I selected eight variables, including lake shape, depth, productivity (nitrogen and phosphorus concentration), water transparency, isolation (distance to river margins and distance to nearest lake), and habitat richness, that may be important because (i) many bird species (ducks, egrets, herons, storks, and allies) rely on lakes with complex shapes, transparent, and shallow waters [3, 5, 6, 13] for foraging, and lakes with more habitats can offer a more complex structure, sustaining more species [6], and may provide more area for feeding and nesting, while providing spots for sexual display for some species; for example, males of Muscovy duck use patches of grass near the shore as arenas where they fight for females [R. Cintra unpublished observations]; (ii) lakes bearing more productive waters may provide higher food availability [2023]; the abundance of small vertebrates such as small fishes and large invertebrates such as large arthropods, crabs, and snails (prime avian food resources) is directly related to water productivity [20, 23]; (iii) degree of isolation is an important component of wetlands; for example, more isolated lakes may increase breeding success because they are located further away from access of terrestrial predators [6].

Although in recent years a great body of knowledge has been gathered on general aspects of the biology and behaviour of many tropical species, little is known about the behaviour and distribution of species in wetland ecosystems, and how a significant part of the aquatic bird community responds to the local variation observed in the physical and biological characteristics of tropical lakes, mainly of black waters. Within- and among-year climatic seasonality can also be determinant in the spatial distribution of bird species richness, so I also evaluated its effects on the bird community.

Since most aquatic birds are good colonizers, opportunistic, and widespread in the Amazonian basin [24], the main prediction of this study is that differences in species composition in relation to lake features is relatively similar among black water lakes throughout the Anavilhanas Archipelago. However, because lakes vary in their features, and because bird species differ in their habitat preferences and ability to occupy the complex environment, some direct relationships are expected, for example, ecological turnover in bird species composition with concomitant variation in lake habitat richness. Because the locally abundant species seem to have ample distribution, I also predict that aquatic bird density in each lake will correlate strongly with the frequency of lakes used by the species.

I verified the effect of the eight variables on bird community using data based on speedboat transect surveys. Specifically, I tested the hypothesis that spatial variation in lake characteristics should influence the richness (total number of bird species/lake), abundance, and variation in bird species’ composition (i.e., how bird community undergoes species substitution across lakes), in a tropical black water wetland system.

The main questions I am asking are the following. (1) Does bird community composition change seasonally (low versus high water period) and among years in the black water lakes of Anavilhanas archipelago? (2) Does bird richness, abundance and species composition vary along a lake gradient that includes lake shape, water transparency, depth, productivity, isolation, and habitat richness?

2. Material and Methods

2.1. Study Area

Bird surveys were conducted on the black water lakes of Anavilhanas Archipelago, state of Amazonas, Brazil (Figures 1 and 2, 2°3′ S and 3°2′ S, 60°22′ W, and 61°12′ W), a National Park administered by Instituto Chico Mendes (ICMBio-MMA, Brazil), located 40 km northwest of Manaus, the Amazonas state capital. Anavilhanas National Park has an area of 350,000 ha, including the entire archipelago, in the Rio Negro, and a continuous and extensive area covered with terra firme (upland) rainforest on its left margin (Figure 1). The Archipelago itself is mainly covered by 100,000 ha of igapó, a forest seasonally inundated by black water relatively poor in nutrients [2527]. The mean annual temperature is 25°C. Rainfall is about 2500 mm/year, most of it between November and May, with a short dry season between June and October [4].

There are three main kinds of Amazonian aquatic systems according to water type: white-, black-, and clear-water systems [28, 29]. White waters originate in the Andes, are turbid and muddy due to high loads of inorganic particles and suspended minerals, with transparency no greater than 20 cm, and more productive because of high nutrient concentration; their pH is almost neutral (6.6 to 7) due to dissolved bicarbonates and high conductivity (above 60 μS/cm). Clear waters originate in the pre-Cambrian crystalline central Brazilian plateau and Guiana shield. They are greenish and transparent, they can have low to high concentration of dissolved and suspended minerals, water transparency is very high (5-6 m), their pH is acid to neutral (5 to 6), and their conductivity (6 to 50 μS/cm) and productivity are low. Black waters originate in washed sands and podzols in Central Amazonian lowlands, in crystalline shields and sedimentary basins. They drain from Campina vegetation (with incomplete decomposition of organic matter) in very porous soils, facilitating the transport of humic acids to the main rivers and lakes. Black waters are acidic, with high concentration of humic acids in a colloidal form and fulvic acids (contents, resp., ten and three times higher than in white waters); consequently, they are dark-brown coffee in color, and their transparency is up to 1 m. The low nutrient concentration probably makes them the poorest in nutrients and electrolytes in Amazonia, and accordingly, conductivity is very low (8 μS/cm) [21, 22, 29, 30].

The lakes of Anavilhanas are delimited by forest, usually in the form of levees. Their sediments originate in the Guianan Shield and in Central Amazonia uplands of the late Cretaceous and the early Tertiary/Quaternary periods [29]. The water is poor in nutrients, with high transparency and low pH (<5). The dark water color is due to dissolved humic and fulvic acids from the forest litter during processes of decomposition [29, 30].

Filoso and Williams [30] provided a characterization of the lakes’ hydrology; for this study, some more data were collected as a complement. The lakes receive water from the Rio Negro, and those located close to river margins, also from terra firme forest streams, locally known as igarapés. The frequency and duration of the floods, depending on the area, vary from 75 to 250 days [31]. They have ecological importance and great impact on the local biota, determining large spatial-temporal variation in the aquatic and terrestrial phases of plant and animal communities [21, 23, 31]. The Rio Negro water level in Anavilhanas varies between 8 and 12 m, with the highest flooding in June-July and the lowest waters in October-November. At the time when the main bird surveys in this study were conducted, lake depth varied from 0.88 to 7 m. In Anavilhanas, during the period of the highest waters, the igapó forest is submerged, with water level fluctuating between 1 and 20 m, depending on the distance to the Rio Negro [18]. During the low water period, macrophytes abound, especially Utricullaria foliosa, Salvinia auriculata, Pistia stratioites, free floating Paspalum repens, Eichhornia crassipes, and Oryza perennis. More information on typical plant species of the area are presented below and in [17]. The vegetation along lake shores allows the coexistence of terrestrial forest and aquatic birds. Long-distance migrant birds from the northern hemisphere annually find along the white sands and muddy beaches of the Anavilhanas lakes a temporary feeding area and refuge [32].

2.2. Bird Surveys

Bird surveys were conducted twice a day, between 6.30 and 10.00 h, and between 16.00 and 18.00 h, in 45 lakes, from 14 November to 17 December 2007. In this period of the annual inundation cycle, water level is low, and starting to rise. At this time, the forest floor of the igapó (black water flooded forest) is not submerged yet, and most Nearctic migrant birds (shorebirds, and allies) are no longer present in the area. Although wading and probably most other aquatic birds often move from a foraging spot to another as the water table drops in the dry season, in response to water recession and prey availability, they were sedentary in the area in this period, using the same lake for weeks (R. Cintra unpublished data). The large size and/or conspicuous foraging behavior of most bird species results in high detection probability. Therefore, for this study, the data were not adjusted for false absences [33]. Water level was typical for the area and the period of the year [29]. For the purpose of investigating the effects of lake features on bird communities (see Section 2.4), the surveys were conducted only once in each of the 45 lakes, in a very short period, but covered a large area; hence, a snapshot of the aquatic bird community in the study area could be obtained. This approach is convenient in studies where birds are not captured and marked and avoids recording the same bird twice, overestimating count numbers.

To investigate potential seasonal effects on the bird community, two additional surveys were done at maximum water level (July 2007 and July 2008) and two at low water level (November-December 2007 and November-December 2009). Because during high flooding bird abundance was very low in both years, and during the shallow waters of 2009 the majority of the lakes were inaccessible, only a subsample of 16 small- and medium-sized lakes was surveyed, between 6.30 and 10.00 h. Also, because bird abundance during the high-water periods is too low, for this part of the study only species presence/absence data was used in the analyses to compare seasonal species richness and community composition between the high- and low-water periods.

All surveys were done aboard a 30 hp outboard speedboat driven by a field assistant, most of the time at a speed of 15–20 km/h. Birds were observed with 10 × 50 and 8–20 × 50 binoculars, sometimes using a portable counter.

The sample unit was the lake, which was completely circumnavigated, travelling at a distance of about 10–15 m from the shore. A Garmin 76 GPS was used to record geographical coordinates in order to obtain lake positions (Table 2); because lakes can have several names, they were identified by numbers instead of names.

Lakes of several sizes were sampled (see Figures 1, 2, Table 2). The influence of lake size was already considered when I estimated and used the lake shape (which includes lake area in the formula). Because of this, I did not evaluate the effects of lake area on birds using rarefaction curves which would be an obvious thing to do, but when various bird surveys are made per lake, which is not the case. However, rarefaction curves were constructed and presented but for the entire lake system.

Although it would not be necessary because none of lakes had an exaggeratedly higher number of birds than the rest (e.g., hundreds of species or thousands of individuals, see Table 2), I followed an earlier protocol of data standardization [4, 5] in order to be able to make regional comparisons with an ongoing study. Therefore, to make bird species richness and abundance, and community composition, comparable among lakes of different sizes, and with other studies, the data were standardized by dividing the number of individuals of each species by the sizes (areas) of the lakes in which the species occurred. The densities obtained were then used in the quantitative statistical analyses. Rarefaction curves of bird richness in relation to number of lakes sampled and to bird number were constructed using “R” program [34] and are also presented.

The density approach assumes that birds potentially could use the entire lake area. However, kingfishers and other birds commonly use trees along lake margins; herons and egrets (Rufescent Tiger-Heron Tigrisoma lineatum; Striated Heron Butorides striatus; Great Egret Ardea alba), rails (Gray-necked Wood-Rail Aramides cajanea), and Sungrebe Heliornis fulica, use floating meadows along the margins and in the central part of the lakes; cormorants (Neotropic Cormorant Phalacrocorax brasilianus), Anhinga Anhinga anhinga, piscivorous hawks (e.g., Black-collared Hawk Busarellus nigricollis and Osprey Pandion haliaetus) are commonly observed perched in dead trees, wherever they are located, and fishing over the entire lake area; and terns (Large-billed Tern Phaetusa simplex and Yellow-billed Tern Sternula superciliaris) use the center and margins of lakes and rivers. Any bird-counting method in Amazonian aquatic systems needs to consider differences in species distribution in relation to habitats. Since the surveys were made in the central and marginal parts of lakes alike, the influence of differences in their habitat use (birds that use trees along lake margins and those that use floating meadows along the margins and in the central part of the lakes) in local species distribution on the results is likely negligible.

The marginal habitat data (see below) were also standardized for the same reasons as above by dividing the number of times each habitat was recorded at lake margins by the perimeter of the lakes in which the habitat occurred. Information on biology of the bird species is available in Del Hoyo et al. [35].

2.3. Lake Characteristics Sampling

Lake depth was recorded at the center of each lake using a 5 kg weight attached to a 20 m rope graduated in meters, with a precision of 0.25 m. Water transparency was measured with a secchi disc. At the same spot, at a depth of 20 cm, a 50 ml water sample was collected in a sterile plastic vial to analyze water productivity. Total nitrogen (N) and total phosphorus (P) concentrations (μmol l−1) in water samples were determined through simultaneous analyses [36] conducted at INPA’s laboratory. A recent study in 20 lakes around Manaus (including three in Anavilhanas) showed significant correlations between both and chlorophyll-a [22], so they were used as indices of lake productivity in this study. Several experienced Amazonian limnologists informed that variation in lake water depth, water transparency, and productivity is very low, and one sample in the centre would represent adequately each lake. Therefore, lakes were sampled for these parameters only once.

The aquatic bird community composition was also compared among lakes of different shapes. The Global Mapper Program [37] was used to estimate lake dimensions (area and perimeter) on a satellite image. To estimate lake shape, or deviation from a circle (a circular lake assuming an shape index, SI = 1.0, and all other shapes assuming higher values), the shape index (SI) of Patton [38] was used as follows: , where SI = lake shape index; = perimeter of the lake in km; = 3.1416; A = area of the lake in km2.

Microhabitat surveys (lake marginal vegetation), hereafter called habitat, were conducted concomitantly with the bird surveys. With the speedboat in motion every 10 s, an assistant (Jessica Cancelli), with extensive experience in the region, looked at the lake margin and recorded habitat type on a form containing habitat names. As the lake littoral strip can vary substantially in area, only the habitats located in the immediate vicinity of water were considered; these are the ones most often used by birds (personal observations). Habitats were named based on physiognomy, complemented with the presence of some dominant plant species. Eleven habitats along lake margins were identified: (1) stands of wild rice (Oryza perennis), sometimes partially submerged and forming small “peninsulas” in lakes; (2) stands of shrubs dominated by araçá-do-igapó (Eugenia inundata) and camu-camu (Myrciaria dubia); (3) stands of carauaçu (Symmeria panicullata); (4) stands of the tree macacarecuia (Eichweilera tenuifolia); (5) groups of snags (standing dead trees); (6) groups of dead logs; (7) secondary forest dominated by the embaúba tree (Cecropia spp.); (8) secondary forests dominated by jauari palms (Astrocaryum jauari); (9) igapó forest (with typical and some of the abundant flood-tolerant trees Alchornea schomburgkiana, Leopoldinia pulchra, arapari (Macrolobium acacifolium), macucu-do-igapó (Aldina latifolia), Mabea nitida, and Hevea spruceana, ipê (Tabebuia barbata) [31, 39]; (10) muddy beaches with no vegetation; (11) muddy beaches covered with scattered patches of grass (Oryza perennis) and rush (Cyperaceae). After standardization (see bird survey data above), habitat richness was used in the analyses to evaluate its influence on the bird community. Piedade and Junk [17, 40] presented information on herbaceous plants of the Amazonian floodplain, some of which occur in Anavilhanas.

2.4. Statistical Analysis

Aquatic bird richness, abundance, and community composition were compared among 45 lakes within the Anavilhanas Archipelago. Of the total community (birds associated with aquatic environments), only a group of 31 species was subjected to statistical analyses: those always associated with the aquatic environment, hereafter called the aquatic bird community, which includes the nine species considered piscivorous by Petermann [11]. They are marked as A (aquatic) and P (piscivorous) in the “Aq./Pisc.” column of Table 2, which lists all the birds associated with aquatic environments in the study area.

Mantel tests, available in PC-ORD [41], were run to verify spatial autocorrelation of the quantitative aquatic bird community (matrix of species abundance by lakes) and lake location (matrix of Euclidian distances with geographical coordinates by sites (lakes), using latitude and longitude metric (UTM) values. The method chosen for the test was Mantel’s asymptotic approximation [41].

To investigate whether the most abundant species were also the most widely distributed in the lakes, a simple linear statistical model was used to verify the relationship between bird abundance (mean bird density, which was log-transformed before analysis), and bird distribution in the Archipelago (number of lakes in which they occur). In this specific analysis, only species occurring in at least five lakes (21 of the 31 bird species) were included.

Multiple linear statistical models were used to look for relationships between qualitative and quantitative bird community composition (NMDS axes, see below) and the independent variables of lake shape, degree of isolation (distances from the centre of a given lake to the centre of the nearest neighbouring lake and to river margins), depth, water transparency, habitat richness, and lake productivity. Simple linear models were used to verify relationships between bird richness and abundance and the same independent variables. Before running the analyses, a Pearson correlation matrix was used to verify the significance of correlations among the independent variables. Only those that were not significantly correlated were considered in the same statistical models.

Qualitative matrices with species presence/absence, and quantitative matrices with species abundance, by lake, were constructed for the aquatic bird communities. In the matrices, the quantitative values (number of birds recorded) for each species were divided by lake area to obtain bird densities, which were then used as dependent variables in models of multiple linear statistical models to test the null hypothesis that bird community composition is similar among lakes.

To compare bird species composition among lakes, all quantitative and qualitative data matrices constructed from the bird surveys were subjected to nonmetric multidimensional scaling (NMDS), an multivariate ordination analysis available in the program PC-ORD [41, 42]. An advantage of using NMDS is that it is based on ranked distances and tends to linearize the relationship between environmental distance and biological distance [43], a procedure that alleviates the “zero truncation problem,” which plagues all ordinations of heterogeneous datasets [41, 44]. According to these authors, NMDS has performed well with both simulated gradients and field data.

The distributions of species densities generally do not conform to the assumptions of multivariate inferential analyses such as multivariate analysis of variance [45]. Therefore, NMDS was performed on the data matrices to obtain linear, orthogonal variables (axes) describing the bird community composition that met the assumptions of the multivariate inferential analyses [46].

Throughout the study, I used two NMDS axes for quantitative data and three axes for qualitative data as dependent variables. For community data from 45 lakes, two NMDS axes captured most of the variance in the original variables for quantitative data in the bird species matrix (cumulative proportion of total variance, CPV = 0.83), whereas three axes were enough for the presence/absence data matrix (CPV = 0.82). Similarly, for the seasonal data, three axes captured most of the variance for presence/absence data (CPV = 0.85).

The Bray-Curtis index was used to describe dissimilarity between lakes. This index is available in PC-ORD program [41], and was used on abundance and presence/absence data. The use of quantitative data may lead to patterns influenced by the higher contribution of the most abundant species to the ordination; however, for qualitative ordination, they contribute less to dissimilarities in the bird communities among lakes because qualitative ordination tends to capture mainly the patterns of the rarest species. The Bray-Curtis coefficient has been recommended and used in ecological gradient studies [42, 43, 47, 48], and, in Amazonia, with birds [4, 5, 49]. Finally, the resulting NMDS axes were used as dependent variables in models of multiple linear statistical models.

A comparison between both bird richness and species composition in the community was made using analysis of variance (ANOVA) and multiple linear statistical models, respectively, to evaluate whether and/or how the bird community changes between high-water and low-water periods. A posteriori Pillai-Trace tests were used to verify whether there were significant differences in bird species composition among the 45 lakes and also when 16 of the 45 lakes during the low-water period (November-December 2007, November-December 2009) were compared with the same 16 lakes for the high-water period (July 2007, July 2008, see bird survey data). The Pillai-Trace statistic has been shown to be less sensitive to deviations from assumptions than other multivariate statistics [50, 51]. To make the water periods comparable, only the morning surveys (6.00 to 10.00 h) were considered. The abundances of most bird species were too low during the peak of flooding; therefore, only a qualitative matrix (presence/absence data) including the two periods of high water and one of low water was used for NMDS. The resulting NMDS axes were also used as dependent variables in models of multiple linear statistical models. The General Linear Model (GLM) procedure available in the program SYSTAT [52] was used to run the analyses.

Multiple linear models followed by Pillai-Trace tests were used to evaluate the effects of lake shape, isolation, and physicochemical characteristics (depth, water transparency, productivity) on the bird qualitative and quantitative community composition, now using data from all 45 lakes (November-December 2007). To verify potential problems of residual analysis in multiple regressions, a graphic method called partial residual plot, available in R (Core Development Team 2008) was used. To verify possible linear relationships among predicting variables, I also used R to estimate the variance inflation factor, which calculates the level of multicolinearity.

3. Results

3.1. Avian Surveys

In this study, the Anavilhanas Archipelago landscape was surveyed along a large gradient of 45 lakes, varying in size (area) from 0.29 km2 (Lake 26) to over 60 km2 (Lake 31), to obtain a nearly complete survey of the aquatic bird community. Lake perimeter ranged from 4 km (Lake 28) to 43.7 km (Lake 31). The lakes were distributed within an archipelago extending 140 km along the Negro River (Figure 1). Since the survey was conducted mostly from a distance of approximately 20 m from lake margins (but most of the time crossing each lake when visiting the central areas to get the geographical coordinates for its position), a minimum of 1,300 km in the 45 lakes was surveyed by speedboat.

A total of 3,626 individuals of 48 species associated with 45 lake environments were recorded at Anavilhanas (Table 2). None are endemic to the Archipelago, and none are listed as threatened. The Fasciated Tiger-Heron (Tigrisoma fasciatum), a bird extinct in southern Brazil, is a new record for the area. During the surveys, three individuals were registered (one in a lake). Most species recorded are permanent residents. One is considered a northern neotropical migrant (the Osprey Pandion haliaetus). The families Ardeidae and Alcedinidae had the highest species richness, with 8 (16.6%) and 5 (10.4%) species, respectively.

The rarefaction curves of aquatic bird species increased until they reached a plateau after 20 lakes (about half of the total) were surveyed, after which they continued to increase, but less than before (Figure 3).

Nine of the ten most abundant aquatic and piscivorous species are among the ten most abundant of the 48 species associated with aquatic environments (Table 2). These species (Neotropic Cormorant Phalacrocorax brasilianus; Black-bellied Whistling-Duck Dendrocygna autumnalis; Greater Ani Crotophaga major; Muscovy Duck Cairina moschata; Cocoi Heron Ardea cocoi; Great Egret Ardea alba; Anhinga; Large-billed Tern Phaetusa simplex; Ringed Kingfisher Megaceryle torquata) summed up 2,567 individuals, representing over four fifths (82.7%) of the 3,105 aquatic birds recorded for the 31 species used in the analyses.

In the 45 lakes, considering the 31 species used in the analyses, the aquatic bird species richness ranged from four (Lakes 19-B and 38) to 19 (Lakes 15 and 38), and the total community abundance ranged from four individuals (Lake 4) to 603 (Lake 12).

3.2. Lake Characteristics

Lake shape (SI) ranged from 1.183 (close to round) to 2.749 (elliptical). Lake distances from river margins ranged from 0.306 km (Lake 40) to 10.4 km (Lakes 5 and 6), and distances to the nearest lake ranged from 0.688 km (Lakes 25 and 26) to 6.3 km (Lakes 42 and 43). Lake depth ranged from 0.88 m (Lake 12) to 7.0 m (Lake 42), and water transparency ranged from 0.47 m (Lake 25) to 1.1 m (Lakes 42 and 46). Lake productivity variation considering nitrogen concentration was higher (from 0.18 in Lakes 21 and 22, to 0.81 mg/L in Lake 16) than for phosphorus (from 0.020 in Lakes 6 and 18 to 0.060 mg/L in Lake 25). Habitat richness (marginal vegetation) ranged from 3 (Lake 12) to 11 habitats (Lake 31). See more detailed information in Table 3.

The correlation between lake habitat richness and distance to nearest lake was significant ( , , ), and so was that between lake shape and phosphorus concentration ( , , ). Therefore, these two variables (habitat richness and lake shape) were not included together with the others in the same statistical models, in which the response variables of bird richness, abundance, and community composition were evaluated. The other independent variables were not significantly correlated among them.

The relationship between distances in the quantitative aquatic bird community composition and distances between lakes was not significant (Mantel test, ; ), indicating absence of spatial autocorrelation, or that bird species composition of the lakes was not grouped by spatial proximity. The relationship between distances in the qualitative aquatic bird community composition and distances between lakes was very weak, although significant (Mantel test, ; ), indicating some degree of spatial autocorrelation.

Because the correlations for quantitative aquatic bird community composition were not significant, and the relationship between the qualitative aquatic bird community composition and distance between lakes were so weak, multivariate analyses were run to compare the quantitative and qualitative bird community composition among lakes, and to evaluate the effects of lake characteristics on them.

There was a positive and highly significant relationship between bird species abundance (mean bird density) and bird species distribution, or the number of lakes in which they occurred (simple linear regression, , , ; Figure 4). The number of lakes in which birds occurred explained 91% of the total variation, indicating that the most abundant species were the ones with the most extensive distribution in the archipelago.

3.3. Effects of Seasonality on Bird Species Composition

The aquatic bird species richness changed significantly among the high-water periods and the low-water periods (ANOVA, ; ), with more species occurring at low water. The species composition of the aquatic bird community during high water (July 2007 and July 2008) was more similar between them than to that during low water (November-December 2007 and November-December 2009, Figure 5) and changed significantly between high- and low- water periods (Pillai-Trace = 0.744, , , Figure 5).

3.4. Bird Species Composition along Gradients Produced by Lake Characteristics

In Figures 6(a) and 6(b), the aquatic bird species are ordered by relative density (bars) along the ecological gradients of lake depth and habitat richness, which influenced significantly bird community composition (see below). Because it depicts a similar pattern to those of Figures 6(a) and 6(b), a figure with the ordination of bird species in relation to that of lake distance to river margins and distance to the nearest neighboring lakes is not shown.

Bird community composition was significantly influenced by lake depth, water transparency, lake isolation, and habitat richness (Table 1). However, in all cases, few species were strongly associated with the limits of the gradients, and most species occurred across the gradients of these variables (Figures 6(a), 6(b)). Differences in species composition (qualitative data) were significantly correlated with lake water transparency, lake isolation and habitat richness (Table 1, Figures 6(a), 6(b)). However, when abundance of individuals was added to the community data (quantitative data), lake depth was also an important determinant and all relations were significant except for water transparency (Table 1).

3.5. Effects of Lake Characteristics on the Aquatic Avifauna

The aquatic bird species richness was not significantly related to lake shape (simple linear regression, ; ; ); however, it was significantly higher in lakes bearing higher habitat richness ( , ; ), and significantly lower in more isolated lakes, or those located farther away in relation to neighbouring lakes ( , ; ). It was not significantly related to any other lake physicochemical properties.

Bird abundance was significantly lower in deeper lakes ( , ; ). However, it was not significantly related to lake shape ( ; ; ), habitat richness of lake margins ( , ; ), or any other lake physicochemical properties or isolation.

For the quantitative data (species abundance), the aquatic bird community composition (NMDS axes) differed significantly with lake isolation (distance to river margins and to the neighbouring lake) and lake depth. Differences in the aquatic bird community composition were not significantly related to water transparency or lake productivity ( and , Table 1).

For the qualitative data (species presence/absence), bird community composition also differed significantly with lake isolation (distance to river margins and to neighbouring lake) and water transparency. However, differences in bird community composition were not significantly related to lake depth and lake productivity ( and , Table 1).

The differences in both quantitative and qualitative aquatic community composition were consistent, and probably, as a result of the degree of lake isolation, decreased significantly with the increase in distance to river margins (Table 1, Figure 7), and increased with the increase in distance to the nearest neighbouring lake (Figure 7). The differences in quantitative community composition increased with lake water depth and decreased with water transparency (Table 1, Figure 7), indicating that aquatic birds use shallow lakes with somewhat turbid waters more often.

Because lake shape and habitat richness were significantly correlated to other variables (see data on correlation between variables), they were included in separate models. For the quantitative data (species abundances), the differences in aquatic bird community composition decreased significantly with the increase in habitat richness. However, differences in the aquatic bird community composition were not significantly related to lake shape (Table 1, Figure 7). For the qualitative data (species presence/absence), the differences in the aquatic bird community composition also decreased significantly as habitat richness increased. However, those differences were not significantly related to lake shape as well (Table 1, Figure 7).

4. Discussion

Several factors have been investigated as determinants in the composition of aquatic bird communities in Brazil. Guadagnin et al. [6] showed the importance of α- and β diversity influencing aquatic bird communities in fragmented wetlands in southern Rio Grande do Sul. They found that some aquatic birds favour small, isolated lakes and suggested among other parameters that waterbody physicochemical conditions, food resources, and habitat, may influence ecological gradients and therefore deserve further investigation.

This study clearly provides further evidence that the natural limnological heterogeneity and flood dynamics of Amazonian black water lakes can affect the local composition of animal and plant communities [5, 39, 53, 54]. In addition to some variables significantly associated with the number of bird species (richness) and individuals (abundance), there were strong and significant relationships between lake depth, isolation, water transparency and bird community composition. Furthermore, I also found that variation in bird community composition can, in part, be attributed to a habitat richness gradient, and the phases of the annual hydrological cycle in Anavilhanas.

Most of the findings were in community composition and fewer in variation in species richness and bird abundance. This can be attributed to species turnover within the archipelago of Anavilhanas, which is not related to regional variation in the aquatic avifauna (unpublished data), but to local differences in the avian communities. The surveys indicate that most bird species occur throughout the lakes of the Anavilhanas and are tracking differences in the archipelago’s landscape and the characteristics of the lakes.

In the study area, there is a strong and significant relationship between the locally most abundant species and the range of their distributions across the Archipelago (Figure 4), which can be related to population dispersion as predicted by the neutral theory of biodiversity [9, 55, 56]. Nevertheless, in the Anavilhanas, the weak correlations between the spatial variation in the bird community composition and geographical distances (distances between lakes, see Section 3 on spatial autocorrelation) can be taken as evidence against the distance decay predicted by Hubbell’s theory [9]. Therefore, long-term studies are needed to specifically test this theory for the lake bird community in the Anavilhanas. In European wetlands, abundance has been considered a good predictor of the ability of species to colonize lakes, with the less abundant being less able to occupy distinct areas and occurring in few large ponds [57], a pattern that seems valid for the Anavilhanas. For example, the Agami Heron, the Sunbittern, and the Horned Screamer, among others, were all very scarce in the surveys (see Table 2) but occurred in a small number of medium to large lakes.

Some lakes were not visited because access was blocked by shallow waters. However, 60–70% of all lakes in the archipelago were surveyed. Therefore, the sampling was extensive and representative of this enormous black water system (Figures 1 and 2). In addition, over 95% of the bird species associated with aquatic environments previously recorded in the archipelago [4] were recorded in this study.

Of the ten most abundant species in lakes, six are piscivorous (Neotropic Cormorant, Cocoi Heron, Great Egret, Anhinga, Large-billed Tern and Ringed Kingfisher), and one is a migrant, piscivorous species (Osprey) whose juveniles can stay in the area for a year or more before they migrate to the Northern Hemisphere (unpublished data). Six also showed ample distribution in the archipelago, with records in more than 80% of the lakes (see Table 2). The Neotropic Cormorant was the most abundant species (nearly 500 individuals). The Black-collared Hawk and the Horned Screamer (abundant in white water Lakes), the Gray-necked Wood-Rail, the Agamia Heron Agamia agami; the Fasciated Tiger-Heron Tigrisoma fasciatum; the Sunbittern Eurypyga helias; were among the ten rarest species (Table 2).

With similar species composition, but lower bird abundance than those from Amazonian white water lakes [5], the communities of birds associated with black water environments at Anavilhanas comprise a few abundant and a higher number of rare species. This is a typical pattern of large animal communities, not only in the tropics, but also in other latitudes [1, 58, 59].

In this study, the aquatic bird species accumulation curve did not reach an asymptote, and more species were being added even after 40 lakes (Figure 3). This indicates that other species that are present in the area could potentially have been added. Indeed, some species unrecorded in the survey have been seen in the area, and also in other conservation units in central and western Brazilian Amazonia [4, 5, 9, 19]. This notwithstanding, most species (30, or 96.8% of the aquatic; 17, or 94.4% of the piscivorous) were found before less than half of the lakes (20, or 44.4%) were surveyed (Figure 3), and several lakes had to be visited before another new species was added. In white water lakes in western Brazilian Amazonia, a similar number of bird species were recorded after surveying only five lakes [5]. This difference may be due to the higher lake productivity of the várzea wetlands.

In the Anavilhanas, there were significant changes in species composition between periods of high and low water, with less variation in low-water than in high-water season (Figure 5). Although species richness is similar between these two periods, only slightly lower during high water, abundance is notably higher during low water. This may result from some species having nonresident individuals which are, therefore, transient in the area, for example, Anhinga, terns, egrets, herons, and ducks (Black-bellied Whistling-Duck). During low waters food resources are more abundant and concentrated because the lakes are drying out, and probably more accessible. In his intensive and detailed study at Marchantaria Island, Central Amazonia, Petermann [11] also found seasonal changes in the composition of the wetland bird fauna, with lower numbers of species occurring during the high water period (flood pulse).

The results of multivariate analyses suggest that, in general, the aquatic bird community composition varies considerably with a gradient in ecological factors across the Archipelago, such as lake depth, lake isolation, and habitat richness (Table 1). Morphometric features of lakes are important determinants of bird and fish community structure in north-central Alberta, Canada [60]; probably, environmental conditions at the margins of round and shallow lakes are more homogeneous, and littoral production higher, than in deeper lakes that have a greater amount of pelagic habitat relative to the littoral zone [55, 61]. Considering the entire bird assemblage, there were no effects of lake shape on the community, contrary to another recent study in western Brazilian Amazonia (Amanã and Mamirauá Reserves), where lake morphometry influenced variation in the avian community [5]. This difference probably is due to the lower heterogeneity in the shapes of the Anavilhanas lakes, which tend to be more elliptical, whereas the ones from Amanã and Mamirauá vary from circular to elongated, and dendritic [5].

There was a trend of higher bird richness and abundance and larger variation in bird community composition in shallower lakes (see Section 3, Figures 6(a) and 7), likely linked to the fact that many piscivorous and omnivorous birds, such as egrets, herons, terns, and ducks (Figure 8), capture prey (aquatic insects, crabs, snails, and small fish) more easily in shallow than in deeper lakes ([55], see also [62]). Higher light penetration in shallow lakes may increase prey detectability, principally for visually oriented predators. In this sense, the lack of relationship observed between bird richness and abundance and water transparency was an unexpected result, but one that can be derived from the fact that some waders and piscivorous birds (egrets, herons, and others) forage more often in relatively turbid waters or shady spots, possibly to avoid being visually detected by their prey (personal observations) and/or because prey biomass tends to be higher in these waters. However, variation in aquatic bird qualitative community composition significantly decreased with increasing water transparency (Table 1, Figure 7) indicating lower species turnover in clearer lakes.

Opposite to studies in North America and Europe [62, 63], in the Anavilhanas, lake productivity did not influence variation in bird community composition, although it may affect bird species individually [64].

Lake isolation was another important determinant of the ecological gradient in the aquatic community of the Anavilhanas. Aquatic bird richness was higher in more isolated lakes. Isolation also influenced bird community composition (Table 1, Figure 7), which showed the highest variation in species in lakes closer to river margins, as expected due to more potential bird colonization from the “mainland” (lakes along river margins). In the Anavilhanas, differences in bird occurrence may stem from differential vagility and ability to colonize lakes farther away from river margins. For example, the Sunbittern, the Sungrebe, the Black-collared Hawk, and the American Pygmy Kingfisher, have small territories; the Striated Heron and the Rufescent Tiger-Heron are quite sedentary.

The distribution of species’ relative densities along, (Figures 6(a), 6(b)) and species variation in relation to, ecological gradients (Figure 7) suggest that some species prefer one side of the gradients (e.g., deeper water, higher habitat richness, and greater isolation), which is one of the reasons why bird richness and abundance were weakly affected by them.

In this study, habitat richness of marginal vegetation emerged as one of the most important lake properties affecting aquatic bird community composition; similar results were found in temperate regions [62, 63, 65]. At the Anavilhanas, aquatic birds forage more frequently in lakes than along river margins (personal observations). This may stem from differences in water flow (slower waters favouring feeding success) and to a greater heterogeneity of lake margin vegetation [66]. Habitat diversity may be especially important in the prebreeding season, since a more complex structure may provide more isolated spots for social reproductive behaviour, such as display and courtship, as Turnbull and Baldassarre [65] found for Mallards and American Wigeons in the Northern Hemisphere.

Higher bird species richness and abundances are consistently found in lakes with higher habitat richness along the margins. Anavilhanas lakes differ in their potential to provide adequate habitat to aquatic birds because species have different biological needs which may influence the way each uses the lakes. For some species (e.g., Horned Screamer, Hoatzin, Rufescent Tiger Heron, Gray-necked Wood-Rail, Green, Black-bellied Whistling-Duck, Muscovy Duck, Greater Ani), habitat preference is likely related to local vegetation type and structure within the lakes; the occurrence of other species (e.g., Anhinga, Neotropic Cormorant, Great Egret, Sunbittern) is related to landscape features at larger spatial scales, while others (Large-billed Tern, Yellow-billed Tern, Ringed Kingfisher, Osprey and swallows) do not depend so much on lake habitat attributes because they use mainly open waters and beaches [66]. These differences in habitat use may contribute to maintain alpha diversity and produce spatial gradients [6].

This study’s findings strongly suggest that the environmental heterogeneities created by differences in lake isolation, physical properties, and habitat richness are important determinants of ecological gradients in aquatic avian richness, abundance, and community composition in the black water systems of Central Amazonia, particularly in the Anavilhanas archipelago.

An enormous archipelago comprised of hundreds of beautiful lakes stretching throughout its area in the middle of a National Park is a unique situation in the entire Amazon region. Since the distribution and abundance of fish is influenced by waterbirds [8], a long-term basis monitoring program of the bird community dynamics will be important to improve the lakes’ conservation value, and to analyze human impacts on them. Future analysis of multiyear data will allow to make temporal comparisons among lakes and to assess if trends reflect local and/or regional patterns of community change in the Amazonian black water system of Anavilhanas.

Acknowledgments

The author thanks the ICMBio Staff, Adenilson Azevedo dos Santos, Sr. Antônio, Jefferson Gomes dos Santos, Antônio Craveiro Martins, Francisco de Assis dos Santos, Raimundo Oliveira da Silva, Janis Monteiro, Marli Gomes da Silva, Alegário Monteiro dos Santos, for their assistance with local logistics. The author thanks the ICMBio researchers at Novo Airão, Bruno Marchena Tardio, Tatiana Alves, Giovanna Palazzi, and Igor Matos Soares for logistical support, hospitality and permission to work in the area. Thanks to Melissa Rosas, Giselle Casetta (in 2007), Adriana Terra, Jessica Cancelli, Pedro Santos (in 2009), and José Antônio Melgueiro dos Santos (Queca) for their diligence and competence during field work. Melissa Rosas helped to draw Figure 1 of the study area. Ítalo Mourthé drew Figure 3, and Helena Aguiar and Victor Landeiro made Figures 6(a) and 6(b) using “R” statistical software. Bruce Nelson (INPA) introduced me to Global Mapper Software. Pedro Santos, José Eugênio Figueiras, Tânia Sanaiotti, and Summer Wilson gave suggestions and revised the text. Andreia Guadalupe (INPA Laboratory of Geosciences) helped with water productivity analysis. Maria Tereza Piedade Fernandez Piedade (Maitê) sent me some articles on Amazonian aquatic plants. Finnaly, I want to thank Luiz Cláudio Marigo who gently make available the excelent quality photos on waterbirds and Anavilhanas Lakes. Financial support was granted by Instituto Nacional de Pesquisas da Amazônia (INPA) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq-Processo: 472276/2008-5), Brazil.