The Central American locust, Schistocerca piceifrons piceifrons Walker, is an economically important agricultural pest in Mexico and Central America. This species is characterized by its ability to gregarize and migrate long distances. High infestations of S. p. piceifrons have occurred in northeastern Mexico since 1998, particularly in the Huasteca region (South Tamaulipas, East San Luis Potosi, North Veracruz and the State of Hidalgo). The aim of this study was to develop a model that relates the life cycle of the Central American locust with the climatic conditions of northeastern Mexico. The model is expressed in terms of a Thermal Time Clock calculated from mean monthly temperatures and photoperiod. Historic temperature data of 20 y, from 20 meteorological stations, were used to calculate average monthly (Ψmonth) and yearly (Ψyear) values for thermal time (units are °C-days). For these calculations we used the Allen sine method, a lower temperature threshold (k1 =15.3 °C) and an upper thermal threshold (k2 = 38.5 °C). An average Thermal Time Locust Development Clock (TTLDC) was calculated using the relationship (Ψmonth/Ψyear)*360, which quantifies the monthly angular contribution on the clock (Ψ = 10.6 α). Essentially, the clock is a translator between thermal time and calendar time (or vice versa) and may be used to forecast locust-stage development dates, starting from a given “Bio-fix” point. This clock can be used to program management activities of this pest. Field observations from 2001 to 2008, on population density and development stages of the Central American locust, were consistent with the development times recorded in the TTLDC.
Introduction
The Central American locust, Schistocerca piceifrons piceifrons Walker, is an economically important agricultural pest in Mexico and Central America, affecting almost all crops. Heavy infestations of this pest have occurred since early in 1998, causing severe damage to a diversity of annual and perennial plants in the Huasteca region (South Tamaulipas, East San Luis Potosi, North Veracruz and the State of Hidalgo) (DGSV 2003, Barrientos-Lozano et al. 2004). Biology, ecology and integrated management of S. p. piceifrons have been the subject of numerous studies (Barrientos-Lozano 2001a, b, c, d; Barrientos-Lozano 2005; Barrientos-Lozano et al. 2004; Poot 2005; Ávila-Valdez et al. 2005; Ávila-Valdez et al. 2006). These works have been crucial in understanding the behavior and the interaction of this insect in the ecosystem.
The species has two generations per year: the first generation occurs from May to August-September (60–80 d) and the second from October-November to May (155–180 d). Second generation, sexually immature adults go through a diapause period from December to April; diapause ends in Spring (March-April) and with the onset of the rainy season adults become sexually mature and begin mating (Barrientos-Lozano 2001b, c).
As poikilotherms, all insects adopt the temperature of the environment and therefore their rates of development are directly related to the air temperature within a particular interval (k1, k2). In contrast to homeotherms, whose rate of development is conveniently expressed as a function of Julian or calendar time (day-1), in poikilotherms, this is not so, unless the temperature remains constant during the whole developmental cycle. Instead, the rate of development of poikilotherms is best expressed as a function of thermal time, whose units are degree-days [°C-days]-1 (Higley et al. 1986). For insect species which enter a diapause period, as is the case of the Central American locust, besides air temperature the day length may also play an important role in regulatory mechanisms. An insect's phenology models, expressed as a function of thermal time, may be useful in predicting pest development dates, monitoring scheduling, risk assessment analysis, predicting migration tracks and planning control actions (Mellors & Bassow 1983, Kauffman et al. 1985, Pinto et al. 2002). Online data bases are available to support insect model designs, for example, temperature thresholds and thermal times for insect development of many crop pests economically important: e.g., (UC-IPM Online): (1) www.ipm.ucdavis.edu/_WEATHER/ddconcepts.html), (2) www.nappfast.org/pdf/IDD%202007.pdf.
Temperature and light are the most important single nonbiotic factors associated with the rate of growth and development of the Central American locust; however, other climatic variables also play significant roles. For example, early rain promotes copulation as well as egg laying, whereas a high relative humidity protects eggs and immature stages from dehydrating, and by so doing improves locust fertility and survival rates. High relative humidity may also increase locust susceptibility to the attack of naturally occurring entomopathogenic fungi, which might result in a decrease in locust population growth. Locust food availability (plant growth) is positively related with rain frequency and total precipitation. Finally, temperature gradients, solar radiation and light intensity may also play a role in migration, inasmuch as locusts show positive thermoand phototropic effects (Porter et al. 1991; Cornford 1991 cited by Retana 2003; Barrientos-Lozano 2001a, b). While studying the effects of temperature on locust swarm migration, Uvarov (1935) reported two types of responses: a) an increase in swarm excitability when in contact with a warm soil surface and b) an increase in reflex actions, as expressed in terms of an increased jump length of individuals belonging to young swarms. More recent work (Retana 2003), establishes the relationship between sea surface temperatures “El Niño Southern Oscillation (ENSO), and the potential of this pest.
The aim of this study was to design an easy-to-use model, a Thermal Time Locust Development Clock (TTLDC), of help to biologists, entomologists and agronomists making decisions in locust management, such as predicting pest development dates, monitoring scheduling, making risk assessment analyses, predicting migration tracks and planning control actions.
Table 1.
Meteorological stations located in south Tamaulipas, México. Comisión Nacional del Agua (CNA-2006).
Materials and methods
The work was conducted in southern Tamaulipas, México, in an area between parallels 22°00′00″ and 23°15′00″ north latitude and meridians 98°00′00″ and 100°00′00 west longitude, at elevations from 0 to 380m a.s.l., from the coastal plain of the Gulf of Mexico, to the foothills of the Sierra Madre Oriental (SMO). Climatic conditions are defined by latitude, elevation and proximity to the Gulf of Mexico. The Tropic of Cancer divides the state into two climatic zones — the south with warm, humid climate and the north, less warm, with its little rain distributed throughout the year (INEGI 1983).
The grand mean monthly thermal time of the study area was calculated from the mean monthly minimum and mean monthly maximum temperatures of each of 20 meteorological stations (Table 1). Thermal time (°C-days) was calculated by the Allen sine method (1976), locating the axis of coordinates at the point (0,Tmed). Thermal thresholds for the Central American locust, k1 and k2, were taken from Barrientos-Lozano et al. (2004). Comparing temperature thermograms Tmax and Tmin with k1 and k2 there were six possible cases, but only three of them occurred in the area. This meant it required only three different equations for calculating thermal time: case1, all thermogram temperatures are within the interval k1≤T≤ k2 (Fig. 1a); case 2, the daily minimum temperature is less than the lower threshold (Tmin<k1) and the maximum is located within the threshold values (k1≤Tmax≤k2) (Fig. 1b); this condition occurs frequently during the winter months, when locusts are in diapause; case 3, the thermogram maximum daily temperature is greater than the upper threshold (Tmax>k2) and the minimum is located within the interval [k1≤Tmin≤k2] (Fig. 1c); this usually occurs during spring and summer months (April-August). Cases 4, 5 and 6 are theoretically possible, but not probable in the area under study: [case 4 (Tmin< k1) and (Tmax> K2), case 5(Tmin≥ K2), case 6 ((Tmax≤ K1),].
Day length (photoperiod), was obtained from latitude values given in Tables 1 and 2, which present monthly average day-length values (hours) for different latitudes (Mexican National Meteorological Service). For latitude values not listed in Table 2, linear extrapolation was used.
To calculate the daily number of hours at which the air temperature is less than 15.3 °C, diapause cool hours (case 2, Fig. 1b), the thermogram was adjusted by the method of Rodríguez-Absi et al. (2007) and Rodríguez-Absi (1997), using two functions: a sine function for the photoperiod and a decreasing exponential function for the night portion, as shown in Figure 2. It is convenient to separately analyze each part of the thermogram for establishing the equations which are needed in each case (Fig. 3).
Microsoft® Office Excel 2003 (Microsoft Corp., 2003) was used for: a) calculating daily thermal time using the derived equations for the three cases (1a,b,c) already discussed; b) calculating the day length (photoperiod) by linear extrapolation; c) calculating the number of hours below 15.3 °C which occur during diapause, using the derived cooling hours equations of Figure 3.
Input data included mean monthly minimum and mean monthly maximum temperatures of the 20 meteorological stations listed in Table 1, latitude values and temperature threshold values (k1 and k2). The program also calculated the monthly angular contribution to the clock and also the relationship between thermal time and a clock's angle. Thus, the following equations were also incorporated into the Excel program:
a) average monthly thermal time:
b) average yearly thermal time:
c) average angular monthly contribution to the thermal clock:
d) relationship between thermal time and a clock's angle:
Once those calculations were performed and the clock was drawn, average locust phenological development data and other relevant information were incorporated into the clock, including dates at which developmental stages do occur, diapause period, etc., (Barrientos-Lozano 2001a,b,c,d; Barrientos-Lozano 2005; Barrientos-Lozano et al. 2004; Ávila-Valdez et al. 2005; Ávila-Valdez et al. 2006).
Results and Discussion
The calculated yearly thermal time (Ψyearly) and Diapause Cool Hours (DCH) for the 20 meteorological stations are shown in Table 2. Maximum (Ψyearly) thermal time variation between the two more discrepant stations (Table 2) was 29%, (i.e., 4272.3/3311.2-1 * 100), while DCH hours variation was about 25% (i.e., 2065/1655-1 * 100). However, a better estimate of variability is obtained by dividing the values by the average value of the group, e.g., when (Ψyearly) thermal time largest and smallest values are compared with the grand mean, the variability drops to 12.5 and 14.7% respectively. Analogous calculations for DCH result in a variation of 8 and 15.6% for extreme values, when such values are compared to the grand mean.
Table 2.
Calculated thermal time (°C-days) and Diapause Cool Hours (DCH) (T≤15.3°C), for 20 meteorological stations located in south Tamaulipas, Mexico
A TTLDC for south Tamaulipas, Mexico is shown in Figure 4. Note that in contrast to a regular calendar time clock in which all month angles would be approximately equal (30°), the monthly angles in a TTLDC may be quite different; for example, May's angle is 50°, whereas December's is only 16°. This means that May offers more than three times as many degree days for locust development as does December, i.e., locust development rate in May is more than three times faster than that in December.
The relationship between (Ψyearly) thermal time (°C-days) and the angular values α is Ψ = 10.6 α, as expressed in the center of the clock (Fig. 4). This means that a unit angle in the clock = 10.6 °C-days of thermal time. For example, an egg laid on 25th April, may hatch on May 15th (first instar nymph); the angle comprised between these two dates is ∼30°, resulting in Ψ = 10.6*30 = 318°C-days, as shown in Fig. 4. The thermal time involved in transition from one stage to another is calculated in an analogous manner, i.e., by multiplying the constant 10.6 by the angle which comprises the change of stage.
The period from November to April is particularly interesting. This is when second generation adult locusts are in diapause. During this period not only do the lowest temperatures of the year occur, but also the photoperiod shortens. The average monthly number of cold days and the average monthly photoperiod for October to April are shown in Fig. 5. It is suggested that the input signal to initiate diapause is controlled by day length (Fp), particularly when Fp≤11h: this occurs November 18th to 30th; the number of cooling hours is increasing at an approximate rate of 288 hours per month as shown in Figure 5.
The diapause stage continues whenever the rate of cooling is greater than 288 hours month-1 (12 to 13 cold-days/month; T≤15.3°C/month) and concludes when Fp12h. The average total number of cooling hours, which occurs during the whole diapause interval, is around 1914 (Table 2).
Immature individuals of the Central American locust, similar to other locust migratory species, show phenotypic plasticity; this is defined as the capacity of a given genotype to produce different phenotypes in response to varying environmental conditions (Pigliucci 2001). For example, diapause entrance of S. p. piceifrons occurs with last instar, green-color nymphs. These individuals are solitary and avoid encounters with other individuals. As the number of solitary nymphs increases, they develop pink to red color with intense black pigments and become gregarious (Fig. 6).
Changes occur not only in color but also in behavior; a density-dependent gregarious mechanism is activated and may result in swarms of millions of individuals. The dorsal black stripes allow the insect to be more efficient at absorbing solar radiation, raising its body temperature and metabolic activity; this is important during the winter (diapause stage) because temperature may drop several degrees below k1 (15.3°C).
To incorporate Locust Development Stages, Diapause Cooling Hours and Day length data into the TTLDC, three circles were added (Fig. 7); the outermost circle comprises life history stages of the Central American locust. This circle moves in phase with the 10°-angle graduated circle. The inner two circles correspond to photoperiod (day length) and cold hours, respectively. For example, when the zero of the angular section (10° degrees graduated circle) is set at the date when last instar nymphs are initially observed (i.e., Nov 15th) (Bio-fix point), the average dates at which other stages of development occur, may be read on the scale dates: since the angle remains constant (it is proportional to thermal time), variation is only in dates. Data on life history, ecology and field observations from 2001 to 2008 on the Central American locust (Barrientos-Lozano 2001 a, b, c; Barrientos-Lozano 2005; Barrientos-Lozano et al. 2004; Avila-Valdez et al. 2005; Avila-Valdez et al. 2006) are consistent with the development time estimated in the average TTLDC. Variation is due to fluctuations of temperature and moisture typical of different years, so the life cycle of this pest is advanced or delayed.
Conclusions
This average Thermal Time Locust Development Clock (TTLDC), obtained for the Central American locust, can be used to predict timing of phenological stages of this pest; this implies a starting date (Bio-fix). This date is experimentally established after obtaining quality data on population dynamics and on the life cycle of the pest. Reliability of predictions depends upon accuracy of the starting date and deviations in weather condition from the average year. When the TTLDC was compared with field observations, the differences were relatively small.
Building a TTLDC involved calculations using average values: temperature, day length, diapause cool hours and calendar dates. Thus it may be expected that the clock provides reasonably good estimates for predicting locust phenological development in south Tamaulipas, México. However, more accurate clocks may be built, incorporating other variables such as relative humidity, rain etc. More specific clocks could be developed for smaller areas and used in local pest management programs. At present, several minimeteorological stations connected to the web have been installed throughout. This may facilitate analysis on pests' life cycle vs climate data using the methods already discussed.
Acknowledgments
Special thanks to the Mexican Council for Science and Technology (Fondo Mixto de Fomento a la Investigación Científica y Tecnológica, CONACYT-Gobierno del Estado de Tamaulipas), for funding the research Project “Desarrollo de un proceso tecnológico para el control biológico de la langosta centroamericana (Schistocerca piceifrons piceifrons Walker) y acridoideos plaga en la region sur de Tamaulipas”. Clave Tamps-2005-C08-26. We acknowledge the facilities provided by the Instituto Tecnológico de Cd. Victoria, and assistance and support during fieldwork provided by Pablo García-Salazar, Junta Local de Sanidad Vegetal, Mante-Comité Estatal de Sanidad Vegetal Tamaulipas. The assistance provided by Jonathan Michel Álvarez-García for editing figures and tables is greatly appreciated.