The dynamics of social networks among female Asian elephants
© de Silva et al; licensee BioMed Central Ltd. 2011
Received: 23 March 2011
Accepted: 27 July 2011
Published: 27 July 2011
Patterns in the association of individuals can shed light on the underlying conditions and processes that shape societies. Here we characterize patterns of association in a population of wild Asian Elephants at Uda Walawe National Park in Sri Lanka. We observed 286 individually-identified adult female elephants over 20 months and examined their social dynamics at three levels of organization: pairs of individuals (dyads), small sets of direct companions (ego-networks), and the population level (complete networks).
Corroborating previous studies of this and other Asian elephant populations, we find that the sizes of elephant groups observed in the field on any particular day are typically small and that rates of association are low. In contrast to earlier studies, our longitudinal observations reveal that individuals form larger social units that can be remarkably stable across years while associations among such units change across seasons. Association rates tend to peak in dry seasons as opposed to wet seasons, with some cyclicity at the level of dyads. In addition, we find that individuals vary substantially in their fidelity to companions. At the ego-network level, we find that despite these fluctuations, individuals associate with a pool of long-term companions. At the population level, social networks do not exhibit any clear seasonal structure or hierarchical stratification.
This detailed longitudinal study reveals different social dynamics at different levels of organization. Taken together, these results demonstrate that low association rates, seemingly small group sizes, and fission-fusion grouping behavior mask hidden stability in the extensive and fluid social affiliations in this population of Asian elephants.
Determining the ecological conditions that shape group formation and social structures is a prerequisite for understanding social evolution [1–4]. Studies in numerous group-living species that follow individuals longitudinally over multiple years find that relationships among individuals change both quantitatively and qualitatively over time [5–9]. This is especially true of societies structured by fission-fusion processes, in which associations among individuals may change over time scales ranging from hours to months [10–13]. Patterns in these dynamics may shed light on the underlying ecological conditions that drive the behavior of individuals [2, 3, 14, 15]. In this paper, we examine social dynamics in a population of Asian elephants.
The longevity and cognitive sophistication of elephants make them potentially capable of maintaining complex social relationships [16–19]. The Asian elephant (Elephas maximus), African savannah elephant (Loxodonta africana) and African forest elephant, (Loxodonta africana cyclotis or Loxodonta cyclotis) are the only living members of the proboscidean clade [20–23]. Among adult African and Asian elephants, females and calves form the basis of social units. In African savannah elephants, relatives associate most closely, but multiple social units associate periodically to form hierarchically stratified 'multi-tiered' societies with fission-fusion dynamics [6, 24, 25]. Relatedness decreases at higher order 'tiers' . Much less is known about the social organization of African forest elephants, but it appears they generally tend to form smaller social groups than savannah elephants [27, 28].
Surprisingly little is known of the social behavior of Asian elephants in the wild, despite their long history of co-habitation with people. While it has been commonly assumed that Asian and African elephants behave similarly , independent field studies in India and Sri Lanka suggest that adult females associate only with maternal relatives [30, 31] and also tend to form smaller groups than African savannah elephants [30, 32]. Early studies report that Asian elephants form 'loose' associations with one another, implying unstable social affiliations [32, 33]. A more recent study of Asian elephants at Uda Walawe and Yala National Parks in southern Sri Lanka found that matrilineal kin in fact associated together only 18-20% of time, based on observation and telemetry data, leading to the conclusion that individuals from different matrilines are unlikely to associate at all, and that inter-group transfer of females most likely does not occur . It was proposed that family groups fission into daughter groups which then become largely independent of one another .
In this paper we investigate how local ecological conditions influence the social dynamics of adult female Asian elephants by tracking associations within a single wild population that inhabits a highly seasonal environment in Sri Lanka. There are at least three distinct temporal patterns of associations that we could potentially observe: random, cyclic, or stable. In the light of previous findings , we expect to be able to reject the null hypothesis of random associations. We expect to observe cyclicity in associations if seasonal resource abundance governs the degree to which individuals associate. Alternately, individuals may form stable associations that are independent of environmental conditions. We take cross-sectional and longitudinal views to determine which patterns emerge at different levels of organization within the society. First, we characterize how the strength of the relationship between pairs of individuals changes over time. We then quantify how close companions of individuals, their putative 'ego-networks,' change over time. Finally, we take a bird's eye view of the population to examine whether this Asian elephant society has a tier structure akin to that of African elephants.
Glossary of terms.
A set of individuals observed in the field moving, resting, or interacting non-aggressively within an approximately 500 m radius of one another.
A set of individuals that repeatedly associate together such that they form a distinct unit which is revealed by an objective analytical clustering method.
The set of individuals that the central individual, called 'ego', is directly connected to.
A general term describing sets of socially affiliated individuals. Social units can form at different levels. For example, an individual's ego-network, a social unit at one level, may be embedded in a larger network, constituting a social unit of a higher level.
Jan - Apr
May - Sept
Oct - Dec
May - Sept
Oct - Dec
No. of individuals
Mean group size1
Max. group size1
% Non-zero SRI
Mean non-zero SRI
Number of clusters2
Mean cluster size2
Max cluster size2
Structure at the level of dyads
Correlation between matched SRI matrices.
T1 vs. D1
T1 vs. W1
D1 vs. W1
D1 vs. D2
D2 vs. W2
W1 vs. W2
Temporal SRI trajectories for the core individuals.
Number of dyads
Number of individuals (percent)
D1 and D2
All ego-network statistics except 'Density' tend to have lower values in wet seasons than in dry seasons (Table 5), indicating that ego-networks are larger in dry seasons than in wet seasons. The measure 'Density' is marginally lower in dry seasons, implying that the ego-networks are less interconnected in dry seasons than in wet seasons. Because ego-network measures for different individuals are non-independent, we cannot directly test for differences between the distributions of ego-network measures in different seasons. Instead, we test for differences in the bootstrap distributions of ego-network statistics . The bootstrap means of ego-network statistics are reported in Table 5. Bootstrap distributions for all ego-network statistics differ significantly between seasons. We note however that these distributions for all statistics in all seasons are systematically shifted with respect to the observed values (Table 5). This suggests that the bootstrap method we used produces biased distributions of ego-network statistics, but we are not aware of a better method for statistical comparison. We therefore must be cautious with the interpretation of Table 5.
The ego-network measures are sensitive to the number and arrangement of companions but ignore their identities. However, visualizations of ego-networks demonstrate that while a subject's direct companions do change over time, she has a few that are almost always present; even those that are not present continuously may have been companions in previous seasons (Figure 5). Thus, individuals maintain long-term relationships with others even though they may be apart for one or several seasons and have low SRI values.
First, we note that the social network formed by individuals in this population is extensive and well-connected (Figure 5). Moreover, in all seasons, a typical cluster detected by the Girvan-Newman algorithm is 2 to 3 times larger than a typical group encountered by observers in the field (Figure 6 and Table 2). These findings are unexpected in the light of previous studies which have suggested that Asian elephants do not display extensive social affiliations [30, 33].
Next, we examine how the structure of the social network changes over time. For each season, we plot the network structure curves (Figures 7 and 8; see Methods for detailed explanation), which show how ties of various strengths are distributed in a manner that is easier to digest than the social network diagrams depicted in Figure 5 (the effect of the thresholding procedure on the T1 network is illustrated in Additional file 5: Figure S5). A network generated from randomized data shows no distinct clusters (see Additional file 6: Figure S6). As a result, the corresponding network structure curve remains flat until the SRI threshold reaches the mean tie strength in the network, at which point the curve rapidly declines to zero (Figure 8). In contrast, the network structure curves for real data have a distinct peak. This is probably a consequence of the fact that both extremely weak ties and extremely strong ties can be influenced by sampling effects, but intermediate values clearly distinguish core social units. Although all seasons have 10-20 well-defined clusters at the SRI threshold of zero, network structure curves for each season behave differently at subsequent SRI threshold values. In particular, the slope of the curve changes only once (from positive to negative) in seasons D1, D2 and W2, but changes twice (from zero to positive and from positive to negative) in seasons T1 and W1 (Figure 8). This suggests that the network structures in each season are different (also visible in Figure 5). Although the network structure curves are different across seasons, there appear to be no patterns characteristic of dry or wet seasons per se. Visual inspection of networks shows that, while many clusters maintain their integrity across seasons, some individuals transfer between clusters, and the connections among clusters change (Figure 5).
We have presented the first detailed quantitative characterization of social organization in an Asian elephant population, at multiple levels of organization and across ecological timescales. We asked whether associations are random, stable, or cyclic. We find that the answer depends on the timescale and level of organization. Most ties are weak (SRI values below 0.3) compared those of African savannah elephants where typical association rates are above 0.6 [6, 13, 25]. Despite the overall weakness of ties, most individuals have a few strong ties (SRI values exceeding 0.3) as well as a few consistent ties (maintained over several seasons) with some of their associates. Individuals do not mix randomly within the population, nor are they always with the same companions, but rather they shuffle amongst a subset of preferred companions. All individuals engage in temporary associations, especially during dry seasons. Some also associate to a greater degree in dry seasons, forming cyclic associations. This cyclicity is evident at the level of dyadic associations (Figure 1), but not at the ego-network or population levels.
Our results suggest a view of Asian elephant social structure that is different from what has been described in the literature before. Earlier studies reported small group sizes typically consisting of less than five adult females [30, 33], a number comparable to the group sizes observed in this study. The number of individuals with the same mitochondrial haplotype ranged from 3 to 12 , which could be taken as matrilineal family sizes. However, families and observer-defined groups are only relevant structural units insofar as they relate to actual social units, expressed through the animals' behavior. The latter are revealed by a quantitative analysis of long-term association data rather than genetics [6, 26]. It took up to two years to obtain repeated observations of all members of some social units in this study (Figure 5 and unpublished data from 2005-2006). In this paper, we used the Girvan-Newman procedure to detect such social units and found that they are in fact much larger than either the group or family sizes reported in previous studies (Table 2, Figures 5 and 6). Such units are more stable across years than individuals' immediate companions. Informal photographic records suggest that some of the individuals in this study have associated since at least 2001 (unpublished data), which supports the existence of long-term stable associations. Since Asian elephants are capable of communicating both chemically [29, 40] and acoustically [41, 42] at distances humans find difficult to observe, they may be aware of the their associates' locations even when the latter are beyond the visual range of human observers. Indeed, vision is not the preferred mode of perception for elephants, as we see individuals track precise paths taken by others using scent even when both parties are plainly visible to humans. Moreover, outside protected areas, elephants are largely nocturnal (personal observations). 'Groups' of Asian elephants are not unlike those formed by African elephants . As with other animal societies that exhibit fission-fusion dynamics, such as that of chimpanzees [43, 44], the social organization of a highly mobile species like Asian elephants is not fully evident without systematic and prolonged observations, particularly in areas where visibility is restricted. Long-term observations of other Asian elephant populations would be extremely useful to corroborate this finding.
Previous authors also concluded that associations among different families were highly unlikely as associations even among family members appeared infrequent . However, these conclusions were drawn from extremely small sample sizes (for instance, only 1 mtDNA haplotype from Uda Walawe and few repeat observations). Our results suggest otherwise. While the association rate of 18-20% reported in a previous study  is roughly analogous to the median SRI value of found in our study (Table 2, Figure 6), we do find reliable SRI values that range as high as 1 (see Additional file 1: Figure S1). Moreover, there is at least some transfer of individuals between social units across seasons (Figure 5). However, the population-level social network structure does not appear to exhibit any clear seasonal patterns (Figure 8). It is not clear whether individuals form hierarchical social 'tiers', such as those observed in African savannah elephants, which form higher-order associations among multiple families in wet seasons [6, 24, 25, 45]. Among savannah elephants, relatedness decreases at higher-order tiers of association, where associations are weaker than 0.6 [13, 26, 46], and are mediated by intra- and inter-group dominance interactions [47, 48]. The network structure curve for Asian elephants peaks near an SRI threshold of 0.3 (Figure 8). Although there is no consistent social stratification in this Asian elephant population, it is possible that the clusters prior to the peak have a lower degree of relatedness than the clusters that follow it. A detailed genetic study of this population examining the hypothesis above, with larger sample sizes than previously obtained, would also be illuminating and is planned in the future.
There is much variation in individuals' long-term fidelity to companions (Figures 3 and 4). For instance, Kamala (KAM) and Kanthi (KAN) were two mature females who appeared close to the same age and were nearly always together (Figure 5). The so-called 'K' unit (Kamala, Kanthi, Karin, Kavitha and Kalyani) almost always contained every member whenever it was seen although they also interacted with others to form a larger cluster. On the other hand, individuals like '471,' also part of a large cluster, had few stable companions (Figure 5). The social placement of a few other females remained unresolved despite numerous repeat observations. The fitness consequences of these different social strategies remain to be seen. Moreover, while it is widely assumed that Asian elephants, like African elephants, form strongly bonded family groups centered around matriarchs [29, 40], the apparent variation and fluidity in social preferences shown in this study would seem to question such a characterization.
It is intriguing that social dynamics differ depending on the level of analysis - the bottom-most (dyadic) and top-most (population) levels of organization exhibit a greater degree of instability than the intermediate level (social units and long-term ego-networks). Uncovering the ecological basis for observed patterns would require separate investigations and hypotheses at each level. Preferences for one another shown by some pairs of individuals might depend on reproductive state (e.g. those with similarly-aged calves), while social units may differ in their strategies depending on whether they are seasonal inhabitants of the park or residents. For instance, killer whales of the same species exhibit different social strategies depending on whether they are resident or transient, in accordance with the associated feeding ecology .
Our analyses are based on association index data, calculated from observations of individuals in the field. There are at least two sources of uncertainty associated with this type of data. First, we expect that some variation comes from the fact that we observe different sets of individuals in different seasons. Our observation area is largely constrained by the road network inside Uda Walawe National Park. While we expect to have a reliable observation record for individuals whose range strongly overlaps with this area, we also observe individuals that presumably move into this area only periodically . Such individuals, being farther away from the centers of their home ranges, might exhibit different behavioral and social patterns than individuals residing more centrally, thus introducing additional noise in our data. The second source of variation stems simply from relatively low counts of events in some seasons for some individuals. To minimize the first type of noise, we have constrained most of our analyses to the so-called resident individuals, i.e. those that we have consistently observed every season. To minimize the second type of noise, we have constrained some of the analyses even further, to individuals that have been seen at least 30 times. Nevertheless, we can extrapolate the conclusions drawn from these analyses to less sampled individuals in the population, since such individuals are not ostensibly different from those sampled thoroughly. One does however need to keep in mind that we describe the behavior of individuals that are close to the center of their home range.
One of our surprising findings is that the elephants at UWNP tend to form a greater proportion of strong ties in dry seasons than in wet seasons. This suggests that aggregation may be more advantageous in the former, perhaps for accessing and protecting scarce resources. This hypothesis remains to be tested with additional seasons of data and behavioral studies. While direct behavioral evidence of resource defense among adult females is rare, we have observed competition over water and mud, dominance interactions when unfamiliar individuals or social units meet, as well as the vocal and physical displacement of one social unit by another . Resource monopolization may more often take the form of competitive exclusion rather than confrontation, in which acoustic and chemical signals facilitate social cohesion as well as avoidance despite the seeming fluidity of associations. Herbivores must balance intraspecific resource competition against potential anti-predator benefits [3, 50, 51]. Among artiodactyles, gregariousness is an anti-predator adaptation seen in species inhabiting open environments . African savannah elephants likewise may be more gregarious than Asian elephants because they typically inhabit more open environments, and also encounter predators other than humans [53–55]. Interestingly, in drier regions of Sri Lanka just a few kilometers east of the study site, elephants are reported to aggregate in wet seasons rather than dry seasons , a similarity to African savannah elephants that could be ecologically driven. Similar longitudinal studies in other Asian elephant populations, especially those in India, where there is likely to be greater variation in habitat quality and home range sizes  would be of great interest. More data are also needed on African elephants occupying various habitats including desert and forest environments, the latter being more similar to those of many Asian populations. It is possible that societies in general and elephant societies in particular are more flexible and responsive to environmental pressures than generally conceded.
Associations among female Asian elephants can be characterized as fission-fusion. Patterns of association differ across levels of organization and ecological timescales. Individuals are found by observers in small groups whose composition changes on the timescale of days. However, most individuals belong to relatively large social units which can be revealed by quantitative analysis of longitudinal data. Such units generally maintain their cohesiveness over longer timescales, although there is some transfer of individuals between units. Social units may fission or fuse without discernible seasonal patterns or clear hierarchical stratification into social 'tiers' at the level of the population. Individuals vary highly in the number of stable companions they maintain over time, and some repeatedly associate in dry seasons. Their companions tend to form a pool of long-term associates. Interesting future directions include examining the relatedness among social units, the ecological basis of observed dynamics at each level of organization, and the fitness consequences of the seemingly different social strategies employed by individuals.
Uda Walawe National Park (UWNP), Sri Lanka, is located between latitudes 6° 25' - 6° 34' N and longitudes 80° 46' - 81° 00' E, at an average altitude of 118 m above sea level. It encompasses 308 km2 around the catchment of the Uda Walawe reservoir. The study area comprises approximately 1/3 of this area, which includes tall grassland, dense scrub, riparian forest, secondary forest, a permanent river, seasonal streams, and other water sources. It has a highly predictable pattern of rainfall with two monsoons per year, which occur in March through April and in October through December [58, 59]. There appear to be no non-human predators of elephants at this location. The Sri Lankan subspecies of Asian lion (Panthera leo sinhaleyus) was extinct prior to the colonization of the island by humans . The leopard (Panthera pardus cotiya) is the current top land predator in Uda Walawe, but there is no evidence that it hunts elephants. The only other large predator is the freshwater crocodile (Crocodylus palustris), but predation on elephants has not been documented. The greatest threat to elephants both historically and currently is human activity, but disturbance within the park is minimal. Tourism has also led to the elephant population becoming well-habituated to people in vehicles.
The data presented here span 259 field days (twenty months) from years 2007 and 2008, two or three days per week on average except between January and April 2008, when UWNP was temporarily closed due to political unrest. We typically entered the park between 0600-0700 h (sunrise), remaining continuously inside until 1730-1830 h (sunset). Driving routes were varied such that all accessible parts of the park were covered in a week. Locations where animals were closest to the road were marked on a hand-held GPS unit. Temperature, humidity and wind were recorded at least three times per day with a Kestrel™ pocket weather station. Rainfall (mm.) was recorded daily using a standard U.S. Weather Bureau rain gauge.
Individuals were identified photographically and catalogued for two years preceding the study period. All individuals were given numbers; the most frequently seen were also given names. A previous study of Asian elephants considered individuals within 100 m of one another to be associated  whereas studies of African elephants have used a distance of up to 500 m to define aggregations [6, 25, 61]. We considered all individuals within visual range of the observer and up to 500 m of one another who moved, rested, shared resources (mud, mineral wells, trees) to be a single aggregation, or group. Occasionally individuals showed affiliative vocal or tactile behavior such as growling and bodily rubbing  but such interactions were uncommon and not required for individuals to be considered part of the same group. Multiple groups occasionally shared water without interaction. The term 'group' here carries no implication of social history or permanence (see also Table 1). Individuals from multiple groups which initially co-occurred in space or even passed through one another, were not counted as associated unless they actively moved together. It was possible to spend several hours with a single group. We recorded identities of known individuals and counted the number of individuals in five size-based age classes . Unidentified individuals were counted, but excluded from analyses.
We examined only relationships among adult females, as most sub-adults and juveniles were not identified individually. The strength of one individual's bond with another individual over the course of a season was quantified in terms of their association index. We used the Simple Ratio Index or SRI [6, 35, 36, 62], which is a symmetric measure that shows the proportion of time two individuals spent with each other. Before computing the SRI, we performed data aggregation which consisted of (a) partitioning the data into day-long sampling intervals; (b) identifying individuals that associated with each other in a given sampling interval by merging those groups that shared at least one identified individual; all individuals within such groups were then, by definition, 'associated'; and (c) excluding individuals that were observed only once in a season and only alone. Step (b) was performed in order to exclude potentially non-independent observations of the same group within the same sampling period. After data aggregation, we compute the association index for each pair of individuals A and B as SRI = X AB /(X t - X n ) where X AB is the number of times individuals A and B were observed together, X t is the total number of observations, and X n is the number of observations in which neither A nor B were observed . We also compute a measure of uncertainty of the SRI value for each dyad (see Additional file 7: Supplementary Text).
Association data were partitioned according to season. Months that had a total rainfall higher than the two year monthly average of 120 cm were designated as 'wet' months and those that had less were designated as 'dry', consistent with the monsoon cycle . May-September constitute the 'Dry season' and October-December constitute the 'Wet season' according to this classification. January-April, with two wet months followed by two dry months, were considered a 'Transitional' period rather than divided into dry and wet periods since two month periods provided insufficient data for analysis. We refer to seasons as T1, D1, W1, D2 and W2 (Table 2).
To reduce variance in our data that arises from differences in the identities of individuals observed in different seasons, we constrain some of the analyses below to the set of 88 individuals that were observed in all seasons (either more than once or in association with other individuals). We refer to these individuals as 'residents'. To further reduce noise in the data due to low number of sightings, we constrain some of the analyses to the 51 residents that were observed at least 30 times throughout the entire study period. We refer to these individuals as the 'core individuals'. The uncertainty in the estimates of association indices for such individuals is generally below 10% for all seasons (see Additional file 1: Figure S1).
To test the data against the null hypothesis that associations within a season are random, we permute the seasonally partitioned datasets so that the number of sightings for each individual and the distribution of group sizes within the season are preserved. We use the 'fill' rather than the 'swap' method to generate 1000 permutations per season [63–65], with the average SRI value as the test statistic. To speed up computations, we partition the dry seasons into two overlapping three-month periods (May-July and July-September). We use some of the random datasets generated by this procedure in other analyses.
We examine the stability of associations among pairs of individuals across seasons in a few different ways. First, we test for correlations between matched SRI matrices across pairs of seasons using the Mantel test [63, 66, 67]. As not all individuals are seen in all seasons, by 'matched' matrices we mean that for each pair of seasons we test correlations only across the subset of individuals seen in both seasons. We used 10,000 permutations per test, with the Pearson product-moment correlation coefficient as the test statistic. This is the most basic way to test whether associations across seasons deviate from random .
Second, we track the SRI for each pair of individuals as a function of time and assess, using K-means clustering, whether such temporal SRI trajectories fall into distinct types. Three temporal patterns are natural to expect. If associations are stable, we expect to see a flat SRI trajectory (type A). If associations are temporary, we expect to see an SRI trajectory with a single peak at a particular season (type B). If associations are cyclic, we expect to see an SRI trajectory with more than one peak, in corresponding seasons across years (type C). In order to minimize noise in the data due to rarely observed individuals, we limit this analysis to the 51 core individuals (see above) that can potentially form 1275 dyads. We exclude those dyads that have never associated during the study period, yielding 478 dyads with at least one non-zero SRI value within the study period. We use the correlation distance between SRI trajectories as the metric for K-means clustering because in this analysis we are interested in similarities in the shape of temporal SRI trajectories rather than in their absolute values. The number of clusters, K, that is most appropriate for the data, is chosen using the Bayesian Information Criterion (BIC) where each K-means cluster is assumed to be generated by a Gaussian distribution . We perform the K-means clustering procedure 100 times for each value of K between 2 and 15, starting with a random initial condition, and choose K at which the expected BIC is maximized as the optimal K. After determining the appropriate value of K, we run the K-means clustering algorithm 1000 times with different initial condition and pick the partition of the SRI trajectories into clusters that minimizes the sum of distances between the SRI trajectories and the cluster centroids to which they belong. In order to avoid confusion with the population-level clustering procedure below, K-means clusters are henceforth characterized by their centroids and are referred to as 'typical SRI trajectories'.
To investigate whether an individual's preferred companions change over time irrespective of the strength of the ties, we determine the top-n associates of each core individual within each season. Top-n associates of a core individual i in the given season are defined as the n core individuals (other than i) who have the highest SRI values with respect to individual i in that season. Then, in each season, a core individual i, by definition, allocates n 'companion slots' of time to spend with her top-n associates. Thus, the total number of companion slots available to each individual in 5 seasons is 5n. If individual j is present in the individual i's top-n over m seasons (1 ≤ m ≤ 5), we say that associate j 'occupies' m companion slots of individual i, or that individual i allocates m of her companion slots to individual j. We then call individual j 'an m-term associate' of individual i. If m = 1, individual j is a short-term associate of individual i, while if m = 5, individual j is a long-term associate of individual i, with obvious gradations in between. We arbitrarily set n = 5. To illustrate these concepts, consider two individuals A and B. If B's SRI index with respect to individual A is ranked 3rd, 2nd, 6th, 5th, and 9th in seasons 1 through 5 respectively, then B is in A's top-5 associates for 3 out of 5 seasons and therefore occupies 3 companion slots. Thus, individual B is a 3-term associate of individual A.
Now, if k i (m) is the number of individual i's m-term associates, then f i (m) = k i (m)m/5n is the fraction of companion slots that individual i allocates to all of her m-term associates (note that because each individual has a total of exactly 5n companion slots). Then the average fraction of companion slots that individuals allocate to their m-term companions is , where N = 51 is the total number of core individuals. indicates the extent to which individuals prefer to associate with the same companions over many seasons or to change companions from season to season.
To quantify individual variation in social behavior, we count, for each of the 51 core individuals, the number of associates that are present in her top-n for the whole observation period for n = 5. This number ranges from 0 to n = 5.
We represent the SRI matrices as weighted social-network graphs [37, 38], where nodes represent individuals, edges connect those individuals who were associated within a season, and edge weights correspond to the SRI values. An ego-network is a social network that consists only of the subject, called 'ego,' and the nodes to which she is directly connected. For each of the 88 residents, we compute five ego-network measures with intuitive biological interpretations: Size, Ties, Pairs, Density, and 2-Step Reach (defined in Table 5). We then compute the corresponding ego-network statistics, i.e. the average values of ego-network measures over all residents. As ego-network measures for different individuals are non-independent, we test for differences in ego-network statistics between seasons using a bootstrap procedure in which association data for each season are re-sampled 1000 times with replacement .
In order to investigate how the social network structure of the whole population changes over time, we construct 'network structure curves' (Figure 7) for each season using the following procedure. First, using the full SRI matrix for each season, we construct the social network graph for the whole population (i.e., we include non-resident individuals in this analysis). We then identify social clusters within the network using the Girvan-Newman algorithm for community detection . This algorithm recursively fragments the network into subnetworks by successively removing the edges with the highest between-ness . For each resulting subdivision of the network, the so-called 'modularity quotient,' Q, is computed. The modularity quotient takes values between 0 and 1, where 0 implies that the number of ties within a cluster does not exceed a random expectation, and values above 0.3 indicate potentially meaningful subdivisions [71–74]. We label the maximum modularity for a given social network as Qmax and take the partition or partitions that yield this maximum value to be the most appropriate way of subdividing the network. The term 'cluster' refers then to the set of individuals belonging to the same subnetwork in an optimal subdivision.
The original Girvan-Newman algorithm was designed to find community structure in unweighted networks. To accounts for edge weights, we add one more step. We apply the Girvan-Newman procedure to a social network that consists only of ties with strength above a particular threshold, and record the number of clusters in the resulting network (Figure 7). If multiple subdivisions of the network at a given threshold yield identical Qmax values, we record the average number or clusters (e.g. if either 14 or 15 clusters have an equivalent Qmax at a threshold of 0.1, we record 14.5 clusters). We perform this procedure iteratively for different SRI thresholds, incrementing by 0.02: 0, 0.02, 0.04, ..., 1. We then plot the average number of clusters against the SRI threshold. We call the resulting plot the 'network structure curve,' as it shows the structure of the network at a glance, being a lower-dimensional visual representation than social network graphs (Figure 7).
In order to test whether network structure curves have significantly different shapes in different seasons we modify the method used by Wittemyer et al. . For each SRI threshold value, we compare the distributions of increments of the network structure curve within a window w before and after this value, using the Mann-Whitney test. We then find points in the curve where these distributions are different from each other at the significance level 0.05. The points with significant P-values are consistent when we vary the size of window w between 0.1 and 0.3 (see Additional file 8: Figure S7 and Additional file 9: Figure S8).
Implementation and ethical statement
Permutations of associations, observed group sizes, and corresponding significance tests within and across time partitions were implemented in the OCaml programming language; code is available upon request. Network visualizations were generated in NetDraw using a graph-theoretic layout with node repulsion . All other statistical procedures and analyses were performed on Matlab v. 7.0, and R v. 2.7. This work was carried out in compliance with requirements of the Institutional Animal Care and Use Committee of the University of Pennsylvania, protocol number 801295.
SdS wishes to thank Drs. Dorothy Cheney, Robert Seyfarth, Marc Schmidt, David White, and Arthur Dunham for comments on the manuscript and research design. The authors thank two anonymous referees for their constructive comments on the manuscript, Dr. Devaka Weerakoon for advice in conducting this study, and the Department of Wildlife Conservation, Sri Lanka for granting permission to work in Uda Walawe National Park. This research was partially supported by funding from the Binns-Williams Fund in the University of Pennsylvania graduate group in Ecology and Evolution, a Doctoral Dissertation Completion grant from the University of Pennsylvania, an Integrative Graduate Education and Research Traineeship grant from the National Science Foundation (NSF-IGERT 0504487), and an Animal Behavior Society student research award to SdS as well as a grant from the US Fish & Wildlife Wildlife Asian Elephant Conservation Fund (grant no. 98210-7-G167) to SdS and ADGR.
- Emlen ST, Oring LW: Ecology, sexual selection, and the evolution of mating systems. Science. 1977, 197: 215-223. 10.1126/science.327542.View ArticlePubMedGoogle Scholar
- Geist V: Relationship of social evolution and ecology in ungulates. Am Zool. 1974, 14: 205-220.View ArticleGoogle Scholar
- Jarman PJ: The social organization of antelope in relation to their ecology. Behaviour. 1974, 48: 215-267. 10.1163/156853974X00345.View ArticleGoogle Scholar
- Wrangham RW: An ecological model of female-bonded primate groups. Behaviour. 1980, 75: 262-300. 10.1163/156853980X00447.View ArticleGoogle Scholar
- Baird RW, Whitehead H: Social organization of mammal-eating killer whales: group stability and dispersal patterns. Can J Zool. 2000, 78: 2096-2105. 10.1139/z00-155.View ArticleGoogle Scholar
- Wittemyer G, Douglas-Hamilton I, Getz WM: The socioecology of elephants: analysis of the processes creating multitiered social structures. Anim Behav. 2005, 69: 1357-1371. 10.1016/j.anbehav.2004.08.018.View ArticleGoogle Scholar
- Silk JB, Alberts SC, Altmann J: Social relationships among adult female baboons (Papio cynocephalus) I. Variation in the strength of social bonds. Behav Ecol Sociobiol. 2006, 61: 183-195. 10.1007/s00265-006-0249-2.View ArticleGoogle Scholar
- Silk JB, Alberts SC, Altmann J: Social relationships among adult female baboons (Papio cynocephalus) II. Variation in the quality and stability of social bonds. Behav Ecol Sociobiol. 2006, 61: 197-204. 10.1007/s00265-006-0250-9.View ArticleGoogle Scholar
- Parsons KM, Balcomb KC, Ford JKB, Durban JW: The social dynamics of southern resident killer whales and conservation implications for this endangered population. Anim Behav. 2009, 77: 963-971. 10.1016/j.anbehav.2009.01.018.View ArticleGoogle Scholar
- Cross PC, Lloyd-Smith JO, Getz WM: Disentangling association patterns in fission-fusion societies using African buffalo as an example. Anim Behav. 2005, 69: 499-506. 10.1016/j.anbehav.2004.08.006.View ArticleGoogle Scholar
- Ramos-Fernández G: Vocal communication in a fission-fusion society: Do spider monkeys stay in touch with close associates?. Int J Primatol. 2005, 26: 1077-1092. 10.1007/s10764-005-6459-z.View ArticleGoogle Scholar
- Sundaresan SR, Fischoff IR, Dushoff J, Rubenstein DI: Network metrics reveal differences in social organization between two fission-fusion species, Grevy's zebra and onager. Oecologia. 2007, 151: 140-149. 10.1007/s00442-006-0553-6.View ArticlePubMedGoogle Scholar
- Archie EA, Moss C, Alberts S: The ties that bind: genetic relatedness predicts the fission and fusion of social groups in wild African elephants. P Roy Soc B Biol. 2006, 273: 513-522. 10.1098/rspb.2005.3361.View ArticleGoogle Scholar
- Ramos-Fernández G, Boyer D: A complex social structure with fission-fusion properties can emerge from a simple foraging model. Behav Ecol Sociobiol. 2006, 60: 536-549. 10.1007/s00265-006-0197-x.View ArticleGoogle Scholar
- Connor R, Whitehead H: Alliances II. Rates of encounter during resource utilization: a general model of intrasexual alliance formation in fission-fusion societies. Anim Behav. 2005, 69: 127-132. 10.1016/j.anbehav.2004.02.022.View ArticleGoogle Scholar
- Hart BL, Hart LA, McCoy M, Sarath CR: Cognitive behaviour in Asian elephants: use and modification of branches for fly switching. Anim Behav. 2001, 62: 839-847. 10.1006/anbe.2001.1815.View ArticleGoogle Scholar
- McComb K, Baker L, Durant SM, Moss C, Sayialel S: Matriarchs as repositories of social knowledge in African elephants. Science. 2001, 292: 491-10.1126/science.1057895.View ArticlePubMedGoogle Scholar
- McComb K, Moss C, Sayialel S, Baker L: Unusually extensive networks of vocal recognition in African elephants. Anim Behav. 2000, 59: 1103-1109. 10.1006/anbe.2000.1406.View ArticlePubMedGoogle Scholar
- Plotnik JM, de Waal FBM, Reiss D: Self-recognition in an Asian elephant. P Natl Acad Sci USA. 2006, 103: 17053-17057. 10.1073/pnas.0608062103.View ArticleGoogle Scholar
- Shoshani J, Tassy P: The Proboscidea: evolution and palaeoecology of elephants and their relatives. 1996, Oxford; New York: Oxford University PressGoogle Scholar
- Fleischer RC, Perry EA, Muralidharan K, Stevens EE, Wemmer CM: Phylogeography of the Asian elephant (Elephas maximus) based on mitochondrial DNA. Evolution. 2001, 55: 1882-1892.View ArticlePubMedGoogle Scholar
- Roca AL, Georgiadis N, O'Brien SJ: Cyto-nuclear dissociation and the African elephant species question. Quatern Int. 2007, 169-170: 4-16.View ArticleGoogle Scholar
- Rohland NA, Malaspinas JL, Pollack M, Slatkin P, Matheus P, Hofreiter M: Proboscidean mitogenomics: chronology and mode of elephant evolution using mastodon as outgroup. Plos Biol. 2007, 5:Google Scholar
- Douglas-Hamilton I: On the ecology and behaviour of the African elephant. 1972, University of OxfordGoogle Scholar
- Moss CJ, Poole JH: Relationships and social structure of African elephants. Primate Social Relationships: An Integrated Approach. Edited by: Hinde RA. 1983, Oxford: Blackwell Scientific Publications, 315-325.Google Scholar
- Wittemyer G, Okello JBA, Rasmussen HB, Arctander P, Nyakaana S, Douglas-Hamilton I, Siegismund HR: Where sociality and relatedness diverge: the genetic basis for hierarchical social organization in African elephants. P Roy Soc B Biol. 2009, 276: 3513-3521. 10.1098/rspb.2009.0941.View ArticleGoogle Scholar
- Turkalo A, Fay JM: Forest elephant behavior and ecology: observations from the Dzanga saline. African Rain Forest Ecology and Conservation. Edited by: Weber W, White LJT, Vedder A, Naughton-Treves L. 2001, New Haven: Yale University PressGoogle Scholar
- White LJT, Tutin CEG, Fernandez M: Group composition and diet of forest elephants, Loxodonta africana cyclotis Matschie 1900, in the Lopé Reserve, Gabon. Afr J Ecol. 1993, 31: 181-199. 10.1111/j.1365-2028.1993.tb00532.x.View ArticleGoogle Scholar
- Rasmussen LEL: Chemical communication: An integral part of functional Asian elephant (Elephas maximus) society. Ecoscience. 1998, 5:Google Scholar
- Fernando P, Lande R: Molecular genetic and behavioral analysis of social organization in the Asian elephant (Elephas maximus). Behav Ecol Sociobiol. 2000, 48: 84-91. 10.1007/s002650000218.View ArticleGoogle Scholar
- Vidya TNC, Sukumar R: Social organization of the Asian elephant (Elephas maximus) in southern India inferred from microsatellite DNA. J Ethol. 2005, 23: 205-210. 10.1007/s10164-005-0144-8.View ArticleGoogle Scholar
- Eisenberg JF, McKay GM, Seidensticker J: Asian Elephants. 1990, Washington, DC: Friends of the National Zoo and National Zoological ParkGoogle Scholar
- McKay GM: The ecology and behaviour of the Ceylon elephant in south-eastern Ceylon. Asian Elephants. Edited by: Eisenberg JF, MGM Seidensticker J. 1973, Washington, DC: Friends of the National Zoo and National Zoological Park (Smithsonian Institute)Google Scholar
- Wittemyer G, Daballen DK, Rasmussen HB, Kahindi O, Douglas-Hamilton I: Demographic Status of elephants in the Samburu and Buffalo Springs National Reserves, Kenya. Afr J of Ecol. 2005, 43: 44-47. 10.1111/j.1365-2028.2004.00543.x.View ArticleGoogle Scholar
- Cairns SJ, Schwager SJ: A comparison of association indexes. Anim Behav. 1987, 35: 1454-1469. 10.1016/S0003-3472(87)80018-0.View ArticleGoogle Scholar
- Gilby IC, Wrangham RW: Association patterns among wild chimpanzees (Pan troglodytes schweinfurthii) reflect sex differences in cooperation. Behav Ecol Sociobiol. 2008, 62: 1831-1842. 10.1007/s00265-008-0612-6.View ArticleGoogle Scholar
- Wasserman S, Faust K: Social network analysis: methods and applications. 1994, Cambridge; New York: Cambridge University PressView ArticleGoogle Scholar
- Borgatti SP, Everett MG, Freeman LC: Ucinet for Windows: Software for Social Network Analysis. Ucinet for Windows: Software for Social Network Analysis. 2002, Analytic TechnologiesGoogle Scholar
- Lusseau D, Whitehead H, Gero S: Incorporating uncertainty in to the study of animal social networks. Anim Behav. 2008, 75: 1809-1815. 10.1016/j.anbehav.2007.10.029.View ArticleGoogle Scholar
- Rasmussen LEL, Krishnamurthy V: How Chemical Signals Integrate Asian Elephant Society: The Known and the Unknown. Zoo Biol. 2000, 19: 405-423. 10.1002/1098-2361(2000)19:5<405::AID-ZOO10>3.0.CO;2-J.View ArticleGoogle Scholar
- de Silva S: Acoustic communication in the Asian elephant, Elephas maximus maximus. Behaviour. 2010, 147: 825-852. 10.1163/000579510X495762.View ArticleGoogle Scholar
- Payne KB, Langbauer WR, Thomas EM: Infrasonic Calls of the Asian Elephant (Elephas Maximus). Behav Ecol Sociobiol. 1986, 18: 297-301. 10.1007/BF00300007.View ArticleGoogle Scholar
- Wakefield ML: Grouping patterns and competition among female Pan troglodytes schweinfurthii at Ngogo, Kibale National Park, Uganda. Int J Primatol. 2008, 29: 907-929. 10.1007/s10764-008-9280-7.View ArticleGoogle Scholar
- Goodall J: The Chimpanzees of Gombe: Patterns of behavior. 1986, Cambridge, MA: Harvard University PressGoogle Scholar
- Thouless CR: Home ranges and social organization of female elephants in northern Kenya. Afr J Ecol. 1996, 34: 284-297. 10.1111/j.1365-2028.1996.tb00623.x.View ArticleGoogle Scholar
- Charif RA, Ramey RR, Langbauer WR, Payne K, Martin RB, Brown LM: Spatial relationships and matrilineal kinship in African savanna elephant (Loxodonta africana) clans. Behav Ecol Sociobiol. 2005, 57: 327-338. 10.1007/s00265-004-0867-5.View ArticleGoogle Scholar
- Archie EA, Morrison TA, Foley CAH, Moss C, Alberts S: Dominance rank relationships among wild female African elephants, Loxodonta africana. Anim Behav. 2006, 71: 117-127. 10.1016/j.anbehav.2005.03.023.View ArticleGoogle Scholar
- Wittemyer G, Getz WM: Hierarchical dominance structure and social organization in African elephants, Loxodonta africana. Anim Behav. 2007, 73: 671-681. 10.1016/j.anbehav.2006.10.008.View ArticleGoogle Scholar
- de Silva S, Rangeewa ADG, Weerakoon D: Demography of Asian elephants from identified individuals at Uda Walawe National Park, Sri Lanka. Biol Conserv. 2011, 144: 1742-1752. 10.1016/j.biocon.2011.03.011.View ArticleGoogle Scholar
- Boinski S, Treves A, Chapman CA: A critical evaluation of the influence of predators on primates: effects on group travel. On The Move. Edited by: Garber SBPA. 2000, Chicago: University of Chicago Press, 43-72.Google Scholar
- Molvar EM, Bowyer RT: Costs and Benefits of Group Living in a Recently Social Ungulate - the Alaskan Moose. J Mammal. 1994, 75: 621-630. 10.2307/1382509.View ArticleGoogle Scholar
- Caro TM, Graham CM, Stoner CJ, Vargas JK: Adaptive significance of anti-predator behavior in artiodactyls. Anim Behav. 2004, 67: 205-228. 10.1016/j.anbehav.2002.12.007.View ArticleGoogle Scholar
- Loveridge AJ, Hunt JE, Murindagomo F, Macdonald DW: Influence of drought on predation of elephant (Loxodonta africana) calves by lions (Panthera leo) in an African wooded savannah. J Zool. 2006, 270: 523-530. 10.1111/j.1469-7998.2006.00181.x.View ArticleGoogle Scholar
- Power RJ, Compion RXS: Lion predation on elephants in the Savuti, Chobe National Park, Botswana. Afr Zool. 2009, 44: 36-44. 10.3377/004.044.0104.View ArticleGoogle Scholar
- McComb K, Shannon G, Durant SM, Sayialel S, Slotow R, Poole J, Moss C: Leadership in elephants: the adaptive value of age. P Roy Soc B Biol. 2011, 1-7.Google Scholar
- Katugaha HIE, de Silva M, Santiapillai C: A long-term study on the dynamics of the elephant (Elephas maximus) population in Ruhuna National Park, Sri Lanka. Biol Conserv. 1999, 89: 51-59. 10.1016/S0006-3207(98)00122-0.View ArticleGoogle Scholar
- Baskaran N, Balasubramanian M, Swamiathan S, Desai A: Home range of elephants in the Nilgiri Biosphere Reserve, South India. A week with elephant: Proceedings of the International Seminar on Asian Elephants. Edited by: Daniel JS, Datye H. 1993, Bombay: Natural History Society, Oxford University Press, 298-313.Google Scholar
- Zubair L, Siriwardhana M, Chandimala J, Yahiya Z: Predictability of Sri Lankan rainfall based on ENSO. Int J Climatol. 2008, 28: 91-101.View ArticleGoogle Scholar
- Puvaneswaran KM, Smithson PA: An Objective Classification of Homogenous Rainfall Regimes in Sri Lanka. Theor Appl Climatol. 1993, 48: 133-145. 10.1007/BF00864920.View ArticleGoogle Scholar
- Manamendra-Arachchi K, Pethiyagoda R, Dissanayake R, Meegaskumbura M: A second extinct big cat from the Late Quaternary of Sri Lanka. Raffles B Zool. 2005, 423-434.Google Scholar
- Moss C: Getting to know a population. Studying Elephants. Edited by: Kangwana K. 1996, Nairobi: African Wildlife Foundation, 58-74. AWF Technical Handbook SeriesGoogle Scholar
- Ginsberg JR, Young TP: Measuring association between individuals or groups in behavioural studies. Anim Behav. 1992, 44: 377-379.View ArticleGoogle Scholar
- Whitehead H, Bejder L, Ottensmeyer A: Testing association patterns: issues arising and extensions. Anim Behav. 2005, 69: e1-e6. 10.1016/j.anbehav.2004.11.004.View ArticleGoogle Scholar
- Sundaresan SR, Fischhoff IR, Dushoff J: Avoiding spurious findings of nonrandom social structure in association data. Anim Behav. 2009, 7: 1-5.Google Scholar
- Bejder L, Fletcher D, Bräger S: A method for testing association patterns of social animals. Anim Behav. 1998, 56: 719-725. 10.1006/anbe.1998.0802.View ArticlePubMedGoogle Scholar
- Mantel N: The detection of disease clustering and a generalized regression approach. Cancer Research. 1967, 27: 209-220.PubMedGoogle Scholar
- Sueur C, Jacobs A, Amblard F, Petit O, King AJ: How can social network analysis improve the study of primate behavior?. Am J Primatol. 2011, 703-719. 73Google Scholar
- Sokal RR, Rohlf FJ: Biometry: the principles and practice of statistics in biological research. 1995, New York: W.H. Freeman, 3Google Scholar
- Pelleg D, Moore A: X means: Extending K-means with Efficient Estimation of the Number of Clusters. Seventeenth International Conference on Machine Learning. 2000Google Scholar
- Girvan M, Newman MEJ: Community structure in social and biological networks. P Natl Acad Sci USA. 2002, 99: 7821-7826. 10.1073/pnas.122653799.View ArticleGoogle Scholar
- Clauset A, Newman MEJ, Moore C: Finding community structure in very large networks. Phys Rev E. 2004, 70: 066111-View ArticleGoogle Scholar
- Newman MEJ: Analysis of weighted networks. Phys Rev E. 2004, 70: 056131-View ArticleGoogle Scholar
- Newman MEJ: Modularity and community structure in networks. P Natl Acad Sci USA. 2006, 103: 8577-8582. 10.1073/pnas.0601602103.View ArticleGoogle Scholar
- Wolf JBW, Mawdsley D, Trillmich F, James R: Social structure in a colonial mammal: unravelling hidden structural layers and their foundations by network analysis. Anim Behav. 2007, 74: 1293-1302. 10.1016/j.anbehav.2007.02.024.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.