Introduction
Alfalfa (Medicago sativa L.) is a summer legume that presents high forage yield potential, good persistence and great tolerance to drought and frost, associated with the accumulated reserves in its crown1)(2. The extensive genetic and phenotypic variation found in alfalfa allows its cultivation in diverse climates3. These characteristics make alfalfa a great alternative among pastures grown in Uruguay4)(5)(6.
The most important climatic variables influencing alfalfa growth are temperature7)(8, solar radiation9, evapotranspiration10)(11 and rainfall regime. The latter two variables are strongly associated with alfalfa productivity concerning deficit and excess of soil water. Uruguay has a sub-humid temperate subtropical climate12, presenting a relatively stable thermal regime between years. These climatic conditions allow the favorable development of alfalfa4, with optimal growth temperatures between 25-30 °C7, considering a wide range of temperatures for growth from 5 to 30 °C, with an optimal daytime temperature of 15 to 25 °C, and 10 to 20 °C optimal night temperature13)(14)(15. On the other hand, frosts in winter and early spring, as well as the intensity and frequency of water deficit in spring and summer, also limit alfalfa growth6. Annual alfalfa production varies widely in Uruguay; the contribution of each growing period to the annual yield is influenced by current climatic conditions and by pasture age, particularly in the first year after installation16. In this regard, climatic variables influence alfalfa growth differently according to the season of the year17.
This study aims to estimate the relative contribution of rainfall, solar radiation, temperature, evapotranspiration and water balance in seasonal alfalfa production in southwestern Uruguayan conditions. These findings will allow to estimate forage production based on climate forecasts and to take appropriate management measures when facing climate change scenarios.
Material and Methods
Alfalfa forage production
Analyzed data correspond to experimental results from rainfed forage alfalfa production in the framework of the National Cultivar Evaluation, INASE-INIA agreement, between 1997 and 2011 at INIA La Estanzuela (34°20’48,46"S; 57°43’48.74"W), Uruguay18)(19)(20)(21)(22)(23)(24)(25)(26)(27)(28)(29)(30)(31)(32. The seeding of the experiments was carried out in April, except in 1997, 1998 and 2007, which took place in May, and in June in 2002. The experimental unit in each trial was a plot of six furrows of 5 m long spaced 0.16 m from each other (total area 4.8 m2). Seeding density was 20 kg seeds ha-1, corrected by the germination percentage. Fertilization was applied at seeding time according to soil analysis; phosphorus was added to a final concentration of 17 ppm of P2O5, and nitrogen to a maximum of 20 kg N ha-1 to ensure implantation. Every time the average of the alfalfa genotypes reached a height of 20 cm and/or regrowths from the crown were observed, the evaluation of forage production was carried out in the three central furrows corresponding to an area of 2.65 m2. Forage harvest was conducted with a mechanical cutter with a collection bag, and field fresh weight was obtained. A forage sample was oven-dried at 60 °C to obtain dry matter percentage, to finally express the yield in kg ha-1 of dry matter. The seeding year was considered the first year of pasture life; and the following, as the second and third year of pasture life. In the latter, the test evaluation finished at the end of November by the protocol; therefore, the production of the third summer of pasture life was not recorded. Alfalfa trials were refertilized with phosphorus at the beginning of their second and third autumn. A randomized incomplete block design was conducted by triplicates. Information was extracted from the database for two cultivars used as test controls, Estanzuela Chaná (E_CHANA) and Crioula (CRIOU), both of intermediate latency. The dry matter production per hectare (kg DM ha-1) was calculated between each season cutting dates, for each cultivar, season and year. The cutting date information was used to estimate the length of the period between cuttings. Daily GR between cuttings (kg DM ha-1 day-1) was calculated from the forage yield (kg DM ha-1) and the duration of each period between cuttings. Cuttings were grouped and averaged according to seasons: autumn (A) from March 1st to May 31st, winter (W) from June 1st to August 31st, spring (Sp) from September 1st to November 30th, summer (Su) from December 1st to February 28th for each particular year, so that the GR of a particular season and year could be constituted by one, two or three cuttings.
The climatic data used, covering the period 1997-2011, were obtained from INIA La Estanzuela meteorological station, 4 km from where the experiments were carried out. The climatic variables analyzed were: minimum (Tmin), medium (Tmed) and maximum (Tmax) daily air temperature (ºC); daily incident radiation (Rs_RADIAT) (MJ m-2 day-1); accumulated incident radiation (RADIAT_Accu), period addition between cuttings, (MJ m-2 period-1); accumulated rainfall in the period between cuttings (RAIN) (mm); effective accumulated rainfall (Pe)(mm)33 and reference evapotranspiration (ETo, mm)34. Maximum pasture evapotranspiration (ETm) corresponded to the reference evapotranspiration (ETo) times the cultivation coefficient (Kc) between cuttings34. Average values of 0.40; 0.95; and 0.90 for the initial Kc, medium Kc, and final Kc were used; when pasture cuttings within the season occurred, values of 0.40; 1.20 and 1.15 were used respectively. These Kc coefficient values were applied immediately after cutting; in full coverage; and immediately before cutting, respectively. The growing season was defined as a set of individual cutting periods34. The actual evapotranspiration (ETa) was estimated through the soil water balance WinIsareg model35. This model, due to its simplicity and versatility, allows estimating, on daily basis, water content of the soil for a particular year or group of years, as well as water requirements and deficits magnitude according to the pre-established cultivation criteria or different irrigation strategies based on soil variables, potential evapotranspiration and rainfall or irrigation, and crop parameters, including simple Kc35. This procedure allowed to obtain the ETa:ETm ratio and the difference between rainfall and maximum evapotranspiration (climatic water balance, cWB, mm). The average or accumulated value of the direct or derived meteorological variables for each period between pasture cuttings was calculated, obtaining for each period the pasture GR associated with an average or accumulated value of the climatic variable or a calculated variable.
Parameters of the WinIsareg model
The daily soil water balance was estimated during the evaluation period 1997-2011. To this goal, the same INIA-LE climate database was used, together with the water retention characteristics of the soil profile where the pasture was planted, as well as the Kc, following the respective phenology and cuttings. The crop Kc was taken from the FAO guide No. 5634, following the cutting dates and growth stages adjusted to each season dates in each particular year. The soil was classified as a Brunosol Eutrico Tipico LAc V (Table 1). The intended effective depth for the root system was 60 cm with a threshold (p) = 0.50, which represents the maximum depletion fraction of available water in the soil at root depth, to which the pasture growth is not affected by the reduction of water content in the soil. The contribution of water coming from lower layers and the effect of soil salinity were considered null due to their low significance in local conditions.
Statistical analysis
To study the magnitude of variability between years and between seasons within the years of climatic variables, regardless of a particular year or season, an analysis of the variance components was carried out, adjusting a Mixed Linear Model MLM36 with a random effect of the years and seasons. The variables used in this analysis were: Tmin, Tmed, Tmax, Rs_RADIAT, RADIAT_Accu, RAIN, Pe, ETa, ETm, ETa:ETm, and cWB. To evaluate the effect of the season of the year, the pasture age, the cultivar and its interactions on the GR, the ANAVA was carried out. The mean comparison was made using Fisher’s least significant difference test (LSD) for a significance level of α = 0.05. To relate the effect of climatic variables on GR, linear regression analysis, two-section non-linear regression and multivariate analysis using regression trees (CART)37 were performed. All analyses were performed using the INFOSTAT38 software and its interface with the R39 software.
Results and Discussion
Climate variables
High dispersion of values found for the main meteorological variables in southwestern Uruguay is mainly explained by the dispersion within the seasons and the differences between seasons (Figure 1). These climatic characteristics are frequent in temperate sub-humid climates of similar latitudes12. This strong climatic seasonality conditions the growth of most pastures, both natural and cultivated (Figure 1). During the annual growth cycle, temperature, radiation and water availability will determine the seasonality of alfalfa forage production7)(13)(14. At the same time we found an important variation in the alfalfa growth in the same growing season among years (Table 2), not only explained by the climatic conditions prevailing in that season but also by the possible influence of the physiological state of alfalfa from the previous growing season40)(41.
The interdependence with the individual meteorological variables and with their associated environmental variables is relatively large42, in such a way that the influence of one of them on the growth factors is almost always associated in greater or lesser degree with another climatic variable8)(42)(43, except for extreme weather events.
Alfalfa growth rate
The alfalfa GR (kg DM ha-1 day-1) changed significantly between seasons of the year (p < 0.0001) and with pasture age (p = 0.0307). No significant differences were found in the GR between cultivars E_CHANA and CRIOU (p = 0.9376), both of intermediate latency, through crop age, seasons and years. The possible combinations of interactions between season, age and cultivar were not significant in any of the cases, allowing the evaluation of the main effects of season and pasture age separately.
Recorded average annual GR between 1997-2011 (37.94 kg DM ha-1 day-1) was similar to that obtained in other areas of similar climatic characteristics44)(45. High seasonality of climatic variables (Figure 1) determines, as in other regions, the seasonality of the GR for an alfalfa cultivar43)(45.
The GR in summer was significantly higher than the GR in autumn and spring, when radiation and temperature are higher and consequently, greater photosynthesis40 and radiation use efficiency (RUE) (g DM MJ- 1 incident radiation)46 are achieved. The lowest GR occurred in winter, influenced by the interactions due to environmental conditions (low temperature and radiation) and defoliation intensity, affecting carbohydrates partition45 and consequently reducing the GR.
The evolution of the GR according to pasture age and season, showed a wide variation of the values within the seasons for all the years analyzed (Table 2). On an annual basis, the calculated variation coefficients (VC) were closely comparable to those obtained in other producing areas in Argentina47.
The average GR during the first year of growth in winter and summer, as well as in the third year in autumn and winter showed variation coefficients greater than 45 %(Table 2). This greater variability between years could indicate a stronger influence of meteorological factors in the GR during these growth periods, given that crop management was similar between years in the plots.
Relation between growth rate and climatic variables by season
Autumn
In autumn, GR was similar to spring (Table 2). As pasture ages, the dispersion of GR mean values increased at this time of the year, shown by VC above 50 % (Table 2). Air temperature in autumn decreased significantly compared to summer (Figure 1), reducing alfalfa GR, however, the temperature range remained within favorable growth thresholds, between 8 and 20 °C43)(46 with no extreme temperatures (Figure 1). Despite this global relationship with temperature, no direct relationship between GR and autumn temperature between years (r2 = 0.05) was found to explain the variability of GR between years. Furthermore, the results of CART analysis allowed us to discriminate four GR response groups concerning cWB (Table 3a), integrating in this index the relationship between rainfall and atmospheric demand.
In autumn, the maximum GR (56.3 kg DM ha-1 day-1) was associated with intermediate water conditions; cWB in autumn ranged between -18 and 140 mm, with an average of 48 mm in 16 out of 46 periods analyzed. Water conditions in autumn higher or lower than these cWB values, significantly reduced alfalfa GR, in 30 out of 46 autumn periods (Table 3a). This relation was found regardless of alfalfa age. This higher GR was associated with cWB values close to 48 mm (Table 3a) when the crop’s accumulated evapotranspiration was slightly lower than rainfall. On the other hand, the lowest GR (23.8 kg DM ha-1 day-1) were found during autumn periods with significant high rainfall compared to ETo, creating excess water situations that did not enhance alfalfa growth. Growth periods with extremely positive cWB would indicate agronomic conditions for possible soil waterlogging, which induce hypoxia, influenced by the internal soil drainage where planting takes place. Waterlogging decreases the growth of shoots and roots, increases mortality and decreases alfalfa vigor48)(49)(50; on the other hand, higher temperatures increment damage produced by root asphyxia49, combining lack of oxygen at root level, modulation reduction and a partial increase in plant metabolism51. In general, these climatic conditions also favor fungal diseases development in foliage and root, which usually reduce yield and compromise its persistence52)(53. The dispersion of the rainfall in autumn was higher in relation to other growth seasons (Figure 1), which could explain the very positive cWB frequency and autumn GR variability.
Winter
Given that E_CHANA and CRIOU cultivars are of intermediate latency, growth stops in response to short and cold days, then lower GR are expected in winter than in other seasons. GR in winter was the lowest among all seasons (Table 2). Analyzing GR in the set of two alfalfa cultivars through CART, showed that high Tmin in winter (> 8.5 °C) produced lower GR (Table 3b), possibly due to the increase in respiration rate and decrease in carbohydrates availability54. This relationship is not linear and enabled the identification of three groups according to the relationship between Tmin and GR in winter. The GR was highest in winter (26.1 kg ha-1 day-1) when air Tmin ranged between 7.8 and 8.5°C; outside this range, GR decreased significantly. Although Tmin showed greater capacity of group separation in GR through CART; Tmin, Tmax, and Tmed were autocorrelated in this period for the same groups (Table 3b).
Daily average radiation during winter was not significantly different between GR groups, or from relative parameters to soil water conditions, shown by high values of ETa:ETm ratio, suggesting GR are independent of environmental factors related to water dynamics in winter (Table 3b).
Spring
The GR in spring was higher than in winter, and there was a relatively low average variation of the GR between years (Table 2). The GR increased when Rs_RADIAT was greater than 16.03 MJ m-2 day-1 and when the air temperature rose (Table 3c); causing a RUE increase46, thus showing better conditions for photosynthesis and nodulation40 and consequently more biomass. As temperature increased the GR increased too, producing a direct effect on alfalfa RUE46. Incident radiation and water well-being (ETa:ETm > 0.6) were significantly associated with GR during spring (Table 3c), where good water availability conditions interacted with radiation and temperature rise, favoring a GR increase, possibly through the rise of the stomatal conductance of the canopy44.
The relation between water well-being and alfalfa GR was different according to average incident radiation during spring (Table 3c). Three GR groups discriminated by CART showed three different relations between GR and ETa:ETm. For Rs_RADIAT lower than 16.03 MJ m-2 day-1, there was no direct relationship between GR and water well-being (r2 = 0.0018, p < 0.6, n = 13); when Rs_RADIAT was between 16.3 and 50.54 MJ m-2 day-1, GR slightly began to increase when ETa:ETm increased (r2= 0.24, p < 0.039, n = 34 ). Finally, when Rs_RADIAT was greater than 50.54 MJ m-2 day-1 the possibility of achieving the highest GR was related to high values of ETa:ETm (r2 = 0.91, p < 0.001, n=9). These three scenarios defined spring periods where radiation and temperature can limit GR. On the other hand, when the highest GR were expected due to maximum radiation (> 50.54 MJ m-2 day-1), GR were limited by rainfall when the maximum atmospheric demand was not satisfied (ETa:ETm). Therefore, GR was constrained in some years, and consequently DM yield.
In Colonia, where these alfalfa trials were planted, periods with water deficit in spring often occurred, shown by the low ratio values between ETa and ETm, especially in those springs where radiation was not limiting. This suggests that some measures, such as supplementary irrigation, could increase alfalfa GR potential through the reduction of water deficit, satisfying the potential demand and achieving the greatest radiation use.
Although alfalfa is considered a relatively drought-tolerant species that can use up to 65-70 % of the available soil water (ASW) before transpiration decreases2, lower soil water content would decrease transpiration and consequently the ETa:ETm ratio34)(55. Alfalfa showed to have a high growth recovering ability under small water deficit conditions (seven days), but a limiting recovery with a deficit of 14 to 21 days56. Equation (1) implies that small variations in ETa:ETm ratio, significantly affect GR in spring.
Air temperature was mainly associated with daily Rs_RADIAT; the lowest average temperature was associated with lower radiations (< 16.03 MJ m-2 day-1), which decreased GR (Table 3c), regardless of the soil hydric state, reducing RUE46. As radiation and temperature increased, GR rise was more related to the hydric regime and water availability, in agreement with Brown and others46 and Ward and Micin10. The effect of temperature increase was reflected in a higher RUE46. In Colonia, growth periods with high daily radiation in spring (> 20.54 MJ m-2 day-1) were directly associated with periods in which air Tmax average was the highest (22.5 °C), and rainfall was low compared to the atmospheric demand (ETo).
Summer
Alfalfa growth was higher during summer (Table 2), when radiation and mostly temperature, were between the optimum range (20-25 °C) for temperate species growth44)(57. These two variables alone did not allow to discriminate GR variability found between years during summer (r2 = 0.20 and 0.01 respectively, p < 0.60, n = 48). However, GR was discriminated in significantly different groups when analyzing the ratio between ETa and ETm (Table 3d). The GR in summer increased (3x) when ETa:ETm ratio was greater than 0.42 (Table 3d), allowing to discriminate GR in three groups according to its relationship with water stress level in summer.
Summer growth periods that had high ETa:ETm values (> 0.88) and the highest GR (101.5 kg ha-1 day-1) were relatively few in comparison to the total number of summer periods analyzed (two periods out of 48). The ETa:ETm average ratio in summer was 0.66 (36 out of 48) (Table 3d)(GR = 60.1 kg ha-1 day-1), value well below the optimum to achieve high potential yields34)(58. This shows suboptimal conditions for alfalfa growth in most evaluated years, in which neither radiation nor temperature were the limiting factor, but water, in maximum yield. During the analyzed period, the value of 0.88 was exceeded only in few years, in which the GR was the highest.
Lower GR in summer is also associated with higher Tmax (Table 3d), which corresponded to lower ETa:ETm ratio periods and greater water deficit. The greater the soil water stress was (ETa:ETm ratio decrease), the lower GR was during summer periods (GR = 84.79 ETa:ETm + 2.75, r = 0.72, p < 0.01, n = 48). This ratio slope showed that small ETa:ETm ratio increase indicates significant alfalfa GR increase. The closer ETa:ETm ratio is to one, the greater transpiration is compared to the atmospheric demand adjusted to the phenological stage (between cuttings) the lower its water stress and the higher its productivity34. The individual analysis of evapotranspiration and rainfall showed no sources of discrimination of productive alfalfa behavior.
During summer in Colonia, variables associated with important GR changes were mainly related to atmospheric demand and water dynamics (Table 3d) under equal solar radiation conditions. Radiation and temperature directly affect RUE, and consequently growth46)(59; however, the effect of these variables on GR is directly limited to plant water balance60)(61. These observations coincide with Henry62 revealing that in our country there are periods of climatic water deficiency, especially from December to March. Depending on the storage capacity of soils and pasture tolerance, this period could entail plant water deficit and could extend in some years from spring (October) to autumn (April). Water stress is the most limiting environmental factor behind productivity and plant growth stability63, and the impact on yield due to water stress is the main limiting factor behind summer crop production in Uruguay64. On the other hand, it is the main reason for crop variability yields between years64)(65)(66.
Conclusions
Alfalfa GR presented significant differences between seasons of the year and according to pasture age in the southwest of Uruguay. Variability of GR according to pasture age, expressed through changes in the variation coefficient, could be associated to a greater meteorological variables influence in the most vulnerable pasture stages: the initial alfalfa installation stage (first year of life), and the stage when it can begin to affect pasture persistence by plant death, root and crown diseases, insect attack, etc. (third year of life).
Climatic variables affected alfalfa GR differentially according to the season of the year considered. In autumn, water deficits and excesses were limiting to GR, showing maximum GR values between -18 and 140 mm of cWB, which correspond to 35 % of the evaluated autumn growth periods. In winter, GR was favored by the increase of minimum temperatures in periods of no water excess. The GR was maximum in winter (26.1 kg ha-1 day-1) when the minimum air temperature ranged between 7.8 and 8.5 °C; outside this range, GR decreased significantly, which only occurred in 36 % of the winters considered. In spring and according to pasture age, the higher the solar radiation and ETa:ETm (lower water deficit), the higher alfalfa GR. Summer was the greatest growth period of this summer legume. The greater the water well-being of the pasture, denoted by the ETa:ETm ratio closer to one, the higher GR, provided that the maximum air temperatures were not limiting. Only in 6 % of the evaluated summers, ETa:ETm ratio exceeded 0.80; showing GR restrictions due to water deficiency during growth. These water restrictions could be alleviated by the incorporation of management measures that improve water use efficiency, more drought tolerant cultivars, or supplementary irrigation technologies during summer; especially during evident water deficit years. Consequently, growth rates would be high and therefore alfalfa forage yield maximized.