elocation-id: e3202
The quantity and quality of potato tubers are related to irrigation. We studied the effect of three levels of water deficit (L= low, M= medium, and H= high) on the yield and quality of three potato varieties (Fiana, Atlantic, and Ágata), induced by varying the parameter α3 of the maximum allowable deficit (MAD) function of a comprehensive irrigation programming model. The work was carried out in the north of Sinaloa, Mexico (autumn-winter cycle 2015-2016), using a design in divided plots with two factors and four repetitions, where the large plots were the levels of water deficit and the sub-plots were the varieties. There was a significant effect (p= 95%) on yield (Y) as it decreased by 3.65 t ha-1 in H compared to L. In the variety factor, Fiana and Ágata were the best in quantity and quality, with Fiana standing out with fewer noncommercial (NCO) tubers. In the water deficit-variety interaction, Fiana was the best, with more tubers of first and second quality (Q1 and Q2) in L. In M, Ágata was the best in Y but the lowest in Q1 and highest in deformed (DEF) and NCO. In H, the behavior was similar; Ágata was the best but with less Q1 and Q2 and a significant presence of DEF and NCO. A 29.5, 28.8, and 28.4 cm water requirement was estimated in L, M, and H, respectively. The behavior of Fiana suggests managing it under drip irrigation, while Ágata under gravity irrigation and limited water conditions.
Solanum tuberosum L., comprehensive model, growing degree-days (GDD), water use efficiency.
In Mexico, 62 000 ha of potato (Solanum tuberosum L.) are established annually, with a production of 1.7 million tonnes and an average yield of 27 t ha-1 (Onamu et al., 2015; SIAP, 2021); in addition, 8 700 producers participate directly, benefiting 20 000 families, generating more than 17 500 direct jobs and 51 600 indirect jobs (ASERCA, 2018). Sinaloa is the largest potato-producing state in Mexico, establishing annually in the autumn-winter cycle (AW), around 14 000 ha under gravity and sprinkler irrigation, which represented 22% of the national area, with an average yield of 30 t ha-1 (CONAGUA, 2018).
The quality and quantity of tubers are related to good irrigation programming due to the high sensitivity of the crop to water stress (Wagg et al., 2021). Steduto et al. (2012) reported that the crop’s accumulated water demand fluctuates from 350 to 650 mm, depending on sowing date, climate, soil moisture, and variety, with water productivity of 4 to 11 kg m-3. The same authors suggested that in order not to affect the yield, water deficit should be avoided during the stages of stolon formation and tuber development, and the usable moisture (UM) of the soil should not be depleted by more than 30 to 50%.
During two agricultural cycles, Camargo et al. (2015) assessed the response of potato crop yield and water productivity to four irrigation treatments (60, 80, 100, and 120% of the water requirement), finding statistically significant differences in the treatment where 60% was applied in the first year, while in the second, the treatments with high levels of water supply (100 and 120%) were also different compared to those of low supply (60 and 80%). Water productivity fluctuated from 8.6 to 11.6 kg m-3 in the first year and 7.1 to 8.4 kg m-3 in the second.
Regarding quality, Steduto et al. (2012); Da Silva et al. (2019) recommended that to reduce the presence of deformed tubers, the crop should not be stressed during the tuber development phase, heavy irrigation after a period of stress at this stage also generates cracking of tubers and black hearts. Jha et al. (2017) evaluated the effect of water stress induced by the concentration of salts in irrigation water on yield and quality under drip and gravity irrigation, finding a gradual reduction of the variables plant height, number of tubers, and yield as stress increased, mainly under gravity irrigation.
This effect was also observed in tuber quality as starch concentration decreased, and tubers of less weight and size increased with increasing salt concentrations. Poor irrigation programming is also related to the presence of diseases that reduce quality, mainly caused by bacteria and fungi (Rich, 1983).
The use of indicators of the level of soil moisture (θ) is one of the most recommended techniques to schedule irrigation; the moment is defined when the moisture reaches a minimum allowable or critical value (θc), dependent on the factor known as maximum allowable deficit (MAD); in this way, soil moisture fluctuates between field capacity (θfc) and θc (Mahima et al., 2018). Another widely used irrigation programming method is water balance, which is based on making a rigorous balance of water inputs and outputs in the root zone that affect soil moisture dynamics.
In the former are irrigation (I); effective precipitation (EP); capillary rise (CR); and runoff (Ru), while in the latter are crop evapotranspiration (ETc); percolation (P); and runoff or drainage (Ru). Because the ETc depends on environmental conditions, soil moisture, species, and stage of crop development, it is difficult to estimate it using fixed values for the duration of each stage, so it is pertinent to incorporate the concept of Growing Degrees-Days (GDD) in the monitoring of programming variables (Ojeda et al., 2006; Ojeda et al., 2015).
This approach is becoming increasingly important due to the possible intensification of climate variability and change in Mexico’s various agricultural regions. The Secretariat of Environment and Natural Resources (SEMARNAT, for its acronym in Spanish), an agency of the Mexican government, reported that according to records from 1971-2008 in the country, the historical average annual temperature was 20.9 °C and since 1990, most of the average annual temperatures have been higher than this value.
Only in the period 1999-2008 the increase was 0.7 °C (SEMARNAT, 2012). Lobato and Altamirano (2017) reported an increase in the average annual temperature greater than 2 °C for northern Sinaloa, where the present study was carried out, taking as reference the base period 1961-1990. For the programming and scheduling of irrigation under the GDD concept, computer systems have been developed, which are applied on a large scale to better couple the water demands of crops to their phenological development, such as the Irrimodel platform (Sifuentes and Macías, 2015).
Nonetheless, most studies using these irrigation programming tools have focused on reporting the effect of water stress on yield (Erdem et al., 2006), and there is limited information on these effects on tuber quality for currently grown varieties.
Given this demand for knowledge, the objective of the present research work was to evaluate the response of the yield and quality of the tubers of the three potato varieties with the largest sown area in northern Sinaloa (Fiana, Ágata, and Atlantic) to the water deficit controlled through parameter ∝3 of the MAD function of a comprehensive-dynamic model of irrigation programming based on the GDD concept, using the Irrimodel software.
The study was conducted in the northern part of the state of Sinaloa, Mexico, during the autumn-winter (AW) cycle 2015-2016, in the central part of the DR-075 Río Fuerte at an altitude of 12 m. The region has an average annual temperature of 24 to 26 °C and an average accumulated annual rainfall of 360 mm, concentrated from July to September. The predominant soils are clayey (C) and clayey loam (CL) in texture, poor in organic matter (OM), flat slope, and without problems of soluble salts (INIFAP, 2010).
An experimental lot of 2 500 m2 was established in the Valle del Fuerte Experimental Field (CEVAF, for its acronym in Spanish) of the National Institute of Forestry, Agricultural, and Livestock Research (INIFAP, for its acronym in Spanish), located at coordinates 25° 45’ 49” north latitude and -108° 51’ 41” west longitude. The soil of the lot was clayey, deep, and without salt problems, typical of the region, and there were no significant rains during the growing period of the crop.
The soil water deficit was managed through parameter ∝3, included in the MAD function of the comprehensive irrigation programming model indicated in Table 1, proposed specifically for this crop by Ojeda et al. (2004) and executed automatically through the Irrimodel software (Sifuentes and Macías, 2015).
Kco= represents the crop coefficient for the first phenological stage of the crop, which depends essentially on soil evaporation; Kmax= is the maximum value of Kc during its development; XKmax= corresponds to an auxiliary variable defined by the accumulated growing degree-days (AGDDs) when the crop reached its maximum; Kc, ∝1= is a regression parameter obtained by adjusting experimental data from the study area; erfc= is the complementary function of the error; and x= is an auxiliary variable calculated with equation (1):
.
Where: AGDD= is the accumulated value of GDD from sowing to the calculation date, and ∝0= is the AGDDs required until the end of the phenological cycle. RDo and RDmax represent the sowing depth and maximum root depth, respectively; the ∝2 value was empirically adjusted from an approximate value to 2/3 of the AGDD value where the potato reaches the maximum root depth.
Three water deficit conditions were generated, using three values for parameter ∝3, which vary according to the irrigation system, as reported by Ojeda et al. (2004): 0.15 for low water deficit (L), 0.3 for medium deficit (M) and 0.45 for high deficit (H), in three commercial varieties currently grown in the study area: Atlantic (V1); Ágata (V2); and Fiana (V3).
The ∝4 value was assumed to be 0.1, suggested for potato crops (Ojeda et al., 2004). The experimental design used was in divided plots (Soltero-Díaz et al., 2010; Badii et al., 2007), with the applied level of moisture or water deficit being the large plot and the three varieties distributed randomly being the sub-plots, generating 36 experimental plots formed of four beds 0.9 m wide by 15 m long (54 m2).
The sowing was carried out dry on October 24, 2015, manually, depositing 4 to 5 seed tubers in the center of each bed at a depth of 0.2 m to have a sowing density of 45 000 seeds ha-1. After sowing, gravity germination irrigation was applied to bring soil moisture to θfc in the upper meter of the profile. Fertilization, crop management, and phytosanitary control were carried out according to the technological package recommended for the area (INIFAP, 2010).
The programming of the irrigations in each treatment, in terms of timeliness (interval) and quantity (sheet), was carried out using the water balance method through the Irrimodel Software (Sifuentes and Macías, 2015), which uses the comprehensive model described in Table 1.
This Software is executed through a computer platform via the internet and also requires the following information: 1) usable moisture by soil stratum (UM) calculated as the difference between θfc and permanent wilting point (θpwp); 2) crop phenology; 3) climatic variables: mean temperature (Tm) and reference evapotranspiration (ETo-FAO) (Allen et al., 1998; Mahima, 2018) and rainfall (PP), acquired in ‘quasi-real’ time from an automated agroclimatic station; and 4) application of irrigation. It was assumed that the root zone reaches moisture at field capacity (θfc) at the end of each irrigation.
To have greater control of the application and evaluation of the irrigations, a drip system was used with a tape of 18 mm in diameter and drippers separated at 0.3 m with a nominal expenditure of 1 L h-1, installed in the upper part to the center of each sowing bed. The irrigations were evaluated in terms of application efficiency (AE) and water productivity (WP), expressed in % and kg m-3, respectively.
The estimation of AE was made with equation (2):
Where: Sn= represents the net sheet or net irrigation requirement (cm) determined with the model, and Sb= represents the applied sheet (cm) measured with the irrigation system (Rendón et al., 2012). The calculation of WP was performed with equation (3):
Where: Y= represents the yield (t ha-1) and VA= the volume of water applied per hectare (m3 ha-1).
The phenology of the crop was monitored weekly using an adaptation of the methodology reported by Jefferies and Lawson (1991), considering the following phases of development: 1) sprouting; 2) emergence; 3) start of stolonization; 4) stolon elongation; 5) start of tuberization (diameter >1 cm); 6) tuber development (diameter from 3 to 5 cm); 7) developed tubers (diameter >5 cm); 8) foliage killing (complete skin) (diameter from 7 to 8 cm); and 9) harvest.
GDDs were accumulated for each development phase using the mean daily air temperature (Ta) from an automated weather station (AWS) located 100 m from the experimental lot. GDDs were calculated with the method reported by Ojeda et al. (2004) in potato crop studies in Sinaloa, using 2 and 29 °C as minimum and maximum development threshold temperatures.
Soil moisture was also monitored biweekly in the two upper strata (0-30 and 30-60 cm) with a Spectrum 300 time-domain reflectometry (TDR)-type portable moisture sensor manufactured in Illinois, USA, and calibrated for the study area.
The killing of foliage was carried out using contact herbicide Velquat 25 (ai. Paraquat) on 10/02/2016 at 1 754 AGDDs, while the harvest was carried out on 11/03/2016 at 2 282 AGDDs once the tuber skin had the appropriate consistency. The tubers were extracted from the soil mechanically by tractor, using 18-inch moldboards mounted on a horizontal support bar. After the killing of foliage, in each experimental plot, the tubers that remained within an area of 9 m2 (useful plot) formed of the two central beds 5 m long (2 x 0.9 x 5 m) were collected and weighed manually.
The samples obtained were determined the variables Y, relative yield (RY), and quality; the first calculated with equation (4):
Where: Ws= represents the weight of the sample expressed in kg, Y= in t ha-1 and the constants 10 and 9 are constants for unit conversion. RY of each variety was calculated with respect to the yield of the treatment with low water deficit, expressed as a percentage.
To analyze the quality of tubers, the different qualities were separated considering the equatorial diameter (mm), according to the criteria of the area (Table 2) and it was expressed in terms of percentage with respect to the total yield. Camargo et al. (2015) used a similar methodology grouping three categories, specifying that tubers of 40 to 80 mm in diameter are the most commercialized, which coincides with the classification of the study area. In addition to the equatorial diameter, deformed (DEF) and diseased (NCO) tubers were separated.
| Quality | Code | Diameter (mm) |
|---|---|---|
| First | Q1 | 71-80 |
| Second | Q2 | 51-70 |
| Third | Q3 | 41-50 |
| Fourth | Q4 | 31-40 |
| Fifth | Q5 | 21-30 |
| Giant | GIA | >80 |
| Deformed | DEF | - |
| Noncommercial | NCO | - |
For the statistical analysis, the SAS online version (www.sas.com/es-mx/software/on-demand-for-academics.htm) computer program was used, with which an analysis of variance of the variables yield (Y) and quality (Q) was carried out, considering the simple effects and the interaction of the factors water stress and varieties, in addition to a comparison of means with Tukey’s test (p≤ 0.05). The model used is shown in equation (6) (Fernández et al., 2010).
Where: Yijr - μ= represents the deviation of an observation from the global mean (μ); ρr= are the deviations caused by the blocks; αi= is the main factor; βj= is the subordinate factor; (αβ)ij= is the interaction between factors; and εjr and εijr= are random errors. For the quality variable, an analysis of the distribution of each classification by variety was carried out to know the effect of water deficit on this important variable.
Table 3 shows the analysis of variance of variables Y and Q of the two factors studied. The analysis shows a significant effect (Tukey, p≤ 0.05) of the water deficit factor only for the variable Y; the results showed a reduction of up to 3.65 t ha-1 of the variable Y in the treatment with high water deficit compared to the low deficit, which is consistent with what was reported by Steduto et al. (2012); Mahima et al. (2018).
Regarding the variety factor, a significant effect was also observed in the variables Y, Q1, Q2, Q4, and Q5; Fiana and Ágata obtained similar yields, with 26.14 and 26.45 t ha-1 respectively, while Atlantic obtained 21.22 t ha-1, which represents a difference of 4.92 t ha-1. In this factor, Q1 and Q2 were higher with Fiana, followed by Atlantic and Ágata; Fiana also presented fewer noncommercial tubers.
In the water deficit-variety interaction, significance was found only for the Q5 variable; however, important effects were observed: at the low level of water deficit, Fiana obtained the highest yield (30.29 t ha-1), with more Q1 and Q2, followed by Ágata and Atlantic, with the latter showing fewer DEF and NCO tubers. At the medium level, Ágata was the best in yield but the lowest in Q1 and highest in DEF and NCO. The behavior was similar at the high level; Ágata was the best in Y but with less Q1 and Q2 and an important presence of DEF and NCO.
Figure 1 shows the effect of the three levels of water deficit on the RY of each variety (stress level x variety interaction). The results indicate that with the high level of water stress (∝3 = 0.45), the Fiana variety had the highest RY drop, with 26.3%, followed by Atlantic with 5.1% and Ágata with only 2.8%, compared to the low level, and tendency not to affect the yield at the medium level. This indicates that the Fiana variety is the most sensitive to water deficit, while Ágata is the most tolerant, Atlantic can be considered of intermediate sensitivity.
Studies such as those by Da Silva et al. (2019) reported a similar behavior regarding the Ágata variety, concluding that the best management conditions were obtained when applying intermediate irrigation sheets (27 mm between irrigations).
An important effect of water deficit on this variable was found in each variety (Figure 2), as reported by Jha et al. (2017). At the low deficit level, the Atlantic variety produced the largest number of Q1 and Q2 tubers, with 15 and 18% of total production, respectively. The high and medium deficits generated the highest values of Q3 tubers, with 30 and 23%, while the highest value of Q4 tubers was obtained with the medium level of water deficit.
The Q5 tubers remained around 10% in the three levels of water deficit due to their little development due to the advance to harvest induced through the defoliant. Similarly, for the Ágata variety, the highest percentages of Q1 and Q2 tubers were generated with the low level of water deficit (5 and 12%, respectively), maintaining similar values in the Q3 tubers in the three levels of deficit (18 to 19%).
A significant increase in Q4 and Q5 tubers was analyzed in the three deficit levels, mainly with the high value, which induced the production of 38% of Q5 tubers. The high values of Q4 and Q5 tubers could also be due to the slower development of this variety compared to the others. This confirms that the values of parameter ∝3 are associated with the irrigation system; the values vary from 0.15 to 0.3 for drip irrigation, from 0.4 to 0.55 for sprinkler irrigation, and they were greater than 0.6 for gravity irrigation (Ojeda et al., 2004).
Regarding the Fiana variety, the highest values of Q1 tubers were similar in the three levels of water deficit, with percentages of 11 to 13%, while the highest values of Q2 tubers occurred with the low level of deficit, with 22%. The Q3, Q4, and Q5 tubers presented similar increases in the three stress levels; although the higher percentage was more frequent with the medium level of water deficit in the three qualities, a behavior similar to the other varieties.
With regard to DEF tubers, a direct relationship with water deficit was examined in all three varieties; that is, when this increases, the proportion of this class of tubers increases, as mentioned by Steduto et al. (2012), it reached values close to 10% of production with the high and medium level of deficit.
In the NCO, a tendency towards an increase with the low level of water deficit was found in the three varieties, which may be associated with a more humid environment generated by the low permeability of the soil, Wagg et al. (2021) suggest high levels of moisture in well-drained soils, since heavy and frequent irrigation in soils with low permeability favors the development of pathogens, mainly fungi and bacteria (Rich, 1983). With the medium and high deficit levels, the tendency was towards decreasing NCO, so in practice, the different stress management scenarios should be analyzed according to the objectives of each producer.
Table 4 presents a summary of applied supplemental irrigations generated with the three levels of water deficit. It was perceived that 14 supplemental irrigations with a maximum interval of 22 days and a minimum of five days were generated for the low level; eight supplemental irrigations with a maximum interval of 30 days and a minimum of eight days were applied with the medium level of deficit.
On the other hand, with the high level, the maximum and minimum intervals were 36 and 13 days, respectively. Regarding the irrigation requirements or Sn, the low level had a value of 29.5 cm, higher than the medium and high levels, due to the lower level of stress to which the crop was subjected, which affects evapotranspiration.
A higher value of AE was also observed at the low and medium levels, with values of 90.6 and 90%, respectively, while at the high level of stress, it was 81.5% because of the reduction of Sn due to the high stress considered in the ETc calculations and the greater Sb derived from a longer irrigation time used in the formation of the wet bulb.
Regarding WP, it was found that with the low level of water deficit, the varieties Atlantic and Fiana obtained the highest values, with 6.7 and 9.32 kg m-3, respectively, while for Ágata, the highest values were similar at the low and medium levels, with 8.03 and 8.7, which indicates some resistance to water stress. In all cases, the values are in the range reported by Erdem et al. (2006); Steduto et al. (2012); Mahima et al. (2018).
The Fiana variety was the most sensitive to water deficit from the point of view of yield, followed by Atlantic and Ágata. Applying irrigation intervals, sheets, and AE recommended by the Irrimodel Software was very convenient in all cases. The effect of water deficit on tuber size was different for each variety due to their specific development pattern and tolerance to water stress.
As the water deficit increased (∝3 high), associated with a longer irrigation interval, the proportion of deformed tubers increased, while diseased and noncommercial tubers increased with the low deficit (∝3 low), so conditions of excess moisture should be avoided, mainly in heavy soils and cold climatic conditions.
The application of the Irrimodel software represents an excellent automated tool for the controlled management of water deficit according to the irrigation system, which translates into the improvement of potato yield and quality for the producing areas of Mexico. For a greater effect of water deficit on yield, it is suggested to use ∝3 values of 0.15, 0.45, and 0.6 for low, medium, and high deficits, which are associated with drip, sprinkler, and gravity irrigation systems, respectively.
ASERCA. 2018. Agencia de Servicios a la Comercialización y Desarrollo de Mercados Agropecuarios. http://www.infoaserca.gob.mx/claridades/revistas/057/ca057.pdfwww.infoaserca.gob.mx/claridades/revistas/057/ca057.pdf.
CONAGUA. 2018. Comisión Nacional del Agua. Estadísticas hidroagrícolas de los distritos de riego. https://www.gob.mx/conagua/documentos/estadisticas-agricolas-de-los-distritos-de-riego.
Da Silva, E. F.; Jadoski, S. O. and Martinkoski, L. 2019. Water depletion depth for irrigation of potato cultivar Ágata. Engenharia agrícola, Jaboticabal. 39(4):434-443. Doi: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v39n4p434-443/2019
INIFAP. 2010. Instituto Nacional de Investigaciones Forestales. Agrícolas y Pecuarias. Informe final del proyecto: evaluación del cambio climático en la productividad de la agricultura de riego y temporal del estado de Sinaloa’. Reporte técnico del proyecto. Instituto Mexicano de Tecnología del Agua (IMTA). 15-35 pp.
Jha, G.; Choudhary, O. P. and Sharda, R. 2017. Comparative effects of saline water on yield and quality of potato under drip and furrow irrigation. Cogent Food & Agriculture 3(1) https://doi.org/10.1080/23311932.2017.1369345.
Ojeda, B. W.; Flores, V. J.; Íñiguez, C. M.; Unland, W. H.; Carrillo, G. M.; Rojano, A. A.; Flores, G. H.; Macías, C. J.; Flores, M. H.; López, L. R.; Sifuentes, I. E. y Ruiz, G. A. 2015. Calendarización del riego: teoría y práctica. Libro, 2da. Ed. Jiutepec, Morelos, México. Instituto Mexicano de Tecnología del Agua. 3-307 pp.
SIAP. 2021. Servicio de Información Agroalimentaria y Pesquera. Avance de siembras y cosechas resumen por estado; producción anual agrícola. https://www.gob. mx/siap/acciones-y-programas/produccion-agricola-33119.
Wagg, C.; Hann, S.; Kupriyanovich, Y. and Li, S. 2021. Timing of short period water stress determines potato plant growth, yield and tuber quality. Agricultural Water Management. 247: https://doi.org/10.1016/j.agwat.2020.106731