Introducción
La fertilización, una de las prácticas agrícolas con mayor impacto en el rendimiento y calidad del melón1. Usualmente, se la utiliza para maximizar su producción, sin embargo, la nutrición mineral también tiene impacto importante en calidad y vida en anaquel de la fruta cosechada2-4. Muchos de estos problemas pueden ser causados por deficiencias, desequilibrios o toxicidad nutricional. Además, afectan la susceptibilidad del cultivo al ataque de plagas y enfermedades, que repercuten en su rendimiento y calidad4.
La fertilización química5-7 está asociada con otros factores, como suelo, clima, riego, prácticas agrícolas, que en conjunto determinan la calidad y cantidad de fruto. Por ello es importante las funciones de los nutrientes, interacciones, fuentes e interacciones, las fuentes y dosis de estos nutrientes, las etapas de crecimiento, desarrollo del cultivo (genética) para lograr que la fertilización cumpla con su papel de alcanzar altos rendimientos y calidad de la fruta1,8.
En el Recinto Puerto la Boca de la Parroquia Puerto Cayo en Manabí, Ecuador, los agricultores de la Asociación Agro-artesanal de Puerto la Boca Cantagallo cuentan con 54 invernaderos de 500 a 1000 m2 (2.74 ha), producen diferentes hortalizas, entre ellas el melón9, desarrollaron estrategias sobre fertilización edáfica básica con macronutrientes para sus cultivos, y generalmente los manejan de forma empírica, aunque ellos están conscientes de la necesidad de lograr cosechas rentables, con calidad en poscosecha y larga vida en anaquel. Por tanto, son conscientes de la necesidad de implementar prácticas de manejo de los cultivos, incluyendo la fertilización, que tiene gran impacto en la calidad del producto cosechado. Esto llevó a que se realice la presente investigación, cuyo objetivo fue determinar una estrategia química de fertilización para la producción y calidad de fruto de melón (Cucumis melo L.) en invernadero.
Materiales y métodos
Ubicación. El experimento fue desarrollado de mayo a agosto 2019, en un invernadero de 1000 m2 con temperatura promedio de 30° C, del Recinto Puerto la Boca perteneciente a la parroquia Puerto Cayo del cantón Jipijapa, ubicado 1°18'20'' latitud Sur, 80°45'42" longitud Oeste, a una altitud aproximada de 53 msnm, temperatura 24.8 °C, precipitación promedio anual 298 mm, concentrándose la mayor cantidad de lluvia en febrero, agosto el mes más seco9.
Tratamientos. Consistieron en 3 dosis de fertilización química (T1 Bajo, T2 Medio y T3 Alto) en 3 estadios de desarrollo (E1 al follaje, E2 a la floración y E3 a fructificación) (Tabla 1).
Diseño experimental. Bloques completamente aleatorios (DBCA) con 4 repeticiones10, al igual para el análisis bromatológico. El diseño permitió bloquear el efecto de la luz solar en el invernadero en las unidades experimentales (UE). Cada UE estuvo constituida por 3 hileras, cada hilera/tratamiento tuvo 41 plantas y la UE para cada tratamiento constaba de 123 plantas11. Para el estudio de las variables respuesta de cada UE se eligieron 10 plantas al azar de cada repetición, con un total de 160 plantas evaluadas al trasplantar al sitio definitivo, el distanciamiento fue 0.20 m entre plantas y 1.20 m entre hileras. La parcela experimental tuvo un total de 1980 plantas11.
Variables de estudio12. Altura de planta (cm) (AP). Se consideró 10 plantas al azar de cada tratamiento Diámetro de tallo (mm) (DT). Se evaluó cuando el cultivo alcanzó 50 % de floración. Número de nudos (NN). Se procedió a contarlos hasta el primer fruto, también se contaron aquellos que se encuentran después del fruto. Número de flores totales (NFT). Se contó el número de flores por planta evaluada. Número de flores fecundadas (NFF). Se contó el número de flores fecundadas por planta evaluada. Número de frutos (NFR). Se evaluó en todas las cosechas y comprendió el total de frutos cosechados. Volumen del fruto (cm 3 ) (VF). Se evaluó en base a largo, ancho y alto del fruto. Peso del fruto (kg) (PF). Se evaluó por fruto cosecha, utilizando una balanza digital.
Tabla 1 Dosis de elementos que se utilizaran en los tratamientos
Tratamiento | Estadios | Dosis (g/parcela) | Dosis (g/planta) |
---|---|---|---|
T1 | E1 | N = 986, P = 223.4: K = 0 | N = 7.95, P = 1.80, K = 0 |
T1 | E2 | N = 986, P = 223, K = 1259 | N = 7.95, P = 1.80, K = 10.15 |
T1 | E3 | N = 0, P = 0, K = 1259 | N = 0, P = 0, K = 10.15 |
T2 | E1 | N = 1391, P = 656, K = 0 | N = 11.22, P = 5.29, K = 0 |
T2 | E2 | N =1391, P = 656, K = 1379 | N = 11.22, P = 5.29, K = 11.12 |
T2 | E3 | N = 0, P = 0, K = 1379 | N = 0, P = 0, K = 11.12 |
T3 | E1 | N =1803, P = 816, K = 0 | N =14.54, P = 6.58, K = 0 |
T3 | E2 | N =1803, P = 816, K = 1512 | N =14.54, P = 6.58, K = 12.19 |
T3 | E3 | N = 0, P = 0, K = 1512 | N = 0, P = 0, K = 12.19 |
T4 | Testigo (sin fertilización) |
E1 Estadio de desarrollo, E2 Estadio de floración, E3 Estadio de fructificación. Para las formulaciones se utilizó: Urea (46 %), Superfosfato (P2O5) y Cloruro de potasio (Cl K)
Variables de estudio en poscosecha. En el Laboratorio de Bromatología de la Universidad Estatal del Sur de Manabí se realizó el análisis de Índice de refracción en pulpa (IRpulpa), Índice de Refacción en núcleo (IRnúcleo), Sólidos Totales Pulpa, medido en °Brix (STpulpa), Sólidos Totales Núcleo, medido en °Brix (STnúcleo), Temperatura Pulpa (Tpulpa), Temperatura Núcleo (Tnúcleo) de los frutos de cada tratamiento. Para el análisis se extrajo una gota de jugo de 2 frutos y se colocó en un refractómetro ABBE marca Reichert Mark II Plus11.
Análisis del benéfico/costo. Permitió definir la rentabilidad o no, si el resultado fuera mayor a 1 fue rentable13.
La preparación del suelo en el invernadero se hizo manualmente, se realizó la remoción, y desterronado de platabandas. Se aplicó materia orgánica (biocompost comprado en casa comercial) para proporcionar un suelo adecuado a las plantas al momento de su trasplante. El biocompost fue aplicado a razón de 75 kg por hilera de 33 m. Luego se medió el terreno con una cinta métrica y estacas de madera, para la formación de las platabandas de 0.50 m de ancho por 33 m de largo. El sustrato se preparó con biocompost, hoja de guaba y tierra del lugar, en proporción 2:1:1. Se puso 10 kg de humus y 1 bolsa (10 g) de micorriza adquiridos de casa comercial, para evitar el ataque de patógenos que causan damping off14. Una vez preparado el sustrato, se procedió a llenar las bandejas germinadoras, teniendo cuidado de humedecerlas, para posteriormente sembrar las semillas de los parentales. El riego de las bandejas se realizó 2 veces por día durante 20 min según la necesidad del cultivo y la capacidad de campo. Se aplicó el desinfectante Captan, en dosis de 30 g/5 L de agua, para evitar damping off14. El trasplante se realizó en las hileras, se cavaron hoyos de 0.15 m a distancia de 0.30 m entre plantas y 1.50 m entre hileras, después se procedió al trasplante.
El control preventivo de acuerdo a los antecedentes de presencia de enfermedades en la zona9. Para el control del mildiu velloso causado por el Oomycetes Pseudoperonospora cubensis y otras manchas foliares, se realizó la aplicación de Metalaxyl + Mancozeb (Ridomil) (2.5g/L) alternado con Clorotalonil (2.5 ml/L), Trichoderma (3 mL/L) y Bacillus subtilis (3 mL/L) a partir del octavo día después del trasplante9.
El control de plagas se realizó de acuerdo al monitoreo y aplicación del umbral de daño para el control de insectos plagas como mosca blanca (Bemisia ta-
baci), negrita (Prodiplosis longifila) y pulgones (Myzus persicae), se utilizó la aplicación de Thiamethoxam (0.25 mL/L), alternando con Abamectina (2.25 mL/L), Confidor (0.60 g/L) y Neem (4 mL/L), a partir de los 10 días después del trasplante15.
La poda se realizó en una sola rama principal y eliminando las ramas restantes, se eliminaron las hojas viejas y los brotes para evitar la formación de otras ramas secundarias. El tutoraje se realizó después de la poda, después de cada poda se trató con un fungicida de contacto para evitar enfermedades causadas por las heridas15. La cosecha se realizó a partir de los 120 días después del trasplante (ddt).
Análisis estadístico. Una vez que los datos satisficieran los supuestos de normalidad y homogeneidad de varianza, fueron analizados bajo un DBCA, de acuerdo al modelo aditivo lineal10.
Se realizaron análisis de varianza (ANVA) para probar hipótesis acerca de los efectos fijos, así como comparaciones de medias de los tratamientos mediante la prueba de Tukey al P<0.05. El ANVA también fue para estimar los componentes de varianza para los efectos aleatorios. Los análisis indicados se realizaron utilizando el Proc GLM del SAS University16. Se realizó análisis de correlación mediante el coeficiente de Pearson entre las variables de respuesta correspondientes10.
Resultados
Análisis de caracteres agronómicos y de rendimiento. La prueba de Kolmogórov-Smirnov P<0.01, no fueron significativos, sugiriendo una distribución normal. Los coeficientes de variación (CV), estuvieron entre 0.72 a 13 %, rango permitido para este tipo de investigaciones (Tabla 2).
Tabla 2 Análisis de normalidad para caracteres agronómicos
AP | DT | NN | NFT | NFF | NFR | VF | PF | |
SD | 3.86 | .51 | 1.27 | 2.55 | .59 | .29 | 1574.24 | .18 |
Varianza | 14.91 | .26 | 1.63 | 6.53 | .0.35 | .08 | 2478235.3 | .03 |
Asimetría | .0.31 | .29 | -.14 | -.34 | .28 | .31 | -0.04 | -.30 |
Curtosis | -1.11 | -1.27 | -.14 | -1.07 | -1.31 | -1.74 | -2.02 | -1.84 |
CV (%) | .72 | 2.63 | 3.98 | 6.42 | 6.77 | 10.90 | 6.20 | 13.107 |
K-S | >.15 ns | .083 ns | >.15 ns | <.15 ns | .096 ns | .021 ns | >.15 ns | <.01 ns |
SD desviacion standart, ns no significativo, CV (%) porcentaje del coeficiente de variación, AP altura de planta, DT diámetro de tallo, NN número de nudos, NFT número de flores totales, NFF número de flores fecundadas, NFR número de frutos, VFR volumen de fruto, PF peso de fruto, K-S Kolmogorov-Smirnov
El análisis de homogeneidad de varianza (Tabla 3), la prueba de Chi-cuadrada no fue estadísticamente significativa P<0.01, que confirmó que las varianzas son comparables.
Con los resultados de ajuste a la distribución normal de las variables y la existencia de la homogeneidad de varianza, se vio conveniente continuar con el proceso de ANVA.
Análisis de varianza para caracteres agronómicos y de rendimiento. La Tabla 4, el CV de las variables evaluadas está entre 0.72 a 13.07 % está dentro los rangos permitidos.
El ANVA para caracteres agronómicos de rendimiento (Tabla 4), fueron altamente significativos P<0.01. Esto estaría indicando que, para AP, DT, NN, NFT, NFF, NFR, VF y PF al menos uno de los tratamientos fue diferente.
Tabla 4 Análisis de varianza para caracteres agronómicos y de rendimiento. Jipijapa 2019
FV | gl | Cuadrados Medios | |||||||
AP | DT | NN | NFT | NFF | NFR | VF | PF | ||
Total | 15 | ||||||||
Repeticiones | 3 | 7.57 | .0001 | 1.03 | 1.79 | .052 | .01 | 231438.19 | .004 |
Tratamientos | 3 | 60.69** | 1.27** | 4.94* | 25.98** | 1.60** | .35** | 11878815.58** | .140** |
Error | 9 | 2.09 | .01 | .74 | 1.64 | .04 | .02 | 93640.91 | .010 |
C.V. (%) | .72 | 2.63 | 3.98 | 6.42 | 6.77 | 10.90 | 6.20 | 13.107 |
AP altura de planta, DT diámetro de tallo, NN número de nudos, NFT número de flores, NFF número de flores fecundadas, NFR número de frutos, VF volumen de fruto, PF peso del fruto, *Significativo P<0.05, ** altamente significativo P<0.01.
Tabla 5 Análisis de medias mediante la prueba múltiple de Tukey al P<0.05 para caracteres agronómicos y de rendimiento
Trat | AP (cm) | DT (mm) | NN | NFT | NFF | NFR | VF (cm3) | PF (kg) |
---|---|---|---|---|---|---|---|---|
T2 | 205.05a | 4.92a | 23.00a | 22.00a | 3.05 | 1.67a | 6616.74a | .97a |
T3 | 200.05bc | 5.37b | 22.00ab | 22.00a | 3.62 | 1.52a | 6181.95a | .89a |
T1 | 109.05bc | 3.90c | 21.00ab | 18.00b | 2.40 | 1.12b | 3747.65b | .66b |
T4 | 196.68c | 3.64d | 20.00b | 18.00b | 2.25 | 1.08b | 3169.89b | .58b |
DSH | .0.05 | .24 | 1.89 | 2.82 | .42 | .32 | 675.00 | .22 |
AP altura de planta, DT diámetro de tallo, NN número de nudos, NFT número de flores totales, NFF número de flores fecundadas, NFR número de frutos, VF volumen de fruto, PF peso de frutos, DSH diferencia significativa honesta.
El análisis de medias (Tabla 5), mediante la comparación múltiple de Tukey para la variable AP hay diferencias significativas P<0.05, T2 tuvo mejor comportamiento (2.05m), respecto T4 (1.96 m).
Para DT se observó T2 (4.92 mm), fue altamente significativo respecto T4 (3.64 mm). En el NN se observó diferencias significativas P<0.05, T2 tuvo 23 nudos/planta, respecto T4 con 20 nudos/planta.
En el NFT se observó diferencias significativas P<0.05. T3 y T2 obtuvieron una media de 22 flores/planta con respecto T4 tuvo 18 flores/planta. Para el NFF, se observó T3 fue altamente significativo, obteniendo una media de 4 flores fecundadas/planta, respecto T4 que tuvo 2 flores fecundadas/planta).
El NFR con diferencias significativas P<0.05, T2 con un promedio de 2 frutos/planta respecto T4 que tuvo 1 fruto/planta. Para el VF, se observó diferencias significativas P<0.05, T2 con un volumen de 6616.74 cm3, respecto T4, que tuvo un volumen de 3169.89 cm3. En el PF, se observó diferencias significativas P<0.05 T2 con una media de 0.97 kg/fruto, en relación al T4 con un peso de 0.58 kg de peso/fruto.
Análisis de fruto en poscosecha.
Análisis de normalidad. La prueba de significancia con Kolmogorov-Smirnov (K-S) al P<0.01 (Tabla 6), no hubo diferencias significativas para IRpulpa, IRnúcleo, STpulpa, STnúcleo, Tpulpa y Tnúcleo, todas tuvieron distribución normal.
Tabla 6 Análisis de normalidad para las variables en poscosecha
IRpulpa | IRnúcleo | STpulpa | STnúcleo | Tpulpa | Tnúcleo | |
DS | .005 | .004 | 1.25 | 1.50 | 1.20 | 4.14 |
Varianza | .000025 | .00002 | 1.57 | 2.26 | 1.46 | 17.17 |
Asimetría | -.57 | -1.27 | -.64 | -.03 | 2.51 | -2.31 |
Curtosis | -1.93 | -.43 | .22 | -.67 | 8.76 | 4.53 |
CV | .37 | .33 | 15.00 | 15.86 | 4.29 | 15.50 |
ns no significativo, IRpulpa índice de refacción en pulpa, IR núcleo índice de refacción en núcleo, STpulpa sólidos totales pulpa, STnúcleo sólidos totales núcleo, Tnúcleo temperatura núcleo, Tpulpa temperatura pulpa. CV porcentaje de coeficiente de variación.
Asimismo, se ratificó mediante la prueba de Chi-cuadrada P<0.05 que no hubo diferencias significativas para las varianzas, estas fueron homogéneas. La no existencia de diferencias significativas del análisis de normalidad y de la homogeneidad de varianzas sugiere que se pueda proseguir con los ANVA.
Análisis de varianza en poscosecha. El ANVA frutos en poscosecha no hubo diferencias significativas para ninguna de las variables evaluadas (Tabla 7).
Los análisis de medias para IRpulpa y IRnúcleo, STpulpa y STnúcleo, Tpulpa y Tnúcleo, no tuvieron diferencias significativas P<0.05 entre los tratamientos (Tabla 8), indicando esto, que todas fueron iguales estadísticamente. Sin embargo, hubo una tendencia de que T2 y T1 fueron las mejores para IRpulpa, IRnúcleo y STnúcleo. En cambio, para STpulpa, Tpulpa y Tnúcleo los mejores fueron T2 y T4.
Tabla 7 Análisis de varianza para cosecha y poscosecha
FV | gl | Cuadrados medios | |||||
IRpulpa | IRnúcleo | STpulpa | STnúcleo | Tpulpa | Tnúcleo | ||
Total | 15 | ||||||
Rep. | 3 | .00 | .00 | 1.54 | 4.42 | 1.41 | 35.16 |
Trat. | 3 | .00 ns | .00 ns | .40 ns | .35 ns | 2.31 ns | 14.96 ns |
Error | 9 | .00 | .00 | 2.09 | 2.18 | 1.19 | 11.92 |
CV | .37 | .30 | 16.82 | 15.57 | 3.87 | 12.91 |
IRpulpa índice de refracción pulpa, IRnúcleo índice de refracción núcleo, STpulpa sólidos totales pulpa, STnúcleo sólidos totales núcleo, Tpulpa temperatura pulpa, Tnúcleo temperatura núcleo, ns no presenta significación, CV porcentaje de coeficiente de variación.
Tabla 8 Análisis de medias de variables de poscosecha mediante la prueba múltiple de tukey al P<.05 de probabilidad
Trat | IRnúcleo | STpulpa | STnúcleo | Tpulpa |
---|---|---|---|---|
T1 | 1.35ns | 8.42ns | 9.85ns | 27.12ns |
T2 | 1.35ns | 8.20ns | 9.45ns | 28.90ns |
T3 | 1.34ns | 8.07ns | 9.50ns | 28.22ns |
T4 | 1.34ns | 8.80ns | 9.12sn | 28.50ns |
DSH | .008 | .009 | 3.25 | 2.00 |
ns no significativo, IRnúcleo índice de refracción del núcleo, STpulpa sólido total pulpa, STnúcleo sólido total núcleo, Tpulpa temperatura pulpa, DHS diferencia significativa honesta.
Análisis del Beneficio/Costo de los tratamientos
Beneficio/Costo (USD) supermercado. En general todos los tratamientos presentaron un B/C>1, cuando el producto fue comercializado a 0.50 dólar/kg de fruto (Tabla 9). Los T2 y T3 fueron los más rentables, con una relación de B/C=3.10 y 2.76 respectivamente. Indicando esto que, por cada dólar invertido en el T2, se logra un beneficio de 3.10 centavos de dólar, asimismo en el caso del T3, por cada dólar invertido se ganaría 2.76 centavos de dólar.
Beneficio costo (USD) mercado de abasto. Cuando el producto fue comercializado en el mercado de abasto local solo T2 y T3 fueron rentables con una relación B/C de 1.46 y 1.26 respectivamente (Tabla 10).
Un análisis de comparación entre el producto comercializado en supermercados y mercado local (Figura 1), las ganancias por la venta en supermercados, son mayores al 50 % en todos los casos
Discusión
Se determinó que la estrategia de fertilización fue T2 aplicado en diferentes estadios del hibrido Edisto, logrando frutos de 0.97 kg, que sería el equivalente a 2.48 t/1000 m2. Al respecto Bazán Quintana17, determino que al aplicar un nivel de fertilización de 150-100-150 kg/ha de NPK, obtuvo el mayor rendimiento comercial de 4.79 t/1000 m2, y observó que conforme se incrementan los niveles de fertilización, la cantidad de frutos no comerciales disminuían18. Pacheco-Avalos & Quiroz19 indicaron que, para lograr 1 t de fruto, el cultivo de melón extrae 1.5 kg N, 0.65 kg P y 3.37 kg K. Mateo Box20, reportó un experimento en campo para producción de melón el nivel fue de 12.5 - 2 - 2.30 kg de NPK, con un rendimiento de 2.50 t /1000 m2. PROAIN21 mencionan que aplicando un nivel de 24 - 6 - 39 kg de NPK se puede producir entre 2.70 a 3.50 t/1000 m2. Todos estos trabajos reportados señalan que la fertilización no ocasiona una respuesta uniforme en el rendimiento y que mucho depende de la composición química del suelo, el tipo de suelo (textura, estructura), el clima, el agua, el pH, la temperatura y otros factores17,21,22.
Tabla 9 Análisis del B/C para los tratamientos de fertilización para un invernadero de 1000 m2, cuyo producto fue comercializado en supermercado
Cultivar | NP (1000 m2) | Peso/cosecha (kg) | Peso/cosecha (kg) | Precio/kg (USD) | Beneficio Bruto(1000 m2) | Costo (1000 m2) | BN (USD) Comisariato | Relación B/C | Rentabilidad |
---|---|---|---|---|---|---|---|---|---|
T1 | 3135 | .65 | 2037.75 | .5 | 1018.88 | 371.15 | 647.73 | 1.75 | Rentable |
T2 | 3135 | .97 | 3040.95 | .5 | 1520.48 | 371.15 | 1149.33 | 3.10 | Rentable |
T3 | 3135 | .89 | 2790.15 | .5 | 1395.08 | 371.15 | 1023.93 | 2.76 | Rentable |
T4 | 3135 | .58 | 1818.30 | .5 | 909.15 | 371.15 | 538.00 | 1.46 | Rentable |
B/C > 1.0 = Rentable.
Tabla 10 Análisis del B/C para los tratamientos de fertilización para un invernadero de 1000 m2, cuyo producto fue comercializado en mercado de abasto
Cultivar | NP (1000 m2) | Peso/cosecha (kg) | Peso/cosecha (kg) | Precio/kg (USD) | Beneficio (1000 m2) | Costo (1000 m2) | BN (USD) Mercado | Relación B/C | Rentabilidad |
---|---|---|---|---|---|---|---|---|---|
T1 | 3135 | .65 | 2037.75 | .3 | 611.33 | 371.15 | 240.18 | .65 | No rentable |
T2 | 3135 | .97 | 3040.95 | .3 | 912.29 | 371.15 | 541.14 | 1.46 | Rentable |
T3 | 3135 | .89 | 2790.15 | .3 | 837.05 | 371.15 | 465.90 | 1.26 | Rentable |
T4 | 3135 | .58 | 1818.30 | .3 | 545.49 | 371.15 | 174.34 | .46 | No rentable |
B/C > 1.0 = Rentable.
En Costa Rica, Barrientos23, para 4 genotipos de melón, entre 8.0 y 11.7 °Brix, y para el genotipo JMX-1006 dicho valor fue 9.5 °Brix, mientras Monges-Pérez24 con valor de solidos solubles de 16.1 °Brix. Charlo et al.18, obtuvieron en frutos de melón entre 9.05 y 13.05 °Brix. Castoldi et al.25, con valores entre 10.85 y 13.59 °Brix, al evaluar 5 genotipos de melón reticulado en Brasil. En nuestro estudio fue notorio observar que T1 también influyó en el contenido de sólidos totales del núcleo (°Brix), cuyos frutos fueron los más dulces con 9.85 °Brix, T2 tuvo 9.45 °Brix.
Finalmente determinamos en general, todos los tratamientos con un B/C > 1, cuando el producto fue comercializado a 0.50 dólar/kg de fruto. Sin embargo, T2 fue más rentable, con B/C=3.10 y 1.46 respectivamente. Indicando esto, que, por cada dólar invertido, ganaría 2.10 y 46 centavos de dólar. Al respecto, Kumar et al.26 observaron que con una fertilización optima de N y P mejoraban el rendimiento y la rentabilidad del cultivo27. Asimismo, pudimos observar mediante un análisis de comparación entre el producto comercializado en supermercados y mercado local que las ganancias por la venta en supermercados, son mayores al 50 % en todos los casos.