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

Background 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 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. Conclusions 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.

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 O 2 , 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 oxygendeficient 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 greenhouseproduced 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.

LAI of muskmelon
Values in the LAI 25 row of Table 1 corresponded to the tendril elongation growth period. The 55 DAT (LAI 55 ) and 75 DAT (LAI 75 ) measurements were taken during the flowering and at fruit harvest. The overall mean of 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  25 and D 10 placement depth gave LAI 25 > 1 indicating that lowering placement depth could offset the negative impact of lowering soil aeration frequency (Fig. 1).
The highest LAI 55 value of 3.60 during the flowering period was obtained for the D 25 25 tended to be higher for combinations with the lower (A 4 and A N ) soil aeration frequency and especially for the higher irrigation levels (I 90 ) ( Table 1). As indicated the LAI 75 values were very close to the LAI 55 for some combinations indicating little leaf growth after the flowering and during the fruit enlargement growth stages. LAI 75 values were markedly lower than LAI 55 for the D 10 A N I 70 and D 40 A N I 90 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 I 70 irrigation level compared to the I 80 and I 90 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 D 40 A 1 I 70 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.
The ternary plot (Fig. 3) show the range of the ratios for root, stem, and leaf DM of muskmelon to total postharvest 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 D 10 A N I 70 > D 10 A 4 I 80 > D 10 A 1 I 80 treatment combinations (Fig. 3). The treatment cluster of low leaf DM percentages was D 10 A N I 70 < D 10 A 4 I 80 < D 10 A 1 I 80 compared with the cluster of high leaf DM percentages D 25 A N I 80 > D 25 A 1 I 90 > D 40 A 1 I 70 treatment combinations. High root DM percentages was obtained for the D 10 A 4 I 80 > D 10 A 1 I 80 cluster of treatment combinations. These clusters tended to confirm the existence of 2-factor interactions ( Table 1).
The overall irrigation use efficiency for marketable fruit yield of muskmelon (IUE f ) 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 IUE f marginal means. The IUE f marginal mean was 21.5 g fresh fruit per plant per liter for the I 70 irrigation level. Corresponding values decreased linearly by 20 and 40% for I 80 (=17.2) and I 90 (=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 (IUE dm ) 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 IUE dm and all 2-factor interactions were highly significant ( Table 1). The IUE dm value for D 40 (=8.94) was significantly higher than for D 10 (= 6.15) by 45%. Similarly, IUE dm progressively increased with increasing aeration frequency from 6.20 for A N to 6.53, for A 4 , to 7.81 for A 2 , and 9.75 for A 1 . The 57% increase between the values for A N and A 1 , and the 49% increase between A 4 and A 1 , were significant (P < 0.05). IUE dm values decreased with increasing irrigation level. The marginal mean of 9.23 for I 70 differed significantly from to the value of 7.29 for I 80 (20% lower), and 6.20 for I 90 (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 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 D 25  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 D 40 A N I 90 . 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, Postharvest stem, leaf, and root DM increased with increasing aeration frequency. Maximal values for stem, leaf, and root DM were obtained in the D 15 A 2 , D 40 A 2 , D 40 A 4 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 D 15 A 2 .

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 D 15 A 2 treatment combination resulted in the maximum tomato yield (Fig. 6). The D 15 A 2 was increased the tomato yield by 27.8% compared with the D 10 A N.

Discussion
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 O 2 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 subsurface 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 O 2 available for root metabolic activities and promote root growth.
Increasing levels of aeration frequency would supply more O 2 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 LAI 55 and LAI 75 , 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 evapotranspiration 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 I 70 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 I 80 and 127 mm for I 90 .
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, LAI 55 , and LAI 75 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 I 80 and I 90 irrigation levels were significantly lower compared to the I 70 level mean (Table 1 and Fig. 2). The irrigation use efficiency for both fruit yield (IUE f ) and total DM accumulation (IUE dm ) decreased significantly with increasing irrigation level (Table 1). This would imply that irrigation water applied above the I 70 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 O 2 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 D 10 > D 25 > D 40 .

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 LAI 25. 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 LAI 25 (Table 1). On the other hand, all the soil aeration treatments positively affected LAI 25 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 LAI 25 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 IUE dm 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 I 80 and I 90 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 IUE f and IUE dm 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 O 2 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.

Conclusions
The results of this experiment show that it is feasible to configure installation and operation of commerciallyavailable, 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 O 2 deficits in the wetted zone around the emitters especially for low-frequency, highvolume irrigation applications. The Total DM of muskmelon at 75 DAT increased with increasing tubing placement depth and aeration frequency, the D 40 was increased the Total DM by 39.3% compared with the D 10 , and the A 1 was increased the Total DM by 46.9% compared with the A N . 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 D 15 A 2 was increased the tomato yield by 27.8% compared with the D 10 A N. 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 , COD Mn 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. 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% P 2 O 5 , and 12% K 2 O), and 1500 kg/ha of diammonium phosphate (18% N and 46% P 2 O e ) 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. Posttransplanting 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 A interval in day (i.e. A 1 , A 2 , A 4 , and A N ). 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 O 2 concentration and low CO 2 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]: where V A 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. I 70 , I 80 , and I 90 ), were based on replenishing the water in the soil volume (V s ) 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] Where V I was the irrigation amount, ρ b the soil bulk density, θ gfc field capacity, q 1 the irrigation level (0.7, 0.8, or 0.9), q 2 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 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.

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).

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)

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.