Revista Mexicana de Ciencias Agrícolas publicación especial número 23 28 de septiembre - 11 de noviembre, 2019
DOI: https://doi.org/10.29312/remexca.v0i23.2019
Artículo
Acumulación de horas frío para la producción de arándano en Nayarit, México
Arturo Álvarez-Bravo1
Rubén Bugarín-Montoya2§
Mairim Elizabeth Arellano-Figueroa3
1Campo Experimental Santiago Ixcuintla-INIFAP. Nayarit, México. CP. 63300. (alvarez.arturo@inifap.gob.mx). 2Universidad Autónoma de Nayarit-Unidad Académica de Agricultura. Amado Nervo s/n, Col. Los Fresnos, Ciudad de la Cultura Tepic, Nayarit, México. CP. 63155. 3Investigador independiente. (maieli51@gmail.com).
§Autor para correspondencia: drbugarin@hotmail.com.
Resumen
El clima es un factor determinante del rendimiento y calidad en sistemas de producción agrícola. El objetivo de este trabajo fue cuantificar las horas frío-acumulables en Nayarit y utilizar dicha información como indicador agroclimático en la regionalización para la producción de arándano (Vaccinium corymbosum L.). En la investigación se emplearon datos de la red de estaciones agrometeorológicas de Nayarit, la cual se compone de 38 estaciones automatizadas que colectaron datos quinceminutales de 11 variables. Con registros de temperatura mínima se calcularon las horas frío, las cuales se organizaron por mes y acumulación por año. Mediante un sistema de información geográfica se representó cartográficamente la distribución de las horas frío anual (HFA), con interpolación y organización de los resultados en cuatro fases (<299, 300-399, 400-499 y 500-699). Para cada clase de HFA, se cuantificó la superficie en cada uno de los 20 municipios. Finalmente se evaluó la precisión del modelo de interpolación contra una base de datos ajena al conjunto de datos principal. Los resultados muestran cómo solamente en el 27.4% de la superficie estatal se acumulan más de 300 HFA que es el requerimiento mínimo para arándano. Las condiciones térmicas para arándano se ubican principalmente en el centro y sur del estado, resaltando los municipios de Compostela, Xalisco, Tepic, Santa María del Oro, Jala, La Yesca e Ixtlán del Río que representan el 81.6% de la superficie apta. Los resultados son un precedente sobre la acumulación de horas de frío en Nayarit y funge como instrumento de planeación para tomadores de decisión, técnicos y productores interesados en el cultivo de arándano.
Palabras clave: agroclimatología, horas frío, temperatura.
Recibido: enero de 2019
Aceptado: abril de 2019
Introducción
En las últimas décadas, el estado de Nayarit ha experimentado un proceso de reconversión productiva, desde cultivos básicos hacia frutales y hortalizas en invernadero (SIAP, 2017). En este proceso de transición se ha incursionado en la adopción de nuevos cultivos y tecnologías de producción, particularmente el de frutillas bajo condiciones de agricultura protegida. Esta última ofrece nueva alternativa de reconversión con alta rentabilidad en poca superficie (SIAP, 2017). A nivel mundial, la producción de frutillas o ‘Berries’ en inglés, como la fresa, arándano, zarzamora o frambuesa, ha cobrado relevancia económica, política, social y se encuentra en pleno crecimiento y desarrollo (Cruz, 2018).
En México, como respuesta a la demanda creciente de frutillas en el mundo, se han desarrollado proyectos de diversificación de cultivos en Michoacán, Baja California, Nayarit y Puebla, entre otros, las cuales cuentan con características climáticas óptimas en diferentes periodos estacionales de producción (FIRA, 2016). En Sinaloa o Sonora, se están buscando variedades que se adapten a condiciones cálidas con el objeto de mejorar la rentabilidad (González, 2013).
El arándano es una de las especies de reciente introducción en la cadena agroalimentaria en México, su producción se remonta a 1996 y en la última década ha tenido un crecimiento mayor a 800%, debido a la demanda del producto en Europa, Asia y Norte América (Pérez, 2018). En el país hay establecidas 2 625 ha de arándano (SIAP, 2017) que generan entre 100 y 110 mil empleos directos e indirectos (FAOSTAT, 2017). El estado de Jalisco se coloca en primer lugar en producción de arándano, con 14 563 t en 1 576 ha y un valor de producción de 524 millones de pesos (SIAP, 2017).
En general, el precio al productor por tonelada de arándano ascendió a 51 966 pesos en ese año y fue más alto comparado con el resto de las frutillas. Este tipo de cultivos, por su alto valor, potencia la rentabilidad de pequeñas superficies y es viable integración de pequeños productores (FIRA, 2016). México produce 8.78 t ha-1, posicionándose con mejor rendimiento, seguido por Italia con 7.5 t ha-1, Rumania 6.66 t ha-1 y Estados Unidos con 6.46 t ha-1 (Hernández y Gutiérrez, 2013).
El arándano se cultiva desde los 600 hasta 2 500 msnm en diversas partes del mundo (Paredes, 2010). En general, la planta de arándano durante el otoño en climas fríos presenta un periodo de dormancia debido a la presencia de bajas temperaturas y fotoperiodo corto, por lo que requiere cierto número de horas frío (HF) para la iniciación floral y crecimiento de hojas en la primavera (Retamales y Hancock, 2012). Las variedades actuales disponibles tienen requerimientos de frío desde 150 a 800 HF en arándano tipo arbusto alto del sur (southern highbush), hasta 800 a 1200 HF para el tipo arbusto alto del norte (northern highbush), y de 300 a 600 HF en variedades tipo ‘ojo de conejo’ (rabbiteye) (Retamales y Hancock, 2012).
Asimismo, las variedades de arándano se clasifican en tres grupos: requerimiento ‘alto’ mayor a 800 horas frío, ‘medio’ de 400 a 600 HF y ‘bajo’ menor a 400 HF (García, 2011). Se ha propuesto considerar que los requerimientos de HF para arándano se satisfacen en el intervalo de temperaturas arriba de 1.4 y menores de 12.4 °C (Retamales y Hancock, 2012). Cesaraccio et al. (2004 y 2006) describieron la importancia de acumulación de frío seguido de un periodo cálido el cual permite la emergencia de los brotes reproductivos en arándano. Por lo tanto, el cálculo de HF es importante, ya que permite definir los sitios potenciales de cultivo (Mainland, 1985).
El arándano a pesar de ser una especie adaptada de climas templados y fríos, la producción en México es posible porque algunas variedades como ‘Biloxy’, ‘Victoria’, ‘Kester’, ‘Rocío’ y ‘Corona’ entre otras, tienen bajo requerimiento de horas frío y se adaptan a la mayoría de los microclimas tropicales y subtropicales presentes (Salgado et al., 2018). En la actualidad, se han introducido a nuestro país, nuevas variedades de arándano que no requieren horas-frío, y que han sido seleccionadas para los sistemas de producción siempre verde, como es el caso de BiancaBlue™ ‘FCM12-087’, AtlasBlue™ ‘FCM12-045’ y Jupiter Blue ‘FCM12-131’ por la empresa Fall Creek Farm & Nursery Inc., lo cual, sin duda, favorecerá una mayor expansión del cultivo de arándano en zonas subtropicales y tropicales.
La precocidad productiva es la principal ventaja comercial de las variedades de arándanos de baja exigencia de horas frío, cuyos frutos alcanzan los precios más altos durante el inicio de la cosecha del hemisferio norte en marzo-abril (Cantuarias-Avilés et al., 2014). Dada la expansión actual y futura que se prevé del arándano en el estado de Nayarit, es importante contar con un instrumento técnico, validado científicamente que permita identificar las zonas con potencial climático que requiere el cultivo. El objetivo del presente estudio fue identificar y cuantificar la superficie con condiciones agroclimáticas óptimas para la producción de arándano en el estado de Nayarit.
Los resultados abonan a un mejor entendimiento de las condiciones ambientales que prevalecen en las zonas agroecológicas y es una herramienta útil para los tomadores de decisión, técnicos y productores.
Materiales y métodos
Área de estudio
El trabajo consideró el territorio del estado de Nayarit, entidad ubicada en el occidente de México. Localizada en, entre los paralelos 23° 05’ 04’’ y 20° 36’ 12’’ latitud norte y entre los meridianos de 103° 43’ 15’’ y -105° 45’ 37’’ longitud oeste (INEGI, 2015).
Datos meteorológicos
Se utilizaron datos de la red de estaciones agrometeorológicas de Nayarit, la cual se compone de equipos automatizados que colectan datos quinceminutales de temperatura, precipitación, radiación solar y viento (Adcon Telemetry, modelo A753, Klosterneuburg, Austria). La red operó durante ocho años desde 2007 con 38 estaciones (Cuadro 1).
Cuadro 1. Descripción de las estaciones meteorológicas.
Nombre | Latitud | Longitud | Altitud (m) | Nombre | Latitud | Longitud | Altitud (m) |
Acaponeta | 22.48 | 105.403 | 15 | Santa María del Oro | 21.341 | 104.635 | 1084 |
Estancia de los López | 20.852 | 104.434 | 892 | Colonia Moderna | 21.467 | 104.66 | 858 |
Valle de Banderas | 20.784 | 105.242 | 62 | Villa Juárez | 21.695 | 105.392 | 57 |
El Capomo | 21.116 | 105.156 | 40 | Santiago Ixcuintla | 21.824 | 105.184 | 10 |
Ixtapa de la Concepción | 21.301 | 105.192 | 20 | El Verdineño | 21.702 | 105.132 | 43 |
Monteón | 20.974 | 105.306 | 21 | Santa Cruz | 21.979 | 105.579 | 1 |
Compostela | 21.231 | 104.884 | 861 | Pozo de Ibarra | 21.872 | 105.276 | 32 |
Mesa del Nayar | 22.214 | 104.647 | 1403 | Quimichis | 22.368 | 105.537 | 9 |
Huajicori | 22.636 | 105.331 | 75 | El Limón | 22.301 | 105.467 | 3 |
Ixtlán del Río | 21.024 | 104.362 | 1131 | San Felipe Aztatán | 22.399 | 105.395 | 22 |
Rosa Blanca | 21.118 | 104.358 | 1936 | Atonalisco | 21.653 | 104.827 | 415 |
Puente de Camotlán | 21.7 | 104.089 | 1113 | V. Carranza | 21.525 | 104.974 | 1063 |
Rosamorada | 22.095 | 105.217 | 25 | Xalisco | 21.425 | 104.892 | 974 |
Guadalupe Victoria | 21.667 | 105.326 | 3 | Malinal | 21.367 | 105.018 | 864 |
Las Palmas | 21.605 | 105.143 | 186 | San Pedro L. | 21.201 | 104.757 | 1271 |
Tequilita | 21.104 | 104.807 | 979 |
Gestión de los datos
Mediante el motor de base de datos Access 2013 (Microsoft Office 2010, Redmond, Washington, Estados Unidos) se integraron los datos de temperatura organizados por estación meteorológica, fecha y hora.
Cálculo de horas frío acumulado
Por cada estación meteorológica, mes y año se cuantificó las horas frío (HF). Las HF son los registros de temperatura menor o igual a 12 °C. Este umbral de temperatura fue sugerido por Norvell & Moore, 1982. El cálculo se realizó mediante la siguiente ecuación:
HFAm = (ΣRT ≤ 12 °C )/ 4
Donde: HFAm son las horas frío acumuladas por mes y RT son los registros de temperatura quinceminutales menor o igual a 12 ºC.
Conjunto de datos de horas frío anual
Para el análisis de los datos anuales (HFA) se promediaron los valores mensuales por la serie de tiempo (8 años) para cada estación. Resultaron dos bases de datos, una para la interpolación constituida por los registros de 31 estaciones y otra, con registros de siete estaciones, que representan una muestra de 20% para la validación de la interpolación (Cuadro 1).
Representación cartográfica
Él análisis espacial, se realizó en el sistema de información geográfica Arcmap Versión 10.1 (ArcGIS ESRI, Redlands, California, Estados Unidos). El método de interpolación utilizado para la variable HFA fue IDW (Distancia Inversa Ponderada) y el resultado de este procedimiento fue clasificado de acuerdo al Cuadro 2.
Cuadro 2. Clasificación de la acumulación de horas frío anual y las variedades de arándano aptas para esa condición.
Umbral | Descripción y condiciones de horas-frío por variedad1 | |
Inferior | Superior | |
0 | 299 | Acumulación de horas frío insuficiente para muchas variedades, excepto variedades para cultivo siempre verde, con bajo o nulo requerimiento de frío. Por ejemplo, Biloxi (150 HF), Victoria (sin requerimiento de frío) |
300 | 399 | Condiciones para variedades Jewel, Rabbiteye y Misty |
400 | 499 | Condiciones para variedades Sharblue y O’Neal |
500 | 699 | Condiciones para variedades Jubilee y Ozarkblue |
1= (Norvell y Moore, 1982; Darnell y Davies, 1990; Darnell, 1992; García, 2011; Salgado et al., 2018).
El mapa de acumulación anual de horas frío en Nayarit se hizo en Arcmap Ver. 10.1, con el método de estadística espacial ArcToolBox se calculó la superficie estatal y municipal para cada clase.
Validación de la interpolación
Se compararon los datos de horas frío acumulado resultado de la interpolación o datos predichos contra los datos observados de una muestra aleatoria de siete estaciones (Cuadro 3). Con el programa estadístico Minitab versión 17 (Universidad Estatal de Pensilvania, Estados Unidos), el análisis se elaboró con el modelo de regresión lineal, calculando el coeficiente de correlación.
Cuadro 3. Descripción de estaciones meteorológicas empleadas en la validación.
Nombre | Latitud | Longitud | Altitud (m) |
San Juan de Abajo | 20.837 | 105.211 | 70 |
Las Varas | 21.192 | 105.147 | 14 |
Jala | 21.081 | 104.438 | 1 045 |
Huajimic | 21.671 | 104.342 | 1 205 |
Amapa | 21.812 | 105.232 | 31 |
Tecuala | 22.4 | 105.475 | 11 |
Tepic | 21.488 | 104.89 | 946 |
Análisis del conjunto de datos mensuales
Se realizó un análisis de estadística descriptiva de los datos diarios organizados por mes, con el apoyo de una gráfica de caja elaborada con el programa estadístico Minitab versión 17.
Resultados y discusión
Horas frío acumulado anual de la zona de estudio
En Nayarit 72.6% (2 038 724 has) del territorio, acumula anualmente menos de 300 horas frío. Las restantes 768 549 ha acumulan desde 300 hasta 700 horas frío anualmente, las cuales según (Norvell y Moore, 1982; INTAGRI, 2017; Salgado et al., 2018) son suficientes para la producción de arándano de variedades tipo ‘Ojo de conejo’, así como tipo arbusto alto del sur Jewel, Misty, Sharpblue, O’Neal, Jubilee o Ozarkblue (Figura 1).
Figura 1. Clasificación de las horas frío acumulado anual en Nayarit.
Las condiciones identificadas entre 300-500 HF anuales, que representan 18.3% de la superficie estatal, pueden ser aptas para la variedad Biloxi lo cual coincide con Salgado et al. (2018). Las HF demandadas por algunas variedades de arándano de bajo y medio requerimiento de frío según García (2011) fueron encontradas en el presente trabajo por lo que en Nayarit es altamente factible su cultivo.
Acumulación anual de horas frío por municipio
En 19 municipios fue posible cuantificar superficie donde se acumula menos de 300 horas frío anualmente. Las 768 549 ha que corresponden a la acumulación de frío mayor a 300 horas anuales fueron ubicadas en los municipios de Ahuacatlán, Amatlán de Cañas, Compostela, Del Nayar, Ixtlán del Río, Jala, La Yesca, San Pedro Lagunillas, Santa María del Oro, Tepic y Xalisco. La superficie en la que se acumulan más de 500 horas frío anual corresponde a los municipios de Ahuacatlán, Del Nayar, Ixtlán del Río, Jala, La Yesca, San Pedro Lagunillas, Santa María del Oro y Tepic, la cual representa 256 125 ha o 9.1% del territorio estatal. De estos, sobresalen los municipios de La Yesca y Santa María del Oro con cerca de 180 000 ha (Cuadro 4).
Cuadro 4. Superficie en hectáreas con horas frío-acumuladas durante el año.
Municipio | Acumulación anual de horas frío | |||
0-299 | 300-399 | 400-499 | 500-699 | |
Acaponeta | 142 610.4 | |||
Ahuacatlán | 25 469 | 13 323.9 | 8 024.3 | 3 623.2 |
Amatlán de Cañas | 50 266 | 1 529.9 | ||
Bahía de Banderas | 77 073 | |||
Compostela | 152 478.1 | 7 785.8 | 27 642.2 | |
Del Nayar | 489 687.3 | 10 008.8 | 6 660.3 | 7 482.5 |
Huajicori | 223 538.3 | |||
Ixtlán del Río | 4 010.6 | 26 943.6 | 3530.5 | 14 781.9 |
Jala | 38.6 | 941.7 | 49 344.4 | |
La Yesca | 6 824.3 | 147 103.3 | 185 904.5 | 91 569.4 |
Rosamorada | 183 914.7 | |||
Ruíz | 52 016.5 | |||
San Blas | 110 328.9 | |||
San Pedro Lagunillas | 44 050.7 | 5 308.4 | 1 928.8 | 237.9 |
Santa María del Oro | 2 275.8 | 8 484.5 | 11 295.3 | 86 998.3 |
Santiago Ixcuintla | 172 747.4 | |||
Tecuala | 104 399.2 | |||
Tepic | 137 124.8 | 8 884 | 15 294.6 | 2 086.7 |
Tuxpan | 31 367.6 | |||
Xalisco | 28 541.4 | 3 088.4 | 18 703.5 |
García (2011), clasificó los requerimientos de HF para arándano en: bajo (<400), medio (400-600) y alto (>800), lo anterior corresponde a nuestro estudio al 8.3% y 19.1% de la superficie estatal del nivel bajo y medio respectivamente.
Distribución espacial de la acumulación anual de horas frío
En la Figura 2 se observó cómo la región donde se acumulan menos de 300 HFA se distribuyen por toda la costa, los municipios desde Tepic hasta Huajicori al norte y en el sur, la región de San Pedro Lagunillas el sur de Ahuacatlán y Amatlán de Cañas. La región con condiciones óptimas se concentra en el centro y sur del estado, principalmente en los municipios de Santa María del Oro, Jala, norte de Ixtlán del Río y sur de La Yesca.
El sistema de información geográfica permitió la representación cartográfica de datos climáticos, así como cuantificar la superficie de las cuatro clases de acumulación anual de HF en Nayarit. Los resultados permitirán localizar puntualmente las zonas con mejores condiciones para el establecimiento de cultivos con requerimientos térmicos específicos, como las horas frío en arándano, similares investigaciones fueron reportadas por Bhatt et al. (2018); Gentilucci et al. (2019) para kiwi (Actinidia deliciosa) en India y uva para vino (Vitis vinífera) en Italia respectivamente. Por su parte, Paredes (2010) señala que el arándano se cultiva con éxito desde los 600 hasta 2 500 msnm lo cual concuerda con nuestros resultados en particular con las zonas con acumulación anual de horas frío >400 que se ubican en altitudes de 900 a 2 000 m.
Figura 2. Distribución espacial de la acumulación anual de horas frío en Nayarit.
Validación del modelo de interpolación
En la Figura 3, se muestra el resultado de relacionar los datos generados por la interpolación del modelo IDW (predichos) con los observados. En 71.4% de los puntos de validación, la interpolación subestima (residuos menores a 25 horas frío), aun así, el ajuste del modelo con un coeficiente de correlación de 0.97 es muy bueno, solamente en dos puntos de validación la interpolación sobrestima con residuos desde 9.9 hasta 16.9%.
Figura 3. Correlación entre los valores observados y predichos por la interpolación de las horas frío acumulados anuales.
Los resultados contrastan con lo reportado por Álzate-Velásquez et al. (2018), donde evaluaron el modelo de interpolación IDW con una R2 de 0.76, en este estudio el modelo de regresión reportó mejores resultados con 0.97 en R2. Yang y colaboradores (2015) reportaron resultados similares a los del presente trabajo, al evaluar la precisión del mismo modelo interpolación (R2 de 0.94).
Distribución mensual de las horas frío
En la región donde se acumula más de 300 HF por año, la mayor variación mensual se presenta en noviembre y diciembre (con valores medios >60 HF), siendo los meses cuando se inicia el periodo de mayor acumulación de HF de noviembre a marzo. En este mismo periodo, 50% de los datos fluctúan entre 40 y 115 HF, lo que indica una alta variabilidad interanual en la región. La menor acumulación se presentó entre junio a septiembre (promedio >15 HF), iniciando un descenso de HF en abril (50 HF promedio) hasta llegar a la menor acumulación en junio con 8.6 HF (Figura 4).
Figura 4. Distribución mensual de horas frío para la región con más de 300 HFA.
Conclusiones
Por primera vez en Nayarit fue posible identificar y cuantificar la superficie con condiciones agroclimáticas óptimas para la producción de arándano. La cantidad de horas frío acumuladas anualmente permiten la producción de variedades comerciales de bajo requerimiento. Considerando la buena precisión del modelo de interpolación, las condiciones óptimas se encuentran en 27% de la superficie estatal (700 mil hectáreas) distribuidas en 11 municipios.
El periodo que acumula la mayor cantidad de horas frío en Nayarit es de noviembre a marzo. Estos resultados servirán para el diseño de una estrategia de reconversión productiva de suelos con baja productividad, nula aptitud agrícola o condición ociosa; para transitar hacia un cultivo de alto valor comercial que demanda mano de obra y de servicios asociados a la cadena de valor.
Literatura citada
Álzate, V. D. F.; Araujo, C. G. A.; Rojas, B. E. O.; Gómez, L. D. A. y Martínez, M. F. E. 2018. Interpolación Regnie para lluvia y temperatura en las regiones andina, caribe y pacífica de Colombia. Colombia Forestal. 21(1):102-118.
Bhatt, C. K.; Nain, A. S.; Bhatt, M. K. and Paliwal, A. 2018. Site suitability analysis for kiwi fruit plantation in Uttarakhand using GIS. IJCS. 6(3):101-106.
Cantuarias-Avilés, T.; Rodríguez-Da Silva, S; Bordignon-Medina, R.; García-Moraes A. F. e Freire-Alberti, M. 2014. Cultivo do mirtilo: atualizações e desempenho inicial de variedades de baixa exigência em frio no estado de São Paulo. Rev. Brasileira de Fruticultura. 36 (1): 139-144.
Cesaraccio, C.; Spano, D.; Snyder, R. L. and Duce, P. 2004. Chilling and forcing model to predict bud-burst of crop and forest species. Agric. For. Meteorol. 126(1-2):1-13.
Cesaraccio, C.; Spano, D.; Snyder, R. L.; Duce, P. and Jones, H. G. 2006. Improvement of chilling and forcing model to predict bud-burst. In: Proceedings of the 17th Conference on Biometeorology and Aerobiology. San Diego, CA. American Meteorological Society. 1-4 pp.
Cruz, O. A. P. 2018. Análisis de la cadena productiva del arándano en México y Chile. PORTES, Revista Mexicana de Estudios sobre la Cuenca del Pacífico. 12(23):31-62.
Darnell, R. L. and Davies, F. S. 1990. Chilling accumulation, budbreak, and fruit set of young rabbiteye blueberry plants. HortSci. 25(6):635-638.
FAOSTAT. 2017. Base de datos estadísticos corporativos de la Organización de Alimentos y Agricultura. http://www.fao.org/faostat/es/#country/138.
FIRA. 2016. Fideicomisos Instituidos en Relación con la Agricultura. Panorama Agroalimentario, Berries 2016. https://www.gob.mx/cms/uploads/attachment/file/200633/Panorama-Agroalimentario-Berries-2016.pdf
García, R. J. C. 2011. El cultivo del arándano en Asturias. Tecnología agroalimentaria. Boletín informativo del SERIDA no. 9:13-20. https://ria.asturias.es/RIA/handle/123456789/1482.
Gentilucci, M.; Barbieri, M. and Burt, P. 2019. Climate and Territorial Suitability for the Vineyards Developed Using GIS Techniques. In: exploring the nexus of geoecology, geography, geoarcheology and geotourism: advances and applications for sustainable development in environmental sciences and agroforestry research. Springer, Cham. 11-13 pp.
González, C. 2013. Alternativas para el cultivo de arándanos. Oficina de estudios y política agraria (ODEPA). Ministerio de Agricultura, Santiago de Chile. 6 p.
Hernández, G. J. A. y Gutiérrez, P. H. 2013. Agenda de innovación tecnológica, Jalisco. 52 p.
INEGI. 2015. Instituto Nacional de Estadística, Geografía e Informática. Anuario Estadístico y Geográfico de Nayarit 2017. Aguascalientes, Ags. 469 p.
INTAGRI. 2017. Variedades comerciales de arándanos en el mundo. Serie frutillas núm. 15. Artículos técnicos de INTAGRI. México. 4 p.
Mainland, C. 1985. Some problems with blueberry leafing, flowering and frutting in a warm climate. III International Symposium on Vaccinium Culture. ISHS Acta Horticulturae. 165:35-46.
Norvell, D. J. and Moore, J. N. 1982. An evaluation of chilling models for estimating rest requirements of highbush blueberries (Vaccinium corymbosum L.). J. Am. Soc. Hortic. Sci. 107:54-56
Paredes, J. I. de S. 2010. Los frutos del bosque o pequeños frutos en la cornisa cantábrica: el arándano. Gobierno de Cantabria. Consejería de Desarrollo Rural, Ganadería, Pesca y Biodiversidad. Cantabria, España. 151 p.
Pérez, O. 2018. Análisis de la cadena productiva del arándano en México y chile. Rev. Mexicana de Estudios sobre la Cuenca del Pacífico. 12(23):31-62.
Retamales, J. and Hancock, J. F. 2012. Blueberries. Crop Production Science in Horticulture Series UK. Núm. 21. 323 p.
Salgado, C.; Sánchez, P.; Volke, V. and Colinas, M. 2018. Respuesta agronómica de arándano (Vaccinium corymbosum L.) al estrés osmótico. Agrociencia. 52(2):231-239.
SIAP. 2017. Atlas agroalimentario 2017. Servicio de Información Agroalimentaria y Pesquera. SIAP-SAGARPA. México, DF. 236 p.
Yang, J. S.; Wang, Y. Q. and August, P. V. 2015. Estimation of land surface temperature using spatial interpolation and satellite-derived surface emissivity. J. Environ. Informatics. 4(1):40-47.