#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Agronomic performance of lettuce cultivars submitted to different irrigation depths


Authors: Claudinei Martins Guimarães aff001;  Fernando França da Cunha aff001;  Francisco Charles dos Santos Silva aff001;  Edcássio Dias Araújo aff001;  Aline Baldez Felismino Guimarães aff002;  Everardo Chartuni Mantovani aff001;  Derly José Henriques da Silva aff003
Authors place of work: Department of Agricultural Engineering, Center of Agricultural Sciences, Federal University of Viçosa, Viçosa-MG, Brazil aff001;  Department of Domestic Economy, Center of Human Sciences, Federal University of Viçosa, Viçosa-MG, Brazil aff002;  Department of Phytotechnics, Center of Agricultural Sciences, Federal University of Viçosa, Viçosa-MG, Brazil aff003
Published in the journal: PLoS ONE 14(12)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0224264

Summary

The use of cultivars adapted to the climate and soil conditions associated with adequate irrigation supply maximizes lettuce agronomic performance. The aim of this work was to evaluate the agronomic performance of four lettuce cultivars submitted to five different drip irrigation depths under a protected environment in Viçosa-MG, Brazil. A randomized block design was applied in a split plot scheme with four replications, and several agronomic characteristics were evaluated by analysis of variance, Tukey range tests, regression and principal component analysis. A higher chlorophyll concentration in the Raider Plus cultivar promoted the production of more leaves, leading to a higher phytomass. The Luara cultivar presented a higher number of commercial leaves per plant than the other cultivars, regardless of the irrigation depth, reflected in a larger diameter and volume of the aerial part of the plants. The Raider Plus and Luara cultivars presented better root development than that of the other cultivars, reducing the effect of plant water stress even under lower irrigation depth conditions. Although higher water productivity (WP) was observed for the lowest irrigation depth (50% of ETc), important variables reached the maximum values at depths higher than 100% water replenishment. Therefore, Raider Plus and Luara cultivars with an irrigation depth of 110% of crop evapotranspiration provided better commercial lettuce quality and are recommended for lettuce cultivation in the research region and under conditions similar to those of the present study.

Keywords:

Leaves – Crops – Chlorophyll – Vegetables – Agricultural irrigation – Lettuce – Vegetable crops – Drip irrigation

Introduction

The cultivated lettuce (Lactuca sativa L.) of the Asteraceae family is one of the most popular and consumed vegetables in the world [1, 2] and the main leafy vegetable marketed in Brazil [3]. Generally, it is consumed in salads and represents a good source of minerals [4]. Grown all over Brazil, traditionally on small family farms, lettuce has great economic and social importance [5], strengthening self-reliance and improving the subsistence of small farmers [6].

The choice of cultivars that are more adapted to the climate and soil type of a region allows an increase in crop yield. Thus, due to its origin of temperate climate, it was necessary to develop new lettuce cultivars that were more adapted to tropical conditions over the years [7]. In addition, each region of tropical climate has very specific local soil and climatic characteristics to be considered, which justifies the study of cultivars adapted to each region. Currently, the main types of lettuce grown in Brazil are crisp, American, smooth and others (red, Mimosa, and Roman, among others), corresponding to 70%, 15%, 10% and 5% of the Brazilian market, respectively [3].

Lettuce cultivation is usually performed in a protected environment or in the open field [8], but the latter does not include the quantification of the necessary irrigation depth [9]. Therefore, the use of a protected environment, considered intensive [10], has become more frequent because it protects crops from climatic adversities, pests and diseases, increasing the overall yield and product quality [11]. The use of irrigation enables successive crops throughout the year [9] and is obligatory in a protected environment to supply plant water needs.

Another important factor is an efficient irrigation method to fully or supplement the water needs of the crops. Second [12], the most commonly used method (conventional sprinkler) does not include aspects that are fundamental to the sustainability of irrigated agricultural activity, such as plant health and multiple uses of water. This fact favors the use of localized drip irrigation, successfully adopted for several crops, with great advantages, such as reduced energy consumption, facilitated and favored plant health control and greater efficiency in water productivity (WP).

Combined with an efficient irrigation method, proper irrigation water management enhances crop quality and productivity. Considering the intensive cultivation in a protected environment, the study of appropriate irrigation depths for each crop and region is even more important. According to [13], water is indispensable for production and therefore its lack or excess significantly affects the productivity of a crop, making its rational management indispensable to achieving the maximization of production.

In this sense, the aim of this work was to evaluate the agronomic performance of lettuce cultivars submitted to different drip irrigation depths, indicating suitable cultivars and soil water replacement depths for the cultivation of lettuce in a protected environment in the region of Viçosa-MG, Brazil.

Materials and methods

Experimental area

The present research was located at the Federal University of Viçosa (UFV)—Viçosa-MG, Brazil, at coordinates 20° 45' 14” S, 42° 52' 55” W and an altitude of 648 m above sea level (S1 Fig). The region climate of the type Cwb [14] presents an annual average temperature of 19.4°C and annual rainfall of approximately 1,200 mm.

Installation and conduct of the research

A protected environment (Fig 1) with a total area of 240 m2 (8 m wide, 30 m long and right foot of 3.2 m) was used, with the sides covered by polyethylene wire mesh (100% polyethylene, 25 mesh—1.0 × 1.0 mm aperture, with 10 wires cm-1 screen, and 25% shading) and a roof constructed of blue plastic (AV Blue, 120 microns, 78% and 67% transmission and diffusion of light, respectively).

Fig. 1.
Research installation location (a, b, c and d) and structure of one of the four blocks (e). C: lettuce cultivars. Viçosa-MG, 2017.

A eutrophic Red-Yellow Latosol (Table 1) was used. Soil preparation was carried out by plowing, harrowing and passing the enchanter. Chemical and physical water analyses (Table 1) were performed, and liming and chemical fertilization were performed based on the results of the soil chemical analysis and recommendations of [15].

Tab. 1. Soil physical-hydrological and chemical analyses, according to [16].
Soil physical-hydrological and chemical analyses, according to [<em class="ref">16</em>].

The sowing of lettuce was carried out in polystyrene trays, with transplanting of the seedlings to the beds (in small depressions spaced 0.25 m between rows and 0.25 m between plants) when they had four definitive leaves (30 days after sowing). Each sampling unit (Fig 2) had an area of 1 m2 (1 m in length and 1 m in width), consisting of 16 plants, and the four central plants were evaluated.

Daily variation of solar radiation, reference evapotranspiration, temperature and relative humidity for the three cycles of lettuce cultivation in Viçosa-MG in 2017.
Fig. 2. Daily variation of solar radiation, reference evapotranspiration, temperature and relative humidity for the three cycles of lettuce cultivation in Viçosa-MG in 2017.

The evaluations were made 39 days after seedling transplantation in three cultivation cycles (01/19 to 02/26/2017—cycle 1, 04/05 to 05/13/2017—cycle 2 and 07/25 to 01/09/2017—cycle 3). A randomized block design (RBD) was used with four replications in a split plot scheme. Five irrigation depths (50, 75, 100, 125 and 150% of the irrigation depth for crop evapotranspiration replacement—ETc) were implemented in the plots, and four lettuce cultivars (Raider Plus—Americana type; Luara—Lisa type; Vanda—Crespa type; Imperial—Roxa type) were planted in the subplots (Fig 1). Therefore, data from 80 observations (5 depths × 4 replicates × 4 plants sampled) were used to evaluate each analyzed characteristic, each cultivar and each cycle.

A drip irrigation system (Fig 1) was used, with lateral lines composed of drip tape (Amanco brand) 16 mm in diameter and 15 thousandths of inch thick. The spacing between the drip tapes was 0.50 m, which allowed for the irrigation of two rows of plants. The emitters (drippers), spaced 20 cm apart, operated at a service pressure of 98 kPa (~ 10 mca), applying an average flow of 1.8 L h-1.

Manual weeding was performed weekly, and there were no pests and diseases capable of causing damage to the lettuce quality or yield.

Irrigation management

The irrigation system was evaluated according to [19], modified by [20], collecting the dripper flows at eight points along the lateral line and in four lateral lines along the line of bypass. The distribution uniformity coefficient (DUC) was calculated by Eq 1:


where q25 is the average of the lowest quartile of flows (L h-1) and qm is the average flow (L h-1).

The application efficiency was considered 100%, and the average DUC was 93.1%. Thus, the irrigation efficiency (Ei) was 93.1% for all three cycles. To maintain pressure uniformity during the irrigation period (beyond this research), two lateral lines were always simultaneously kept open, such as the practical example outlined in S1 Table.

The real irrigation required (RIRLOC) and total irrigation required (TIRLOC) in the localized systems (mm) for 100% ETc treatment were estimated using Eqs 2 and 3, respectively, adapted from [21].



where ET0 is the reference evapotranspiration (mm day-1), KC is the coefficient of culture (dimensionless), KS is the coefficient of soil water depletion (dimensionless), KL is the location coefficient (dimensionless), and C is the constant relative to the elevation of the water table (mm).

ET0 was estimated according to Penman-Monteith FAO-56 [22] by Eq 4, with a wind speed equal to 0.2 m s-1 (recommended for protecting the environment).


where s is the slope of the saturation pressure curve versus temperature (kPa°C-1), RN is the radiation balance (MJ m2 day-1), G is the soil heat flux (MJ m2 day-1), γ is the psychrometric constant (kPa°C-1), t is the average air temperature (°C), U2 is the wind speed at 2 m in height (m s-1), es is the water vapor saturation pressure (hPa), and e is the current water vapor pressure (hPa).

Kc values were adopted from [23] and [24] for each stage (phase) of crop development: phase I (Kc = 0.65)—from transplant until the crop covers 10% of the land surface; phase II (Kc = 0.85)—from the end of the first phase to the effective full coverage of the ground surface; and phase III (Kc = 1.00)—from the end of the second phase until the maximum vegetative development, when lettuce is harvested. KS and KL were calculated according to [21] and [25] by Eqs 5 and 6, respectively:



where CWB is the current water depth (mm), TWC is the total water capacity (mm), and P is the highest percentage value between the wet and shaded area (%).

Over five days prior to beginning climate irrigation management, soil water content was monitored at a depth of 0–0.20 m using tensiometers based on the soil water retention curve θ = 0.1813 + (0,3443 − 0.1813)/[1 + (0.01|Ψm|)2]0.49 and by the direct stove method to ensure soil moisture at field capacity. Irrigation depths were differentiated for 7 days after seedling transplantation (acclimatization period of the seedlings to the soil).

Weather data

Meteorological data (Fig 2) were obtained through the automatic meteorological station E 4000 (IRRIPLUS) centrally installed in the protected environment (Fig 1). Daily reference evapotranspiration (ET0) was estimated from the data of solar radiation (Rs), maximum and minimum temperatures [and with these values the average air temperature (Taverage)] and relative humidity (RH).

The average daily Rs decreased over the crop cycles, ranging from 4.7 to 11.3 MJ m-2 in the first cycle, from 1.4 to 8.6 MJ m-2 in the second, and from 2.0 to 7.7 MJ m-2 in the third (Fig 2). This behavior influenced the decrease in Taverage and ET0. The Taverage varied from 22.7 to 26.1°C in the first cycle, from 16.5 to 24.3°C in the second cycle, and from 15.4 to 20.1°C in the third cycle. The average RH ranged from 36.1 to 62.6%, from 46.1 to 57.0% and from 33.6 to 61.6% in cycles 1, 2 and 3, respectively.

The total accumulated degrees per day for each crop cycle were obtained according to [26] and [27] and were determined as 546.9, 423.5 and 287.6°C day-1 for cycles 1, 2 and 3, respectively.

ET0 presented average values of 2.1 mm day-1 (ranging from 1.6 to 2.6 mm day-1) in the first cycle, 1.5 mm day-1 (from 1.0 to 1.8 mm day-1) in the second cycle, and 1.3 mm day-1 (from 1.0 to 1.7 mm day-1) in the third cycle. The ET0 that occurred during each cycle was used to determine the irrigation depth of 100% ETc applied in the treatments (Table 2).

Tab. 2. Real and total irrigation required in each treatment and lettuce-growing cycle.
Real and total irrigation required in each treatment and lettuce-growing cycle.

Due to the absence of rainfall inside the protected environment, a high frequency of irrigation (2-day irrigation shifts) was required, maintaining the soil with moisture near the field capacity and with little water requirement to reach the total storage at each irrigation period (Table 2).

Characteristics analyzed

  • Leaves: relative content of chlorophylls (SPAD index) A (Chl-A) and B (Chl-B) and total chlorophyll (Chl-T) were expressed as indirect in vivo measurements with a portable chlorophyll meter SPAD-502 (Soil Plant Analysis Development, Minolta Camera Co., Osaka, Japan) in SPAD values (higher values of SPAD indicate a higher content of chlorophyll [28]) from in intact and freshly mature leaves at four points on the leaf and two centimeters from the margin, always at 12 o'clock; the number of commercial leaves (NCL, units per plant) were expressed as the number of marketable leaves per plant; the number of noncommercial leaves (NNCL, units per plant) were expressed as the number of damaged or unlabeled leaves per plant; leaf length and width (LL and LW, cm) were measured as the greatest length and width of the more developed leaves; the leaf area index (LAI, dimensionless) was measured by the leaf discs method, using 30 discs of 2.0 cm in diameter (3.14 cm2) per plant taken from the basal, median and apical portions of three leaves, avoiding the sampling of the central vein, similar to [29]; the fresh leaf mass (FLM, g) was expressed as the mass of all leaves per plant; the dry leaf mass (DLM g) was expressed as the mass of all leaves per plant, dried in a forced ventilation oven at 65°C until reaching constant mass.

  • Roots: the fresh root mass (FRM, g) and dry root mass (DRM, g) per plant were obtained as outlined above for the fresh and dry leaf masses; the root length (RL, cm) was obtained for the longest root per plant, measured with a ruler; the root volume (RV, cm3) per plant was measured by the displacement volume of water inside a graduated cylinder.

  • Stems: the fresh stem mass (FSM, g) was determined as the fresh shoot mass without leaves; the stem length (SL, cm) was determined using a ruler; the stem average diameter (SMD, mm) was expressed as the average of three measurements (base, middle and tip) using a pachymeter; the stem volume (SV, cm3) was measured by the displacement volume of water inside a graduated cylinder.

  • Aerial Part: the aerial part height (APH, cm) was measured as the distance from the soil surface to the top of the shoot, measured just before harvest; the aerial part diameter (APD, cm) was obtained by taking two measurements of the diameter of each plant with a ruler, one in the horizontal plane and one in the perpendicular direction; the aerial part volume (APV, cm3) was measured by the displacement volume of water inside a graduated cylinder of known dimensions; the aerial part total mass (APTM, g) was expressed as the sum of the commercial and noncommercial masses per plant; the aerial part commercial mass (APCM, g) was measured as the mass of the marketable aerial part; the mass of the aerial part loss at harvest (MLH, g) per plant was determined by subtracting the commercial mass from the total aerial part mass; the crop yield (Y, kg m-2) was an estimation of the ratio between the production of the useful plot (based on APCM) and the area occupied by it, extrapolating to an area of 1 m2; the water productivity (WP, kg m-3) (water use efficiency by area) was expressed as the ratio between Y (kg m-2) and the amount of water applied (mm) per treatment.

Statistical analysis

The obtained data were subjected to analysis of variance at 5% and 1% probability by the F test. When significant at 5%, the effects of depth were subjected to regression analysis, and the effects of cultivars were compared by the Tukey test at 5% probability. For regression analysis, the linear, quadratic, cubic, square root, potential, exponential, hyperbolic, logarithmic, cubic root, log-log, Ln-Ln and Exp(x) models were tested. The best model was obtained considering the significance of the F test for the regression equation, up to a 5% probability, the coefficient of determination (R2), and the representativeness of biological behavior by the equation. The curve inflection points were obtained by the first-order partial derivative of the respective regression equations. Data were subjected to principal component analysis (PCA). We used GENES and SISVAR software for statistical analysis.

Results and discussion

The coefficients of variation (CV%) from the analysis of variance (Tables 3, 4 and 5) for the three crop cycles ranged from 5.4% to 45.0% in the plot (average of 25%) and from 4.7% to 5.0% in the subplot (average of 27%). The highest CV values occurred for MLH, APV and LAI, which are directly linked to fresh phytomass, especially of the leaves. Lettuce has, on average, 94% water in its constitution, and therefore, part of the inherent variability of this culture occurs due to variations of the water content from one plant to another, increasing the CV of the experiment [30].

Tab. 3. Average squares, significance of the F test (ANOVA) and average values for leaf characteristics: relative content of chlorophylls A (Chl-A) and B (Chl-B) and of total chlorophyll (Chl-T), number of commercial leaves (NCL), number of noncommercial leaves (NNCL), leaf length and width (LL and LW), leaf area index (LAI), fresh leaf mass (FLM), and dry leaf mass (DLM g) as a function of different cultivars and irrigation depths in three lettuce crop cycles.
Average squares, significance of the F test (ANOVA) and average values for leaf characteristics: relative content of chlorophylls A (Chl-A) and B (Chl-B) and of total chlorophyll (Chl-T), number of commercial leaves (NCL), number of noncommercial leaves (NNCL), leaf length and width (LL and LW), leaf area index (LAI), fresh leaf mass (FLM), and dry leaf mass (DLM g) as a function of different cultivars and irrigation depths in three lettuce crop cycles.
Tab. 4. Average square significance of the F test (ANOVA) and average values for characteristics of roots and stems: fresh root mass (FRM), dry root mass (DRM), root length (RL), root volume (RV), fresh stem mass (FSM), stem length (SL), stem average diameter (SMD) and stem volume (SV) as a function of different cultivars and irrigation depths in three lettuce crop cycles.
Average square significance of the F test (ANOVA) and average values for characteristics of roots and stems: fresh root mass (FRM), dry root mass (DRM), root length (RL), root volume (RV), fresh stem mass (FSM), stem length (SL), stem average diameter (SMD) and stem volume (SV) as a function of different cultivars and irrigation depths in three lettuce crop cycles.
Tab. 5. Average square significance of the F test (ANOVA) and average values for aerial part characteristics: aerial part height (APH), diameter (APD) and volume (APV), aerial part total mass (APTM) and commercial mass (APCM), mass of aerial part loss at harvest (MLH), crop yield (Y) and water productivity (WP) as a function of different cultivars and irrigation depths in three lettuce crop cycles.
Average square significance of the F test (ANOVA) and average values for aerial part characteristics: aerial part height (APH), diameter (APD) and volume (APV), aerial part total mass (APTM) and commercial mass (APCM), mass of aerial part loss at harvest (MLH), crop yield (Y) and water productivity (WP) as a function of different cultivars and irrigation depths in three lettuce crop cycles.

Tables 3, 4 and 5 show that there were interactions between irrigation depth and lettuce cultivars only for cycle 2 (APCm2, Y2, WP2, WP3, APD2, NCL2, FLm2, LAI2, FSm2, SMD2, and SV2) and for cycle 3 (WP3). A significant effect for lettuce cultivars, as a variation factor, was observed in almost all evaluated characteristics for the three crop cycles. Irrigation depth had a significant effect on APTMT2, WP1, APH1, APH2, DLm2 and SL2.

The existence of a linear relationship between the number of leaves and the air temperature, accumulated through the degrees per day (GD) thermal summation, has confirmed that temperature is the main factor controlling the rate of leaf emission [27]. However, in the present study, the highest number of leaves per plant for all cultivars tested (Table 3) occurred in cycle 2, in which the average temperature was lower than that of cycle 1 (Fig 2). The lower incidence of radiation in cycle 2 (Fig 2) induced an increase in the number of leaves in the plant, increasing the amount of chlorophyll and consequent maintenance of the photosynthetic rate required for the plant.

The Raider Plus cultivar better adapt to environments with a lower incidence of light, similar to the protected environment used in this study, since it presented a higher percentage of chlorophylls A and B and total chlorophyll (Table 3) than that of the other cultivars and, consequently, carried out more photosynthesis. A higher concentration of chlorophyll increased the generation of leaf quantity (higher NCL and NNCL), leading to relatively high phytomass values (FLM and DLM). Leaf quality (higher LL and LW) was also influenced when the plants had higher leaf areas (LL and LW) and, consequently, a higher specific leaf area [31] increased the LAI of the Raider Plus cultivar, which also propagated greater phytomass of the plants [32].

The quality and freshness of commercially sold vegetables can be monitored by measuring the chlorophyll content [28]. According to [33], changes in total chlorophyll concentration relative to the photosynthetic rate may follow an opposite-to-expected pattern, exhibiting higher values under reduced radiation intensities. This increase in pigment content in shaded leaves has been attributed to the increase in the number and size of the chloroplasts, the amount of chlorophyll produced by chloroplasts and/or better development of the Grana within the chloroplasts. This effect is considered an efficient adaptive process of the acclimatization of plants to environments with low light intensity [34].

This behavior benefited the development of the Raider Plus cultivar, intensifying the green coloration and increasing its visual quality, which promotes greater attractiveness and commercialization in relation to the other cultivars.

The Luara cultivar also showed high values of chlorophyll (Table 3), but its NCL values were the most prominent in all crop cycles, independent of the applied irrigation depth, which reflected in larger aerial parts (APD and APV) (Table 5). According to [35], a higher NCL increases the quality of lettuce, adding value to the crop during its commercialization, and this characteristic is often linked to genetic factors.

Better photosynthetic characteristics, for both the Raider Plus and Luara cultivars, promoted better root contribution (greater FRM, DRM, and RL), with a higher volume of explored soil (RV), reducing water stress of the plant even under smaller applied irrigation depths (50% and 75% of ETc) (Table 4). The optimal root characteristics possibly generated higher nutrient absorption and better soil exploration, increasing the amount and quality of leaves in the plant and therefore facilitating photosynthesis [36]. The increase in leaf quantity and quality provided higher LAI values (Table 3), with fast closure of the plant stand and less competition with unwanted plants, with further development of the Luara cultivar [37].

Although the Vanda cultivar had higher values of some of the characteristics, such as LAI, LL, FRM, DRM, RL, RV and APD (Tables 3, 4 and 5), higher values of stem characteristics [FSM, SMD, SV, and mainly SL (plant weeding)] and lower values of leaf characteristics reduce the visual and commercial quality of this cultivar compared to the others. The Imperial Roxa cultivar was less adapted to the conditions studied since it presented lower values in almost all characteristics, except for SL, an undesirable characteristic in lettuce commercialization. Different species, under water stress, may present survival mechanisms, including variations in the stomatal responses, osmotic adjustment and greater movement of photoassimilates to the roots [38], with metabolic energy consumption, which may affect the plant’s production and commercial quality [9].

As described previously, the plant physiological factors of the Raider Plus and Luara cultivars were better than those of the Vanda and Imperial Roxa cultivars. The plant physiological factors were influenced in a directly proportional way in production (APTM and APCM), as well as in Y. The Raider Plus and Luara cultivars, with larger numbers of leaves (NCL and NNCL), smaller stems (lower SL, SMD and SV) and higher volumes (higher APV), presented higher APTM and APCM values, which are directly proportional to Y.

Although the Raider Plus cultivar presented high crop loss indexes (MLH), APTM was enough to confer higher APCM and higher Y for this crop. The highest yield average of the three cycles was observed for the Raider Plus cultivar (3.0 kg m-2), followed by the Luara (2.8 kg m-2), Vanda (2.7 kg m-2) and the Imperial Roxa (1.5 kg m-2) cultivars.

The highest values of Y observed in the present study for the Raider Plus, Luara and Vanda cultivars were similar to the maximum value (3.6 kg m-2) presented by [39] for conventional lettuce cultivation in soil.

Fig 3 presents the average behavior of the agronomic characteristics of the four cultivars that did not show significant interaction at a 5% probability by the F test between the cultivars studied and the irrigation depths, with variable maximum values and an irrigation depth corresponding to the curve inflexion point (IDIP).

Average values of the aerial part height (APH) and water productivity (WP) of cycle 1, and the aerial part total mass (APTM), aerial part height (APH), dry leaf mass (DLM) and stem length (SL) of cycle 2 as a function of irrigation depths (ID) for lettuce cultivars. The number after the variable abbreviation indicates the cultivation cycle.
Fig. 3. Average values of the aerial part height (APH) and water productivity (WP) of cycle 1, and the aerial part total mass (APTM), aerial part height (APH), dry leaf mass (DLM) and stem length (SL) of cycle 2 as a function of irrigation depths (ID) for lettuce cultivars. The number after the variable abbreviation indicates the cultivation cycle.
IDIP: irrigation depth referring to the inflection point of the curve.

An approximate average value of water replacement depth equal to 110% ETc provides water conditions for the lettuce crop to express its greatest potential in important variables (APTM and DLM), with PA near to the average (Fig 3).

Interactions between the studied cultivars, as a function of the irrigation depth and referring to the maximum point of the curve inflection point (IDIP) obtained by deriving the respective equations, are shown in Fig 4. Note that trend lines were used to highlight only the cultivars that were statistically influenced (5% probability by the F test) by the variation depths, but the averages of the other cultivars were also inserted (Fig 4).

Average values of the number of commercial leaves (NCL), aerial part diameter (APD), fresh leaf mass (FLM), leaf area index (LAI), fresh stem mass (FSM), stem average diameter (SMD), stem volume (SV), aerial part total mass (APTM), aerial part commercial mass (APCM), and crop yield (Y) of cycle 2 and the water productivity (WP) of cycles 2 and 3 as a function of irrigation depths (ID) for lettuce cultivars.
Fig. 4. Average values of the number of commercial leaves (NCL), aerial part diameter (APD), fresh leaf mass (FLM), leaf area index (LAI), fresh stem mass (FSM), stem average diameter (SMD), stem volume (SV), aerial part total mass (APTM), aerial part commercial mass (APCM), and crop yield (Y) of cycle 2 and the water productivity (WP) of cycles 2 and 3 as a function of irrigation depths (ID) for lettuce cultivars.
The number after the variable abbreviation indicates the cultivation cycle. IDIP: irrigation depth referring to the inflection point of the curve.

The Luara cultivar was the most influenced by irrigation depth (Fig 4), explained by second-degree polynomial equations for the characteristics NCL2, FLm2, LAI2, APTm2, APCm2, Y2, FSm2, SMD2, and SV2 and by linear equations for WP in cycles 2 and 3. The Vanda cultivar was the second most influenced by irrigation depth, adjusted by linear equations for FSm2, SMD2, SV2, and WP in cycles 2 and 3. The Raider Plus cultivar was slightly affected by irrigation depth, as explained by the second-degree polynomial equation for APD2, LAI2 and WP in the last two cycles. Finally, the Imperial Roxa cultivar showed an irrigation depth influence only for WP in the third cycle.

Although higher water productivity (WP) was observed for the lowest irrigation depth used (50% of ETc), similar to [9], important variables, such as NCL and FLM, reached maximum values by applying an average depth of approximately 110% of ETc (Fig 4).

Note that the optimal aerial part characteristics, which directly influence Y, were obtained by irrigation depths higher than 100% water replenishment. Possibly, the global efficiency of the system did not achieve 100% efficiency in water absorption because of losses by percolation, redistribution of water in the soil and areas with a water deficit [40]. Another explanation for this fact is the potential soil salinity caused by the high salt levels (Table 1) accumulated in the soil used in the present research.

The principal component analysis (PCA) is represented in Fig 5. The first major component (PC1) explains a large portion of the parameter variances (agronomic characteristics). The second main component (PC2) explains a smaller part of said variance but more than that of the last component (coincident with the number of parameters evaluated), which explains little or almost nothing of the existing variance. Usually, only the first two components are graphically represented (Fig 5) because together, they generally represent more than 80% of this variance [41, 42, 43].

Principal component analysis (PCA) for the aerial part total mass (APTM), aerial part commercial mass (APCM); mass of aerial part loss upon harvesting (MLH), crop yield (Y), water productivity (WP), aerial part height (APH), aerial part diameter (APD), aerial part volume (APV), leaf length (LL), leaf width (LW), fresh leaf mass (FLM), dry leaf mass (DLM), leaf area index (LAI), fresh stem mass (FSM), stem length (SL), stem average diameter (SMD), stem volume (SV) and content of total chlorophyll (Chl-T) in cycles 1, 2 and 3 under different irrigation depths for lettuce cultivars.
Fig. 5. Principal component analysis (PCA) for the aerial part total mass (APTM), aerial part commercial mass (APCM); mass of aerial part loss upon harvesting (MLH), crop yield (Y), water productivity (WP), aerial part height (APH), aerial part diameter (APD), aerial part volume (APV), leaf length (LL), leaf width (LW), fresh leaf mass (FLM), dry leaf mass (DLM), leaf area index (LAI), fresh stem mass (FSM), stem length (SL), stem average diameter (SMD), stem volume (SV) and content of total chlorophyll (Chl-T) in cycles 1, 2 and 3 under different irrigation depths for lettuce cultivars.
Trat: Treatment. ID: irrigation depth. Circles: green (Raider Plus); blue (Luara); red (Vanda); and purple (Imperial Purple). Viçosa-MG, DEA-UFV, 2017.

Comparatively, the closer that a treatment is plotted to a given characteristic, the more suitable this treatment will be to obtain higher values of that characteristic.

Treatments involving the Raider Plus cultivar (1, 5, 9, 13 and 17), regardless of the irrigation depth, were concentrated on the right of the central axis (Fig 5). This behavior of Raider Plus occurred mainly in the first and third cycles, with intense association with variables that positively influence the quality and productivity of the lettuce, such as Chl-T, LAI and Y (first cycle), Chl-T, LAI and WP (second cycle), and APV, Y and WP (third cycle). Treatments involving the Luara cultivar (Fig 5) favored the NCL and APD variables in the first cycle (treatments 2, 10 and 14), NCL, APV and Y in the second cycle (treatments 8 and 11), and APD, APH and NCL in the third cycle (treatments 2, 6, 10, 14 and 18).

Unlike the Raider Plus and Luara cultivars, the Vanda and Imperial Roxa cultivars concentrated their treatments on the negative side of the CP1 axis (Fig 5), in the opposite direction of the important variables for lettuce quality, confirming their nonpreference for the important variables. The situation of the Imperial Roxa cultivar was even more removed, with reduced association with the variables of interest for lettuce production.

The angle between the direction of any two vectors (Fig 5) indicates the association amplitude between the variables represented by these vectors, and it is inversely proportional with the correlation; that is, the smaller the angle between vectors referring to the characteristics, the greater the correlation degree between them. Thus, Y, as the main characteristic, was positively influenced by the APV, APCM, APTM, LAI, and Chl-T characteristics in the first cycle, by APV, APCM and APTM in the second cycle, and by LAI, APD, APH, NCL, APTM and APCM in the third cycle. In contrast, Y was negatively influenced by APH and APD in cycle 1 and by Chl-T in cycles 2 and 3.

Conclusions

The Raider Plus and Luara cultivars can be recommended for the cultivation of lettuce in protected environments, with blue plastic covering, in the Zona da Mata of Minas Gerais State and in regions with similar conditions.

Lettuce irrigation should be performed with a depth to restore 110% of the crop evapotranspiration to the region under conditions similar to those of the current study.

Future research comparing the same cultivars and irrigation depths in protected environments with different coverings is encouraged by the authors and may provide different conclusions than that of the current research.

Supporting information

S1 Fig [tif]
Research localization (Viçosa-MG, Brazil).

S1 Table [tif]
Practical example of opening sideline valves.


Zdroje

1. Figueiredo CC, Ramos MLG, Mcmanus CM, Menezes AM. (2012). Mineralização de esterco de ovinos e sua influência na produção de alface. Horticultura Brasileira 30(1): 175–179. Portuguese.

2. Koefender J, Schoffel A, Manfio CE, Golle DP, Silva NA, Horn RC. (2016). Consorciação entre alface e cebola em diferentes espaçamentos. Horticultura Brasileira, 34: 580–583. Portuguese.

3. Sala FC, Costa CP. (2012). Retrospectiva e tendência da alfacicultura brasileira. Horticultura Brasileira 30(2): 187–194. Portuguese.

4. Caliskan S, Yetisir H, Karanlik S. (2014). Combined use of green manure and farmyard manure allows better nutrition of organic lettuce. Notulae Botanicae Horti Agrobotanici Cluj-Napoca 42(1): 248–254.

5. Kapoulas N, Koukounaras A, Ilić ZS. (2017). Nutritional quality of lettuce and onion as companion plants from organic and conventional production in north Greece. Scientia Horticulturae 219: 310–318.

6. Nandi R, Gowdru N, Bokelmann W, Dias G. (2015). Smallholder organic farmer’s attitudes, objectives and barriers towards production of organic fruits and vegetables in India: a multivariate analysis. Emirates Journal of Food and Agriculture 27(5): 396–406.

7. Arantes CRDA Junior SS, Camili EC Diamante MS, Pinto ESC. (2014). Produção e tolerância ao pendoamento de alface-romana em diferentes ambientes. Revista Ceres 61(5): 558–566. Portuguese.

8. Foteinis S, Chatzisymeon E. (2016). Life cycle assessment of organic versus conventional agriculture. A case study of lettuce cultivation in Greece. Journal of Cleaner Production 112(4): 2462–2471.

9. Bandeira GRL, Pinto HCS, Magalhães OS, Aragão CA, Queiroz SOP, Souza ER, et al. (2011). Manejo de irrigação para cultivo de alface em ambiente protegido. Horticultura Brasileira 29: 237–241. Portuguese.

10. Tringovska I, Yankova V, Markova D, Mihov M. (2015). Effect of companion plants on tomato greenhouse production. Scientia Horticulturae 186(21): 31–37.

11. Carvalho J de A, Rezende FC, Aquino RF, Freitas WA de; Oliveira EC. (2012). Produção da ervilha cultivada em ambiente protegido sob diferentes tensões de água no solo. Revista Brasileira de Engenharia Agrícola e Ambiental 16(1): 44–50. Portuguese.

12. Silva JA, Pires RCM, Sakai E, Silva TJA, Andrade JE, Arruda FB, et al. (2007). Desenvolvimento e produtividade da cultura da batata irrigada por gotejamento em dois sistemas de cultivo. Engenharia Agrícola 27(2): 354–362. Portuguese.

13. Reichardt K. (1978). A água na produção agrícola. São Paulo: McGraw-Hill do Brasil. Portuguese.

14. Alvares CA, Stape JL, Sentelhas PC, Goncalves JLM, Sparovek G. (2013). Köppen's climate classification map for Brazil. Meteorologische Zeitschrift 22(6): 711–728.

15. Ribeiro AC, Guimarães PTG, Alvarez VVH. (1999). Recomendações para o uso de corretivos e fertilizantes em Minas Gerais– 5ª aproximação. Viçosa: CFSEMG, p.177. Portuguese.

16. Embrapa—Empresa Brasileira de Pesquisa Agropecuária. (2011). Manual de métodos de análise de solo. 2. ed. Rio de Janeiro: Embrapa Solos, 230 p. Portuguese.

17. Richards LA. (1949). Methods of measuring soil moisture tension. Soil Science of American Journal 68(1): 95–112.

18. Embrapa—Empresa Brasileira de Pesquisa Agropecuária. (2013). Centro Nacional de Pesquisa de Solos. Sistema brasileiro de classificação de solos. 3.ed. Rio de Janeiro, Embrapa Solos, 353 p. Portuguese.

19. Keller J, Karmeli D. (1975). Trickle irrigation design. Glendora: Rain Bird Sprinkler Manufacturing, 133 p.

20. Denículi W, Bernardo S, Thiébaut JTL, Sediyama GC. (1980). Uniformidade distribuição de água, em condições de campo num sistema de irrigação por gotejamento. Revista Ceres 27(150): 155–162. Portuguese.

21. Bernardo S, Soares AA, Mantovani EC. (2006). Manual de Irrigação. 8. ed. Viçosa: Editora UFV, 625 p. Portuguese.

22. Allen RG, Pereira LS, Raes D, Smith M. (1998). Crop evapotranspiration: guidelines for computing crop water requirements. Rome: FAO, 297 p. (Irrigation and Drainage Paper, 56).

23. Doorenbos J, Pruit WO. (1977). Guidelines for predicting crop water requirements. Rome: FAO, 179p. (FAO: Irrigation and Drainage Paper, 24).

24. Marouelli WA, Silva WLC, Silva HR. (2008). Irrigação por aspersão em hortaliças: qualidade da água, aspectos do sistema e método prático de manejo. 2 ed. Brasília: Embrapa Informação Tecnológica, 150p. Portuguese.

25. Keller J. (1978). Trickle irrigation. Colorado: En Soil Conservation Service National Engineering Handbook, 129 p.

26. Brunini GA, Lisbão RS, Bernardi JB. (1976). Temperatura base para a Alface (Lactuca sativa L.) cultivar “White Boston”, em sistema de unidades térmicas. Revista de Olericultura 35(19): 214–219. Portuguese.

27. Araújo TS, Fideles Filho J, Kumar KK, Rao TVR. (2010). Crescimento da alface-americana em função dos ambientes, épocas e graus-dias. Revista Brasileira de Ciências Agrárias 5(4): 441–449. Portuguese.

28. Limantara L, Dettling M, Indrawati R, Brotosudarmo THP. (2015). Analysis on the chlorophyll content of commercial green leafy vegetables. Procedia Chemistry, 14: 225.

29. Lucena RRM, Batista TMV, Dombroski JLD, Lopes WAR, Rodrigues GSO. (2011). Medição de área foliar de aceloreira. Revista Caatinga 24: 40–45. Portuguese.

30. Lúcio AD, Haesbaert FM, Santos D, Benz V. (2011). Estimativa do tamanho de parcela para experimentos com alface. Horticultura Brasileira 29: 510–515. Portuguese.

31. Puiatti M, Finger FL. (2005). Fatores climáticos. In: FONTES, P.C.R. Olericultura: teoria e prática. Viçosa: Ed. da UFV, p.17–30. Portuguese.

32. Radin B, Reisser Júnior C, Matzenauer R, Bergamashi H. (2004). Crescimento de cultivares de alface conduzidas em estufa e a campo. Horticultura Brasileira 22: 178–181. Portuguese.

33. Vandana V, Bhatt RK. (1999). Physiological changes in Sesbania species to reducing light intensities. Journal of Agronomy and Crop Science 182:43–47.

34. Kosma C, Triantafyllidis V, Papasavvas A, Salahas G, Patakas A. (2013). Yield and nutritional quality of greenhouse lettuce as affected by shading and cultivation season. Emirates Journal of Food and Agriculture 25(12): 974–979.

35. Jordan RA, Geisenhoff LO1, Oliveira FC de, Santos RC, Martins EAS. (2018). Yield of lettuce grown in aquaponic system using different substrates. Revista Brasileira de Engenharia Agrícola e Ambiental 22(1): 27–31.

36. Taiz L, Zeiger E, Moller IM, Murphy A. (2017). Fisiologia e desenvolvimento vegetal. 6. ed.–Porto Alegre: Artmed. Portuguese.

37. Oliveira F de A de, Carrilho MJS de O, Medeiros JF de, Maracajá PB, Oliveira MKT de. (2011). Desempenho de cultivares de alface submetidas a diferentes níveis de salinidade da água de irrigação. Revista Brasileira de Engenharia Agrícola e Ambiental 15(8): 771–777. Portuguese.

38. Jaimez RE, Rada F, Garcia-Nuñez C. (1999). The effect of irrigation frequency on water and carbon relations in three cultivars of sweet pepper (Capsicum chinense Jacq), in a tropical semiarid region. Scientia Horticulturae 81: 301–305.

39. Peixoto Filho JU, Freire MBG dos S, Freire FJ, Miranda MFA, Pessoa LGM, Kamimura KM. (2013). Produtividade de alface com doses de esterco de frango, bovino e ovino em cultivos sucessivos. Revista Brasileira de Engenharia Agrícola e Ambiental 17: 419–424. Portuguese.

40. Vilas Boas RC, Carvalho JA, Gomes LAA, Sousa AMG, Rodrigues RC, Souza KJ. (2008). Avaliação técnica e econômica da produção de duas cultivatires de alface tipo crespa em função de lâminas de irrigação. Ciência e Agrotecnologia 32(2): 525–531. Portuguese.

41. Pearson K. (1901). On lines and planes of closest fit to systems of points in space. Philosophical Magazine 2(6): 559–572.

42. Johnson RA, Wichern DW. (1999). Applied multivariate statistical analysis. 4th ed. Upper Saddle River, New Jersey: Prentice-Hall, 815 p.23.

43. Gotelli NJ, Ellison AM. (2010). Princípios de Estatística em Ecologia. Artmed; Edição 1. Portuguese.


Článek vyšel v časopise

PLOS One


2019 Číslo 12
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#