Acessibilidade / Reportar erro

Mathematical models for Isoptera (Insecta) mound growth

Modelos matemáticos para crescimento dos ninhos de Isoptera (Insecta)

Abstracts

In this research we proposed two mathematical models for Isoptera mound growth derived from the Von Bertalanffy growth curve, one appropriated for Nasutitermes coxipoensis, and a more general formulation. The mean height and the mean diameter of ten small colonies were measured each month for twelve months, from April, 1995 to April, 1996. Through these data, the monthly volumes were calculated for each of them. Then the growth in height and in volume was estimated and the models proposed.

Nasutitermes coxipoensis; Isoptera; Insecta; mound growth


Nessa pesquisa, propusemos dois modelos matemáticos derivados da curva de Von Bertalanffy para o crescimento dos ninhos de Isoptera; um específico para Nasutitermes coxipoensis; e outro mais geral. A altura média e o diâmetro médio de dez colônias pequenas foram medidas a cada mês, de abril de 1995 a abril de 1996. Através desses dados, foi calculado o volume mensal de cada colônia. Posteriormente os crescimentos em altura e volume foram estimados e os modelos propostos.

Nasutiternes coxipoensis; Isoptera; Insecta; crescimento de ninhos


ECOLOGY

Mathematical models for Isoptera (Insecta) mound growth

Modelos matemáticos para crescimento dos ninhos de Isoptera (Insecta)

Buschini, MLT.I, * * e-mail: isatunes@yahoo.com.br ; Abuabara, MAP.II; Petrere-Jr, MiguelIII

IDepartamento de Biologia, Universidade Estadual do Centro-Oeste – UNICENTRO, Rua Presidente Zacarias, 875, CEP 85010-990, Guarapuava, PR, Brazil

IIDepartamento de Matemática da Universidade Estadual Paulista – UNESP, CP 199, CEP 13506-900, Rio Claro, SP, Brazil

IIIDepartamento de Ecologia, Universidade Estadual Paulista – UNESP, CP 199, CEP 13506-900, Rio Claro, SP, Brazil

ABSTRACT

In this research we proposed two mathematical models for Isoptera mound growth derived from the Von Bertalanffy growth curve, one appropriated for Nasutitermes coxipoensis, and a more general formulation. The mean height and the mean diameter of ten small colonies were measured each month for twelve months, from April, 1995 to April, 1996. Through these data, the monthly volumes were calculated for each of them. Then the growth in height and in volume was estimated and the models proposed.

Keywords:Nasutitermes coxipoensis, Isoptera, Insecta, mound growth.

RESUMO

Nessa pesquisa, propusemos dois modelos matemáticos derivados da curva de Von Bertalanffy para o crescimento dos ninhos de Isoptera; um específico para Nasutitermes coxipoensis; e outro mais geral. A altura média e o diâmetro médio de dez colônias pequenas foram medidas a cada mês, de abril de 1995 a abril de 1996. Através desses dados, foi calculado o volume mensal de cada colônia. Posteriormente os crescimentos em altura e volume foram estimados e os modelos propostos.

Palavras-chave:Nasutiternes coxipoensis, Isoptera, Insecta, crescimento de ninhos.

1. Introduction

All termites are social insects and live in communities, large or small, within the limits of a nest-system. The nest and associated structures, such as mounds, subterranean galleries and covered runways, comprise a closed system largely isolated from the external environment but allowing for the egress of foraging parties and flight of alates (Lee and Wood, 1971). In this system, within which the micro-climate can be controlled within certain limits, food can be stored, and they obtain some protection from natural enemies (Lee and Wood, 1971; Banerjee, 1975). The nest-system of these insects is not a static structure and is enlarged as the colony grows (Noirot, 1970).

Three recognisable developmental phases exist in termite societies. The first corresponds to a juvenile phase, during which there is just differentiation of the neuter castes. Soon after, the colony passes to the mature or adult phase, period in which the winged reproductives are produced. Finally, a phase of senescence occurs, when there is a decline in the nest population (Noirot, 1969). Collins (1981) deduced that the populations in the nests of Macrotermes bellicosus Smeathman 1781 increase exponentially for 4-6 years reaching the maximum number of individuals in its neuter population at the end of this period. Resources are subsequently channelled into the production of winged reproductives and, years later, its population declines and the colony dies.

Considering all three development stages of the Isoptera colonies, recognized by Noirot (1969), one can consider the growth of the colony analogous to that of an organism. Thus, in this study we proposed two mathematical models in order to describe the growth of Nasutitermes coxipoensis Holmgren 1910 (Model 1) mounds derived from the Von Bertalanffy (1938) growth curve and a more general formulation (Model 2) for other termite species. Von Bertalanffy developed one of the first models for organic growth, particularly for fishes.

2. Material and Methods

2.1. Study area

Mounds sampled in this study were located in a Cerrado area in the Municipality of Itirapina, São Paulo, Brazil (22° 15' S and 47° 49' W, altitude of 765 m). The word "Cerrado" is a Portuguese term meaning "half-closed" or "dense". In Brazilian terminology it also describes a particular kind of vegetation similar to savannah, although having much broader physiogonomic variation in size and density of trees. That is, the word Cerrado describes a gradient of vegetation comprising "savannah grassland", "low trees and shrub savannah" and "savannah woodland". It is important to point out that the similarities to savannah are only physiogonomic and not floristic. It covers about 2 million km2, or 25% of the whole Brazilian territory (Ferri, 1976). This supports private farming enterprises where one of the main activities is cattle rearing. Therefore, at times of drought when fire becomes frequent in this locality, the impact of the cattle also becomes more apparent.

The study area is an alluvial plain covered by sandy sediment, with a deep, quartzes' sand soil type (Oliveira and Prado, 1984) and has 12,750 m2.

Mean annual rainfall is 1,425 mm, with the rainy season extending from October to March, when 84% (1,199 mm) of the precipitation occurs. The most rainy months are December, January and February, with precipitation average values equal to 288, 266 and 262 mm, respectively. The driest months are July and August, with 16 and 19 mm of precipitation, respectively. The mean annual temperature is 19.7 °C, with January and February being the hottest months with respective mean values of 22.2 and 22.3 °C. The coldest months are June and July with respective mean temperatures of 16.4 and 16.2 °C.

2.2. Sampling program

The mean height (n = 3) and the mean diameter (n = 3) of ten previously marked mounds, whose initial volumes varied from 3.88 dm3 to 19.38 dm3, were measured each month for twelve months, from April, 1995 to April, 1996, with the aid of two wood stakes approximately 1 metre in length. Considering that the biggest volume found was equal to 192.46 dm3 (Buschini, 1996), this research was been accomplished with small colonies. During the measurements the stakes were vertically placed on each side of the mound and the diameter measured at this distance. To measure the height, a single stake was placed horizontally on the apex of the mound and, with the aid of the metric tape, the distance between the stake and the soil was measured as closely as possible to the base of the mound.

Although some of the mounds present deformations and considering that they are entirely epigean, it was assumed that they have approximately a spherical calotte shape (Figure 1). Thus, the volume of the mounds was evaluated by the formula: V = π/6 .h{3r2 + h2} where r = d/2 = radius and h = height of the mound (Gieck, 1979).


3. Results and Discussion

3.1. Model 1 for Nasutitermes coxipoensis mounds

3.1.1. Growth in height

The growth in height can be expressed by the differential equation:

where:

ht = height of the mounds in the instant t; k = constant of growth; hmax = maximum height, that is, the maximum height reached by an adult mound.

Let us suppose that in the instant t = 0 ⇒ ht = 0

Then,

integrating:

When t = 0 ⇒ -ln(hmax)= c, thus,

The graph of this function is given in Figure 2. Notice that hmax is an asymptotic value that needs to be determined together with k, through a non-linear regression.


A more flexible model was proposed by Richards (1959) to describe different patterns of organic growth in which the growth rate, in length, is determined by the differential equation:

where Dmax is the maximum size reached by the organism and k is a growth constant. The solution of this equation is the growth curve defined by:

If m = 2/3 then we have Von Bertalanffy's model (Brown and Rothery, 1993).

3.1.2. Growth in volume

Considering that the nests of Nasutitermes coxipoensis possess the volume of a spherical shape with height h and radius r, we have:

A dispersion diagram of r vs. h of each N. coxipoensis mound was drawn. From this graph (Figure 3) we verified that a linear relationship exists, approximately, among them. Then, a model of linear regression r = βo + β1h was used, and we observed that the slope coefficient of the adjusted straight line is b1 = 0.569.


A test of the hypothesis {Ho1 = 0.5 and H11 diferente_symbol 0.5} was accomplished and it showed that we should accept the null hypothesis Ho, that is, there is strong evidence that β1 = 0.5. Thus, firstly, we suppose that:

Replacing (2) in (6) then:

we have:

Studying the function: Vt = Vmax (1 - e-kt)3

that is, Vt is crescent

Thus, in t = we have one inflexion point;

for t < , Vt has an upward concavity;

for t < , Vt has a downward concavity.

Now, limt→∞ Vt = limt→∞ Vmax(1 - e-kt)3 =

limt→∞ Vmax = Vmax1 = Vmax

Figure 4 illustrates the behaviour of this function, also asymptotic for Vmax.


3.2. Model 2 for other termite species

We propose below a more general formulation, without the supposition that rt = 2ht.

The graph of Vt will not be sketched here because, as it can be observed, the sign of the second derivate d2V ⁄dt2 depends on the specifications of rt and ht, for each species of termites to be studied.

Faced with these models proposed for N. coxipoensis in particular, and for other termites species, it can be thought that the juvenile phase would be equivalent to the initial stages of development of the colony with rates of relative growth {(w2-w1)/(w1)} and absolute {(w2-w1)/(t2-t1)} height. The mature phase, for its time, would correspond to the inflexion point with a decrease in the growth rates due to a larger reproductive investment (beginning of the production of winged reproductives in termites). And, finally, the older phase would be the stage during which the colonies reach an asymptotic weight (volume). However, other studies will be necessary to test and to discuss the models presented in this study.

This model can be tested simply by measuring height and diameter of the studied species nests. As mentioned above, the Vt graphs have to be sketched for each termite species with their specific values of rt and ht. It is important to investigate the relationship between these two variables by initially plotting both variables on the graph and verifying the existence of linearity between them. If rt = 2ht, the species in question will present the same growth pattern as N. coxipoensis, equivalent to that proposed by Von Bertalanffy (1938).

Acknowledgements – We would like to thank Prof. Dr. Ana Maria Costa Leonardo for her trust in us. CNPq supplied financial support for the accomplishment of this study.

Received December 18, 2006

Accepted April 26, 2007

Distributed August 31, 2008

  • BANERJEE, B., 1975. Growth of mounds and foraging territories in Odontotermes redemanni (Wasmann) (Isoptera: Termitidae). Insects Sociaux, vol. 22, no 2., p. 207-212.
  • BROWN, D. and ROTHERY, P., 1993. Models in biology: mathematics, statistics and computing. New York: John Wiley & Sons.
  • BUSCHINI, MLT., 1996. Ecologia e biologia de uma espécie de Nasutitermes (Isoptera : Termitidae) no cerrado de ItirapinaSP Rio Claro : Instituto de Biociências, UNESP. 129 p. [PhD Thesis]
  • COLLINS, NM., 1981. Populations, age structure and survivorship of colonies of Macrotermes bellicosus (Isoptera: Macrotermitinae). Journal of Animal Ecology, vol. 50, no. 1, p. 293-311.
  • FERRI, MG., 1976. Ecologia dos cerrados. In Simpósio sobre o cerrado São Paulo: Itatiaia. p. 15-36.
  • GIECK, K., 1979. Manual de fórmulas técnicas São Paulo: Hemus.
  • LEE, KE. and WOOD, TG., 1971. Termites and soils London: Academic Press.
  • NOIROT, C., 1969. Formation of castes in the higher termites. In KRISHNA, K. and WEESNER, FM. (Eds.). Biology of termites New York: Academic Press. p. 311-350
  • -, 1970. The nests of termites. In KRISHNA, K. and WEESNER, FM. (Eds.). Biology of termites New York: Academic Press. p. 73-126
  • OLIVEIRA, JB. and PRADO, H., 1984. Levantamento pedológico semidetalhado do Estado de São Paulo: quadrícula de São Carlos II. Campinas: Instituto Agronômico. Memorial Descritivo.
  • RICHARDS, FJ., 1959. A flexible growth function for empirical use. Journal of Experimental Botany vol. 10, no. 10, p. 290-300.
  • VON-BERTALANFFY, LV., 1938. A quantitative theory of organic growth (inquiries on growth laws II). Hum. Biol, vol. 10, no 2., p. 181-213.
  • *
    e-mail:
  • Publication Dates

    • Publication in this collection
      23 Sept 2008
    • Date of issue
      Aug 2008

    History

    • Received
      18 Dec 2006
    • Accepted
      26 Apr 2007
    Instituto Internacional de Ecologia R. Bento Carlos, 750, 13560-660 São Carlos SP - Brasil, Tel. e Fax: (55 16) 3362-5400 - São Carlos - SP - Brazil
    E-mail: bjb@bjb.com.br