Layout 1


INTRODUCTION

Cladocera are one of the most important components
of zooplankton in lacustrine ecosystems. This group is also
well-known to promptly respond to environmental changes
(Jeppesen et al., 2001). Previous studies have highlighted
that temperature and altitude are among the main factors
determining Cladocera communities (Campanelli, Mortari
and Henry, 2016; Green, 1995; Padhye et al., 2016; Sweet-
man, 2010; Zawiska et al., 2015). In addition, these studies
have shown comparable sensitivities to changes in lake
water trophic state, water depth, conductivity and pH (Ko-
rhola and Rautio, 2001). These responses to environmental
variables makes Cladocera one of the most powerful pale-
olimnological indicators and for this reason they have been
widely used for the reconstruction of past lacustrine envi-
ronments (Bjerring, 2007; Kienast et al., 2011; Korosi,
2012; Luoto et al., 2009; Mirosław-Grabowska and Zaw-
isza, 2014; Nováková et al., 2013 Sarmaja-Korjonen and
Hyvärinen, 2008; Schmidt, 2000; Szeroczyńska, 1991; Sze-
roczyńska et al., 2007; Szeoczyńska and Zawisza, 2011
a,b). Subfossil remains of Cladocera are commonly used in

paleolimnological studies in Eurasia and America (Bjer-
ring, 2007; Mirosław-Grabowska and Zawisza, 2014; Pa-
terson, 1994; Szeroczyńska et al., 2007). However, only a
few paleolimnological studies based on water fleas were
conducted in Central America (Cuna et al., 2014; Sze-
roczyńska et al., 2015; Zawisza et al., 2012, 2014, 2016).
So far, the majority of studies focused on the living Clado-
cera and their taxonomy (Elías-Gutiérrez et al., 1999, 2006,
2008; Sinev, 2015; Sinev and Silva-Briano, 2012; Sinev
and Zawisza, 2013). 

The knowledge of Cladocera species composition, dis-
tribution and ecology in freshwater lakes of Central Amer-
ica is crucial for the development of reliable bioindicators,
and for paleolimnological investigations. The main objec-
tive of the present study was to analyze and identify sub-
fossil Cladocera species in northern Central America, in
particular in Guatemala, El Salvador and Honduras, and
to explore their relationship to lake environmental char-
acteristics such as altitude, lake area, depth, pH, trans-
parency, conductivity and dissolved oxygen. The obtained
results will provide the basis for the reconstruction of re-
gional and global climatic and envoronemntal changes at
long-term scale. 

Advances in Oceanography and Limnology, 2016; 7(2): 151-162 ARTICLE
DOI: 10.4081/aiol.2016.6266
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0).

Ecology of Cladocera species from Central America based on subfossil assemblages

Marta Wojewódka,1* Edyta Zawisza,1 Sergio Cohuo,2 Laura Macario-González,2 Antje Schwalb,2 Izabela Zawiska,3
Liseth Pérez4

1Institute of Geological Sciences, Research Centre in Warsaw, Polish Academy of Sciences, Twarda 51/55, 00818 Warsaw, Poland;
2Institut für Geosysteme und Bioindikation, Technische Universität Braunschweig, Langer Kamp 19c, 38106 Braunschweig, Germany;
3Institute of Geography and Spatial Organization, Polish Academy of Sciences, Twarda 51/55, 00818 Warsaw, Poland; 4Instituto de
Geología, Universidad Nacional Autónoma de México (UNAM), Ciudad Universitaria, 04510, Ciudad de México, México

*Corresponding author: m.wojed@twarda.pan.pl

ABSTRACT
Cladocera species composition was analyzed in surface sediments of 29 lakes in Central America (Guatemala, El Salvador and Hon-

duras). The material studied was collected with an Ekman grab in autumn 2013 from lakes located in lowland, highland and mountain re-
gions. The study revealed high variability in qualitative and quantitative composition of subfossil Cladocera. A total of 31 Cladocera species
(5 planktonic and 26 littoral) were identified, as well as 4 morphotypes that could not be identified (NRR 1-4). Planktonic Bosminidae and
Daphniidae were the most abundant families. Daphniidae were restricted to water bodies in mountain regions, whereas Bosminidae were
widely distributed in lakes with different abiotic conditions. Moreover, Bosminidae species also occurred in highly mineralized waters (>
900 µS cm–1). The great majority of the identified Cladocera species belonged to the littoral family Chydoridae. Chydorus cf. sphaericus
was the most common species (found in 20 lakes), which probably reflects its tolerance to a wide spectrum of habitat conditions. Cluster
analysis discriminated 6 groups of Cladocera species with a high correlation level within groups (≥0.8), which showed different types of
correlation with lake characteristics and environmental variables. Canonical correspondence analysis (CCA) showed that altitude and sec-
ondly water electrical conductivity were the most important drivers of Cladocera species composition in the region studied. Furthermore,
CCA analysis indicated lowland lakes with low water transparency were also characterized by peculiar species assemblages.

Key words: Subfossil Cladocera; Central America; waterfleas ecology; Guatemala; El Salvador; Honduras.

Received: September 2016. Accepted: November 2016.

No
n c

om
me

rci
al 

us
e o

nly



152 M. Wojewódka et al.

METHODS

Study sites

The study included 29 waterbodies located between 13°
and 17° N geographical latitude and from 87° to 91° W lon-
gitude in Guatemala, El Salvador and Honduras (Fig. 1).
The study area is located in the Central American Bioregion
of the American Neotropical Ecozone, which is character-
ized by two distinct seasons: a dry winter and a moist sum-
mer season (Taylor and Alfaro, 2005). The subtropical ridge
(belt of high atmospheric pressure) from the northern At-
lantic is controlling the climate of the region, whereas the
Intertropical Convergence Zone (ITCZ) and polar fronts of
mid-latitude origin have a secondary impact (Taylor and
Alfaro, 2005). One of the main factors affecting the micro-
climate of the area is topography (Karmalkar et al., 2011;
Taylor and Alfaro, 2005). The considerable differences in

mean annual temperature and precipitation which charac-
terize the entire region mainly depend on altitude. The
Guatemalan lowlands are characterized by a mean annual
temperature of 25-28°C, whereas temperatures in highland
and mountain regions are around 12-15°C (Atlas Climato-
logico for 1928-2003, the National Institute for Seismology,
Volcanology, Meteorology and Hydrology of Guatemala;
www.insivumeh.gob.gt). Annual precipitation in the area
varies considerably during the year and ranges from <1000
mm in the plains to >2500 mm in the mountains (Taylor
and Alfaro, 2005). 

The lakes studied are located at altitudes ranging from
sea level (El Muchacho at 3 m asl) to high altitude lakes
(Magdalena and Chicabal at 2863 and 2726 m asl, respec-
tively, Tab. 1). Lake area and depth also were quite vari-
able. The largest lake were Atitlán and Yojoa (125 and
~79 km2, respectively, Tab. 1), while several lakes have
an area <0.10 km2, such as Madre Vieja (0.10 km2), Verde

Fig. 1. Map of the study area showing the location of lakes sampled in Central America, i.e. Guatemala, Honduras and El Salvador.
Numbers indicate the names of the lakes: 1, Atitlán; 2, Calderas; 3, Chicabal; 4, Comandador; 5, El Muchacho; 6, El Pino; 7, El Rosario;
8, Grande; 9, Ipala; 10, Lachuá; 11, Las Pozas; 12, Magdalena; 13, Sacnab; 14, Salpetén; 15, Quexil; 16, Chiligatoro; 17, Jucutuma; 18,
Madre Vieja; 19, Ticamaya; 20, Yojoa; 21, Apastepeque; 22, Aramuaca; 23, Chanmico; 24, El Espino; 25, Jocotal; 26, Los Negritos;
27, Metapan; 28, Olomega; 29, Verde.

No
n c

om
me

rci
al 

us
e o

nly



Subfossil Cladocera of Central America 153

(0.10 km2), Chiligatoro (0.04 km2), El Rosario (0.02 km2),
Magdalena (0.01 km2) (Tab. 1). The study covered both
very shallow (≤3 m depth, such as Comandador, El
Muchacho, El Rosario, Grande, Magdalena, Jucutuma,
Ticamaya, Jocotal, Los Negritos, Olomega), and very
deep lakes (>300 m), such as Lachuá (a 378 m deep
karstic lake), Atitlán (a 340 m deep crater lake) and Chi-
cabal (330 m, Tab. 1).

Sampling and analyses

Sediments were collected in autumn 2013 using an
Ekman sediment sampler. The geographical location of the
lakes was determined by a handheld navigator (GPSmap
60c). Water transparency was determined using a Secchi
disk, and the maximum lake depth was identified using a

portable depth sounder (Echosounder Eagle Mach 1).
Physical and chemical parameters of the surface waters
(i.e., pH, conductivity, dissolved oxygen) were measured
with a WTW multi set 350i multiparametric probe. The
area of each lake (Tab. 1) was calculated by the measure
tool of Google Earth. Lake sediment samples were col-
lected both from the littoral as well as from the open-water
(pelagic) zone. In eight shallow lakes, where there was no
pelagic zone (i.e., Comandador, El Muchacho, El Rosario,
Grande, Jucutuma, Ticamaya, Jocotal, Los Negritos), sam-
ples were collected only from the littoral zone (Tab. 1).
Lakes El Pino, Ipala, Quexil, Madre Vieja, Chanmico, El
Espino, Metapan and Verde, were sampled only in the pro-
fundal zone due to the shape of the lake basin. Crater lake
Aramuaca (107 m depth) was sampled only in the littoral
zone (Tab. 1). 

Tab. 1. Location, morphometry and selected limnological variables measured at water surface at the waterbodies sampled in Central
America. 

ID  Name                Country Geographic       Altitude     Area Samples           Zmax   Secchi         pH     Conductivity    O2      HCl
                                             coordinates
                                                               N              W        (m asl)     (km2)          L              P          (m)      (m)              -          (µS cm–1)   (mgL–1)     

1     Atitlán               Guatemala         14.6837     91.2239      1556      125.00         +              +          340      4.2            8.9             442            5.7        +
2     Calderas                                      14.4117     90.5913      1790        0.35           +              +           26       2.9            9.2             100            4.9        -
3     Chicabal                                      14.7875     91.6561      2726        0.21           +              +          330      1.6            9.0              12             5.4        -
4     Comandador                               13.9600     90.2544        20          0.65           +                         1.7       0.5            7.4             251            5.6        -
5     El Muchacho                              13.8892     90.1918         3           0.36           +                           2        0.4            9.1             439            5.9        -
6     El Pino                                        14.3447     90.3941      1038        0.64                         +           18       2.8            8.3             100            2.0        -
7     El Rosario                                   16.5255     90.1601       126         0.02           +                           3        0.2            7.2             987            1.7        +
8     Grande                                        13.8903     90.1703         5           0.95           +                           2        0.2            7.4             110            3.7        -
9     Ipala                                            14.5571     89.6394      1495        0.56                         +           25       3.6            8.0             100            2.8        -
10   Lachuá                                        15.9184     90.6732       170         4.00           +                         378      4.4            7.9             906            4.8        +
11   Las Pozas                                    16.3432     90.1660       152         2.16           +              +           35       3.1            8.3             277            5.4        +
12   Magdalena                                  15.5426     91.3956      2863        0.01           +              +            3        2.8            8.8             331            6.2        +
13   Sacnab                                        17.0583     89.3725       170         4.28           +              +            9        1.8            9.0             412            6.1        +
14   Salpetén                                      16.9815     89.6755       105         2.77           +              +           32       1.7            7.0            4520           5.8        +
15   Quexil                                         16.9231     89.8099       120         2.20                         +           32       2.7            8.5             204            4.5        -
16   Chiligatoro        Honduras           14.3756     88.1830      1925        0.04           +              +           5.5       0.9            7.4             100            1.9        -
17   Jucutuma                                    15.5123     87.9028        27          4.34           +                           2        2.0            7.3             100            1.0        -
18   Madre Vieja                                14.3569     88.1376      1866        0.10                         +           3.4       0.9            8.5             100            2.7        -
19   Ticamaya                                    15.5506     87.8897        17          2.91           +                           2        0.8            7.2             100            1.6        -
20   Yojoa                                          14.8606     87.9847       639        79.70          +              +           22       3.3            8.3             100            2.7        -
21   Apastepeque     Salvador            13.6925     88.7448       509         0.38           +              +           47       6.1            8.6             100            2.8        -
22   Aramuaca                                   13.4294     88.1065        76          0.40           +                         107      6.7            8.4             100            4.1        -
23   Chanmico                                   13.7786     89.3541       477         0.78                         +           51       0.9            9.2             100            3.4        +
24   El Espino                                    13.9530     89.8652       689         0.99                         +           5.5       0.4            8.5              85             6.6        -
25   Jocotal                                        13.3371     88.2519        26          8.70           +                           3        1.5            8.0             595            3.0        -
26   Los Negritos                               13.2831     87.9370       102         0.69           +                           2        0.3            9.2              40             5.2        -
27   Metapan                                      14.3094     89.4655       450        16.00                        +            6        0.2            8.4             255            3.1        -
28   Olomega                                     13.3072     88.0551        66         25.20          +              +           2.5       0.9            7.7             105            2.5        -

29   Verde                                          13.8915     89.7872      1609        0.10                         +           12       2.7            7.5              83             4.5        -
L, littoral zone; P, pelagic zone; Zmax, maximum lake depth; HCl, reaction with acid during preparation. The World Geodetic System of 1984 (WGS84)
datum was used for the Cladocera samples.

No
n c

om
me

rci
al 

us
e o

nly



154 M. Wojewódka et al.

Sediment samples were placed in plastic bottles
(100 mL) immediately after collection, and stored under
refrigeration. Samples were then transported to the Pale-
olimnology Laboratory of Universidad National Au-
tonomous de Mexico (UNAM), where subsamples for
subfossil Cladocera analysis were obtained and preserved
in cooling condition. 

Subfossil Cladocera were analyzed at the Bioindicator
Laboratory at the Warsaw Research Center of the Institute
of Geological Sciences, the Polish Academy of Sciences.
Each sample consisting of 1 cm3 of wet sediment was pre-
pared according to standard methods (Frey, 1986). Sedi-
ments were treated with 10% HCl to eliminate carbonate
and then heated in 10% KOH for 30 min. The residue was
centrifuged and sieved through a 38 µm sieve, transferred
into a test tube and filled up to 5 cm3 with distilled water.
Samples were stained with safranin. Three to four slides
obtained from of 0.1 cm3 liquid suspension each were ex-
amined at the microscope for each sample. The identifi-
cation of Cladocera taxa followed Elías-Gutiérrez et al.
(2006, 2008), Hudec (2000), Kotov (2009), Kotov et al.
(2003), Rajapaksa and Fernando (1987), Rey and Vasquez
(1986), Sinev (1998, 2001, 2015a, 2015b), Sinev and
Dumont (2016), Sinev and Zawisza (2013), Van Damme
et al. (2011). The Cladocera relative abundance diagram
was prepared using the C2 program (Juggins 2005, 2007).

In order to determine the relation between species dis-
tribution and lake characteristics (area, altitude, maximum
depth) and water variables (pH, Secchi transparency, dis-
solved oxygen), cluster analysis and canonical correspon-
dence analysis (CCA; ter Braak, 1986) were performed.
Cluster analysis is based on the similarity measure per-
formed on qualitative data (Hammer et al., 2011), while
CCA is a combination of mutual averaging algorithm with
a multiple regression (ter Braak, 1986) which performed
on combined quatitative Cladocera and environmental
data. Statistical analyses were conducted using the PAST
software (Hammer et al., 2001).

RESULTS

Environmental variables

In general, the surveyed lakes were characterized by
relatively low water transparency, with exception of lakes
Aramuaca and Apastepeque, which showed Secchi depth
>6 m (Tab. 1). Secchi depth values ranging between 2.7-
4.4 m were measured in 10 lakes, i.e. Lachuá, Atitlán,
Ipala, Yojoa , Las Pozas , Calderas, El Pino, Magdalena,
Quexil and Verde (Tab. 1). The waters of Lake Jucutuma
(2.0 m depth) were characterized by visibility reaching
the bottom, while four lakes were characterized by clarity
of 1.5-1.8 m (Sacnab, Salpetén, Chicabal, Jocotal). The
other lakes (Comandador, El Muchacho, El Rosario,

Grande, Chiligatoro, Madre Vieja, Ticamaya, Chanmico,
El Espino, Los Negritos, Metapan, Olomega) were char-
acterized by Secchi disk visibility <1.0 m (Tab. 1). 

The lakes investigated did not show considerable dif-
ferences in pH values, which ranged from neutral (7.0 in
Lake Salpetén) to alkaline (9.2 in Lake Chanmico and
Lake Los Negritos, Tab. 1). Circumneutral water pH was
recorded for lakes El Rosario, Ticamaya, Jucutuma,
Comandador, Chiligatoro, Grande, Olomega, Verde and
Lachuá. Waters of five lakes were strongly alkaline (pH
≥9.0), i.e., Sacnab, El Muchacho, Calderas, Los Negritos,
and Chanmico. The other 14 lakes were characterized by
slightly alkaline waters, with pH values ranging from 8.0
to 9.0 (Tab. 1). Electric conductivity ranged from 12 µS
cm–1, measured at Lake Chicabal, to 4520 µS cm–1 in Lake
Salpetén (Tab. 1). High conductivity was recorded in the
waters of El Rosario and Lachuá (Tab. 1). Values between
~400 and 600 µS cm–1 were determined in lakes Jocotal,
Atitlan, El Muchacho, and Sacnab (Tab. 1). In the remain-
ing lakes this parameter was <350 µS cm–1 (Tab. 1). Ten
of the sampled lakes were characterized by the dissolved
oxygen content of over 5.0 mg L–¹, i.e., El Espino, Mag-
dalena, Sacnab, El Muchacho, Salpetén, Atitlán, Coman-
dador, Las Pozas, Chicabal and Los Negritos (Tab. 1). The
amount of dissolved oxygen in nineteen lakes ranged from
1.0 mg L–¹ (Jucutuma) to 4.9 mg L–¹ (Calderas). The low-
est dissolved oxygen concentration (<2.0 mg L–¹) were
recorded in lakes Chiligatoro, El Rosario, Ticamaya and
Jucutuma (Tab. 1). 

Subfossil Cladocera

A total of 31 Cladocera taxa belonging to three fami-
lies – Daphniidae, Bosminidae and Chydoridae – were
found in the sediments of the lakes studied. Some Clado-
cera remains that could not be identified to species level
belong to four different morphotypes and were referred
to as “Not Recognized Remains” (NRR) 1, 2, 3 and 4, re-
spectively (Fig. 2). Species relative abundances are pre-
sented in Fig. 3. Frequency (n), mean and maximum
relative abundance of each Cladocera species are collated
in Tab. 2.

The species richness per lake ranged from one in Lake
Calderas to 12 species (lakes Atitlán and Verde), and in-
cluded both planktonic and littoral forms. Only in the sed-
iments of the deep Aramuaca crater lake no subfossil
cladocerans remains were found. Planktonic species were
represented by five taxa, i.e., Bosmina longirostris
(O.F.Müller, 1785), Bosmina (E.) longispina (Leydig,
1860), Bosmina (E.) coregoni (Baird, 1857), group of
Daphnia longispina (O.F. Müller, 1785), group of Daph-
nia pulex (Leydig, 1860), which occurred in 17 lakes
(Fig. 3). Cladocera communities from lakes El Rosario
and Calderas were entirely dominated by planktonic
species (Fig. 3), while the share of Bosminidae and/or

No
n c

om
me

rci
al 

us
e o

nly



Subfossil Cladocera of Central America 155

Daphniidae exceeded 90% in lakes Ipala, El Pino and
Apastepeque (Fig. 3). Pelagic species were dominant (60-
88% of total Cladocera abundance) in lakes Salpetén, Yojoa,
El Espino, Atitlán, Sacnab, Chicabal, Lachuá (Fig. 3). 

The remaining lakes (16) were dominated by littoral
species of Chydoridae, which also represented the most
species-rich group. In fact, 26 littoral Cladocera species
and four NRR (Fig. 3 and Tab. 2) were identified in the
sediments studied, with the NRR most likely belonging
to also to littoral Chydoridae. Chydorus cf. sphaericus
(O.F. Müller, 1776) (20 lakes), Alona glabra (Sars, 1901)
(Ovalona glabra, following Sinev, 2015a) (13 lakes) and
Cladocera comparable to (13 lakes) and Alona quadran-
gularis type (O.F. Müller, 1776) (12 lakes) were the most
common taxa (Fig. 3 and Tab. 2). In five lakes (Las Pozas,
Quexil, Verde, Chiligatoro, Magdalena), both littoral and
planktonic species were present. The share of littoral
species ranged from 60% to 80% in lakes Las Pozas,

Quexil, Verde and Chiligatoro. In the sediments of lakes
Las Pozas and Quexil, eight and nine littoral species were
identified, respectively, and two planktonic species were
detected. In Las Pozas, Alona ossiani (Sinev, 1998), was
the most abundant littoral species (approx. 20%), whereas
in Quexil – Chydorus cf. sphaericus (ca. 15%) and Alona
quadrangularis type (ca. 15%, Fig. 3) were dominant.
Lake Verde and Lake Chiligatoro were characterized by
the presence of 12 and 6 species, respectively. Chydorus
cf. sphaericus (ca. 25%) and group of Anthalona verru-
cosa (Sars, 1901) (ca. 15%) were the dominant species of
the littoral zone in Lake Verde, while Alona ossiani (ca.
35%) and Simocephalus sp. (Schoedler, 1858) (ca. 20%)
were the most abundant species in Lake Chiligatoro (Fig. 3
and Tab. 2). In Lake Magdalena, the deep-water species
accounted for a minor portion of cladocerans, and Chy-
dorus cf. sphaericus was the dominant species (ca. 90%).
Planktonic taxa were completely absent in the sediments

Fig. 2. Unidentified subfossil remains (NRR) of Cladocera. A) Probably headshield of Euryalona sp. B) Headshield NRR1. C) Postab-
domen of Leydigiopsis ornata. D) Probably headshield of Leberis sp. E) NRR3, probably shell of Cladocera from group of Coronatella
monacantha. F) shell NRR4. Scale bars: 500 µm.

No
n c

om
me

rci
al 

us
e o

nly



156 M. Wojewódka et al.

of eleven lakes (Fig. 3). Lakes Comandador, Los Negritos,
Ticamaya and Jocotal were dominated by Chydorus cf.
sphaericus, which represented from 36% to 50% of the
Cladocera communities. This taxon was also common in
Madre Vieja (15%), Metapan (23%) and Olomega (16%).
In Lake Madre Vieja, Alona ossiani was the most abun-
dant species (ca. 42%), Camptocercus dadayi (Stingelin,
1913) was the dominant species in Metapan (ca. 31%),
and NRR1 (ca. 32%) (Figs. 2 and 3) were dominant in
Olomega. In Lake Jucutuma, Acroperus sp. (Baird, 1843)
(ca. 22%) and Camptocercus dadayi (ca. 22%) were dom-
inant, whereas Alona quadrangularis type (ca. 43%) and
Alona glabra (ca. 29%) dominated in Lake Grande. Lake
El Muchacho and Lake Chanmico were characterized by
low species diversity (3) and low frequency of individu-
als. Group of Leydigia acanthocercoides (Fischer, 1854)
was dominant in Lake El Muchacho (ca. 52%), whereas
Leberis sp. (Smirnov, 1989) was the most prominent
species in Lake Chanmico (50%, Fig. 3). The largest num-

ber of subfossil Cladocera remains was determined for
Lakes Sacnab (22,000 ind. cm–3) and Ipala (21,000 ind.
cm–3), where open-water species dominated. The maxi-
mum number of water fleas in the waterbodies dominated
by littoral forms was ca. 10,000-11,000 ind. cm–3, as
recorded in Madre Vieja, Quexil and Verde. 

Statistical analysis

The “cophenetic coefficient” in the cluster analysis was
0.8009 and well reflected the similarities within the data.
The dendrogram shown in Fig. 4 (see also Tab. 3) defined
six Cladocera groups. The highest correlation coefficient
between species was determined for Alona rustica (Scott,
1895) (group of Flavalona rustica, following Sinev and
Dumont, 2016) and Simocephalus sp. (0.99, Group 6), and
for Daphnia pulex-group and Bosmina longirostris (0.94,
Group 1, Fig. 4a and Tab. 3). The species Bosmina (E.)
coregoni, Bosmina (E.) longispina, Alona quadrangularis

Fig. 3. Relative abundances (%) diagram for the subfossil Cladocera species identified in surface sediments from freshwater ecosystems
of Central America. Lakes are ordered by decreasing contribution of planktonic forms. NRR1, unidentified remains type 1; NRR2,
unidentified remains type 2, probably headshield of Leberis sp.; NRR3, unidentified remains type 3, probably of Coronatella of monacan-
tha-group); NRR4, unidentified remains type 4.

No
n c

om
me

rci
al 

us
e o

nly



Subfossil Cladocera of Central America 157

type, Alona glabra and Leberis sp. (Group 2) resulted to
be correlated at the similarity level of 0.92. In the three
other groups the correlation coefficient between species
ranges between 0.8 and 0.86 (Fig. 4a). 

The CCA analysis was conducted using the species in-
cluded the six distinguished clusters and an additional set
of other Cladocera species that were not included in any
cluster. Due to their rarity and low contribution to total
Cladocera remains, NRR were not included in the analysis.
The CCA ordination biplot is presented in Fig. 4b. Arrows
show environmental variables, and their length and direc-
tion indicate the strength and direction of their influence on
the species composition (ter Braak, 1986). 

Axis 1 and 2 explained almost 80% of the variability
(axis 1=49.98%, axis 2=29.29%). Altitude, transparency, pH
and the maximum depth were positively correlated with axis
1, while conductivity and lake area were negatively corre-
lated with the first canonical axis (Fig. 4b). The location of
Group 1 in quadrant I, and of Group 3, Group 6 and the
Daphnia longispina-group in quadrant IV underlines the
positive correlation of the abundances of these species
groups with altitude and negative correlation with conduc-
tivity. Moreover, the occurrence of Graptoleberis testudi-
naria (Fischer, 1848) and Group 1 also seem to be strongly
determined by the Secchi disk visibility (Fig. 4b). Species and
groups of Cladocera located in quadrants II and III were cor-
related negatively with altitude (Fig. 4b) and positively with
conductivity and lake area (Fig. 4b and Tab. 3). Group 2 and
Alonella excisa (Fischer, 1854) were placed in the 2nd quar-
ter. Species of Group 4 and 5 and Epheromorphus sp. (Frey,
1982), Leydigia louisi louisi (Jenkin, 1934), Leydigiopsis
ornata (Daday, 1905), Notalona sculpta (Sars, 1901) were
located in the 3rd quarter. According to Fig. 4b, these species
dominated in lowlands and in the waters with low visibility
(negative correlation with Secchi disk visibility and alti-
tude). Chydorus cf. sphaericus was located between quad-
rants III and IV, which may indicate a weak relation with all
of the considered environmental variables. 

DISCUSSION

Thirty-one Cladocera species were identified in surface
sediments of the lakes studied. The number of species
found in individual lakes was generally small. In fact, the
maximum number of species found in lakes Atitlán and
Verde was only 12. Compared to European lakes, which
are normally inhabited by over 30 species (Dumont, 1994),
the species diversity of Cladocera was low in the consid-
ered region of Central America. These findings seems to
confirm that the number of Cladocera species found in one
waterbody in the Neotropics is in general low and usually
ranges from several to over a dozen species (Cuna et al.,
2014; Peréz et al., 2013; Zawisza et al., 2012). The differ-
ences in the number of cladoceran species between the con-

tinents may result from species competition (Feniova et al.,
2011; Shurin, 2000), fish or invertebrate predation (Korhola
and Rautio, 2001), as well as from different geological his-
tory of the area (Korhola and Rautio, 2001). 

The identified Cladocera included both littoral and
pelagic species. Planktonic species of Daphniidae and
Bosminidae were the most abundant. Remains of the
Daphnia pulex-group and the Daphnia longispina-group
were mainly associated to lakes located at altitudes above

Tab. 2. Number of occurrences, maximum and mean percentage
contribution of cladoceran species identified in the sediments of
the surveyed waterbodies. 

Group or species                                 n                  Max           Mean
                                                                                  (%)             (%) 

Bosmina longirostris                             9                 100.0           43.8
Bosmina (E.) coregoni                          8                  31.3            17.0
Bosmina (E.) longispina                      10                 80.0            36.2
Daphnia longispina-group                   8                  67.9            25.6
Daphnia pulex-group                            2                   7.0              7.0
Chydorus cf. sphaericus                      20                 89.5            19.7
Acroperus sp.                                        3                  22.2            12.2
Alona glabra                                        13                 28.6             6.7
Alona quadrangularis type                 12                 42.9            10.7
Alona cf. intermedia                             1                   5.6                -
Alona manueli                                      6                  40.0             9.6
Alona ossiani                                        8                  42.4            16.4
Alona rustica                                        1                   4.3                -
Alonella excisa                                      3                   4.3              2.5
Alonella nana                                        3                  13.6             7.1
Alonella pulchella                                4                   6.1              3.7
Anthalona verrucosa-group                  5                  14.3             5.9
Camptocercus dadayi                           6                  30.8            13.7
Disparalona dadayi                              2                  11.9              9.0
Dunhevedia sp.                                     3                  14.3             8.4
Ephemeroporus sp.                               2                  10.5             9.0
Euryalona sp.                                      11                 25.0             11.3
Graptoleberis testudinaria                   3                   3.6              2.5
Kurzia sp.                                             4                  13.3             7.2
Leberis sp.                                             5                  50.0            16.6
Leydigia louisi louisi                           10                 34.4             7.8
Leydigia acanthocercoides-group         6                  51.7            15.9
Leydigiopsis ornata                              2                  11.1              8.1
Notalona sculpta                                   7                  15.9             7.1
Pleuroxus denticulatus                         4                  11.1              6.4
Simocephalus sp.                                  2                  22.4            19.0
NRR1                                                    4                  31.8             11.4
NRR2                                                    1                  25.0               -
NRR3                                                    1                   4.3                -

NRR4                                                    1                   4.3                -
n, occurrence; NRR1-4, not recognize remains type 1 to 4; NRR2, prob-
ably headshield of Leberis sp.; NRR3, probably Coronatella monacan-
tha (Sars 1901) group. 

No
n c

om
me

rci
al 

us
e o

nly



158 M. Wojewódka et al.

Fig. 4. Results of multivariate statistical analysis applied to qualitative and quantitative results of species composition of subfossil
Cladocera, and environmental variables. A) Cluster analysis showing species grouping according to reciprocal correlation level; the
groups were distinguished for species that showed a high level of similarity >0.8; species composition of each group is presented in
Tab. 3; B) Canonical correspondence analysis (CCA) for species relative abundances and environmental variables; length and direction
of arrows show respectively the strength and direction of a certain environmental variable.

No
n c

om
me

rci
al 

us
e o

nly



Subfossil Cladocera of Central America 159

1000 m asl. CCA analysis showed the preference of the
Daphnia pulex-group for clear waters, such as lakes Ati-
tlán and Ipala, that showed transparency of 4.2 m and 3.6
m, respectively. Moreover, the contribution of Daphniidae
was negatively correlated with conductivity and was
higher in waters with conductivity ≤100 µS cm–1. The re-
sults show the preference of Daphniidae for highlands and
mountain areas and clean, low-mineralized and nutrient-
poor waters. Similar environmental preferences by Daph-
niidae were also suggest by Hart (2004) and by Mergeay
et al. (2005) for African lakes. The occurrence of Daph-
niidae species in lakes of Central American highlands was
also observed by Peréz et al. (2013). 

Bosminidae species were observed along the whole al-
titudinal gradient of the lake studied, from lowland, to high-
land and mountain lakes. CCA analysis revealed a negative
correlation of Eubosmina (Bosmina (E.) coregoni, Bosmina
(E.) longispina) and a positive correlation of Bosmina lon-
girostris with altitude. Eubosmina was found in six lowland
lakes, but also in two lakes located at >1000 m asl (Atitlán
and Ipala), in Lake Yojoa (639 m asl) and in Lake Apaste-
peque (509 m asl). This suggests that even though Eu-
bosmina prefered lakes located at a lower altitude, this is
probably not the main factor determining the distribution
of this genus. On the other hand, Eubosmina mainly oc-
curred in waters with high conductivity. These species were
dominant in highly mineralized (4520 µS cm–1) and brack-
ish (salinity=2.5‰) waters of Lakes Salpetén, as well as in
El Rosario and Lachuá, which showed conductivity values
>900 µS cm–1. The presence of Eubosmina, along with
Bosmina longirostris in Lake Salpetén, confirms the toler-
ance of Bosminidae to waters with high mineral content
and/or salinity (Aladin, 1991). The obtained results also
confirm the wide distribution of planktonic Daphniidae and
Bosminidae in Central America (Elías-Gutiérrez et al.,
1999; Korovchinsky, 2006). 

Chydoridae was the most species rich group of Clado-
cera in the sediment studied. Chydorus cf. sphaericus was
the most common species, as it was present in 20 out of
the 29 surveyed lakes located both in lowland and moun-
tain areas. The relative abundance of Chydorus cf. sphaer-
icus in the lake located at the highest altitude (Lake
Magdalena, 2863 m asl) was almost 90%. CCA analysis

showed that none of the measured environmental factors
had a significant effect on the distribution of the species
in the region studied of Central America. This confirmed
the ubiquity of Chydorus cf. sphaericus and its wide range
of tolerance to environmental and ecological conditions
(Flössner, 2000; Fryer, 1968; Korhola and Rautio, 2001;
Zawisza and Szeroczyńska, 2011). 

The genus Alona was mostly represented by Alona
glabra and the Alona quadrangularis type. These species
were found in 13 and 12 lakes, respectively, which were
characterized by different conductivity, oxygenation,
water visibility and altitude. The Alona quadrangularis
type was the unique littoral species found in Lake
Salpetén, which may indicate its tolerance to waters with
higher mineralization. The genus Alona was also repre-
sented by Alona manueli (Sinev and Zawisza, 2013) and
Alona ossiani, which, based on the cluster analysis, were
grouped together, along with Alonella nana (Baird, 1843),
Alonella pulchella (Herrick, 1884), the Anthalona verru-
cosa-group and Disparalona dadayi (Birge, 1910). The
multivariate analysis showed a negative correlation be-
tween Group 3 and lake conductivity, which most likely re-
sulted from the presence of species included in Group 3 in
lakes with conductivity ≤100 µS cm–1 and with medium
conductivity (204-331 µS cm–1). Lake Sacnab (412 µS
cm–1) and Lake Lachuá (906 µS/cm) are exceptions, and
two species from Group 3 (i.e., Alona ossiani and Alona
manueli, respectively) were present. 

Sinev (1998) showed that Alona ossiani belongs to the
Alona affinis complex, which is considered an olig-
otrophic species in Eurasia (Kamenik et al., 2007; Ko-
rhola and Rautio, 2001). On the other hand, Alona
manueli (an endemic species of Central America) and
Alonella pulchella were identified in oligotrophic lakes in
Mexico (Cuna et al., 2014; Sinev and Zawisza, 2013).
Moreover, this studied provided the first record of
Alonella pulchella in the tropical region, as previously this
species was recorded only in Canada (Korosi and Smol,
2012). The obtained results seem to confirm the prefer-
ence of these species for waters of low trophic status. The
negative correlation of Group 3 with the lake area may in-
dicate the preference of these species for smaller lakes
with well-developed littoral zone. In particular, Alonella

Tab. 3. Groups of species correlated with each other at the minimum level of 0.8 and distinguished based on the cluster analysis. 

Group Species composition

1 Daphnia pulex-group, Bosmina longirostris
2 Bosmina (E.) coregoni, Bosmina (E.) longispina, Alona glabra, Alona quadrangularis type, Leberis sp.
3 Alona manueli, Alona ossiani, Alonella nana, Alonella pulchella, Anthalona verrucosa-group, Disparalona dadayi
4 Euryalona sp., Kurzia sp., Leydigia acanthocercoides-group
5 Acroperus sp., Alona cf. intermedia, Camptocercus dadayi, Dunhevedia sp., Pleuroxus denticulatus
6 Alona rustica, Simocephalus sp.

No
n c

om
me

rci
al 

us
e o

nly



160 M. Wojewódka et al.

nana, was present in only three small lakes with area
≤0.1 km2. This suggests that, similarly to what observe in
Europe, Alonella nana may prefer small lakes (Fryer,
1968; Korosi and Smol, 2012). Furthermore, the presence
of Alona ossiani in mountain lakes may indicate a pro-
nounced tolerance toward the harsher climatic conditions
determined by altitude. The Eurasian species Alona affi-
nis, is considered to tolerate lower temperatures (Kamenik
et al., 2007; Locke and Sprules, 2000) and the results of
the present study indicated that Allona ossiani shows sim-
ilar characteristics in Central America. 

Leydigia louisi louisi and Leydigia acanthocercoides-
group played a significant role in Lake El Muchacho (al-
most 90%), which was characterized by a water
transparency of only 0.4 m. Moreover, these species were
relatively abundant (>10%) in other lakes with water trans-
parency below 0.9 m. The CCA analysis confirmed the
preference of these species for waters with lower Secchi
visibility. The presence of a species characteristic of Central
America such as Leydigiopsis ornata, which is a typical
benthic species (Van Damme and Dumont, 2010), was
recorded only in 2 shallow (<2.5 m water depth) and
scarcely transparent lakes, which were characterized
(Secchi depth <1 m). Graptoleberis testudinaria, which
was recorded in the sediments of Lakes Atitlán, Quexil and
Verde, preferred waters with higher Secchi visibility
(≥2.7 m), as showed also by the CCA analysis. However,
the results of the multivariate statistical analysis are highly
controversial, because Garptoleberis testudinaria is com-
monly considered to be associated with aquatic vegetation
(Fryer, 1968; Rybak and Błędzki, 2016). This suggests that
other factors not included in our analysis control the distri-
bution of this species. 

Pleuroxus denticulatus (Birge, 1879) remains were
identified in the sediments from four lakes (Chicabal,
Jucutuma, Ticamaya, Los Negritos) and its maximum
contribution to the Cladocera communities did not exceed
12% (average=6.4%). The presence of only one Pleuroxus
species in the lakes of the study area and its local distri-
bution confirms the rarity of this Cladocera group in the
tropical region (Korovchinsky, 2006).

Canonical correspondence analysis outlined altitude
as a key driver of the cladoceran fauna of in Central
America. Other significant environmental factors affect-
ing the Cladocera distribution were water transparency
and conductivity. Among the environmental variables
considered in this study, the less significant in affecting
the waterflea population in the study area was the concen-
tration of dissolved oxygen.

CONCLUSIONS

This study showed a relatively typical distribution of
Cladocera in lakes of the Neotropic region. In summary,

our results indicate that: i) the most common Cladocera
species of the study area was Chydorus cf. sphaericus,
likely in relation to the ability of this species to adapt to dif-
ferent ecological conditions; ii) planktonic species of Daph-
niidae and Bosminidae were the most abundant; iii)
Daphniidae species were recorded primarily in highland
lakes; iv) Eubosmina and the Alona quadrangularis type
showed pronounced tolerance to waters with high miner-
alizatiovn level; e) Alona ossiani, Alona manueli, Alonella
nana and Alonella pulchella were negatively correlated
with water conductivity; vi) the presence of Alonella nana
was restricted to lakes with a smaller surface area; vii)
Alona ossiani was an important subdominant species in
shallow lakes located at high altitudes; viii) among the en-
vironmental variables considered, altitude was a key driv-
ing factor for distribution of the Cladocera fauna whereas
dissolved water oxygenation was almost insignificant. 

These results provide basic information on Cladocera
communities in freshwater ecosystems of Central Amer-
ica, a region where research carried out so far is insuffi-
cient to provide robust taxonomical and ecological
information. 

ACKNOWLEDGMENTS

The study was funded by the Polish Ministry of Science
(Grant NCN 2014/13/B/ST10/02534) and the German Re-
search Foundation (DFG, SCHW 671/16-1). Furthermore,
scientific cooperation was supported by the Polish - Ger-
man governments (MNiSW-DAAD, 2016-2017). We
would like to thank Prof. A.Y. Sinev and Prof. A.A. Kotov
for their help with identification of specimens. Special
thanks is due to to Cuauhtémoc Ruiz (Instituto Tecnológico
de Chetumal), Ramón Beltrán (Centro Interdisciplinario de
Ciencias Marinas, Mexico), and Lisa Heise (Universidad
Autónoma de San Luis Potosí, Mexico) for their excellent
work on field. We also like to thank all people involved in
this work: Margarita Caballero (Instituto de Geofísica,
UNAM), Alexis Oliva and the team from the Asociación
de Municipios del Lago de Yojoa y su área de influencia
(AMUPROLAGO, Honduras), María Reneé Alvarez, Mar-
garita Palmieri, Eleonor de Tott (Universidad del Valle de
Guatemala, Guatemala), Consejo Nacional de Áreas Pro-
tegidas (CONAP, Guatemala), Néstor Herrera and Minis-
terio de Medio Ambiente (San Salvador, El Salvador).
CONACYT (Mexico) provided fellowship (218604,
405326) for the third and fourth authors. 

REFERENCES

Aladin NV, 1991. Salinity tolerance and morphology of the os-
moregulation organs in Cladocera with special reference to
Cladocera from the Aral Sea. Hydrobiologia 225:291-299.

Bigler C, Heiri C, Krskova R, Lotter AF, Sturm M, 2006. Dis-

No
n c

om
me

rci
al 

us
e o

nly



Subfossil Cladocera of Central America 161

tribution of diatoms, chironomids and Cladocera in surface
sediments of thirty mountain lakes in south-eastern Switzer-
land. Aquat. Sci. 68:154-171.

Bjerring R, 2007. Lake response to global change: nutrient and
climate effects using cladoceran (Crustacea) subfossils as
proxies. PhD Thesis, National Environmental Research In-
stitute University of Aarhus, Denmark.

Campanelli Mortari R, Henry R, 2016. Horizontal distribution
of Cladocera in a subtropical lake marginal to a river. J. Pa-
leolimnol. 75:109-120.

Cuna E, Zawisza E, Caballero M, Ruiz-Fernández AC, Lozano-
García S, Alcocer J, 2014. Environmental impacts of Little
Ice Age cooling in central Mexico recorded in the sediments
of a tropical alpine lake. J. Paleolimnol. 51:1-14.

Dumont HJ, 1994. On the diversity of the Cladocera in the trop-
ics. Hydrobiologia 272:27-38.

Elías-Gutiérrez M, Ciros-Pérez J, Suárez-Morales E, Silva-Bri-
ano M, 1999. The freshwater Cladocera orders Ctenopoda
and Anomopoda) of Mexico, with comments on selected
taxa. Crustacean 72:171-186.

Elías-Gutiérrez M, Kotov AA, Garfias-Espejo T, 2006. Clado-
cera (Crustacea: Ctenopoda: Anomopoda) from southern
Mexico, Belize and Northern Guatemala.Zootaxa 119:
1-27.

Elías-Gutiérrez M, Suárez-Morales E, Gutiérrez-Aguirre MA,
Silva-Briano M, Granados-Ramírez JG, Garfias-Espejo T,
2008. [Cladocera y Copepoda de las aguas continentales de
México].[Book in Spanish]. Universidad Nacional Autónoma
de México: 322 pp.

Feniova IY, Razlutsky VI, Palash AL, 2011. Temperature effects
of interspecies competition between Cladoceran species in
experimental conditions. Inland Water. Biol. 4:65-71.

Flössner D, 2000. [Die Haplopoda und Cladocera (ohne
Bosminidae) Mitteleuropas].[Book in German]. Backhuys
Publisher, Leiden: 428 pp.

Frey DG, 1986. Cladocera analysis, p. 667-692. In: B.E.
Berglund (ed.), Handbook of Holocene palaeoecology and
palaeohydrology. J. Wiley & Sons, Chichester.

Frolova L, Nazarova L, Pestryakova L, Herzschuh U, 2014.
Subfossil Cladocera from surface sediment in thermokarst
lakes in northern Siberia, Russia, in relation to limnological
and climatic variables, J. Paleolimnol. 52:107-119.

Fryer G, 1968. Evolution and adaptive radiation in the Chydoridae
(Crustacea: Cladocera): a study in comparative functional mor-
phology and ecology. Philos. T. R. Soc. B 254:221-385.

Green J, 1995. Altitudinal distribution of tropical planktonic
Cladocera. Hydrobiologia 307:75-84.

Hammer O, Harper DAT, Ryan PD, 2001. Past Paleontological
Statistics Software Package for Education and Data Analy-
sis. Palaeontol. Electron. 4:4-9.

Hart RC, 2004. Cladoceran periodicity patterns in relation to se-
lected environmental factors in two cascading warm-water
reservoirs over a decade, Hydrobiologia 526:99-117.

Hudec I, 2000. Subgeneric differentiation within Kurzia (Crus-
tacea: Anomopoda: Chydoridae) and a new species from
Central America. Hydrobiologia 421:165-178.

Jeppesen E, Leavitt P, De Meester L, Peder Jensen J, 2001.
Functional ecology and paleolimnology: using Cladocera
remains to reconstruction anthropogenic impact. Trends
Ecol. Evol. 16:191-198.

Juggins S, 2005. New features in C2 version 1.4. University of
Newcastle.

Juggins S, 2007. User guide C2 Software for ecological and pa-
leoecological data analysis and visualisation User guide Ver-
sion 1.5. University of Newcastle.

Kamenik C, Szeroczyńska K, Schmidt R, 2007. Relationships
among recent Alpine Cladocera remains and their environ-
ment: implications for climate-changes studies. Hydrobiolo-
gia 594:33-46.

Karmalkar AV, Bradley RS, Diaz HF, 2011. Climate change in
Central America and Mexico: regional climate model valida-
tion and climate change projections. Clim. Dyn. 37:605-629.

Kienast F, Wetterich S, Kuzmina S, Schirrmeister L, Andreev
AA, Tarasov P, Nazarova L, Kossler A, Frolova L, Kunitsky
VV, 2011. Paleontological records prove boreal woodland
under dry inland climate at today’s Arctic coast in Beringia
during the last interglacial. Quat. Sci. Rev. 30:2134-2159.

Kotov AA, Elías-Gutiérrez M, Guadalupe Nieto M, 2003. Ley-
digia louisi louisi Jenkin, 1934 in the Neotropics, L. louisi
mexicana n. subsp. in Central Mexican highlands. Hydrobi-
ologia 510:239-255.

Kotov AA, 2009. A revision of Leydigia Kurz, 1875
(Anomopoda, Cladocera, Branchiopoda), and subgeneric
differentiation within the genus. Magnolia Press: 84 pp.

Korhola A, 1999. Distribution patterns of Cladocera in subarctic
Fennoscandian lakes and their potential in environmental re-
construction. Ecography 22:357-373.

Korhola A, Rautio M, 2001. Cladocera and other Branchiopod
crustaceans, p. 5-41. In: J.P. Smol, J.B. Birks and W.M. Last
(eds.), Tracking environmental change using lake sediments.
Zoological indicators. Kluwer Academic, Dordrecht.

Korosi JB, Smol JP, 2012. An illustrated guide to the identifica-
tion of cladoceran subfossil from lake sediments in north-
eastern North America: part 2 - the Chydoridae, J.
Paleolimnol. 48:587-622.

Korovchinsky NM, 2006. The cladocera (Crustacea: Bra-
chiopoda) as a relict group. Zool. J. Linn. Soc. 147:109-124.

Locke A, Sprules WG, 2000. Effects of acidic pH and phy-
tolankton on survival and condition of Bosmina longirostris
and Daphnia pulex. Hydrobiologia. 437:187-196.

Luoto T, Sarmaja-Korjonen K, Nevalainen L, Kauppila T, 2009. A
700 year record of temperature and nutrient changes in a small
eutropied lake in southern Finland. Holocene 19:1063-1072.

Mergeay J, Verschuren D, De Meester L, 2005. Daphnia species
diversity in Kenya, and a key to the identification of their
ephippia. Hydrobiologia 542:261-274.

Mirosław-Grabowska J, Zawisza E, 2014. Late Glaciale early
Holocene environmental changes in Charzykowskie Lake
(northern Poland) based on oxygen and carbon isotopes and
Cladocera data. Quat. Int. 328-329:156-166.

Nováková K, Van Hardenbroek M, van der Knaap WO, 2013.
Response of subfossil Cladocera in Gerzensee (Swiss
Plateau) to early Late Glacial environmental change. Palaeo-
geogr. Palaeoclimatol. Palaeoecol. 391:84-89.

Padhye SM, Kotov AA, Dahanukar N, Dumont HJ, 2016. Bio-
geography of the ‘water flea’ Daphnia O.F. Müller, (Crustacea:
Branchiopoda: Anomopoda) on the Indian subcontinent. J.
Limnol. 75:1476.

Paterson MJ, 1994. Paleolimnological reconstruction of recent
changes in assemblages of Cladocera from acidified lakes

No
n c

om
me

rci
al 

us
e o

nly



162 M. Wojewódka et al.

in the Adirondack Mountains (New York). J. Paleolimnol.
11:189-200.

Peréz L, Lorenschat J, Massaferro J, Pailles Ch, Sylvestre F, Holl-
wedel W, Brandorff GO, Brenner M, Islebe G, Lozano MS,
Scharf B, Schwalb A, 2013. Bioindicators of climate and
trophic state in lowland and highland aquatic ecosystems of
the Northern Neotropics. Rev. Biol. Trop. 61:603-644.

Rajapaksa R, Fernando CH, 1987. Redescription and assignment
of Alona globulosa Daday, 1898 to a new genus Notoalona
and a description of Notoalona freyi sp.nov. Hydrobiologia
144:131-153.

Rey J, Vasquez E, 1986. [Contribution à la connaissance des
Cladocères néotropicaux: redescription de Leydigiopsis or-
nata Daday, 1905. (Crustacea, Cladocera)].[Article in
Fench]. Ann. Limnol. 22:169-176.

Rybak JI, Błędzik LA, 2016. Freshwater Crustacean zooplank-
ton of Europe. Springer: 918 pp.

Schmidt R, Müller J, Drescher-Schneider R, Krisai R, Sze-
roczyńska K, Barić A, 2000. Changes in lake level and tro-
phy at Lake Vrana, a large karstic lake on the Islanad of Cres
(Croatia), with respect to palaeoclimate and anthropogenic
impacts during the last approx. 16,000 years. J. Limnol.
59:113-130.

Shurin JB, 2000. Dispersal limitation, invasion resistance, and
the structure of pond zooplankton comunities. Ecology
81:3074-3086.

Sinev AY, 1998. Alona ossiani sp. n., a new species of the Alona
affinis complex from Brazil, deriving from the collection of
G.O. Sars (Anomopoda Chydoridae). Arthropoda Sel.
7:103-110.

Sinev AY, 2001. Redescripton of Alona glabra Sars, 1901, a
South American species of the pulchella group (Bran-
chiopoda: Anomopoda: Chydoridae). Arthropoda Sel.
10:273-280.

Sinev A, Zawisza E, 2013. Comments on cladocerans of crater
lakes of the Nevado de Toluca Volcano (Central Mexico),
with the description of a new species, Alona manueli sp.
Nov. Zootaxa 3647:390-400.

Sinev A Y, 2015a. Revision of the pulchella-group of Alona s.
lato leads to its translocation to Ovalona Van Damme et Du-
mont, 2008 (Branchiopoda: Anomopoda: Chydoridae).
Zootaxa 4044:451-492.

Sinev AY, 2015b. Morphology and phylogenetic position of
three species of genus Camptocercus Baird, 1843 (Clado-
cera: Anomopoda: Chydoridae). Zootaxa 4040:169-186.

Sinev AY, Dumont HG, 2016. Revision of the costata-group of
Alona s. lato (Cladocera: Anomopoda: Chydoridae) con-
firms its generic status. Eur. J.Taxon. 223:1-38.

Sinev AY, Silva-Briano M, 2012. Cladocerans of genus Alona
Baird, 1843 (Cladocera: Snomopoda: Chydoridae) and re-
lated genera from Aguascalientes State, Mexico. Zootaxa.
3569:1-24.

Sweetman JN, Rühland KM, Smol JP, 2010. Environmental and
spatial factors influencing the distribution of cladocerans in
lakes across the central Canadian Arctic treeline region. J.
Limnol. 69:76-87.

Szeroczyńska K,1991. Impact of prehistoric settlements on the
Cladocera in the sediments of Lakes Suszek, Błędowo, and
Skrzetuszewskie. Hydrobiologia 225:105-114.

Szeroczyńska K, Tatur A, Weckstrom J, Gąsiorowski M, No-
ryśkiewicz A, Sienkiewicz E, 2007. Holocene environmen-
tal history in northwest Finnish Lapland reflected in the
multi-proxy record of small subartic lake. J. Paleolimn.
38:25-47.

Szeroczyńska K, Zawisza E, 2011a. Subfossil faunal and floral
remains (Cladocera, Pediastrum) in two northern Lobelia
lakes in Finland. Knowl. Manag. Aquat. Ect. 403:09.

Szeroczyńska K, Zawisza E, 2011b. Records of the 8200 cal BP
cold event reflected in the compositionof subfossil Clado-
cera in the sediments of three lakes in Poland. Quat. Int.
233:185-193.

Szeroczyńska K, Zawisza E, Wojewódka M, 2015. Initial time
of two high altitude crater lakes (Nevado De Toluca, Central
Mexico) recorded in subfossil Cladocera. Stud. Quatern.
32:109-116.

Taylor MA, Alfaro EJ, 2005. Climate of Central America and
the Caribbean, p. 181-189. In: J.E. Oliver (ed.), Encyclope-
dia of world climatology. Springer, Dordrecht.

ter Braak CJF, 1986. Canonical correspondence analysis: a new
eigenvector technique for multivariate direct gradient analy-
sis. Ecology 67:1167-1179.

Van Damme KV, Dumont HJ, 2010. Cladocera of the Lençóis
Maranhenses (NE - Brazil): faunal composition and a reap-
praisal of Sars’ method. Braz. J. Biol. 70:755-779.

Van Damme KV, Sinev AY, Dumont HJ, 2011. Separation of An-
thalona gen.n. from Alona Baird, 1843 (Branchiopoda:
Cladocera: Anomopoda): morphology and evolution of
scraping stenothermic alonines. Zootaxa 2875:1-64.

Zawiska I, Słowiński M, Correa-Metrio A, Obremska M, Luoto
T, Nevalainen L, Woszczyk M, Milecka K, 2015. The re-
sponse of a shallow lake and its catchment to Late Glacial
climate changes - A case study from eastern Poland. Catena
126:1-10.

Zawisza E, Szeroczyńska K, 2011. Cladocera species composi-
tion in lakes in the area of the Hornsund Fjord (Southern
Spitsbergen)-Preliminary results, Knowl. Manag. Aquat. Ec.
402:04.

Zawisza E, Caballero M, Ruiz-Fernandez C, 2012. 500 years of
ecological changes recorded in subfossil Cladocera in a
high-altitude tropical lake Lago de La Luna, Central Mex-
ico. Stud. Quatern. 29:23-29.

Zawisza E, Correa-Metrio A, Caballero M, Lozano S, Sze-
roczyńska K, 2014. Paleoecology of tropical Lake Zirhauen
(Western Mexico) recorded in Cladocera remains. Abstract
book of the 8th Shallow lakes Conf., Antalya, Turkey.

Zawisza E, Cuna E, Caballero M, Ruiz-Fernandez AC, Sze-
roczyńska K, Woszczyk M, Zawiska I, 2016. Environmental
changes during the last millennium recorded in subfossil
Cladocera, diatoms and sediment geochemistry from Lake
El Sol (Central Mexico). Geol. Q. (in press).

No
n c

om
me

rci
al 

us
e o

nly