Revista Mexicana Ciencias Agrícolas   volume 10  number 4   May 16 - June 29, 2019

DOI: https://doi.org/10.29312/remexca.v10i4.1661

Article

Production potential of Vanilla planifolia Jacks in Totonacapan,
Mexico, using geographic techniques

Alberto Santillán-Fernández1

Maryjose Trejo Cabrera1

Araceli Martínez Sánchez1

Luis Martínez Ángel1

Nehemías Vásquez Bautista2

Salomón Luis Mejía1§

1Forestry Engineering-Higher Technological Institute of Venustiano Carranza. Av. Tecnológico s/n, Col. El Huasteco, Lázaro Cárdenas, Venustiano Carranza, Puebla, Mexico. CP. 73049. (asantillanq@itsvc.edu.mx; marytreca@gmail.com; chelyforestal@itsvc.edu.mx; luismartinezangel99@gmail.com;). 2Trusts Instituted in Relation to Agriculture. Av. Tulum 238 and 240, lot 81 and 82, super block 4 and block 12, Cancún, Quintana Roo, Mexico. CP. 77500. (nvasquez@fira.gob.mx).

§Corresponding author: forestal@itsvc.edu.mx.

Abstract

The Totonacapan region, between the northern limits of Puebla and Veracruz, produces 80% of the vanilla in Mexico, and is considered the center of origin of the species Vanilla planifolia Jacks. However, the area sown has decreased, which affects the competitiveness of the sector. In the present work a socio-economic characterization of the producers of Vanilla planifolia Jacks was carried out during 2017 and 2018 and the productive potential for the cultivation in the Totonacapan region was determined through the use of geographic information systems (GIS) based on factors edaphoclimatic and physiographic. Results show that production is a secondary economic activity, among producers over 55 years of age and in a smallholder situation. The highest productive potential is located northwest of the Totonacapan region and is conditioned by the increase in rainfall in the driest month (April: <60 mm year-1) and temperatures no higher than 30 °C in the warmest quarter (June-July-August). These areas are mostly destined to livestock and agriculture of more profitable cyclical crops, so in the expansion of the crop should consider the socioeconomic aspects of producers.

Keywords: maxent, maximum entropy, multicriteria evaluation, potential distribution.

Reception date: March 2019

Acceptance date: June 2019

Introduction

The cultivation of vanilla (Vanilla planifolia Jacks) had its origin in Mexico, particularly in the northern area of Veracruz, where the Totonac culture settled (Lubinsky et al., 2008), which used and domesticated vanilla for the first time (Sinha et al., 2008). In the Totonacapan region, the cultivation of vanilla has been determinant for the socioeconomic development of producers (Herrera-Cabrera et al., 2012). The vanilla produced there is considered the best quality and largest aromatic compounds worldwide, these characteristics are attributable to the geographical environment in which it develops, natural factors and crop management practices (Lubinsky et al., 2008; Flores et al., 2017).

Despite being the center of origin, as a result of a reduction in the planted area, historically Mexico went from being the largest center of vanilla production at the international level, to occupy the fifth place, with 4.97% of the annual net volume, behind from Indonesia (34.93%), Madagascar (31.81%), China (11.63%) and Papua New Guinea (6.97%) (FAOSTAT, 2018). Particularly in Mexico, Veracruz is the main producer with 70% of the volume of production, in order of importance continue Puebla and Oaxaca, which together contribute 29% of total production and in smaller quantities, also occurs in San Luis Potosí, Hidalgo, Chiapas and Quintana Roo.

The most important production area in Mexico is located north of Puebla and in north-central Veracruz and is known as the Totonacapan region with more than 80% of the national production (SIAP, 2018), where germplasm is also located of greater biological and commercial importance (Herrera-Cabrera et al., 2012). However, the regional and, in general, the national sector, faces low productivity, attributed to the reduction of the sown area and the abandonment of the crop, as a consequence of a fragile organization of producers and problems associated with crop management. The good field yields presented by Mexico (0.41 t ha-1) are insufficient to reverse this situation (Barrera-Rodríguez et al., 2011; Jaramillo et al., 2012; Barrera et al., 2014; Santillán-Fernández et al., 2018).

Given this scenario, the vanilla production system in the Totonacapan region seems destined to disappear; however, according to Lubinsky et al. (2008); Barrera-Rodríguez et al. (2011) the production of vanilla in the region has an ancestral connotation that can revalue the product as a center of origin and have the slogan of producing the best quality vanilla in the world. For this purpose, INIFAP (2011); SAGARPA (2012); Santillán-Fernández et al. (2018) propose that at the same time as solving the socio-economic and crop management problems associated with production, areas with productive potential should be identified that consider edaphoclimatic aspects, to increase the sown areas, which is the problem of international competitiveness of the national vanilla producer sector.

In this sense, geographic information system (GIS) techniques such as multicriteria evaluation (MCE) and species distribution models are useful tools to characterize the agroecological potential of the area destined for a crop, even more so if in the area of there is no similar background in the study (Jiménez et al., 2004; Kazuya et al., 2006; Aguilar et al., 2010). Within the species distribution models, the Maximum Entropy algorithm (MaxEnt) has the best evaluations in comparison with other programs that model the probability of occurrence and adaptation of a species from edaphoclimatic variables (Hernández et al., 2006; Phillips et al., 2006; Elith et al., 2011; Navarro-Cerrillo et al., 2011).

This tool has been applied successfully in the determination of the productive potential of potato crops (Solanum tuberosum) (Hijmans et al., 2000), papaya (Carica papaya) (Scheldeman et al., 2007) and vanilla (Vanilla planifolia) (Hernández-Ruiz et al., 2016; Flores et al., 2017). Phillips et al. (2006) and Elith et al. (2011) mention that many authors prefer the MaxEnt algorithm due to parsimony and its proper adjustment to the distribution of species when only presence information is used.

For its part, the MCE associated with GIS is an exact system when information of an environmental, social and economic nature is combined, which based on the weighting and compensation of determining variables or fitness factors, generates cartography that expresses the productive potential of a region with respect to the specific requirements of a crop (Gómez and Barredo, 2005) and has been applied in the determination of suitable areas for the cultivation of Agave durangensis (Olivas et al., 2007), of sugar cane (Saccharum officinarum) (Aguilar et al., 2010) and alternative crops such as amaranth (Amaranthus spp.) and nopal (Opuntia ficus-indica) (Ceballos-Silva and López-Blanco, 2010).

The cultivation of vanilla relates technical, social, ecological and climatic factors that condition its development, yield, conservation and sown areas, by virtue of which it is necessary to know its current state of production and characterize potential areas with higher yields in the field; knowledge about these factors that influence the management of the crop are determining factors in its expansion (Flores et al., 2017; Rocha et al., 2018).

In this context, the objectives of this study were to a) characterize socioeconomically the vanilla producers in the Totonacapan region; and b) determine the areas with the greatest productive potential to establish the vanilla crop through MaxEnt and MCE models. This information will allow to make reliable decisions for the establishment of this crop in areas with greater potential, that have the best edaphoclimatic conditions for its development and thus improve yields in the field, which help to improve the international competitiveness of the sector.

Materials and methods

Study area

The Totonacapan region is the most important production area in Mexico and is located north of Puebla (it includes 61 municipalities) and in north-central Veracruz (20 municipalities) (Figure 1). Due to its physiographic and edaphoclimatic characteristics, it is an optimal region for the development of vanilla cultivation, since subtropical, warm and humid climates are located in altitudes below 600 m, with temperatures ranging from 20 °C to 32 °C and average annual rainfall ranging from 1 500 to 3 000 mm, in addition to having land with excellent drainage, organic matter greater than 3.5%, with a pH of 5.5 to 7 and slopes of 10% to 40% (INIFAP, 2011; SAGARPA, 2012).

The Totonacapan region collects the waters of the Panuco, Tuxpan, Cazones and Tecolutla rivers, which allows the development of livestock, forestry and agriculture, with crops of corn (Zea mays), beans (Phaseolus vulgaris), squash (Cucurbita spp.) and pepper (Capsicum annuum) as main for self-consumption. In addition to commercial species such as coffee (Coffea spp.), sugar cane (Saccharum officinarum), tobacco (Nicotiana tabacum), pepper (Pimenta dioica), banana (Musa paradisiaca), sesame (Sesamum indicum), ginger (Zingiber officinale), litchi (Litchi chinensis), citrus (Citrus sinensis) and vanilla (SIAP, 2018).

Figure 1. Spatial location of the Totonacapan region composed of municipalities of Puebla and Veracruz.

Socioeconomic characterization of vanilla production in the Totonacapan region

To determine the socioeconomic factors that explain the changes in area planted and yield in the field in Totonacapan, a survey was conducted between November 2017 and February 2018 to 175 vanilla producers (20% of the total, statistically sufficient sample according to the theory of sampling) based on their availability and references in the region. In addition, specialists from the Autonomous University of Chapingo (2), the Postgraduate School (1) and representatives of the Product System Committee of Puebla (2) were interviewed. The interview guide considered phytosanitary, climatic, social and market aspects.

The survey included questions about age, schooling, cultivated area, land tenure and percentage that represents the activity in the family income. According to the CIMMYT survey design guide (1993), aspects regarding the knowledge that the producer has about irrigation, fertilization, pest and disease control, as well as other factors that could condition the expansion of the crop and increase its productivity, such as the climate, the physiography of the plot, production costs, management of subsidies and credits, management and parcel equipment, among others.

Vanilla productive potential using the maximum entropy algorithm (MaxEnt)

Through in situ tours in 2016 and with the help of the directors of the Association of Vanilla producers of the State of Puebla, a total of 167 points with presence of vanilla in the field of the municipalities with production were georeferenced for their accessibility; this information was combined with climatic and edaphic variables of the area, with a resolution of 30 s (1 km2), the first correspond to the 19 bioclimatic variables and annual evaporation, and the second are: calcium, electrical conductivity, organic carbon, potassium , magnesium, organic matter, sodium, pH and sodium adsorption ratio, available in the global database of climatic surfaces WorldClim.

The points with presence of vanilla in field (167), were introduced in the MaxEnt algorithm (MaxEnt, 2016) in comma delimited format (.csv) and the climatic and edaphic variables in ASCII format (.asc), following the methodology de Hijmans et al. (2005) was determined the productive potential of the region for the cultivation of vanilla. The results were imported into DivaGis v7.5 (DivaGis, 2016) and converted into grig format (.grd) to be viewed as images. The probability of occurrence of the species was categorized into areas with potential i) very low [0.0 - 0.2]; ii) low (0.2-0.4); iii) medium (0.4-0.6); iv) high (0.6-0.8); and v) very high (0.8-1), according to a Likert scale of homogeneous intervals.

Vanilla productive potential using multicriteria evaluation

From the geoportal of the National Commission for the Knowledge and Use of Biodiversity (CONABIO), the georeferenced data were extracted in vector format 1: 1 000 000 in altitude, precipitation, temperature, pH and soils and in raster format with a resolution of 5 min the digital elevation model that served to delimit the slopes of the region. The climatic and edaphic ranges were used where the highest field yield is obtained: altitude (less than 600 m), precipitation (1 500 to 3 000 mm), temperature (20 to 32 °C), type of soil (with organic matter) greater than 3.5% and pH higher than 5.5) and pending (10 to 40%) (INIFAP, 2011; SAGARPA, 2012), to discriminate those areas that do not meet these specific needs.

For the procedure of weighting and comparison of said aptitude factors, the Weighted Overlay module of ArcMap® (ArcMap, 2010) was used, the evaluation scale (very low, low, medium, high and very high) was determined by routes in field and according to the opinion of specialists in the production of vanilla of the Autonomous University Chapingo (UACH), College of Postgraduates (CP) and representatives of the Committee System Product of Puebla AC.

Discrimination of the productive potential by land use and current vegetation

Finally, to delimit those regions suitable for the expansion of the crop, areas were determined by productive potential (very low, medium, high and very high) based on the use of soil and vegetation reported by the INEGI (2016) scale 1: 250 000 of the VI series. Through field trips with directors of the Association of Vanilla Producers of the State of Puebla AC, the information obtained by the two methods (MaxEnt and MCE) was verified, which served as the basis for the proposal to establish new areas for the cultivation of Vanilla planifolia with greater productive potential.

Results and discussion

Socioeconomic characterization of vanilla producers in the Totonacapan region

The Product System Committee of Puebla AC indicates that in the region of Totonacapan there are a total of 871 producers who cultivate 2 653 ha of vanilla, of which 80% are men and 20% women, belong mostly to the Otomi ethnic groups and totonaco. Of the total of producers, 487 have this crop in land of small property with a planted surface of 1 486 ha, and 384 producers develop the activity in lands of common property with a total of 1 167 ha.

95% of the producers surveyed produce vanilla as a secondary activity, and have as main sources of income, the production of corn, beans, coffee, pepper and livestock. In most cases the production of vanilla is practiced more by a cultural factor than by the profitability of the crop, since this activity is inherited from generation to generation and the pods are used for the manufacture of crafts that are occupied in the festivities regional

Experts from the UACH, CP and representatives of the Product System Committee of Puebla AC, consider the smallholding (<0.5 ha in 79% of producers), the low level of study (<8 years), the age (from 55 to 62) years) and the fact that the production of vanilla means a secondary activity that complements the family income (Table 1), impedes the development of the vanilla sector in the region, by limiting the area sown and in many cases converting it to more profitable crops.

Table 1. Socioeconomic characterization of the surveyed producers of vanilla in the Totonacapan region according to the sown area and net production.

IF= family income that represents the activity.

Alcon et al. (2008) point out that older producers are more reluctant to change their production systems, considering that they have always done so and for them it is well done. On the other hand, Soto et al. (2006) establish that it is more complicated for producers with small plots to strengthen their production systems with new infrastructure, since in many cases they represent a secondary activity in their family income. In this regard, Álvarez and Crecente (2000) describe the benefits of land consolidation (through rights-of-use session) as a measure of adaptation to smallholdings and crop spraying.

Factors of crop management and access to credit

Added to the problems of smallholding, the vanilla producer of the Totonacapan region faces a scarce development of technological packages for the management of the crop, which reduce production by not considering programs for the control of pests, diseases, fertilization and irrigation. 98% of the planted area is managed under a rainfed water regime, only 2% uses irrigation. 90% of the producers surveyed mentioned that they have noticed variations in rainfall patterns, so they are aware that in the future watering may be required to maintain the activity. However, they do not consider the climatic factors, nor the physiography of the plots as limiting aspects for production.

95% of the producers surveyed consider that the expansion of the crop is limited by phytosanitary aspects related to the control of pests and diseases such as red bug (T. confusus), hairy worm (P. aunifera), root and stem rot (F. oxysporum) and anthracnose (C. gloeosporioides). In addition to the zero support for access to credit and subsidies that facilitate the parcel equipment and incorporate technological packages (especially fertilization) that improve productivity in the region. Practically 100% of producer’s self-finance production costs, which in an unstable market in prices causes many of them to abandon the activity by converting their plots to more profitable crops.

INIFAP (2011); SAGARPA (2012), consider that the adequate management of the problems of fertilization, pests and diseases associated with the production of vanilla in the region, substantially improve the productivity of the sector. By virtue of this, as of 2011 they implemented technological packages that included these aspects, with good results in their beginning when increasing the yields in the field. However, due to lack of investment, continuity has been scarce. To solve this problem, the Product System Committee of Puebla AC, currently manages trusts instituted in relation to agriculture (FIRA) credits to the producer to strengthen the production infrastructure.

The market factor: lack of added value

The average yields in field of the region are of 0.449 t ha-1 with production costs of $60 000.00 year ha-1 and with gains superior to the $88 000.00 year ha-1, these gains correspond to the sale of the pod in green whose market value ranges from $180 to $200 kg. In the region 94% of the producers do not give the beneficiary service where the selling costs are over $3 000 kg. This aspect has been addressed by Jaramillo et al. (2012), Jaramillo et al. (2013) and Xochipa-Morante et al. (2016) who consider that this factor is what limits the expansion of the crop, since the producer considers the activity to be unprofitable and prefers to channel its inputs to economically more attractive crops.

Barrera et al. (2014) mention that the international vanilla market imposes new and strict quality and safety standards on producers that they must abide by in order to be competitive; so the international market is a factor that directly affects the productivity of vanilla. By virtue of this, the Product System Committee of Puebla AC, implemented a series of workshops to promote among the producers of the region the vanilla beneficiation service, with which it is expected that they will substantially improve their income and the competitiveness of the vanilla sector.

Vanilla productive potential using the maximum entropy algorithm (MaxEnt)

The results obtained by the MaxEnt algorithm are highly significant with a reliability level of 95% (p-value <0.0001) and a coefficient of determination (R2) of 99.6%. Of the 29 edaphoclimatic variables used in the MaxEnt probabilistic model, four contributed with 88.1% of the prediction of zones with the highest productive potential for the vanilla crop. The most important were: increase in rainfall in the driest month (April: <60 mm year-1) that contributed 56.3% of the prediction, followed by a low content of sodium in the soil with 13.1%, temperatures no higher than 30 °C in the hottest quarter (June-July-August) with 10.7% and high organic carbon content in the soil with 8%.

Therefore, according to the model, the factor that most affects the optimal development of the vanilla crop is the increase in precipitation in the driest month. This aspect was documented by Barrera-Rodríguez et al. (2009); INIFAP (2011); SAGARPA (2012); Flores et al. (2017); Santillán-Fernández et al. (2018) who agree that the adoption of irrigation during periods of low water substantially improves yields in the field, even without fertilization. On the other hand, the surveyed producers mentioned that they have perceived variations in rain patterns, so they do not rule out irrigation in the future.

Figure 2 shows the productive potential of the region classified in: i) very low; ii) low; iii) medium; iv) high; and v) very high. The current production areas and the areas planted during the year 2017 are also located, which were proposed for this crop before the Association of Vanilla Producers of the State of Puebla AC, and which were selected based on the results obtained and availability in the field. , when the producers voluntarily adopt the vanilla crop.

Figure 2. Spatial representation of the estimation of the productive potential for vanilla in the Totonacapan region obtained by means of the maximum entropy algorithm.

For the MCE an evaluation scale of five classes was built, with areas of potential i) very low; ii) low; iii) medium; iv) high; and v) very high (Table 2). The MCE turned out to be a more robust methodology, which considers larger areas compared to MaxEnt.

Table 2. Criteria for evaluating climatic and edaphic variables to establish the productive potential of vanilla in the Totonacapan region, built from field trips and expert opinion.

*= FAO classification system (1999).

Figure 3 shows those areas in the Totonacapan region that meet the agroclimatic conditions of the crop, specific areas with medium, high and very high aptitude for the establishment of the vanilla, where to improve the management of the crop with aspects such as control phytosanitary, irrigation and fertilization, it is expected to obtain the best yields in the field of the region. The results were taken as a reference by the Association of Vanilla Producers of the State of Puebla AC, to manage the sowing of 100 ha with vanilla between May and June 2017.

When considering edaphoclimatic factors (annual precipitation, annual temperature, pH, type of soil and organic matter) and physiographic factors (altitude and slope) in the MCE, a greater area with medium to very high productive potential was delimited with respect to the results obtained with the MaxEnt algorithm, which did not include physiographic variables. However, both methods agree that the greatest productive potential is located northwest of the Totonacapan region, in the border areas of Puebla and Veracruz, which is where the current producing areas are located, with medium to very high productive potential.

Figure 3. Spatial representation of the estimation of the productive potential for vanilla in the Totonacapan region obtained through multicriteria evaluation.

Discrimination of the productive potential by land use and current vegetation

By using land use and vegetation cartography (INEGI, 2016), it was possible to identify the current land use of the areas with productive potential for the cultivation of vanilla obtained by the MaxEnt and MCE algorithms (Table 3). Coverages where the expansion of the crop can be potentiated are located in those areas where the continuity of an ecosystem is not put at risk; therefore, agriculture and grassland areas are the most optimal, as long as socioeconomic factors manage to reconvert them.

Table 3. Coverage for use of soil and vegetation of the productive potential of vanilla in the Totonacapan region.

*= includes mangrove, aquatic vegetation and halophyte. MA= very high; A= high; M= medium; B= low; MB= very low.

The results obtained by MaxEnt and MCE indicate that the areas with the greatest productive potential are located within the grassland areas that are destined for livestock, this information was validated through on-site tours with directors of the Association of Vanilla Producers of the State of Puebla AC and served as the basis so that between May and June 2017, vanilla is planted 100 ha to the northwest of the Totonacapan region, considering that in these areas with greater productive potential (high to very high) the probability of adaptation is increased, development and growth of vanilla.

Both methodologies agree that the areas with the greatest productive potential are located on agricultural, pasture and forest areas. It is proposed the cultivation of vanilla in those areas whose current use is destined to livestock (pasture) and to the agriculture of corn, beans, squash and pepper mainly. Barrera-Rodríguez et al. (2011); Jaramillo et al. (2012 and 2013) establish that to increase the areas planted with vanilla in the Totonacapan region, the profitability of the crop must be improved in the face of other activities, giving an added value to the product that is marketed in green through the beneficiation process.

Conclusions

The production of vanilla in the Totonacapan region is a secondary activity, among producers of advanced age (>55 years), low level of studies (<8 years) and in a smallholding (<0.5 ha). In addition to these factors, phytosanitary problems, fertilization, access to credit and market, have also contributed to a decline in the productivity of vanilla in the region, marked by the reduction in the area sown.

The highest productive potential for the cultivation of vanilla is found in the northwest of the Totonacapan region, in areas that are mostly destined for cattle raising and agriculture with more profitable cyclical crops, between the limits of Puebla and Veracruz. The productive potential is conditioned by the increase in rainfall in the driest month (April: <60 mm year-1) and temperatures no higher than 30 °C in the warmest quarter (June-July-August), so that the irrigation can substantially help improve yields in the field and consequently productivity in the region.

The productive potential for the sowing of new areas with vanilla in the Totonacapan region, obtained by means of geographic techniques, helped to make decisions in relation to the expansion of this crop, towards those surfaces with greater probability of adaptation for the development and optimal growth of the vanilla, with what is expected to increase the productivity of the national vanilla producer sector and therefore its international competitiveness. Both methods should not be considered different but complementary; however, the limitation of these techniques to socioeconomic, market and crop management factors should be considered, in virtue of which it is necessary to complement them with fieldwork.

Acknowledgments

To C. Miguel Ángel García Cabrera, technical secretary of the Association of Vanilla Producers of the State of Puebla, for the facilities provided for the collection of information in the field. To the anonymous reviewers of the paper, for their comments, which helped to enrich the investigation.

Cited literature

Aguilar, R. N.; Galindo, M. G.; Fortanelli, M. J. y Contreras, S. C. 2010. Evaluación multicriterio y aptitud agroclimática del cultivo de caña de azúcar en la región de Huasteca. Ciencia y Tecnología Agropecuaria. 11(2):144-154.

Alcón, F.; De Miguel, M. D. y Burton, M. 2008. Adopción de tecnología de distribución y control del agua en las comunidades de regantes de la región de Murcia. Economía Agraria y Recursos Naturales. 8(1):83-102.

Álvarez, L. C. J. y Crecente, M. R. 2000. Una revisión de la concentración parcelaria en Europa. Estudios Agrosociales y Pesqueros. 187:221-274.

ArcMap. 2010. Environmental Scientific Research Institute. Arc Map 10.0. Redlands, California, USA. https://www.esri.com/es-es/home.

Barrera, R. A.; Santoyo-Cortés, V.; Baca, D. e Altamirano, C. J. 2014. Perspectives d’avenir et compétitivité des organisations de producteurs de vanille du Mexique. Cahiers Agricultures. 23(6):374-381.

Barrera-Rodríguez, A. I.; Herrera-Cabrera, B. E.; Jaramillo-Villanueva, J. L.; Escobedo-Garrido, J. S. y Bustamante-González, Á. 2009. Caracterización de los sistemas de producción de vainilla (Vanilla planifolia A.) bajo naranjo y en malla sombra en el Totonacapan. Trop. Sub. Agroecosys. 10(2):199-212.

Barrera-Rodríguez, A. I.; Jaramillo-Villanueva, J. L.; Escobedo-Garrido, J. S. y Herrera-Cabrera, B. E. 2011. Rentabilidad y competitividad de los sistemas de producción de vainilla (Vanilla planifolia J.) en la región del Totonacapan, México. Agrociencia. 45(5):625-638.

Ceballos-Silva, A. P. y López-Blanco, J. 2010. Delimitación de áreas adecuadas para cultivos de alternativa: una evaluación multicriterio-sig. Terra Latinoam. 28(2):109-118.

CIMMYT. 1993. Centro Internacional de Mejoramiento de Maíz y Trigo. La adopción de las tecnologías agrícolas: Guía para el diseño de encuestas. http://libcatalog.cimmyt.org/ download/cim/42408.pdf.

DivaGis. 2016. DIVA-GIS. Version 7.5. A geographic information system for the analysis of species distribution data. http://www.diva-gis.org/download.

Elith, J.; Phillips, S.; Hastie, T.; Dudík, M.; Chee, Y. y Yates, C. 2011. A statistical explanation of MaxEnt for ecologists. Div. Distrib. 17(1):43-57.

FAO. 1999. Organización de las Naciones Unidas para la Agricultura y la Alimentación. Centro Internacional de Referencia e Información en Suelos y Sociedad Internacional de las Ciencias del Suelo (FAO-ISRIC y SICS). Base Referencial Mundial del Recurso Suelo. Reporte No. 84. Roma, Italia. 93 p.

FAOSTAT. 2018. The Food and Agriculture Organization Corporate Statistical. Production. http://faostat3.fao.org/download/Q/QC/S.

Flores, J. A.; Reyes, L. D.; Jiménez, G. D.; Romero, A. O.; Rivera, T. J. A.; Huerta, L. M. y Pérez, S. A. 2017. Diversidad de Vanilla spp (Orchidaceae) y sus perfiles bioclimáticos en México. Rev. Biol. Trop. 65(3):975-988.

Gómez, M. y Barredo, J. I. 2005. Sistemas de información geográfica y evaluación multicriterio en la ordenación del territorio. Ra-Ma. Madrid. 279 p.

Hernández, P. A.; Graham, H. C.; Master, L. L. and Albert, L. D. 2006. The effect of sample size and species characteristics on performance of different species distribution modeling methods. Ecography. 29(5):773-785.

Hernández-Ruíz, J.; Herrera-Cabrera, B. E.; Delgado-Alvarado, A.; Salazar-Rojas, V. M.; Bustamante-González, Á.; Campos-Contreras, J. E. y Ramírez-Juárez, J. 2016. Distribución potencial y características geográficas de poblaciones silvestres de Vanilla planifolia (Orchidaceae) en Oaxaca, México. Rev. Biol. Trop. 64(1):251-262.

Herrera-Cabrera, B. E.; Salazar-Rojas, V. M.; Delgado-Alvarado, A.; Contreras, J.; Contreras, C. y Cervantes-Vargas, J. 2012. Use and conservation of Vanilla planifolia J. in the Totonacapan Region, México. Eur. J. Environ. Sci. 2(1):43-50.

Hijmans, R. J.; Cameron, E. S.; Parra, L. J.; Jones, G. P. and Jarvis, A. 2005. Very high-resolution interpolated climate surfaces for global land areas. Inter. J. Climatol. 25(15):1965-1978.

Hijmans, R. J.; Garrett, K. A.; Huaman, Z.; Zhang, D. P.; Schreuder, M. and Bonierbale, M. 2000. Assessing the geographic representativeness of genebank collections: the case of Bolivian wild potatoes. Conservation Biol. 14(6):1755-1765.

INEGI. 2016. Instituto Nacional de Estadística, Geografía e Informática. Conjunto de datos vectoriales de uso de suelo y vegetación. Escala 1:250 000. Serie VI. http://www.conabio.gob.mx/informacion/metadata/gis/usv250s6gw.xml?-httpcache=yes& -xsl=/db/metadata/xsl/fgdc-html.xsl&-indent=no.

INIFAP. 2011. Instituto Nacional de Investigación Forestal, Agrícola y Pesquera. Paquete tecnológico vainilla (Vainilla planifolia Jackson) establecimiento y mantenimiento. www.inifap.gob.mx/Documents/inicio/paquetes/vainilla-establecimiento.pdf.

Jaramillo, V. J. L.; Escobedo, G. J. D. y Barrera, R. A. 2012. Competitividad de sistemas de beneficiado de vainilla (Vanilla Planifolia J.) en la región del Totonacapan, México. Panorama Socioeconómico. 30(45):80-93.

Jaramillo, V. J. L.; Garrido, J. S. E.; Rodríguez, A. B. y Cabrera, B. E. H. 2013. Eficiencia económica en el beneficiado de vainilla (Vanilla planifolia J.) en la región del Totonacapán, México. Rev. Mex. Cienc. Agríc. 4(3):477-483.

Jiménez, C. A.; Vargas, T. V.; Salinas, C. W. E.; Aguirre, B. M. de J. y Rodríguez, C. D. 2004. Aptitud agroecológica para el cultivo de caña de azúcar en el sur de Tamaulipas, México. Investigaciones Geográficas. 53:58-74.

Kazuya, N.; Gómez, M.; López, R.; Meneses, R. y Vargas, J. 2006. Comparación de modelos de distribución de especies para predecir la distribución potencial de vida silvestre en Bolivia. Ecol. Bolivia. 41(1):65-78.

Lubinsky, P.; Bory, J. S.; Hernández, H. S.; Kim, C. and Gómez, P. G. 2008. Origins and dispersal of cultivated Vanilla (Vanilla planifolia Jacks. [Orchidaceae]). Econ. Bot. 62(2):127-138.

MaxEnt. 2016. Maxent software for modeling species niches and distributions (Version 3.4.1). http://biodiversityinformatics.amnh.org/open-source/maxent/.

Navarro-Cerrillo, R. M.; Hernández-Bermejo, J. E. and Hernández-Clemente, R. 2011. Evaluating models to assess the distribution of Buxus balearica in southern Spain. Appl. Vegetation Sci. 14(2):256-267.

Olivas, G. U.; Valdez, L. J.; Aldrete, A.; González, G. M. y Vera, C. G. 2007. Áreas con aptitud para establecer plantaciones de maguey cenizo: definición mediante análisis multicriterio y SIG. Rev. Fitotec. Mex. 30(4):411-419.

Phillips, S.; Anderson, R. and Schapire, R. 2006. Maximum entropy modeling of species geographic distributions. Ecol. Modelling. 190(3-4):231-259.

Rocha, F. R. G.; Herrera, C. B. E.; Velasco, V. J.; Salazar, R. V. M.; Delgado, A. A. y Mendoza, C. M. C. 2018. Determinación preliminar de componentes de rendimiento para el cultivo de vainilla (Vanilla planifolia Jacks. ex Andrews) en la región Totonacapan, México. Agroproductividad. 11(3):9-14.

SAGARPA. 2012. Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación. Plan rector sistema producto vainilla. http://dev.pue.itesm.mx/sagarpa/estatales/ept% 20comite%20sistema%20producto%20vainilla%20puebla/plan%20rector%20que%20contiene%20programa%20de%20trabajo%202012/pr-vainilla-puebla-%202012.pdf.

Salazar-Rojas, V. M.; Herrera-Cabrera, B. E.; Delgado-Alvarado, A.; Soto-Hernández, M.; Castillo-González, F. and Cobos-Peralta, M. 2012. Chemotypical variation in Vanilla planifolia Jack (Orchidaceae) from the Puebla-Veracruz Totonacapan region. Gen. Res. Crop Evol. 59(5):875-887.

Santillán-Fernández, A.; Salas-Zúñiga, A. y Vásquez-Bautista, N. 2018. La producción de la vainilla (Vainilla planifolia Jacks. Ex Andrews) en México de 2003 a 2014. Rev. Mex. Cienc. Fores. 9(47):050-060.

Scheldeman, X.; Willemen, L.; d’Eeckenbrugge, G. C.; Romeijn, P. E.; Restrepo, M. T.; Motoche, J. R.; Jiménez, D.; Lobo, M.; Medina, C. I.; Reyes, C.; Rodríguez, D.; Ocampo, J. A.; Damme, V. P. and Goetgebeur, P. 2007. Distribution, diversity and environmental adaptation of highland papayas (Vasconcellea spp.) in tropical and subtropical America. Bio. Conserv. 16(6):1867-1884.

SIAP. 2018. Sistema de Información Agrícola y Pesquera. Producción Agrícola. http://www.siap.gob.mx/.

Sinha, A. K.; Sharma, U. K. and Sharma, N. 2008. A comprehensive review on Vanilla flavor: extraction, isolation and quantification of vanillin and other constituents. Inter. J. Food Sci. Nutr. 59(4):299-326.

Soto, J. L.; Hartwich, F.; Arispe, T.; Monge, M. y Ampuero, L. 2006. Innovación en el cultivo de Quinua en Bolivia: efectos de la interacción social y de las capacidades de absorción de los pequeños productores. Rev. Análisis. 1(3):17-20.

Xochipa-Morante, R. C.; Delgado-Alvarado, A.; Herrera-Cabrera, B. E.; Escobedo-Garrido, J. S. y Arévalo-Galarza, L. 2016. Influencia del proceso de beneficiado tradicional mexicano en los compuestos del aroma de Vanilla planifolia Jacks. ex Andrews. Agroproductividad. 9(1)55-62.