Región del estudio

El estudio se realizó en los años 2017 y 2018, bajo condiciones de secano, en el municipio de Pocyaxum, principal zona productora de soya en Campeche, que se localiza a 19° 37’ de LN y 90° 16’ de LO, con una altitud media de 30 m. El suelo del experimento es clasificado como Nitisol ródico, con textura arcillosa caolinita ^{(Palma-López et al., 2017)}. La precipitación acumulada en 2017 y en 2018 durante la estación de crecimiento fue de 1,282 y 1,290 mm, respectivamente y una temperatura media de 32.6 y 33.3 °C, respectivamente ^{(CONAGUA, 2019)}.

Diseño, tratamientos, genotipos y parcela experimental

En cada año de estudio, el experimento se estableció en un diseño de bloques completos al azar en parcelas divididas, con tres repeticiones. Las parcelas grandes fueron cuatro fechas de siembra (FS): 19 de junio (FS1), 10 julio (FS2), 2 de agosto (FS3) y 20 de agosto (FS4), mientras que las parcelas chicas fueron nueve genotipos: Huasteca 100 (H100), Huasteca 200 (H200), Huasteca 300 (H300), Huasteca 400 (H400), Huasteca 600 (H600), Tamesí, Vernal, Luziania y Mariana; estos dos últimos como testigos. La unidad experimental fueron cuatro surcos de 3 m de longitud, separados a 0.8 m; la parcela útil fueron los dos surcos centrales. La densidad de población fue 250,000 plantas ha^-1.

Manejo agronómico

La semilla se trató con una mezcla fungicida de Carboxín y Thiram a una dosis de 130 mL 100 kg^-1 de semilla y se inoculó con Bradyrhrizobium japonicum a una dosis de 365 mL 100 kg^-1 de semilla. Para la fertilización se usó la fórmula 60-80-00; el fertilizante se aplicó en banda al momento de la siembra. El control de maleza preemergente se realizó con 2 L ha^-1 de Glifosato, mientras que en postemergencia se controló la de hoja angosta y ancha con 1 L ha^-1 de Fluzifop-p-buthyl y Fomesafen. El control de plagas se realizó alternando en tres aplicaciones 100 mL ha^-1 de Novalurón y 200 mL ha^-1 de Cipermetrina.

Variables evaluadas

Se registraron datos de: altura de planta, medida en cm desde la superficie del suelo y hasta la yema terminal superior de la planta; ramas por planta, contadas en 10 plantas con competencia completa, tomadas al azar dentro de cada unidad experimental; vainas por planta, se contaron en 10 plantas tomadas al azar con competencia completa dentro de cada unidad experimental; número de semillas por planta, se cuantificaron las semillas por vaina de cada una de las 10 plantas tomadas al azar en cada unidad experimental y se multiplicó por el número de vainas por planta; peso de 100 granos, de cada unidad experimental se tomaron 100 granos al azar y se pesaron en g en una balanza de precisión; rendimiento de grano, que fue el peso experimental en g en cada parcela útil transformado a t ha^-1, previo ajuste al 14 % de humedad del grano.

Análisis estadístico

Para considerar el efecto de los años de estudio se realizó el análisis de varianza combinado de un diseño experimental de bloques completos a azar con arreglo en parcelas divididas. También se hicieron una prueba de comparación de medias mediante la diferencia mínima significativa (DMS), un análisis de componentes principales y un análisis con el modelo SREG ^{(Crossa y Cornelius, 1997)} para la interacción genotipo x ambiente, donde se consideró cada fecha de siembra como un ambiente de evaluación. Todos los análisis se hicieron con SAS® ^{(SAS, 2003)}.

Análisis de varianza

El efecto de las fechas de siembra (FS) fue el de mayor importancia en la expresión del rendimiento de grano, pues aportó el 75.3 % de la variabilidad total del mismo, mientras que la altura de planta aportó 59.6 %, el número de ramas por planta 52.3 % y el peso de cien granos 47.5 %. Además, entre FS hubo diferencias estadísticas significativas (Cuadro 1) para todas las variables (P ≤ 0.01). Al respecto, ^{Yari et al. (2013)}, en un trabajo similar, encontraron que las fechas de siembra aportaron la mayor variabilidad en la expresión del rendimiento (62.9 %) y de la altura de planta (54.1 %); esto en un ambiente con menor precipitación (720 mm anuales) que Campeche.

La variación debida a los genotipos (G) tuvo un impacto menor al de FS, pero estadísticamente importante (P ≤ 0.01), pues generaron 19.7 % de la variabilidad en el número de semillas, el 16 % en ramas por planta y el 11.8 % en la expresión del rendimiento de grano. En el número de vainas y de semillas por planta se observó un impacto similar al de FS, pero en esta última variable el genotipo mostró una mayor influencia, pues superó al efecto de fechas de siembra (16.5 %) y al de años (17 %). Al igual que entre FS, entre genotipos hubo alta significancia (P ≤ 0.01) para todas las variables.

El factor años (A) generó el 49.9 % de la variabilidad total en la expresión del número de vainas por planta y el 22.5 % del peso de cien granos y únicamente el número de semillas y el de ramas por planta mostraron significancia (P ≤ 0.01) entre años. Con respecto a las interacciones entre factores, la de mayor importancia fue genotipos x fechas de siembra; en ésta se generó variabilidad estadísticamente significativa para todas las variables evaluadas, lo que coincide con lo reportado en trabajos similares por ^{Robinson et al. (2009)} y ^{Matsuo et al. (2016)}. La interacción de fechas de siembra x años (FS x A) y la de genotipos x años (G x A) no resultaron significativas para rendimiento. No obstante, en esta última variable, los efectos de FS, genotipos (G) y su interacción (G x FS) aportaron la mayor variabilidad, lo que dio lugar a diferencias estadísticas significativas (P ≤ 0.01). Por otro lado, con base en los valores de los coeficientes de variación (CV) de fechas de siembra y genotipos se asume una adecuada precisión en la comparación de dichos factores de estudio ^{(Gomez y Gomez, 1984)}.

Año de producción

Solamente se observaron diferencias estadísticas (Cuadro 2) entre años para el número de semillas (8.5 semillas) y de ramas por planta (0.3 ramas). La anterior respuesta se derivó probablemente de la homogeneidad de la precipitación acumulada durante ambos años (Figura 1).

Figura 1 Precipitación registrada durante los dos ciclos de cultivo de soya en Pocyaxum, Campeche, México. 

Fechas de siembra

En el Cuadro 3 se puede observar que la evaluación de los genotipos en los ensayos establecidos el 19 de junio y 10 de julio presentó la mejor respuesta agronómica, al superar estadísticamente (P ≤ 0.05) a la de las fechas de agosto (2 y 20) en rendimiento promedio con 650 kg ha^-1, 1.5 g en el peso de 100 granos, 3.6 cm en altura de planta, 9.5 vainas, así como con 14 semillas y 1 rama por planta. Además, en la medida en que se retrasó la fecha de siembra, los genotipos redujeron su rendimiento en 5.7 g dia^-1, entre el 19 de junio y 19 de julio y durante agosto en promedio 10 g dia^-1 (Cuadro 3). Una respuesta similar en fechas de siembra es la reportada por ^{Ascencio y Maldonado (1998)}, en Tamaulipas, en 16 fechas de siembra, donde los rendimientos de las fechas de junio y julio superaron al de las de agosto. Un comportamiento semejante lo encontraron ^{Matsuo et al. (2015}, ²⁰¹⁶⁾ con siembras de junio en el suroeste de Japón.

Genotipos

En las variedades se observó que, a una mayor altura de planta y número de semillas por planta, además de un mayor peso de cien granos, estas presentaron un mayor rendimiento, y entre los componentes más asociados a este estuvieron el número y peso de semillas (Cuadro 4). De igual manera, se observó que entre más ramas tuvieron las variedades hubo mayor número de vainas, aunque estas variables fueron las que menor relación presentaron con el rendimiento de grano. Entre los genotipos, las variedades H200, H600, H100, H400 y H300 tuvieron el mejor comportamiento agronómico, pues superaron al resto en 394 kg ha^-1 de grano, 0.5 g en el peso de cien granos y en mayor altura de planta (7 cm); mientras que en el número de vainas, semillas y ramas por planta presentaron respuesta similar ambos grupos (Cuadro 4). Entre las mejores variedades destacaron por su rendimiento H200, H600, H100, H400 y H300, ya que superaron las 2.6 t ha^-1 en promedio, mientras que las de menor rendimiento fueron Vernal y Tamesí, que superaron a las variedades testigo Luziania y Mariana.

Componentes principales

Mediante este análisis (Cuadro 5) se encontró que los dos primeros componentes principales (CP) explicaron el 83.3 % de la varianza total, donde sobresalieron el rendimiento de grano, número de ramas por planta y el peso de cien granos dentro del CP1, mientras que en el CP2 destacaron el número de semillas y vainas por planta. En el Cuadro 5 también se observa que el rendimiento se asocia principalmente con el peso de semillas y número de ramas por planta; es decir, que estas mantuvieron una respuesta similar al rendimiento a través de las fechas de siembra, lo que coincide con lo reportado por ^{Van Roekel et al. (2015)}, quienes encontraron que los componentes del rendimiento de mayor importancia fueron el número y masa de los granos.

En la Figura 2 se observa que las FS de junio y julio de ambos años (2017 y 2018) tendieron a agruparse por una respuesta similar del rendimiento de grano, número de ramas, semillas y vainas por planta. Hubo una respuesta diferencial a favor del 2017 por un mayor rendimiento, pero con menor número de vainas y semillas. También destaca que en la medida en que se retrasó la siembra del 19 de junio al 20 de agosto de cada año, los genotipos evaluados redujeron su rendimiento y el número de ramas por planta. Si bien las fechas de siembra temprana (19 de junio y 10 de julio) de 2018 tuvieron más alto rendimiento, con base en un mayor número de semillas y vainas por planta, las mismas fechas tempranas en 2017 superaron en rendimiento a la de 2018, con un menor número de semillas y vainas; es decir, se produjo un mayor peso de granos por unidad de área.

Figura 2 Dispersión de los ocho ambientes de evaluación (dos años y cuatro fechas de siembra), para seis variables de nueve materiales genéticos de soya evaluados en Pocyaxum, Campeche, México. 

Las siembras tardías del 2 de agosto, que presentaron rendimientos bajos, se asociaron con un menor número de vainas y semillas por planta, mientras que las siembras del 20 de agosto presentaron un mayor número de éstas, por lo que su menor rendimiento estuvo asociado con un menor peso de granos por área.

Interacción Genotipo Ambiente

El análisis de la interacción genotipo x ambiente para el rendimiento de las variedades de soya (Figura 3) indica que H200 tuvo el mayor promedio (3.57 t ha^-1) en los ambientes con mayor potencial, los cuales fueron las FS del 10 de julio de ambos años (2017 y 2018). Este resultado indica que esta variedad tuvo un comportamiento superior al de 2.7 t ha^-1, encontrado por ^{Maldonado y Ascencio (2010)}; mientras que H100 sobresalió con un promedio de 3.12 t ha^-1 pues en las FS del 19 de junio rindió 3.401 t ha^-1 y en las del 2 de agosto 2.84 t ha^-1, en ambos años (2017 y 2018). Las FS del 02 y 20 de agosto fueron superadas en rendimiento promedio por las FS del 10 de julio y 19 de junio con 0.56 y 0.74 t ha^-1, respectivamente. La variedad H300 sobresalió únicamente en el ambiente del 20 de agosto de 2018, donde registró un rendimiento de 2.31 t ha^-1, mientras que las variedades H400 y H600, con un promedio de 2.77 y 2.92 t ha^-1, respectivamente, no destacaron en ninguno de los ambientes de evaluación. El rendimiento de las variedades testigo comerciales osciló entre 2.3 y 2.4 t ha^-1 en promedio, por lo que no presentaron un comportamiento sobresaliente en ninguno de los ambientes de evaluación.

Figura 3 Interacción genotipo x ambiente en rendimiento de nueve variedades de soya evaluadas en dos años y cuatro fechas de siembra en Pocyaxum, Campeche, México.