RESUMEN
ABSTRACT Population dynamics of aphids have been studied in sole and intercropping systems. These studies have required the use of more precise analytical tools in order to better understand patterns in quantitative data. Mathematical models are among the most important tools to explain the dynamics of insect populations. This study investigated the population dynamics of aphids Aphis gossypii and Aphis craccivora over time, using mathematical models composed of a set of differential equations as a helpful analytical tool to understand the population dynamics of aphids in arrangements of cotton and cowpea. The treatments were sole cotton, sole cowpea, and three arrangements of cotton intercropped with cowpea (t1, t2 and t3). The plants were infested with two aphid species and were evaluated at 7, 14, 28, 35, 42, and 49 days after the infestations. Mathematical models were used to fit the population dynamics of two aphid species. There were good fits for aphid dynamics by mathematical model over time. The highest population peak of both species A. gossypii and A. craccivora was found in the sole crops, and the lowest population peak was found in crop system t2. These results are important for integrated management programs of aphids in cotton and cowpea.