Crecimiento, índice de cosecha y rendimiento de frijol (Phaseolus vulgaris L.) en unicultivo y asociado con girasol (Helianthus annuus L.)

Growth, harvest index and yield of common bean (Phaseolus vulgaris L.) in sole crop and in intercropped with sunflower (Helianthus annuus L.)

EJ Morales-Rosales¹, JAS Escalante-Estrada¹, JA López-Sandoval²

¹ Facultad de Ciencias Agrícolas Universidad Autónoma del Estado de México Campus Universitario El Cerrillo. El Cerrillo, Piedras Blancas, Toluca 50200 Edo. México, México ejmorales@colpos.mx.

² Departamento de Botánica. Colegio de Postgraduados. Montecillo. Edo. México, Méx.

Artículo recibido: 15 de diciembre de 2006
Aceptado: 31 de marzo de 2008

Resumen

El objetivo de la presente investigación fue estimar el crecimiento, la producción de biomasa, el rendimiento y sus componentes principales en tres cultivares de frijol (Phaseolus vulgaris L.) en unicultivo y asociados con girasol (Helianthus annuus L.). La investigación se realizó durante la primavera de 2002 bajo condiciones de temporal en Montecillo, Estado de México. Las variables evaluadas fueron: fenología del cultivo, número de hojas m^-2 (NHV), índice de área foliar (IAF), tasa de asimilación neta (TAN), biomasa total (BT), índice de cosecha (IC), número de racimos m^-2 (NR), número de vainas m^-2 (NV), número de semillas vaina^-1 (NSV), número de semillas m^-2 (NS), peso de cien semillas (PCS) y rendimiento de semilla m^-2 (REN). La evaluación se realizó bajo un diseño experimental de bloques completos al azar con cuatro repeticiones en arreglo factorial, donde los factores de estudio fueron: cultivares de frijol y sistemas de siembra. Entre el rendimiento y sus principales componentes se realizó un análisis de correlación lineal simple. En el factor cultivares se detectaron diferencias estadísticas (p = 0.001) en todas las variables evaluadas, mientras que el factor sistemas sólo presentó diferencias significativas para las variables NR, NV, NS y PCS. Aunque la interacción de cultivares por sistema no fue significativa (p = 0.05), la producción de biomasa y semilla en los cultivares Bayomex y Michoacán se incrementó cuando se asociaron con girasol. Las variables NR, NV, NSV, NS y PCS se correlacionaron de forma positiva y significativa (p = 0.01) con el REN. Los resultados indicaron que Bayomex y Michoacán respondieron positivamente a la asociación con girasol Victoria, lo que constituye una opción para los productores del valle de Texcoco.

Palabras clave: Phaseolus vulgaris L., Helianthus annuus L., cultivos múltiples, rendimiento de semilla, tasa de asimilación neta.

Abstract

The purpose of this study was to estimate growth, biomass production, yield and yield components of three cultivars of kidney bean (Phaseolus vulgaris L.) in sole crop and intercropped with sunflower (Helianthus annuus L.). The study was carried out during the Spring of 2002 under rainy conditions in Montecillo, State of Mexico. The evaluated variables included: crop phenology, number of leaves m^-2 (NL), leaf area index (LAI), net assimilation rate (NAR), total biomass (TB), harvest index (HI), number of racemes m^-2 (NR), number of pods m^-2 (NP), number of seeds per pod (NSP), number of seeds m^-2 (NS), weight of one hundred seeds (WHS) and seed yield m^-2 (SY). The evaluation took place with an experimental design of randomised complete blocks with four replications in a factorial arrangement, where the study factors were cultivars of kidney bean and planting systems. A simple linear correlation analysis was carried out between seed yield and its principal components. Statistical differences (p = 0.001) were detected for all evaluated variables in the cultivar factor, while the systems factor presented significant differences for the variables NR, NP, NS and WHS. Although the interaction among cultivars per system was not significant (p = 0.05), biomass and seed production increased in the Bayomex and Michoacan cultivars when they were intercropped with sunflower. The variables NR, NP, NSP, NS and WHS were positively and significantly correlated (p = 0.01) with seed yield. Results indicate that Bayomex and Michoacan responded positively to being associated with the Victoria sunflower, and this represents an option for the Texcoco valley farmers.

INTRODUCCIÓN

La productividad en los cultivos asociados depende de las relaciones de complementación y de competencia entre las especies componentes del agrosistema. Algunos factores de competencia pueden ser controlados agronómicamente a través del uso de fertilizantes, densidades de población y elección adecuada de las especies en la asociación de cultivos (Sarandón & Chamorro 2003). La reducción en el rendimiento de frijol cuando se siembra asociado con otras especies, está en función de la habilidad competitiva del cultivar utilizado (Davis & García 1983). Generalmente, el frijol de crecimiento determinado (tipo I) es menos "competitivo" que los de crecimiento indeterminado (tipo II, III y IV) cuando son asociados con diferentes especies. Por esta razón, la reducción en el rendimiento en estos tipos de frijol es menos severa (Escalante & Kohashi 1993). Típicamente, los cereales como el maíz, sorgo, mijo y otros cultivos como el girasol, son especies que comúnmente se asocian con leguminosas como frijol, cacahuate y soya. Tsubo & Walker (2004) indicaron que las diferencias entre estos cultivos y las leguminosas en características morfológicas, como la estructura del dosel y sistema radical (mayor área foliar y profundidad de la raíz en los cereales y el girasol), provocaron que los cultivos componentes de la asociación tengan una "diferencia espacial" (cuando las raíces o doseles de los cultivos captan recursos o insumos en diferentes zonas) y hagan un "uso temporal" de los recursos ambientales (cuando las mayores demandas de los recursos o insumos de los cultivos se producen en periodos diferentes). Por lo que se piensa que el aprovechamiento de la radiación solar, agua y nutrimentos es mejor en las especies en asociación (Tsubo et al. 2003).

Cuando el frijol se siembra asociado, la producción de biomasa y rendimiento de semilla ha sido afectada por la competencia ejercida por la otra especie. Stern (1993), Kohashi (1996) y Santalla et al. (1999) indicaron que el rendimiento de semilla y el índice de cosecha de frijol se redujeron significativamente cuando fue sembrado asociado con maíz, debido principalmente a la competencia por insumos de cada especie. En ese sentido, Tsubo & Walker (2004) mencionaron que el maíz ejerce sobre el frijol de tipo arbustivo una fuerte competencia que puede reducir en 67% la producción de biomasa total, asimismo el frijol acumuló 40% menos materia seca en la vaina.

Ann & Francis (1985) al evaluar diferentes cultivares de frijol de hábito de crecimiento determinado e indeterminado asociados con maíz observaron reducciones de 25 y 8 % en el número de nudos, 46 y 39 % en el número de racimos, y 53 y 44 % en el número de vainas, respectivamente, con relación a sus unicultivos. En contraste, Tsubo et al. (2003) no encontraron reducciones significativas en la producción de materia seca y rendimiento de semilla en cultivares de frijol de crecimiento indeterminado cuando fueron asociados con maíz, lo anterior fue debido al menor sombreado del maíz sobre estos cultivares de frijol.

Dentro de los estudios más relevantes donde se combinan leguminosas con girasol están los siguientes: Stern (1993) afirmó que el crecimiento de las leguminosas se redujo cuando fueron sombreadas por el girasol. Sin embargo, algunas especies continúan fijando nitrógeno aún bajo estas condiciones. Asimismo, Stern (1993) registró reducciones en el rendimiento de soya y frijol cuando son sembrados en asociación con girasol. Kandel et al. (1997) encontraron una alta producción de materia seca de frijol (4098 kg ha^-1) cuando este se sembró asociado y en la misma fecha que el girasol. Ibrair et al. (2002) al evaluar diferentes leguminosas asociadas con girasol, encontraron que la competencia ejercida por esta oleaginosa, redujo el número de vainas por planta, la longitud de la vaina y el número de semillas por vaina. Morales et al. (2006) afirmaron que el frijol y el girasol se "complementan en el uso de los recursos" cuando son sembrados en asociación, en ese sentido, mencionaron que la raíz del girasol tiene la capacidad de profundizar en el suelo, lo que facilita la penetración y arraigo del sistema radical del frijol. En suelos alcalinos, el frijol sufre severos daños debido a la baja disponibilidad de nutrimentos por el pH elevado, al sembrarlo junto al girasol, la raíz de esta especie excreta iones H+ que bajan el pH y favorecen el aprovechamiento de los nutrimentos por parte del frijol (Escalante 1999). No obstante, estos antecedentes, los estudios sobre la asociación frijol - girasol son limitados por lo que se requiere más investigación sobre un número mayor de cultivares de frijol para determinar el potencial que puede ofrecer este sistema de producción. El objetivo del presente trabajo fue estudiar el crecimiento, producción de biomasa, rendimiento de semilla y sus principales componentes de tres cultivares de frijol (Phaseolus vulgaris L.) en unicultivo y asociados con girasol (Helianthus annuus L.).

MATERIALES Y MÉTODOS

La investigación se realizó en condiciones de temporal en un lote experimental ubicado en Montecillo, Estado de México (19° 29' N - 98° 54' O), a una altitud de 2 250 metros sobre el nivel del mar, con temperatura y precipitación media anual de 14.6 °C y 558 mm, respectivamente (García 1988). El suelo es de textura arcillosa con pH de 7.8 y contenido de materia orgánica de 3.8%. La siembra se realizó el 25 de mayo de 2002.

Durante el desarrollo del experimento se registró, según los criterios de Escalante & Kohashi (1993), la fenología del cultivo: días a emergencia (E), inicio de floración (R6) y madurez fisiológica (R9). Para evaluar el crecimiento de los tres cultivares de frijol en ambos sistemas de producción, se realizaron muestreos destructivos de tres plantas de la parcela útil a los 20, 34, 49, 67 y 85 días y a la madurez fisiológica (97, 103 y 120 días para Canario, Bayomex y Michoacán, respectivamente) en cada uno de los muestreos se contaron el número de hojas verdes, se midió el área foliar (sin incluir peciolo) con un integrador de área foliar (Licor 3100) y el índice de área foliar (IAF). El IAF se estimó según Escalante & Kohashi (1993) mediante la relación siguiente:

donde:

IAF = índice de area foliar, AF = área foliar por planta (cm²) y DP = número de plantas m^-2.

A la madurez fisiológica se estimó la tasa de asimilación neta (TAN), biomasa total (BT), índice de cosecha (IC), así como el rendimiento de semilla de frijol y sus principales componentes tales como: número de racimos m^-2 (NR), número de vainas m^-2 (NV), número de semillas por vaina (NSV), número de semillas m^-2 (NS) y peso de cien semillas (PCS). La TAN se estimó según Escalante & Kohashi (1993) mediante la siguiente relación:

donde:

TAN = tasa de asimilación neta; PS₂ y PSi = Peso seco final y peso seco inicial; AF₂ y AF₁ = área foliar final y área foliar inicial; ln = logaritmo natural de AF₂ y AF₁; t₂ y t₁ = tiempo final y tiempo inicial.

La BT es la resultante de la suma del peso seco de cada órgano de la planta sin considerar la raíz y se estimó después de someter tres plantas de frijol (sin considerar las hojas desprendidas) de cada tratamiento a 80 ^°C en una estufa de aire forzado hasta peso constante. El IC se estimó mediante la relación de rendimiento agronómico/biomasa total.

RESULTADOS

Variables fisiológicas

Las variables inicio de floración y madurez fisiológica presentaron diferencias estadísticas (ANDEVA; p < 0.001) (Tabla 1). La emergencia en todos los tratamientos se presentó a los nueve días después de la siembra (dds), con el 95 % en todos los genotipos. El cultivar Canario 107 fue el más precoz, seguido por Bayomex y Michoacán con días a floración y madurez fisiológica de 43 y 97; 50 y 103; y 66 y 120, respectivamente (Tabla 2).

El NHV en ambos sistemas de siembra se incrementó hasta alcanzar su máximo a los 67 dds para frijol Canario 107 y a los 85 dds para Bayomex y Michoacán (Figura 1). Los cultivares de frijol Michoacán en unicultivo y asociado con girasol fueron significativamente mayores con valores de 412 y 444 hojas verdes m^-2. Los frijoles Bayomex y Michoacán presentaron mayor NHV en asociación y Canario 107 formó más NHV en unicultivo. El IAF tuvo una tendencia similar, los mayores valores correspondieron al frijol Michoacán asociado (4.3) y en unicultivo (4.2) seguido por Bayomex (3.7 y 3.4) y Canario 107 (0.7 y 0.6) (Figura 2).

Los valores de F correspondientes a la TAN mostraron diferencias estadísticas (p < 0.001) para cultivares y en relación a sistemas e interacción cultivares por sistema no existieron diferencias estadísticas (Tabla 1). La TAN para los tratamientos en unicultivo y asociados con girasol fue de 3.5 y 3.4 g m^-2 día^-1 (Canario 107); 3.3 y 3.4 g m^-2 día^-1 (Bayomex); 4.3 y 4.31 g m^-2 día^-1 (Michoacán), lo cual indicó que entre cultivares de frijol la velocidad promedio de fotosíntesis neta durante todo el ciclo de cultivo fue mayor en Michoacán. La asociación de cultivos no afectó la TAN (Tabla 2).

Rendimiento de grano y sus componentes

El rendimiento de grano y otras características evaluadas reveló para el factor cultivares de frijol diferencias estadísticas altamente significativas (ANDEVA; p < 0.001) en todas las variables (Tabla 3). En el factor sistemas de siembra, las variables NV, NS y PCS mostraron diferencias estadísticas altamente significativas (ANDEVA; p < 0.001) y el NR fue significativo (ANDEVA; p < 0.05). La interacción cultivares de frijol por sistema de siembra no fue significativa para las variables REN, IC, NR, NV, NSV, NS y PCS. Los valores mínimos y máximos de los coeficientes de variación fueron de 9 a 21 % que correspondieron a PCS y para IAF, respectivamente.

La producción de materia seca y rendimiento de semilla fue mayor en el frijol Michoacán que superó en 50.8 y 69.2 % a Bayomex y en 78.3 y 87.8 % a Canario 107. En consecuencia, el IC fue mayor en el cultivar Michoacán, cuyo valor fue de 22.4% (Tablas 2 y 4). En los componentes de rendimiento, Michoacán fue superior al resto de los cultivares de frijol, excepto en la variable PCS donde Canario 107 mostró el mayor valor y superó a Bayomex y Michoacán en 27.2 y 27.5%, respectivamente (Tabla 4). En el sistema asociado, el NV, NS y el PCS presentaron valores superiores respecto al registrado en el sistema unicultivo, sin embargo, en este sistema de siembra se obtuvo el mayor NR.

El análisis de correlación lineal simple, mostró que los caracteres NR, NV, NSV y NS presentaron correlación positiva altamente significativa con el rendimiento de grano. El PCS presentó correlación significativa (p < 0.05) y tendencia positiva con el rendimiento de semilla (Tabla 5).

El rendimiento de semilla del girasol Victoria en unicultivo y asociado con los frijoles Canario 107, Bayomex y Michoacán reveló que no existieron diferencias estadísticas significativas en esta variable. Sin embargo, el girasol Victoria en unicultivo superó numéricamente en 11.2, 0.26 y 7.87 % a los sistemas asociados Victoria + Canario, Victoria + Bayomex y Victoria + Michoacán, respectivamente (Figura 5).

DISCUSIÓN

Variables fisiológicas

Durante el desarrollo del experimento la precipitación pluvial fue 272 mm y la temperatura media 16 °C. Con estas condiciones ambientales, los resultados de esta investigación indicaron que la producción de biomasa y rendimiento de semilla en el cultivar Michoacán de hábito de crecimiento indeterminado en ambos sistemas de siembra fue superior respecto a Bayomex (57.1, 50.8% en biomasa; 79.2, 40.2% en rendimiento) y Canario 107 (79.3, 80.3% en biomasa; 85.5, 90.2% en rendimiento) para unicultivo y asociado, respectivamente (Figuras 3 y 4). Estos resultados coinciden con Tsubo et al. (2003), quienes en condiciones similares de temperatura y precipitación no estimaron reducciones significativas en la producción de materia seca y rendimiento en cultivares de frijol de crecimiento indeterminado (tipo III y IV) cuando fueron asociados con maíz. Esta investigación confirmó lo mencionado por Santalla et al. (1999) y Tsubo et al. (2003), quienes indicaron que los cultivares de frijol de tipo III y IV, cuando son asociados con maíz y girasol, resultaron más competitivos que los de crecimiento determinado de tipo I.

Aunque no existieron diferencias estadísticas significativas para ninguna variable evaluada en la interacción cultivar por sistema de siembra, los genotipos Bayomex y Michoacán respondieron de forma positiva cuando se establecieron asociados con girasol. Estos resultados confirman lo publicado por Kandel et al. (1997), quienes afirmaron que la producción de biomasa y rendimiento de semilla es semejante en frijol cuando se establece en condiciones de unicultivo y asociado con girasol. Asimismo, el cultivar Bayomex se adaptó de forma positiva a la asociación y produjo mayor materia seca y rendimiento de grano con respecto a su unicultivo, lo que confirmó lo publicado por Morales & Escalante (2004). Ellos indicaron que el frijol Bayomex y el girasol Victoria se complementaron mutuamente en el uso de los insumos para la producción agrícola. Esto se puede deber al mayor NHV que presentó este frijol con relación a su unicultivo (Figura 1) y a que el girasol le permitió mejor explotación del medio, al profundizar más el sistema radical de la oleaginosa, facilitó a Bayomex mejor penetración y arraigo del sistema radical al suelo (Morales et al. 2006). El frijol sufre severos daños debido a la baja disponibilidad de nutrimentos en suelos alcalinos por el pH elevado, sin embargo, al sembrarlo junto al girasol, la raíz de esta especie excreta iones H+, los cuales bajan el pH y proporcionan un microambiente alrededor del área radical del frijol. Ambas condiciones permiten una adecuada asimilación de nutrimentos (Escalante 1999).

En ambos sistemas de siembra, el frijol Canario 107 presentó los promedios más bajos en la producción de biomasa y rendimiento de semilla. Lo anterior se debe a que este genotipo es susceptible al pH alcalino del terreno (el pH óptimo de este cultivar es entre 6 y 6.5) y no se adaptó a la asociación (Morales et al. 2006). Cuando este cultivar se asoció con girasol, se observó un efecto de competencia, expresado en el menor NHV e IAF, que limitó su crecimiento y desarrollo, y el cual fue provocado por el pH alcalino del suelo (7.8) y por el sombreado del girasol, que a su vez redujo la radiación interceptada y como consecuencia de esto la baja producción de biomasa y distribución de fotoasimilados hacia el rendimiento agronómico (Hawkins 1982; Tsubo & Walker 2004).

Durante el crecimiento del cultivo, el NHV y el IAF fue mayor en Michoacán lo que permite afirmar que este cultivar tuvo un mayor porcentaje de radiación interceptada y mejor eficiencia en el uso de este recurso, por lo que la producción de materia seca y rendimiento de grano se incrementó en más de 50% (Morales et al. 2006). La mayor TAN indica que el follaje del frijol Michoacán fue más eficaz en la producción de fotoasimilados y mostró mayor velocidad de fotosíntesis neta (Escalante & Kohashi 1993). La asociación no modificó la TAN de cada uno de los materiales de frijol empleados en la investigación (Tabla 2). En contraste, en un experimento donde se evaluó la asociación maíz - haba, existió un 18% más en la TAN en el haba asociada con respecto a su unicultivo (Manu 1994).

La mayor producción de biomasa y rendimiento de semilla alcanzado por el frijol Michoacán en ambos sistemas de siembra se debe a su mayor ciclo de cultivo, y a la expresión adecuada de sus características fisiológicas, tales como mayor NHV, IAF y TAN (Morales & Escalante 2004), lo que le permitió competir con el girasol.

Rendimiento y sus componentes

El sistema de siembra asociado no afectó el rendimiento de semilla y sus principales componentes (excepto en Canario 107), sino que en algunas características como NV, NS y PCS fue mayor, lo cual indica que no existió una relación de competencia entre ambas especies y sí una relación de complementación en el uso de los recursos. Lo anterior fue expresado por Fukai & Trenbath (1993), quienes señalaron que al asociar dos especies, los cultivos interactúan de tal modo que cuando uno ejerce un efecto negativo sobre el otro se establece el principio de competencia. La relación de complementación se origina cuando al menos una de las especies componentes de la asociación ejerce un efecto positivo sobre la otra.

En esta investigación, el crecimiento, producción de biomasa, rendimiento y sus principales componentes de los cultivares de frijol Bayomex y Michoacán no fueron afectados por la siembra asociada con girasol. Además en las asociaciones se obtuvo el rendimiento adicional de girasol (380 g m^-2 y 351 g m^-2, para Victoria + Bayomex y Victoria + Mi-choacán, respectivamente) (Figura 5) lo que hace atractiva esta asociación de especies y por lo tanto representa una buena alternativa para los pequeños productores de temporal del valle de Texcoco, Estado de México.

LITERATURA CITADA

Ann CE, Francis CA (1985) Bean - maize intercrops: A comparison of bush and climbing bean growth habits. Field Crops Res. 1: 319-335.         [ Links ]

Cárcova JG, Abeledo G, López PM (2003) Análisis de la generación del rendimiento: Crecimiento, partición y componentes. En: Satorre E, Benech AR, Slafer GA, de la Fuente E, Miralles D, Otegui ME, y Savin R (eds) Producción de cultivos de granos. Bases funcionales para su manejo. Editorial Facultad de Agronomía, Buenos Aires. 783 pp.         [ Links ]

Davis JHC, García S (1983) Competitive ability and growth habit of indeterminate beans and maize for intercropping. Field Crops Res. 6: 59-75.         [ Links ]

Escalante EJA (1999) Área foliar, senescencia y rendimiento del girasol de humedad residual en función del nitrógeno. Terra 17: 149-157.         [ Links ]

Escalante EJA, Kohashi SJ (1993) El rendimiento y crecimiento del frijol. Manual para la toma de datos. Colegio de Postgraduados, Montecillo. 84 pp.         [ Links ]

Fukai S, Trenbath BR (1993) Processes determining intercrop productivity and yields of component crops. Field Crops Res. 34: 247-271.         [ Links ]

García E (1988) Modificaciones al sistema de clasificación climática de Kóppen. Universidad Nacional Autónoma de México, Distrito Federal.         [ Links ]

Gomez KA, Gomez AA (1978) Statistical procedures for agricultural research. John Wiley & Sons, Manila. 653 pp.         [ Links ]

Hawkins AF (1982) Light interception, photosynthesis and crop productivity. Outlook Agric. 11: 104-113.         [ Links ]

Ibrair R, Ahmad S, Azim MM (2002) Sunflower summer legumes intercropping systems under rainfed conditions: yield and yield components. Pakistan J. Agric. Res. 17(3): 6-11.         [ Links ]

Kandel HJ, Schneiter AA, Johnson BL (1997) Intercropping legumes into sunflower at different growth stages. Crop Sci. 37: 1532-1537.         [ Links ]

Kohashi SJ (1996) Aspectos de la morfología y fisiología del frijol Phaseolus vulgaris L. y su relación con el rendimiento. Instituto de Recursos Naturales, Programa de Botánica. Colegio de Postgraduados. Montecillo. 42 pp.         [ Links ]

Manu AJ A (1994) Effects of resource on growth and development of maize and faba bean intercrop. University Nottingham Press. Nottingham. 88 pp.         [ Links ]

Morales REJ, Escalante EJA (2004) Biomass and seed yield of common bean in sole crop and intercrop. Annual report of the bean improvement cooperative 47: 207-208.         [ Links ]

Morales REJ, Escalante EJA, Tijerina CL, Volke HVH, Sosa ME (2006) Biomasa, rendimiento, eficiencia en el uso del agua y de la radiación solar del agrosistema girasol-frijol. Terra 24: 55-64.         [ Links ]

Santalla M, Casquero PA, de Ron AM (1999) Yield and yield components from intercropping bush bean cultivars with maize. J. Agron. Crop Sci. 183: 263-269.         [ Links ]

Sarandón JS, Chamorro AM (2003) Policultivos en los sistemas de producción de granos. En: Satorre E, Benech AR, Slafer GA, de la Fuente E, Miralles D, Otegui ME, y Savin R (eds) Producción de cultivos de granos. Bases funcionales para su manejo. Editorial Facultad de Agronomía, Buenos Aires: 783 pp.         [ Links ]

Stern WR (1993) Nitrogen fixation and transfer in intercrop systems. Field Crops Res. 34: 335-356.         [ Links ]

Tsubo M, Walker S (2004) Shade effects on Phaseolus vulgaris L. intercropped with Zea mays L. under well-watered conditions. J. Agron. Crop Sci. 190 (3): 168-176.         [ Links ]

Tsubo M, Mukhala E, Ogindo HO, Walker S (2003) Productivity of maize-bean intercropping in semi-arid region of South Africa. Water 29 (4): 381-388.         [ Links ]