Skip to main content

Growth response of greenhouse-produced muskmelon and tomato to sub-surface drip irrigation and soil aeration management factors



Hypoxia causes injury and yield loss. Soil aeration has been reported to accelerate the growth of plants and increase crop yield. The aim of this study was to examine growth response of greenhouse-produced muskmelon to 3 levels of sub-surface drip irrigation (I), 3 different installation depths of drip laterals in the soil (D), and 4 levels of supplemental soil aeration frequency (A). A fractional factorial experiment was designed to examine these treatment effects on marketable fresh fruit yield, leaf area index during 3 growth stages, and dry matter partitioning at harvest. In addition, we studied the response of fruit yield and dry matter of tomato to 2 levels of burial depths of subsurface tubing in combination with 3 frequency levels of soil aeration.


Results showed that soil aeration can positively influence the yield, leaf area index, dry matter and irrigation use efficiency of the muskmelon (p < 0.05). The fruit yield of muskmelon and tomato were increased by 21.5 and 30.8% respectively with 1-d and 2-d aeration intervals compared with the no aeration treatment.


The results suggest that soil aeration can positively impact the plant root zone environment and more benefits can be obtained with aeration for both muskmelon and tomato plants.


Use of sub-surface drip irrigation under protective structures to produce high-value vegetables for the fresh food market is increasing in the semi-arid areas of Northwest China. These production technologies significantly improve water use efficiency compared to traditional spray or furrow irrigation systems [1,2,3]. The topsoil water content depends on the sub-surface drip application frequency, rate, depths, leakage and evapotranspiration. Manufacturer can vary the sub-surface depth of the buried drip tubing to facilitate cultivation along with the rate, amount, and frequency of irrigation to control the wetted portion of the root zone and optimize soil water availability and accessibility to the plant root system [4,5,6].

In addition to soil water, rhizosphere O2 content is one of the most important factor affecting respire, grow, develop, and function normally of plant roots [5, 7]. Most of O2 supply is obtained directly by continuous diffusive air exchange between the atmosphere and soil [8]. Previous studies have shown that hypoxia stress limited leaf growth, the rate of leaf expansion declined, and the size of leaf reduced [9]. High CO2 concentrations (at 2500 μL L− 1) or low (10% by volume) O2 in the root-zone significantly increased the activity of lactate dehydrogenase (LDH), alcohol dehydrogenase (ADH), and pyruvate decarboxylase (PDC) in muskmelon plants, resulting in inhibition of growth, and decrease in fruit yield and quality [10,11,12].

Subsurface drip irrigation allows the application of air and soluble materials (such as O2, nutrients, pesticide) directly to the root zone. Studies have shown that these practices can be very useful in overcoming problems associated with hypoxia in the root zone of irrigated crops over a range of crops, soil water contents and soil types, and for improving crop performance under oxygen-deficient conditions [13,14,15]. Aerated irrigation positively affected yield and improved quality for cotton, beans, squash, pumpkin, cucumber, tomato, rice and other crops especially on heavy and saline soils [16,17,18,19,20].

These demonstrated positive effects of aerated irrigation and supplemental aeration on growth and yield have not been completely explained in the literature. It is quite possible that improved LAI of the leaves may be a key mechanism for the increased yield and quality. For given ambient conditions (light, temperature, water, nutrition) crop growth and yield is directly related to LAI. Higher LAI, photosynthetic rate and chlorophyll content are key indicators and predictors of higher crop yields [21]. However, it is not known how changing sub-surface irrigation practices (tubing placement depth, irrigation level and frequency) and artificial aeration frequency would affect LAI.

We suspect that varying the aeration frequency, aeration volume, aeration position and irrigation amount would result in a possible change in rhizosphere soil environment, LAI, and irrigation use efficiency of plants, thus changing the plant growth and fruit output. The specific objective of this study was to determine the response of greenhouse-produced muskmelon and tomato plants on clay loam soils of the semi-arid area of Northwest China under various combinations of sub-surface irrigation volumes, tubing placement depth, and supplemental aeration frequency. Fresh fruit yield, leaf area index, and dry matter partitioning were taken as measures of response to these treatment factors.


Muskmelon results

LAI of muskmelon

Values in the LAI25 row of Table 1 corresponded to the tendril elongation growth period. The 55 DAT (LAI55) and 75 DAT (LAI75) measurements were taken during the flowering and at fruit harvest. The overall mean of 2.69 at 40 DAT was close to double the value of 1.31 at 25 DAT. At 75 DAT the LAI of 2.86 was only about 4% higher than at 40 DAT. The 10 cm tubing placement depth tended to be significantly lower compared to 25 and 40 cm depth for the 55 and 75 DAT measurements. Differences between the LAI marginal means for irrigation level were non-significant for all 3 sampling dates (Table 1). Among the single-factor analyses, tubing placement depth had a significant impact on both LAI55 and LAI75, aeration frequency had an extremely significant impact on LAI25. The identified significant D x A interaction for LAI25, and A x I interaction for LAI during the tendril elongation growth period and fruit harvest period.

Table 1 Marginal means and 2-factor interactions for marketable fruit yield, LAI, post-harvest dry matter partitioning, and irrigation use efficiency response of greenhouse-produced muskmelon plants to sub-surface drip-irrigation tubing placement depth x soil aeration frequency x irrigation level treatment factors. In each row, the means for the levels of a given factor not followed by the same letter are significantly different at the 5% level. Two-factor interaction is significant at the 1% (**), 5% (*), or not significant (ns) level

The 1-way ANOVA showed that LAI25 values were all < 2.00 for all treatment combinations (Fig. 1). The highest LAI25 value of 1.87 for the D10A1I80 treatment combination was significantly higher (P<0.01) than all the non-aerated combinations with values < 1.00 (namely D10ANI70 = 0.81, D25ANI80 = 0.85, D40ANI90 = 0.85) regardless of depth and irrigation level. Treatment combinations with daily soil aeration tended to be higher or similar to those with the 2-day and 4-day soil aeration frequency. Combination with D40 and A4 gave LAI25 < 1. Combinations with higher A2 and A1 aeration frequency and shallower D25 and D10 placement depth gave LAI25 > 1 indicating that lowering placement depth could offset the negative impact of lowering soil aeration frequency (Fig. 1).

Fig. 1

LAI of muskmelon plants measured by canopy solar radiation transmittance at 25, 55, and 75 days after transplanting for each combination of drip-tubing placement depth, soil aeration frequency, and irrigation level treatment factors. Combination treatment means are annotated along polar axis. Ducan 0.05 for mean comparison at 25, 55, and 75 DAT were 0.50, 1.16, and 1.12

The highest LAI55 value of 3.60 during the flowering period was obtained for the D25A1I90. LAI55 values > 3.00 were obtained only for the D25A4I80 (= 3.22) and D40A1I70 (= 3.32) treatment combinations. These LAI75 values for these combinations were practically the same as those for LAI55. It would appear that 3.00 to 3.60 may represent the maximum LAI range for the Shantian No. 1 muskmelon cultivar as measured under the experimental conditions.

Except for the D10ANI70 treatment combination with LAI55 = 1.87, all the 8 remaining combination LAI55 values were between 2.00 and 3.00. For the D25A1I90 combination, the LAI55 = 3.60 represented a 116% increase on the LAI25 = 1.66 value. Corresponding increases were 112 and 154% for the D25A4I80 and D40A1I70 treatment combinations. Increases for LAI55 relative to LAI25 tended to be higher for combinations with the lower (A4 and AN) soil aeration frequency and especially for the higher irrigation levels (I90) (Table 1). As indicated the LAI75 values were very close to the LAI55 for some combinations indicating little leaf growth after the flowering and during the fruit enlargement growth stages. LAI75 values were markedly lower than LAI55 for the D10ANI70 and D40ANI90 treatment combinations indicating possible earlier senescence.

Dry matter partitioning and irrigation use efficiency of muskmelon

Total dry matter of muskmelon progressively increased with increasing tubing placement depth and progressively decreased with decreasing soil aeration frequency (Table 1). Decreased soil aeration frequency tended to progressively reduce all 3 components of the total DM. The effects were more marked for the leaf and root contributions. The marginal mean for leaf and root DM with daily soil aeration was 50 and 71% higher than the values for no aeration. The corresponding value for stem DM was 32%. The marginal mean for root DM was significantly higher for the I70 irrigation level compared to the I80 and I90 levels. In addition, interaction analysis found that D x A and D x I significantly affected the root DM, D x I and A x I significantly affected the Total DM and Stem DM. Overall mean total DM was 52.2 g per plant partitioned as 11.9 g in stem+ 38.3 g in leaf + 2.0 g in root representing 22.7, 73.4, and 3.9% respectively. These DM partitioning percentages were similar for each of the tubing placement depths, soil aeration frequencies, and irrigation levels. Values ranged between 21 and 26%, 68–75%, and 3–4% respectively for the stem, leaf, and root contribution to the total DM for all levels of the treatment factors (Table 1). The D40A1I70 treatment combination resulted in the maximum plant root, stem, leaf and total DM values (Fig. 2). As shown, the plant root, stem, leaf and total DM values for treatment combinations with daily soil aeration tended to be higher than those with the no aeration.

Fig. 2

Post-harvest stem, leaf, and root dry matter of muskmelon plants measured at 75 days after transplanting for each combination of drip-tubing placement depth, soil aeration frequency, and irrigation level treatment factors. Bars on the stacked sections for each component represent the standard error of the mean. Means not followed by the same letter are significantly different at the 5% level

The ternary plot (Fig. 3) show the range of the ratios for root, stem, and leaf DM of muskmelon to total post-harvest DM at 75 days after transplanting for each treatment combination. The root DM percentages fell between 2.8 to 5.4%. The leaf DM range was 67.5 to 75.7% and 20.6 to 28.5% for the stem DM. These ranges were quite narrow indicating that the treatment combinations did not markedly affect DM partitioning in the muskmelon plants. Nevertheless, it would be interesting to examine the clusters of treatment combinations at the high and low ends of these ranges as shown in Fig. 3. The cluster with high stem DM percentages was D10ANI70 > D10A4I80 > D10A1I80 treatment combinations (Fig. 3). The treatment cluster of low leaf DM percentages was D10ANI70 < D10A4I80 < D10A1I80 compared with the cluster of high leaf DM percentages D25ANI80 > D25A1I90 > D40A1I70 treatment combinations. High root DM percentages was obtained for the D10A4I80 > D10A1I80 cluster of treatment combinations. These clusters tended to confirm the existence of 2-factor interactions (Table 1).

Fig. 3

Ternary plot of root, stem, and leaf dry weight of muskmelon expressed as % of total post-harvest plant dry weight measured at 75 days after transplanting foreach combination of drip tubing placement depth, soil aeration frequency, and irrigation level treatment factors

The overall irrigation use efficiency for marketable fruit yield of muskmelon (IUEf) across each treatment factor in Table 1 was 17.2 g per plant per liter of total irrigation water applied (subsurface drip + initial 70 mm surface irrigation immediately at prior to transplanting). Neither tubing placement depth nor soil aeration resulted in significant differences in the IUEf marginal means. The IUEf marginal mean was 21.5 g fresh fruit per plant per liter for the I70 irrigation level. Corresponding values decreased linearly by 20 and 40% for I80 (=17.2) and I90 (=12.9). These decreases were significantly (p < 0.05) different. The 3-way ANOVA indicated a highly significant 2-factor A x I interaction effect (Table 1).

The overall irrigation use efficiency (IUEdm) for total dry matter (stem + leaf + root) of muskmelon across each treatment factor in Table 1 was 7.58 g per plant per liter of total irrigation water applied. All treatment factors significantly affected IUEdm and all 2-factor interactions were highly significant (Table 1). The IUEdm value for D40 (=8.94) was significantly higher than for D10 (=6.15) by 45%. Similarly, IUEdm progressively increased with increasing aeration frequency from 6.20 for AN to 6.53, for A4, to 7.81 for A2, and 9.75 for A1. The 57% increase between the values for AN and A1, and the 49% increase between A4 and A1, were significant (P < 0.05). IUEdm values decreased with increasing irrigation level. The marginal mean of 9.23 for I70 differed significantly from to the value of 7.29 for I80 (20% lower), and 6.20 for I90 (33% decrease).

Marketable fresh fruit yield of muskmelon

The results in Table 1 show significant treatment main effects for the drip tubing placement depth and soil aeration frequency experimental factors. Increasing depth of tubing tended to decrease the marketable fresh fruit yield per plant in a non-linear manner. Values increased non-significantly from 1.48 to 1.59 kg per plant between the 10 cm to 25 cm tubing placement depth. However, the mean of 1.38 kg per plant for the 40 cm placement depth was significantly lower than for the 25 cm depth (Table 1). Only the D x I two-factor interaction on fresh fruit yield was significant (Table 1) indicating that the tubing placement depth effect would be different for the various irrigation levels. This was due to tubing placement depth determine the position of the soil wetted. Uneven redistribution of soil water in different longitudinal position result in different water use efficiency and fruit yield. It seems that the optimal installation depth of drip pipes is around 25 cm. However, the fractional factorial design does not permit examination of the simple effects of placement depth at each irrigation level.

Yield progressively decreased with decreasing soil aeration frequency from 1.64 kg/plant with daily aeration to 1.51, 1.43 1.35 with 2-day, 4- day, and no aeration, respectively (Table 1). A comparison of means showed that although the daily value was not significantly better the 2-day frequency, it was significantly better than the 4-day and no aeration treatment (Table 1). Irrigation means were not significantly different. Values were 1.53 kg/plant for irrigating to 70% of field capacity and 1.51 and 1.40 kg/plant for the 80 and 90% of field capacity levels (Table 1).

The one-way ANOVA for marketable yield supported the foregoing results. When arranged in descending order, the marketable yield in kg per plant tended to fall into several groups for the 12 treatment combinations (Fig. 4). A high-value group comprising in order D25A2I70 (1.67 kg per plant) > D10A1I80 (1.65) = D40A1I70 (1.65) > D25A1I90 (1.63); an intermediate group D25A4I80 (1.54) > D25ANI80 (1.50) > D10A2I90 (1.48); a low yield group D10A4I70 (1.40) > D10ANI70 (1.38) > D40A2I80 (1.37) > D40A4I90 (1.34). The lowest yield of 1.16 kg per plant was obtained for D40ANI90. The high values group included all the high frequency (A1) soil aeration. The lower frequency aeration (A2, A4, AN) and the higher irrigation level means (I80 and I90) were associated with the low value group. These results tend to confirm a positive effect of soil aeration.

Fig. 4

Marketable fruit yield of muskmelon measured for each combination of tubing placement depth (Dcm), aeration frequency (A interval in day), and irrigation level (I to% field capacity) treatment factors. Ducan0.05 for combination treatment mean comparison = 0.40 kg

Given a measured 4 to 5 fruits per plant, the foregoing results imply increasing fruit size for the grouping of the treatment combinations. Average fresh fruit weight ranged from 232 to 290 g for the lowest yield of 1.16 ± 0.08 kg for D40ANI90. The range for the low value group average of 1.37 ± 0.20 kg was 275 to 343 g, for the intermediate group with 1.41 ± 0.29 kg was 300 to 376 g, and for the high value group with 1.65 ± 0.43 kg was 330 to 412 g. The overall mean across all treatment factors in Table 1 and treatment combinations was 1.48 kg per plant giving an overall average fresh fruit weight ranging between 300 to 360 g for the Shantian No. 1 cultivar.

Tomato results

Dry matter partitioning of tomato

Figure 5 shows that Post-harvest stem, leaf, and root dry matter of tomato was greater with aeration than without aeration treatment. For both tube burial depths, Post-harvest stem, leaf, and root DM increased with increasing aeration frequency. Maximal values for stem, leaf, and root DM were obtained in the D15A2, D40A2, D40A4 treatments, respectively. Our investigation found that root DM of tomato with a 40-cm burial of the drip irrigation tube was higher than 15-cm treatment. The highest Total dry matter (stem + leaf + root DM) value of 98.4 g was obtained for the D15A2.

Fig. 5

Post-harvest stem, leaf, and root dry matter of tomato measured for each combination of tubing placement depth (Dcm) and aeration frequency (A interval in day) treatment factors. The differences were compared using the Duncan’s test with a significance level of 0.05

Marketable fresh fruit yield of tomato

Soil aeration improved plant root growth environment, thus increasing yield. Fig. 6 shows that the tomato yield was significantly different between aeration treatments. The results are found for the yield as aeration treatments had significantly higher yield than CK. The D15A2 treatment combination resulted in the maximum tomato yield (Fig. 6). The D15A2 was increased the tomato yield by 27.8% compared with the D10 AN.

Fig. 6

Fruit yield of tomato for each combination of tubing placement depth (Dcm) and aeration frequency (A interval in day) treatment factors. The differences were compared using the Duncan’s test with a significance level of 0.05


Sub-surface drip irrigated systems are managed to optimize soil water availability, accessibility, and use efficiency for maximum growth and marketable yield. On the other hand, irrigation always results in displacement and/or redistribution of soil air. This can reduce the level and availability of O2 which is continuously needed to act as the electron acceptor in the tricarboxylic acid metabolic cycle, ATP production, and normal root cell activity [22, 23]. Oxygen-deprivation destroys mitochondria and proteins in root cells leading to cytoplasm acidosis and plant growth inhibition [24]. In this context, it would be expected that any negative growth impact of increasing irrigation application can be offset by soil aeration delivered via the drip-irrigation tubing when placed at a given depth in the soil.

Effect of soil aeration on muskmelon and tomato growth

Hypoxia induces a shift from normal respiration to anaerobic respiration. Previous studies have shown that artificial aeration significantly improved dry matter accumulation, yield, and water use efficiency of muskmelon cotton, wheat, tomato and pineapple [5, 25, 26]. Our results showed that soil aeration had a positive effect on both muskmelon and tomato plant growth. Both muskmelon and tomato DM increased with increasing aeration frequency indicating that rhizospheric oxygen stress do exist, and soil aeration would more likely ameliorate hypoxic conditions.

Applying 405 L of air through the average of 36 emitters would imply soil ventilation at a rate of about 11 L per emitter. It would be expected that most of this air would move upwards by convection displacing existing gases in the air-filled porosity above the plane of the drip lines. In short, the aeration treatments would in effect flush this air-filled porosity. This flushing effect would impact a larger volume of soil with increasing tubing placement depth. Permanent or semi-permanent sub-surface drip irrigation systems need to be placed at depths that permit mechanical surface cultivation. This would reduce the potential of hypoxia. The air-filled porosity in a given soil volume would increase as water is removed by evapo-transpiration. More frequent flushing would increase O2 available for root metabolic activities and promote root growth.

Increasing levels of aeration frequency would supply more O2 for plant aerobic respiration and generally this would positively increase growth and yield. Any positive effect of aeration increasing soil ventilation on root growth would translate into improved yield and growth as reflected in the consistently higher values of yield, LAI, and DM with increasing soil aeration frequency. Aeration frequency significantly increased tomato growth performance and fruit yield (Table 1). As expected, increasing aeration frequency progressively increased dry matter and marketable fresh fruit yield of tomato (Fig. 5, Fig. 6). This study indicate that soil aeration can significantly increase the plant growth performance. It is consistent with previous studies of tomato, pepper and potato [27,28,29,30]. In addition, the aeration treatments tended to significantly impact most responses measured after 25 DAT. Although soil aeration did not significantly affect the LAI55 and LAI75, the leaf DM response was significant. The reason for this is most likely due to uncontrollable factors (such as solar angle, cloudiness, etc.) that increase the variability in the measured LAI by the canopy solar radiation transmittance method.

Effect of irrigation on muskmelon growth

The plots were drip-irrigated at 23 and 60 DAT based on the measured gravimetric water content averaged over the 0–60 cm depth of the soil profile on these dates. Based on the volumetric water content of 0.24 measured at transplanting there would be an additional 100 mm x (0.24–0.13) = 11 mm of plant available water in the 50–60 cm profile depth. This would imply a total of 136 mm when added with the 125 mm made available in the 0–50 cm from the surface irrigation at the time of transplanting. This stored water would support an evapo-transpiration rate of close to 6 mm per day over the 0–23 DAT period. This implies there was very likely no soil water deficits during this period. Any effect of the deficit drip irrigation levels at the different tubing placement depths would therefore be more apparent for responses measured after 25 DAT.

At 23 and 60 DAT the plots were drip irrigated with a volume equivalent to the amount required to replenish the measured soil water storage in the 0–60 cm profile depth to 70, 80, or 90% of field capacity. For the I70 treatment, this would imply making plant available water in the soil wetted volume around the emitters equal to 600 x (0.7 × 0.38–0.13) mm = 82 mm for both the 37-day period between 23 and 60 DAT, and for the 15-day period between 60 DAT and harvest at 75 DAT. Corresponding values were 104 mm for I80 and 127 mm for I90.

The physical and hydraulic properties of the deep, loess-derived clay loam soil used in this study were quite uniform with depth. Therefore, it would be expected that the shape and volume of the wetted soil around each emitter would be similar for all tubing placement depths. The shape and dimensions of the wetted soil around the buried emitters would be important in trying to explaining the growth responses to the subsurface drip irrigation treatments for different tubing placement depths (Table 1 and Figs. 1, 2, 3, 4).

In this study, the emitters were spaced 30 cm along the tubing and the tubing was placed 50 cm apart. The foregoing review indicates that the entire thickness of the soil profile for at most 13 cm above and 13 cm below the plane of the drip line placement depths would be wetted by drip irrigation. This thickness would increase with increasing irrigation level. These literature results relating to a soil with the same texture as that used in this study, along with the non-significant effects of irrigation level on most of the post-25 DAT responses (Table 1), tend to confirm that soil water was not a limiting factor in this experiment. The marketable fruit yield, LAI55, and LAI75 remained unaffected by irrigation level (Table 1, Figs. 1 and 2). For all the measured dry matter accumulation responses at 75 DAT, only the root DM means for the I80 and I90 irrigation levels were significantly lower compared to the I70 level mean (Table 1 and Fig. 2). The irrigation use efficiency for both fruit yield (IUEf) and total DM accumulation (IUEdm) decreased significantly with increasing irrigation level (Table 1). This would imply that irrigation water applied above the I70 level did not contribute to growth and yield since there was no effect of irrigation level on fruit yield and total DM accumulation.

Effect of burial depths of subsurface tubing on muskmelon and tomato growth

Increasing the placement depth for a given irrigation/aeration level would correspondingly change the position of the soil moisture or soil O2 around the emitter. The shape of the wetted volume is expected to change with depth also in homogenous soils due to the effects of transpiration, evaporation and/or ground water. In addition, in soils with marked changes in physical and hydraulic properties with depth increasing the emitter depth would also induce changes in the shape of the wetted volume. Both muskmelon and tomato are moderate-rooted plant, the roots are mainly distributed within 40 cm of the soil surface. The results showed that muskmelon plant DM increased with increasing tubing placement depth (Fig. 2, Table 1). We speculate that lower dry matter accumulation at a shallow depth (10–25 cm) of aeration/ irrigation lines due to the chimney effect. For the root zone soil texture used in this study, the experimental results, in their entirety, tend to support this reasoning. In addition, Wetting the profile 13 cm above and 13 cm below the plane of the drip line placement would imply increasing potential water losses due to soil evaporation for the shallow placement depth with D10 > D25 > D40.

Interactions of soil aeration, burial depths and irrigation on plant growth

Tomato is a moderate rooted plant, the roots are mainly distributed within 40 cm of the soil surface. Until 23 DAT only the soil aeration treatments were applied via the 10, 25, and 40 cm tubing placement depths. Neither irrigation levels or tubing placement depth had any effect on the marginal mean for LAI25. Any effect of these 2 treatment factors would be apparent in the LAI measured at 25 DAT. Nevertheless, the 2-factor D x A interaction was significant affected LAI25 (Table 1). On the other hand, all the soil aeration treatments positively affected LAI25 compared to the no aeration treatment. This suggests that low natural soil aeration may be a limiting factor in the clay loam soil used in this study. It confirms that artificial ventilation would positively affect growth, and provides more support for studies showing beneficial effects of this practice as was reported by Bhattarai et al. [17] in his comprehensive review of soil aeration research. According to our studies we found that only the D x I two-factor interaction was significant affected muskmelon fruit yield (Table 1). This suggests that tubing placement depth determine the position of the soil wetted. The difference of rhizosphere soil moisture results in the plants grow differently. Soil aeration would more likely ameliorate hypoxic conditions above, rather than below, the emitter. If tubing is placed at a shallow depth (10 cm), plant growth and fruit yield might not be improved significantly because of the chimney effect (Table 1). Also, the 2-factor D x A interaction on LAI25 and Root DM were significant indicating that the aeration effect may not be independent of application depth in an initially wet soil profile. On the other hand, the significant 2 factor D x I and A x I interactions on IUEdm indicated that the irrigation level effects varied differently for the tubing placement and aeration treatment factors.

The foregoing discussion raises the question as to whether the roots were able to fully access all the plant available water. Also, how could soil aeration influence the IUE and dry matter partitioning? Unfortunately, there are no recent field investigations on the root system development of muskmelon plants and how the roots explore the soil profile. Assuming a similar overall root architecture and distribution for the Hami melon, it was quite possible that some of the total applied irrigation at the I80 and I90 may not have been accessible to the plants. If so, this would explain the lack of irrigation effect on the growth responses and the significant decreases in IUEf and IUEdm with increasing irrigation level (Table 1, Figs. 1, 2, 3).

In general, all management practices to optimize the root environment can only indirectly influence crop growth and yield through their direct effect on the root systems with limits dictated by the species genetics [31]. Muskmelon and tomato plant roots are especially vulnerable to Hypoxia stress O2 deficiency in the soil reduces adenosine triphosphate (ATP) production leading to lower uptake and transport of nutrients to the shoot leading to a reduction in leaf growth and photosynthesis [32]. The results of this study on muskmelon add new information to the findings of Bhattarai et al. [18, 33,34,35] and Sharma et al. [36] that show positive effects of soil aeration on growth and yield of other crops.


The results of this experiment show that it is feasible to configure installation and operation of commercially-available, permanent or semi-permanent, sub-surface drip irrigation systems to manipulate the root environment and indirectly influence growth and yield of greenhouse-produced muskmelon and tomato. For muskmelon experiment, deeper drip tubing placement permits a wider range of mechanical cultivation practices but can induce O2 deficits in the wetted zone around the emitters especially for low-frequency, high-volume irrigation applications. The Total DM of muskmelon at 75 DAT increased with increasing tubing placement depth and aeration frequency, the D40 was increased the Total DM by 39.3% compared with the D10, and the A1 was increased the Total DM by 46.9% compared with the AN. The fruit yield of muskmelon and tomato increased with increasing aeration frequency, and 1-d aeration intervals was increased the yield of muskmelon by 21.5% compared with the control, and the D15A2 was increased the tomato yield by 27.8% compared with the D10AN. Our results demonstrated that the negative impact of any such deficit on the measured responses can be offset by soil aeration via the drip tubing.


Experimental conditions, setup, and treatments

The muskmelon experiment was conducted in a 108 m long and 8 m wide greenhouse located in Yangling, Shaanxi Province, Northwest China (34°17′N latitude, 108°02′E longitude) from April 24 to July 12, 2014. Tomato experiment was conducted from Oct 182,014 to May 20, 2015 in the same greenhouse. The climate is semi-arid with a long-term average annual precipitation of 550–650 mm, with the average annual sunshine of 2163.8 h and 210 frost-free days. The soil texture of root growth area was clay loam with 25.4% sand, 44.1% silt, and 30.5% clay. pH was 7.82, dry bulk density 1.35 g cm− 3, porosity 49.4%, and gravimetric field capacity 28.2% (volumetric field capacity 38%). The physical and chemical characteristics of the irrigation water: pH was 7.90, EC was 2.71 ds m− 1, total suspended solids were 15 mg L− 1, CODMn was 1.2 mg L− 1.

The main drip irrigation tube was connected to an air pump, both water and air were supplied to the soil through the tube (Fig. 7). Drip irrigation tubes (φ16 subsurface drip irrigation pipe, Qinchuan water-saving irrigation equipment engineering Co. Ltd., Yangling, China), 16 mm diameter with emitter spacing of 0.30 m, were buried to the appropriate depth spaced 0.50 m apart along the ridges. The subsurface drip irrigation tubes had a peculiar labyrinth channel structure, so it can promote the uniformity discharge rate of air/water as much as possible. The drip irrigation tubes were connected to a distribution system designed to supply both water and air. Each plot had 35 to 36 emitters.

Fig. 7

experimental arrangement of an example block

Muskmelon cultivar Shantian No.1 (Northwest New Horizon Facilities Agriculture Development Co. Ltd., Yangling Shaanxi, China) is one of over 100 cultivars of the Chinese ‘Hami’ reticulated melon group (Cucumis melo L.). It is an early maturing, monoecious cultivar with separate male and female flowers on the same plant. Each plant produces 4–5 small (300–350 g), fragrant, sweet, crisp-fleshed, fruits that are very popular with consumers. The tested cultivar of tomato (Lycopersicon esculentum Mill.) was Fen-Yu-Yang-Gang (Northwest New Horizon Facilities Agriculture Development Co. Ltd., Yangling Shaanxi, China).

Before transplanting, the soil was rototilled, and 120 t/ha of decomposed organic manure (pig and sheep manure), 400 kg/ha of compound fertilizer (18% N, 15% P2O5, and 12% K2O), and 1500 kg/ha of diammonium phosphate (18% N and 46% P2Oe) were broadcast uniformly in the soil as the basal fertilizer. After 20 days, 26 plugs of muskmelon (or tomato) were transplanted to the experimental plots spaced 0.40 m apart within 2 rows spaced 0.5 m apart on the ridges. To prevent the lateral spread of air and water into adjacent treatments, the plots were separated from each other by a 1.5 m wide empty space. Post-transplanting management practices (i.e. fertilization, agricultural chemicals spraying, fruit pruning etc.) for all plots were consistent with local production practice. Fruiting was not restricted, although some farmers prune the flowers.

Based on the balance between soil air update rate and labor cost, artificial aeration treatments of tomato were none or aeration at 2-day, and 4-day intervals, and interval of aeration for muskmelon were none or aeration at daily, 2-day, and 4-day intervals beginning the first day after transplanting and designated as Ainterval in day (i.e. A1, A2, A4, and AN). For each treatment, 405 L of air was applied to each plot via the drip tubing using a manifold connected to an air compressor. The flow rate for each plot was about 10.2 L min− 1. The air which injected via subsurface drip irrigation tubes have a characteristic of high O2 concentration and low CO2 concentration compared with air which originally stored in soil pore space. The injected air/oxygen consumed by soil microorganisms, soil animals, crop roots. Moreover, a fraction of oxygen diffuses to the atmosphere because of the chimney effect. This volume was calculated as [5, 37]:

$$ {\mathrm{V}}_{\mathrm{A}}=\left[ SL\left(1-{\rho}_{\mathrm{b}}/{\rho}_s\right)\right]/1000\kern0.5em \eta . $$

where VA was the amount of air injected, S the area of a cross-section between rows, L the length of plot, ρb the soil bulk density, and ρs the soil particle density. η an application efficiency coefficient for aeration system.

Irrigation timing and amount demand of tomato experiment is mainly driven by farmers’ perceptions and climatic conditions. For the muskmelon experiment, following recommended production practices the plots were surface-irrigated at the time of transplanting. Soil water content was measured and controlled using a Field TDR 200 soil moisture meter (Spectrum, Aurora, IL, USA). A 60-cm deep probe was installed in the center of each plot. Soil water content was measured at 10-cm intervals down to a depth of 60 cm. The gravimetric water content (θg) averaged over the 0–60 cm at the time of transplanting was measured. Subsurface drip irrigation treatments were based on the measured gravimetric water content averaged over the 0–60 cm depth of the soil profile on these dates. The location of soil sampling was between two rows of muskmelon (or tomato), each soil sampling with three replicates. Muskmelon treatments, designated as I to% field capacity (i.e. I70, I80, and I90), were based on replenishing the water in the soil volume (Vs) in 60 cm of the soil profile to 70, 80, and 90% of the gravimetric field capacity. The irrigation treatment amount in litre of muskmelon were calculated as: [5, 38]

$$ {V}_I={V}_s\times \uprho b\times \left({\theta}_{gfc}\times {q}_1-{q}_2\right)/\left(\eta \times 1000\right). $$

Where VI was the irrigation amount, ρb the soil bulk density, θgfc field capacity, q1 the irrigation level (0.7, 0.8, or 0.9), q2 the measured soil moisture content, and η an application efficiency coefficient for irrigation system. Irrigation timing demand is mainly driven by farmers’ perceptions and climatic conditions. Details of irrigation timing and amounts during muskmelon growth periods are shown in Table 2. Pre-irrigation soil moisture of muskmelon for different treatments and soil depth at 22 and 59 DAT are shown in Fig. 8.

Table 2 Details of irrigation timing and amounts during muskmelon growth period
Fig. 8

Soil moisture for different treatments and soil depth at (a) 22 days and (b) 59 days after transplanting

Experimental design

Muskmelon experimental design

A full 41 × 32 factorial of the 36 experimental treatment combinations and would require 108 plots for 3 replicates (Table 3). As with all fractional factorial designs, some effects are not estimable and confounding of lower and higher order effects is unavoidable. In this case, the three-factor interaction effect was not estimable. The main effects are independently estimable and it was possible to determine the existence of two-factor interactions. Even with these limitations, the design provided good stability and permitted estimation of the effects of primary interest comparable to the output of the much more expensive complete factorial experiment requiring 108 experimental plots.

Table 3 Muskmelon experimental design

Tomato experimental design

The experiment was limited to 6 treatments arranged as a randomised complete block design with drip irrigation tube burial depth and aeration frequency (Table 4).

Table 4 Tomato experimental design

Plant measurements

At 75 days after transplanting (DAT), all marketable fruits of muskmelon were harvested from each plot. And 214 days after transplanting (DAT), all marketable fruits of tomato (3 trusses) were harvested from each plot. The fruits were weighed and the fresh fruit weight per plant calculated. At 25, 55, and 75 DAT leaf area index (LAI) was measured by using AccuPARLP-80 canopy analyzer (Decagon Devices, Pullman, Washington 99,163, USA) based on solar radiation transmittance. After harvesting, the remaining above-ground portion of all plants in each experimental plot were collected. The roots were removed and washed through a sieve. Stem and leaves and root were oven dried at 60 °C for 48 h according to Shao et al., 2008; Ahmadi et al., 2014; Ramírez et al., 2014 [39,40,41]. Root dry weight was then measured on 3 muskmelon (tomato) plants from each plot.

Data analyses and statistics

The experimental data were organized in Microsoft Excel 2016. The statistical analyses were performed with the SPSS 22 software package (IBM, Armonk, New York). All figures were constructed using the graphing software OriginPro 9.0 (Origin Lab Corporation, One Roundhouse Plaza, Suite 303, Northampton, MA 01060, USA).

Muskmelon data

Data are expressed as the mean. Three replicates were used for each experimental determination. Data were analyzed using a residual test method before statistical analysis, and the data met the assumption of homogeneity of variances and followed normal distribution. Mean differences between treatments were assessed by analysis of variance (ANOVA). Post-hoc pairwise comparisons of the treatment means were performed using Duncan’s multiple range test. Differences were considered significant at the level of 0.05. The main purpose of the ANOVA was to assess the main effects and the 2-way interactions between drip irrigation tube burial depth, irrigation level and artificial aeration frequency on tomato yield, LAI, post-harvest dry matter partitioning, and irrigation use efficiency. The significances were defined at the level of 0.05 and 0.01.

Tomato data

The experimental design was taken as a 2 × 3 factorial with 3 replicates. Data were analyzed using a residual test method before statistical analysis, and the data met the assumption of homogeneity of variances and followed normal distribution, and all treatment means were compared for significant differences using the Duncan’s new multiple-range test at level of P = 0.05.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.



Electrical conductivity


Days after transplanting


Leaf area index


Dry matter

Dcm (i.e. D10, D25, and D40):

Irrigation tubing placement depths were 10, 25, and 40 cm below the surface of the ridge

IUEf or IUEdm :

Irrigation use efficiency as (marketable fresh fruit yield per plant) or as (g total DM per plant) divided by (liter total irrigation water applied)

Ainterval in day (i.e. A1, A2, A4, and AN):

Interval of aeration for muskmelon/tomato were none or aeration at daily, 2-day, and 4-day intervals beginning the first day after transplanting

Ito% field capacity (i.e. I70, I80, and I90):

Muskmelon treatments, designated as, were based on replenishing the water in the soil volume (Vs) in 60 cm of the soil profile (Vs = 5.5 m2 × 0.6 m) to 70, 80, and 90% of the gravimetric field capacity

LAI days (i.e. LAI25, LAI55, and LAI75):

LAI of muskmelon plants measured by canopy solar radiation transmittance at 25, 55, and 75 days after transplanting.


  1. 1.

    Chang J, Wu X, Wang Y, Meyerson LA, Gu BJ, Min Y, Xue H, Peng CH, Ge Y. Does growing vegetables in plastic greenhouses enhance regional ecosystem services beyond the food supply? Front Ecol Environ. 2013;11:43–9.

    Article  Google Scholar 

  2. 2.

    Jiang W, Qu D, Mu D, Wang L. Protected cultivation of horticultural crops in China, In Horticultural Reviews, Volume 30 (ed J. Janick). John Wiley & Sons. Oxford, UK: Inc; 2003. Chapter 4.

    Google Scholar 

  3. 3.

    Wittwer SH, Castilla N. Protective cultivation of horticultural crops worldwide. Hort Tech. 1995;5:6–23.

    Article  Google Scholar 

  4. 4.

    Mostaghami S, Mitchell KJ, Lembke WD. Effect of discharge rate on distribution of moisture in heavy soils irrigated from a trickle source. Trans ASAE. 1981;25:2081–5.

    Google Scholar 

  5. 5.

    Li Y, Niu WQ, Xu J, Wang JW. Zhang MZ, lv W. Root morphology of greenhouse produced muskmelon under sub-surface drip irrigation with supplemental soil aeration Sci Hortic. 2016;201:287–94.

    Google Scholar 

  6. 6.

    Michelakis N, Vougioucalou E, Clapaki G. Water use, wetted soil volume, root distribution and yield of avocado under drip irrigation. Agr Water Manage. 1993;24:119–31.

    Article  Google Scholar 

  7. 7.

    Ben-Noah I, Friedman SP. Review and evaluation of root respiration and of natural and agricultural processes of soil aeration. Vadose Zone J. 2018;17:170119.

    Article  Google Scholar 

  8. 8.

    Cook FJ, Thorburn PJ, Fitch P, Bristow KL. WetUp: a software tool to display approximate wetting patterns from drippers. Irrig Sci. 2003;22:129–34.

    Article  Google Scholar 

  9. 9.

    Smit B, Stachowiak M, Volkenburgh V. Cellular processes limiting leaf growth in plants under hypoxic root stress. J Exp Bot. 1989;40:89–94.

    Article  Google Scholar 

  10. 10.

    Li TL, Chen YD, Liu YL, Qin X. Effects of rhizosphere CO2 concentration on root growth and activity of netted muskmelon. Transactions of the CSAE. 2009;25:210–5.

    Google Scholar 

  11. 11.

    Liu YL, Li TL, Sun ZP, Chen YD. Effects of Rhizosphere CO2 concentration on plant growth and root nitrogen metabolism of muskmelon. Sci Agric Sin. 2010;43:2315–24.

    CAS  Google Scholar 

  12. 12.

    Liu YL, Li TL, Sun ZP, Chen YD. Impacts of root-zone hypoxia stress on muskmelon growth, its root respiratory metabolism, and anti oxidative enzyme activities. Chin J Appl Ecol. 2010;21:1439–45.

    CAS  Google Scholar 

  13. 13.

    Bhattarai SP, Dhungel J, Midmore DJ. Oxygation improves yield and qualityand minimizes internal fruit crack of cucurbits on a heavy clay soil in the semi-arid tropics. J Agric Sci. 2010;2:17–25.

    Google Scholar 

  14. 14.

    Liu Y, Zhou Y, Wang T, Pan J, Zhou B, Muhammad T, Zhou C, Li Y. Micro-nano bubble water oxygation: synergistically improving irrigation water use efficiency, crop yield and quality. J Clean Prod. 2019;222:835–43.

    CAS  Article  Google Scholar 

  15. 15.

    Pendergast L, Bhattarai SP, Midmore DJ. Benefits of oxygation of subsurface drip-irrigation water for cotton in a vertosol. Crop Pasture Sci. 2013;64:1171–81.

    CAS  Article  Google Scholar 

  16. 16.

    Bhattarai SP, Huber S, Midmore DJ. Aerated subsurface irrigation water gives growth and yield benefits to zucchini, vegetable soybean and cotton in heavy clay soils. Ann Appl Biol. 2004;144:285–98.

    Article  Google Scholar 

  17. 17.

    Bhattarai SP, Midmore DJ. Influence of soil moisture on yield and quality of tomato on a heavy clay soil. Proc Int Symp Harnessing Potential Horticulture Asian-Pac Region. 2005:451–4.

  18. 18.

    Bhattarai SP, Pendergast L, Midmore DJ. Root aeration improves yield and water use efficiency of tomato in heavy clay and saline soils. Sci Hortic. 2006;108:278–88.

    Article  Google Scholar 

  19. 19.

    Bhattarai SP, Midmore DJ, Pendergast L. Yield, water-use efficiencies and root distribution of soybean, chickpea and pumpkin under different subsurface drip irrigation depths and oxygation treatments in vertisols. Irrig Sci. 2008;26:439–50.

    Article  Google Scholar 

  20. 20.

    Zhu LF, Yu SM, Jin QY. Effects of aerated irrigation on leaf senescence at late growth stage and grain yield of rice. Rice Sci. 2012;19:44–8.

    Article  Google Scholar 

  21. 21.

    Murchie EH, Pinto M, Horton P. Agriculture and the new challenges for photosynthesis research. Tansley Rev. 2009;181:532–52.

    CAS  Google Scholar 

  22. 22.

    Rieu M, Sposito G. Fractal fragmentation, soil porosity, and soil water properties: I. Theory Soil Sci Soc Am J. 1991;55:1231–8.

    Article  Google Scholar 

  23. 23.

    Biemelt S, Keetman U, Albrecht G. Re-aeration following hypoxia or anoxia leads to activation of the antioxidative defense system in roots of wheat seedlings1. Plant Physiol. 1998;116:651–8.

    CAS  Article  Google Scholar 

  24. 24.

    Drew MC. Oxygen deficiency and root metabolism: injury and acclimation under hypoxia and anoxia. Annu. Rev. plant. Physiol. Plant Mol. Biol. 1997;48:223–50.

    CAS  Google Scholar 

  25. 25.

    Munns R, Sharp RE. Involvement of abscisic acid in controlling plant growth in soils of low water potential. Aust J Plant Physiol. 1993;20:425–37.

    CAS  Google Scholar 

  26. 26.

    Chen XM, Dhungel J, Bhattarai SP, Torabi M, Pendergast L, Midmore DJ. Impact of oxygation on soil respiration, yield and water use efficiency of three crop species. J Plant Ecol. 2011;4:236–48.

    Article  Google Scholar 

  27. 27.

    Bonachela S, Quesada J, Acuna RA, Magan JJ, Marfa O. Oxyfertigation of a greenhouse tomato crop grown on rockwool slabs and irrigated with treated wastewater: oxygen content dynamics and crop response. Agric Water Manag. 2010;97:433–8.

    Article  Google Scholar 

  28. 28.

    Ityel E, Ben-Gal A, Silberbush M, Lazarovitch N. Increased root zone oxygen by a capillary barrier is benefcial to bell pepper irrigated with brackish water in an arid region. Agric Water Manag. 2014;131:108–14.

    Article  Google Scholar 

  29. 29.

    Shahien MM, Abuarab ME, Magdy E. Root aeration improves yield and water use efficiency of irrigated potato in sandy clay loam soil. Int J Adv Res. 2014;2(10):310–20.

    Google Scholar 

  30. 30.

    Li Y, Niu WQ, Wang JW, Liu L, Zhang MZ. Effects of artificial soil aeration volume and frequency on soil enzyme activity and microbial abundance when cultivating greenhouse tomatoes. Soil Sci Soc Am J. 2016;80:1208–21.

    CAS  Article  Google Scholar 

  31. 31.

    Boyer JS. Plant productivity and environment. Science. 1982;218:443–8.

    CAS  Article  Google Scholar 

  32. 32.

    Barrett-Lennard EG. The interaction between waterlogging and salinity in higher plants: causes, consequences and implications. Plant Soil. 2003;253:35–54.

    CAS  Article  Google Scholar 

  33. 33.

    Bhattarai SP, Midmore DJ. Oxygation enhances growth, gas exchange and salt tolerance of vegetable soybean and cotton in a saline Vertisol. J Integr Plant Biol. 2009;51:675–88.

  34. 34.

    Bhattarai SP, Su N, Midmore DJ. Oxygation unlocks yield potentials of crops in oxygen-limited soil environments. Adv Agron. 2005;88:313–77.

    CAS  Article  Google Scholar 

  35. 35.

    Bhattarai SP, Salvaudon C, Midmore DJ. Oxygation of the Rockwool substrate for hydroponics. Aquaponics Journal. 2008;49:29–33.

    Google Scholar 

  36. 36.

    Sharma SP, Leskovar DI, Crosby KM, Volder A, Ibrahim AMH. Root growth, yield, and fruit quality responses of reticulatus and inodorus melons (Cucumis melo L.) to deficit subsurface drip irrigation. Agr Water Manage. 2014;136:75–85.

    Article  Google Scholar 

  37. 37.

    Niu WQ, Fan WT, Persaud N, Zhou XB. Effect of post-irrigation aeration on growth and quality of greenhouse cucumber. Pedosphere. 2013;23:790–8.

    Article  Google Scholar 

  38. 38.

    Pei Y, Bie ZL. Efects of diferent irrigation minima on the growth and physiologica1 characteristics of lettuce under plastic greenhouse. Trans CSAE. 2008;24:207–11.

    Google Scholar 

  39. 39.

    Shao G, Zhang Z, Liu N, Yu S, Xing W. Comparative effects of deficit irrigation (DI) and partial rootzone drying (PRD) on soil water distribution, water use, growth and yield in greenhouse grown hot pepper. Sci Hortic. 2008;119:11–6.

    Article  Google Scholar 

  40. 40.

    Ahmadi SH, Agharezaee M, Kamgar-Haghighi AA, Sepaskhah AR. Effects of dynamic and static deficit and partial root zone drying irrigation strategies on yield, tuber sizes distribution, and water productivity of two field grown potato cultivars. Agric Water Manag. 2014;134:126–36.

    Article  Google Scholar 

  41. 41.

    Ramírez DA, Yactayo W, Gutiérrez R, Mares V, De Mendiburu F, Posadas A, Quiroz R. Chlorophyll concentration in leaves is an indicator of potato tuber yield in water-shortage conditions. Sci Hortic. 2014;168:202–9.

    Article  Google Scholar 

Download references


We acknowledge the support of Jian Xu, Jun Gu and Yahui Guan (Institute of Soil and Water Conservation, Northwest A&F University) for their assistance in the field and laboratory work.


This work is supported jointly by the National Key Research Project of China “13th Five Year Plan” (No. 2016YFC0400202), Natural Science Foundation of China (No. 41807041; No. 41831284; No. 51679205), Natural Science Foundation of Guangdong Province (No. 2018A0303130149), Science and Technology Program of Guangzhou (No.20181002SF0530), Science and Technology Program of Xi ‘an (No. 20193052YF040NS040), Fundamental Research Funds for the Central Universities (GK201903115; 2412018JC002), and the Program of Introducing Talents of Discipline to Universities (B16011). The Funding bodies were not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




Conceived and designed the experiments: LY and NWQ. Performed the experiments: LY and ZMZ. Analyzed the data: LY and WJW. Contributed reagents/materials/analysis tools: LY CXS WJW ZMZ ZZX. Wrote the paper: LY. Revised the manuscript: LY, NWQ, CXS and ZZX. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Wenquan Niu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Niu, W., Cao, X. et al. Growth response of greenhouse-produced muskmelon and tomato to sub-surface drip irrigation and soil aeration management factors. BMC Plant Biol 20, 141 (2020).

Download citation


  • Supplemental soil aeration
  • Leaf area index
  • Dry matter partitioning
  • Muskmelon
  • Tomato