INTRODUCCIÓN

La cebada ( Hordeum vulgare L.) es el quinto cereal de mayor producción a nivel mundial ( ^{INEC, 2010} ); con el 50% del área y 63% del volumen de producción concentrados en Europa, donde se produce noventa millones de t/año ( ^{Abbassian, 2010} ), con una productividad promedio de 4,00 t/ha ( ^{INEC, 2010} ). A pesar de que hay una tendencia leve en la reducción de la demanda mundial de este cereal (Abbassian, 2010), ya sea por limitaciones agronómicas y económicas que controlan el mercado o por el repunte en la demanda de arroz ( Oryza sativa L.) y maíz ( Zea mays L.) ( ^{Castro et al., 2011} ), se mantiene como un insumo importante para la industria alimentaria, en especial para la industria cervecera.

En Ecuador, las condiciones agroclimáticas para la producción de cebada incluyen zonas de 2400 - 3300 msnm, precipitaciones de 400 - 600 mm durante el ciclo de cultivo, suelos franco arenoso y profundos con buen drenaje, y con un pH de 6,5 - 7,5 ( ^{Coronel y Jiménez, 2011} ). El 40% de la producción ecuatoriana se usa para producir cerveza, mientras que los excedentes se comercializan en mercados locales y sirven para generar subproductos para la alimentación animal y humana ( ^{Garófalo et al., 2010} ). A pesar de la importancia de la cebada para la economía agrícola, en Ecuador se producen solo 24 000 t/año, con una productividad promedio de 0,60 t/ha, y con costos de producción de hasta US$ 700 por hectárea ( ^{INEC, 2010} ). Es por ello que el país importa hasta 40 000 t/año, por un valor superior a US$ 10 millones, para suplir la demanda de la industria cervecera ( ^{BCE, 2014} ).

El N es el elemento más limitante para la producción de cereales ( ^{Arnall et al., 2009} ; ^{Scharf, 2015} ), porque contribuye a la acumulación de biomasa, forma la molécula de clorofila y participa del balance nutricional de la planta al facilitar la absorción y asimilación de otros nutrientes como el K+ ( ^{López et al., 2009} ). En la cebada, la disponibilidad de N es importante, sobre todo en las etapas de mayor demanda que son el ahijado y el encañado ( ^{López-Bellido, 2013} ). Una adecuada disponibilidad de N puede incrementar no solo la calidad del grano, en especial el contenido de proteína, sino también el nivel de producción hasta un 45% ( ^{Reussi-Calvo et al., 2006} ; ^{Ferraris et al., 2008} ; ^{Blackshaw et al., 2011} ; ^{Morojele y Kilian, 2015} ). Aunque para la industria cervecera se prefieren variedades de cebada con contenidos de proteína menores al 14,00% y alto contenido de almidón ( ^{CN y MAGAP, s.f.} ).

La cebada requiere 25 kg N/t de grano producido ( ^{García-Serrano et al., 2009} ), pero ^{Mahler y Guy (2005)} reportaron requerimientos de hasta 35 kg N/t de grano. Sin embargo, con un manejo intensivo y alta producción, la cebada absorbe 35 kg/ha solo después de la floración ( ^{HGCA y SEERAD, 2006} ), por lo que el requerimiento total puede ser mayor a este valor. ^{Moreno et al. (2003)} indicaron que una producción óptima de cebada requiere 120 kg N/ha como fertilizante, y que aplicaciones superiores a 150 kg N/ha resultan en excesivo costo de producción y disminución en la productividad. En cambio, si se desea maximizar la rentabilidad del cultivo, ^{Jankovic et al. (2011)} sugirieron fertilizaciones entre 50 y 110 kg N/ha, y ^{Hajighasemi et al. (2016)} indicaron que solo se debe aplicar 100 kg N/ha.

Debido a que la mayoría de suelos de Ecuador no proporcionan todo el N que el cultivo requiere, debido a la lenta liberación de N a partir de la materia orgánica, se necesita una adecuada fertilización nitrogenada (FN) para corregir y completar la nutrición del cultivo ( ^{López et al., 2009} ). Para una producción de 3,00 - 4,00 t/ha de grano de cebada se deberían aplicar entre 45 - 60 kg N/ha ( ^{Falconí et al., 2010} ), aunque las compañías cerveceras que conforman el consorcio de Cervecería Nacional y el Ministerio de Agricultura, Ganadería, Acuacultura y Pesca (MAGAP) recomiendan aplicaciones que varían de 0 (cuando se rota con el cultivo de leguminosas) a 90 kg N/ha (en suelos con bajo contenido de materia orgánica) ( ^{CN y MAGAP, s.f.} ). Las variaciones con respecto a la dosis de FN, se deben a que la extracción de N por parte del cultivo depende de la rotación de cultivos, tipo de suelo y contenido de materia orgánica, y a la humedad disponible para transportar el N en el sistema suelo-planta ( ^{Garófalo y Calvache, 2012} ; ^{Stevens et al., 2015} ). Sin embargo, aún con una adecuada fertilización, el país no alcanza la productividad de latitudes templadas donde se reporta producciones de 4,00 - 7,00 t/ha ( ^{Moreno et al., 2003} ; ^{Halvorson y Reule, 2007} ; ^{Malecka y Blechareczyk, 2008} ; ^{Alley et al., 2009} ; ^{Blackshaw et al., 2011} ).

La urea es el fertilizante nitrogenado de mayor uso a nivel mundial ( ^{Aguiar y Herrera, 2011} ; ^{Frame et al., 2014} ), y en Ecuador es la fuente de N más utilizada para la producción de cebada; sin embargo, no hay una industria productora de urea que pueda abastecer el mercado local, por lo que este fertilizante se importa por un valor cercano a los US$ 200 millones al año ( ^{BCE, 2015} ).

Cuando la urea entra en contacto con el suelo, su humedad se hidroliza (amonifica) para convertirse en NH₄ ⁺; luego el NH₄ ⁺ se transforma en NO₂ ^- y NO₃ ^- (nitrificación) ( ^{Tenuta y Beauchamp, 2000} ; ^{Frame et al., 2014} ). Sin embargo, se ha reportado una baja eficiencia en el uso del N por parte de los cereales cuando se utiliza urea como fertilizante ( ^{Muurinen et al., 2007} ; ^{Manlla et al., 2012} ). La eficiencia se reduce por la volatilización del NH₃ ( ^{Frame et al., 2012} ), desnitrificación (que genera N₂, NO, y N₂O) (Wilson et al., 2015), lixiviación y escorrentía de NO₂ ^- y NO₃, fijación de NH₄ ⁺, e inmovilización ( ^{Arregui et al., 2005} ; ^{García-Serrano et al., 2009} ; ^{Gao et al., 2015} ; ^{Thapa et al., 2015} ). Además, la aplicación excesiva y frecuente de N resulta en menor disponibilidad natural del N del suelo en el largo plazo, porque se acelera la descomposición de la materia orgánica, lo que puede resultar en mayor requerimiento de FN ( ^{Mulvaney et al., 2009} ).

En Ecuador la fertilización de la cebada es deficiente porque muchos agricultores no aplican fertilizante, mientras que otros aplican sin hacer un manejo adecuado del suelo y el cultivo, por lo que, se exceden en los requerimientos de fertilización. Esto resulta en baja eficiencia en el uso del fertilizante e incrementa la contaminación ambiental ( ^{Brinkman, 1999} ; ^{Moreno et al., 2003} ; ^{Prasad y Shivay, 2015} ). De hecho, en la mayoría de países en vías de desarrollo, como Ecuador, hay una tendencia creciente en el uso de fertilizantes nitrogenados para incrementar la producción de cereales, aunque poco se conoce sobre la eficiencia en el uso del fertilizante y sobre la magnitud de las pérdidas de N por procesos naturales ( ^{Prasad y Shivay, 2015} ).

Una de las alternativas para reducir la pérdida del N proveniente de la urea y mejorar la eficiencia en el uso de N, es el uso de fertilizantes nitrogenados de lenta liberación ( ^{IDEA, 2007} ; ^{Prasad y Shivay, 2015} ) como la urea polimerizada, que es la urea normal revestida con un polímero ( ^{Bierman et al., 2015} ; ^{Gao et al., 2015} ). La urea polimerizada que se comercializa en Ecuador está revestida con tiofosfato de N-n-butiltriamida (NBPT) ( ^{Garzón y Cárdenas, 2013} ). Esta cobertura retrasa la hidrólisis, lo que resulta en una liberación más lenta de N al sistema suelo-planta, mejora la eficiencia en el uso de N, y reduce las aplicaciones de fertilizante ( ^{Hauck, 1985} ). Al usar urea polimerizada, la dosis óptima de N (DON), desde un punto de vista económico, puede bajar hasta un 50% ( ^{Patil et al., 2010} ), y por eso el 20% de la urea que se importa es del tipo polimerizada ( ^{BCE, 2015} ). Además, el uso de urea polimerizada puede incrementar la producción en un 53% en comparación con la urea normal ( ^{Fontanetto et al., 2011} ), y en algunas casos se observa incluso un aumento en el contenido de proteína en el grano ( ^{Blackshaw et al., 2011} ).

A pesar de las ventajas que se reportan con el uso de urea polimerizada, en el país todavía hay escepticismo en su uso por parte de los agricultores, porque no creen que el costo adicional en el que se incurre para comprar la urea polimerizada se compense con una mayor producción y mejor eficiencia en el uso de N. Además, su precio es 20% más alto que el de la urea normal ( ^{SINAGAP, 2015} ), no hay suficientes investigaciones de FN en cebada, y las existentes se han limitado a la evaluación de resultados mediante un análisis de varianza (ANDEVA). Aunque esta metodología es útil, no permite identificar de forma clara la respuesta del cultivo a la fertilización (ej., lineal, lineal platea, cuadrática, cuadrática platea, exponencial, o la ausencia de respuesta).

El uso de modelos matemáticos (regresiones) puede ser una mejor alternativa para explicar la respuesta de los cultivos a la fertilización ( ^{Cerrato y Blackmer, 1990} ). A partir de esta metodología se puede determinar la dosis óptima de fertilización (DOF) y la producción que se alcanza con la DOF (PDOF). Esto es importante si se desea conocer la DON y la producción que se puede alcanzar con la DON (PDON) para aumentar la eficiencia de la FN. El uso de regresiones también permite incluir otros tratamientos (e.g., tipo de fertilizante, tipo de labranza, cultivar o variedad, otro nutriente). Al usar una metodología que permite determinar la dosis que maximiza la producción y reduce los problemas ambientales, también se optimiza la rentabilidad del cultivo ( ^{Pagani et al., 2009} ). Al determinar la DON en la cebada se puede mejorar la producción, reducir los costos de esta al incrementar la eficiencia en el uso del N, mejorar el valor nutricional del grano y mitigar el impacto ambiental por el uso de fertilizantes nitrogenados. Por eso, esta investigación tuvo como objetivo evaluar la respuesta de la producción de cebada a la fertilización nitrogenada (FN) y el impacto de esta en el contenido proteico del grano y la disponibilidad de nutrientes del suelo.

MATERIALES Y MÉTODOS

Localización de la investigación

Esta investigación se desarrolló en dos sitios experimentales de la Hacienda Valencia, que administra la carrera de Ingeniería Agropecuaria de la Universidad de las Fuerzas Armadas (ESPE). Esta hacienda se encuentra en la parroquia Pintag, cantón Quito, provincia de Pichincha, Ecuador. Los sitios se localizaron en 0°22’36”S y 78°23’34”O (sitio 1), y 0°23’29”S y 78°23’14”O (sitio 2). Estos sitios estaban separados entre sí a una distancia aproximada de 2 km. El sitio 1 estaba a 2820 msnm y el sitio 2 a 2940 msnm. En el sitio 1 el cultivo se sembró el 10 de marzo de 2014 y se cosechó el 31 de agosto de 2014 (desde finales de la época lluviosa hasta mediados de la época seca), y en el sitio 2 la siembra se hizo el 16 de septiembre de 2014 y la cosecha se hizo el 21 de febrero de 2015 (desde finales de la época seca hasta mediados de la época lluviosa). De acuerdo con el Instituto Nacional de Meteorología e Hidrología (INAMHI), la zona se caracteriza por tener un clima de tipo templado frio, con una temperatura media de 14 °C y una precipitación promedio anual de 1300 mm. La mayor parte de las precipitaciones ocurren entre octubre y mayo.

Manejo del cultivo

Para esta investigación se seleccionó la variedad INIAP Cañicapa 2003, debido a que cubre cerca del 60% de la demanda para la producción de cerveza en Ecuador ( ^{Grey, 2013} ). El cultivo antecesor fue papa ( Solanum tuberosum ), con producción de tipo convencional. La siembra de la cebada se realizó dos semanas después de la cosecha de la papa en cada sitio; durante ese periodo el crecimiento de malezas fue bajo y cubrió menos del 5% del suelo. Por esa razón, no se aplicaron herbicidas previo a la siembra, y la preparación presiembra del suelo fue suficiente con dos pases de rastra pesada a una profundidad media de 0,20 m. La siembra se realizó después de la labranza, al voleo, con una densidad de 110 kg semilla/ha. Luego se necesitó un pase de rastra liviana, a una profundidad media de 0,05 m, para cubrir la semilla. Por prevención, y con base en antecedentes históricos de la zona de cultivo, dos meses después de la siembra, en los dos sitios, se fumigó con el fungicida oxicarboxin con una dosis de 300 g/100 l de agua para el control de roya ( Puccinia hordei ) de la cebada, y 2, 4-D éster en dosis de 400 ml/ha para el control de malezas de hoja ancha.

Los datos de precipitación y temperatura se recopilaron en la estación meteorológica de la carrera de Ingeniería Agropecuaria de la ESPE, que se localiza a 3 km de cada sitio experimental, y que es monitoreada por el INAMHI. La variedad INIAP Cañicapa 2003 requiere precipitaciones de 400 mm durante el ciclo del cultivo ( ^{Rivadeneira et al., 2003} ). Con respecto a la temperatura, durante la germinación la cebada requiere una temperatura mínima de 6 °C, mientras que la floración y maduración requieren una temperatura de 15 - 20 °C ( ^{Castañeda-Saucedo et al., 2009} ). Con temperaturas mayores a 20 °C disminuye la germinación, y si no hay suficiente humedad en el suelo también disminuye la producción ( ^{Rivadeneira et al., 2003} ).

Muestreo inicial de suelos

Se realizó un muestreo inicial del suelo para conocer sus condiciones nutricionales antes de establecer la investigación. Se recolectaron de forma aleatoria veinte submuestras de toda el área en cada sitio experimental, a una profundidad de 0,20 m. Se mezcló todo el suelo recolectado y 1 kg de esa mezcla se usó para análisis en el laboratorio de suelos de la Agencia Ecuatoriana de Aseguramiento de la Calidad del Agro (AGROCALIDAD). Ahí la muestra se secó a 40 °C con circulación continua de aire por 48 h, luego se trituró y pasó por un tamiz de 2 mm para eliminar impurezas, raíces, y rocas. Los análisis se realizaron siguiendo la metodología descrita por ^{Carrera (2008)} , estos incluyeron la determinación de la textura del suelo mediante el método del hidrómetro de Bouyoucos ( ^{Bouyoucos, 1962} ); el pH se evaluó con el potenciómetro al utilizar una relación 1:2 suelo:agua; la materia orgánica y el N se evaluaron mediante titulación con el método volumétrico. Para determinar la disponibilidad del P, K⁺, Ca²⁺, Mg²⁺, Fe, Mn, Cu, y Zn se realizó la extracción con el método Olsen modificado; luego el P se analizó mediante colorimetría y los demás nutrientes se analizaron con el equipo de absorción atómica.

Diseño experimental y aplicación de tratamientos

En cada sitio se utilizó un diseño de bloques completos al azar (DBCA) en un arreglo de parcelas divididas con cuatro repeticiones. Los tratamientos incluyeron el tipo de urea (normal y polimerizada) como parcela principal y los seis niveles de FN (0, 30, 60, 90, 120, y 150 kg/ha) como subparcela. Cada subparcela fue de 64 m² (8×8 m), con un total de 48 parcelas por sitio. La FN se hizo al voleo, aplicando el 50% de cada dosis tres semanas después de la siembra, y el 50% restante se aplicó ocho semanas después de la siembra.

Variables evaluadas

En cada sitio experimental se evaluaron tres variables: a) la producción de grano de cebada con respecto a la FN, b) el contenido de proteína en el grano, y c) el efecto de la FN sobre las propiedades químicas del suelo al final del ciclo de cultivo. El análisis del contenido de proteína en el grano y el análisis de suelo durante la cosecha del cultivo se realizaron solo para las subparcelas que recibieron 0, 90, y 150 kg N/ha.

Para evaluar la producción se recolectaron de forma manual tres submuestras de 1,00 m² en cada subparcela (3,00 m² en total). Para ello se colocó de forma aleatoria un marco de plástico de tipo PVC de 1 m² en tres lugares diferentes dentro de cada subparcela y se cortó la cebada a nivel del suelo. Las tres submuestras se juntaron para formar una muestra homogénea. Estas se secaron a temperatura ambiente por diez días, para luego separar el grano de la espiga con una trilladora. El grano recolectado se volvió a secar a temperatura ambiente hasta alcanzar una humedad promedio y estable de aproximadamente 10%. El peso del grano que se obtuvo de cada subparcela se corrigió para obtener la producción en t/ha.

Para el análisis del contenido de proteína se tomaron submuestras de 20 g de grano de las subparcelas seleccionadas, las cuales se enviaron para análisis al laboratorio de AGROCALIDAD. Las muestras se secaron a 60 °C por 48 h con circulación permanente de aire y luego se molieron ( ^{Shirin et al., 2008} ). En el análisis se determinó el contenido de N mediante el método Kjeldahl ( ^{García y Fernández, s.f.} ), el valor que se obtuvo, se multiplicó por el factor de 6,25 para obtener el contenido de proteína ( ^{De García y Gallardo, 2011} ).

Para el muestreo del suelo al final del ciclo de cultivo, se recolectaron cuatro submuestras en forma aleatoria de cada subparcela a 0,20 m de profundidad. Estas se mezclaron para obtener una muestra homogénea, y 1,00 kg de esa mezcla se analizó en AGROCALIDAD con los mismos métodos que se utilizó para el análisis inicial, pero sin incluir el análisis de textura.

Análisis estadístico

Análisis de varianza

En virtud de que la cebada se estableció en épocas diferentes para cada sitio, los ANDEVA para las variables evaluadas se hicieron por separado para cada sitio. Para ello, los tratamientos de tipo de urea y niveles de FN se consideraron factores fijos, mientras que los bloques se consideraron aleatorios. Para realizar los ANDEVA se utilizó el procedimiento PROC MIXED del programa estadístico SAS ( ^{SAS, 2009} ). Las diferencias entre los tratamientos se determinaron con la opción DIFF en PROC MIXED, con el uso del procedimiento de diferencia mínima significativa (DMS), y fueron consideradas significativas con un p≤0,10.

Respuesta del cultivo a la fertilización nitrogenada

Para evaluar la respuesta de la producción de grano de cebada a la FN (lineal, lineal platea, cuadrática, cuadrática platea, exponencial), se utilizaron los procedimientos PROG REG o PROC NLIN de SAS. Esta evaluación se hizo primero para cada sitio, y luego combinando los resultados de producción de los dos sitios para obtener la respuesta promedio y compararla con las recomendaciones de FN que se manejan en Ecuador. Las regresiones que se obtuvieron se consideraron significativas cuando p≤0,10, pero cuando varias regresiones fueron significativas, se seleccionó aquella con el coeficiente de determinación (r²) más alto con base en la recomendación de ^{Cerrato y Blackmer (1990)} . En los resultados, para el sitio 1 y para el promedio de los dos sitios, se identificó una respuesta de tipo cuadrática (ecuación 2), mientras que al haber una interacción significativa entre el tipo de urea y la FN con respecto a la producción de grano para el sitio 2 (ver más adelante), se identificó una respuesta de tipo lineal (ecuación 1) para la urea normal y de tipo cuadrática platea (ecuaciones 2 y 3) para la urea polimerizada.

y = a + b (1)

y = a + bx + cx² (2)

y = a + bx+ cx² si x < x_o (3)

y = a + bxo + cxo² si x ≥ x_o (4)

En estas regresiones, y representa la respuesta de la producción de cebada (t/ha) a la FN, x la dosis de FN (kg N/ha), y a (intercepto), b (coeficiente lineal), c (coeficiente cuadrático) y xo (dosis de FN en el punto de unión) son las constantes de las regresiones. La dosis óptima de N (DON) y la producción que se puede alcanzar con esa DON (PDON), se determinaron al resolver para x la regresión correspondiente, utilizando una relación de 0,002706 que resulta al dividir el precio de 1 kg de N (en US$) entre el precio de una tonelada de grano ( ^{Cerrato y Blackmer, 1990} ). Para ello se utilizaron precios referenciales de US$ 30 por un saco de 50 kg de urea y US$ 482 por tonelada de grano (US$ 22 por un saco de 45 kg) ( ^{MAGAP, 2015} ).

Correlación entre producción y contenido de proteína

También se evaluó la relación entre el nivel de producción de grano y el contenido de proteína al calcular el coeficiente de correlación de Pearson con el procedimiento PROC CORR de SAS, y utilizando los datos de los dos sitios.

RESULTADOS Y DISCUSIÓN

Análisis inicial de suelos

El análisis inicial de suelos mostró que el sitio 1 tenía una textura franca con una composición de 44, 44, y 12% de arena, limo y arcilla, respectivamente; mientras que el suelo del sitio 2 tenía una textura franco arenosa con una composición de 60, 32, y 8% de arena, limo y arcilla, respectivamente ( Cuadro 1 ). Los suelos con textura franca y sin problemas de compactación, por lo general, tienen una alta capacidad para retener agua y nutrientes, pero altas precipitaciones durante el ciclo de cultivo pueden resultar en pérdidas de nutrientes por lixiviación; mientras que los suelos con textura franco arenosa retienen poca humedad, y por su naturaleza poseen fertilidad media, alta porosidad, y percolación considerable que puede resultar en lixiviación ( ^{Schjonning et al., 1994} ). El pH en los dos sitios fue ligeramente ácido, pero es considerado óptimo para el crecimiento y desarrollo de la cebada ( ^{Garófalo et al., 2010} ). La materia orgánica se encontró en niveles altos (> 5%), pero se considera que la liberación de N es lenta, debido a que la zona tiene un clima templado frío, por lo que las aplicaciones históricas de FN para producir cebada en los dos sitios han sido cercanas a los 90 kg N/ha, en una mezcla que incluye urea y el fertilizante 15-15-15. El contenido de elementos minerales reportó concentraciones óptimas o altas de K⁺, Ca²⁺, Mg²⁺, Fe, Mn, Zn y Cu, por lo que estos nutrientes no se consideraron en esta investigación y no se necesitó realizar fertilizaciones correctivas para suplirlos.

Cuadro 1 Análisis inicial de suelo previo al establecimiento de la cebada, variedad INIAP Cañicapa 2003, en dos sitios de la Hacienda Valencia de la Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014. Table 1. Initial soil’s analysis before establishment of barley, variety INIAP Cañicapa 2003, in two sites of the Valencia Farm. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014.  

^*M.O.: material orgánica/ M.O.: organic matter

^**Rango óptimo que maneja la Agencia Ecuatoriana de Aseguramiento de la Calidad del Agro (AGROCALIDAD), institución rectora de la Red de Laboratorios de Suelos del Ecuador (RELASE) / Optimum range according to the Agencia Ecuatoriana de Aseguramiento de la Calidad del Agro (AGROCALIDAD), leading institution of the Red de Laboratorios de Suelos del Ecuador (RELASE).

Clima

En los dos sitios la precipitación fue mayor a los requerimientos del cultivo (600 y 650 mm durante el ciclo de cultivo para los sitios 1 y 2, respectivamente). En el sitio 1 se observó que en el mes de marzo (luego de la siembra) y en mayo las precipitaciones fueron altas, llegando a 260 mm/mes ( Figura 1A ), mientras que en el sitio 2 las precipitaciones fueron más frecuentes y abundantes hacia el final del ciclo de cultivo, en especial en los meses de diciembre y febrero ( Figura 1C ). Debido a que en los dos sitios el nivel de precipitación superó los requerimientos del cultivo, esto pudo afectar la respuesta del cultivo a la FN. Las excesivas precipitaciones provocan lixiviación ( ^{Perdomo et al., 2011} ) o desnitrificación ( ^{Wilson et al., 2015} ), lo que puede resultar en menor disponibilidad de N y menor producción del cultivo. Los suelos con cantidades considerables de arena tienden a ser porosos y esto facilita la lixiviación de los nutrientes en cantidades mayores a los suelos donde predomina la arcilla ( ^{Mengel y Kirkby, 2001} ), algo que pudo influir en los resultados de producción del cultivo y la respuesta a la FN en el sitio 2 (ver más adelante).

Figura 1  Datos climatológicos de la Hacienda Valencia del sitio 1 (A y B) y sitio 2 (C y D) donde se cultivó la cebada, variedad INIAP Cañicapa 2003. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015. Figure 1. Weather data of the Valencia Farm, site 1 (A y B) and site 2 (C y D), where barley, variety INIAP Cañicapa 2003, was grown. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015. 

La temperatura promedio en los dos sitios fue de 14 °C ( Figuras 1B y 1D ), excepto en octubre en el sitio 2, donde la temperatura promedio fue ½ °C inferior al resto del ciclo de cultivo. Esta temperatura es ligeramente menor a lo que se recomienda para el cultivo, en especial para la etapa de maduración ( ^{Castañeda-Saucedo et al., 2009} ), aunque no se espera un efecto significativo en la producción o contenido de proteína en el grano ( ^{CN y MAGAP, s.f.} ).

Producción

Se obtuvo una producción promedio de 4,33 t/ha en el sitio 1 y de 2,08 t/ha en el sitio 2 ( Cuadro 2 ). En Ecuador la producción de cebada, variedad INIAP Cañicapa 2003, es por lo general de 3,00 t/ha bajo un sistema de tipo convencional, por lo que en el sitio 1 se obtuvo 45% más de esa producción y en el sitio 2 un 30% menos. La diferencia en producción entre los dos sitios pudo deberse a la excesiva precipitación que se tuvo durante la investigación. Esto pudo generar procesos temporales de desnitrificación en los meses más lluviosos ( ^{Wilson et al., 2015} ), en especial en el sitio 2, donde también la textura más arenosa pudo resultar en mayor lixiviación de N en la etapa de maduración del cultivo. Esto coincide con una investigación realizada con siete cultivares de cebada en Victoria, Australia, donde se reportó que con precipitaciones mayores a 170 mm/mes y bajas temperaturas la producción disminuyó hasta un 12% ( ^{Eagles et al., 1995} ).

Cuadro 2 Producción de cebada, variedad INIAP Cañicapa 2003, en dos sitios de la Hacienda Valencia. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015. Table 2. Yield of barley, variety INIAP Cañicapa 2003, in two sites of the Valencia Farm. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014 – 2015. 

^*Valores con la misma letra y dentro de un mismo sitio (urea normal y polimerizada) no fueron diferentes / Values with the same letter within a site (normal and coated urea) were not different.

^**Los valores de p≤0,10 fueron significativos / Values of p≤0,10 were significant.

El tipo de urea no aumentó la producción en ninguno de los dos sitios (p = 0,48 y 0,34 para los sitios 1 y 2, respectivamente) ( Cuadro 2 ). La limitada respuesta de la cebada al tipo de urea, podría atribuirse a que las altas precipitaciones provocaron lixiviación de N al inicio del ciclo de cultivo en el sitio 1 y al final del ciclo en el sitio 2. Similares resultados se observaron en una investigación realizada en Colorado, USA, donde la urea polimerizada resultó en menores pérdidas de N y mostró potencial para incrementar la eficiencia en el uso de este nutriente, pero no se reflejó en un aumento en la producción de cebada ( ^{Shoji et al., 2001} ). En cambio, en una investigación sobre FN en veinte sitios-ciclos diferentes de Canadá, solo se encontró mayor producción de cebada en tres sitios-ciclos con el uso de urea polimerizada ( ^{Blackshaw et al., 2011} ); en esta investigación se concluyó que si el suelo no se satura y la temperatura es relativamente baja, la liberación de N de la urea polimerizada es lenta y eso resulta en una baja disponibilidad de N en etapas tempranas del cultivo, lo que puede reflejarse en una baja producción de cebada ( ^{Blackshaw et al., 2011} ).

La FN aumentó la producción de cebada en el sitio 1 (p<0,01), pero no lo hizo en el sitio 2 (p= 0,33) ( Cuadro 2 ). Para ^{Prystupa et al. (2008)} , una adecuada FN en las primeras etapas del cultivo que garantice la disponibilidad de 60 - 90 kg N/ha, puede mejorar el crecimiento y desarrollo de la planta, lo que resulta en mayores producciones. Pero esta respuesta puede verse afectada por el exceso de precipitación o riego que ocasiona pérdidas de N ( ^{Thapa et al., 2015} ). Para ^{Stevens et al. (2015)} , quienes realizaron investigaciones de FN en cebada en doce sitios desde el año 2005 hasta el 2011 en North Dakota, USA, el nivel de producción no solo resulta de la FN, sino también del manejo del riego, la labranza, y la rotación de cultivos. Según estos autores, la rotación con leguminosas reduce la FN y mejora la producción de grano.

Respuesta a la fertilización nitrogenada

En el sitio 1 no hubo interacción entre el tipo de urea y la FN (p= 0,16), por lo que al promediar los resultados entre las dos ureas se obtuvo una respuesta de tipo cuadrática a la FN (p= 0,03; r2= 0,91), con la que se pudo determinar una DON de 90 kg N/ha y una PDON de 4,57 t/ha ( Figura 2 ). Aunque la producción es más alta que la reportada para la variedad INIAP Cañicapa 2003 ( ^{CN y MAGAP, s.f.} ), la DON fue similar a la dosis recomendada por el Instituto Nacional de Investigaciones Agropecuarias (INIAP), ente regulador del manejo de cultivos en Ecuador, y la dosis más alta recomendada por el consorcio de Cervecería Nacional ( ^{CN y MAGAP, s.f.} ), pero superior a las dosis de 45 - 60 kg N/ha que aplican muchos agricultores ( ^{Rivadeneira et al., 2003} ).

Figura 2 Respuesta de la producción de cebada, variedad INIAP Cañicapa 2003, a la fertilización nitrogenada en el sitio 1 de la Hacienda Valencia. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014 – 2015. La línea entrecortada indica la dosis óptima de N y la producción que se obtiene con esa dosis. Figure 2. Barley, variety INIAP Cañicapa 2003, yield response to N fertilization in site 1 of the Valencia Farm. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015. The discontinue line indicates the optimum N rate and the production obtained with that rate. 

En el sitio 2 hubo una interacción significativa entre el tipo de urea y la FN (p=0,03), aunque las respuestas a la FN para cada tipo de urea no fueron significativas (p>0,10), pero mostraron tendencias diferentes entre sí ( Figura 3 ). Para la urea normal se obtuvo una tendencia en la respuesta de tipo lineal creciente (p=0,15; r2= 0,44), misma que indicó que las dosis de FN utilizadas en este sitio fueron bajas con respecto a los requerimientos del cultivo, y no permitieron alcanzar la DON o calcular la PDON, resultado que refleja la alta pérdida de N, probablemente por la excesiva precipitación. Una respuesta de tipo lineal en un estudio realizado en Karaj, Irán, la obtuvieron ^{Hajighasemi et al. (2016)} . En cambio, la urea polimerizada mostró una tendencia en la respuesta de tipo cuadrática platea con una DON de 25 kg N/ha y una PDON de 2,35 t/ha. Los resultados para la urea polimerizada en este sitio fueron inconsistentes con las recomendaciones de fertilización que maneja Cervecería Nacional. En algunos casos, la variabilidad climática y de manejo del cultivo tienen alta influencia en la producción del cultivo y esto puede resultar en respuestas no significativas a la FN ( ^{Hajighasemi et al., 2016} ). La alta variabilidad en la respuesta a la FN en el sitio 2 fue un reflejo de la excesiva precipitación en este sitio durante la evaluación.

Figura 3 Respuesta de la producción de cebada, variedad INIAP Cañicapa 2003, a la fertilización nitrogenada en el sitio 2 de la Hacienda Valencia. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014 – 2015. La línea entrecortada indica la dosis óptima de N y la producción que se obtiene con esa dosis. Figure 3. Barley, variety INIAP Cañicapa 2003, yield response to N fertilization in site 2 of the Valencia Farm. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015. The discontinue line indicates the optimum N rate and the production obtained with that rate. 

Diferentes respuestas de la cebada a la FN se han reportado en Montana, USA, en una evaluación que incluyó cuatro dosis de FN (0, 40, 80, y 120 kg N/ha), en cuatro sitios y con diferentes rotaciones de cultivos; en este caso, la DON varió de 75 a 130 kg N/ha según el sitio, y la PDON fue menor a 3,00 t/ha en todos los casos ( ^{Sainju et al., 2013} ). En este estudio, las diferencias climáticas y de textura del suelo, que afectan la disponibilidad de N en el corto y mediano plazo ( ^{Cárcova et al., 2004} ; ^{Wilson et al., 2015} ), pudieron causar la diferencia en la respuesta de la producción de cebada a la FN entre los dos sitios, y entre los dos tipos de urea para el sitio 2.

A pesar de que cada sitio presenta una respuesta diferente de la producción de cebada a la FN, se promediaron los resultados para comparar la respuesta global de esta investigación con los resultados de otros trabajos. Este promedio dio como resultado una respuesta de tipo cuadrática (p<0,01; r² = 0,96), con la que se determinó una DON de 90 kg N/ha y una PDON de 3,41 t/ha ( Figura 4 ). Estas DON y PDON fueron similares a las reportadas por ^{Rivadeneira et al. (2003)} a través del INIAP en Ecuador. ^{Alley et al. (2009)} mencionaron que la FN debe ser mayor a 30 kg N/ha, pero no debe exceder los 100 kg N/ha, ya que podría resultar en una alta producción de biomasa que, bajo condiciones de alta precipitación y velocidad del viento, puede ocasionar el acame de las plantas. Además, debido al alto impacto del fertilizante en el costo de producción, en los sistemas productivos de cereales como la cebada se deberían considerar rotaciones de cultivos con leguminosas para disminuir la DON ( ^{Mulvaney et al., 2009} ). Esto ayudaría a reducir los problemas ambientales que ocasiona la FN, en especial la que utiliza fertilizantes sintéticos como la urea.

Figura 4  Respuesta promedio de la producción de cebada, variedad INIAP Cañicapa 2003, a la fertilización nitrogenada en la Hacienda Valencia. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015. La línea entrecortada indica la dosis óptima de N y la producción que se obtiene con esa dosis. Figure 4. Barley, variety INIAP Cañicapa 2003, average yield response to N fertilization in the Valencia Farm. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015. The discontinue line indicates the optimum N rate and the production obtained with that rate. 

Contenido de proteína en el grano

En Ecuador hay reportes que indican que la cebada, variedad INIAP Cañicapa 2003, puede alcanzar hasta un 17,00% de contenido proteico, pero por lo general, esta variedad tiene entre 11,00 - 14,00% de proteína, con un promedio de 13,00% ( ^{Rivadeneira et al., 2003} ). En esta investigación el tipo de urea no afectó el contenido de proteína en el grano (p=0,65 y 0,27 para los sitios 1 y 2, respectivamente). Esto difiere de lo que encontraron otros investigadores ( ^{Blackshaw et al., 2011} ), quienes obtuvieron mayor contenido proteico con el uso de urea polimerizada en ocho de los veinte sitios que se evaluaron en Canadá. Sin embargo, en la actual investigación la FN resultó en mayor contenido proteico en los dos sitios (p≤0,10), llegando a 14,50% con la dosis de 150 kg N/ha ( Figura 5 ). En el sitio 1 las dosis de 90 y 150 kg N/ha, resultaron en contenidos de proteína que no mostraron diferencia estadística, mientras que sí mostraron diferencia en el sitio 2. Las diferencias de clima, textura de suelo, y nivel de producción entre los dos sitios pudieron ocasionar estas variaciones del contenido de proteína. Sin embargo, los resultados de contenido proteico de esta investigación están dentro de los rangos presentados por ^{Blackshaw et al. (2011)} y de ^{Stevens et al. (2015)} , quienes reportaron contenidos promedio de proteína de 11,00 - 17,00% con una FN entre 70 - 150 kg N/ha.

Figura 5 Contenido de proteína en el grano de cebada, variedad INIAP Cañicapa 2003, en dos sitios de la Hacienda Valencia. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015.* Barras con la misma letra no son diferentes (P≤0,10). Figure 5. Grain’s protein content of barley, variety INIAP Cañicapa 2003, in two sites of the Valencia Farm. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014 – 2015. * Bars with the same letter aren’t different (p≤0,10). 

Altos niveles de FN producen un grano con mayor contenido proteico, lo cual es deseable en la industria de procesamiento de alimentos para animales y humanos ( ^{López et al., 2007} ). Sin embargo, contenidos altos de proteína no son deseables en la industria cervecera, porque bajan el contenido de glúcidos, y la fracción insoluble (hordeína y glutelina) puede producir turbidez en la cerveza ( ^{Giménez et al., 2007} ). De hecho, la industria cervecera prefiere altos contenidos de almidón y bajos contenidos de proteína (menores al 14,00%, ideales entre el 10,00 - 12,00%) para asegurar la producción de espuma ( ^{CN y MAGAP, s.f.} ; ^{Echagüe et al., 2001} ; ^{Jankovic et al., 2011} ). Los resultados indican que si el destino de la cebada es la producción de cerveza, no se debería exceder la FN sobre el nivel de la DON, puesto que esto no solo resulta en el uso ineficiente del fertilizante, sino también en pérdidas económicas por el elevado contenido de proteína y menor contenido de almidón en el grano.

Al correlacionar el nivel de producción con el contenido de proteína en el grano, se encontró una relación positiva de tipo lineal (p<0,01; r2 = 0,27; coeficiente de Pearson = 0,52) ( Figura 6 ). La relación fue baja, probablemente porque la cantidad de datos (48 en total) no fue lo suficientemente alta como para obtener una respuesta más significativa, y también por la variabilidad en la producción entre los sitios 1 y 2. Sin embargo, relaciones similares se han reportado en cebada y trigo ( Triticum aestivum ) en suelos deficientes en cuanto a la disponibilidad natural de N ( ^{Eagles et al., 1995} ; ^{Jankovic et al., 2011} ; ^{Hafeez, 2012} ; ^{Sainju et al., 2013} ; ^{Stevens et al., 2015} ). Pero esta respuesta contradice de lo expuesto por ^{Blackshaw et al. (2011)} , quienes mencionaron que altos niveles de producción crean un efecto de dilución que resulta en menor contenido de N en el grano, lo que se refleja en menor contenido de proteína. Al parecer, cuando la disponibilidad natural de N en el suelo es baja, la FN no solo promueve la producción de grano, sino que también incrementa su valor nutricional.

Figura 6 Correlación entre el contenido de proteína en el grano y la producción de cebada, variedad INIAP Cañicapa 2003, obtenida en la Hacienda Valencia. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015. Figure 6. Correlation between grain’s protein content and yield of barley, variety INIAP Cañicapa 2003, in the Valencia Farm. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015. 

Análisis de suelo al final del ciclo del cultivo

Al final del ciclo del cultivo el tipo de urea no resultó en cambios significativos en el pH del suelo (p=0,38 en ambos sitios) ( Cuadro 3 ). Esto es contradictorio con los resultados que reportan ^{Bierman et al. (2015)} , quienes mencionaron que al haber una nitrificación más lenta de la urea polimerizada, los cambios de pH del suelo son menores en comparación con los que causa la urea normal. La leve, pero significativa disminución del pH del suelo con el aumento de la FN (p<0,10) en los dos sitios, pudo deberse a que los fertilizantes amoniacales como la urea generan acidez en el suelo, ya que en el proceso de nitrificación se libera H+ ( ^{Bierman et al., 2015} ; ^{Gao et al., 2015} ). Además, debido a su baja capacidad buffer, los suelos con texturas de tipo arenosa y franca tienden a acidificarse más rápido que los suelos arcillosos con la aplicación de fertilizantes ( ^{Perdomo et al., 2011} ), lo que explica el menor grado de acidificación en el suelo del sitio 1 (suelo franco) en comparación con el sitio 2 (suelo franco arenoso), en el periodo de evaluación con la dosis más alta de FN (reducción de 0,10 en el sitio 1 y de 0,50 en el sitio 2). La acumulación de H+ desplaza los cationes de las cargas del suelo, lo que intensifica la potencial pérdida por lixiviación de cationes básicos como el Ca²⁺ ( ^{Campillo y Sadzawka, 2010} ; ^{Sierra, 2010} ). Sin embargo, la intensidad de la acidificación y su efecto sobre la disponibilidad de nutrientes dependen de factores como la cantidad de fertilizante que aplica el agricultor, el periodo de aplicación, el manejo del suelo, y el nivel de precipitación ( ^{Campillo y Sadzawka, 2010} ).

Cuadro 3 Análisis final de suelo durante la cosecha de la cebada, variedad INIAP Cañicapa 2003, en dos sitios de la Hacienda Valencia. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014-2015. Table 3. Final soil’s analysis during harvest of barley, variety INIAP Cañicapa 2003, in two sites of the Valencia Farm. Universidad de las Fuerzas Armadas. Pintag, Quito, Pichincha, Ecuador. 2014 – 2015. 

^*M.O.: materia orgánica/ M.O.: organic matter.

^**Valores seguidos de la misma letra en una misma columna no fueron diferentes / Values with the same letter within a column were not different.

^***Los valores de p≤0,10 fueron significativos / Values of p ≤ 0,10 were significant.

No se observó un efecto significativo del tipo de urea, o una interacción del tipo de urea con la FN en la mayoría de contenido de nutrientes del suelo, al final del ciclo del cultivo en ninguno de los dos sitios (p>0,10). Esto pudo deberse a que el periodo de evaluación fue relativamente corto (solo un ciclo de cultivo), y se necesitarían evaluaciones a largo plazo con continuas aplicaciones de FN para encontrar diferencias significativas ( ^{Pantoja et al., 2015} ). Sin embargo, se observó una leve disminución en la disponibilidad de Mg²⁺ en el sitio 2 al aumentar la FN (p=0,07), pero no ocurrió lo mismo con otros nutrientes. Al no haber un efecto significativo de la FN en la producción de cebada para este sitio, se puede inferir que este resultado refleja cierta variabilidad espacial de la concentración de Mg²⁺ en el sitio 2 ( ^{Pantoja et al., 2015} ).