Next Article in Journal
A Limiting Factor of Sex Attractants of Bactrocera dorsalis (Diptera: Tephritidae), Verified under Laboratory Conditions
Previous Article in Journal
A New Approach for Detecting Sublethal Effects of Neonicotinoids on Bumblebees Using Optical Sensor Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of the Effect of the ENSO Cycle on the Distribution Potential of the Genus Anastrepha of Horticultural Importance in the Neotropics and Panama

by
Arturo Batista Degracia
1,*,
Julián Ávila Jiménez
2,
Anovel Barba Alvarado
3,
Randy Atencio Valdespino
3 and
Mariano Altamiranda-Saavedra
4
1
Instituto de Innovación Agropecuaria de Panamá (IDIAP), C. Carlos Lara 157, Ciudad del Saber 0843-03081, Panama
2
Facultad de Ciencias, Universidad Pedagógica y Tecnológica de Colombia, Avenida Central del Norte 39-115, Tunja 150003, Colombia
3
Institute of Agricultural Innovation of Panama/National Research System of Senacyt-Panama, Panama City 0816-02852, Panama
4
Grupo de Investigación Bioforense, Tecnológico de Antioquia Institución Universitaria, Medellín 050005, Colombia
*
Author to whom correspondence should be addressed.
Insects 2023, 14(8), 714; https://doi.org/10.3390/insects14080714
Submission received: 14 July 2023 / Revised: 2 August 2023 / Accepted: 15 August 2023 / Published: 18 August 2023
(This article belongs to the Section Insect Pest and Vector Management)

Abstract

:

Simple Summary

The genus Anastrepha (Diptera: Tephritidae) is a group of fruit flies that contain significant insects (pests) of agricultural importance since they alter the quality and reduce the yield of fruit and vegetable production. This in turn projects great uncertainty regarding safety due to quarantine restrictions around the world. In this work, potential distribution models are presented for four selected species of this genus of fruit flies in different scenarios of climate variability; these models can constitute an input for the design of preventive control plans and early warning systems for these pest insects. The models were designed using the MaxEnt Maximum Entropy algorithm and showed a wide suitable area for the potential distribution of each species in the neotropics, as promoted by climate variability. This work is important both in the scientific and technical communities since it allows for the direction of economic and environmental policies, creating early warning systems, mitigating the impacts of these pest species in the fruit and vegetable economies of Latin America and Panama, and contributing to improving the social economy and food security.

Abstract

Climate variability has made us change our perspective on the study of insect pests and pest insects, focusing on preserving or maintaining efficient production systems in the world economy. The four species of the genus Anastrepha were selected for this study due to their colonization and expansion characteristics. Models of the potential distribution of these species are scarce in most neotropical countries, and there is a current and pressing demand to carry out this type of analysis in the face of the common scenarios of climate variability. We analyzed 370 presence records with statistical metrics and 16 bioclimatic variables. The MaxEnt method was used to evaluate the effect of the ENSO cycle on the potential distribution of the species Anastrepha grandis (Macquart), Anastrepha serpetina (Wiedemann), Anastrepha obliqua (Macquart), and Anastrepha striata (Schiner) as imported horticultural pests in the neotropics and Panama. A total of 3472 candidate models were obtained for each species, and the environmental variables with the greatest contribution to the final models were LST range and LST min for A. grandis, PRECIP range and PRECIP min for A. serpentina, LST range and LST min for A. obliqua, and LST min and LST max for A. striata. The percentage expansion of the range of A. grandis in all environmental scenarios was 26.46 and the contraction of the range was 30.80; the percentage expansion of the range of A. serpentina in all environmental scenarios was 3.15 and the contraction of the range was 28.49; the percentage expansion of the range of A. obliqua in all environmental scenarios was 5.71 and the contraction of the range was 3.40; and the percentage expansion of the range of A. striata in all environmental scenarios was 41.08 and the contraction of the range was 7.30, and we selected the best model, resulting in a wide distribution (suitable areas) of these species in the neotropics that was influenced by the variability of climatic events (El Niño, Neutral, and La Niña). Information is provided on the phytosanitary surveillance systems of the countries in areas where these species could be established, which is useful for defining policies and making decisions on integrated management plans according to sustainable agriculture.

Graphical Abstract

1. Introduction

The fluctuations that are perceived in the climate (temperature, precipitation, and humidity) over relatively short periods are the cause of climate variability [1,2,3,4,5]. Previous studies have shown that climate variability affects physical and biological systems, and is of particular importance in cultivation areas, with changes in the frequency, intensity, spatial extent, duration, and timing of heavy rains and periods of drought during El Niño and La Niña phenomena. These variations can alter the behavior of ecosystems and cause the displacement and extinction of species [6,7,8,9]. El Niño–Southern Oscillation (ENSO) is a natural phenomenon involving fluctuating ocean temperatures in the equatorial Pacific [10]. Varying from 2 to 7 years, the ENSO cycle is the prevailing form of variability in the Pacific Ocean [6,11]. This variability causes uncertainties in productive agroecosystems (fruit and vegetables), in turn causing social and environmental stress that puts food security on the planet at risk [6,12]. It is expected that, in the next few years, climatic variability will cause changes in the geographic distribution of insects because of the reordering of climatic zones [13,14]. Therefore, the magnitudes of the impacts will be associated with climate variability. Regions of Latin America, South America, and the Caribbean are among the most vulnerable areas to climate variability, since most of the species that live there are endemic or restricted to a specific tropical ecosystem [15,16,17,18]. Most living beings are susceptible to the effects of climate variability due to specific physiologies and phenological qualities [19]. Poikilothermic organisms (ectotherms), such as insects, whose body temperatures vary according to the surrounding climate, are affected by unstable weather [20]. Climatic variations in temperature and precipitation can alter the ecological interactions of Anastrepha species, allowing species to colonize higher elevations and extend the ranges of their geographical distributions [21,22,23]. Two of the most commonly used methods to assess the effects of climatic variations in species distribution patterns have been the ecological niche (ENM) and potential distribution (SDM) models [24,25,26,27]. In addition to characterizing the environmental and geographic patterns of the species, these methods have been applied to multiple research questions in various areas of knowledge, including conservation, epidemiology and management of invasive species, eradication programs, and identification of cultivable areas [28,29,30].
The Tephritidae (Diptera) family consists of about 4700 species grouped into approximately 500 genera [31,32], distributed in temperate, tropical, and subtropical regions worldwide [33]. Anastrepha Schiner (Diptera: Tephritidae: Toxotrypa-nini) is the largest genus of Tephritidae from the Americas and includes more than 250 species distributed from the southern United States (Texas and Florida) to northern Argentina. [34,35]. At least seven species of the genus Anastrepha are considered economically important insects because they directly affect the production of cultivated fruit (Mangifera indica L. and Citrus sp.) and have a wide range of hosts [36]. Among them, Anastrepha fraterculus (Wiedemann 1830), Anastrepha grandis (Macquart 1846), Anastrepha ludens (Loew 1873), Anastrepha obliqua (Macquart 1835), Anastrepha serpentina (Wiedemann 1830), Anastrepha striata (Schiner 1868), and Anastrepha suspensa (Loew 1862) stand out. The species A. fraterculus is recognized as a cryptic species complex [37,38,39]. The damage caused by some species of the genus Anastrepha can be devastating to fruit and its commercialization; total losses of up to 90% of the crop have been documented in recent decades [40,41,42,43,44]. The potential distribution of fruit flies as part of the management of pest insects has been considered important in recent years, with the management of information bases, the use of models, and the implementation of new technologies in the fields of climate, species biodiversity, and better control planning in productive zones in regions of the United States [45], Europe [46,47], and Latin America [48,49,50].
Over the past decade, the world has seen at least 690 million people undernourished, 750 million people suffering from food insecurity, 2 billion people lacking access to safe and nutritious food, and 3 billion people unable to afford to pay the cost of a healthy diet [15]. Within the agricultural production sector, tropical fruit constitutes a relatively new group in the world trade of basic products, having acquired importance in the international market since 1970 due to advances in transportation, trade agreements, and changes in consumer fruit preferences [51]. With the declaration of the International Year of Fruits and Vegetables in 2021, healthy diets and lifestyles have been promoted among consumers in an environmentally sustainable way, recommending a minimum intake of 400 g of fruits and vegetables per person per day, which is only met in some parts of Asia and in countries with high incomes and good fruit availability [52,53]. In Panama, growing fruit has enormous potential for development due to the climate and soil conditions in various regions of the country, as well as the opportunities offered by international markets [54]. Of the total planted area (227,551 ha−1), as coverage of interest to the producer, 13% (29,581 ha) belongs to fruit trees [55].
Evaluating the effect of the ENSO cycle on the potential distribution of four species of the genus Anastrepha, of economic importance in horticulture, is of great relevance for Central America and Panama, and the scientific, economic, social, and environmental fields; it can provide information for decision-making in countries of the American continent that present the greatest vulnerability (exports and quarantines) for the establishment and distribution of each species. In addition, the information generated constitutes the basis for risk analysis and the implementation of preventive control and mitigation plans, which will reduce socioeconomic impacts in the potentially affected sectors of Latin America, Panama, and the world.

2. Materials and Methods

2.1. Study Area

The American neotropics, which extend between the Tropics of Cancer and Capricorn, was used as the study area [56,57]. This included the south of Mexico, all of Central America, the Caribbean, and a large part of South America, including the great Amazon jungle [58]. The natural landscape of these bioregions is made up of tropical forests, wetlands, savannahs, temperate grasslands, deserts, and Andean herbaceous formations [59,60,61,62]. Subsequently, layers were cut for the Republic of Panama, which is geographically located in the low northern latitudes (7°12′07″ and 9°38′46″ north latitude), between 77°09′24″ and 83°03′07″ western longitude, and characterized as the narrowest and longest country of the Central American isthmus. It has a land surface of 74,177.3 km2 (not including areas of continental water masses of 1,142,506.9 km2), is arranged in a west–east direction, and is bordered to the north by the Caribbean Sea, to the south by the Pacific Ocean, to the east by Colombia, and to the west by Costa Rica [63,64] (Figure 1).

2.2. Presence Records

As a model for this study, four species of the genus Anastrepha that are of horticultural importance in the neotropics and Panama were selected: A. grandis, A. serpentina, A. obliqua, and A. striata [65,66]. The various hosts, direct fruit damage, some polyphagous species, and distribution within the region were recorded [67]. Species presence records from the following sources were used: the Global Biodiversity Information Facility (GBIF; https://www.gbif.org/es/, accessed on 5 March 2022, Species Link (https://specieslink.net, accessed on 16 June 2022, Center for International Agricultural Bioscience (CABI; https://www.cabi.org/, accessed on 3 March 2022 and the National Plant Health Directorate of the Ministry of Agricultural Development (MIDA; https://mida.gob.pa/sanidad-vegetal, accessed on 10 March 2022. (Table S1). All presence points were reviewed in Excel, eliminating duplicate records and those with coordinate errors, and using a shapefile of the neotropical polygon [68]. Once the database was formed, columns with the scientific name and longitude and latitude coordinates were ordered for each species. The data were refined using the thinning function with the spThin package in R [69], using a 30 km buffer to carry out the spatial thinning of records and taking the dispersal capacity reported for the genus Anastrepha as a reference [70,71,72]. Each final dataset was divided into presence records for calibration and evaluation (80% and 20%, respectively).

2.3. Climate Data

To characterize this phenomenon, the following specialized agencies were consulted for their records of climatic episodes occurring in the Pacific Ocean: NOAA in the United States (National Weather Service, Los Angeles, CA, USA, 2018), the Australian Government Meteorological Office (Bureau of Meteorology, Melbourne, Australia, 2018), and the Tokyo Climate Center in Japan (Japan Meteorological Agency, Tokyo, Japan, 2019). The information available from 2000 to 2019 was integrated; from the consensus of these three agencies, six episodes of La Niña, five of El Niño, and three Neutral episodes were obtained. Each of these episodes was characterized by four rasters (minimum, maximum, mean, and range), for which each of the environmental layers included the improved vegetation index (EVI, monthly Modis-Terra MOD 11C2v006), the temperature of the land surface (LST, monthly Modis-Terra MOD 11C3v006), near-real-time precipitation rate (NRTPR, 3 h TRMM 3B42RTv7), and normalized difference vegetation index (NDVI, monthly Modis-Terra MOD 11C2v006). From these, 16 environmental layers were obtained for each episode at a spatial resolution of 0.25° or 25 km at the equator [73,74,75].
With these variables, 8 datasets were built using the following criteria for their selection: all variables (set 1); analysis of the Pearson correlation coefficient to reduce collinearity (set 2) [76], where variables with a correlation value >|0.8| were removed using the corrplot package of the statistical software R 3.6.0 [77,78]; Jackknife analysis in MaxEnt to evaluate the individual contribution of variables without spatial autocorrelation to the models, which included variables that contribute ≥80% (set 3) [78]; variables with a variance inflation factor (VIF) <10 (set 4); all variables related to NDVI (set 5); all variables related to EVI (set 6); all variables related to LST (set 7); and all variables related to precipitation (set 8) [78,79,80,81] (Table S2).

2.4. Construction of the Models

For each species, a calibration area was proposed for the model based on the spatial distribution of occurrences to define the environmental sampling area for the algorithm [82,83]. These areas were demarcated by taking the records of the presence of the species at the intersections of the polygons of the ecoregions proposed by the World Wide Fund for Nature (WWF) [56,84]. Using the MaxEnt program, which uses an automatic learning algorithm that combines Bayesian, maximum entropy, and statistical models, all possible models resulting from combinations of the variables were built, estimating the probability distribution of the maximum entropy of presences in a cell, determined according to the setting [85,86]. This probability distribution is the result of calculating the response curves, where suitability is described based on each of the model’s environmental variables [87]. All models were generated with the kuenm package in R version 1.1.5 using the method of exhaustive selection of variables [29,88,89], a sequence of 17 regularization multipliers (0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 2, 3, 4, 5, 6, 8, and 10), 8 sets of variables, and a combination of entity classes (linear, quadratic, product, and their combination). Subsequently, the performance of the candidate models was evaluated with the following metrics: (a) statistical significance by partial receiver operating characteristic curve (pROC with 500 iterations and 50 percent of the data for bootstrapping); (b) evaluation of the omission rate in the test presence records, selecting models with results lower than 5%; and (c) selection of models with less complexity to the Akaike information criterion corrected for small samples, and the difference between the pretender model and the outstanding model (AICc and ∆AICc, respectively). The models were evaluated by selecting those that were statistically significant (criterion a) and reducing the number of models by incorporating only those that executed the omission rate (criterion b). Among the significant and low omission candidate models, models with a delta less than 2 (criterion c) were selected [76]. Based on the results, each final model was built from 10 repetitions, and individual response curves were obtained. The type of output of the final model was selected (extrapolation, free extrapolation, and without extrapolation and subjection) using Jackknife analysis of the three variables with the greatest contribution to the models and graphing them in a three-dimensional space. In this environmental space, the location of the points of occurrence is compared with the environmental space of area M [77,81]. We reclassified the maps to binary by a threshold that omits all regions with habitat suitability below the suitability values for the lowest 10% of occurrence records [50] using the reclassify to binary function in SDMToolbox V. 2.4.

2.5. Analysis of Geographic Space

Each final model by species was projected to the geographic space of the area of the neotropics, and later the map was cut to the area of Panama, estimating the potential distribution of each species in the different average conditions investigated (Neutral, El Niño, and La Niña). These projections were reclassified as binary (1 = suitable; 0 = not suitable) [90,91,92], calculating a threshold for each species that omitted all regions with habitat suitability below the suitability values for the lower 10% of occurrence records. Additionally, we implemented a mobility-oriented parity (MOP) analysis that compared the environmental similarity of the projection area with the model calibration area to identify areas where there was strict extrapolation [93,94,95,96,97,98].

3. Results

A total of 370 presence records were obtained for the four Anastrepha species, in which A. grandis had the lowest number of records (39), distributed between Panama, Colombia, Ecuador, Venezuela, Peru, Bolivia, Brazil, Paraguay, Uruguay, and northern Argentina. Anastrepha serpentina had 88 records distributed among Mexico, Belize, El Salvador, Guatemala, Costa Rica, Panama, Colombia, Ecuador, Venezuela, French Guyana, Peru, Brazil, Bolivia, Paraguay, and northern Argentina. Anastrepha obliqua had 93 records distributed in Mexico, Guatemala, Belize, El Salvador, Nicaragua, Costa Rica, Panama, Colombia, Venezuela, Peru, Ecuador, Suriname, French Guiana, Brazil, Bolivia, Paraguay, and northern Argentina. Finally, A. striata had 150 records, mainly distributed in Mexico, Belize, Honduras, Guatemala, Nicaragua, Costa Rica, Panama, Colombia, Ecuador, Peru, French Guyana, Suriname, Bolivia, Brazil, and Paraguay. A total of 3472 candidate models were obtained for each Anastrepha sp. Comparing the locations of the presence records in the environmental space and the environmental space available in the calibration area (M) and the projection area for the four species, the records were dispersedly distributed in the environmental space of the M. This result allowed the selection of all species in the models with output and extrapolation, considering the methodology proposed by Owens et al. [93] (Figure S1). After analyzing all candidate models, the best of each species was selected according to the criteria outlined in the Materials and Methods section (Table 1).
The environmental variables with the greatest contributions to the final models were LST range (32%) and LST min (28%) for A. grandis, PRECIP range (58%) and PRECIP min (18%) for A. serpentina, LST range (45%) and LST min (28%) for A. obliqua, and LST min (47%) and LST max (24%) for A. striata.

3.1. Potential Distribution in the Neotropics

Under El Niño, Neutral, and La Niña conditions, the models for the four species showed wide areas with environmentally suitable conditions (adequate for their establishment) in coastal, savannah, montane, temperate, and tropical forest zones. For A. grandis, the potential distribution models in El Niño (Figure 2A) and La Niña (Figure 2C) events were very similar, with suitable areas in northern and southeastern to western Mexico, Guatemala, Belize, northwestern Honduras and Nicaragua, southwestern Guatemala, Costa Rica, the Atlantic zone and western Panama, northeastern Colombia and Ecuador, southeastern Venezuela, northeastern Guyana, northeastern Peru, northern Bolivia, much of the northern and southern Amazon in Brazil, southeastern Paraguay, Uruguay, northern and southern Argentina, and southern Chile. For Neutral events (Figure 2B), a larger suitable area was evident in northern Mexico, the northern Amazon in Brazil, central Uruguay and Paraguay, northern Argentina, and several countries in Central America. For A. serpentina, the El Niño (Figure 3A) and La Niña (Figure 3C) episodes have almost the same pattern, where an increase in the suitable areas was observed in the northeastern areas of Mexico during the La Niña event and a decrease from the ideal area in southern Guatemala, northern Honduras, and El Salvador, and an increase in Costa Rica and Panama. In South America, there was an increase in the ideal areas in a large part of Colombia, Venezuela, Ecuador, Brazil, Bolivia, Paraguay, and northern Argentina and a reduction in the suitable area in La Niña in southern Brazil, northern Argentina, and northwestern Chile. In the Neutral event, there was a larger suitable area in the entire neotropical region except in the central and western areas of Mexico, a small area of northwestern and southern Colombia, east-central Brazil, the coastal area of Peru and Chile, and northeastern and southern Argentina (Figure 3B).
For A. obliqua, in the El Niño and La Niña events, the ideal area increases in a large part of the neotropics, in Mexico, Central America, and a large part of the South American territory (Figure 4A,C). In Neutral events, it maintains a pattern similar to that of El Niño and La Niña, except for the increase in the suitable area that is observed in central Mexico (Figure 4B). Lastly, for A. striata, the El Niño and La Niña events (Figure 5A,C) presented a more suitable area in Mexico, Guatemala, Belize, El Salvador, Nicaragua, Costa Rica, Panama, Colombia, Venezuela, central Brazil, Argentina, Paraguay, Uruguay, and part of Chile. In Neutral events, the increase was maintained in central Brazil and southern Argentina (Figure 5B).

3.2. Potential Distribution in Panama

In the case of Panama, the potential distribution with suitable areas for A. grandis in the El Niño event was in the eastern (Darién, Panama) and northern parts of the country (Colon, Coclé, Veraguas, and Bocas del Toro). In Neutral events, the ideal area included the provinces of Darién and Panama to the east, Panama West, the Atlantic area of Colon, Veraguas, and Bocas del Toro, and the central areas of Herrera, Coclé, Veraguas, and Chiriquí. The ideal areas remained in eastern and western Panama, Darién, the Atlantic zone of Colon, Veraguas, and Bocas del Toro, and the Pacific zone of Veraguas and central Chiriquí (Figure 6A). For A. serpentina and A. obliqua, the ideal areas covered all of Panama in all events (Figure 6B,C). For A. striata, in El Niño, the ideal area covered almost all of Panama, except for isolated points in Colon and Bocas del Toro. In Neutral events, the ideal areas were located in the central regions of the provinces of Darién, Panama, Coclé, Veraguas, Chiriquí, and Bocas del Toro, and to the south in Los Santos. In La Niña events, the ideal areas were eastern Darién, central Panama, and the provinces of Coclé, Herrera, Los Santos, Veraguas, Chiriquí, and Bocas del Toro (Figure 6D).

3.3. Movement-Oriented Parity Analysis (MOP)

When identifying the risk of extrapolation of El Niño, Neutral, and La Niña conditions, all species presented similar patterns, where most of the projection area presented environmental similarities with the calibration area (Figure 2D–F). Specifically, for A. grandis, in El Niño events, there was a risk of extrapolation in northern and central Mexico, southern Peru, northern Chile, northeastern Argentina, and east-central Brazil. In Neutral events, there was a risk of extrapolation in central Mexico, southern Peru, northern Chile, northeastern Argentina, and Brazil. For La Niña, we identified areas of extrapolation risk in western and eastern Mexico, east-central Brazil, southern Peru, northern Chile, and northeastern Argentina. For A. serpentina, in El Niño events, there was no risk of extrapolation. For Neutral events, there was a risk of extrapolation in west-central Mexico, southern Peru, and parts of southern Argentina. For La Niña events, there was a risk of extrapolation in southwestern Mexico, northern and southern Guatemala, northern El Salvador, northern Belize, northern Colombia and Venezuela, and east-central Brazil (Figure 3D–F). For A. obliqua, in El Niño events, there was a risk of extrapolation in west-central Mexico, southern Peru, northern Chile, northern Argentina, and east-central Brazil. For Neutral events, there was a risk of extrapolation in central Mexico and northern Chile, and for La Niña events, there was a risk of extrapolation in a large part of Mexico, northern Venezuela, a large part of east-central and southern Brazil, Argentina, Paraguay, Uruguay, southern Peru and Bolivia, and Chile (Figure 4D–F). For A. striata, in El Niño events, there was a risk of extrapolation in northwestern Mexico, northern Venezuela, east-central Brazil, southern Peru, and on the coasts of Chile and Argentina. In Neutral events, there was a risk of extrapolation at only a few points in northern Mexico, northern Chile, and northern Argentina. In La Niña events, there was a risk of extrapolation in a large part of Mexico, northern Venezuela, east-central and southern Brazil, Paraguay, Uruguay, Argentina, southern Peru and Bolivia, and Chile. In the remaining points not mentioned by events, there was no risk of extrapolation (Figure 5D–F).
In the MOP analysis, in Panama, there was no risk of extrapolation in most of the projected models of the species, except for A. grandis in Neutral events in Darién Province (Figure 2E), A. serpentina in El Niño events in Darién Province (Figure 3D), and La Niña in the provinces of Herrera and Los Santos (Figure 3E).

3.4. Analysis of Changes in Environmental Suitability in Different Scenarios

For A. grandis, when comparing the Neutral vs. El Niño scenarios, there was a 14.72% expansion range to northern Mexico, southwestern Brazil and Colombia, southern Peru, and a large part of Chile. In addition, 61.53% of the area was shown as an area of no change (occupied area in both scenarios), mainly for the Pacific zone, southern Peru, Chile, and Argentina. In the comparison of Neutral vs. La Niña events, there was a 16.47% contraction in northern Mexico, southeastern Venezuela, east-central Brazil, southern Bolivia, and northern Paraguay and Argentina, with a 59.38% chance of no change in the occupied area (Figure 7A). For A. serpentina (Figure 7B), a 14.11% contraction range was observed in Neutral vs. El Niño events in northern Mexico, northeastern Venezuela, east-central Brazil, central Bolivia, Paraguay, and Argentina, with a 74.57% chance of no change in the occupied area. In Neutral vs. La Niña events, there was a 14.38% contraction range in northern Mexico, central Colombia and Peru, and southeastern Brazil, Uruguay, and Argentina, with a 74.29% chance of no change in the occupied area. In A. obliqua (Figure 7C), in Neutral vs. El Niño events, we observed a 9.77% unoccupied area in parts of Colombia, Ecuador, Peru, Bolivia, Chile, and Argentina, and an 86.4% chance of no change in the occupied area. In Neutral vs. La Niña events, there was an 8.04% unoccupied area in Colombia, Ecuador, Peru, Bolivia, Chile, Argentina, and southern Brazil, with an 86.59% chance of no change in the occupied area. For A. striata (Figure 7D), in Neutral vs. El Niño events, we observed an 18.91% range expansion in a large part of Mexico, northern Guatemala and Belize, and central Brazil, in addition to a 66.84% chance of no change in the occupied area. In Neutral vs. La Niña events, there was a 22.17% expansion range in Mexico, northern Guatemala and Belize, and central Brazil, and a 68.44% chance of no change in the occupied area (Table 2).

4. Discussion

In recent years, the study of new technologies has allowed strategic plans to be made within productive systems, where the management of information on biodiversity is important, since these elements allow prevention and control plans to be organized in various countries, improving food security in society [99,100,101]. The results of this study indicate that the studied species of the genus Anastrepha can increase in the native areas at the neotropical level. In ENSO events (El Niño, Neutral, and La Niña), the information generated is important due to the limited research on species of this genus [102,103]. For A. grandis, there are reports of its presence from eastern Panama to a large part of South America, for which it is important to know the places with ideal environmental conditions for its distribution to the rest of Central America [104]. In El Niño and La Niña events, the geographic spaces with ideal environmental conditions for potential distribution were similar, where the LST variable (range and min) was shown to be the most important in the models. This result could be related to the stability of soil moisture and temperature in tropical forested areas at the time of oviposition, hatching, and the availability of fruit hosts and wild host plants [104,105]. The results of this study showed that areas with adequate conditions for these species of fruit flies corresponded to localities with temperate and sub-humid climates in Central America and northern Mexico. Therefore, the government entities in charge of the control of fruit flies and other pest insect species should take measures that allow the early arrest of these species and thus reduce the risk of the establishment of reproductive populations [104]. During Neutral events, the model showed an increase in geographic spaces with adequate conditions in the center of the Amazon. This could be related to fluctuations in average and low temperatures in this region, which benefit the biological processes of hatching and larval development of A. grandis [105]. In Panama, this species could inhabit tropical forests, where they could expand toward suitable productive zones, according to the results [104,105,106].
The results of the MOP analysis for A. grandis showed a high environmental similarity between the calibration and projection areas in the different episodes that could favor the latitudinal and altitudinal expansion of the geographic distribution of the species. This pattern is potentially due to high climate variability in the average temperatures in the regions of Central America and central South America where there is a greater concentration of tropical forests, and with little extrapolation irrigation in central Mexico and Panama due to their arid zones [106,107,108]. When comparing Neutral vs. El Niño events, the potential distribution of A. grandis could expand because the biology of the species responds well to average temperatures in northern Central America, where there are savannahs and temperate forests, similar to South American regions with temperate and semi-arid savannah forests in which mild temperatures and little rainfall favor the expansion of this species [109,110]. The changes observed in Neutral vs. La Niña events would cause a contraction in the species possibly associated with high temperatures in northern Mexico, due to its dry areas and low humidity, affecting the species’ reproduction in the same way as in South America due to the semi-arid conditions with high temperatures in western Amazonia and semi-arid conditions in Argentina [111,112].
For A. serpentina, the potential distribution in El Niño and La Niña episodes could be interpreted according to the model based on the climate variable PRECIP (range and min), where in arid areas with little rainfall, such as northern Mexico and the Southern Cone countries, there is little probability of establishing this species [113,114,115]. In Neutral episodes, the adequate area could increase due to the average rainfall in the neotropics since the biological processes associated with the species adapt better with adequate field humidity [114]. When projecting MOP analyses in Neutral and El Niño events, the environmental conditions could be similar since by maintaining the average rainfall regime in the neotropics, the species increases its distribution [113]. In La Niña events, the environmental similarity is lower since the increase or variation in the amount of rain in east-central Brazil affects the reproduction of the species due to the death of males [114,115,116]. In Panama, suitable zones with environmental similarities are high, which would not affect the distribution patterns of the species, where the natural landscape and orography probably maintain good humidity conditions in the environment [117]. When comparing the Neutral vs. El Niño events, the potential distribution could contract in the north of Mexico, part of the coastal areas of the central Pacific, and the Southern Cone due to their arid conditions or lack of rainfall, affecting the hatching and lifespan of the species [117,118]. When comparing Neutral and La Niña events, there could be a contraction in the distribution due to the low rainfall produced on the Mexican Atlantic coast and a large part of the Southern Cone, which affects egg oviposition and pupal development due to a lack of humidity [40,118,119]. For A. obliqua, one of the most polyphagous species of this genus, the climate variable LST (range and min), according to the model, could influence the potential distribution in El Niño and La Niña events, since it adapts to a wide range of temperatures, colonizing almost the entire neotropics, with countless hosts, limited only by high temperatures and competition between species and predators, which is why their development can be delayed or accelerated [120]. At an ecological level, the areas with high elevation, as well as cold climates, reduce their distribution, such as the Cordillera of Mexico, the Andes in South America, and the coldest parts of the Southern Cone [121,122]. In Neutral events, the possible extension of suitable distribution areas in northern Central America could be due to the stability in temperatures and the greater range of fruit hosts, which benefit the reproduction of this species [122]. Observing the MOP analysis, in El Niño, Neutral, and La Niña events, the areas with the greatest environmental similarity between the calibration and projection areas were from Neutral events, which favor the greater distribution of the species due to its characteristics of adaptation to variations in temperatures [120]. In El Niño and La Niña events, the environmental similarity is reduced in parts of Mexico due to their hot, dry, and very humid climates, in eastern and southern Brazil due to dry climates, and in Argentina due to its cold climate, which harms the oviposition of the eggs of this species [120,122]. In Panama, the area of potential distribution and environmental similarity presents adequate zones in its geography, which could promote reproduction (variety of hosts and ideal conditions) in the country due to its suitable temperature (30 °C) and humidity (98%) [123,124]. The changes observed for Neutral vs. El Niño events showed that the potential distribution of the species would not include the Andes or Argentina; this could be related to low temperatures and the few natural hosts present in these areas [40,45]. In the comparison of Neutral vs. La Niña, both presented the same pattern as El Niño, with their potential distribution areas not including those with cold climates, affecting all reproduction patterns [120,121,122,123,124,125].
For A. striata, better known as the guava fly, the potential distribution could be influenced by the climate variable LST (min and max), according to the projected model, since El Niño and La Niña events present similar distribution patterns. We estimate that the species adapts to various ranges of temperatures between 18 °C and 30 °C and good soil moisture throughout the neotropics, except for the mountainous regions of the Andes and the wetlands of Brazil [126,127,128]. In Neutral events, the areas of potential distribution are reduced in central Mexico due to its desert areas with high temperatures and in northern Brazil, in the jungle areas of the Amazon with high temperatures and high humidity, which affect the hatching and death of adults [126]. In Panama, for El Niño events, the area of potential distribution is wide, which tells us that the species adapts to several not very humid tropical temperature ranges. In Neutral and La Niña events, the areas of potential distribution are limited only to coastal areas due to their high temperatures and low humidity, which harms the reproduction of the species [128]. When reviewing the MOP maps in Neutral events, there is a high environmental similarity compared to the area of distribution. The species may adapt to variations in climate, temperatures, and tropical humidity, except in very high areas [126,127]. During El Niño and La Niña, there is a high risk of extrapolation of the species in northwestern Mexico and southern Brazil, Uruguay, Paraguay, Chile, and Argentina, since the areas of environmental similarity are low, and the species cannot be present due to the climatic conditions of adverse temperatures for its reproduction [102]. The changes observed for Neutral vs. El Niño and Neutral vs. La Niña events could present an increase in the expansion areas in central Mexico due to its semi-warm and tropical conditions, places where temperatures are high, such as desert areas and tropical coastal areas, dry areas in Venezuela, and tropical forests in east-central Brazil, which demonstrates the great adaptability of A. striata in its adult phase at high temperatures [126,127]. The same happens in Panama, where the range of expansion is toward the tropical mountain areas on the Atlantic and Pacific coasts, where temperatures and humidity are high, which benefits the reproduction of the species during its dispersal phase in search of hosts [20,126,127,128,129,130].

5. Conclusions

There is evidence that indicates that changes in temperature and rainfall induced by the ENSO cycle cause changes in the potential distribution of insect species by influencing the community structure and population dynamics. Our results show that the climatic variability caused by the event in turn causes changes in the expansion of Anastrepha spp. of horticultural importance, occupying an extensive area where favorable and unfavorable zones can be altered in a varied and complex way in the neotropical region. The anomalies in precipitation and temperature induced by the El Niño and La Niña phenomena generate changes in the potential distribution of species in the region, and an increase or reduction in the environmentally suitable zone is projected. The applied modeling methodology is efficient and makes the simulation of the system more accessible, indicating areas of greater risk. Information is provided for use in phytosanitary surveillance systems of the countries in areas where these fruit fly pests could be established, to define policy and decision-making on integrated management plans according to sustainable agriculture principles.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects14080714/s1, Table S1: Anastrepha spp. species database, Table S2: The main environmental variables that contribute to the model for A. grandis, A. serpentina, A. obliqua, and A. striata, Figure S1: Analysis of the environmental space represented by presence records, calibration area, and projection area.

Author Contributions

Conceptualization, A.B.D., R.A.V., A.B.A., J.Á.J. and M.A.-S.; methodology, A.B.D., R.A.V., A.B.A., J.Á.J. and M.A.-S.; validation, A.B.D., R.A.V., A.B.A., J.Á.J. and M.A.-S.; formal analysis, A.B.D., R.A.V., A.B.A., J.Á.J. and M.A.-S.; investigation, A.B.D., R.A.V., A.B.A., J.Á.J. and M.A.-S.; data review, A.B.D., R.A.V., A.B.A., J.Á.J. and M.A.-S.; writing—original draft preparation, A.B.D.; writing—proofreading and editing, A.B.D., R.A.V., A.B.A., J.Á.J. and M.A.-S.; supervision, R.A.V., A.B.A., J.Á.J. and M.A.-S.; activity management, A.B.D., R.A.V., A.B.A. and M.A.-S.; funding acquisition, A.B.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financed by SENACYT—IFARHU-Panama, as part of the doctoral scholarships project for public officials in the agricultural sector. Contract number 270-2021-117.

Data Availability Statement

All data are contained in the article and the Supplementary Materials.

Acknowledgments

First of all, to God for letting me live this moment and thanking everyone who contributed to this work; the group of colleagues from the IDIAP-Panama Agroclimatic Unit Project; Rubiela García, MIDA-Panama; Milagros Castillo, University of Panama; Agustin Merino, University of Santiago de Compostela; and the group of directors of the research work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Montealegre, J.; Pabon, J. La Variabilidad Climática Interanual asociada al ciclo El Niño-La Niña–Oscilación del Sury su efecto en el patrón pluviométrico de Colombia. Meteorol. Colombina 2000, 2, 7–21. [Google Scholar]
  2. Trenberth, K.E.; Fasullo, J.T.; Kiehl, J. Earth’s global energy budget. Bull. Am. Meteorol. Soc. 2009, 90, 311–324. [Google Scholar] [CrossRef]
  3. Organización de las Naciones Unidas. Convención Marco De Las Naciones Unidas Sobre El Cambio Climático; Organización de las Naciones Unidas: NewYork, NY, USA, 2015; p. 40. Available online: https://unfccc.int/resource/docs/2015/cop21/spa/l09s.pdf (accessed on 1 March 2023).
  4. D’Antoni, H.L. Cambio global. Procesos naturales e intervención humana. Acta Bioquímica Clínica Latinoam. 2012, 46, 1–76. [Google Scholar]
  5. Moyano, E.; Paniagua, Á.; Lafuente, R. Políticas ambientales, cambio climático y opinión pública en escenarios regionales. El caso de Andalucía. Rev. Int. Sociol. 2009, 67, 681–699. [Google Scholar] [CrossRef]
  6. Intergovernmental Panel on Climate Change. Climate Change, Impact, Adaptation and Vulnerability, Summary for Policymakers; Intergovernmental Panel on Climate Change: Geneva, Switzerland, 2014; p. 40. Available online: https://www.ipcc.ch/site/assets/uploads/2018/02/SYR_AR5_FINAL_full_es.pdf (accessed on 3 March 2023).
  7. Gonzales, G.F.; Zevallos, A.; Gonzales, C.C.; Nuñez, D.; Gastañaga, C.; Cabezas, C.; Steenland, K. Environmental pollution, climate variability and climate change: A review of health impacts on the peruvian population. Rev. Peru. Med. Exp. Salud Publica 2014, 31, 547–556. [Google Scholar]
  8. Zotelo, C. Variabilidad climática y ciclos naturales. In Jornada sobre “Evolución y Futuro del Desarrollo de Producciones Agrícola-Ganaderas en el SO Bonaerense”; Universidad Nacional de la Plata: Bahía Blanca, Argentina, 2011; pp. 374–381. Available online: http://sedici.unlp.edu.ar/handle/10915/27824 (accessed on 5 March 2023).
  9. Tobón, C.J. Evaluación de los Impactos Potenciales de la Variabilidad Climática y el Cambio Climático en Algunos Indicadores para Seguridad Alimentaria en Zonas Productoras de Mercados Campesinos; Trabajo de Grado-Mestria, Universidad de Colombia: Bogotá, Colombia, 2014. [Google Scholar]
  10. Krishnamurthy, L.; Vecchi, G.A.; Msadek, R.; Murakami, H.; Wittenberg, A.; Zeng, F. Impact of strong ENSO on regional tropical cyclone activity in a high-resolution climate model in the North Pacific and North Atlantic Oceans. J. Clim. 2016, 29, 2375–2394. [Google Scholar] [CrossRef]
  11. Collins, M.; CMIP Modeling Groups. El Niño- or La Niña-like climate change? Clim. Dyn. 2005, 24, 89–104. [Google Scholar] [CrossRef]
  12. Balvanera, P.; Astier, M.; Gurri, F.D.; Zermeño, H.I. Resiliencia, vulnerabilidad y sustentabilidad de sistemas socioecológicos en México. Rev. Mex. Biodivers. 2017, 88, 141–149. [Google Scholar] [CrossRef]
  13. Gouveia, S.F.; Hortal, J.; Tejedo, M.; Duarte, H.; Cassemiro, F.A.; Navas, C.A.; Diniz, F.J.A. Climatic niche at physiological and macroecological scales: The thermal tolerance-geographical range interface and niche dimensionality. Blackwell Publishing. Glob. Ecol. Biogeogr. 2014, 23, 446–456. [Google Scholar] [CrossRef]
  14. Lehmann, P.; Ammunét, T.; Barton, M.; Battisti, A.; Eigenbrode, S.D.; Jepsen, J.U.; Kalinkat, G.; Neuvonen, S.; Niemelä, P.; Terblanche, J.S.; et al. Complex responses of global insect pestss to climate warming. Front. Ecol. Environ. 2020, 18, 141–150. [Google Scholar] [CrossRef]
  15. Comisión Económica para América Latina. La economía del Cambio Climático en América Latina y el Caribe: Paradojas y Desafíos del Desarrollo Sostenible. Programa EUROCLIMA Cambio Climático, Componente Socioeconómico 2015. (CEC/10/001). Available online: https://hdl.handle.net/11362/37311 (accessed on 10 March 2023).
  16. Sanchez, R.R. Respuestas Urbanas al Cambio Climático en América Latina. Documentos de Proyectos e Investigación. CEPAL-Naciones Unidas. Instituto Interamericano para la Investigación del Cambio Global No. 563. 2013. Available online: https://hdl.handle.net/11362/36622 (accessed on 15 March 2023).
  17. Magrin, G. Adaptación al Cambio Climático en América Latina y el Caribe. Documentos de Proyectos e Investigación. CEPAL-Naciones Unidas. Instituto Interamericano para la Investigación del Cambio Global No. 692. 2015. Available online: https://hdl.handle.net/11362/39842 (accessed on 15 March 2023).
  18. Rascón, V.A.E.; Cervantes, R.E. Vulnerabilidad social y clima extremo en estudios de América Latina. 2000–2019. Tlalli. Rev. Investig. Geogr. 2022, 8, 6–32. [Google Scholar] [CrossRef]
  19. Sheldon, K.S. Climate change in the tropics: Ecological and evolutionary responses at low latitudes. Annu. Rev. Ecol. Evol. Syst. 2019, 50, 303–333. [Google Scholar] [CrossRef]
  20. Régnière, J.; Powell, J.; Bentz, B.; Nealis, V. Effects of temperature on development, survival and reproduction of insects: Experimental design, data analysis and modeling. J. Insect Physiol. 2012, 58, 634–647. [Google Scholar] [CrossRef]
  21. Freeman, B.G.; Scholer, M.N.; Ruiz, V.; Fitzpatrick, J.W. Climate change causes upslope shifts and mountaintop extirpations in a tropical bird community. Proc. Natl. Acad. Sci. USA 2018, 115, 11982–11987. [Google Scholar] [CrossRef]
  22. Soberón, J. Grinnellian and Eltonian niches and geographic distributions of species. Ecol. Lett. 2007, 10, 1115–1123. [Google Scholar] [CrossRef]
  23. Powell, A.J.; Logan, A.J. Insect seasonality: Circle map analysis of temperature-driven life cycles. Theor. Popul. Biol. 2005, 3, 161–179. [Google Scholar] [CrossRef]
  24. Peterson, A.T.; Papeş, M.; Soberón, J. Rethinking receiver operating characteristic analysis applications in ecological niche modeling. Ecol. Modell. 2008, 213, 63–72. [Google Scholar] [CrossRef]
  25. Peterson, A.; Soberón, J.; Pearson, R.; Anderson, R.; Martínez, M.E.; Nakamura, M.; Araújo, M. Ecological Niches and Geographic Distributions (MPB-49); Princeton University Press: Princeton, NJ, USA, 2011. [Google Scholar] [CrossRef]
  26. Beckler, A.A.; French, B.W.; Chandler, L.D. Using GIS in areawide pests management: A case study in South Dakota. Trans. GIS 2005, 9, 109–127. [Google Scholar] [CrossRef]
  27. Aluja, M.; Birke, A.; Ceymann, M.; Guillén, L.; Arrigoni, E.; Baumgartner, D.; Pascacio, C.; Villafán, J. Agroecosystem resilience to an invasive insect species that could expand its geographical range in response to global climate change. Agric. Ecosyst. Environ. 2014, 186, 54–63. [Google Scholar] [CrossRef]
  28. Peterson, A.T. Uses and requirements of ecological niche models and related distributional models. Biodivers. Inform. 2006, 3, 59–72. [Google Scholar] [CrossRef]
  29. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Modell. 2006, 190, 231–259. [Google Scholar] [CrossRef]
  30. Aluja, M.; Mangan, R.L. Fruit fly Diptera: Tephritidae host status determination: Critical conceptual, methodological, and regulatory considerations. Annu. Rev. Entomol. 2008, 53, 473–502. [Google Scholar] [CrossRef]
  31. Norrbom, A.L. Host Plant Database for anastrepha and toxotrypana (Diptera: Tephritidae: Toxotrypanini). Diptera Data Dissemination Disk. CD—Not a Journal. 2004. Available online: https://www.ars.usda.gov/research/publications/publication/?seqNo115=108854. (accessed on 17 March 2023).
  32. Norrbom, A.L.; Neder, L.E. New neotropical species of Trupanea Diptera: Tephritidae with unusual wing patterns. Zootaxa 2014, 3821, 443–456. [Google Scholar] [CrossRef] [PubMed]
  33. Hernandez, O.V. El género Anastrepha Schiner en México Diptera: Tephritidae, taxonomía, distribución y sus plantas huéspedes. La Cienc. Y El Hombre 1992, 12, 190–191. [Google Scholar]
  34. Hernández, O.V.; Bartolucci, A.F.; Morales, V.P.; Frías, D.; Selivon, D. Cryptic species of the Anastrepha fraterculus complex: A multivariate approach for the recognition of South American morphotypes. Ann. Entomol. Soc. 2012, 105, 305–318. [Google Scholar] [CrossRef]
  35. Norrbom, A.L.; Korytkowski, C.A. New species of Anastrepha Diptera: Tephritidae, with a key for the species of the megacantha clade. Systematic Entomology Lab., USDA, ARS, c/o Smithsonian Institution. Zootaxa 2012, 34781, 11. [Google Scholar] [CrossRef]
  36. Cruz, B.; Bacca, M.L.; Nelson, A. Diversidad de las Moscas de las frutas Diptera: Tephritidae y sus parasitoides en siete minicipios del departamento de Narino. Boletín Científico. Cent. Museos. Mus. Hist. Nat. 2017, 21, 81–98. [Google Scholar] [CrossRef]
  37. Hernández, O.V.; Gómez, A.J.; Sánchez, A.; McPheron, B.; Aluja, M. Morphometric analysis of Mexican and South American populations of the Anastrepha fraterculus complex (Diptera: Tephritidae) and recognition of a distinct Mexican morphotype. Bull. Entomol. Res. 2004, 94, 487–499. [Google Scholar] [CrossRef]
  38. Selivon, D.; Perondini, A.L.; Morgante, J.S. A Genetic–Morphological Characterization of Two Cryptic Species of the Anastrepha fraterculus Complex (Diptera: Tephritidae). Ann. Entomol. Soc. 2005, 98, 367–381. [Google Scholar] [CrossRef]
  39. Silva, J.C.; Dutra, S.V.; Santos, M.S.; Silva, N.M.O.; Vidal, D.V.; Nink, R.V.; Guimarães, J.G.; Araujo, E.L. Diversity of Anastrepha spp. (Diptera: Tephritidae) and Associated Braconid Parasitoids from Native and Exotic Hosts in Southeastern Bahia, Brazil. Environ. Entomol. 2021, 39, 1457–1465. [Google Scholar] [CrossRef]
  40. Nolasco, N.; Iannacone, J. Fluctuación estacional de moscas de la fruta Anastrepha spp. y Ceratitis capitata (Wiedemann, 1824) (Diptera: Tephritidae) en trampas McPhail en Piura y en Ica, Perú. Acta Zool. Mex. 2008, 24, 33–44. [Google Scholar] [CrossRef]
  41. Aluja, M.; Rull, J.; Sivinski, J.; Allen, L.; Norrbom, A.; Wharton, R.M.; Díaz, F.; López, M. Fruit Flies of the Genus Anastrepha (Diptera: Tephritidae) and Associated Native Parasitoids (Hymenoptera) in the Tropical Rainforest Biosphere Reserve of Montes Azules, Chiapas, Mexico. Environ. Entomol. 2003, 32, 1377–1385. [Google Scholar] [CrossRef]
  42. Diznarda, S.B.; Flores, R.; Terrazas, G.G.; Leyva, R.E. Evaluación Económica de la Campaña Nacional Contra las Moscas de la Fruta en los Estados de Baja California, Guerrero Nuevo León, Sinaloa, Sonora y Tamaulipas; IICA-Ciudad de Mexico D.F: México city, Mexico, 2010; Available online: http://repiica.iica.int/docs/B2041e/B2041e.pdf (accessed on 17 March 2023).
  43. Valderrama, J.K.; Serrano, M.S.; Fischer, G. Mortalidad de larvas de Anastrepha fraterculus (Wiedemann) (Diptera: Tephritidae) en frutos de feijoa (Acca sellowiana [O. Berg] Burret) sometidos a un tratamiento cuarentenario de frio. Rev. Colomb. Entomol. 2005, 31, 171–176. [Google Scholar] [CrossRef]
  44. Saavedra, D.J.; Galeano, O.P.; Canal, N. Ecological relationships between host fruits, frugivorous flies and parasitoids in a fragment of tropical dry forest. Rev. Sci. Agric. 2017, 34, 32–49. [Google Scholar] [CrossRef]
  45. Sequeira, R.; Millar, L.; Bartels, D. Identification of Susceptible Areas for the Establishment of Anastrepha spp. Fruit Flies in the United States and Analysis of Selected Pathways; USDA-APHISPPQ Center for Plant Health Science and Technology: Raleigh, North Carolina, 2001; p. 47. Available online: https://www.aphis.usda.gov/plant_health/plant_pests_info/fruit_flies/downloads/isa.pdf (accessed on 18 March 2023).
  46. Godefroid, M.; Cruaud, A.; Rossi, J.P.; Rasplus, J.Y. Assessing the Risk of Invasion by Tephritid Fruit Flies: Intraspecific Divergence Matters. PLoS ONE 2015, 10, 8. [Google Scholar] [CrossRef]
  47. Fu, L.Z.H.; Huang, G.S.; Wu, X.X.; Ni, W.L.; Qü, W.W. The current and future potential geographic range of West Indian fruit fly, Anastrepha obliqua (Diptera: Tephritidae). Insect Sci. 2014, 21, 234–244. [Google Scholar] [CrossRef]
  48. Vázquez, P.; Escalona, A.H.; Segura, G.; Esparza, O.L.G. Modelación de la distribución geográfica potencial de dos especies de psitácidos neotropicales utilizando variables climáticas y topográficas. Acta Zool. Mex. 2014, 30, 471–490. [Google Scholar] [CrossRef]
  49. Jiménez, M.E.; Núñez, M.R.G.; Maradiaga, B.E.J. Distribución temporal de insectos asociados a maracuyá (Passiflora edulis Sims) en Matagalpa, Nicaragua. La Calera 2020, 20, 10–19. [Google Scholar] [CrossRef]
  50. Altamiranda, S.M.; Gutiérrez, D.J.; Araque, A.; Valencia, D.J.; Gutiérrez, R.; Martínez, R.A. Effect of El Niño Southern Oscillation cycle on the potential distribution of cutaneous leishmaniasis vector species in Colombia. PLoS Negl. Trop. Diseases 2020, 14, e0008324. [Google Scholar] [CrossRef]
  51. Food and Agriculture Organization of the Unites Nations, Fondo Internacional de Desarrollo Agricola, Organization Mundial de la Salud, Programa Mundial de Alimentos y Fondo de las Naciones Unidas por la Infancia. Versión Resumida de El Estado de la Seguridad Alimentaria y la Nutrición en el Mundo. Transformación de los Sistemas Alimentarios para que Promuevan Dietas Asequibles y Saludables. Roma. 2020. Available online: https://www.fao.org/documents/card/en/c/ca9699es (accessed on 18 March 2023).
  52. Food and Agriculture Organization of the Unites Nations; Centro de Cooperacion Internacional en Investigación Agrícola para el Desarrollo; Frutas y Hortalizas—Oportunidades y Desafíos para la Agricultura Sostenible a Pequeña Escala. Roma. 2021. Available online: https://www.fao.org/documents/card/en/c/cb4173es (accessed on 19 March 2023).
  53. Food and Agriculture Organization of the Unites Nations. Análisis del Mercado de las Principales Frutas Tropicales en Roma. 2022. Available online: https://www.fao.org/3/cb6897es/cb6897es.pdf (accessed on 19 March 2023).
  54. Instituto Interamericano de Cooperación para la Agricultura. La Fruticultura en Panamá: Su Potencial Socioeconómico e Iniciativas para su Desarrollo; Ministerio de Desarrollo Agropecuario, Instituto de Innovación Agropecuaria de Panamá, Ciudad de Panamá: Ancon, Panamá, 2008; Volume 167, p. 23. Available online: http://repiica.iica.int/docs/b0760e/b0760e.pdf (accessed on 19 March 2023).
  55. Ministerio de Desarrollo Agropecuario. Panamá. Cierre Agrícola Nacional. Dirección de Agricultura. Documento Panamá. 2020–21. Available online: https://mida.gob.pa/ (accessed on 20 March 2023).
  56. Olson, D.; Dinerstein, E.; Wikram, E.D.; Burgess, N.D.; Powell, G.V.N.; Underwood, E.C.; D’amico, J.A.; Itoua, I.; Strand, H.E.; Morrison, J.C.; et al. Terrestrial Ecoregions of the World: A New Map of Life on Earth. Biosci. J. 2001, 51, 933. [Google Scholar] [CrossRef]
  57. Dinerstein, E.; Olson, D.; Joshi, A.; Vynne, C.; Burgess, N.B.; Wikramanayake, E.; Hahn, N.; Palminteri, S.; Hansen, M.; Locke, M.; et al. An Ecoregion-Based Approach to Protecting Half the Terrestrial Realm. Biosci. J. 2017, 67, 534–545. [Google Scholar] [CrossRef]
  58. Myers, N.; Mittermeier, R.; Mittermeier, C. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853–858. [Google Scholar] [CrossRef] [PubMed]
  59. Metz, B.; Davidson, O.; Bosch, P.; Dave, R.; Meyer, L. Climate Change-Mitigation of Climate Change; OSTI Identifier: 21017235; Intergovernmental Panel on Climate Change, Geneva (Switzerland). Working Group III; U.S. Department of Energy Office of Scientific and Technical Information: Oak Ridge, TN, USA, 2007. Available online: https://www.osti.gov/biblio/21017235 (accessed on 3 June 2023).
  60. Jarvis, A.; Hijmans, R. The effect of climate change on crop wild relatives. Agric. Ecosyst. Environ. 2008, 126, 13–23. [Google Scholar] [CrossRef]
  61. Ramírez, O.; Ruiz, G.; Corral, J.A.; Pérez, M.C.; Villavicencio, G.R.; Mena, M.; Puga, N. Impactos del cambio climático en la distribución geográfica de Gossypium hirsutum L. en México. Rev. Mex. Cienc. 2014, 5, 1885–1895. [Google Scholar]
  62. Morrone, J.J. Biogeographical regionalisation of the Neotropical region. Zootaxa 2014, 3782, 1–110. [Google Scholar] [CrossRef] [PubMed]
  63. Jaén, S.O. Geografía de Panamá. Estudio Introductorio y Antología; Universidad de Panamá, Cuidad de Panama, Campus Central: Panama City, Panama, 1985; p. 472. Available online: https://books.google.com.pa/books?id=WlBrAAAAMAAJ&source=gbs_navlinks_s (accessed on 5 March 2023).
  64. Makay, A. Cien años de Geografía en Panamá. Universidad de Panamá. Articulo. 2003. Available online: http://bdigital.binal.ac.pa/bdp/artpma/cienanosdegeografia.pdf (accessed on 25 March 2023).
  65. Uchoa, A.M. Fruit Flies (Diptera: Tephritoidea): Biology, Host Plants, Natural Enemies, and the Implications to Their Natural Control; IntechOpen: London, UK, 2012; Volume 12, pp. 271–300. [Google Scholar] [CrossRef]
  66. Alvarado, G.L.; Medianero, E. Especies de parasitoides asociados a moscas de la fruta del género Anastrepha (Diptera: Tephritidae) en Panamá, República de Panamá. Scientia 2021, 25, 47–62. [Google Scholar]
  67. Uchôa, M.A.; Nicácio, J. New records of Neotropical fruit flies (Tephritidae), lance flies (Lonchaeidae) (Diptera: Tephritoidea), and their host plants in the South Pantanal and adjacent areas, Brazil. Ann. Entomol. Soc. 2010, 103, 723–733. [Google Scholar] [CrossRef]
  68. Cobos, M.E.; Jiménez, L.; Nuñez, P.C.; Romero, A.D.; Simoes, M. Sample data and training modules for cleaning biodiversity information. Biodiversity 2018, 13, 49–50. [Google Scholar] [CrossRef]
  69. Lammens, A.; Boria, R.A.; Radosavljevic, A.; Vilela, B.; Anderson, R.P. spThin: An R package for spatial thinning of species occurrence records for use in ecological niche models. Ecography 2015, 38, 541–545. [Google Scholar] [CrossRef]
  70. Shaw, J.G.; Sanchez, M.; Spishakoff, L.M.; Trujillo, G.F.; Loppez, D. Dispersal and Migration of Tepa-Sterilized Mexican Fruit Flies. J. Econ. Entomol. 1967, 60, 992–994. [Google Scholar] [CrossRef]
  71. Chambers, D.L.; OConnell, T.B. A Flight Mill for Studies with the Mexican Fruit Fly. Ann. Entomol. Soc. 1969, 62, 917–920. [Google Scholar] [CrossRef]
  72. Mayara, R.; Dos Santos, M.D.; Martins, M.; Fornazier, J.M.; Uramoto, K.; Ferreira, F.; Zucchi, R.A.; Conde, W.A. Aggregation and spatio-temporal dynamics of fruit flies (Diptera, Tephritidae) in papaya orchards associated with different area delimitations in Brazil. Acta Sci. 2020, 44. [Google Scholar] [CrossRef]
  73. Japan Meteorological Agency. El Nino Monitoring and Outlook /TCC; Japan Meteorological Agency: Tokyo, Japan, 2019; Available online: https://www.data.jma.go.jp/multi/index.html?lang=es (accessed on 13 September 2019).
  74. Bureau of Meteorology. Climate Influences Timeline. Australian Government. 2018. Available online: http://www.bom.gov.au/ (accessed on 13 November 2018).
  75. National Weather Service. National Oceanic and Atmospheric Administration; NOAA’s Climate Prediction Center, National Weather Service: Silver Spring, MD, USA, 2018. Available online: https://www.weather.gov/ (accessed on 18 November 2018).
  76. Dupin, J.; Smith, S.D. Integrating historical biogeography and environmental niche evolution to understand the geographic distribution of Datureae. Am. J. Bot. 2019, 106, 667–678. [Google Scholar] [CrossRef]
  77. Wei, R.; Chan, K.W.; So, W.W.M. A systematic review of remote laboratory work in science education with the support of visualizing its structure through the Hist Cite and Cite Space software. Int. J. Sci. Math. Educ. 2017, 15, 1217–1236. [Google Scholar] [CrossRef]
  78. Moo-Llanes, D.A.; Arenas, C.Y.; Baak, B.C. Shifts ecological niche of Lutzomyia peruensis under climate change scenarios in Peru. Med. Vet. Entomol. 2017, 31, 123–131. [Google Scholar] [CrossRef]
  79. NASA, MO DIS Web 2018. Available online: https://modis.gsfc.nasa.gov/ (accessed on 19 December 2018).
  80. NOAA National Oceanic and Atmospheric Administration, What Are El Nino and La Nina? 2020. Available online: https://oceanservice.noaa.gov/facts/ninonina.html (accessed on 20 December 2020).
  81. Acker, J.; Leptoukh, G. Online analysis improves the use of NASA earth science data. Eos Trans. Am. Geophys. Union 2007, 88, 14–17. [Google Scholar] [CrossRef]
  82. Lobo, J.; Jiménez, V.A.; Real, R. AUC: Erratum: Predicting species distribution: Offering more than simple habitat models. Glob. Ecol. Biogeogr. 2008, 17, 145–151. [Google Scholar] [CrossRef]
  83. Carrington, A.M.; Manuel, D.G.; Fieguth, P.W.; Ramsay, T.; Osmani, V.; Wernly, B.; Bennett, C.; Hawken, S.; Magwood, O.; Sheikh, Y.; et al. Deep ROC Analysis and AUC as Balanced Average Accuracy, for Improved Classifier Selection, Audit and Explanation. IEEE Trans. Pattern Anal. Mach. Intell. 2023, 45, 329–341. [Google Scholar] [CrossRef]
  84. Barve, N.; Barve, V.; Jiménez, V.A.; Lira, N.A.; Maher, S.P.; Peterson, A.T.; Soberón, J.; Villalobos, F. The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecol. Modell. 2011, 222, 1810–1819. [Google Scholar] [CrossRef]
  85. Cruz, C.G.; Villaseñor, J.L.; López, M.L.; Martínez, M.E.; Ortiz, E. Selección de predictores ambientales para el modelado de la distribución de especies en Maxent. Rev. Chapingo. Ser. 2014, 20, 187–201. [Google Scholar]
  86. Díaz, A.C.J.; Romero, A.L.V.; Miranda, E.D.R. Neotropical páramos as biogeographic units. Biol. Trop. 2020, 68, 503–516. [Google Scholar] [CrossRef]
  87. Anderson, R.P.; Peterson, T. Evaluating predictive models of species’ distributions: Criteria for selecting optimal models. Ecol. Modell. 2003, 162, 211–232. [Google Scholar] [CrossRef]
  88. Elith, J.; Phillips, S.J.; Hastie, T.; Dudík, M.; Chee, Y.E.; Yates, C.J. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 2011, 17, 43–57. [Google Scholar] [CrossRef]
  89. Mateo, R.; Felcisimo, Á.M.; Muñoz, J. Species distributions models: A synthetic revisión. Rev. Chil. Hist. Nat. 2011, 84, 217–240. [Google Scholar] [CrossRef]
  90. Moo, L.D.A.; López, O.T.; Torres, M.J.A.; Mosso, G.C.; Casas, M.M.; Samy, A.M. Assessing the Potential Distributions of the Invasive Mosquito Vector Aedes albopictus and Its Natural Wolbachia Infections in México. Insects 2021, 12, 143. [Google Scholar] [CrossRef]
  91. Elith, J.; Leathwick, J.R. Species Distribution Models: Ecological Explanation and Prediction Across Space and Time. Annu. Rev. Ecol. Evol. Syst. 2009, 40, 677–697. [Google Scholar] [CrossRef]
  92. Cobos, M.E.; Peterson, A.T.; Osorio, O.L.; Jiménez, G.D. An exhaustive analysis of heuristic methods for variable selection in ecological niche modeling and species distribution modeling. Ecol. Inform. 2019, 53. [Google Scholar] [CrossRef]
  93. Pliscoff, P.; Fuentes, T. Modelación de la distribución de especies y ecosistemas en el tiempo y en el espacio: Una revisión de las nuevas herramientas y enfoques disponibles. Rev. Geogr. Norte Gd. 2011, 48, 61–79. [Google Scholar] [CrossRef]
  94. Owens, H.L.; Campbell, L.P.; Dornak, L.L.; Saupe, E.E.; Barve, N.; Soberón, J.; Ingenloff, K.; Lira, N.A.; Hensz, C.M.; Myers, C.E.; et al. Constraints on interpretation of ecological niche models by limited environmental ranges on calibration areas. Ecol. Modell. 2013, 263, 10–18. [Google Scholar] [CrossRef]
  95. Brown, J.L. SDM toolbox: A python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods Ecol. Evol. 2014, 5, 694–700. [Google Scholar] [CrossRef]
  96. Brown, J.L.; Bennett, J.R.; Connor, M. SDMtoolbox 2.0: The next generation Python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. PeerJ 2017, 5, e4095. [Google Scholar] [CrossRef]
  97. Pearson, R.G.; Raxworthy, C.J.; Nakamura, M.; Peterson, A.T. Predicting species distributions from small numbers of occurrence records: A test case using cryptic geckos in Madagascar. J. Biogeogr. 2006, 34, 102–117. [Google Scholar] [CrossRef]
  98. Pinto, J.N.; Bares, C.J. Predicting species distributions and community composition using satellite remote sensing predictors. Sci. Rep. 2021, 11, 16448. [Google Scholar] [CrossRef] [PubMed]
  99. Paini, D.R.; Sheppard, A.W.; Cook, D.C.; De Barro, P.J.; Worner, S.P.; Thomas, M.B. Global threat to agriculture from invasive species. Proc. Natl. Acad. Sci. USA 2016, 113, 7575–7579. [Google Scholar] [CrossRef]
  100. Lódola, A. Contratistas, Cambios Tecnológicos y Organizacionales en el Agro Argentino, Documentos de Proyectos No. 176. Available online: https://hdl.handle.net/11362/36772 (accessed on 25 April 2023).
  101. Atencio, V.R.; Collantes, G.R.; Caballero, E.M.; Hernández, A.P.; Vaña, H.M. Impacto de los Insectos en la Seguridad Alimentaria en Panamá. Sci. Agropecu 2021, 36, 139–165. [Google Scholar]
  102. Amat, E.; Altamiranda, S.M.; Canal, N.; Gómez, P.L. Changes in the potential distribution of the guava fruit fly Anastrepha striata (Diptera, Tephritidae) under current and possible future climate scenarios in Colombia. Bull. Entomol. Res. 2022, 112, 469–480. [Google Scholar] [CrossRef]
  103. Machado, T.C.; Krüger, P.; Edson, N.D.; Mello, G.F. Potential global distribution of the south American cucurbit fruit fly Anastrepha grandis (Diptera: Tephritidae). Crop Prot. 2021, 45, 0261–2194. [Google Scholar] [CrossRef]
  104. Silva, J.G. Biologia e Comportamento de Anastrepha Grandis (Macquart, 1846) (Diptera: Tephritidae); Dissertação Mestrado. Universidade de São Paulo: São Paulo, Brazil, 1991; Available online: https://repositorio.usp.br/item/000733447. (accessed on 1 March 2023).
  105. Topón, R.L.M. Ciclo Biológico de la Mosca de la Fruta del Género (anastrepha spp.) a dos temperaturas, Universidad Tecníca de Cotopaxi. Latacunga Tesis en Ingieneria Agronomica, Av. Simón Rodríguez, Latacunga, Ecuador. 2020, p. 93. Available online: http://repositorio.utc.edu.ec/handle/27000/7050. (accessed on 15 April 2023).
  106. Fúnez, X.I. Características Biológicas de Anastrepha Grandis (Macuart,1846) en Relación con su Hospedero Natural, Fevillea Cordifolia en Darién, Panamá. Maestría Thesis, Universidad de Panamá, Ciudad de Panamá, Panamá, 2014. Available online: http://up-rid.up.ac.pa/id/eprint/348. (accessed on 15 April 2023).
  107. Bonebrake, T.C.; Deutsch, C.A. Climate heterogeneity modulates impact of warming on tropical insects. Ecology 2012, 93, 449–455. [Google Scholar] [CrossRef]
  108. De La Vega, G.J.; Schilman, P.E. La importancia de la fisiología en la distribución geográfica de los insectos. Rev. Soc. Entomol. Arg. 2015, 74, 101–108. [Google Scholar]
  109. Addo, B.A.; Chown, S.L.; Gaston, K.J. Thermal tolerance, climatic variability and latitude. Proc. R. Soc. Lond. Ser. B Biol. Sci. 2000, 267, 739–745. [Google Scholar] [CrossRef]
  110. Esparza, M. La sequía y la escasez de agua en México: Situación actual y perspectivas futuras. Secuencia 2014, 89, 193–219. [Google Scholar] [CrossRef]
  111. González, O.D.; Córdoba, A.A.; Dáttilo, W.; Noriega, L.A.; Sánchez, G.R.A.; Villalobos, F. Insect responses to heat: Physiological mechanisms, evolution and ecological implications in a warming world. Biol. Rev. 2020, 95, 802–821. [Google Scholar] [CrossRef] [PubMed]
  112. Bruniard, E.D. La diagonal árida Argentina: Un límite climático real. Rev. Geográfica 1982, 95, 5–20. [Google Scholar]
  113. Mason, S.J.; Goddard, L. Probabilistic Precipitation Anomalies Associated with EN SO. Research Paper. Bull. Am. Meteorol. Soc. 2001, 82, 619–638. [Google Scholar] [CrossRef]
  114. Karlin, M.S. Cambios temporales del clima en la subregión del Chaco Árido. Multequina 2012, 21, 3–16. [Google Scholar]
  115. Fava, G.A.; Acosta, J.C.; Blanco, G.M. The effects of seasonality and precipitation in the avifauna of the Argentine Southern Chaco Serrano. Rev. Biol. Trop. 2017, 65, 953–961. [Google Scholar] [CrossRef]
  116. Flores, A.R.; Kazuz, Y.E.M.; García, V.E.; Ayala, B.A.; Garrido, R.E.R.; Aceves, M.A.C.; Sánchez, O.M.A.; Tejacal, A.I. Control de Anastrepha serpentina (Wiedemann) y calidad de los frutos de zapote mamey Pouteria sapota (Jacq) Moore & Stearn tratados con aire caliente forzado. Rev. Chapingo. Ser. Cienc. 2009, 15, 9–15. [Google Scholar]
  117. Pinson, E.P.; Tejada, L.O.; Toledo, J.; Enkerlin, W.; Hurtado, C.H.; Valle, J.; Pérez, J.N.; Liedo, P. Caracterización de la adaptación de anastrepha serpentina (wied.) (diptera: Tephritidae) a condiciones de cría masiva. Fol. Entomol. Mex. 2006, 45, 97–112. [Google Scholar]
  118. Alexander, M.; Kulminski, F.M.; Irina, V.; Culminskaya, K.G.; Arbeev, S.V.; Ukraintseva, J.R.; Carey, A.I. Date of eclosion modulates longevity: Insights across dietary-restriction gradients and female reproduction in the mexfly Anastrepha ludens. Exp. Gerontol. 2009, 44, 718–726. [Google Scholar] [CrossRef]
  119. Alonso, E.H. Sincronía Biológica, Relación Interespecífica y Análisis de Calidad Hospedera de Pouteria Buenaventurensis (Sapotacea) con Anastrepha Serpentina y Anastrepha Intermedia, n.sp. en Altos de Pacora. Maestría Thesis, Universidad de Panamá, Ciudad de Panamá, Panamá, 2000. Available online: http://up-rid.up.ac.pa/id/eprint/3949 (accessed on 15 April 2023).
  120. Omkar, R.S.; Pandey, P. Effect of temperature on development and immature survival of Zygogramma bicolorata (Coleoptera: Chrysomelidae) under laboratory conditions. Int. J. Trop. Insect Sci. 2008, 28, 130–135. [Google Scholar] [CrossRef]
  121. Chaverri, L.G.; Soto, M.J.; Jirón, L.F. Biology and ecology of Anastrepha obliqua (diptera: Tephritidae), plague of Anacardiaceae plants in tropical America. II. Mature stages. Agron. Mesoam. 2006, 10, 99–102. [Google Scholar] [CrossRef]
  122. Soto, M.J.; Chaverri, L.G.; Jirón, L.F. Notes on the biology and ecology of Anastrepha obliqua (Diptera: Tephritidae), pests of plants in Tropical América. I. Imnature forms. Agron. Mesoam. 2016, 8, 116–120. [Google Scholar] [CrossRef]
  123. Hernández, E.; Ruiz, M.L.; Toledo, J.; Montoya, P.; Liedo, P.; Aceituno, M.M.; Perales, H. A comparison of sexual competitiveness and demographic traits of Anastrepha obliqua (Macquart) (Diptera: Tephritidae) among fruit-associated populations. Bull. Entomol. Res. 2019, 109, 333–341. [Google Scholar] [CrossRef]
  124. Candanedo, M.; Villarreal, D.; Bernal, S. Uso de registros de temperatura máxima promedio de las estaciones meteorológicas de ETESA, para la creación de mapas de temperatura mediante el uso de programa ArcGIS. Rev. De Iniciación Científica 2020, 6, 9–14. [Google Scholar] [CrossRef]
  125. Kemp, W.P.; Bosch, J. Effect of Temperature on Osmia lignaria (Hymenoptera: Megachilidae) Prepupa–Adult Development, Survival, and Emergence. J. Econ. Entomol. 2005, 98, 1917–1923. [Google Scholar] [CrossRef]
  126. Bateman, M.A. The Ecology of Fruit Flies. Annu. Rev. Entomo. 1972, 17, 493–518. [Google Scholar] [CrossRef]
  127. Cruz, L.L.; Malo, E.A.; Rojas, J.C. Sex Pheromone of Anastrepha striata. J. Chem. Ecol. 2015, 41, 458–464. [Google Scholar] [CrossRef] [PubMed]
  128. Péres, D.; Aluja, M. Anastrepha striata (Diptera: Tephritidae) Females That Mate with Virgin Males Live Longer. Ann. Entomol. Soc. Am. 2004, 97, 1336–1341. [Google Scholar] [CrossRef]
  129. García, A.E.C.; Martínez, A.J.O.; Gómez, P.L.M. Distribución Geográfica Potencial de Anastrepha Striata (Schiner 1868)(Diptera: Tephritidae) en Colombia; Documento de conferencia, III Congreso Colombiano de Zoología: Medellín, Colombia, 2010; Available online: https://dspace.tdea.edu.co/handle/tdea/1462 (accessed on 12 May 2023).
  130. Ferrer, S.Y.; Jacho, S.A.A.; Wilmer, R.; Zambrano, U.; Absalo, P.J.A.; Vázquez, P.F.; Herminio, A.; Zambrano, M.G.J.; Castillo, M.M.J.; Rosado, C.A.; et al. Invasiones Biológicas en Agroecosistemas de Ecuador Continental: Nicho Ecológico de Especies Exóticas y Cultivos Agrícolas Bajo Riesgo. Acta Biol. Colomb. 2021, 26, 352–364. [Google Scholar] [CrossRef]
Figure 1. Study area map of potential distribution models.
Figure 1. Study area map of potential distribution models.
Insects 14 00714 g001
Figure 2. Potential distribution of Anastrepha grandis in the neotropics for El Niño, Neutral, and La Niña events ((AC), respectively), where high suitability values (suitable areas) are represented in gray and low suitability (unsuitable areas) are in green. In addition, various maps are shown analyzing the MOP in the different El Niño, Neutral, and La Niña events ((DF), respectively), where high environmental similarity values are represented in blue and low environmental similarity in brown.
Figure 2. Potential distribution of Anastrepha grandis in the neotropics for El Niño, Neutral, and La Niña events ((AC), respectively), where high suitability values (suitable areas) are represented in gray and low suitability (unsuitable areas) are in green. In addition, various maps are shown analyzing the MOP in the different El Niño, Neutral, and La Niña events ((DF), respectively), where high environmental similarity values are represented in blue and low environmental similarity in brown.
Insects 14 00714 g002
Figure 3. Potential distribution of Anastrepha serpentina in the neotropics for El Niño, Neutral, and La Niña events ((AC), respectively), where high suitability values (suitable areas) are represented in gray and low suitability (unsuitable areas) are in orange. In addition, various maps are shown analyzing the MOP in the different El Niño, Neutral, and La Niña events ((DF), respectively), where high environmental similarity values are represented in blue and low environmental similarity in brown.
Figure 3. Potential distribution of Anastrepha serpentina in the neotropics for El Niño, Neutral, and La Niña events ((AC), respectively), where high suitability values (suitable areas) are represented in gray and low suitability (unsuitable areas) are in orange. In addition, various maps are shown analyzing the MOP in the different El Niño, Neutral, and La Niña events ((DF), respectively), where high environmental similarity values are represented in blue and low environmental similarity in brown.
Insects 14 00714 g003
Figure 4. Potential distribution of Anastrepha obliqua in the neotropics for El Niño, Neutral, and La Niña events ((AC), respectively), where high suitability values (suitable areas) are represented in gray and low suitability (unsuitable areas) are in yellow. In addition, various maps are shown analyzing the MOP in the different El Niño, Neutral, and La Niña events ((DF), respectively), where high environmental similarity values are represented in blue and low environmental similarity in brown.
Figure 4. Potential distribution of Anastrepha obliqua in the neotropics for El Niño, Neutral, and La Niña events ((AC), respectively), where high suitability values (suitable areas) are represented in gray and low suitability (unsuitable areas) are in yellow. In addition, various maps are shown analyzing the MOP in the different El Niño, Neutral, and La Niña events ((DF), respectively), where high environmental similarity values are represented in blue and low environmental similarity in brown.
Insects 14 00714 g004
Figure 5. Potential distribution of Anastrepha striata in the neotropics for El Niño, Neutral, and La Niña events ((AC), respectively), where high suitability values (suitable areas) are represented in gray and low suitability (unsuitable areas) are in green. In addition, various maps are shown analyzing the MOP in the different El Niño, Neutral, and La Niña events ((DF), respectively), where high environmental similarity values are represented in blue and low environmental similarity in brown.
Figure 5. Potential distribution of Anastrepha striata in the neotropics for El Niño, Neutral, and La Niña events ((AC), respectively), where high suitability values (suitable areas) are represented in gray and low suitability (unsuitable areas) are in green. In addition, various maps are shown analyzing the MOP in the different El Niño, Neutral, and La Niña events ((DF), respectively), where high environmental similarity values are represented in blue and low environmental similarity in brown.
Insects 14 00714 g005
Figure 6. Potential distribution of Anastrepha spp. in Panama for the El Niño, Neutral, and La Niña events: A. grandis (A), A. serpentina (B), A. obliqua (C), and A. striata (D). The suitable areas are shown in green, and the unsuitable areas are represented in gray.
Figure 6. Potential distribution of Anastrepha spp. in Panama for the El Niño, Neutral, and La Niña events: A. grandis (A), A. serpentina (B), A. obliqua (C), and A. striata (D). The suitable areas are shown in green, and the unsuitable areas are represented in gray.
Insects 14 00714 g006
Figure 7. Analysis of changes in the ranges of expansion, no occupation, no change, and contraction range in the different scenarios for each species ((A)–A. grandis, (B)–A. serpentina, (C)–A. obliqua and (D)–A. striata) in the comparison of Neutral vs. El Niño and Neutral vs. La Niña events.
Figure 7. Analysis of changes in the ranges of expansion, no occupation, no change, and contraction range in the different scenarios for each species ((A)–A. grandis, (B)–A. serpentina, (C)–A. obliqua and (D)–A. striata) in the comparison of Neutral vs. El Niño and Neutral vs. La Niña events.
Insects 14 00714 g007
Table 1. Models selected based on defined criteria with partial ROC statistical significance, skip rate performance, and AICc complexity.
Table 1. Models selected based on defined criteria with partial ROC statistical significance, skip rate performance, and AICc complexity.
SpeciesRecordsTrainSelected ModelPartial ROCOmission Rate (<5%)AICc∆AICc
Anastrepha grandis3932M_0.9_F_t_Set700392.8820.000
Anastrepha serpentina8870M_0.4_F_lq_Set8001017.8490.000
Anastrepha obliqua9374M_4_F_l_Set7001080.4750.000
Anastrepha striata150125M_3_F_lt_Set7001854.4800.000
Table 2. Suitability percentages in the neotropics for each species of Anastrepha in the comparison of Neutral vs. El Niño and Neutral vs. La Niña events.
Table 2. Suitability percentages in the neotropics for each species of Anastrepha in the comparison of Neutral vs. El Niño and Neutral vs. La Niña events.
A. grandis Neutral vs. NiñoA. grandis Neutral vs. NiñaA. serpentina Neutral vs. NiñoA. serpentina Neutral vs. NiñaA. obliqua Neutral vs. NiñoA. obliqua Neutral vs. NiñaA. striata Neutral vs. NiñoA. striata Neutral vs. Niña
Range expansion14.7211.741.171.981.993.7218.9122.17
No occupancy9.4112.3810.149.339.778.049.826.56
No change61.5359.3874.5774.2986.486.5966.8468.44
Range contraction14.3316.4714.1114.381.771.634.412.89
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Degracia, A.B.; Jiménez, J.Á.; Alvarado, A.B.; Valdespino, R.A.; Altamiranda-Saavedra, M. Evaluation of the Effect of the ENSO Cycle on the Distribution Potential of the Genus Anastrepha of Horticultural Importance in the Neotropics and Panama. Insects 2023, 14, 714. https://doi.org/10.3390/insects14080714

AMA Style

Degracia AB, Jiménez JÁ, Alvarado AB, Valdespino RA, Altamiranda-Saavedra M. Evaluation of the Effect of the ENSO Cycle on the Distribution Potential of the Genus Anastrepha of Horticultural Importance in the Neotropics and Panama. Insects. 2023; 14(8):714. https://doi.org/10.3390/insects14080714

Chicago/Turabian Style

Degracia, Arturo Batista, Julián Ávila Jiménez, Anovel Barba Alvarado, Randy Atencio Valdespino, and Mariano Altamiranda-Saavedra. 2023. "Evaluation of the Effect of the ENSO Cycle on the Distribution Potential of the Genus Anastrepha of Horticultural Importance in the Neotropics and Panama" Insects 14, no. 8: 714. https://doi.org/10.3390/insects14080714

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop