La enfermedad conocida como Huanglongbing o HLB, es una plaga cuarentenaria ausente de la mayoría de la regiones citrícolas de Sudamérica, a excepción de Brasil, por lo que representa una de las mayores amenazas para la citricultura de la región, ya que es considerada la enfermedad más seria de los cítricos a nivel mundial (Halbert y Manjunath, 2004). Diversos autores estiman que las plagas agrícolas serían capaces de responder a los cambios climáticos mejor que los cultivos agrícolas que afectan y las comunidades de controladores biológicos nativos, por los que el cambio climático podría potencialmente generar condiciones propicias para la ocurrencia de un mayor número de invasiones biológicas, así como un incremento de la severidad de los daños ocasionados por las plagas agrícolas y/o forestales (Aurambout et al., 2006).

Dado que los factores climáticos son críticos para el desarrollo de cítricos y de los insectos vectores del HLB (Hall et al., 2008), el objetivo del trabajo fue estimar los efectos del cambio climático en la distribución geográfica potencial de HLB en Sudamérica.

El análisis de riesgo climático de HLB en Sudamérica, se realizó mediante el empleo de una técnica de comparación climática basada en el Análisis de Componentes Principales (PCA), descrita por Robertson et al. (2001). Esta técnica predice la aptitud ambiental para el establecimiento de una especie, a partir del empleo de variables ambientales predictoras, por medio de la construcción de un hiperespacio para la especie objetivo. Las variables ambientales predictoras utilizadas fueron a) temperatura máxima media mensual, b) temperatura mínima media mensual y c) precipitación mensual promedio, considerando escenarios históricos del Worldclim (1950-2000) y del Intergovernmental Panel on Climate Change (IPCC) vía de concentración representativa 2,6, modelo HadGEM2-ES para las décadas de 2050 y 2070, con una resolución de 1 km (Hijmans et al. 2005).

Como primer paso del análisis propuesto por Robertson et al. (2001), se obtuvieron los valores de las variables climáticas asociados al Estado de San Pablo y algunos municipios de los Estados de Paraná y Minas Gerais (Brasil), a partir de 450 puntos elegidos al azar, en celdas que tuvieran un NDVI mayor a 0,3 en todos los mese del año (Lopes et al., 2009).

Este set de datos fue considerado como valores de entrenamiento. Se generó una matriz de valores estandarizados, denominada matriz (U), cuyos valores medios y desvíos estándar fueron conservados para el tercer paso del análisis (ver figura 1). Posteriormente se realizó un PCA de la matriz U, dando como resultado la matriz (V). En el tercer paso, las variables de predicción fueron estandarizadas por medio de la media y el desvío estándar calculados a partir de los datos de entrenamiento, generando la Matriz (W). Esta matriz fue multiplicada por la matriz V, para producir la matriz (Z) de puntuación de los componentes para todo Sudamérica. A continuación, se seleccionaron aquellos componentes cuyos eigen valores excedieran por la media más dos desvíos estándar de las estimaciones. A continuación, las varianzas de cada componente fueron divididas por sus respectivos eigen valores, para finalmente calcular la probabilidad de cada píxel, al sustituir por una función de probabilidades de Chi-cuadrado (α=0,1), la suma de los cuadrados de la matriz Z estandarizadas. Se tomaron estos valores de probabilidad para indicar la idoneidad ambiental para la especie objetivo, en cada píxel del área bajo estudio.

Figura 1. Figura 1. Esquema general de implementación de análisis de componentes principales modificado de Robertson et al. (2001).

Para estimar los cambios en la susceptibilidad ambiental para el establecimiento de HLB en nuevas áreas, ante escenarios futuros de cambio climático, se procedió a estandarizar las variables de predicción pertenecientes los escenarios futuros estimados por el IPCC para las décadas de 2050 y 2070, por medio de la media y el desvió estándar calculadas a partir de los datos de entrenamiento calculados en base a las estadísticas climáticas históricas de Sudamérica. Se calculó el estadístico kappa de las capas resultantes contra una matriz de sitios con presencia de HLB en Brasil. Para el tratamiento de la información se utilizaron los softwares QGIS 2.12 (Quantum GIS Development Team, 2015), IDRISI Selva e Infostat (Di Rienzo et al., 2013).

Considerando las estadísticas climáticas históricas mensuales de Sudamérica, las áreas de similar riesgo climático a los valores de referencia (Estado de San Pablo y algunos municipios de los Estados de Paraná y Minas Gerais, Brasil), representarían el 19,32 % de la superficie de Sudamérica, distribuidas entre el sur de Brasil, el sudeste de Paraguay y Noreste de Argentina. También se puede destacar el centro de Chile, Bolivia y Perú (ver Figura 2).

Figura 2. Figura 2. Mapa de probabilidad de establecimiento de HLB en Sudamérica, en comparación al riesgo climático observado en San Pablo (Brasil).

Considerando el escenario climático estimado para Sudamérica por IPCC5 para la década del 2050, en base al escenario RCP 2,6 y el modelo HadGEM2-ES, las áreas de riesgo climático para el establecimiento de HLB representarían el 22,66 % de la superficie del subcontinente. Estas áreas se desplazarían hacia el sudoeste del área climáticamente riesgosa en la actualidad. Por ende se incluiría a Uruguay, en Argentina el noreste de la región pampeana y el área oriental de la región del Noroeste Argentino. Se incrementarían las superficies de riesgo en el centro de Chile, Bolivia y Perú. Se reducirían las áreas de riesgo en Brasil. Las áreas de riesgo climático relativo de establecimiento de HLB hacia la década del 2070, representarían el 19,98 % de la superficie del sub-continente y se desplazarían aún más hacia el sudoeste del área considerada climáticamente predisponente en el 2050. Incluiría una mayor superficie del noroeste de la Región Pampeana y el oeste de la Región Chaqueña. En Chile se reduciría el área afectada en comparación con la década del 2050 y se desplazarían más al sur. Se mantendría en categorización de riesgo climático en el área oriental de la Región del Noroeste Argentino, en Uruguay y Paraguay. Se reducirían las superficies de riesgo en el centro de Bolivia y Brasil. En base a los criterios de Viera y Garrett (2005), el modelo logró un ajuste sustancial al considerar las estadísticas históricas Kappa: 0,69; un ajuste moderado para la década del 2050, Kappa: 0,44 y un valores de concordancia justo para la década de 2070, Kappa: 0,36.

Los resultados hallados en este trabajo coinciden con la distribución potencial de HLB obtenida por Narouei-Khandan et al. (2015), mediante el algoritmo de máxima entropía (MAXENT), pero difieren significativamente del estimado mediante la utilización de Máquinas de soporte de vectores (SVM). Por tal motivo el modelo de consenso de ambos algoritmos, difiere de todas las estimaciones halladas en este trabajo.

Una la fortaleza de este trabajo radica en que las variables ambientales procesadas mediante PCA, resultaron en la incorporación a los modelos de componentes principales que constituyen las proyecciones ortogonales de las variables transformadas y, en consecuencia, libres de autocorrelación (Al-Kandari y Jolliffe, 2001). Por otro lado, los modelos de distribución generados con predictores ambientales en bruto generalmente predicen una superficie de ocupación más pequeña, es decir se encuentran sobreajustados (Phillips et al., 2006)

Aurambout et al. (2009), demostraron al menos teóricamente, que las variaciones estacionales de temperatura esperadas ante diversos escenarios de cambio climático, afectarán la velocidad de desarrollo y la supervivencia de D. citri, en sincronía con la de una especie cítrica hospedera, la naranja Valencia. Sin embargo, argumentan que el calentamiento global tendrá repercusiones complejas y espacialmente heterogéneas sobre la dinámica poblacional y la supervivencia de D. citri y por ende sobre el potencial establecimiento y dispersión del HLB (Figura 3) . El cambio climático tendrá un impacto significativo sobre la fecha de ocurrencia de la brotación de los cítricos y sobre la duración de las mismas. En este último caso, al acortarse el tiempo necesario para que los tejidos blandos de los brotes de los cítricos se endurezcan, se reducirá el tiempo disponible para que D. citri se reproduzca en los mismos, durante esta fase tan crítica para el desarrollo poblacional de las ninfas de esta especie (Pearson y Dawson, 2004).

Figura 3. Figura 3. Mapa de riesgo relativo de establecimiento de HLB en Sudamérica hacia el 2050 (izq.) y 2070 (der.), en comparación al riesgo climático históricamente observado en San Pablo (Brasil).

La naturaleza de las interacciones existentes entre el hospedero vegetal, el HLB, D. citri y el clima es muy compleja y heterogénea, no es el objetivo de este trabajo evaluar un modelo que refleje con precisión los mecanismos inherentes a un proceso espacio-temporal epidemiológicamente tan complejo. Por lo que las evaluaciones propuestas en este trabajo deben considerarse como un acercamiento útil y mejorable y como una herramienta aplicable a la estimación de los potenciales impactos del cambio climático en la distribución de las especies.