MATERIALES Y MÉTODOS

El área de estudio se sitúa en el Complejo Valle del Río Dulce dentro de la Subregión del Chaco Semiárido de la República Argentina. Está comprendida dentro del área de riego del Río Dulce en la provincia de Santiago del Estero, a 27° 39' S, 64° 14' O (Morello, 2012) (Figura 1). El clima se caracteriza por ser subtropical con estación seca, con déficit hídrico y variabilidad climática dentro del ciclo anual y entre períodos anuales. El área en estudio presenta una temperatura promedio anual de 20,4 °C y precipitación promedio anual de 624 mm. Según datos de las estaciones meteorológicas situadas en el Aeropuerto Santiago del Estero y en el INTA La Abrita, ambas ubicadas en la zona de estudio, se corrobora la información presentada, en lo que respecta a precipitación promedio anual con valores entre 500 a 700 mm aproximadamente a partir de datos históricos (1961-2017), distribuida a través del año bajo un régimen estival. La estación lluviosa se extiende de octubre a marzo y el resto de los meses comprende la estación seca; se clasifica el área como subhúmeda-seca, la temperatura del mes más cálido es de 26,3 °C y la del mes más frío 11,9 °C (Morello, 2012).

El área donde se ubican actualmente los rodales de algarrobo blanco en la provincia corresponde a la Llanura aluvial del río Dulce, Paleo Llanura de Inundación, con material originario fluvio-loéssico, que posee gran heterogeneidad espacial en cuanto a las propiedades químicas de los suelos principalmente. Según la información registrada a través de GeoINTA (2017) a partir del mapa de suelos de la República Argentina basada en la clasificación taxonómica SoilTaxonomy (Soil Survey Staff, 2014), la mayor parte de las áreas con plantaciones se ubican en suelos MNai-7 y 10 (Haplustoles arídicos), los que se caracterizan por presentar perfiles de textura franca a franco-arenosa, con baja humedad, no sódicos. Algunos rodales se ubican en suelos Ektc-12 (Haplustoles salortídicos) de texturas arenosas con algún tipo de limitante de salinidad en los 50 cm de profundidad y generalmente ubicados en vías de escurrimiento plano cóncavas, y en AEtc-34 (Natracualfes típicos) de texturas franco-arenosas, de drenaje deficiente, de plano aluvial. El área de estudio en general se caracteriza por registrar valores de conductividad eléctrica que varían de 5 a 20 dSm^-1 y concentración de sodio del extracto de saturación en un rango de 21 a 137 meq/l. La presencia de la capa freática se encuentra cercana entre los 2 y 5 m de profundidad, aproximadamente.

Las áreas de influencia de este estudio comprenden bosques y vegetación dentro de la cuenca del Río Dulce, que adopta características distintivas con respecto a la subregión del semiárido. Se presentan comunidades riparias, con la presencia de los géneros y especies Aspidosperma quebracho-blanco Schltdl., Lycium spp., Schinopsis spp., Prosopis sericantha Gill. Ex Hook &Arn. y P alba Griseb. Se trata de un fachinal abierto dominado por Vachellia aroma ex Hook. & Arn. (tusca), Celtis iguanaea (Jacq.) Sarg. (tala), con árboles aislados de Schinopsis lorentzii Griseb (quebracho colorado) y Sarcomphalus mistol Griseb (mistol) (Sarmiento, 1963; Red Agroforestal Chaco Argentina [RedAF], 1999).

En este estudio el criterio de selección de las plantaciones fue una combinación de la edad y la variabilidad en los tipos de suelos, teniendo en cuenta principalmente la salinidad y las características propias de los rodales, buscando maximizar la variabilidad entre muestras y minimizar la variabilidad dentro de la unidad de muestreo. Cabe aclarar que se relevaron todas las plantaciones que superaron la edad mínima establecida de 9 años. Para evaluar la variabilidad se calculó el coeficiente de variación (CV) de la altura dominante (Hd), que resultó en un valor máximo de 25 % de variación dentro de la unidad muestral (Herrera, 1996). Se trabajó con un total de 26 parcelas ubicadas en las diferentes plantaciones del área de riego. Cada parcela posee un número fijo de 48 árboles y por consiguiente representan un área variable entre 576 a 1728 m², según el espaciamiento inicial de la plantación. A la fecha de plantación de estos rodales no se contaba con material genético mejorado, por lo cual el origen del germoplasma empleado provino en su totalidad de poblaciones locales. La base de datos a utilizar provino de 11 parcelas temporales (PT), o sea aquellas parcelas donde se ha medido la variable de crecimiento solo una vez; 15 parcelas permanentes (PP), o sea parcelas donde se ha medido la variable de crecimiento en medidas sucesivas en el tiempo, en algunos casos 2, 3 o 4 mediciones; y se han obtenido datos de crecimiento provenientes del análisis epidométrico de 11 árboles dominantes extraídos en dichas plantaciones. Se dispuso de 91 pares de datos altura-edad, incluyendo datos del análisis fustal de los árboles seleccionados para el análisis epidométrico. Las muestras para el análisis fustal se extrajeron de cada individuo seleccionando tortas de 5 cm de espesor a las siguientes alturas desde la base: 0,30 m, 1,3 m, 2,0 m y 3,0 m. Estas fueron preparadas en la forma descripta por Giménez et al., (1998) y luego digitalizadas con escáner HP Scanjet G2410, a una resolución de 1200 dpi para la medición de los anillos de crecimiento mediante el software IPWIN 4. Se trabajó sobre cuatro radios perpendiculares a partir de los que se consignó el número de años y el espesor del anillo con una precisión de centésima de milímetro.

Figura 1 Mapa de la provincia de Santiago del Estero, localizado en la ecorregión del Chaco Seco, República Argentina. Fuente: Morello et al. (2012) y Unidad Ejecutora del servicio de riego del río Dulce. Elaboración propia 

Para la caracterización del sitio se hizo un relevamiento preliminar del suelo en cada una de las parcelas. El estudio del perfil se realizó con calicatas situadas en el centro de cada parcela, extrayendo muestras de 0 a 90 cm de profundidad, discriminando por horizontes. Los análisis de las muestras de suelo se realizaron en el laboratorio del INSIMA, de la Facultad de Ciencias Forestales en la Universidad Nacional de Santiago del Estero (FCF-UNSE) y en el Laboratorio de Suelo, Agua y Fertilizantes (LabSA) de la Estación Experimental Agropecuaria Salta, Grupo Recursos Naturales. Los parámetros edáficos físicos determinados fueron densidad aparente (DA) (g/cm³) por el método del cilindro y porosidad de la capa arable (PTot) por el método de gravimetría de muestras volumétricas (100 cm³) (% v/v) (Klute, 1986). Los parámetros edáficos químicos se determinaron de acuerdo a las siguientes técnicas: pH por potenciometría, conductividad eléctrica (CE) estandarizada a 25 °C (dSm^-1), Sodio (Na+) por fotometría de llama y Calcio más Magnesio (Ca+Mg) por volumetría con EDTA (meq/l), todos sobre el extracto de saturación (Sparks, 1996) y presencia de carbonatos (%) mediante método descriptivo de reacción con HCl al 10 % según la guía de reconocimiento de suelos de la Organización de las Naciones Unidas para la Alimentación y la Agricultura (FAO, 2009). La presencia de carbonatos se categorizó en una escala ordinal donde 0 indicó la ausencia de carbonatos en todo el perfil; 1 representó la presencia de carbonatos en concentraciones entre 0-2 % en la profundidad del perfil a partir de los 70 cm de profundidad; 2 representó la presencia de carbonatos concentraciones entre 2 a 10 % también a partir de los 70 cm de profundidad, y 3 la presencia de carbonatos en todo el perfil con concentraciones entre 10 a 25 %. Se calculó la relación adsorción sodio (RAS) tomando la base general de datos con la siguiente formula RAS= Na+/raíz ((Ca+Mg)/2)).

En la evaluación del desarrollo de los rodales las variables medidas fueron: el diámetro a la altura de pecho (DAP) expresado en centímetros (cm) y medido con cinta diamétrica; la altura total del árbol (HT), expresada en metros (m) y medida con vara telescópica, y la edad (E) que fue expresada en años. En el análisis epidométrico, el espesor de anillos de crecimiento fue expresado en milímetros (mm). Para la determinación del índice de sitio se consideró la altura dominante (Hd) definida como la altura promedio de los 100 árboles más gruesos por hectárea (Assmann, 1970) expresada en metros (m).

Para la determinación de las curvas de altura dominante, se generó la serie edad-altura dominante construida a partir de registros dasométricos de las parcelas temporales y permanentes de medición y de los datos del análisis epidométrico de los fustes de árboles seleccionados. Se procedió según el método de curva guía que consistió en determinar una única función promedio a partir de la serie de datos altura dominante-edad y sobre la base de esta, derivar las demás curvas por un principio de proporcionalidad constante (Barrero-Medel et al. 2011; Hernández-Ramos et al., 2014).

Se utilizó el modelo matemático de Gompertz sugerido por Senilliani et al. (2019) como el más adecuado para la especie, un modelo de regresión no lineal mixto de Gompertz, con correlación por simetría compuesta y varianza de los errores con función de varianza exponencial. A partir de la curva guía y considerando la edad base determinada se trazaron las curvas de igual forma o curvas anamórficas, que pasaron por diferentes valores de IS; a este proceso se lo denomina armonización (Alder, 1980; Prodan et al, 1997). Se generó la familia de curvas anamórficas haciendo variar la edad y el IS deseado para la edad base (Eb) predeterminada (Andenmatten y Letourneau, 1998; Mares Arreola et al., 2004; Barrero-Medel et al., 2011; Hernández-Ramos et al., 2014). El procedimiento de armonización se realizó teniendo en cuenta la fórmula de Gompertz [Ecuación 1] para modelar la Hd en función de la edad base y en función de la definición de índice de sitio [Ecuación 2]. Las ecuaciones del proceso de armonización son las siguientes:

En resumen, la familia de curvas anamórficas, se obtuvo despejando a en la ecuación de índice de sitio y sustituyendo la expresión resultante en la ecuación original [1], con lo que el valor asintótico se consideró implícito y se asumió que los parámetros de forma son comunes para todos los sitios (Hernandez-Ramos et al., 2014). A partir de la dispersión de los datos de altura dominante obtenidos a la edad base se determina la equidistancia entre las curvas de IS. La validación del modelo de predicción del índice de sitio se realizó a través de remuestreos con reposición técnica denominada Bootstrap (Efron y Tibshirani, 1993). A partir de los remuestreos se estimó el promedio de la raíz cuadrada del error (RECM) calculado como los desvíos al cuadrado. En cada clase de sitio establecida se calculó el promedio de las alturas observadas para cada año de edad y con estos datos se estimaron los índices de sitio (ISest). Posteriormente se calculó el error respecto al índice de sitio para cada clase, (ISobs) (Barret, 1978; Fassola y Wabo, 1993 y Andenmatten y Letourneau, 1998).

Para realizar la caracterización de los factores edáficos que determinan la calidad de sitio se hizo una evaluación preliminar de la relación de la Hd y las variables edáficas relevadas en las parcelas. Se asignó a cada una de las parcelas de estudio su correspondiente clase de calidad a través del cálculo del IS, utilizando la altura dominante y la edad actual con la función del IS descripta en el ítem anterior. Se realizó un análisis de la varianza multivariado (MANOVA) para evaluar diferencias entre clases de calidad de sitio considerando las propiedades físicas y químicas. La variable Carbonatos fue analizada mediante tablas de contingencia según las categorías descriptas anteriormente; se trabajó con el estadístico Chi cuadrado para analizar las diferencias estadísticas entre clases de calidad.

La determinación de diferencias entre las clases de calidad de sitio para la variable IS y Hd se realizó mediante ANOVA. Si bien se analizaron las muestras de suelo de cada horizonte del perfil, se tomó el primer horizonte (A) de referencia para los análisis estadísticos por ser el que mejor permitió diferenciar calidades de sitio. Los datos fueron acondicionados, depurados y analizados estadísticamente con InfoStat (Di Rienzo et al., 2020) y su interface con el programa R (R CoreTeam, 2020).

RESULTADOS Y DISCUSION

Se registraron los valores promedios de altura dominante según clase de calidad e IS con diferencias estadísticamente significativas entre clases en la Hd (F = 19,5; p < 0,0001) e IS (F = 53,95; p < 0,0001). Los resultados del análisis de la varianza confirman que se encuentran correctamente representadas las tres calidades como resultado del crecimiento en altura de los rodales en estudio. Las muestras responden en 26,6 % a la mejor calidad de sitio (calidad I), 55,5 % a la calidad regular (calidad II) y 17,7 % a la menor calidad (calidad III).

Los valores obtenidos indicaron que la Hd se incrementa con la edad a una tasa mayor en la etapa inicial de la vida del rodal que en la madurez y se registraron mayores coeficientes de variación en los rodales de menor edad con respecto a los adultos (71 % vs. 13 %). El proceso de crecimiento se describe en tres etapas: una etapa inicial o juvenil del rodal, denominada temprana, que se estableció hasta los 8 años, donde hay una marcada aceleración del incremento de la Hd en el cual se observa un crecimiento de 0,8 a 1 m/año; luego una etapa de edades intermedias, denominada media, entre los 8 a los 15 años, donde el incremento disminuye respecto a las edades anteriores, observando valores de incremento desde 0,4 a 0,8 m/año; y una etapa posterior, denominada avanzada, en la cual el rodal alcanza la altura final de aproximadamente 10 a 12 m, en esta etapa se incluyen los rodales mayores a 15 años (Tabla 1).

La variabilidad de la Hd fue mayor en los rodales en la primera y segunda etapa de crecimiento (rango etario denominado temprano y medio) y, a medida que aumenta la edad del mismo, dicha variabilidad disminuye (DE = 1,6< 15 años y 1,1 > 15 años) (Tabla 1). En concordancia con lo expresado por Herrera y Alvarado (1998), la mayor variabilidad que se presenta a edades tempranas responde al comportamiento de las masas influenciadas por factores como la variabilidad genética, las técnicas de plantación y la densidad inicial, que pueden ser clasificados como factores genéticos y factores de manejo silvícola.

Curvas de crecimiento en altura dominante e índice de sitio

La muestra de árboles evaluados en este trabajo se encuentra entre los 9 a 21 años de edad y un rango de alturas entre 6 a 10 m (Figura 2), coincidiendo con las máximas alturas registradas para la especie en bosques nativos, en un rango de 10 a 12 m de altura en edades avanzadas de 90 a 100 años (Juárez de Galindez et al., 2005; Cisneros y Moglia, 2017). El patrón de comportamiento de la altura dominante (Hd) en función de la edad de los individuos no refleja un patrón polimórfico que indique que poseen diferente forma para cada calidad de sitio, demostrando una tendencia de carácter anamórfico (Figura 2). Estos perfiles anamórficos llevaron a la decisión de la construcción de curvas anamórficas a partir de la curva guía.

El criterio para definir la edad base habitualmente es arbitrario y elegido en función de la situación (Thrower, 1986). Sin embargo, el mismo autor expresa que en especies del este y oeste de Norte América usualmente se definió la edad de base con relación a la longevidad de la especie. Así se estableció en 100 años para el caso de especies de larga vida, 50 años en especies de corta longevidad y 25 años o menos para plantaciones jóvenes o especies manejadas con rotaciones cortas (Thrower, 1986). Otro criterio es el expresado por Ortega y Montero (1988) quienes establecen que la edad base se encuentra próxima a la mitad del turno de corta final. En forestaciones de P alba en la región del Chaco húmedo provincia de Chaco, Kees et al. (2017) proponen una edad base de 17 años y Pérez et al., (2011) para forestaciones con esta especie en la provincia de Formosa proponen como edad base 15 años. En este estudio la selección de la edad base en 15 años se considera adecuada, ya que según el registro de los datos, es la edad a la cual el potencial del sitio se ha expresado en el crecimiento en altura, coincidiendo aproximadamente con la mitad del turno de corta.

Figura 2: Patrón de comportamiento de la altura dominante en función de la edad en forestaciones de Prosopis alba, ubicadas en Santiago del Estero, Argentina. Línea continua: arboles evaluados a través de análisis de fuste; línea entrecortada: arboles evaluados en parcelas permanentes; puntos llenos: arboles relevados en parcelas temporales 

El modelo para obtener la curva guía y determinar el índice de sitio fue el modelo de regresión no lineal mixto de Gompertz con estructura de correlación simetría compuesta y función de varianza de los errores exponencial (Senilliani et al., 2019). El patrón de evolución en altura se ajusta satisfactoriamente al modelo Gompertz, adecuándose a los criterios biológicos de crecimiento. La aplicación de los modelos de efectos mixtos permitió modelar la estructura de los errores en datos provenientes de medidas repetidas y la presencia de heterocedasticidad. Los parámetros estimados en la función fueron estadísticamente significativos (p < 0,0001).

Tabla 1: Altura dominante (Hd) promedio y variabilidad (DE, mínimo y máximo) con relación a la edad de los rodales, definida por rangos etarios en edades tempranas (de 1 a 8 años), edades medias (de 8 a 15 años) y edades avanzadas (más de 15 años) en rodales de Prosopis alba en Santiago del Estero 

La función es la siguiente:

donde, Hd= Altura dominante (m) y E = edad (años). A partir de la ecuación [5] se graficó la curva guía que establece un IS de 8,2 metros a los 15 años (Figura 3a) y se pudo derivar la función que permitió expresar la familia de curvas anamórficas [6] que cubrían el espacio de crecimiento en altura dominante (Figura 3b).

Figura 3: 3a) Perfiles de la Hd en función de la edad medida en años para estimar la curva guía a partir de la ecuación de Hd en función de la edad base de 15 años. La curva de crecimiento fue contrastada con series longitudinales de parcelas permanentes (líneas segmentadas), análisis de fustes (línea sostenida) y datos de parcelas temporales (líneas de puntos). 3b) Familia de curvas anamórficas derivadas de la primera (línea continua), describiendo el índice de sitio y los límites de cada clase (línea entrecortada) 

La ecuación armonizada es:

A partir de los datos del estudio se generó la familia de curvas anamórficas mediante la ecuación [6], coincidiendo con el procedimiento empleado por Barrero et al. (2011) y Hernández-Ramos et al. (2014). A partir de la curva guía obtenidas en este trabajo se establecieron tres clases de calidad de sitio que fueron establecidas en función al rango de alturas encontradas a la edad base, con una dispersión de los datos de Hd entre 6,5 a 10,5 m. Los rangos de IS definidos cubren el espacio de crecimiento en altura dominante que tiene la especie. Los mismos se definieron de la siguiente forma: Clase I, correspondiente a la mejor calidad de sitio con 11,2 >IS> 9,2; Clase II, correspondiente a una calidad de sitio intermedia con 9,2 > IS > 7,2; y Clase III, que se corresponde con la calidad de sitio más baja con 7,2 > IS > 5,2 (Figura 3b). A partir de la dispersión de los datos de altura dominante descripta anteriormente se estableció la equidistancia entre las curvas de IS de 2 m entre clases, determinándose de esta manera los valores de altura en función de la edad para cada una de las curvas de IS (Tabla 2). En el manejo forestal de la especie P alba, conocer las clases de calidad en la zona de riego de la provincia es valioso para predecir o modelar la producción, como así también para hacer una correcta selección de los mejores sitios para plantación una vez determinados los factores del sitio que más influyan en el crecimiento.

Tabla 2: Valores predichos de altura dominante e índice de sitio según la clase de calidad de sitio (CI, CII y CIII) en función de la edad en plantaciones de Prosopis alba en la Región Chaqueña de la República Argentina 

Relación crecimiento-suelo

La capacidad productiva de la especie está influenciada por el sitio de crecimiento; así lo expresa la Hd. El resultado del análisis de la varianza multivariado en relación con las propiedades físicas en el horizonte A, PTot (porosidad de la capa arable) y DA (densidad aparente) no registraron variaciones estadísticamente significativas a través de las clases de sitio. Si bien la porosidad (PTot) es una variable de importancia entre las propiedades físicas, registra para el horizonte analizado un rango de 49,9 a 53,6 % sin diferencias estadísticamente significativas entre clases de calidad de sitio (F = 0,9 p = 0,4). Los valores de PTot obtenidos en este trabajo fueron menores a los reportados por Meza Pérez et al. (2003) para bosques templados con comunidades de Pinus y Quercus en México (60 a 62 %), similares a los reportados por González-Barrios et al. (2011) y Mellado (2012) para bosques húmedos de coníferas en España (50 %) y superiores a los encontrados por Acuña (2012) en suelos de plantaciones con coníferas en áreas subtropicales de Perú (39 a 40 %). Se considera relevante poder determinar en futuros estudios la proporción según el tipo de poros, ya sea macroporos, mesoporos o microporos con funciones importantes en el balance de agua y aire y en el desarrollo de organismos (Lorenz, 2005).

El análisis de la varianza multivariado con las propiedades químicas pH, CE, RAS y cationes solubles (Na+ y Ca+Mg) presentó diferencias estadísticamente significativas entre clases de calidad de sitio en base al estadístico de Wilks (0,5) (F = 2,8; p = 0,004). En la Tabla 3 se presentan los valores promedios para la Hd de los rodales y los respectivos valores de las variables edáficas evaluadas según la clase de sitio.

Las propiedades químicas a diferencia de las físicas, son determinantes, presentando diferencias estadísticamente significativas entre clases como se expresó anteriormente. Sin embargo, en comparaciones a posteriori entre clases, se observa que las clases I y II no difieren estadísticamente entre sí (Clase I vs Clase II, valor p = 0,67) mientras que sí difieren ambas con respecto a la clase III (Clase I vs Clase III, valor p < 0,001 y Clase II vs Clase III, valor p = 0,001). Estos resultados sugieren que si bien los cambios registrados en las propiedades químicas entre las clases I y II no son estadísticamente significativos, la Hd de los rodales estudiados difiere entre estas clases, siendo sus diferencias estadísticamente significativas. Esta diferencia registrada en la Hd pero no en las propiedades químicas podría deberse a la sensibilidad de la especie al crecimiento ante pequeños incrementos en la salinidad, resultando en una disminución de su crecimiento en altura. También sería posible considerar que los cambios en Hd que experimenta la especie de la clase I a la II se deban a otros factores no considerados en este estudio.

Considerando la clase I y II se propone como valores indicadores para la CE valores de 5,4 a 7 dSm^-1 al igual que la concentración de sodio de hasta aproximadamente 39 meq/l. Esta estimación se realiza considerando que estas clases, donde se registra el valor medio de IS de 10,2 y 8,2, son consideradas las clases de mayor producción maderera. El RAS no supera el valor de 8,8 para la clase I y II; sin embargo, en la clase III los valores llegan a 22. Esto indica una creciente relación de adsorción de sodio en la solución del suelo de la clase I a la III con una consecuente disminución de la capacidad potencial de crecimiento de los

Tabla 3: Valores promedios y desviación estándar (DE) para la variable dasométrica Hd y propiedades químicas y físicas del suelo por clase de calidad de sitio (la muestra por clase es igual a n= 12 en la clase I, n=25 en clase II y n=8 en clase III). Los valores corresponden al horizonte A en el perfil del suelo 

Letras iguales para cada fila simbolizan que las medias no son significativamente diferentes (p > 0,05) rodales. En la variable concentración de Ca+Mg, el comportamiento es similar; sin embargo, el pH no muestra diferencias entre clases (Tabla 3). En la variable presencia de carbonatos (p < 0,1) no se observaron valores que difieran entre las clases. Ante las características geomorfológicas de los suelos de la región es notable la relevancia que tienen las variables químicas del suelo en contraposición a otros sitios de estudio como la provincia del Chaco donde el 58 % de las causas de variación de Hd en plantaciones de P alba se relaciona a variables fisiográficas (relieve, drenaje) y edáficas (textura) mientras que las variables químicas del suelo no explicaron diferencias en la calidad de sitio (Kees et al., 2017).

A partir de los resultados encontrados en este trabajo, se puede inferir que el crecimiento está condicionado negativamente por parámetros característicos de las condiciones de salinidad en los suelos. Las propiedades químicas evaluadas en los suelos de las plantaciones de P alba en estudio demuestran ser un factor de importancia e impacto en la productividad de la especie.

Validación del modelo

El modelo ajustado para la altura dominante (Hd) indicó que en cada una de las clases de calidad de sitio y a edades tempranas se pueden presentar algunos valores extremos de diferencia de altura que pueden llegar a más de 1 m, pero en general a partir de los 10 años el error es menor, decreciendo el sesgo de la estimación a edades intermedias para las tres clases. La variabilidad en la estimación a edades medias para las clases presenta un rango de error de 0,2 a 1,1 y la tendencia es a sobreestimar en la clase I y a subestimar en la clase II. La clase III a partir de los 11 a 12 años muestra una tendencia a sobreestimar el IS dentro del rango de error nombrado. En general se puede concluir que el error estaría en un rango del 5 al 10% para la edad base de 15 años, estimada en este trabajo. En concordancia con estudios realizados sobre otras especies, diferentes autores (Amico et al., 2010; Ribeiro et al., 2016) evaluaron la precisión de los modelos de índice de sitio como satisfactoria con valores similares que no superan el 10% RECM. Se visualizaron diferencias en la estimación de las alturas, enunciadas anteriormente, de hasta 1 m y un alejamiento máximo de un valor observado en relación a la curva de clase I de 2 m a la edad de 9 años (Figura 4).

Figura 4: Variabilidad de la estimación entre los valores estimados (curva guía y curvas correspondientes a clase II y III) y los valores observados (rombos, triángulos y asteriscos) de las alturas dominantes promedio en relación a la edad, correspondientes a cada clase de calidad de sitio (C I, C II y CIII) respectivamente. 

CONCLUSIONES

El patrón de evolución en altura muestra que la especie culmina su crecimiento a edades tempranas en plantaciones del área de estudio, cercanas a los 20 años, si bien esto deberá corroborarse con estudios a edades más avanzadas. Una primera aproximación en la estimación de la capacidad productiva de los sitios para P alba Griseb se expresa a partir de la familia de curvas anamórficas a la edad base de 15 años, representativas de tres calidades de sitio, I (buena), II (regular) y III (mala). El proceso de validación determinó que el error estaría en un rango del 5 y 10 % para la edad base (15 años). Los factores del sitio que son determinantes en el comportamiento de las clases de calidad son los factores edáficos a través de las propiedades químicas del suelo que a su vez describen las condiciones de salinidad.

La caracterización de los sitios de plantación de P. alba y la estimación del índice de sitio en el área de riego de la provincia es una herramienta importante para poder realizar una correcta selección de un sitio para plantación. Ante lo cual se recomienda para algarrobo blanco, tener en consideración los valores indicadores de CE y concentración de cationes solubles de sodio, los cuales no deberían superar los 5 a 7 dS m^-1 y 39 meq/l respectivamente. Se considera necesario profundizar el estudio con la ampliación de la base de datos, aumentar la muestra de estudio para la evaluación potencial de los sitios incorporando mayor número de mediciones en parcelas permanentes de diferentes edades y sitios e incorporar, si se justificara, el análisis de patrones polimórficos, que podría conferir mayor precisión en la estimación de las calidades de sitio.

Para avanzar en una adecuada zonificación de la actividad de plantación con P alba se recomiendan las calidades de sitio I y II para forestación con fines de obtención de servicios de provisión (madera y frutos) y la calidad III, debido a su baja productividad, podría destinarse a la producción de frutos y servicios ecosistémicos de regulación, entre ellos la recuperación de áreas degradadas.