INTRODUCTION
In tropical ecosystems, termites are among the most important soil invertebrate groups for playing an essential role in decomposing processes, nutrient cycling and changes in soil properties (Lee and Wood, 1971; Holt and Lepage, 2000; Korb and Linsenmair, 2001). Many other organisms depend on termites as a source of food and organic matter, or use their mounds for shelter (Constantino and Schlemmermeyer, 2000). In the Brazilian Cerrado, termite mounds are abundant, and some types have a high population density (Lima-Ribeiro et al., 2006). In cultivated pastures of this region, termite populations can be very high as indicated by large densities of mounds, but the diversity of species is greatly reduced.
Species that feed on soil can either build their mounds on the surface or underground (Oliveira et al., 2011). These mounds have different shapes and characteristics, depending on the termite species. Constantino (2005) classified these mounds as epigeous, subterranean, arboreal and within a single piece of wood. According to Bandeira and Martius (2009), one mound usually represents one colony. There are, however, species whose colonies may have two or more sub-mounds spread over large areas, that is, polycalic or compound mounds (Martius, 2001).
Termite colonies can contain thousands and even millions of individuals (Bandeira and Vasconcellos, 2004). A mature colony produces winged termites that can establish new colonies. Termite flights usually occur once a year, when mature nymphs develop wings and transform themselves into adult imagoes or ‘alates’ (Bandeira and Martius, 2009). This phenomenon is primarily seasonal and associated with climatic variations, especially heat and relative humidity (Rebello and Martius 1994; Medeiros et al., 1999; Martius, 2003). Flights of the same species in the same area are synchronized (Cancello and Schlemmermeyer, 1999).
The spatial distribution of termites mounds is related to environmental characteristics such as altitude, temperature, rainfall, local availability of food resources and nesting; as well as intra- and inter-specific relationships and degree of environmental disturbance (Eggleton et al., 2002; Florencio and Diehl, 2006). Species that build epigeous mounds, such as Cornitermes cumulans, can reach high population densities, especially in open environments (Constantino, 2005), such as pastures.
According to Cunha and Morais (2010), the implantation of pastures often increases mound density because of the environmental homogeneity and absence of competitors. Despite the beneficial importance of termites to soil processes the useful pasture area can be reduced by very dense populations of epigeous mounds (Valério, 1995), which has led cattle farmers to find ways to destroy the mounds (Ávila and Goulart, 1995). Termite mounds have been commonly associated with pasture degradation, but this hypothesis has been recently refuted by Lima et al. (2011) in an evaluation of soil and plant parameters in Cerrado areas.
The spatial distribution of mounds may be an indicator of underlying mechanisms regulating population dynamics (Korb and Linsenmair, 2001), For instance, regular spacing may indicate intra-specific competition or density-dependent predation. Clumped distribution patterns may be the product of incomplete dispersion (i.e. when new termite colonies are founded by alates whose ‘nuptial flights’ end not very far from the mother colony). The spatial distribution of termites has been studied using various methods. Clark and Evans (1954) and later, Korb and Linsenmair (2001), evaluated spatial variability by a distance index using the nearest-neighbor technique. Geostatistical techniques were applied by Crist (1998) to assess the spatial distribution of subterranean termites, and more recently, Grohmann et al. (2010) analyzed mound distribution using Ripley’s K function, pair correlation function and similar techniques, to analyze mound distribution patterns of a fungus-growing termite species.
In Brazil, studies on the spatial distribution of termites were conducted by Bruschini (2006), who applied the Clark and Evans (1954) distance index, and by Cunha (2011), who analyzed the spatial distribution of mounds by dispersion and clustering. Martius (2003) mapped termite distribution in a varzea floodplain forest and associated it to environmental parameters such as flooding. Understanding the processes that regulate the population dynamics of termites is particularly important since the distribution is regulated at two stages in the lifetime of a termite colony, with different factors affecting the founding alates than the established colony. Mounds can be prominent landscape features even after the colony that built them died. Therefore, counting mounds without assessing their vitality can lead to an overestimation of the termite population. Only assessing the vitality allows a true diagnosis of termite occupation; this information has been ignored in many surveys conducted in Brazilian pastures with epigeous mounds.
The semivariogram analysis can show whether the spatial distribution of vitality, height and circumference of termite mounds is random or spatially dependent. We hypothesized that if the spatial distribution is not random, some (biotic and/or abiotic) factors are controlling the dynamics of termite mounds. Therefore, the purpose of this study was to apply geostatistical methods to assess the spatial distribution of epigeous termite mounds in pasture areas and to evaluate their vitality.
MATERIAL AND METHODS
The study was carried out in a pasture area of the Fazenda Campana in Rio Brilhante, Mato Grosso do Sul (22° 16’ 29’’ S; 54º 48’ 57” W and altitude of 353 m asl). This area, covered with Cerradão vegetation, characterized by a continuous vegetation canopy and tree cover that can range from 50 to 90 %, is part of the Cerrado biome, where the climate is tropical with a dry season (Köppen-Geiger classification Aw). The mean annual temperature is 20 °C, ranging from 18 to 22 °C, and annual rainfall between 1,500 and 1,750 mm. The predominant soil type in the region is Latossolo Vermelho (Oxisol) (Embrapa, 2006).
Two cultivated pastures with Brachiaria brizantha (cv. Marundu) consisting of two paddocks with a total area of nearly 60 ha were selected for their management history, proximity, topographic similarity, age and a contrast, first assessed visually only, in terms of termite mound density. In each area, cattle grazed in a rotational system, moving from one paddock to another every 7 to 8 days. After pasture establishment 1992, the areas had been maintained without fertilizer. One area had more epigeous mounds and a lower number of trees and shrubs (Pasture 1); the other area (Pasture 2) had lower nest density and a higher number of trees and shrubs.
In each pasture, four continuous 1-ha paddocks were demarcated. A large paddock area was used to enhance the accuracy of mound distribution assessment and because it was suitable for geostatistical methods, which require sampling point insertion into a grid. The mean number of mounds ha-1 was calculated based on the nests counted per plot.
Mound evaluation and termite collection were performed in May 2010 (Pasture 1) and May 2011 (Pasture 2). One 1.0-ha plot per pasture was randomly selected to study spatial variability. Plot 1 was chosen in Pasture 1 and plot 3 in Pasture 2 (highlighted plots in figure 1). The epigeous mounds in each plot were georeferenced using a GPS (Trimble Pathfinder - ProXT) with horizontal precision of ± 0.8 m and vertical precision of ± 1.3 m. Next, the height, circumference and vitality (termite presence) of each mound were evaluated. Mound vitality was assessed in the early morning and late afternoon, by observations of the external aspect of the mound, and detection of termite presence or absence when perforating the mounds with an iron shovel. In some cases, the mounds were destroyed to certify that they had been completely abandoned.

Figure 1 Semivariogram of termite nest circumference and height evaluated in Plot 1 from Pasture 1, on the Fazenda Campana, Rio Brilhante, MS.
The area occupied by the mounds in pastures 1 and 2 was calculated from the relationship between the sum of areas occupied by the mounds in the chosen plots and the total area of the plots (1 ha), as shown in equation 1:
where: Atn (%) is the area (%) occupied by epigeous mounds in a chosen plot; and r the radius of each epigeous mound.
The radius of each mound was determined from its circumference, according to equation 2:
where: c is the circumference at the base of the mound, measured in the field.
For species identification, 10 % of the mounds were separated into their epigeous and underground portions to sample all castes of the colony. Samples were stored in flasks with 80 % alcohol for further taxonomic identification. Identification of the species and morphospecies level was performed according to Cancello (1989) and Constantino (1999), at the Zoology Museum of the University of São Paulo (MZUSP). After identification, part of the material was added to the MZUSP collection.
Descriptive analysis and geostatistics were performed using GEOESTAT software, according to Vieira et al. (1983). Descriptive analysis was performed using the STAT.EXE module of the software to check for normal distribution of the dataset and detect outliers. The AVARIO.EXE module was then used to analyze the spatial dependence of termite mounds based on mound height and circumference. The AVARIO.EXE module calculates the semivariance of “mound height” and “circumference”, by the following equation:
where γ(h): semivariance estimated at distance h; N(h): number of pairs of values [Z(x i), Z(x i + h)] separated by the h vector; xi: spatial position of the Z variable; and Z: values of the variables (mound height and circumference).
Semivariance graphs (semivariograms) were constructed and adjusted to numerical models using Microsoft Excel. Maps of the spatial variability of termite mound properties were produced using Surfer 7.02 software (Golden Software, 2000).
RESULTS AND DISCUSSION
All mounds in both pastures had been constructed by a single termite species (Cornitermes cumulans). In Pasture 1, the mean number of mounds ha-1 was 68, and 71 were selected in a plot for spatial analysis. Pasture 2 contained 127 mounds, with 115 in the selected plot (Table 1). The selected mounds occupied 0.4 and 1 % of the total area of Pastures 1 and 2, respectively.
Table 1 Number of termite mounds in the plots and mean number of mounds per pasture on the Fazenda Campana, Rio Brilhante, MS
Plot 1 | Plot 2 | Plot 3 | Plot 4 | Mean | Atn(1) | |
---|---|---|---|---|---|---|
% | ||||||
Pasture 1 | 71 | 62 | 77 | 63 | 68.3 | 0.4 |
Pasture 2 | 90 | 195 | 115 | 108 | 127.0 | 1.0 |
(1) % Atn: percent of the area occupied by epigeous mounds in the plot selected.
Altitude, mound heights and circumferences were similar in both areas (Table 2). The coefficients of variation (CV) for altitude were consistently low due to the flat topography in both areas. For mound circumference and height, CVs were high. The distribution of data frequency was asymmetric and in general close to zero, particularly in Pasture 1. The highest asymmetry value was found for altitude data in Pasture 2. However, data distribution can be assumed to be close to normal. The average mound height observed in the two areas was 65 and 59 cm in Pasture 1 and 2, respectively, and the maximum height was 130 cm. The basal circumference of the mounds ranged from 0.1 to 3.9 m in Pasture 1 and from 0.6 to 3.4 m in Pasture 2, with means of 2.58 and 2.26 m, respectively.
Table 2 Descriptive statistical analysis of data of epigeous termite mounds
Parameter | Altitude | Circumference | Height |
---|---|---|---|
m asl | m | cm | |
Pasture 1 | |||
Minimum value | 330.6 | 1.00 | 20.0 |
Maximum value | 334.6 | 3.90 | 130.0 |
Mean | 332.7 | 2.58 | 65.7 |
Asymmetry | -0.001 | -0.51 | 0.19 |
Kurtosis | -1.210 | -0.75 | -0.88 |
CV (%) | 0.330 | 28.82 | 42.60 |
Standard deviation | 1.11 | 0.74 | 27.96 |
Pasture 2 | |||
Minimum value | 323.2 | 0.60 | 10.0 |
Maximum value | 329.7 | 3.44 | 128.0 |
Mean | 325.6 | 2.26 | 59.5 |
Asymmetry | 0.70 | -0.38 | 0.21 |
Kurtosis | -0.77 | -1.03 | -0.72 |
CV (%) | 0.54 | 31.62 | 46.40 |
Standard deviation | 1.76 | 0.71 | 27.58 |
The correlation coefficients between the circumference and height of epigeous mounds were high and significant (p<0.05; Table 3), but these variables were not associated to altitude, as expected, given the low altitude variation. Mound circumference was highly correlated to height, as the highest mounds tended to have a higher base circumference.
Table 3 Matrix of correlation between sizes of epigeous termite mounds
Altitude | Circumference | Height | |
---|---|---|---|
Pasture 1 | |||
Altitude | 1.0000 | ||
Circumference | 0.2963 | 1.000 | |
Height | 0.2863 | 0.803 | 1.000 |
Pasture 2 | |||
Altitude | 1.0000 | ||
Height | 0.0659 | 1.0000 | |
Circumference | 0.0849 | 0.8163 | 1.000 |
In both pasture areas, data were recorded in Plot 1.Values without bold indicate no statistical significance (Pearson’s correlation, p<0.05).
The semivariograms of mound measurements in Pasture 1 (Figures 1 and 2), show a “pure nugget effect” for height and circumference. According to Vieira (2000), an experimental semivariogram with a pure nugget effect for a property means that its spatial distribution is random. This pattern can be visualized using experimental semivariance values (γ(h)) at short distances (5 to 15 m). In Pasture 1, height and circumference measures of the closest mounds (approximately 15 m) were more divergent than those of mounds separated by larger distances. In Pasture 2, γ(h) values were practically the same at any distance, confirming the random distribution pattern.

Figure 2 Semivariogram of termite nest circumference and height evaluated in Plot 3 from Pasture 2, on the Fazenda Campana, Rio Brilhante, MS.
In Plot 1 of Pasture 1, 91 % of the mounds were active, and 84 % were active in Plot 3 of Pasture 2. Spatial variability maps show vitality distribution (Figure 3a,c) and its relationship to mound size (Figure 3b,d). Vitality distribution was particularly associated with mound height, indicating that the largest mounds were usually active, mainly in Pasture 2, where most abandoned mounds were found.

Figure 3 Spatial distribution of the vitality (a, c) and height measurements and vitality (b, d) of epigeous termite mounds in Pasture 1 and 2, respectively, on the Fazenda Campana, Rio Brilhante, MS. (+ = abandoned mounds without termite activity; • = active mounds). The circles in the graph explain the random distribution of mounds.
The occurrence of a single termite species (C. cumulans) likely results from the absence of old epigeous mounds, which are commonly cohabited by more than one termite species. Based on the pattern of mounds produced by C. cumulans, records on recent termite activity, mound integrity and the low mound height, the mounds analyzed were assumed to be new. Cornitermes cumulans is the most common termite species found in South-Central Brazil (Cancello, 1989). In addition, Valério et al. (2006) reported that Cornitermes is the predominant species, occupying 94 % of the mounds recorded in a pasture area covered with different species of forage grass.
Our findings on the occupation of pasture areas by termite mounds are similar to those reported in previous studies conducted in the Cerrado (Czepak et al., 2003; Cunha and Morais, 2010). These authors found a prevalence of C. cumulans mounds, ranging from 0.59 to 0.65 m in height, and, as in this study, concluded that the mounds do not significantly compromise the area of useful pasture. The higher occurrence of mounds with low mound height and circumference in Pasture 2 (Figure 4c,d) explains the high coefficient of variation detected. Fadini et al. (2002) reported that termite mounds frequently reach a height of 1.5 m in pasture areas. However, Cunha (2011) found that epigeous C. cumulans mounds are on average 0.42 m high, which was is in line with our study. The positive correlation between circumference and height is related to mound development, and increases with colony size. Therefore, considering the shape of C. cumulans mounds, higher mounds have a greater circumference. In older mounds, where the upper portion deteriorates due to rainfall, height commonly decreases in relation to circumference.

Figure 4 Spatial distribution of height (a, c) and circumference (b, d) measurements of epigeous termite mounds in Pasture 1 and 2, respectively, on the Fazenda Campana, Rio Brilhante, MS. The circles in the graph explain the random distribution of mounds.
The results of spatial distribution agree with Grohmann et al. (2010) and Korb and Linsenmair (2001). These authors analyzed mound distribution patterns as a function of size and vitality, and although they used other methods and studied other species than in the present study, they also observed no spatial dependence of mound distribution in most plots. Crist (1998) also applied geostatistical techniques to study termites, but focused on individuals attracted by bait rather than epigeous mounds. As in the present study, the author also reported random distribution in termite mounds.
The random behavior detected in the semivariograms can be inferred in figures 3a to 4b, where the symbols on the graphs indicate mounds of closest proximity and differences in height and circumference. These symbols show high γ(h) values at short distances. In addition, the spatial distribution of both height and circumference at these points is similar, which is expected, given the high correlation between these properties. Small mounds are usually new constructions that grow near old, large mounds. In certain cases, they can also be ramifications of a larger mound, that is, a polycalic mound.
Other studies investigated the spatial distribution of epigeous Cornitermes spp mounds in Brazil using different procedures. According to Cunha (2011), random distribution of mounds can be explained by the absence of intraspecific competition. Using a distance index, Bruschini (2006) found that low mound density (65) and absence of intraspecific competition and food resources affected the spatial distribution of mounds. Comparing these studies, Cunha (2011) took all mounds in the plots into consideration, whereas Bruschini (2006) analyzed only the inhabited.
According to Bruschini (2006), the area size may affect mound distribution and should therefore be considered in comparative results. Cunha (2011), and Korb and Linsenmair (2001) hypothesized that intraspecific competition can explain the spatial distribution of termite mounds. In addition, according to these authors, colony growth and competition for food seems to change random distribution to a more regular pattern. Moreover, Jiménez et al. (2006) reported that unequal distribution of resources can determine the spatial distribution of soil organisms, which can alternate between high and low density in a same area.
Lima et al. (2011) recorded high productivity and consequent homogeneity of food distribution for termites in both pasture areas studied. Owing to these characteristics, the hypothesis of intraspecific competition cannot explain the random spatial distribution of termite mounds. It is important to underscore that mound distribution may be directly related to flight distance and climatic conditions at flight time, most likely affecting the location of new colonies formed by breeding pairs. Since the flight distance cannot be accurately evaluated, randomness is an acceptable pattern of termite mound distribution.
The high proportion of active mounds indicates the potential for an increase in the number of mounds in the area. Oliveira et al. (2011) found higher mound occurrence (408 mounds ha-1) on the Central Cerrado Plateau. Confirming findings of Lima et al. (2011), the pasture areas studied were considered homogeneous, productive and undegraded, when assessing indicators of pasture degradation (soil fertility and microbial biomass, light fraction of organic matter, litter and existing regrowth) in the areas of this study, differing from many pastures areas in Brazil that are densely occupied by epigeous termite mounds. Ecological studies on the spatial distribution pattern and other termite mound parameters, such as vitality, are important in explaining the higher or lower formation intensity of epigeous mounds in pasture areas. Further studies should investigate this issue, testing different conditions of pasture productivity and ages.
CONCLUSIONS
In the pasture areas, the number of mounds without epigeous termite activity was much lower than that of active or inhabited mounds.
The spatial distribution of mounds was epigeous, confirming reports in the literature and suggesting the use of geostatistics as important tool for ecological studies.
The random distribution of the mounds indicates that more research is needed to explain the dynamics of termites in the soil and their relationship with biotic and abiotic factors.