Next Article in Journal
Exogenously Applied Cyclitols and Biosynthesized Silver Nanoparticles Affect the Soluble Carbohydrate Profiles of Wheat (Triticum aestivum L.) Seedling
Previous Article in Journal
In Vitro Cytotoxic Potential of Selected Jordanian Flora and Their Associated Phytochemical Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Root Distribution of Tomato Cultivated in Greenhouse under Different Ventilation and Water Conditions

Henan Key Laboratory of Crop Water Use, School of Water Conservancy, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
*
Author to whom correspondence should be addressed.
Plants 2023, 12(8), 1625; https://doi.org/10.3390/plants12081625
Submission received: 16 January 2023 / Revised: 20 March 2023 / Accepted: 30 March 2023 / Published: 12 April 2023
(This article belongs to the Section Horticultural Science and Ornamental Plants)

Abstract

:
Mastering root distribution is essential for optimizing the root zone environment and for improving water use efficiency, especially for crops cultivated in greenhouses. Here, we set up two irrigation amount levels based on measurements of the cumulative 20 cm pan evaporation (Ep) (i.e., K0.9: 0.9 Ep; K0.5: 0.5 Ep), and three ventilation modes through opening the greenhouse vents at different locations (TR: open the roof vents only; TRS: open both the roof and south vents; TS: open the south vents only) to reveal the effects of the ventilation mode and irrigation amount on the root distribution of greenhouse tomato. Six treatments were designed in blocks with the ventilation mode as the main treatment and the irrigation amount as the vice treatment. On this basis, the normalized root length density (NRLD) model of six treatments was developed by considering air environment, soil water and temperature conditions, root length density (RLD) and yield. The results showed that air speed of the TRS was significantly higher than TR and TS (p < 0.01), and the air temperature and relative humidity under different ventilation showed the rule: TR > TS > TRS. There was a significant third-order polynomial function relationship between NRLD and soil depth, and the coefficient of the cubic term (R0) had a bivariate quadratic polynomial function relationship with irrigation amount and air speed (determination coefficient, R2 = 0.86). Root mean square errors of the simulated and measured value of NRLD under TR, TRS and TS were 0.20, 0.23 and 0.27 in 2020, and 0.31, 0.23 and 0.28 in 2021, respectively, normalized root mean squared errors were 15%, 17%, 20% in 2020, and 23%, 18% and 21% in 2021. The RLD distribution ratio from the ground surface to a one-quarter relative root depth was 74.1%, and 88.0% from the surface to a one-half relative root depth. The results of the yield showed that a better combination of ventilation and irrigation was recommended as TRS combined with K0.9.

1. Introduction

With the rapid development of modern agriculture, the cultivation technology of solar greenhouses was used to improve vegetable yield, quality and water use efficiency, and had gradually attracted much attention from farmers, especially in northern China [1]. Many greenhouses have a semi-enclosed structure, in which crops are greatly influenced by meteorological factors and the soil environment [2,3]. The root system is an important organ for the transfer of soil moisture and crop water; mastering the root distribution characteristics and water absorption law can well explain the response mechanism of crops to water and environmental stress.
Crop root distribution characteristics in greenhouses have been reported in some previous studies. For example, the root distribution characteristics of pepper and tomato grown in greenhouse were studied by Zapata-Sierra et al. (2021) [4], who indicated that 90% of the root system was concentrated in the range of 0–9 cm (pepper) and 0–11 cm (tomato). Different irrigation amounts would lead to different water conditions around the roots [5,6], thus affecting the growth and distribution of crop roots. Insufficient irrigation would inhibit tomato root growth, while excessive irrigation may lead to a waste of water resources [7,8]. Wu et al. (1999) [9] studied various crop root distribution characteristics and found that the relative root length density and relative soil depth can be described as probabilistic functions. Zuo et al. (2006) [10] built a relative root length density distribution model for winter wheat by using the normalization method. Thereafter, Ning et al. (2019) [11] compared the normalized model with the exponential distribution model, linear distribution model and piecewise distribution model, and indicated that the root length density calculated by the normalized model was in good agreement with the measurements, with an error of 4.0%. Some studies have also reported that the normalized method can better express the spatial distribution characteristics of roots, and calculate root length density based on these theories. However, these studies mainly focus on crops grown under field conditions, and few reports focus on vegetables grown in greenhouse conditions.
In recent years, some studies have tried to analyze the root distribution under different ventilation modes, and have concluded that ventilation is also one of the important factors affecting root growth in greenhouse conditions. It is known that crop roots need soil oxygen to maintain respiration and growth properly [12] and that most oxygen (O2) is obtained directly through diffusive gas exchange from the atmosphere to the inter-root soil. Ventilation can increase the concentration of oxygen in greenhouse. Ventilation can also reduce the physiological stress of the crop by changing the temperature and humidity in the greenhouse [13]. At the same time, suitable temperature and humidity would reduce the production of crop pests and diseases [14], hence promoting crop growth. In addition, the effect of different drip irrigation amounts on tomato root distribution was also studied [15,16], and the results indicated that mild water stress favored the growth of deep root. Ullah et al. (2021) [17] reported that root length and root surface area can significantly increase by reducing the irrigation amount. Mild water stress could promote root growth and plant development, which is beneficial for improving water use efficiency [18].
One study found that ventilation and irrigation affected crop root growth and distribution by altering the greenhouse environment and soil structure, and suitable ventilation mode promoted root growth [19], while another study found that excessive irrigation led to poor root zone permeability, thus limiting root development [20]. In addition, crop root growth indicators, such as root length, root diameter, root surface area and root volume, also changed with ventilation and irrigation management. Previous studies on root systems mainly focused on the distribution patterns of crop root length density in the profile [21]. Root length density represents the total length of the root system per unit volume of soil and reflects the number of roots. Relative root length density represents the proportion of root length distributed at different relative depths and it is used to describe the relative distribution of root length density. Novak (1994) [22] built an exponential distribution model of root length density by analyzing the horizontal and vertical distribution of the underground root system. Yang et al. (2009) [23] established the relative root length density distribution model of wheat by polynomial fitting and verified it through measured and simulated water content with high accuracy. Zuo et al. (2004) [24] established a nonlinear equation with four parameters to simulate the root distribution. Wu et al. (1999) [9] expressed the plant root distribution in terms of a third-order polynomial equation. However, these models only described the distribution of root systems in the horizontal and vertical directions [25,26] and had not considered the test conditions. Therefore, test variables were added to the simulation of root distribution to obtain a more accurate model of root distribution.
From the above research results, we know that ventilation and irrigation amount affect root growth simultaneously in greenhouse conditions. Crop root growth and distribution also affect fruit yield and water use efficiency [27]. However, few studies have focused on the interaction effect of ventilation and irrigation on root distribution, and root length density models under interactive conditions are also lacking. Therefore, two years of studies were conducted in a solar greenhouse to achieve the following objectives: (1) analyze the variations in soil water content, temperature and meteorological factors under different ventilation and irrigation amounts; (2) explore the root distribution and root length density of tomato with drip irrigation under different combinations; (3) establish the relative root length density distribution model and evaluate its performance; and (4) explore the effect of an optimal combination of ventilation and water conditions on tomato yield.

2. Results

2.1. Meteorological, Soil Water and Temperature Condition

2.1.1. Variations in Meteorological Factors

Variations in air speed inside the greenhouse are shown in Figure 1; the right panel indicates the daily average air speed during the whole growth stages, the left panel indicates the variation in air speed with the number of days after transplanting. We found that the magnitude of air speed was TRS > TS > TR, and it varied between 0.06 and 0.27 m s−1 in the TR treatment, 0.06 and 0.72 m s−1 in the TRS treatment, and 0.06 and 0.36 m s−1 in the TS treatment in the two study years. The difference in air speed for the three ventilation treatments began to increase significantly 60 days after transplanting, and also significant differences between TRS and the other two treatments were found in our study (Table 1), which was entirely due to the different opening states of these vents.
Variations in Ta and RH under the three ventilation treatments are shown in Figure 2; the right and middle panels indicate the average Ta and RH during the whole growth stages, and the left panel indicates the cumulative Ta and RH with the number of days after transplanting, for the whole growth stages, and the opposite trends for daily Ta and RH are obvious [28]. The average daily Ta and RH of the TRS treatment were highly significantly different from the TR and TS treatments (p < 0.01) in the middle of the day during the ventilation period (9:00–17:00) for the two seasons. The daily Ta and RH of the TRS were 28.06 °C and 59.36%, which were 6.4%, 5.1% and 7.4%, 5.5% lower than the TR and TS, respectively. The cumulative temperature was closer between TR and TS but ~6.48% and ~5.38% higher than TRS in 2020 and 2021, respectively. The daily Ta and RH of TRS were 28.06 °C and 59.36%, which were 6.4% and 5.1% and 7.4% and 5.5% lower than those of TR and TS, respectively. The maximum RH (96.61%) during both growth periods occurred in TR, which was 1.07% and 1.12% higher than that TRS and TS, respectively.

2.1.2. Variations in Soil Water and Temperature

Crop root development was largely affected by soil water and heat conditions [29,30]. Soil temperature (ST), soil water content (SWC) and the ratio of water to heat ratio (SWC/ST) are shown in Figure 3, where SWC/ST was calculated by the ratio of the average volumetric water content to soil temperature within the 0–20 cm soil layer. A general opposite trend between ST and SWC was observed during the whole growth stage. The average SWC/ST of TRS was 0.87% °C−1 in K0.9, which was 11.2% and 13.3% lower than TR and TS, respectively. However, for TRS, the ST and SWC of the two irrigation treatments (K0.9 and K0.5) were significantly different, with an overall higher soil temperature and water content in K0.9 treatment than K0.5 treatment, while the difference of SWC/ST between K0.5 and K0.9 was not significant.

2.2. Distribution of Root Length Density

The distribution of RLD was directly affected by the soil environment. The tomato RLD distribution of each treatment is shown in Figure 4. RLD under the same ventilation treatment showed that K0.9 > K0.5, and RLD within rows was higher than between rows. K0.9 and K0.5 within rows increased by 19.9% and 18.5% (TR), 15.0% and 21.3% (TRS) and 19.0% and 14.9% (TS) compared to between rows. The RLD under TR and TS varied from 10.44 to 15.20 cm cm−3 and 5.00 to 9.94 cm cm−3, respectively, while TRS was between 6.55 and 7.58 cm cm−3 under the K0.9 and K0.5 treatments, respectively. There was a larger change between TR and TS, while the RLD of TRS was similar under two different irrigation treatments. A similar phenomena can be found where the RLD under the three ventilation treatments decreased gradually with the deepening of the soil layer below 10 cm soil depth.

2.3. Development of the NRLD Model

Equation (7) was used to fit the relationship between relative RLD and soil depth under different treatments in 2021, resulting in correlation coefficients shown in Table 2. R2 was greater than 0.96 for all treatments, and the F−statistics reached a highly significant level (p < 0.01). There was a highly significant correlation between soil depth and RLD, indicating that the function was good. Correlation coefficients for each treatment (R0, R1 and R2) can represent the trend of curve and varied from −97.16 to −79.01, 128.15 to 152.77 and −72.73 to −66.67, respectively. Under the same ventilation mode, R0 and R2 decreased and then increased, while R1 increased then decreased with the increase in irrigation amount. Under the same irrigation amount, R0 and R2 increased then decreased, while R1 decreased then increased with the increase in air speed. R3 represented the relative RLD value at the relative depth Zr = 0. As no root system existed at soil surface (0 cm), it could be regarded as a theoretical value and did not have practical significance.
The ventilation effect from 10:00 to 16:00 was used as the research object. A generalized NRLD distribution model was established according to the irrigation amount (K) and air speed (S), and then the related parameters were obtained. We found a highly significant correlation between R1, R2, R3 and R0 (p < 0.01) (Figure 5). With K and S as independent variables and parameter R0 as the dependent variable over the whole growth stages, the binary quadratic polynomial function could be fitted by using the Levenberg–Marquardt method combined with the general global optimization algorithm by using 1stOpt software (R2 of 0.86), which could be derived under different irrigation amounts and air speeds in the equation for parameter R0.
R 0 = 174.502 1.343 K + 30.188 S + 3.411   ×   10 3 K 2 0.924 S 2 0.00717 KS
where K is the irrigation amount, mm; and S is the average air speed over the whole growth stage, m s−1.
Based on the relationship equation between parameters R1, R2, R3 and R0, and the calculated Equation (1), expressions of the one-dimensional NRLD distribution model for drip-irrigated greenhouse tomatoes were obtained as follows:
N R L D = R 0 Z r 3 + 1.3968 R 0 + 19.053 Z r 2 + 0.546 R 0 24.91 Z r 0.0509 R 0 + 7.7351 R 0 = 174.502 1.343 K + 30.188 S + 3.411 × 10 3 K 2 0.924 S 2 0.00717 K S

2.4. Validation and Application of the NRLD Distribution Model

2.4.1. Validation

The NRLD distribution model was validated by using measured data in 2020 and 2021. The simulated values of relative RLD at each relative depth (Zr) were obtained by substituting the irrigation amount and average air speed into Equation (2). The measured and simulated values were compared to obtain a 1:1 histogram at three ventilation treatments (data points = 72) (Figure 6) and a plot of NRMSE statistics (data points = 12) (Table 3). Under different irrigation amounts, the simulated values of three ventilation treatments agreed well with the measured values, and the RMSE under TR, TRS and TS was 0.20, 0.23 and 0.27 in 2020, and 0.31, 0.23 and 0.28 in 2021, respectively. Table 3 shows the NRMSE between the simulated and measured values of the relative RLD in the two study years. The NRMSE between the simulated and measured values under TR, TRS and TS was 19%, 18% and 21%, respectively, indicating that the performance of the model was perfect.

2.4.2. Application

The relative RLD can be simulated by using Equation (2) when the maximum root depth (Zmax) and RLD are known. A relative RLD distribution model was used to accurately estimate the ratio of root length to total root length (RL/TRL) under different relative sampling depths by using two irrigation amounts and three ventilation treatments in 2020 (irrigation amount: K0.9 = 247.5 mm, K0.5 = 137.5 mm; average air speed: TR = 0.092 m s−1, TRS = 0.152 m s−1, TS = 0.115 m s−1). The relative RLD and RL/TRL were estimated and compared with the measured data (Table 4). The drip-irrigated tomato roots were mainly distributed in the upper soil layer (0–20 cm), which was similar to the research results of Li et al. (2020) [31]. Root length from the surface to a one-quarter depth of relative root system accounted for 74.1% of total root length, while that from the surface to a one-half depth of relative root system accounted for 88.0%. Under the same irrigation amount, RL/TRL increased then decreased with the increase in air speed. Under the same ventilation, RL/TRL increased gradually with the increase in irrigation amount, the difference between the simulated and measured values was greater in the high irrigation amount than in the low irrigation amount, and measured values were lower than the simulated values. Using 2020 as an example, the root distribution under different ventilation and irrigation amount with pan evaporation coefficients of 0.6, 0.8 and 1.0 was predicted (Table 5). The one-dimensional NRLD distribution model can clarify the proportion of root distribution at different soil layers under different ventilation and irrigation amounts, which is beneficial to guide the middle and late ventilation and irrigation management of greenhouse tomatoes.

2.5. Root System and Yield

The yield and RLD for different treatments in 2020 and 2021 are shown in Table 6. The TRSK0.9 treatment has the highest yield (147.6 t ha−1 and 148.1 t ha−1), followed by the TSK0.9 treatment (143.1 t ha−1 and 144.5 t ha−1). The lowest is the TRSK0.5 treatment (119.4 t ha−1 and 126.3 t ha−1). At the highest yield, the root length density is not the highest, but it is at the middle level of the six treatments. Considering the tomato yield, root length density, irrigation amount and ventilation, the combination of TRS and K0.9 is beneficial to increase tomato yield.

3. Discussion

3.1. Effect of Ventilation Mode and Irrigation Amount on Greenhouse Environment

The greenhouse environment was significantly affected by the ventilation mode and irrigation amount [32,33,34]. There was a highly significant difference (p < 0.01) in air speed under different ventilation conditions. In our study, the air speed under TRS was significantly higher than TR and TS. The magnitude of the indoor air speed was a good measure of the water vapor diffusion capacity of the greenhouse [35]. Ta and RH were also influenced by the ventilation mode, which showed that TR > TS > TRS in our study. The air speed regulated the Ta and RH around the leaves, as well as the rate of water vapor diffusion, which had a significant impact on the internal environment of the greenhouse [36]. For instance, Ta and RH of TRS were significantly lower than TR (Table 1). ST and SWC, as important indicators of the soil environment, were influenced by the irrigation amount. The ASWC/ST in the 0–20 cm soil layer of TR was significantly higher than TRS, because air speed had a greater effect on crop evapotranspiration and low air speed slowed the soil water depletion process.

3.2. Effect of Ventilation Mode and Irrigation Amount on RLD Distribution

The tomato roots were dominantly concentrated in shallow soils (0–20 cm) in the vertical direction, which is the same as Liu et al. [37], who indicated that the tomato root system had an initial growth at a faster rate than the infiltration stress at the soil surface, and the root system extended deeper and further, with a maximum root length density in the range of 10–20 cm. However, Zapata-Sierra et al. (2021) [4] showed that tomato roots were mainly distributed in the 0–11 cm or 15.5 cm soil layer, which may be due to different sampling accuracies. Most studies have found that the maximum root depth was in the range of 60–80 cm [8,38], which was similar to the present study, i.e., the maximum rooting depth of greenhouse tomatoes could be up to 80 cm [39]. Wang et al. (2021) [40] showed that for drip-irrigated tomatoes, the RLD decreased with decreasing irrigation, and the RLD under the same ventilation in this study was K0.9 > K0.5. In this study, root length increased then decreased from the surface to one-half of the relative root depth, with 85.8%, 88.5% and 89.7% for TR, TRS and TS, respectively, indicating that excessive air speed did not increase root amount at one-half, and suitable ventilation is conducive to the growth of tomato roots. Similar findings have been proven by Ge et al. (2019) [41]. The reason was primarily due to excessive air speed causing the crop stomata to close, affecting transpiration and preventing the roots from absorbing water and nutrients properly. A higher proportion of deep soil root distribution in better plants would be able to absorb nutrients from deep soil to improve yield [42,43]. Shu et al. (2020) [44] analyzed the effect of different drip irrigation rates on tomato roots and found that moderate deficit irrigation was beneficial to root “active rooting”, which promoted the development of aboveground parts by absorbing water and nutrients from soil, increasing yields and harvests and improving water use efficiency.

3.3. Performance of NRLD Distribution Models in Different Ventilation Mode and Irrigation Amount

The NRLD of greenhouse tomato under different ventilation and irrigation amounts was constructed in this study, and the model parameters R0, R1 and R2 can reflect the simulation accuracy. R3 is a theoretical value of the NRLD at the surface, and other parameters accurately describe the accuracy of root development fit in soil. Wu et al. (1999) [9] used the normalization method to establish a third-order polynomial function for the relative RLD of wheat and maize, and the function was able to accurately simulate root variation in different soil layers, which was similar to our study. Use of normalization had a good applicability for root modeling. Compared with the maize root system model established by Zou et al. (2018) [28], the parameters R0 and R2 fitted in this study were small, while the parameters R1 and R3 were large. The mean value of R3 was 12.19, which was higher than the value of 5.51 fitted by Zou et al. (2018) [28]. The reason for the difference were that the root system data collected by Zou et al. were all roots, while the tomato root system in our study was partial roots, resulting in the differences in model parameters. However, the R2 of the binary quadratic polynomial fit (Equation (1)) under different ventilation modes and irrigation amounts was 0.86, higher than the fitting accuracy of Zou et al. (R2 of 0.84), which may be caused by differences in crop type, root sampling depths and sample numbers.

3.4. Application of the NRLD Distribution Model

The RLD distribution characteristics of greenhouse tomatoes and the estimated proportion of root distribution in different soil layers were obtained by using the NRLD distribution model, quantifying the influence of ventilation and irrigation on root growth. Therefore, we provide a theoretical basis for the study of root growth models, soil water and solute transport. Previous studies indicated that relative RLD distribution models were better applied to other models for simulating soil water transport and other processes to elucidate crop water uptake patterns. For example, Ning et al. (2015) [45] developed an NRLD distribution model for sunflower to validate and estimate the RLD distribution, which was applied to the Hydrus model to simulate water transport with good results and improve the accuracy of soil water and salt transport simulation. Zuo et al. (2004) [24] developed a model for the NRLD distribution of wheat and applied the measured water content profile to estimate the relative RLD distribution, which provided a basis for the continuous simulation of soil water transport. Compared to the traditional method of obtaining RLD, this model only needs the maximum rooting depth and average RLD at any soil layer. Using the NRLD distribution model, RLD data could be obtained, and real-time information on the proportion of roots in different soil layers is important for accurate ventilation and irrigation management. In this study, NRLD distribution models were established using different ventilation and water treatments to estimate the proportion of root distribution under different treatments. The results showed that the proportion of root length in the upper half of the root zone for K0.9 and K0.5 averaged 88.9% and 87.1%, respectively, and the simulated values were close to the measured values. Combining the 2020 ventilation data and setting evaporation pan coefficients of 1.0, 0.8 and 0.6 for these three irrigation treatments, the results show that reducing irrigation leads to an increase in the proportion of deep roots, and conversely, an increase in the proportion of shallow roots. Combined with the NRLD distribution model, the proportion of greenhouse drip-irrigated tomato roots distributed in different soil layers could be determined in real time, providing a reference for optimizing water and environmental management strategies for crops in middle and late periods. For instance, when the proportion of shallow soil roots is high, the irrigation amount could be reduced, and ventilation could be controlled in an appropriate amount to save water and increase yields.

3.5. Effect of Ventilation and Irrigation Amount on Tomato Yield

Ventilation changes the micro-environment in which crops grow, and also affects the temperature and humidity of the air layer and soil layer, and micro-environment varies depending on the location and size of the vents in the greenhouse (Figure 1, Figure 2 and Figure 3). The differences in tomato yields under the three ventilation treatments are not significant (Table 6), which could be due to the high temperature in the greenhouse affecting colored fruits of the crops, and did not have an effect on crop yield [46]. Irrigation treatments had a significant effect on the yield for the greenhouse tomatoes. Under the same ventilation treatment, the yield for the high irrigation treatment (K0.9) was significantly higher than the low irrigation treatment (K0.5), which is the same as the results of Li et al. (2020) [47], in which a water deficit would be caused under low irrigation treatment, leading to difficulties in water uptake by plant roots, limiting nutrient translocation, affecting carbohydrate synthesis and fruit quantity and quality, as well as affecting plant physiology and reducing plant yield. In addition, a low irrigation amount would lead to a reduced stomatal closure and photosynthetic capacity of plants, and then result in a reduced substance synthesis and lower yields [48]. In our study, a combination of TRS and K0.9 treatments was the most beneficial for greenhouse tomato yield.

4. Materials and Methods

4.1. Experimental Site and Design

The experiment was conducted from March to July in 2020 and 2021 in a solar greenhouse at the Xinxiang Integrated Experimental Base of the Chinese Academy of Agricultural Sciences (35°9′ N, 113°5′ E, 78.7 m above sea level). The greenhouse walls are made of brick and concrete construction, with an area of 510 m2. The direction of the greenhouse is east–west, and it sinks 0.5 m. A steel frame construction is used to support the roof of the greenhouse and is covered with a 0.2 mm thick polyethylene drip-free film to maintain the air temperature inside. Meanwhile, 5 cm thick insulation quilts are used to maintain warmth at the seedling stage. There are three vents, one on the roof (60 m × 30 cm) and another on the bottom of the south side (60 m × 1.5 m). The soil in the greenhouse at a depth of 0–100 cm is a silt loam, including 16.3% clay, 77.1% silt and 6.6% sand. Mean field water capacity and wilting water content at a soil depth of 0–100 cm are 0.31 and 0.11 cm3 cm−3, respectively, with an average bulk density of 1.59 g cm−3.
Tomato seedlings “Solanum lycopersicum L. c.v. Jinpeng M6” were transplanted on 4 March 2020 and 7 March 2021. The size of the experimental plot was 8.0 m long and 2.2 m wide, and a wide (65 cm) and narrow (45 cm) row planting mode was used, with an interval of 30 cm. The planting density was 5.7 trees m−2. Drip irrigation was used, with a drip head flow rate of 1.1 L h−1. Each plot was replicated three times, and six treatments were designed in blocks with the ventilation mode as the main treatment and the irrigation amount as the vice treatment. Here, three ventilation treatments were set: TR (open the roof vents only), TRS (open both the roof and south vents), TS (open the south vents only), and two irrigation treatments were set according to the cumulative water evaporation (Ep) from a standard 20 cm evaporation pan (20 cm diameter and 11 cm deep): K0.9 (0.9 Ep) and K0.5 (0.5 Ep). The irrigation events were performed based on the average surface evaporation of the three ventilation treatments. According to the results of our previous research [49], soil water content of 0–60 cm accounts for 80–90% and 60–65% of field water capacity, respectively, under an irrigation amount of 0.9 Ep and 0.5 Ep. The evaporation pan was placed 30 cm above the crop canopy and adjusted according to tomato plants’ development. The evaporation water amount was measured every day at 8:00 using a measuring cylinder with an accuracy of 0.1 mm, and 20 mm of distilled water was refilled after the measurement. When the accumulated water evaporation reaches 20 ± 2 mm, irrigation events was conducted [50]. The irrigation amount (Ir) was calculated according to Equation (3).
I r = E p × φ
where Ir is the irrigation amount, mm; Ep is the accumulated evaporation, mm; and φ is the water surface evaporation coefficient.
A water meter with an accuracy of 0.001 m3 was installed at the head of each plot to precisely control the irrigation amount. A supplementary irrigation amount of 20 mm was performed by drip irrigation after transplanting to maintain the seedlings alive. In this study, 112 kg hm−2 urea (containing 46% N), 150 kg hm−2 potassium sulfate (containing 50% K2O), and 120 kg hm−2 superphosphate (containing 14% P2O5) were used as base fertilizers and ploughed to a depth of 16 cm with a rotary spade. Thereafter, differential pressure fertilizer tanks were used for topdressing urea at 18.8 kg hm−2 and potassium sulfate at 25 kg hm−2. The fertilizer was applied four times when the first, second, third and fourth truss fruits began to expand. The agronomic practices (e.g., topping, spraying, fruit thinning) were the same as those used locally. The irrigation amounts were 247.5 and 245.7 mm (K0.9) and 137.5 and 136.5 mm (K0.5) in 2020 and 2021, respectively.

4.2. Measurement

4.2.1. Environmental Factors

Air speed in the greenhouse was monitored by using an anemometer (Wind Sonic, Gill, UK) at the vents with an accuracy of ±0.02 m s−1. Data were collected every 5 s, and the 15 min average was recorded in a CR1000 data logger (Campbell Scientific Inc., Logan, UT, USA). The air temperature and relative humidity were measured by using an automatic climate station (CS215, Campbell Scientific, Inc, Monterrey, CA, USA) with accuracies of 0.02 °C and 0.18 °C, and 15 min averages were calculated and stored.

4.2.2. Soil Water Content

A ZL6 cloud data collector (METER Group, USA) with an accuracy of 1 × 10−6% was used to determine the water content of the soil layer at 0, 10, 20, 30, 40 and 60 cm in the middle of two drip heads of the same drip tape [51], with data collected automatically every 15 min.

4.2.3. Soil Temperature

An eight-channel soil temperature logger (JL-04, Ningbo, China) with an accuracy of 0.1 °C was used to monitor soil temperatures at 0, 10, 20, 30, 40 and 60 cm, with data collected automatically at 30 min intervals.

4.2.4. Root Distribution

At the end of picking, root systems were removed from the soil depth between 0 and 60 cm at different orientations by using a layered segmental soil auger with an auger diameter of 7 cm. The soil enclosing the roots was removed at intervals of 10 cm within and between rows (Figure 7). Tomato roots were placed in mesh bags, rinsed and scanned into JPG files using a 4800 (H) × 9600 (V) dpi (MRS-9600TFU2L, WANSHEN, China) scanner and analyzed for morphological characteristics, root length and surface area image were analyzed by software (Win RHIZO Pro2004 b, Canada). Root length density (RLD, cm cm−3) was calculated according to Equation (4):
R L D = R L / V
where RLD is the root length density, cm cm−3; RL is the root length in different soil layers, cm and V is the root auger volume (384.85 cm3).

4.2.5. Yield

Twenty plants were selected in the middle of each plot to measure yield, and this was repeated 3 times. When the tomato fruits were picked, an electronic balance with an accuracy of 0.005 kg was used to weigh the tomato and record the number to calculate tomato yield.

4.3. Root Length Density Distribution Model

To facilitate the modeling of tomato root length density under different treatments, the normalization method [9] was used to convert the penetration depth of tomato roots at different soil layers to a standardized root depth in the range of 0 to 1. The root length density was expressed as a generalization of the relative root length density by using the following equation:
Z r = Z i / Z m a x
N R L D ( Z r ) = R L D ( Z r ) 0 1 R L D ( Z r ) d Z r
where Zr is the standardized root depth, between 0 and 1, dimensionless; Zi is the depth of the rooted soil layer, cm; Zmax is the maximum rooting depth, cm, obtained for the root-free soil layer, and the maximum rooting depth in this experiment is 80 cm; NRLD (Zr) is the relative root length density value, dimensionless; and RLD (Zr) is the root length density value at Zr, cm cm−3.
A third-order polynomial was used to fit the mean values of NRLD (Zr) at each lateral position at the relative sampling depth (Zr) based on previous studies [9,52]. The equation is as follows:
NRLD ( Z r ) = R 0 Z r 3 + R 1 Z r 2 + R 2 Z r + R 3
where R0, R1 and R2 are polynomial parameters and R3 represents the theoretical value of the NRLD at the surface.

4.4. Model Evaluation

The regression fitting method was used to model the relative root length density distribution for the 2021 root data, and the model was validated using the 2020 measured data. Microsoft Excel 2010 and SPSS 26.0 were used for data processing and graphing. Model simulations were carried out using 1stOpt software, and the model was evaluated using the coefficient of determination (R2), root mean square error (RMSE) and normalized root mean squared error (NRMSE).
R 2 = i = 1 N Y i Y ¯ Y ^ i Y ^ ¯ 2 i = 1 N Y i Y ¯ 2 i = 1 N Y ^ i Y ^ ¯ 2
RMSE = i = 1 N ( Y i Y ^ i ) 2 / N 1 2
NRMSE = RMSE Y
where N is the number of samples; Yi is the measured values; Y ^ i is the simulated values; Y ¯ is the measured mean value; and Y ^ ¯ is the simulated mean value. R2 is close to 1, the better the correlation is; RMSE can indicate the average difference between the simulated and observed values, and the closer it is to 0, the smaller the deviation is. NRMSE indicates how good the model simulation performance is; when NRMSE < 10%, the model simulation performance is excellent; when 10% ≤ NRMSE < 20%, the model simulation performance is good; when 20% ≤ NRMSE < 30%, the model simulation performance is average; and when NRMSE ≥ 30%, the simulation performance is considered poor.

5. Conclusions

In this study, we investigated the effects of different ventilation and irrigation amounts on air environment, soil water and temperature conditions and the root distribution of greenhouse tomatoes under drip irrigation. The main conclusions were as follows.
The ventilation and irrigation amount had a significant effect on soil water, temperature and meteorological factors. Ventilation mainly had a significant effect on meteorological factors. TRS was more beneficial to greenhouse gas exchange than TR and TS. Ta and RH were affected by the ventilation treatments in the following manner: TR > TS > TRS. The ventilation and irrigation amount had an interaction effect on soil water and temperature.
The root distribution differences were substantially influenced by both the irrigation amount and ventilation. However, consistent trends of decreasing root length density with deepening of the soil layer were observed below the 10 cm level. Thereafter, a relative root length density distribution model was established according to the following principles: firstly, the relationship between NRLD and relative depth of soil profile was a third-order polynomial; secondly, the coefficient of cubic term (R0) had a bivariate quadratic polynomial relationship with the irrigation amount and air speed. The RMSE between the simulated and measured values of NRLD at K0.9 and K0.5 were 0.20, 0.23 and 0.27 in 2020 and 0.31, 0.23 and 0.28 in 2021, respectively, indicating that performance of the model was perfect.
The irrigation amount and ventilation had significant effects on RLD and yield. The combined treatment of TRS and K0.9 was the most beneficial to increase tomato yield.

Author Contributions

J.G. conducted the experiments and contributed to writing of the manuscript. H.L. designed the study, analyzed the data and contributed to writing of the manuscript. X.G. and Z.Y. contributed to preparation of figures and tables. L.L. and Y.L. reviewed and edited the manuscript. J.G., X.G. and L.L. supervised the research project. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (51709110, 51809094, 51779093), Henan Provincial Higher Education Institutions Young Backbone Teacher Training Program Project (2020GGJS100) and Henan Provincial Key Science and Technology R&D and Promotion Program (202102110128).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhen, H.; Gao, W.; Jia, L.; Qiao, Y.; Ju, X. Environmental and economic life cycle assessment of alternative greenhouse vegetable production farms in peri-urban Beijing, China. J. Clean. Prod. 2020, 269, 122380. [Google Scholar] [CrossRef]
  2. Wei, Z.; Du, T.; Li, X.; Fang, L.; Liu, F. Interactive Effects of Elevated CO2 and N Fertilization on Yield and Quality of Tomato Grown under Reduced Irrigation Regimes. Front. Plant Sci. 2018, 9, 328. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ghannem, A.; Ben Aissa, I.; Majdoub, R. Effects of regulated deficit irrigation applied at different growth stages of greenhouse grown tomato on substrate moisture, yield, fruit quality, and physiological traits. Environ. Sci. Pollut. Res. 2021, 28, 46553–46564. [Google Scholar] [CrossRef] [PubMed]
  4. Zapata-Sierra, A.; Moreno-Pérez, M.; Reyes-Requena, R.; Manzano-Agugliaro, F. Root distribution with the use of drip irrigation on layered soils at greenhouses crops. Sci. Total. Environ. 2021, 768, 144944. [Google Scholar] [CrossRef] [PubMed]
  5. Ji, T.; Guo, X.; Wu, F.; Wei, M.; Li, J.; Ji, P.; Wang, N.; Yang, F. Proper irrigation amount for eggplant cultivation in a solar greenhouse improved plant growth, fruit quality and yield by influencing the soil microbial community and rhizosphere environment. Front. Microbiol. 2022, 13, 981288. [Google Scholar] [CrossRef] [PubMed]
  6. Gong, X.; Li, X.; Li, Y.; Bo, G.; Qiu, R.; Huang, Z.; Gao, S.; Wang, S. An improved model to simulate soil water and heat: A case study for drip-irrigated tomato grown in a greenhouse. Agric. Water Manag. 2023, 277, 108121. [Google Scholar] [CrossRef]
  7. Hu, J.; Gettel, G.; Fan, Z.; Lv, H.; Zhao, Y.; Yu, Y.; Wang, J.; Butterbach-Bahl, K.; Li, G.; Lin, S. Drip fertigation promotes water and nitrogen use efficiency and yield stability through improved root growth for tomatoes in plastic greenhouse production. Agric. Ecosyst. Environ. 2021, 313, 107379. [Google Scholar] [CrossRef]
  8. Shabbir, A.; Mao, H.; Ullah, I.; Buttar, N.A.; Ajmal, M.; Solangi, K.A. Improving Water Use Efficiency by Optimizing the Root Distribution Patterns under Varying Drip Emitter Density and Drought Stress for Cherry Tomato. Agronomy 2021, 11, 3. [Google Scholar] [CrossRef]
  9. Wu, J.; Zhang, R.; Gui, S. Modeling soil water movement with water uptake by roots. Plant Soil 1999, 215, 7–17. [Google Scholar] [CrossRef]
  10. Zuo, Q.; Shi, J.; Li, Y.; Zhang, R. Root length density and water uptake distributions of winter wheat under sub-irrigation. Plant Soil 2006, 285, 45–55. [Google Scholar] [CrossRef]
  11. Ning, S.; Chen, C.; Zhou, B.; Wang, Q. Evaluation of normalized root length density distribution models. Field Crops Res. 2019, 242, 107604. [Google Scholar] [CrossRef]
  12. Li, Y.; Niu, W.; Dyck, M.; Wang, J.; Zou, X. Yields and Nutritional of Greenhouse Tomato in Response to Different Soil Aeration Volume at two depths of Subsurface drip irrigation. Sci. Rep. 2016, 6, 39307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Gázquez, J.; Lopez, J.; Pérez-Parra, J.; Baeza, E.; Sáez, M.; Parra, A. Greenhouse Cooling Strategies for Mediterranean Climate Areas. Acta Hortic. 2008, 45, 425–432. [Google Scholar] [CrossRef]
  14. Tantau, H.-J.; Lange, D. Greenhouse climate control: An approach for integrated pest management. Comput. Electron. Agric. 2003, 40, 141–152. [Google Scholar] [CrossRef]
  15. Babalola, O.; Fawusi, M.O.A. Drought susceptibility of two tomato (Lycopersicum esculentum) varieties. Plant Soil 1980, 55, 205–214. [Google Scholar] [CrossRef]
  16. Sacco, A.; Greco, B.; Di Matteo, A.; De Stefano, R.; Barone, A. Evaluation of Tomato Genetic Resources for Response to Water Deficit. Am. J. Plant Sci. 2013, 04, 131–145. [Google Scholar] [CrossRef] [Green Version]
  17. Ullah, I.; Mao, H.; Rasool, G.; Gao, H.; Javed, Q.; Sarwar, A.; Khan, M. Effect of Deficit Irrigation and Reduced N Fertilization on Plant Growth, Root Morphology, and Water Use Efficiency of Tomato Grown in Soilless Culture. Agronomy 2021, 11, 228. [Google Scholar] [CrossRef]
  18. Mardani, S.; Tabatabaei, S.H.; Pessarakli, M.; Zareabyaneh, H. Physiological responses of pepper plant (Capsicum annuum L.) to drought stress. J. Plant Nutr. 2017, 40, 1453–1464. [Google Scholar] [CrossRef] [Green Version]
  19. Niu, W.-Q.; Jia, Z.-X.; Zhang, X.; Shao, H.-B. Effects of Soil Rhizosphere Aeration on the Root Growth and Water Absorption of Tomato. CLEAN—Soil Air Water 2012, 40, 1364–1371. [Google Scholar] [CrossRef]
  20. Howell, T.; Moriconi, J.I.; Zhao, X.; Hegarty, J.; Fahima, T.; Santa-Maria, G.E.; Dubcovsky, J. A wheat/rye polymorphism affects seminal root length and yield across different irrigation regimes. J. Exp. Bot. 2019, 70, 4027–4037. [Google Scholar] [CrossRef] [Green Version]
  21. Zubaidi, A.; McDonald, G.K.; Hollamby, G.J. Shoot growth, root growth and grain yield of bread and durum wheat in South Australia. Aust. J. Exp. Agric. 1999, 39, 709–720. [Google Scholar] [CrossRef]
  22. Novák, V. Water uptake of maize roots under conditions of non-limiting soil water content. Soil Technol. 1994, 7, 37–45. [Google Scholar] [CrossRef]
  23. Yang, D.; Zhang, T.; Zhang, K.; Greenwood, D.J.; Hammond, J.P.; White, P.J. An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop–soil system: Validation and application. J. Hydrol. 2009, 370, 177–184. [Google Scholar] [CrossRef] [Green Version]
  24. Zuo, Q.; Jie, F.; Zhang, R.; Meng, L. A Generalized Function of Wheat’s Root Length Density Distributions. Vadose Zone J. 2004, 3, 271–277. [Google Scholar] [CrossRef]
  25. Yadav, B.K.; Mathur, S. Modeling Soil Water Uptake by Plants Using Nonlinear Dynamic Root Density Distribution Function. J. Irrig. Drain. Eng. 2008, 134, 430–436. [Google Scholar] [CrossRef]
  26. Novák, V. Estimation of soil-water extraction patterns by roots. Agric. Water Manag. 1987, 12, 271–278. [Google Scholar] [CrossRef]
  27. Metselaar, K.; Pinheiro, E.A.R.; Lier, Q.D.J.V. Mathematical Description of Rooting Profiles of Agricultural Crops and its Effect on Transpiration Prediction by a Hydrological Model. Soil Syst. 2019, 3, 44. [Google Scholar] [CrossRef] [Green Version]
  28. Gong, X.W.; Qiu, R.J.; Sun, J.S.; Ge, J.K.; Li, Y.B.; Wang, S.S. Evapotranspiration and crop coefficient of tomato grown in a solar greenhouse under full and deficit irrigation. Agric. Water Manag. 2020, 235, 106154. [Google Scholar] [CrossRef]
  29. Chandran, U.S.S.; Tarafdar, A.; Mahesha, H.S.; Sharma, M. Temperature and Soil Moisture Stress Modulate the Host Defense Response in Chickpea during Dry Root Rot Incidence. Front. Plant Sci. 2021, 12, 653265. [Google Scholar] [CrossRef]
  30. Lahti, M.; Aphalo, P.J.; Finer, L.; Ryyppo, A.; Lehto, T.; Mannerkoski, H. Effects of soil temperature on shoot and root growth and nutrient uptake of 5-year-old Norway spruce seedlings. Tree Physiol. 2005, 25, 115–122. [Google Scholar] [CrossRef] [Green Version]
  31. Li, M.; Schmidt, J.E.; Lahue, D.G.; Lazicki, P.; Kent, A.; Machmuller, M.B.; Scow, K.M.; Gaudin, A.C.M. Impact of Irrigation Strategies on Tomato Root Distribution and Rhizosphere Processes in an Organic System. Front. Plant Sci. 2020, 11, 360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Liu, R.; Wang, Z.; Ye, H.; Li, W.; Zong, R.; Tian, X. Effects of Different Plastic Mulching Film on Soil Hydrothermal Status and Water Utilization by Spring Maize in Northwest China. Front. Mater. 2021, 8, 774833. [Google Scholar] [CrossRef]
  33. Boulard, T.; Fatnassi, H.; Roy, J.; Lagier, J.; Fargues, J.; Smits, N.; Rougier, M.; Jeannequin, B. Effect of greenhouse ventilation on humidity of inside air and in leaf boundary-layer. Agric. For. Meteorol. 2004, 125, 225–239. [Google Scholar] [CrossRef]
  34. Gong, X.; Li, X.; Qiu, R.; Bo, G.; Ping, Y.; Xin, Q.; Ge, J. Ventilation and irrigation management strategy for tomato cultivated in greenhouses. Agric. Water Manag. 2022, 273, 107908. [Google Scholar] [CrossRef]
  35. López, A.; Valera, D.L.; Molina-Aiz, F. Sonic Anemometry to Measure Natural Ventilation in Greenhouses. Sensors 2011, 11, 9820–9838. [Google Scholar] [CrossRef] [PubMed]
  36. He, X.; Wang, J.; Guo, S.; Zhang, J.; Wei, B.; Sun, J.; Shu, S. Ventilation optimization of solar greenhouse with removable back walls based on CFD. Comput. Electron. Agric. 2018, 149, 16–25. [Google Scholar] [CrossRef]
  37. Liu, A.; Qu, Z.; Nachshon, U. On the potential impact of root system size and density on salt distribution in the root zone. Agric. Water Manag. 2020, 234, 106118. [Google Scholar] [CrossRef]
  38. Machado, R.M.A.; Do Rosário, M.; Oliveira, G.; Portas, C.A.M. Tomato root distribution, yield and fruit quality under subsurface drip irrigation. Plant Soil 2003, 255, 333–341. [Google Scholar] [CrossRef] [Green Version]
  39. Hou, M.; Jin, Q.; Lu, X.; Li, J.; Zhong, H.; Gao, Y. Growth, Water Use, and Nitrate-15N Uptake of Greenhouse Tomato as Influenced by Different Irrigation Patterns, 15N Labeled Depths, and Transplant Times. Front. Plant Sci. 2017, 8, 666. [Google Scholar] [CrossRef]
  40. Wang, L.; Ning, S.; Chen, X.; Li, Y.; Guo, W.; Ben-Gal, A. Modeling tomato root water uptake influenced by soil salinity under drip irrigation with an inverse method. Agric. Water Manag. 2021, 255, 106975. [Google Scholar] [CrossRef]
  41. Ge, J.K.; Liu, Y.F.; Gong, X.W.; Liu, Z.J.; Li, Y.B.; Xu, C.D. Response of Greenhouse Crop Ecophysiology, Water Consumption and Yield to Ventilation Environment Regulation. J. Inst. Eng. Ser. A 2019, 100, 743–752. [Google Scholar] [CrossRef]
  42. Wu, Y.; Yan, S.; Fan, J.; Zhang, F.; Xiang, Y.; Zheng, J.; Guo, J. Responses of growth, fruit yield, quality and water productivity of greenhouse tomato to deficit drip irrigation. Sci. Hortic. 2021, 275, 109710. [Google Scholar] [CrossRef]
  43. Fang, Y.; Du, Y.; Wang, J.; Wu, A.; Qiao, S.; Xu, B.; Zhang, S.; Siddique, K.; Chen, Y. Moderate Drought Stress Affected Root Growth and Grain Yield in Old, Modern and Newly Released Cultivars of Winter Wheat. Front. Plant Sci. 2017, 8, 672. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Shu, L.-Z.; Liu, R.; Min, W.; Wang, Y.-S.; Hong-Mei, Y.; Zhu, P.-F.; Zhu, J.-R. Regulation of soil water threshold on tomato plant growth and fruit quality under alternate partial root-zone drip irrigation. Agric. Water Manag. 2020, 238, 106200. [Google Scholar] [CrossRef]
  45. Ning, S.; Shi, J.; Zuo, Q.; Wang, S.; Ben-Gal, A. Generalization of the root length density distribution of cotton under film mulched drip irrigation. Field Crops Res. 2015, 177, 125–136. [Google Scholar] [CrossRef]
  46. Hernández, V.; Hellín, P.; Fenoll, J.; Garrido, I.; Cava, J.; Flores, P. Increasing Yield and Quality of Tomato Cultivated under High Temperature Conditions through the Use of Elicitors. Procedia Environ. Sci. 2015, 29, 184. [Google Scholar] [CrossRef] [Green Version]
  47. Li, H.; Liu, H.; Gong, X.; Li, S.; Pang, J.; Chen, Z.; Sun, J. Optimizing irrigation and nitrogen management strategy to trade off yield, crop water productivity, nitrogen use efficiency and fruit quality of greenhouse grown tomato. Agric. Water Manag. 2020, 245, 106570. [Google Scholar] [CrossRef]
  48. Talbi, S.; Romero-Puertas, M.C.; Hernandez, A.; Terron, L.; Ferchichi, A.; Sandalio, L.M. Drought tolerance in a saharian plant oudneya Africana: Role of antioxidant defences. Environ. Exp. Bot. 2015, 111, 114–126. [Google Scholar] [CrossRef]
  49. Gong, X.; Qiu, R.; Ge, J.; Bo, G.; Ping, Y.; Xin, Q.; Wang, S. Evapotranspiration partitioning of greenhouse grown tomato using a modified Priestley–Taylor model. Agric. Water Manag. 2021, 247, 106709. [Google Scholar] [CrossRef]
  50. Gong, X.; Liu, H.; Sun, J.; Gao, Y.; Zhang, H. Comparison of Shuttleworth-Wallace model and dual crop coefficient method for estimating evapotranspiration of tomato cultivated in a solar greenhouse. Agric. Water Manag. 2019, 217, 141–153. [Google Scholar] [CrossRef]
  51. Liu, H.; Duan, A.-W.; Li, F.-S.; Sun, J.-S.; Wang, Y.-C.; Sun, C.-T. Drip Irrigation Scheduling for Tomato Grown in Solar Greenhouse Based on Pan Evaporation in North China Plain. J. Integr. Agric. 2013, 12, 520–531. [Google Scholar] [CrossRef]
  52. Zou, H.Y.; Zhang, F.C.; Wu, L.F.; Xiang, Y.Z.; Fan, J.L.; Li, Z.J.; Li, S.E. Normalized root length density distribution model for spring maize under different water and fertilizer combinations. Trans. Chin. Soc. Agric. Eng. 2018, 34, 133–142. [Google Scholar] [CrossRef]
Figure 1. Variations in air speed inside greenhouses under different ventilation treatments in 2020 and 2021.
Figure 1. Variations in air speed inside greenhouses under different ventilation treatments in 2020 and 2021.
Plants 12 01625 g001
Figure 2. Characteristics of air temperature (Ta) and relative humidity (RH) under different ventilation treatments in 2020 and 2021.
Figure 2. Characteristics of air temperature (Ta) and relative humidity (RH) under different ventilation treatments in 2020 and 2021.
Plants 12 01625 g002
Figure 3. Change curves of soil temperature (ST), soil water content (SWC) and ratio of soil water content to temperature (SWC/ST) in 2021. ASWC is the average soil water content, AST is the average soil temperature in the 0–20 cm soil layer and ASWC/ST is the average soil water-heat ratio.
Figure 3. Change curves of soil temperature (ST), soil water content (SWC) and ratio of soil water content to temperature (SWC/ST) in 2021. ASWC is the average soil water content, AST is the average soil temperature in the 0–20 cm soil layer and ASWC/ST is the average soil water-heat ratio.
Plants 12 01625 g003
Figure 4. Tomato RLD distribution under different ventilation and irrigation amounts in 2021.
Figure 4. Tomato RLD distribution under different ventilation and irrigation amounts in 2021.
Plants 12 01625 g004
Figure 5. Relationship between parameters R1, R2, R3 and R0. Relationship between parameters R1 and R0, relationship between parameters R2 and R0 and relationship between parameters R3 and R0.
Figure 5. Relationship between parameters R1, R2, R3 and R0. Relationship between parameters R1 and R0, relationship between parameters R2 and R0 and relationship between parameters R3 and R0.
Plants 12 01625 g005
Figure 6. Comparison of simulated and measured NRLD values for greenhouse tomato under different ventilation and irrigation amounts.
Figure 6. Comparison of simulated and measured NRLD values for greenhouse tomato under different ventilation and irrigation amounts.
Plants 12 01625 g006
Figure 7. Diagram of the greenhouse tomato growing pattern.
Figure 7. Diagram of the greenhouse tomato growing pattern.
Plants 12 01625 g007
Table 1. Air speed (m s−1), air temperature (Ta, °C) and relative humidity (RH, %) inside the greenhouse under different ventilation treatments in 2020 and 2021.
Table 1. Air speed (m s−1), air temperature (Ta, °C) and relative humidity (RH, %) inside the greenhouse under different ventilation treatments in 2020 and 2021.
YearTreatmentCropping SeasonsTime (9:00–17:00)
TRTRSTSTRTRSTS
2020Air speed0.085 ± 0.00 c0.159 ± 0.01 a0.113 ± 0.01 b0.086 ± 0.00 c0.148 ± 0.01 a0.110 ± 0.01 b
Ta30.40 ± 0.62 a28.55 ± 0.54 b30.22 ± 0.62 a29.76 ± 0.25 a28.03 ± 0.22 b29.59 ± 0.24 a
RH64.31 ± 1.39 a61.12 ± 1.42 b62.41 ± 1.44 ab65.40 ± 0.57 a62.44 ± 0.67 b64.56 ± 0.69 a
2021Air speed0.121 ± 0.01 b0.253 ± 0.02 a0.143 ± 0.01 b0.091 ± 0.00 c0.163 ± 0.01 a0.109 ± 0.00 b
Ta30.98 ± 0.59 a28.61 ± 0.58 b30.14 ± 0.52 ab30.23 ± 0.32 a28.09 ± 0.26 b29.52 ± 0.26 a
RH61.52 ± 1.47 a54.90 ± 1.73 b61.18 ± 1.52 a62.80 ± 0.83 a56.28 ± 0.89 b62.05 ± 1.00 a
Note: Data in the table are means ± standard errors, and different lowercase letters after peer data indicate significant differences at the level of 0.05.
Table 2. Experiential parametric simulation results of the NRLD fit function for greenhouse tomato at the late growth stage under different ventilation and irrigation amounts in 2021.
Table 2. Experiential parametric simulation results of the NRLD fit function for greenhouse tomato at the late growth stage under different ventilation and irrigation amounts in 2021.
TreatmentR0R1R2R3R2F
TRK0.9−82.33134.68−70.7812.190.98519.29 **
TRK0.5−97.16152.77−75.6812.150.962119.88 **
TRSK0.9−80.79130.72−67.7511.570.968620.23 **
TRSK0.5−86.21140.18−72.7312.280.987620.2 **
TSK0.9−79.01128.15−66.6711.450.977519.94 **
TSK0.5−80.80131.13−68.3611.740.980320.43 **
Note: ** denotes reaching a highly significant level (p < 0.01).
Table 3. Results of NRMSE statistics for simulated and measured values of NRLD.
Table 3. Results of NRMSE statistics for simulated and measured values of NRLD.
TreatmentTRK0.9TRK0.5TRSK0.9TRSK0.5TSK0.9TSK0.5
20200.1290.1750.1600.1870.2010.203
20210.1930.2730.1590.1950.2090.220
Table 4. Estimating the proportion of greenhouse tomato roots allocated in 2020.
Table 4. Estimating the proportion of greenhouse tomato roots allocated in 2020.
LocationIrrigation TreatmentVentilation TreatmentAverage Values
TRTRSTS
SimulatedMeasuredSimulatedMeasuredSimulatedMeasuredSimulatedMeasured
Zr = 1/4K0.973.9%71.8%74.3%68.4%74.7%74.5%74.3%71.6%
K0.573.4%73.8%74.1%71.2%74.3%75.1%73.9%73.4%
Average73.7%72.8%74.2%69.8%74.5%74.8%74.1%72.5%
Zr = 1/2K0.987.0%86.4%89.1%86.1%90.6%87.8%88.9%86.8%
K0.584.6%88.4%88.0%83.9%88.7%88.3%87.1%86.9%
Average85.8%87.4%88.5%85.0%89.7%88.1%88.0%86.8%
Table 6. RLD and yield under different treatments in 2020 and 2021.
Table 6. RLD and yield under different treatments in 2020 and 2021.
YearTreatmentTRK0.9TRK0.5TRSK0.9TRSK0.5TSK0.9TSK0.5
2020Yield (t ha−1)139.5 ± 10.06 a124.2 ± 2.47 b147.6 ± 4.39 a119.4 ± 7.28 b143.1 ± 4.87 a125.1 ± 6.53 b
RLD (cm cm−3)4.85 ± 0.72 a4.35 ± 0.40 ab3.82 ± 0.28 bc3.02 ± 0.57 cd3.15 ± 0.43 cd2.39 ± 0.26 d
2021Yield (t ha−1)143.7 ± 9.23 ab128.5 ± 4.65 c148.1 ± 7.49 a126.3 ± 4.76 c144.5 ± 5.58 ab130.1 ± 7.24 bc
RLD (cm cm−3)6.58 ± 1.23 a5.22 ± 0.67 ab3.79 ± 0.81 bc3.27 ± 0.79 c4.97 ± 0.58 b2.50 ± 0.14 c
Note: Data in the table are means ± standard errors, and different lowercase letters after peer data indicate significant differences at the level of 0.05.
Table 5. Predicted root distribution in 2020 by setting different evaporation coefficient (1.0, 0.8 and 0.6).
Table 5. Predicted root distribution in 2020 by setting different evaporation coefficient (1.0, 0.8 and 0.6).
Ventilation TreatmentIrrigation TreatmentSoil Depth (cm)Zr = 1/4Zr = 1/2
102030405060
TRK1.05.131.520.220.220.550.2184.77%90.32%
K0.85.141.460.170.230.600.2384.26%89.36%
K0.65.171.330.080.250.700.2883.24%87.42%
TRSK1.05.111.630.290.200.460.1785.63%91.96%
K0.85.121.570.250.210.510.1985.18%91.11%
K0.65.131.520.210.220.550.2184.72%90.23%
TSK1.05.101.720.360.190.390.1486.33%93.29%
K0.85.111.660.310.200.440.1685.85%92.38%
K0.65.131.560.240.220.520.2085.04%90.85%
Note: Table shows NRLD model simulation values for different soil layers (10, 20, 30, 40, 50, 60 cm), with reference to the ventilation of TR, TRS and TS in 2020, for model application.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ge, J.; Liu, H.; Gong, X.; Yu, Z.; Li, L.; Li, Y. Root Distribution of Tomato Cultivated in Greenhouse under Different Ventilation and Water Conditions. Plants 2023, 12, 1625. https://doi.org/10.3390/plants12081625

AMA Style

Ge J, Liu H, Gong X, Yu Z, Li L, Li Y. Root Distribution of Tomato Cultivated in Greenhouse under Different Ventilation and Water Conditions. Plants. 2023; 12(8):1625. https://doi.org/10.3390/plants12081625

Chicago/Turabian Style

Ge, Jiankun, Huanhuan Liu, Xuewen Gong, Zihui Yu, Lusheng Li, and Yanbin Li. 2023. "Root Distribution of Tomato Cultivated in Greenhouse under Different Ventilation and Water Conditions" Plants 12, no. 8: 1625. https://doi.org/10.3390/plants12081625

APA Style

Ge, J., Liu, H., Gong, X., Yu, Z., Li, L., & Li, Y. (2023). Root Distribution of Tomato Cultivated in Greenhouse under Different Ventilation and Water Conditions. Plants, 12(8), 1625. https://doi.org/10.3390/plants12081625

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop