Características del sitio experimental

El estudio se realizó en el Campo Experimental Valle de México del INIFAP, ubicado, en la localidad de Santa Lucía, Coatlinchán, Estado de México. Este sitio se caracteriza por su clima templado frío con lluvias en verano, como se puede observar en el Cuadro 1, donde también se indican la localización geográfica y las características climáticas y edáficas.

Determinación de COS

La determinación del COS se llevó a cabo antes del establecimiento del experimento en el mes de junio de 2016. Se tomaron seis muestras de suelo distribuidas en la superficie experimental de 240 m², en la profundidad de suelo (Pm) de 0-30 cm con el siguiente procedimiento: las muestras de suelo se secaron al aire, se registró su peso; con auxilio de pinzas se retiraron raíces, restos vegetales y animales visibles; se tamizaron a 2 mm; y se tomaron submuestras. La determinación de la concentración de COS se realizó con el método de ^{Walkley y Black (1934)} que reporta la materia orgánica del suelo (MOS) y se utilizó el factor 0.58 (1/1.724) para convertir MOS a COS.

Para el cálculo de COS (Mg ha^-1 año^-1), se obtuvo la densidad aparente (Da) (g cm^-3) por el método del cilindro. La Da fue la relación entre el peso del suelo seco y el volumen del suelo: el peso de suelo mineral más el humificado, sin considerar el peso de las raíces, residuos y piedras y el volumen del suelo, se calculó al descontar del volumen total, el volumen de raíces, residuos y piedras. La cantidad de COS (Mg ha^-1) fue el producto de Pm, Da y el% COS.

Quinua

Al considerar que cerca de 80% de la MST de quinua queda en el terreno de cultivo después de la cosecha, se evaluó la acumulación de carbono (C) en la materia seca por estructura vegetal y total en tres en tratamientos de dosis de fertilización y un testigo. Según ^{Bazile et al. (2014)}, en la cosecha en forma manual o mecánica, después de cortar las panojas o inflorescencias, estas se dejan en el campo agrupadas en montones para que se sequen completamente.

Posteriormente se procede a la trilla, sobre una lona se colocan las panojas, con ayuda de un implemento agrícola o golpeándolas con una vara se desprende el grano, de esta manera, la mayor parte de los residuos de cosecha quedan en los terrenos de cultivo. Los tratamientos de dosis de fertilización fueron los siguientes: T1= 40-40-00; T2= 80-40-00; T3= 120-40-00 y T0= 00-00-00. Las unidades experimentales fueron de cinco surcos de 0.8 m de ancho por 5 m de largo (20 m²).

Estos tratamientos se evaluaron en un diseño experimental en bloques al azar con tres repeticiones en arreglo de parcelas divididas. En las parcelas se establecieron las variedades de quinua Amarilla Maranganí y Blanca y en los bloques los tratamientos de dosis de fertilización, en total hubo 24 unidades experimentales (240 m²). Además de los RV, el aporte de C al suelo se dio por el abono que se colocó sobre el suelo junto con la semilla al momento de la siembra.

La materia seca se midió cuando el cultivo alcanzó la madurez fisiológica a los 147 días después de la siembra y consistió de las siguientes actividades: se cortaron cuatro plantas al ras del suelo, se separaron las estructuras vegetales: hojas, tallos, inflorescencias y grano, se obtuvo el peso fresco total y el peso en seco de una submuestra de 100-200 g, la submuestra fue puesta a secar en una estufa a una temperatura de 65 ^oC. Con estos últimos datos se obtuvo la materia seca por órgano al substraer el porcentaje de humedad del peso de la biomasa total en fresco.

El porcentaje de humedad fue el peso húmedo menos el peso seco entre el peso húmedo de las submuestras. La biomasa de las raíces fue estimada como un 10% de la biomasa aérea total, al considerar que en otros cereales como trigo y cebada es 10 a 15% según el estudio de ^{Kuzyakov y Domanski (2000)}. El carbono por órgano fue el producto de la materia seca (g) y la concentración de carbono promedio de tres repeticiones en el tejido vegetal (%). La concentración de carbono en el tejido vegetal se determinó por combustión seca en un analizador automático de carbono orgánico total (Shimadzu TOC 5000-A).

Variables de crecimiento y rendimiento

Para comparar el crecimiento y rendimiento de las variedades de quinua en estudio, en la etapa de madurez fisiológica se midieron las siguientes variables: altura de planta, medida desde la base del tallo a la punta de la panoja principal, diámetro de tallo, medida a 10 cm de la base de la planta, longitud de panoja, materia seca total, rendimiento de grano por planta y por hectárea e índice de cosecha (IC).

Maíz

En la estimación de los RV que aporta el maíz al suelo, se consideró la información para la zona de mediana productividad bajo temporal estricto de Texcoco que es similar a las condiciones bajo las cuales se evaluó el cultivo de quinua. Se consideró un rendimiento de 3 Mg ha^-1 según la ^{SAGARPA (2015)} y un índice de cosecha de 0.34. Esta información dio una entrada anual de RV de 3.5 Mg C ha^-1 año^-1 que corresponde a 40% de la materia seca total, 10% lo aportan las raíces de acuerdo con las estimaciones de ^{Kuzyakov y Domanski (2000)} para cereales y 30% los tallos de maíz (^{Pérez et al., 2000}).

Modelo RothC

El modelo RothC de ^{Coleman y Jenkinson (2005)} divide el COS en compartimentos que contienen materiales con diferentes tasas de descomposición, cuatro son activos y uno es pasivo: (i) material vegetal de fácil descomposición (MVF); (ii) material vegetal resistente (MVR); (iii) biomasa microbiana (BIO); (iv) materia orgánica humificada (HUM); y (v) materia orgánica inerte (MOI).

Los compartimentos activos experimentan descomposición por cinética de primer orden, según la siguiente expresión: Y = Y₀ (1 ‑ e^-abckt). Donde: Y₀ es el C inicial del compartimento activo; k es la tasa constante de descomposición anual; t es 1/12 para obtener la tasa de descomposición al final de cada mes; y a, b y c son los factores que modifican a k y son temperatura, humedad y cobertura del suelo, respectivamente. El valor de k para cada compartimento presenta los siguientes valores: MVD (10), MVR (0.3), BIO (0.66) y HUM (0.2).

El compartimento pasivo, MOI (Mg ha^‑1) se obtuvo con la ecuación de ^{Falloon et al. (1998)} expresada por MOI= 0.049 × COT^1.139, donde COT es el carbono orgánico (Mg ha^‑1). Los datos de entrada del modelo RothC son: 1) factores climáticos: promedio mensual de temperatura del aire (^oC), precipitación (mm) y evaporación (mm), TPM, PPM y EPM, respectivamente; 2) edáficos: contenido de COS (Mg ha^-1), contenido de arcilla (%) y profundidad de muestreo del suelo; 3) entrada mensual de C al suelo de residuos vegetales (RV) y abonos orgánicos (AO); 4) cobertura vegetal mensual del suelo; es decir, si el suelo está desnudo o con cobertura vegetal; 5) la relación MVD/MVR, que de acuerdo con ^{Coleman y Jenkinson (2005)} fue de 1.44 (59% son para MVD y 41% son para MVR); y 6) el valor de la MOI obtenido de la ecuación de ^{Falloon et al. (1998)}.

Simulaciones con el modelo RothC

La simulación de los cambios de COS incluyó la inicialización y elaboración de escenarios. En la inicialización, se obtuvo el contenido de C inicial de los compartimentos activos en la condición del suelo en equilibrio. Esto se logró ejecutando el RothC iterativamente 10 000 años con la información de clima, suelo, relación MVD/MVR, MOI y entrada de C de los RV al suelo obtenidos por el modelo. Los escenarios, se ejecutaron con el valor de C de los compartimentos activos y pasivo, que se obtuvieron en la inicialización y con la información de clima, suelo, MVD/MVR, cobertura del suelo y MOI.

Las simulaciones de los escenarios de la dinámica del COS con el RothC se efectuaron con la información: tres periodos de tiempo: 20, 60 y 100 años; los sistemas anuales: monocultivo de quinua (MQ), monocultivo de maíz (MM), rotación de quinua-maíz (RQM) y el uso de aporte de C al suelo proveniente de los RV: 60, 70 y 80% de la MST de quinua y 40% de la MST de maíz.

El aporte de C por hectárea (Mg ha^-1 año^-1) en quinua se calculó como el producto de la proporción de la MST (%) y la densidad de población (plantas ha^-1): 129, 166 y 136, 458 para Amarilla Maranganí y Blanca, respectivamente. Para determinar la tasa de cambio (TC) del COS en los periodos de tiempo 20, 40 y 60 años, se usó la siguiente expresión: TC COS_{20 años}= COS_{año 2000} - COS_{año 2020}. Un resumen de la información de entrada al modelo asumida y medida del sitio de estudio se muestra en los Cuadros 1, 2 y 3.

Análisis estadístico de las variables de estudio

La acumulación de C en la materia seca por estructura vegetal y total, y las variables de crecimiento y de rendimiento se sometieron a un análisis de varianza con el procedimiento LGM de SAS. Para comparación de medias se usó la prueba de Tukey (p= 0.05).

Distribución de carbono en la biomasa

Los resultados de la comparación de medias de las variables de crecimiento y rendimiento para las quinuas Amarilla Maranganí y Blanca (Cuadro 4), indicaron diferencias estadísticas para altura de planta, diámetro de tallo, materia seca total y rendimiento de grano por hectárea, en todos los casos la Amarilla Maranganí fue superior a la Blanca.

El rendimiento por hectárea en Amarilla Mananganí fue cercano al promedio nacional de 2.8 Mg ha^-1 y en ambas fue alto en el tratamiento de fertilización T3 (120-40-00), aunque no presentó diferencias significativas entre tratamientos de fertilización. La longitud de panoja estuvo dentro del rango de 30 a 80 cm, reportado por ^{Mujica et. al. (2003)} para las quinuas peruanas. El IC obtenido en el presente estudio fue menor al rango promedio reportado en otros estudios, ^{Mujica et al. (2003)} quien reportó un rango en quinua de 0.24 a 0.45.

El rendimiento (kg ha^-1) del presente estudio fue superior al promedio reportado experimentalmente por ^{Mujica et al. (2003)} en las mismas variedades cultivadas de Perú (0.93 kg ha^-1) y también fue superior para el caso de Amarilla cultivada en Viacollo, Chile donde presentó una producción de 1.35 kg ha^-1 (^{Delatorrre et al., 2008}).

Aporte de carbono vegetal al suelo

Las medias de la acumulación de C en hoja, inflorescencia, y grano, MST e IC fueron estadísticamente iguales en ambas quinuas, la excepción se presentó en los tallos, con mayor valor en Amarilla Maranganí (Cuadro 5).

El Cuadro 6, muestra el aporte de C al suelo por las quinuas y el maíz. La distribución de C por órgano: tallo, hoja, inflorescencia, grano y raíz, en los tratamientos de fertilización en quinua no presentó diferencias significativas a pesar de esto, el T3 (120-40-00) tuvo la mayor cantidad de C por órgano y MST total.

Con base en estos resultados, para efectuar las simulaciones con el RothC, se usó el valor promedio de los tratamientos (g de C planta^-1) y con este, se obtuvieron los porcentajes de entrada C de RV de 60, 70 y 80% del a MST, estos últimos, fueron cercanos a los reportados por ^{González et al. (2017)} para maíz bajo labranza de conservación en la región de Bajío (2.6 a 5.9 Mg de C ha^-1 año^-1).

En la comparación de medias de las variables de crecimiento y rendimiento por variedad de quinua en los tratamientos de fertilización, no hubo diferencias significativas (Cuadro 5), lo que se explica por el nivel de fertilidad favorable del suelo donde la quinua se estableció, con un contenido de materia orgánica del suelo de 2.3% y de nitrógeno mineral de 91 ppm, ya que la quinua es un cultivo que mantiene su productividad aun en suelos con baja fertilidad, poca humedad e incluso alta salinidad (^{Ruiz et al., 2014}). Sobre la distribución del C por órgano en los tratamientos de fertilización (Figura 1).

Figura 1 Distribución de carbono por estructura vegetal en los tratamientos de fertilización en madurez fisiológica en las quinuas; a) Amarilla Maranganí y b) Blanca. T1: 40-40-00; T2: 80-40-00; T3: 120-40-00; y T0: 00-00-00.  

Se encontró que en ambas quinuas los tallos y las inflorescencias son cerca de un 80% de la MST. En Amarilla Maranganí es mayor el C en los tallos, 41-46% de la MST, mientras que en Blanca es mayor en las inflorescencias, 44-46% de MST. En ambas quinuas, las hojas y grano representan menor proporción de la MST, de un 0.8 a 2.8% y de 16 a 22% de la MST, respectivamente. Esto dio como resultado, bajo IC (0.17-0.18) que se refleja en alta producción de MST en hojas, tallos e inflorescencias y baja producción de MST en grano. Resultados que estuvieron cercanos con lo reportado por ^{Mujica et al. (2003)} para quinua, con 45% para los tallos, 29% para las hojas e inflorescencias y 26% para el grano.

Simulación de escenarios de COS con diferentes aportes de RV

Los resultados de simulación de escenarios con diferentes aportes C y en diferentes periodos de tiempo se presentan en el Cuadro 7 y Figura 2. En los periodos de tiempo evaluados, hubo secuestro de COS solo en MQ con el aporte de C al suelo de RV de 80% de MST (0.02-0.1 Mg C ha^-1año^-1), mientras que, el orden de pérdidas de COS de menor a mayor fue el siguiente: MQ con RV de 70% MTS< MQ con RV de 60% MST (entre 0 a -0.16)< RQM (-0.02 a -0.18); y MM (entre -0.11 a -0.19).

Figura 2 Escenarios de cambios de carbono orgánico del suelo con diferentes aportes de carbono al suelo en las quinuas: a) Blanca y b) Amarilla Maranganí. 

Esto puede explicarse por el tiempo de conversión o cambio del COS, que de acuerdo con ^{Jenkinson y Rayner (1977)}, es el movimiento del C orgánico, a través de un volumen dado de suelo, y se obtiene de la relación entre el COS y la entrada anual de RV. El tiempo de conversión promedio (años) para los diferentes periodos de simulación con el RothC en los sistemas estudiados tuvo el orden inverso al COS secuestrado: MM (9.3)> RQM (8.5)> MQ con 40% de MST (8.9)> MQ con 60% de MST (7.7)> MQ con 80% de MST (6.7).

Los valores más altos de conversión de COS: MM y RQM, indicaron la estabilización del C en el suelo, y los más bajos: MQ con 80% de MST, indicaron que hubo una migración más rápida del C y por lo tanto, indicaron que el suelo está más lejano del estado de equilibrio por una mayor entrada de RV (Figura 1, Cuadro 4). El COS almacenado o secuestrado estuvo dentro de lo reportado a escala mundial por ^{Lal (2004)}, quien reporta de 0.02 a 0.76 Mg ha^-1 año^-1 para el caso de sistemas agrícolas que adoptan sistemas de manejo mejorado.

En general, la quinua Amarilla Maranganí en monocultivo y bajo rotación con maíz superó las tasas de cambio de COS de la Blanca en monocultivo y de monocultivo de maíz con las entradas de C al suelo evaluadas con el RorhC, resultados que tuvieron una relación lineal con su mejor crecimiento, rendimiento y entradas de C al suelo con relación a Blanca. Las entradas de C al suelo por los RV evaluados tuvieron relación lineal con las tasas de secuestro (Cuadro 5 y Cuadro 6).

Similarmente a lo encontrado por ^{Wang et al. (2017)} quienes, a nivel global en cereales con diferentes entradas de C al suelo, tuvieron una correlación alta con los cambios del COS. Particularmente, las tasas de cambio de COS encontradas en quinua de nuestro trabajo respondieron a las condiciones particulares de la productividad de la materia seca influenciada por el clima y retención de residuos, aplicación de abonos y las salidas de MST, como reporta ^{Smith et al. (2008)}; ^{Luo et al. (2010)} en los sistemas agrícolas. Las ejecuciones de escenarios con el RothC detectaron que ambas quinuas responden a las necesidades actuales del uso de cultivos con potencial de secuestro de COS y para mantener el suelo en una condición de equilibrio entre las ganancias y pérdidas de este elemento (Figura 2).