# Effects of Micropore Group Spacing and Irrigation Amount on Soil Respiration and Yield of Tomato with Microsprinkler Irrigation under Plastic Film in Greenhouse

Authors: Mingzhi Zhang; Xiaoqun Yan; Zhenguang Lu; Qingjun Bai; Yushun Zhang; Donglin Wang; Yuangang Zhou; Yuqing Yin
Journal: Journal of Sensors
 (2021)

---

## Abstract

Microsprinkler irrigation under a plastic film in the greenhouse (MSPF) is a water-saving way which adopts the porous discharge form of a pipe under the plastic film. The effects of different micropore group spacings (L1:30 cm, L2: 50 cm) and irrigation amounts [I1: 0.7 Epan; I2: 1.0 Epan; and I3: 1.2 Epan (Epan is the diameter of 20 cm standard pan evaporation, mm)] of the MSPF on the soil respiration and yield of tomato were studied. A completely randomized trial design was used, and drip irrigation under plastic film (CK1) and microsprinkler irrigation (CK2) were also used as controls. The results showed that under the same irrigation amount, the soil respiration rate, tomato yield, and water use efficiency (WUE) of MSPF in spring and autumn are 8.09% and 6.74%, 19.39% and 4.54%, and 10.03% and 2.32% higher than those of CK1, respectively; they are significantly increased by 31.02% and 20.46%, 49.22% and 38.38%, and 58.05% and 34.66% compared with those of CK2, respectively, indicating that MSPF increased the amount of CO2 emission, but tomato yield and WUE were effectively improved, and a dynamic balance was reached among them. Compared with the 50 cm micropore group spacing, the spring and autumn tomato yields and WUE under the 30 cm micropore group spacing were significantly increased by 16.00% and 13.01% and 20.85% and 14.25%, respectively, and the micropore group spacing had no significant effect on the soil respiration rate in both root and nonroot zones. When the I increased from 0.7 Epan to 1.2 Epan, the soil respiration rate and yield in the root and nonroot zones of the spring and autumn tomatoes increased at first and then decreased, and the WUE showed a decreasing trend. The relationship of soil respiration rate between the nonroot and root zones obeys a logarithmic function, and the soil respiration rate in the nonroot zone has a quadratic curve relationship with the yield of tomato. This study can provide data support for the development of water-saving irrigation and yield increase of facility agricultural tomato and the analysis of the soil carbon cycling mechanism.


---

## Body

## 1. Introduction

With global warming, research on global carbon cycling has become the focus of the scientific community. Soil respiration is an important circulation pathway of the global carbon cycle, and it plays a key role in regulating atmospheric CO2 concentration. According to the IPCC2012 report, agricultural soil is a major contributor to carbon emissions in global carbon cycling, accounting for about 19-29% of global cumulative carbon emissions [1]. Water is one of the five elements in agricultural production, which plays an important role in the carbon cycling of agricultural soil. Rainfall is scarce in arid and semiarid areas of northwest China [2]. Irrigation water accounts for more than 60% of the total water consumption in this area, which is the main way for local crops to grow, and it is also the focus of this paper [3, 4]. The development of facility agriculture provides a strong guarantee for vegetable production in northwest China. However, the irrigation water for facility agriculture in this area mainly comes from groundwater, and the progress in the exploitation of groundwater resources aggravates the water resource crisis in this region [5–7]. Therefore, under the background of the sharp increase in atmospheric carbon emissions and the shortage of water resources [8], it is of great significance to explore the water-saving technology of facility agriculture for crop growth and soil carbon cycling.As one of the main microirrigation methods in facility agriculture, microsprinkler irrigation under plastic film (MSPF, see Figure1; Figure 1 is reproduced from Zhang et al. [9]) is used to irrigate tomato with multiple groups of small holes under the plastic film [10]. The MSPF has obtained good application effects in the facility agriculture tomato irrigation. Compared with the traditional drip irrigation in facility agriculture, under the same working pressure, the single micropore velocity of MSPF is about 15 times that of labyrinth drip irrigation, which has stronger sediment carrying capacity and anticlogging performance, and the MSPF can solve the clogging problem of some drip irrigation emitters [11, 12]. When the irrigation amount is the same, the single group flow of MSPF is about 40 times that of drip irrigation, and the irrigation duration is short, so it is easy to increase the ratio of horizontal and vertical migration distances of the soil water wetting peak. As a result, the soil wetting body per unit area of the tillage layer can be increased, which can reduce the deep transport of soil water and limit the lateral development of roots under water stress [13, 14]. The soil carbon cycle was also changed due to the increase in the soil dry-wet cycle [15, 16]. Compared with the field microsprinkler pipe, the MSPF in facility agriculture can reduce the influence of the external environment on the soil wetting body [17–19]. The MSPF can solve the spray atomization of microsprinkler irrigation, which is easy to increase air humidity. At the same time, it can solve the problems of crop diseases and insect pests caused by high temperature and humidity in the facility’s agricultural environment [20–22]. Therefore, the exploration of MSPF is of great significance for enriching the technical system of greenhouse microirrigation, expanding the scope of application of microsprinkler irrigation, saving water and increasing production of crops, and emitting greenhouse gas from the soil.Figure 1
Microsprinkler irrigation under plastic film (MSPF).As one of the main vegetables planted in facility agriculture, the tomato has rich nutrition and health care value. Tomato in the greenhouse belongs to sparse planting crops. The selection of the distance between orifices in the capillary and the irrigation amount has a direct impact on soil water and heat, ventilation, mineralization and decomposition rate of organic matter, nitrogen transformation, microbial biomass, activity, etc. [23, 24]. Therefore, the study on the emission of greenhouse gases in greenhouse tomato farmland and the formation mechanism of tomato yield is of great significance to obtain the best balance point of water saving and yield increase of greenhouse tomato and to formulate measures to reduce greenhouse gas CO2 emissions in vegetable fields of MSPF. In the practice of facility agricultural production, the capillary distance between orifice outflow and the irrigation amount are often selected according to their own experience. Studies have shown that under the same emitter flow, the smaller the emitter spacing is, the closer the horizontal wetting shape is, the shorter the confluence time of wetting peaks between emitters is, and the larger the soil surface wetting area is [14, 25]. In addition, there is a positive correlation between soil moisture and soil respiration rate, which leads to an increase of CO2 emission per unit area of soil respiration [1]. The research of Enciso et al. showed that the onion yield increased first and then decreased as the distance between the emitters decreased from 30 to 15 cm [26]. It should be noted that when the emitter spacing is small, the investment and operation cost of drip irrigation equipment will also increase. In practical application, the amount of empirical irrigation by farmers is often much higher than the actual water demand of crops, and the deep leakage of irrigation water is serious, which leads to the waste of water resources [27, 28]. Studies have shown that under the same emitter spacing, the greater the emitter flow rate is, the higher the ratio of the distance between the horizontal and vertical transports of soil water content is [13]. The larger the amount of drip irrigation is, the larger the soil wetting volume is and the lower the ratio of horizontal and vertical migration distance of soil water is. There is a positive correlation between the amount of irrigation and soil water content [14, 29, 30]. Suitable soil water content can enhance root respiration and microbial activity, thus accelerating the decomposition of soil organic matter and increasing yield, as well as increasing CO2 emission from soil respiration [31–33]. Previous studies have found that reducing the irrigation amount can reduce the soil respiration rate, and excessive irrigation will also inhibit the increase of the soil respiration rate [34, 35]. In tomato and maize under drip irrigation, the irrigation amount was positively correlated with yield and soil CO2 emission [36–38] and negatively correlated with water use efficiency [23].At present, the research on soil respiration of facility agriculture is mainly focused on drip irrigation with small flow, but there are few studies on water management of MSPF on soil respiration of facility agriculture. Furthermore, there is a lack of qualitative and quantitative descriptions of the correlation between nonroot zone soil respiration and root zone soil respiration in tomato soil respiration, and there is also a lack of quantitative analysis of the relationship between soil respiration rate and yield. Therefore, it is necessary to study the soil respiration and yield of tomato by changing the distance between capillary orifice outflow (micropore group spacing) and irrigation water quantity of microsprinkler irrigation under plastic film. In this study, the experiment of MSPF was conducted to explore the effects of different micropore group spacing and irrigation amounts on soil respiration, root zone microenvironment, and yield of greenhouse tomato. The method of regression analysis was used to make a quantitative study in order to obtain the best combination model of micropore group spacing and irrigation amount for reducing soil respiration rate (CO2 emission) and increasing yield and water use efficiency of tomato in greenhouse under MSPF in northwest China. Through greenhouse experiment and regression analysis, this paper is aimed at enriching the greenhouse tomato microirrigation technology system and providing data support for water saving and high yield of protected agricultural crops and soil carbon cycle in this area.


## 2. Materials and Methods

### 2.1. Experimental Site and Management

This study was carried out in the greenhouse of Modern Agricultural Science and Technology Convention and Exhibition Center in Xi’an, Shaanxi Province (108°52°E, 34°03°N, 435 m above sea level). The greenhouse is 85 m long and 15 m wide (see Figure2; Figure 2 is reproduced from Zhang et al. [9]). The average annual temperature is 13.3°C, and the rainfall is about 613.75 mm. The soil type is sandy loam in this area, and the mass fractions of sand, silt, and clay are 63.9%, 29.63%, and 6.47%, respectively. The average bulk density of a 1.0 m soil layer is 1.48 g/cm3, the field capacity is 27.40%, and the buried depth of the groundwater table is more than 30 m [9]. The content of organic matter, namely, total phosphorus (P), total potassium (K), total nitrogen, available nitrogen, available phosphorus, and available potassium, in the plough layer before sowing was 15.53 g/kg, 10.12 g/kg, 2.01 g/kg, 1.36 g/kg, 0.70 g/kg, 0.11 g/kg, and 0.08 g/kg, respectively. The irrigation water originated from groundwater, the pH of which was 6.8, the chemical oxygen demand (COD) was 53.2 mg/L, the anionic surfactant content was 3.2 mg/L, and the chloride content was 0.48 mg/L [9].Figure 2
Schematic diagram of the greenhouse plot layout. Note: 0: water tank; 1: the pump (WQD10-12-0.75S, People Pump, Corp., Shanghai, China); 2: filter; 3: backwater valve; 4: electromagnetic flowmeter; 5: pressure gauge; 6: capillary; 7: tomato; 8: capillary valve; 9: plastic screens.The Jingfan 401 of tomato (Jingyan Yinong seed Technology Co., Ltd., Beijing) is planted on a ridge, in which the ridge is 3.4 m long and 1.2 m wide, and the irrigated plot is shown in Figure2. The row spacing of the tomato is 50 cm and the plant spacing is 40 cm. The distance between each plot is 4 m, and a layer of styrene-butadiene-styrene block copolymer building waterproof film with a depth of 1.0 m is embedded in the middle to prevent the horizontal infiltration and movement of soil moisture and avoid the influence on other plot tests. The pipe of MSPF (Hebei Plentirain Irrigation Equipment Technology Co., Ltd., Hebei, China) adopted thin-walled oblique 3 micropores; the pipe of MSPF structure parameters is shown in Figure 3 (Figure 3 is reproduced from Zhang et al. [9]) and Table 1. The control drip irrigation under plastic film (CK1, Hebei Plentirain Irrigation Equipment Technology Co., Ltd., Hebei, China) with a thin-wall labyrinth tooth channel was selected. The geometric parameters of the channel were 54.3mm×1.1mm×0.83mm, the pipe diameter was 16 mm, the distance between drippers was 30 cm, and the dripper flow rate was 2 L/h. The control of the pipe of the microsprinkler irrigation (CK2, Hebei Plentirain Irrigation Equipment Technology Co. Ltd, Hebei, China) adopted thin-walled oblique 3 micropores in which the pipe diameter was 32 mm and the micropore diameter was 0.8 mm. The micropore group spacing was 10 cm [9].Figure 3
Schematic diagram of microporous group (inside) spacing structure parameters. Note: diameter of micropore isd=0.8mm; the internal spacing of the microporous group was I=0.4cm; the angle of micropores is equal to 68; the micropore group spacing is L.Table 1
Experimental factor and design.
No.TreatmentMicropore group spacing cmkcpIrrigation amount (mm)SpringAutumn1L1I1300.7247.12152.732L1I2301.0353.03218.193L1I3301.2423.64261.834L2I1500.7247.12152.735L2I2501.0353.03218.196L2I3501.2423.64261.837CK1301.0353.03218.198CK2101.0353.03218.19Note: L: micropore group spacing of MSPF; I: irrigation amount of MSPF; CK1: drip irrigation under plastic film; CK2: microsprinkler irrigation.The experiment of planting tomatoes in a greenhouse during spring and autumn was conducted. Spring tomato and autumn tomato were planted on March 27, 2019, and August 23, 2019; irrigation began on April 4, 2019, and August 30, 2019; irrigation was stopped on July 15, 2019, and January 17, 2019; and the tomatoes were harvested on July 25, 2019, and January 30, 2020. The field management measures of all treatments of greenhouse tomato were consistent. In order to ensure the survival rate of tomato, irrigation was done uniformly after planting with reference to the planting experience of local farmers [9].


### 2.2. Experimental Design

Factors including the micropore group spacing (see Figure3) and irrigation amount (see Figure 4 and see Table 1) were set up in this study. Among them, the micropore group spacing (L) is set to 2 levels: 30 cm (L1) and 50 cm (L2); the irrigation amount was controlled based on the cumulative evaporation from a 20 cm diameter standard pan (Epan) [39], which was realized by a control coefficient (kcp), The kcp (the crop-pan coefficient) was set to 3 levels: 0.7 (I1), 1.0 (I2), and 1.2 (I3) Epan. The CK1 and CK2 were used as the control treatment. There were a total of 8 treatments, each of which was repeated 3 times, for a total of 24 test areas.Figure 4
Irrigation records.The evaporation amount was measured at 08:00 AM for each irrigation frequency; the irrigation was initiated after the measurement. The irrigation amount (W) was calculated according to formula (1) [23, 40]. The irrigation times and amounts were recorded (see Table 1), taking kcp=1.0 as an example to draw the irrigation time and irrigation amount (see Figures 4 and 5].
(1)W=A×Epan×kcp.Figure 5
Regression analysis of nonroot zone soil respiration rate (xn), root zone soil respiration rate (xr), and yield (y) of spring tomato (a) and autumn tomato (b).
(a)(b)In the given formula,Epan represents the evaporation within the interval between 2 irrigations (mm), A is the capillary control area, and kcp is the crop pan coefficient. The kcp of control CK1 and CK2 was 1.0 [9].


## 2.1. Experimental Site and Management

This study was carried out in the greenhouse of Modern Agricultural Science and Technology Convention and Exhibition Center in Xi’an, Shaanxi Province (108°52°E, 34°03°N, 435 m above sea level). The greenhouse is 85 m long and 15 m wide (see Figure2; Figure 2 is reproduced from Zhang et al. [9]). The average annual temperature is 13.3°C, and the rainfall is about 613.75 mm. The soil type is sandy loam in this area, and the mass fractions of sand, silt, and clay are 63.9%, 29.63%, and 6.47%, respectively. The average bulk density of a 1.0 m soil layer is 1.48 g/cm3, the field capacity is 27.40%, and the buried depth of the groundwater table is more than 30 m [9]. The content of organic matter, namely, total phosphorus (P), total potassium (K), total nitrogen, available nitrogen, available phosphorus, and available potassium, in the plough layer before sowing was 15.53 g/kg, 10.12 g/kg, 2.01 g/kg, 1.36 g/kg, 0.70 g/kg, 0.11 g/kg, and 0.08 g/kg, respectively. The irrigation water originated from groundwater, the pH of which was 6.8, the chemical oxygen demand (COD) was 53.2 mg/L, the anionic surfactant content was 3.2 mg/L, and the chloride content was 0.48 mg/L [9].Figure 2
Schematic diagram of the greenhouse plot layout. Note: 0: water tank; 1: the pump (WQD10-12-0.75S, People Pump, Corp., Shanghai, China); 2: filter; 3: backwater valve; 4: electromagnetic flowmeter; 5: pressure gauge; 6: capillary; 7: tomato; 8: capillary valve; 9: plastic screens.The Jingfan 401 of tomato (Jingyan Yinong seed Technology Co., Ltd., Beijing) is planted on a ridge, in which the ridge is 3.4 m long and 1.2 m wide, and the irrigated plot is shown in Figure2. The row spacing of the tomato is 50 cm and the plant spacing is 40 cm. The distance between each plot is 4 m, and a layer of styrene-butadiene-styrene block copolymer building waterproof film with a depth of 1.0 m is embedded in the middle to prevent the horizontal infiltration and movement of soil moisture and avoid the influence on other plot tests. The pipe of MSPF (Hebei Plentirain Irrigation Equipment Technology Co., Ltd., Hebei, China) adopted thin-walled oblique 3 micropores; the pipe of MSPF structure parameters is shown in Figure 3 (Figure 3 is reproduced from Zhang et al. [9]) and Table 1. The control drip irrigation under plastic film (CK1, Hebei Plentirain Irrigation Equipment Technology Co., Ltd., Hebei, China) with a thin-wall labyrinth tooth channel was selected. The geometric parameters of the channel were 54.3mm×1.1mm×0.83mm, the pipe diameter was 16 mm, the distance between drippers was 30 cm, and the dripper flow rate was 2 L/h. The control of the pipe of the microsprinkler irrigation (CK2, Hebei Plentirain Irrigation Equipment Technology Co. Ltd, Hebei, China) adopted thin-walled oblique 3 micropores in which the pipe diameter was 32 mm and the micropore diameter was 0.8 mm. The micropore group spacing was 10 cm [9].Figure 3
Schematic diagram of microporous group (inside) spacing structure parameters. Note: diameter of micropore isd=0.8mm; the internal spacing of the microporous group was I=0.4cm; the angle of micropores is equal to 68; the micropore group spacing is L.Table 1
Experimental factor and design.
No.TreatmentMicropore group spacing cmkcpIrrigation amount (mm)SpringAutumn1L1I1300.7247.12152.732L1I2301.0353.03218.193L1I3301.2423.64261.834L2I1500.7247.12152.735L2I2501.0353.03218.196L2I3501.2423.64261.837CK1301.0353.03218.198CK2101.0353.03218.19Note: L: micropore group spacing of MSPF; I: irrigation amount of MSPF; CK1: drip irrigation under plastic film; CK2: microsprinkler irrigation.The experiment of planting tomatoes in a greenhouse during spring and autumn was conducted. Spring tomato and autumn tomato were planted on March 27, 2019, and August 23, 2019; irrigation began on April 4, 2019, and August 30, 2019; irrigation was stopped on July 15, 2019, and January 17, 2019; and the tomatoes were harvested on July 25, 2019, and January 30, 2020. The field management measures of all treatments of greenhouse tomato were consistent. In order to ensure the survival rate of tomato, irrigation was done uniformly after planting with reference to the planting experience of local farmers [9].


## 2.2. Experimental Design

Factors including the micropore group spacing (see Figure3) and irrigation amount (see Figure 4 and see Table 1) were set up in this study. Among them, the micropore group spacing (L) is set to 2 levels: 30 cm (L1) and 50 cm (L2); the irrigation amount was controlled based on the cumulative evaporation from a 20 cm diameter standard pan (Epan) [39], which was realized by a control coefficient (kcp), The kcp (the crop-pan coefficient) was set to 3 levels: 0.7 (I1), 1.0 (I2), and 1.2 (I3) Epan. The CK1 and CK2 were used as the control treatment. There were a total of 8 treatments, each of which was repeated 3 times, for a total of 24 test areas.Figure 4
Irrigation records.The evaporation amount was measured at 08:00 AM for each irrigation frequency; the irrigation was initiated after the measurement. The irrigation amount (W) was calculated according to formula (1) [23, 40]. The irrigation times and amounts were recorded (see Table 1), taking kcp=1.0 as an example to draw the irrigation time and irrigation amount (see Figures 4 and 5].
(1)W=A×Epan×kcp.Figure 5
Regression analysis of nonroot zone soil respiration rate (xn), root zone soil respiration rate (xr), and yield (y) of spring tomato (a) and autumn tomato (b).
(a)(b)In the given formula,Epan represents the evaporation within the interval between 2 irrigations (mm), A is the capillary control area, and kcp is the crop pan coefficient. The kcp of control CK1 and CK2 was 1.0 [9].


## 3. Measurements and Computational Methods

### 3.1. Soil Respiration Rate (CO2 Emission)

The soil respiration rate under different treatments was measured by an Li-8100 (LI-COR, Inc., USA) infrared gas analyzer. The PVC soil survey ring of diameter 20 cm was injected into the soil 2 days before the measurement and did not destroy the undisturbed soil on the surface. During the measurement, each measuring ring was measured twice to monitor the respiration rate of nonroot zone soil (soil without plant main roots) and root zone soil (soil containing plant main roots) (see Figure6). When measuring root zone soil, the aboveground plants were cut. The soil respiration rate was measured at 11:00 AM on the 98th day of the planting of spring tomato and 97th day of the planting of autumn tomato.Figure 6
The distribution map of soil respiration monitoring points for greenhouse tomato (unit: cm). Note: 1: tomato; 2: TRIME 1; 3: TRIME 2; 4: nonroot zone soil respiration rate PVC measurement ring; 5: root zone soil respiration rate PVC measurement ring; 6: capillary; 7: single group outflow location.


### 3.2. Soil Microenvironment

#### 3.2.1. Soil CO2 Concentration

The soil CO2 concentration in the nonroot zone and root zone of tomato was measured by the injector method (see Figure 7), and 0 (CO2-0) and 20 cm (CO2-20) samples were taken from the surface. The CO2 concentration in the sample was determined by gas chromatography, and the determination time was the same as the soil respiration rate.Figure 7
Schematic diagram of the position of sampling points (unit: cm). Note: 1: pipe; 2: tomato; 3: plastic film; 4: surface soil; 5: root zone of soil sampling point; 6: nonroot zone of soil sampling points; 7: the CO2 concentration of 20 cm measurement point; 8: temperature measurement point; 9: the CO2 concentration of 0 cm measurement point; 10: soil volume moisture content measurement point.
(a)(b)


#### 3.2.2. Soil Temperature and Soil Volume Moisture Content

The soil temperature of the nonroot zone and root zone of tomato was measured by a 15 cm geothermometer. TZS-W (Shanghai Heyi Instruments & Meters Co., Ltd., China) soil moisture meter was used to measure the soil moisture on the same day as the soil respiration rate. Firstly, the soil plane surface of the nonroot zone and root zone of 0~20 cm tomato was dug out, and then the probe was inserted into the 15 cm area away from the surface to determine the soil volume moisture content (see Figure7) in the nonroot zone and root zone of tomato.


#### 3.2.3. The Water-Filled Pore Space of Soil (WFPS)

The nonroot zone and root zone of WFPS were calculated by dividing volumetric water content by total soil porosity. Total soil porosity was calculated by measuring the bulk density of the soil according to the relationship: soil porosity 1/4 (1 soil bulk density/2.65), assuming a particle density of 2.65 mg/m3 [39].


#### 3.2.4. The Tomato Nonrhizosphere and Rhizosphere Soil Total Organic Carbon (TOC), Microbial Biomass Carbon (MBC), Humic Acid Carbon (HA-C), and Fulvic Acid Carbon (FA-C) in Greenhouse

The nonrhizosphere soil adopts the conventional soil drilling method, which is to take out the soil in Figure7 and mix it well; the rhizosphere soil was collected by the shaking soil method (randomly selecting a plant, digging out the 10~20 cm root system from the soil (see Figure 7), shaking off the soil loosely combined with the root system, and then brushing the root topsoil with a soft brush as the rhizosphere soil sample). The nonrhizosphere and rhizosphere soils were randomly sampled 3 times, and the samples were taken back to the room immediately. The plant residues were removed from the fresh soil and stored in the refrigerator at 20°C after 2 mm sieving. The contents of TOC, MBC, HA-C, and FA-C in the soil were measured within 10 days, and the sampling time was the same as the soil respiration rate. The TOC and MBC were determined by the fumigation extraction method as outlined by Chatterjee et al. [40]. The HA-C and FA-C were extracted from soil samples according to the procedure of Sun et al. [41].


### 3.3. Yield and Water Use Efficiency

Four tomatoes were randomly selected from each plot, and the yield was measured by 0.01 g precision electronic scale after ripening. After the yield per plant was obtained, the tomato yield per hectare was calculated [9].The volumetric soil water content of 0~10 cm, 10~20 cm, 20~30 cm, 30~40 cm, 40~50 cm, 50~60 cm, 60~70 cm, and 70~80 cm was measured by a soil moisture sensor (TRIME-PICO-IPH, IMKO, Inc., Ettlingen, Germany). Monitoring before and after planting, two monitoring points were selected in each district (monitoring point 1 was arranged at the outflow micropore, and monitoring point 2 was arranged at the distance m between the two groups of micropores in the vertical flow direction, in whichm=25cm, Figure 6). Crop water consumption (ETa) and crop water use efficiency (WUE) were obtained by formula (2) and formula (3) [9, 42]:
(2)ETa=I±1000×H×θt1−θt2.In the formula:ETa is the crop water consumption during the growth period, mm; I is the irrigation quota of crop growth period, mm; H is the depth of the wetting layer with plan, H=0.8m; and θt1 and θt2 are the 80 cm average soil volumetric water contents at times t1 and t2, respectively (cm3/cm3).
(3)WUE=1000∗YETa.In the formula,WUE is the crop water use efficiency, kg/m3; Y is the crop grain yield, kg/hm2.


## 3.1. Soil Respiration Rate (CO2 Emission)

The soil respiration rate under different treatments was measured by an Li-8100 (LI-COR, Inc., USA) infrared gas analyzer. The PVC soil survey ring of diameter 20 cm was injected into the soil 2 days before the measurement and did not destroy the undisturbed soil on the surface. During the measurement, each measuring ring was measured twice to monitor the respiration rate of nonroot zone soil (soil without plant main roots) and root zone soil (soil containing plant main roots) (see Figure6). When measuring root zone soil, the aboveground plants were cut. The soil respiration rate was measured at 11:00 AM on the 98th day of the planting of spring tomato and 97th day of the planting of autumn tomato.Figure 6
The distribution map of soil respiration monitoring points for greenhouse tomato (unit: cm). Note: 1: tomato; 2: TRIME 1; 3: TRIME 2; 4: nonroot zone soil respiration rate PVC measurement ring; 5: root zone soil respiration rate PVC measurement ring; 6: capillary; 7: single group outflow location.


## 3.2. Soil Microenvironment

### 3.2.1. Soil CO2 Concentration

The soil CO2 concentration in the nonroot zone and root zone of tomato was measured by the injector method (see Figure 7), and 0 (CO2-0) and 20 cm (CO2-20) samples were taken from the surface. The CO2 concentration in the sample was determined by gas chromatography, and the determination time was the same as the soil respiration rate.Figure 7
Schematic diagram of the position of sampling points (unit: cm). Note: 1: pipe; 2: tomato; 3: plastic film; 4: surface soil; 5: root zone of soil sampling point; 6: nonroot zone of soil sampling points; 7: the CO2 concentration of 20 cm measurement point; 8: temperature measurement point; 9: the CO2 concentration of 0 cm measurement point; 10: soil volume moisture content measurement point.
(a)(b)


### 3.2.2. Soil Temperature and Soil Volume Moisture Content

The soil temperature of the nonroot zone and root zone of tomato was measured by a 15 cm geothermometer. TZS-W (Shanghai Heyi Instruments & Meters Co., Ltd., China) soil moisture meter was used to measure the soil moisture on the same day as the soil respiration rate. Firstly, the soil plane surface of the nonroot zone and root zone of 0~20 cm tomato was dug out, and then the probe was inserted into the 15 cm area away from the surface to determine the soil volume moisture content (see Figure7) in the nonroot zone and root zone of tomato.


### 3.2.3. The Water-Filled Pore Space of Soil (WFPS)

The nonroot zone and root zone of WFPS were calculated by dividing volumetric water content by total soil porosity. Total soil porosity was calculated by measuring the bulk density of the soil according to the relationship: soil porosity 1/4 (1 soil bulk density/2.65), assuming a particle density of 2.65 mg/m3 [39].


### 3.2.4. The Tomato Nonrhizosphere and Rhizosphere Soil Total Organic Carbon (TOC), Microbial Biomass Carbon (MBC), Humic Acid Carbon (HA-C), and Fulvic Acid Carbon (FA-C) in Greenhouse

The nonrhizosphere soil adopts the conventional soil drilling method, which is to take out the soil in Figure7 and mix it well; the rhizosphere soil was collected by the shaking soil method (randomly selecting a plant, digging out the 10~20 cm root system from the soil (see Figure 7), shaking off the soil loosely combined with the root system, and then brushing the root topsoil with a soft brush as the rhizosphere soil sample). The nonrhizosphere and rhizosphere soils were randomly sampled 3 times, and the samples were taken back to the room immediately. The plant residues were removed from the fresh soil and stored in the refrigerator at 20°C after 2 mm sieving. The contents of TOC, MBC, HA-C, and FA-C in the soil were measured within 10 days, and the sampling time was the same as the soil respiration rate. The TOC and MBC were determined by the fumigation extraction method as outlined by Chatterjee et al. [40]. The HA-C and FA-C were extracted from soil samples according to the procedure of Sun et al. [41].


## 3.2.1. Soil CO2 Concentration

The soil CO2 concentration in the nonroot zone and root zone of tomato was measured by the injector method (see Figure 7), and 0 (CO2-0) and 20 cm (CO2-20) samples were taken from the surface. The CO2 concentration in the sample was determined by gas chromatography, and the determination time was the same as the soil respiration rate.Figure 7
Schematic diagram of the position of sampling points (unit: cm). Note: 1: pipe; 2: tomato; 3: plastic film; 4: surface soil; 5: root zone of soil sampling point; 6: nonroot zone of soil sampling points; 7: the CO2 concentration of 20 cm measurement point; 8: temperature measurement point; 9: the CO2 concentration of 0 cm measurement point; 10: soil volume moisture content measurement point.
(a)(b)


## 3.2.2. Soil Temperature and Soil Volume Moisture Content

The soil temperature of the nonroot zone and root zone of tomato was measured by a 15 cm geothermometer. TZS-W (Shanghai Heyi Instruments & Meters Co., Ltd., China) soil moisture meter was used to measure the soil moisture on the same day as the soil respiration rate. Firstly, the soil plane surface of the nonroot zone and root zone of 0~20 cm tomato was dug out, and then the probe was inserted into the 15 cm area away from the surface to determine the soil volume moisture content (see Figure7) in the nonroot zone and root zone of tomato.


## 3.2.3. The Water-Filled Pore Space of Soil (WFPS)

The nonroot zone and root zone of WFPS were calculated by dividing volumetric water content by total soil porosity. Total soil porosity was calculated by measuring the bulk density of the soil according to the relationship: soil porosity 1/4 (1 soil bulk density/2.65), assuming a particle density of 2.65 mg/m3 [39].


## 3.2.4. The Tomato Nonrhizosphere and Rhizosphere Soil Total Organic Carbon (TOC), Microbial Biomass Carbon (MBC), Humic Acid Carbon (HA-C), and Fulvic Acid Carbon (FA-C) in Greenhouse

The nonrhizosphere soil adopts the conventional soil drilling method, which is to take out the soil in Figure7 and mix it well; the rhizosphere soil was collected by the shaking soil method (randomly selecting a plant, digging out the 10~20 cm root system from the soil (see Figure 7), shaking off the soil loosely combined with the root system, and then brushing the root topsoil with a soft brush as the rhizosphere soil sample). The nonrhizosphere and rhizosphere soils were randomly sampled 3 times, and the samples were taken back to the room immediately. The plant residues were removed from the fresh soil and stored in the refrigerator at 20°C after 2 mm sieving. The contents of TOC, MBC, HA-C, and FA-C in the soil were measured within 10 days, and the sampling time was the same as the soil respiration rate. The TOC and MBC were determined by the fumigation extraction method as outlined by Chatterjee et al. [40]. The HA-C and FA-C were extracted from soil samples according to the procedure of Sun et al. [41].


## 3.3. Yield and Water Use Efficiency

Four tomatoes were randomly selected from each plot, and the yield was measured by 0.01 g precision electronic scale after ripening. After the yield per plant was obtained, the tomato yield per hectare was calculated [9].The volumetric soil water content of 0~10 cm, 10~20 cm, 20~30 cm, 30~40 cm, 40~50 cm, 50~60 cm, 60~70 cm, and 70~80 cm was measured by a soil moisture sensor (TRIME-PICO-IPH, IMKO, Inc., Ettlingen, Germany). Monitoring before and after planting, two monitoring points were selected in each district (monitoring point 1 was arranged at the outflow micropore, and monitoring point 2 was arranged at the distance m between the two groups of micropores in the vertical flow direction, in whichm=25cm, Figure 6). Crop water consumption (ETa) and crop water use efficiency (WUE) were obtained by formula (2) and formula (3) [9, 42]:
(2)ETa=I±1000×H×θt1−θt2.In the formula:ETa is the crop water consumption during the growth period, mm; I is the irrigation quota of crop growth period, mm; H is the depth of the wetting layer with plan, H=0.8m; and θt1 and θt2 are the 80 cm average soil volumetric water contents at times t1 and t2, respectively (cm3/cm3).
(3)WUE=1000∗YETa.In the formula,WUE is the crop water use efficiency, kg/m3; Y is the crop grain yield, kg/hm2.


## 4. Data Analysis

The significant difference of SPSS22.0 (IBM Crop., Armonk, New York, NY, USA) was analyzed by theF test, and the significance level was set to P<0.05. OriginPro2019 (OriginLab Corporation, Northampton, MA, USA) was used to draw the picture. Mathematica12.0 (Wolfram Research, New York, NY, USA) was used for regression analysis.


## 5. Results

### 5.1. Effects of Different Treatments on Soil Water and Heat of Tomato in Greenhouse

The micropore group spacing (L) had no significant effect on soil water and heat in the nonroot zone and root zone of tomato (P>0.05, see Table 2). The irrigation amount (I) also exhibited a significant effect on soil water and heat in the nonrhizosphere and root zone of tomato (P≤0.05).Table 2
Effects of different treatments on soil moisture and temperature in nonroot zone and root zone of tomato in greenhouse.
TreatmentSpringAutumnTemperature (°C)Soil volume moisture content (%)Temperature (°C)Soil volume moisture content (%)Nonroot zone soilRoot zone soilNonroot zone soilRoot zone soilNonroot zone soilRoot zone soilNonroot zone soilRoot zone soilL1I126.1±0.71abc25.07±0.87ab20.55±5.97c19.56±3.18c20.14±0.52a18.69±0.5ab22.04±3.92bc20.55±3.44cdL1I226.99±0.22a25.57±1.18a27.73±7.56ab27.48±3.84ab20.15±0.93a19.61±2.75a28.22±3.98a28.22±3.39abL1I325.32±1.13bc23.31±2.49bc31.44±3.31a29.21±5.12a17.57±1.66c16.91±1.48bc28.72±4.67a28.72±2.77aL2I126.37±2.34abc24.93±1.8ab20.06±4.28c18.32±3.34c19.67±1.53ab19.25±0.61a19.81±3.83ab18.32±2.92dL2I226.77±0.81ab24.84±1.97ab27.73±4.06ab26.91±3.35ab19.34±1.14ab18.49±1.07ab25.75±2.24c25.62±3.37abL2I325.12±1.36c22.23±1.37c30.2±3.07ab26.99±3.69ab17.41±0.98c16.37±1.98c28.47±2.18a28.22±3.51abCK125.43±1.13abc25.12±1.66ab25.5±1.74bc27.32±2.35ab18.94±0.6abc18.09±1.01abc27.23±5.12a27.23±3.47abCK225.09±1.14c23.55±0.93abc25.24±2.28bc23.93±2.25b18.41±1.91bc17.08±0.31bc25.81±2.02ab24.16±2.01bcF valueL0.014ns1.311ns0.091ns1.124ns1.456ns0.464ns1.588ns2.697nsI5.168∗7.567∗∗13.453∗∗21.070∗∗15.497∗∗8.708∗∗17.142∗∗26.607∗∗L∗I0.145ns0.230ns0.043ns0.143ns0.223ns0.851ns0.173ns0.362nsNote: L: micropore group spacing; I: irrigation amount. The data are allaverage±standarddeviation in the chart; different letters in the same column mean significant difference at 0.05 level, ∗P<0.05; ∗∗P<0.01; ns: P>0.05, the same as below.The nonroot zone soil temperature and root zone soil temperature of spring and autumn tomatoes treated with L1I2 increased by 6.16% and 6.38% and 1.78% and 8.36%, respectively, compared with those of CK1 and 7.58% and 9.45% and 8.58% and 14.78% higher than those of CK2, respectively. The volume water content of the nonroot zone and root zone of tomato treated with L1I2 increased by 8.75% and 3.64% and 0.57% and 3.64%, compared with that of CK1, respectively, and increased by 9.85% and 9.34% and 14.83% and 16.80% compared with that of CK2, respectively. With an increase in irrigation amount, the volume water content of the nonroot zone and root zone of tomato first increased, while the soil nonroot zone temperature and root zone temperature of tomato decreased. The soil nonroot zone temperature, root zone temperature, nonroot zone volume water content, and root zone volume water content of tomato with the micropore group spacing of 30 cm were slightly higher than those of 50 cm by about 0.19% and 2.55%, 2.71% and 2.02%, 2.22% and 6.69%, and 5.59% and 7.38%, respectively.


### 5.2. Effects of Different Treatments on the Soil CO2 Concentration in Nonroot Zone and Root Zone of Tomato

The micropore group spacing (L) had no significant effect on CO2 concentration in the nonroot zone and root zone of tomato soil at different depths (see Table 3). The irrigation amount (I) also exhibited a significant effect on CO2 concentration in the nonroot zone and root zone of tomato soil at different depths. The interaction between them had no significant effect on the concentration of CO2 in the nonroot zone and root zone of soil at different depths.Table 3
Effects of different treatments on soil CO2 concentration of greenhouse tomato at different depths.
TreatmentSpringAutumnThe CO2 concentration of 0 cm (μmol/mol)The CO2 concentration of 20 cm (μmol/mol)The CO2 concentration of 0 cm (μmol/mol)The CO2 concentration of 20 cm (μmol/mol)Nonroot zone soilRoot zone soilNonroot zone soilRoot zone soilNonroot zone soilRoot zone soilNonroot zone soilRoot zone soilL1I1430.93±22.38b448.28±37.33ab11549.36±3088.94b15387.65±3353.14b354.37±17.44b362.06±33.02b9911.93±1542.99b12220.01±2087.97bL1I2460.74±14.86a475.31±32.63a13909.31±2817.28ab18369.31±3474.77ab409.53±25.36a457.31±53.35ab11910.05±2638.63ab14832.84±3354.07abL1I3445.51±23.51ab468.43±19.07ab15238.13±4111.09a19529.57±3560.48a406.69±42.87a424.605±45.46ab13106.25±3051.99a15800.78±2645.28aL2I1429.97±29.42b433.85±32.33b11116.69±1199.53b15209.09±3858.32b347.53±27.89b356.79±47.27b9894.33±2011.61b12087.22±2232.93bL2I2450.31±19.7ab476.66±24.27a13466.39±4681.6ab17815.02±5512.68ab402.3±32.05a435.43±32.39ab11663.94±3104.59ab14547.37±3466.27abL2I3446.74±23.26ab474.06±30.25a14979.31±2576.25a19353.99±2662.52a402.22±47.14a421.88±63.02ab12829.23±2418.95a15614.11±3259.88aCK1448.66±13.2ab473.11±17.92a13366.68±1382.84ab18065.81±2218.28ab402.34±18.45a451.73±44.67ab12069.23±2366.66ab14276.91±3646.29abCK2432.44±5.09b441.15±17.82ab12527.66±2109.63ab16997.02±2780.95ab390.18±12.52a406.75±27.51ab11653.99±2744.25ab12261.35±1285.67bF valueL0.202ns0.062ns0.180ns0.084ns0.303ns0.405ns0.069ns0.066nsI3.766∗4.810∗6.101∗∗5.458∗10.369∗∗11.013∗∗6.756∗∗7.202∗∗L∗I0.225ns0.374ns0.004ns0.014ns0.006ns0.146ns0.014ns0.003nsThe soil CO2 concentration of 0 cm of tomato was significantly lower than the depth of the 20 cm soil layer, indicating that as the depth of the soil layer increased, the soil CO2 concentration increased. With the seasonal change, the soil CO2 concentration of autumn tomato was lower than that of spring tomato. The CO2 concentration in the nonroot zone of the 0 cm soil layer and the CO2 concentration in the root zone of 0 cm soil treated with L1I2 were 2.69% and 1.79% and 0.47%and 1.24% higher than those of CK1 and 6.55% and 4.96% and 7.74% and 12.43% higher than those of CK2, respectively. The CO2 concentration in the nonroot zone of the 20 cm soil layer and the CO2 concentration in the root zone of the 20 cm soil treated with L1I2 were 4.06% and -1.32% and 1.68% and 3.89% higher than those of CK1 and 11.03% and 2.20% and 8.07% and 20.97% higher than those of CK2, respectively.With the increase of irrigation amount, the soil CO2 concentration in the nonroot zone and root zone of tomato at different depths increased first and then decreased. The CO2 concentration in the nonroot zone of 0 cm soil, the CO2 concentration in the root zone of 0 cm soil, the CO2 concentration in the nonroot zone of 20 cm soil, and the CO2 concentration in the root zone of 20 cm soil were higher than those of 50 cm in tomato soil with a micropore group spacing of 30 cm.


### 5.3. Effects of Different Treatments on Soil Respiration Rate of Greenhouse Tomato

Figure8 shows that the micropore group spacing had no significant effect on the respiration rate of the nonroot zone and root zone of tomato soil. The irrigation amount had a significant effect on the respiration rate of the nonroot zone and the root zone of tomato soil (P≤0.05). The soil nonroot zone and root zone respiration rates of autumn tomato were about 2.81% and 4.14% lower than those of spring tomato, respectively. The soil respiration rate in the root zone of spring tomato and autumn tomato was about 14.66% and 16.25% higher than that in the nonroot zone, respectively.Figure 8
Effects of different treatments on soil respiration rate of tomato in greenhouse. Note: the data are allaverage±standarddeviation in the figure; different letters in the same column meant significant difference at 0.05 level, the same as below.Compared with CK1, the soil nonroot zone respiration rate and root zone respiration rate of L1I2 spring tomato and an autumn tomato increased by 11.39% and 6.22% and 4.80% and 7.26%, respectively, and increased by 27.08% and 39.70% and 42.25% and 37.06%, respectively, compared with CK2. With the increase of irrigation amount, the respiration rate of tomato soil nonroot zone and root zone increased first and then showed a slight decreasing trend. The soil respiration rate of 30 cm with the micropore group spacing is higher than 50 cm.


### 5.4. Effects of Different Treatments on the Related Carbon Content of Soil in Greenhouse Tomato

It can be seen from Table4 that the micropore group spacing has no significant effect on TOC, MBC, HA-C, and FA-C in the nonrhizosphere and rhizosphere of soil with spring tomato and autumn tomato; irrigation amount has a significant effect on TOC, MBC, HA-C, and FA-C in the nonrhizosphere and rhizosphere of soil with spring tomato and autumn tomato. The TOC, MBC, HA-C, and FA-C of the rhizosphere soil of spring tomato and autumn tomato were 27.94% and 26.97%, 19.47% and 25.23%, 19.47% and 22.20%, and 8.15% and 2.78% higher than those of the nonrhizosphere soil, respectively.Table 4
Effects of different treatments on soil related carbon content of greenhouse tomato.
TreatmentF valueL1I1L1I2L1I3L2I1L2I2L2I3CK1CK2LIL∗ISpringTOC (%)Nonrhizosphere0.96±0.12c1.18±0.11a1.16±0.09a0.95±0.17c1.17±0.07a1.16±0.09a1.09±0.04ab1.04±0.07bc0.032ns21.533∗∗0.016nsRhizosphere1.31±0.09b1.46±0.14a1.45±0.11a1.3±0.17ba1.46±0.09a1.46±0.1a1.37±0.05ab1.33±0.09b0.001ns9.964∗∗0.029nsMBC (mg/kg)Nonrhizosphere108.27±17.47bc121.53±8.38a121.31±10.36a107.35±13.48c120.78±6.35a121.55±11.58a108.17±4.7ab115.34±7.53abc0.022ns7.804∗∗0.012nsRhizosphere128.11±14.12c147.64±10.26a145.27±15.19a125.88±10.7c142.4±12.98ab145.56±15.5a144.21±8.09ab132.55±8.21bc0.439ns11.295∗∗0.196nsHA-C (g/kg)Nonrhizosphere1.63±0.29b2.07±0.14a2.04±0.24a1.64±0.5b2.05±0.24a2.07±0.36a2.03±0.3a1.89±0.12ab0.001ns10.806∗∗0.021nsRhizosphere2.07±0.6b2.47±0.18a2.39±0.26ab2.08±0.43b2.4±0.23ab2.37±0.15ab2.39±0.31ab2.1±0.37b0.073ns5.721∗0.060nsFA-C (g/kg)Nonrhizosphere2.15±0.45bc2.66±0.31a2.59±0.4a2.06±0.34c2.57±0.37a2.57±0.34a2.65±0.38a2.48±0.37ab0.386ns10.462∗∗0.053nsRhizosphere2.46±0.27ab2.82±0.26a2.72±0.46ab2.36±0.43b2.77±0.39a2.76±0.37a2.79±0.26a2.65±0.42ab0.156ns5.760∗0.162nsAutumnTOC (%)Nonrhizosphere0.92±0.07b1.07±0.05a1.07±0.07a0.9±0.07b1.07±0.05a1.05±0.05a1.03±0.07a1.01±0.11a0.657ns39.2936∗∗0.225nsRhizosphere1.18±0.08b1.34±0.06a1.35±0.07a1.16±0.09b1.34±0.07a1.35±0.07a1.32±0.08a1.22±0.11b0.034ns33.320∗∗0.045nsMBC (mg/kg)Nonrhizosphere92.31±6.01b107.48±7.51a107.14±8.85a94.5±17.2b107.52±10.74a106.96±7.42a101.37±11.89ab97.1±7.65ab0.060ns10.881∗∗0.073nsRhizosphere119.39±7.6c134.61±9.5a134.31±11.2a118.83±20.62c133.68±13.96a134.08±9.39a132.74±16.58ab120.94±6.66bc0.027ns8.320∗∗0.003nsHA-C (g/kg)Nonrhizosphere1.7±0.31bc2.05±0.12a1.94±0.26ab1.63±0.31c1.95±0.19ab1.96±0.26ab1.93±0.3ab1.76±0.36abc0.530ns9.415∗∗0.240nsRhizosphere2.13±0.39c2.49±0.16a2.46±0.21a2.15±0.27bc2.29±0.18abc2.36±0.27abc2.41±0.18ab2.32±0.29abc1.6767ns6.337∗∗0.871nsFA-C (g/kg)Nonrhizosphere2.14±0.21b2.58±0.46a2.57±0.39a2.12±0.36b2.52±0.44ab2.55±0.39a2.51±0.43ab2.32±0.39ab0.090ns7.410∗∗0.017nsRhizosphere2.3±0.25ab2.65±0.4a2.57±0.38ab2.21±0.26b2.59±0.4ab2.59±0.36ab2.58±0.42ab2.35±0.3ab0.188ns6.041∗0.096nsNote: TOC: soil total organic carbon; MBC: soil microbial biomass carbon; HA-C: soil humic acid carbon; FA-C: soil fulvic acid carbon.In terms of nonrhizosphere of soil with spring tomato and autumn tomato, the TOC, MBC, HA-C, and FA-C of L1I2 treatments were higher than those of CK1 by about 8.24% and 4.44%, 2.36% and 6.03%, 2.09%and 6.23%, and 0.36% and 2.72%, respectively, and also higher than those of CK2 by about 13.19% and 5.59%, 7.22% and 10.69%, 9.72% and 16.15%, and 6.89% and 11.24%, respectively. For the rhizosphere soil of spring tomato and autumn tomato, the TOC, MBC, HA-C, and FA-C of L1I2 treatments were higher than those of CK1 by about 6.64% and 1.60%, 2.38% and 1.14%, 3.37% and 3.49%, and 1.32% and 2.57%, respectively, and also higher than those of CK2 by about 10.12% and 10.31%, 11.39% and 11.30%, 17.55% and 7.55%, and 6.39% and 12.86%, respectively. With the increase of irrigation amount, the TOC, MBC, HA-C, and FA-C of nonrhizosphere and rhizosphere soils of spring tomato and autumn tomato increased at first and then decreased. The TOC, MBC, HA-C, and FA-C of nonrhizosphere and rhizosphere soils of spring tomato and autumn tomato in the micropore group spacing of 30 cm was higher than 50 cm.


### 5.5. Effects of Different Treatments on Tomato Yield and Water Use Efficiency

As can be seen from Figure9, the micropore group spacing had a significant effect on the yield and water use efficiency of spring and autumn tomatoes (P≤0.05). The irrigation amount had a significant effect on the yield and water use efficiency of spring and autumn tomatoes.Figure 9
Effects of different treatments on tomato yield and water use efficiency.Compared with CK1, the yield and water use efficiency of spring tomato and autumn tomato treated with L1I2 increased by 19.39% and 4.54% and 10.03% and 2.32%, respectively. Compared with CK2, the yield and water use efficiency of spring tomato and autumn tomato treated with L1I2 increased by 20.46% and 49.22% and 31.02% and 58.05%, respectively. The yield of spring tomato and autumn tomato treated with L1I2 was about 31.88% and 28.03%, 0.96% and 1.43%, 47.73% and 44.49%, 24.74% and 21.09%, and 12.21% and 7.48% higher than that of L1I1, L1I3, L2I1, L2I2, and L2I3, respectively. With the increase of irrigation amount, the yield of spring tomato and autumn tomato increased, while the water use efficiency of spring tomato and autumn tomato decreased. The yield and water use efficiency of spring tomato and autumn tomato with micropore group spacing 30 cm were higher than those of 50 cm by about 16.00% and 13.01% and 20.85% and 14.25%, respectively.


### 5.6. Regression Analysis of Soil Respiration Rate and Yield of Microsprinkler Irrigation under Plastic Film

As can be seen from Figure5, there is a logarithmic relationship between xn and xr (formula (4)) and a conic relationship between xn and y (formula (5)).
(4)xr=2.6749lnxn+0.6565,(5)y=−0.425.90xn2+41378xn+11143.Through the simultaneous establishment of the above two equations, the relationship betweenxr and y can be obtained, that is, y=2602.13e0.7478xr+32370.01e0.3739xr+11143. In the regression model (4), the determination coefficient of xn and xr regression is R2=0.8237, indicating that the degree to which xn can explain xr in this model is up to 82%, and the coefficient xr can be predicted by xn. In the regression model (5), the determination coefficient of xn and y regression is R2=6992, indicating that the degree to which xn can explain y in this model is up to 69%, and the coefficient y can be predicted by xn.


## 5.1. Effects of Different Treatments on Soil Water and Heat of Tomato in Greenhouse

The micropore group spacing (L) had no significant effect on soil water and heat in the nonroot zone and root zone of tomato (P>0.05, see Table 2). The irrigation amount (I) also exhibited a significant effect on soil water and heat in the nonrhizosphere and root zone of tomato (P≤0.05).Table 2
Effects of different treatments on soil moisture and temperature in nonroot zone and root zone of tomato in greenhouse.
TreatmentSpringAutumnTemperature (°C)Soil volume moisture content (%)Temperature (°C)Soil volume moisture content (%)Nonroot zone soilRoot zone soilNonroot zone soilRoot zone soilNonroot zone soilRoot zone soilNonroot zone soilRoot zone soilL1I126.1±0.71abc25.07±0.87ab20.55±5.97c19.56±3.18c20.14±0.52a18.69±0.5ab22.04±3.92bc20.55±3.44cdL1I226.99±0.22a25.57±1.18a27.73±7.56ab27.48±3.84ab20.15±0.93a19.61±2.75a28.22±3.98a28.22±3.39abL1I325.32±1.13bc23.31±2.49bc31.44±3.31a29.21±5.12a17.57±1.66c16.91±1.48bc28.72±4.67a28.72±2.77aL2I126.37±2.34abc24.93±1.8ab20.06±4.28c18.32±3.34c19.67±1.53ab19.25±0.61a19.81±3.83ab18.32±2.92dL2I226.77±0.81ab24.84±1.97ab27.73±4.06ab26.91±3.35ab19.34±1.14ab18.49±1.07ab25.75±2.24c25.62±3.37abL2I325.12±1.36c22.23±1.37c30.2±3.07ab26.99±3.69ab17.41±0.98c16.37±1.98c28.47±2.18a28.22±3.51abCK125.43±1.13abc25.12±1.66ab25.5±1.74bc27.32±2.35ab18.94±0.6abc18.09±1.01abc27.23±5.12a27.23±3.47abCK225.09±1.14c23.55±0.93abc25.24±2.28bc23.93±2.25b18.41±1.91bc17.08±0.31bc25.81±2.02ab24.16±2.01bcF valueL0.014ns1.311ns0.091ns1.124ns1.456ns0.464ns1.588ns2.697nsI5.168∗7.567∗∗13.453∗∗21.070∗∗15.497∗∗8.708∗∗17.142∗∗26.607∗∗L∗I0.145ns0.230ns0.043ns0.143ns0.223ns0.851ns0.173ns0.362nsNote: L: micropore group spacing; I: irrigation amount. The data are allaverage±standarddeviation in the chart; different letters in the same column mean significant difference at 0.05 level, ∗P<0.05; ∗∗P<0.01; ns: P>0.05, the same as below.The nonroot zone soil temperature and root zone soil temperature of spring and autumn tomatoes treated with L1I2 increased by 6.16% and 6.38% and 1.78% and 8.36%, respectively, compared with those of CK1 and 7.58% and 9.45% and 8.58% and 14.78% higher than those of CK2, respectively. The volume water content of the nonroot zone and root zone of tomato treated with L1I2 increased by 8.75% and 3.64% and 0.57% and 3.64%, compared with that of CK1, respectively, and increased by 9.85% and 9.34% and 14.83% and 16.80% compared with that of CK2, respectively. With an increase in irrigation amount, the volume water content of the nonroot zone and root zone of tomato first increased, while the soil nonroot zone temperature and root zone temperature of tomato decreased. The soil nonroot zone temperature, root zone temperature, nonroot zone volume water content, and root zone volume water content of tomato with the micropore group spacing of 30 cm were slightly higher than those of 50 cm by about 0.19% and 2.55%, 2.71% and 2.02%, 2.22% and 6.69%, and 5.59% and 7.38%, respectively.


## 5.2. Effects of Different Treatments on the Soil CO2 Concentration in Nonroot Zone and Root Zone of Tomato

The micropore group spacing (L) had no significant effect on CO2 concentration in the nonroot zone and root zone of tomato soil at different depths (see Table 3). The irrigation amount (I) also exhibited a significant effect on CO2 concentration in the nonroot zone and root zone of tomato soil at different depths. The interaction between them had no significant effect on the concentration of CO2 in the nonroot zone and root zone of soil at different depths.Table 3
Effects of different treatments on soil CO2 concentration of greenhouse tomato at different depths.
TreatmentSpringAutumnThe CO2 concentration of 0 cm (μmol/mol)The CO2 concentration of 20 cm (μmol/mol)The CO2 concentration of 0 cm (μmol/mol)The CO2 concentration of 20 cm (μmol/mol)Nonroot zone soilRoot zone soilNonroot zone soilRoot zone soilNonroot zone soilRoot zone soilNonroot zone soilRoot zone soilL1I1430.93±22.38b448.28±37.33ab11549.36±3088.94b15387.65±3353.14b354.37±17.44b362.06±33.02b9911.93±1542.99b12220.01±2087.97bL1I2460.74±14.86a475.31±32.63a13909.31±2817.28ab18369.31±3474.77ab409.53±25.36a457.31±53.35ab11910.05±2638.63ab14832.84±3354.07abL1I3445.51±23.51ab468.43±19.07ab15238.13±4111.09a19529.57±3560.48a406.69±42.87a424.605±45.46ab13106.25±3051.99a15800.78±2645.28aL2I1429.97±29.42b433.85±32.33b11116.69±1199.53b15209.09±3858.32b347.53±27.89b356.79±47.27b9894.33±2011.61b12087.22±2232.93bL2I2450.31±19.7ab476.66±24.27a13466.39±4681.6ab17815.02±5512.68ab402.3±32.05a435.43±32.39ab11663.94±3104.59ab14547.37±3466.27abL2I3446.74±23.26ab474.06±30.25a14979.31±2576.25a19353.99±2662.52a402.22±47.14a421.88±63.02ab12829.23±2418.95a15614.11±3259.88aCK1448.66±13.2ab473.11±17.92a13366.68±1382.84ab18065.81±2218.28ab402.34±18.45a451.73±44.67ab12069.23±2366.66ab14276.91±3646.29abCK2432.44±5.09b441.15±17.82ab12527.66±2109.63ab16997.02±2780.95ab390.18±12.52a406.75±27.51ab11653.99±2744.25ab12261.35±1285.67bF valueL0.202ns0.062ns0.180ns0.084ns0.303ns0.405ns0.069ns0.066nsI3.766∗4.810∗6.101∗∗5.458∗10.369∗∗11.013∗∗6.756∗∗7.202∗∗L∗I0.225ns0.374ns0.004ns0.014ns0.006ns0.146ns0.014ns0.003nsThe soil CO2 concentration of 0 cm of tomato was significantly lower than the depth of the 20 cm soil layer, indicating that as the depth of the soil layer increased, the soil CO2 concentration increased. With the seasonal change, the soil CO2 concentration of autumn tomato was lower than that of spring tomato. The CO2 concentration in the nonroot zone of the 0 cm soil layer and the CO2 concentration in the root zone of 0 cm soil treated with L1I2 were 2.69% and 1.79% and 0.47%and 1.24% higher than those of CK1 and 6.55% and 4.96% and 7.74% and 12.43% higher than those of CK2, respectively. The CO2 concentration in the nonroot zone of the 20 cm soil layer and the CO2 concentration in the root zone of the 20 cm soil treated with L1I2 were 4.06% and -1.32% and 1.68% and 3.89% higher than those of CK1 and 11.03% and 2.20% and 8.07% and 20.97% higher than those of CK2, respectively.With the increase of irrigation amount, the soil CO2 concentration in the nonroot zone and root zone of tomato at different depths increased first and then decreased. The CO2 concentration in the nonroot zone of 0 cm soil, the CO2 concentration in the root zone of 0 cm soil, the CO2 concentration in the nonroot zone of 20 cm soil, and the CO2 concentration in the root zone of 20 cm soil were higher than those of 50 cm in tomato soil with a micropore group spacing of 30 cm.


## 5.3. Effects of Different Treatments on Soil Respiration Rate of Greenhouse Tomato

Figure8 shows that the micropore group spacing had no significant effect on the respiration rate of the nonroot zone and root zone of tomato soil. The irrigation amount had a significant effect on the respiration rate of the nonroot zone and the root zone of tomato soil (P≤0.05). The soil nonroot zone and root zone respiration rates of autumn tomato were about 2.81% and 4.14% lower than those of spring tomato, respectively. The soil respiration rate in the root zone of spring tomato and autumn tomato was about 14.66% and 16.25% higher than that in the nonroot zone, respectively.Figure 8
Effects of different treatments on soil respiration rate of tomato in greenhouse. Note: the data are allaverage±standarddeviation in the figure; different letters in the same column meant significant difference at 0.05 level, the same as below.Compared with CK1, the soil nonroot zone respiration rate and root zone respiration rate of L1I2 spring tomato and an autumn tomato increased by 11.39% and 6.22% and 4.80% and 7.26%, respectively, and increased by 27.08% and 39.70% and 42.25% and 37.06%, respectively, compared with CK2. With the increase of irrigation amount, the respiration rate of tomato soil nonroot zone and root zone increased first and then showed a slight decreasing trend. The soil respiration rate of 30 cm with the micropore group spacing is higher than 50 cm.


## 5.4. Effects of Different Treatments on the Related Carbon Content of Soil in Greenhouse Tomato

It can be seen from Table4 that the micropore group spacing has no significant effect on TOC, MBC, HA-C, and FA-C in the nonrhizosphere and rhizosphere of soil with spring tomato and autumn tomato; irrigation amount has a significant effect on TOC, MBC, HA-C, and FA-C in the nonrhizosphere and rhizosphere of soil with spring tomato and autumn tomato. The TOC, MBC, HA-C, and FA-C of the rhizosphere soil of spring tomato and autumn tomato were 27.94% and 26.97%, 19.47% and 25.23%, 19.47% and 22.20%, and 8.15% and 2.78% higher than those of the nonrhizosphere soil, respectively.Table 4
Effects of different treatments on soil related carbon content of greenhouse tomato.
TreatmentF valueL1I1L1I2L1I3L2I1L2I2L2I3CK1CK2LIL∗ISpringTOC (%)Nonrhizosphere0.96±0.12c1.18±0.11a1.16±0.09a0.95±0.17c1.17±0.07a1.16±0.09a1.09±0.04ab1.04±0.07bc0.032ns21.533∗∗0.016nsRhizosphere1.31±0.09b1.46±0.14a1.45±0.11a1.3±0.17ba1.46±0.09a1.46±0.1a1.37±0.05ab1.33±0.09b0.001ns9.964∗∗0.029nsMBC (mg/kg)Nonrhizosphere108.27±17.47bc121.53±8.38a121.31±10.36a107.35±13.48c120.78±6.35a121.55±11.58a108.17±4.7ab115.34±7.53abc0.022ns7.804∗∗0.012nsRhizosphere128.11±14.12c147.64±10.26a145.27±15.19a125.88±10.7c142.4±12.98ab145.56±15.5a144.21±8.09ab132.55±8.21bc0.439ns11.295∗∗0.196nsHA-C (g/kg)Nonrhizosphere1.63±0.29b2.07±0.14a2.04±0.24a1.64±0.5b2.05±0.24a2.07±0.36a2.03±0.3a1.89±0.12ab0.001ns10.806∗∗0.021nsRhizosphere2.07±0.6b2.47±0.18a2.39±0.26ab2.08±0.43b2.4±0.23ab2.37±0.15ab2.39±0.31ab2.1±0.37b0.073ns5.721∗0.060nsFA-C (g/kg)Nonrhizosphere2.15±0.45bc2.66±0.31a2.59±0.4a2.06±0.34c2.57±0.37a2.57±0.34a2.65±0.38a2.48±0.37ab0.386ns10.462∗∗0.053nsRhizosphere2.46±0.27ab2.82±0.26a2.72±0.46ab2.36±0.43b2.77±0.39a2.76±0.37a2.79±0.26a2.65±0.42ab0.156ns5.760∗0.162nsAutumnTOC (%)Nonrhizosphere0.92±0.07b1.07±0.05a1.07±0.07a0.9±0.07b1.07±0.05a1.05±0.05a1.03±0.07a1.01±0.11a0.657ns39.2936∗∗0.225nsRhizosphere1.18±0.08b1.34±0.06a1.35±0.07a1.16±0.09b1.34±0.07a1.35±0.07a1.32±0.08a1.22±0.11b0.034ns33.320∗∗0.045nsMBC (mg/kg)Nonrhizosphere92.31±6.01b107.48±7.51a107.14±8.85a94.5±17.2b107.52±10.74a106.96±7.42a101.37±11.89ab97.1±7.65ab0.060ns10.881∗∗0.073nsRhizosphere119.39±7.6c134.61±9.5a134.31±11.2a118.83±20.62c133.68±13.96a134.08±9.39a132.74±16.58ab120.94±6.66bc0.027ns8.320∗∗0.003nsHA-C (g/kg)Nonrhizosphere1.7±0.31bc2.05±0.12a1.94±0.26ab1.63±0.31c1.95±0.19ab1.96±0.26ab1.93±0.3ab1.76±0.36abc0.530ns9.415∗∗0.240nsRhizosphere2.13±0.39c2.49±0.16a2.46±0.21a2.15±0.27bc2.29±0.18abc2.36±0.27abc2.41±0.18ab2.32±0.29abc1.6767ns6.337∗∗0.871nsFA-C (g/kg)Nonrhizosphere2.14±0.21b2.58±0.46a2.57±0.39a2.12±0.36b2.52±0.44ab2.55±0.39a2.51±0.43ab2.32±0.39ab0.090ns7.410∗∗0.017nsRhizosphere2.3±0.25ab2.65±0.4a2.57±0.38ab2.21±0.26b2.59±0.4ab2.59±0.36ab2.58±0.42ab2.35±0.3ab0.188ns6.041∗0.096nsNote: TOC: soil total organic carbon; MBC: soil microbial biomass carbon; HA-C: soil humic acid carbon; FA-C: soil fulvic acid carbon.In terms of nonrhizosphere of soil with spring tomato and autumn tomato, the TOC, MBC, HA-C, and FA-C of L1I2 treatments were higher than those of CK1 by about 8.24% and 4.44%, 2.36% and 6.03%, 2.09%and 6.23%, and 0.36% and 2.72%, respectively, and also higher than those of CK2 by about 13.19% and 5.59%, 7.22% and 10.69%, 9.72% and 16.15%, and 6.89% and 11.24%, respectively. For the rhizosphere soil of spring tomato and autumn tomato, the TOC, MBC, HA-C, and FA-C of L1I2 treatments were higher than those of CK1 by about 6.64% and 1.60%, 2.38% and 1.14%, 3.37% and 3.49%, and 1.32% and 2.57%, respectively, and also higher than those of CK2 by about 10.12% and 10.31%, 11.39% and 11.30%, 17.55% and 7.55%, and 6.39% and 12.86%, respectively. With the increase of irrigation amount, the TOC, MBC, HA-C, and FA-C of nonrhizosphere and rhizosphere soils of spring tomato and autumn tomato increased at first and then decreased. The TOC, MBC, HA-C, and FA-C of nonrhizosphere and rhizosphere soils of spring tomato and autumn tomato in the micropore group spacing of 30 cm was higher than 50 cm.


## 5.5. Effects of Different Treatments on Tomato Yield and Water Use Efficiency

As can be seen from Figure9, the micropore group spacing had a significant effect on the yield and water use efficiency of spring and autumn tomatoes (P≤0.05). The irrigation amount had a significant effect on the yield and water use efficiency of spring and autumn tomatoes.Figure 9
Effects of different treatments on tomato yield and water use efficiency.Compared with CK1, the yield and water use efficiency of spring tomato and autumn tomato treated with L1I2 increased by 19.39% and 4.54% and 10.03% and 2.32%, respectively. Compared with CK2, the yield and water use efficiency of spring tomato and autumn tomato treated with L1I2 increased by 20.46% and 49.22% and 31.02% and 58.05%, respectively. The yield of spring tomato and autumn tomato treated with L1I2 was about 31.88% and 28.03%, 0.96% and 1.43%, 47.73% and 44.49%, 24.74% and 21.09%, and 12.21% and 7.48% higher than that of L1I1, L1I3, L2I1, L2I2, and L2I3, respectively. With the increase of irrigation amount, the yield of spring tomato and autumn tomato increased, while the water use efficiency of spring tomato and autumn tomato decreased. The yield and water use efficiency of spring tomato and autumn tomato with micropore group spacing 30 cm were higher than those of 50 cm by about 16.00% and 13.01% and 20.85% and 14.25%, respectively.


## 5.6. Regression Analysis of Soil Respiration Rate and Yield of Microsprinkler Irrigation under Plastic Film

As can be seen from Figure5, there is a logarithmic relationship between xn and xr (formula (4)) and a conic relationship between xn and y (formula (5)).
(4)xr=2.6749lnxn+0.6565,(5)y=−0.425.90xn2+41378xn+11143.Through the simultaneous establishment of the above two equations, the relationship betweenxr and y can be obtained, that is, y=2602.13e0.7478xr+32370.01e0.3739xr+11143. In the regression model (4), the determination coefficient of xn and xr regression is R2=0.8237, indicating that the degree to which xn can explain xr in this model is up to 82%, and the coefficient xr can be predicted by xn. In the regression model (5), the determination coefficient of xn and y regression is R2=6992, indicating that the degree to which xn can explain y in this model is up to 69%, and the coefficient y can be predicted by xn.


## 6. Discussion

### 6.1. Effects of Irrigation Methods on Soil Microenvironment and Soil Respiration Rate

Previous studies have shown that different irrigation methods lead to different forms of water entering the soil. It can change the shape of the soil wetting body, the size of the dry and wet zones, and the number of cycles and affect the mineralization decomposition rate of soil organic matter, microbial biomass, and soil enzyme activity and then increase or decrease the concentration of CO2 in soil. Under the concentration gradient drive, the diffusion rate of soil CO2 to the surface increases or slows down [38, 43]. There is a positive correlation between soil respiration rate and soil moisture in arid and semiarid areas, in which the limitation of soil moisture is mainly reflected in reducing the availability of soil organic matter and the death of soil microorganisms caused by low water potential [44, 45]. In this study, the soil respiration rate of tomato with MSPF was higher than that of drip irrigation under plastic film (see Figure 8). It is possible that under the same working pressure and irrigation amount, the effluent velocity of drip irrigation under plastic film is about 0.067 times that of MSPF, and the irrigation duration increases, so that the ratio of horizontal to vertical migration distance of soil water in drip irrigation under plastic film is smaller than that of MSPF, and the higher horizontal movement of soil water is easy to increase the wetting volume of the tillage layer (0~40 cm) per unit area [46]. The soil moisture of tomato with drip irrigation under a plastic film is mainly concentrated under the emitter, and there is an obvious dry-wet zone and a long dry-wet duration, which restricts the activities of microorganisms and enzymes in the dry area, and the soil water content is too high in the wet area, which reduces the soil gas exchange. As a result, the soil respiration rate of drip irrigation under plastic film is lower than that of MSPF. Burger et al. found that the peak of CO2 emission from tomato soil respiration occurred when WFPS was about 60% [47]. It is mainly due to the fact that soil moisture is the main factor limiting soil respiration when WFPS<60%. Too high or low soil moisture will reduce the stability of soil aggregates. The destruction of soil aggregates is usually accompanied by the increase of the utilization rate of organic carbon sources, which can increase the activities of microorganisms and enzymes and increase the mineralization of soil carbon [48, 49]. At the same time, the increase of soil moisture increased microbial activity and plant physiological activities; increased soil TOC, MBC, HA-C, and FA-C; and further supplied energy sources for microbial activities and increased soil microbial respiration. When WFPS>60%, higher soil moisture occupies a large number of soil pores, reducing soil oxygen content, and limiting soil gas diffusion rate [39]. In this study, the WFPS of MSPF is close to 60%, and the moisture and temperature of MSPF are higher than 5.77% and 4.30% and 3.99% and 13.47% (see Table 5) of drip irrigation under plastic film, respectively, which may also be one of the reasons for the higher soil respiration of MSPF.Table 5
Effects of different treatments on soil water-filled porosity of tomato.
TreatmentSpringAutumnNonroot zone soil (%)Root zone soil (%)Nonroot zone soil (%)Root zone soil (%)L1I155.69±10.04ab54.41±10.89ab54.25±9.78ab49.85±12.37bL1I262.55±9.26ab61.76±5.81a60.52±8.52ab58.05±7.4abL1I363.2±5.35a62.29±8.05a63.49±10.09a63.09±8.88aL2I153.72±8.72b50.82±7.17b50.47±8.79b47.85±13.4bL2I256.77±6.99ab57.17±6.54ab55.39±12.51ab54.34±8.03abL2I363.92±9.08a58.56±8.79ab63.55±9.63a62.34±8.91aCK162.63±6.31ab61.84±4.79a60.13±7.72ab57.82±13.7abCK257.68±11.19ab55.18±9.51ab56.82±7.96ab54.69±7.44abF valueL1.052ns3.289ns1.181ns0.617nsI5.039∗5.042∗5.635∗∗8.516∗∗L∗I0.682ns0.021ns0.328ns0.097nsIt was found that the soil respiration rate of the tomato root zone was higher than that of the nonroot zone. It may be due to the fact that the soil in the crop root zone contains plant roots. Previous studies have shown that the diversity, abundance, and activity of microorganisms in the plant rhizosphere soil are high [50], and the respiration of microorganisms and plant roots is easy to increase soil respiration rate (CO2 emission) [51, 52]. In this study, it was also found that the soil respiration rate of spring tomato was higher than that of autumn tomato. Mainly because the soil temperature of spring tomato is 32.19% higher than that of autumn tomato (L1I2, see Table 2), lower temperature can reduce cell activity, which is not conducive to microorganisms and crop respiration [53]. Previous studies have found that the effect of irrigation on the soil respiration rate is attributed to the increase of crop yield to stimulate root autotrophic respiration and the change of soil water content to stimulate soil heterotrophic respiration [54]. In this study, the yield of tomato and soil moisture of MSPF were higher than those of drip irrigation under plastic film by 19.39 and 4.54% (see Figure 9) and 5.77 and 4.30% (see Table 2), respectively, indicating that MSPF can increase soil respiration rate and lead to the increase of the soil CO2 emission rate. The increase of the soil respiration rate in greenhouse agriculture contributes to the formation of greenhouse tomato yield, which can achieve a new dynamic balance among tomato crop water use efficiency, yield, and CO2 emission from greenhouse vegetable fields.This study found that the humidity of spring tomato and autumn tomato canopy under MSPF was about 36.03% and 25.20% lower than that of microsprinkler irrigation (see Figure10), respectively. This may be due to the large number of micropores per unit length of microsprinkler irrigation, so the water spray is easy to atomize, and the increase of surface wetting area also easily increases ineffective water transpiration, which leads to a significant increase in air humidity in the microenvironment of tomato growth, which is consistent with the finding of Man et al. [15] who explored the mechanism of action of microsprinkler irrigation. The results also showed that the soil respiration rate of tomato under MSPF was higher than that of microsprinkler irrigation, and the yield of tomato under MSPF was 4.69% and 49.22% higher than that of microsprinkler irrigation (see Figure 9), which was consistent with the conclusion of Scheer et al. [54]. Mainly due to the low oxygen content of the air humidity, it will reduce the soil aerobic respiration microbial activity, and its byproduct CO2 emissions will also be reduced. The increase in crop yield will further stimulate the autotrophic respiration of roots, resulting in a higher soil respiration rate of MSPF. This study also found that the soil WFPS of 0~10 cm under MSPF was about 60%, which was higher than 10.15% and 6.33% of microsprinkler irrigation. At the same time, the soil moisture volume content of 20 cm under MSPF was 7.96% and 7.22% (see Table 2) higher than that of microsprinkler irrigation, indicating that soil moisture was also one of the main reasons affecting the soil respiration rate of tomato under MSPF.Figure 10
Effects of different treatments on air humidity in tomato canopy with spring and autumn.


### 6.2. Effects of Micropore Group Spacing and Irrigation Amount on Soil Respiration Rate

In this study, the micropore group spacing has little effect on the soil respiration rate, which may be due to the large flow rate and short irrigation duration of the single group of MSPF. When the flow of WFPS is greater than the soil infiltration rate, it is easy to produce local microrunoff, which increases the soil water horizontal transport distance and reduces the effect of relatively small micropore group spacing. There is no significant difference in 0~10 cm soil WFPS between the two groups (P>0.05). It was also found that the micropore group spacing had no significant effect on soil moisture and temperature. Previous studies have found that water and heat changes directly affect soil microorganisms and enzyme activities, thus changing soil TOC, MBC, HA-C, and FA-C cycles [38, 55]. There was no significant difference in soil TOC, MBC, HA-C, and FA-C among different micropore group spacings, which further indicated that the micropore group spacing had little effect on soil respiration rate.Previous studies have found that properly reducing the irrigation amount can reduce the soil respiration rate. With the increase of irrigation amount, soil water stimulates root respiration and soil microbial activity, resulting in an increase in the soil respiration rate [34, 35]. However, when the irrigation amount increased to a certain extent, crop root aerobic respiration and soil microbial activity were gradually inhibited and further inhibited soil respiration rate [56, 57]. This may also be one of the reasons why the soil respiration rate increased at first and then decreased with the increase of irrigation amount in this study. Chen et al.’s [39] study shows that the soil respiration rate increases with the increase of the irrigation amount, which is inconsistent with the conclusion of this paper. It may be caused by the difference in irrigation amount. The highest irrigation amount in Chen et al.’s study was 1.0 Epan, while the highest in this study was 1.2 Epan. Agbna et al. [58] and Zhu et al. [59] found that when the irrigation amount of drip irrigation with tomato increased from 0.5 Epan to 1.0 Epan, the yield increased while the WUE decreased. It shows that the effect of MSPF on tomato yield and water use efficiency was similar to drip irrigation, and the MSPF was suitable for greenhouse tomato irrigation.


### 6.3. Correlation between Soil Microenvironment, Soil Respiration Rate, and Yield of Tomato with Microsprinkler Irrigation under Plastic Film

Previous studies have shown that when the protected agricultural space is relatively closed and the irrigation water is limited, reducing the irrigation amount and increasing greenhouse CO2 concentration can improve the yield and quality of greenhouse tomato [60, 61]. In this study, it was found that there was a significant positive correlation between the soil respiration rate and yield (see Figure 11), which was consistent with the conclusion that Mancinelli et al. [62] found that there is a positive correlation between the CO2 emission rate of tomato soil and yield under drip irrigation. It may be because the greenhouse microenvironment is relatively closed and the CO2 content is limited. For tomatoes with a high photorespiration rate, increasing atmospheric CO2 concentration in a certain range can not only provide more sufficient raw materials for photosynthesis but also increase the activity of (RuBP) carboxylase of ribulose 5-diphosphate, which is beneficial to accelerate the binding of RuBP in chloroplasts to CO2 entering chloroplasts, thus enhancing the ability of photosynthesis to fix CO2, formatting 3-phosphoglyceric acid (PGA), and further synthesizing photosynthetic carbohydrate through C3 cycling [63, 64]. In this study, there was a positive correlation between soil moisture and soil respiration rate (see Figure 11), which was consistent with that obtained by Wei et al. [65] when soil moisture content in maize planting under drip irrigation was lower than the field water holding capacity. This study also found that there was a synergistic change of TOC and MBC in the soil of tomato under MSPF, which was positively correlated with the soil respiration rate (see Figure 11). It is consistent with the finding of Han et al. [66] to study the changes of soil TOC and MBC under drip irrigation, which shows that the effect of MSPF on soil TOC and MBC is similar to that of drip irrigation under plastic film.Figure 11
Correlation between soil microenvironment, soil respiration rate, and yield of tomato. Note: WFPS: water-filled pore space of soil; CO2-20: soil CO2 concentration of 20 cm; CO2-0: soil CO2 concentration of 0 cm; M: soil moisture; T: soil temperature; NSCF: soil respiration rate in nonroot zone; RSCF: soil respiration rate in root zone. Different color depth means correlation difference; correlation value is showed by the value of the number from -0.1 to 1. ∗ means significance correlation (P<0.05); ∗∗ means significance correlation (P<0.01), the same below.Previous studies have found that there is a significant positive correlation between the soil temperature and soil respiration rate [34, 65, 67]. However, there was no significant difference in the positive effect of temperature on the soil respiration rate in this study. It may be that the soil water distribution of MSPF is relatively uniform, the soil moisture in the root zone and the nonroot zone of tomato is not too high or too low, and the soil microenvironment is relatively stable, which lead to the decrease of the effect of temperature on the soil respiration rate. Wang et al. [68] believe that when the WFPS of the soil is less than 60%, there is a positive correlation between the soil respiration rate and soil WFPS; otherwise, there is a negative correlation between them. In this study, there was a positive correlation between soil WFPS and the soil respiration rate of greenhouse tomato, which is inconsistent with the conclusion that there was a negative correlation between the soil respiration rate and soil WFPS under drip irrigation [69]. It is possibly due to the differences in soil types, irrigation methods, and the irrigation amount controlled by the experiment. Chen used drip irrigation of tomato in Yangling soil (brown loam), and the irrigation amount increased from 0.6 Epan to 1.0 Epan. However, in this study, tomato was irrigated with MSPF on sandy loam (yellow cinnamon soil), and the irrigation amount increased from 0.7 Epan to 1.2 Epan.


## 6.1. Effects of Irrigation Methods on Soil Microenvironment and Soil Respiration Rate

Previous studies have shown that different irrigation methods lead to different forms of water entering the soil. It can change the shape of the soil wetting body, the size of the dry and wet zones, and the number of cycles and affect the mineralization decomposition rate of soil organic matter, microbial biomass, and soil enzyme activity and then increase or decrease the concentration of CO2 in soil. Under the concentration gradient drive, the diffusion rate of soil CO2 to the surface increases or slows down [38, 43]. There is a positive correlation between soil respiration rate and soil moisture in arid and semiarid areas, in which the limitation of soil moisture is mainly reflected in reducing the availability of soil organic matter and the death of soil microorganisms caused by low water potential [44, 45]. In this study, the soil respiration rate of tomato with MSPF was higher than that of drip irrigation under plastic film (see Figure 8). It is possible that under the same working pressure and irrigation amount, the effluent velocity of drip irrigation under plastic film is about 0.067 times that of MSPF, and the irrigation duration increases, so that the ratio of horizontal to vertical migration distance of soil water in drip irrigation under plastic film is smaller than that of MSPF, and the higher horizontal movement of soil water is easy to increase the wetting volume of the tillage layer (0~40 cm) per unit area [46]. The soil moisture of tomato with drip irrigation under a plastic film is mainly concentrated under the emitter, and there is an obvious dry-wet zone and a long dry-wet duration, which restricts the activities of microorganisms and enzymes in the dry area, and the soil water content is too high in the wet area, which reduces the soil gas exchange. As a result, the soil respiration rate of drip irrigation under plastic film is lower than that of MSPF. Burger et al. found that the peak of CO2 emission from tomato soil respiration occurred when WFPS was about 60% [47]. It is mainly due to the fact that soil moisture is the main factor limiting soil respiration when WFPS<60%. Too high or low soil moisture will reduce the stability of soil aggregates. The destruction of soil aggregates is usually accompanied by the increase of the utilization rate of organic carbon sources, which can increase the activities of microorganisms and enzymes and increase the mineralization of soil carbon [48, 49]. At the same time, the increase of soil moisture increased microbial activity and plant physiological activities; increased soil TOC, MBC, HA-C, and FA-C; and further supplied energy sources for microbial activities and increased soil microbial respiration. When WFPS>60%, higher soil moisture occupies a large number of soil pores, reducing soil oxygen content, and limiting soil gas diffusion rate [39]. In this study, the WFPS of MSPF is close to 60%, and the moisture and temperature of MSPF are higher than 5.77% and 4.30% and 3.99% and 13.47% (see Table 5) of drip irrigation under plastic film, respectively, which may also be one of the reasons for the higher soil respiration of MSPF.Table 5
Effects of different treatments on soil water-filled porosity of tomato.
TreatmentSpringAutumnNonroot zone soil (%)Root zone soil (%)Nonroot zone soil (%)Root zone soil (%)L1I155.69±10.04ab54.41±10.89ab54.25±9.78ab49.85±12.37bL1I262.55±9.26ab61.76±5.81a60.52±8.52ab58.05±7.4abL1I363.2±5.35a62.29±8.05a63.49±10.09a63.09±8.88aL2I153.72±8.72b50.82±7.17b50.47±8.79b47.85±13.4bL2I256.77±6.99ab57.17±6.54ab55.39±12.51ab54.34±8.03abL2I363.92±9.08a58.56±8.79ab63.55±9.63a62.34±8.91aCK162.63±6.31ab61.84±4.79a60.13±7.72ab57.82±13.7abCK257.68±11.19ab55.18±9.51ab56.82±7.96ab54.69±7.44abF valueL1.052ns3.289ns1.181ns0.617nsI5.039∗5.042∗5.635∗∗8.516∗∗L∗I0.682ns0.021ns0.328ns0.097nsIt was found that the soil respiration rate of the tomato root zone was higher than that of the nonroot zone. It may be due to the fact that the soil in the crop root zone contains plant roots. Previous studies have shown that the diversity, abundance, and activity of microorganisms in the plant rhizosphere soil are high [50], and the respiration of microorganisms and plant roots is easy to increase soil respiration rate (CO2 emission) [51, 52]. In this study, it was also found that the soil respiration rate of spring tomato was higher than that of autumn tomato. Mainly because the soil temperature of spring tomato is 32.19% higher than that of autumn tomato (L1I2, see Table 2), lower temperature can reduce cell activity, which is not conducive to microorganisms and crop respiration [53]. Previous studies have found that the effect of irrigation on the soil respiration rate is attributed to the increase of crop yield to stimulate root autotrophic respiration and the change of soil water content to stimulate soil heterotrophic respiration [54]. In this study, the yield of tomato and soil moisture of MSPF were higher than those of drip irrigation under plastic film by 19.39 and 4.54% (see Figure 9) and 5.77 and 4.30% (see Table 2), respectively, indicating that MSPF can increase soil respiration rate and lead to the increase of the soil CO2 emission rate. The increase of the soil respiration rate in greenhouse agriculture contributes to the formation of greenhouse tomato yield, which can achieve a new dynamic balance among tomato crop water use efficiency, yield, and CO2 emission from greenhouse vegetable fields.This study found that the humidity of spring tomato and autumn tomato canopy under MSPF was about 36.03% and 25.20% lower than that of microsprinkler irrigation (see Figure10), respectively. This may be due to the large number of micropores per unit length of microsprinkler irrigation, so the water spray is easy to atomize, and the increase of surface wetting area also easily increases ineffective water transpiration, which leads to a significant increase in air humidity in the microenvironment of tomato growth, which is consistent with the finding of Man et al. [15] who explored the mechanism of action of microsprinkler irrigation. The results also showed that the soil respiration rate of tomato under MSPF was higher than that of microsprinkler irrigation, and the yield of tomato under MSPF was 4.69% and 49.22% higher than that of microsprinkler irrigation (see Figure 9), which was consistent with the conclusion of Scheer et al. [54]. Mainly due to the low oxygen content of the air humidity, it will reduce the soil aerobic respiration microbial activity, and its byproduct CO2 emissions will also be reduced. The increase in crop yield will further stimulate the autotrophic respiration of roots, resulting in a higher soil respiration rate of MSPF. This study also found that the soil WFPS of 0~10 cm under MSPF was about 60%, which was higher than 10.15% and 6.33% of microsprinkler irrigation. At the same time, the soil moisture volume content of 20 cm under MSPF was 7.96% and 7.22% (see Table 2) higher than that of microsprinkler irrigation, indicating that soil moisture was also one of the main reasons affecting the soil respiration rate of tomato under MSPF.Figure 10
Effects of different treatments on air humidity in tomato canopy with spring and autumn.


## 6.2. Effects of Micropore Group Spacing and Irrigation Amount on Soil Respiration Rate

In this study, the micropore group spacing has little effect on the soil respiration rate, which may be due to the large flow rate and short irrigation duration of the single group of MSPF. When the flow of WFPS is greater than the soil infiltration rate, it is easy to produce local microrunoff, which increases the soil water horizontal transport distance and reduces the effect of relatively small micropore group spacing. There is no significant difference in 0~10 cm soil WFPS between the two groups (P>0.05). It was also found that the micropore group spacing had no significant effect on soil moisture and temperature. Previous studies have found that water and heat changes directly affect soil microorganisms and enzyme activities, thus changing soil TOC, MBC, HA-C, and FA-C cycles [38, 55]. There was no significant difference in soil TOC, MBC, HA-C, and FA-C among different micropore group spacings, which further indicated that the micropore group spacing had little effect on soil respiration rate.Previous studies have found that properly reducing the irrigation amount can reduce the soil respiration rate. With the increase of irrigation amount, soil water stimulates root respiration and soil microbial activity, resulting in an increase in the soil respiration rate [34, 35]. However, when the irrigation amount increased to a certain extent, crop root aerobic respiration and soil microbial activity were gradually inhibited and further inhibited soil respiration rate [56, 57]. This may also be one of the reasons why the soil respiration rate increased at first and then decreased with the increase of irrigation amount in this study. Chen et al.’s [39] study shows that the soil respiration rate increases with the increase of the irrigation amount, which is inconsistent with the conclusion of this paper. It may be caused by the difference in irrigation amount. The highest irrigation amount in Chen et al.’s study was 1.0 Epan, while the highest in this study was 1.2 Epan. Agbna et al. [58] and Zhu et al. [59] found that when the irrigation amount of drip irrigation with tomato increased from 0.5 Epan to 1.0 Epan, the yield increased while the WUE decreased. It shows that the effect of MSPF on tomato yield and water use efficiency was similar to drip irrigation, and the MSPF was suitable for greenhouse tomato irrigation.


## 6.3. Correlation between Soil Microenvironment, Soil Respiration Rate, and Yield of Tomato with Microsprinkler Irrigation under Plastic Film

Previous studies have shown that when the protected agricultural space is relatively closed and the irrigation water is limited, reducing the irrigation amount and increasing greenhouse CO2 concentration can improve the yield and quality of greenhouse tomato [60, 61]. In this study, it was found that there was a significant positive correlation between the soil respiration rate and yield (see Figure 11), which was consistent with the conclusion that Mancinelli et al. [62] found that there is a positive correlation between the CO2 emission rate of tomato soil and yield under drip irrigation. It may be because the greenhouse microenvironment is relatively closed and the CO2 content is limited. For tomatoes with a high photorespiration rate, increasing atmospheric CO2 concentration in a certain range can not only provide more sufficient raw materials for photosynthesis but also increase the activity of (RuBP) carboxylase of ribulose 5-diphosphate, which is beneficial to accelerate the binding of RuBP in chloroplasts to CO2 entering chloroplasts, thus enhancing the ability of photosynthesis to fix CO2, formatting 3-phosphoglyceric acid (PGA), and further synthesizing photosynthetic carbohydrate through C3 cycling [63, 64]. In this study, there was a positive correlation between soil moisture and soil respiration rate (see Figure 11), which was consistent with that obtained by Wei et al. [65] when soil moisture content in maize planting under drip irrigation was lower than the field water holding capacity. This study also found that there was a synergistic change of TOC and MBC in the soil of tomato under MSPF, which was positively correlated with the soil respiration rate (see Figure 11). It is consistent with the finding of Han et al. [66] to study the changes of soil TOC and MBC under drip irrigation, which shows that the effect of MSPF on soil TOC and MBC is similar to that of drip irrigation under plastic film.Figure 11
Correlation between soil microenvironment, soil respiration rate, and yield of tomato. Note: WFPS: water-filled pore space of soil; CO2-20: soil CO2 concentration of 20 cm; CO2-0: soil CO2 concentration of 0 cm; M: soil moisture; T: soil temperature; NSCF: soil respiration rate in nonroot zone; RSCF: soil respiration rate in root zone. Different color depth means correlation difference; correlation value is showed by the value of the number from -0.1 to 1. ∗ means significance correlation (P<0.05); ∗∗ means significance correlation (P<0.01), the same below.Previous studies have found that there is a significant positive correlation between the soil temperature and soil respiration rate [34, 65, 67]. However, there was no significant difference in the positive effect of temperature on the soil respiration rate in this study. It may be that the soil water distribution of MSPF is relatively uniform, the soil moisture in the root zone and the nonroot zone of tomato is not too high or too low, and the soil microenvironment is relatively stable, which lead to the decrease of the effect of temperature on the soil respiration rate. Wang et al. [68] believe that when the WFPS of the soil is less than 60%, there is a positive correlation between the soil respiration rate and soil WFPS; otherwise, there is a negative correlation between them. In this study, there was a positive correlation between soil WFPS and the soil respiration rate of greenhouse tomato, which is inconsistent with the conclusion that there was a negative correlation between the soil respiration rate and soil WFPS under drip irrigation [69]. It is possibly due to the differences in soil types, irrigation methods, and the irrigation amount controlled by the experiment. Chen used drip irrigation of tomato in Yangling soil (brown loam), and the irrigation amount increased from 0.6 Epan to 1.0 Epan. However, in this study, tomato was irrigated with MSPF on sandy loam (yellow cinnamon soil), and the irrigation amount increased from 0.7 Epan to 1.2 Epan.


## 7. Conclusion

In this paper, by exploring the response mechanism of greenhouse soil microenvironment, soil respiration rate, and yield of tomato under different micropore group spacing and irrigation amount of MSPF, the suitable micropore group spacing and irrigation amount of MSPF for greenhouse tomato growth were obtained. The results showed that compared with drip irrigation under a plastic film, and the yield of spring tomato and autumn tomato increased by 19.39% and 4.54%, respectively. Compared with microsprinkler irrigation, tomato yield, water use efficiency, and soil respiration rate under MSPF for spring and autumn were increased by 20.46% and 49.22%, 31.02% and 58.05%, and 34.66% and 38.38%, respectively. In this study, greenhouse tomato under MSPF increased the soil respiration rate (CO2 emission); however, the tomato yield and WUE were also improved, indicating that the MSPF could achieve an effective dynamic balance among greenhouse CO2, yields and WUE. The micropore group spacing had no significant effect on the soil respiration rate. Compared with 50 cm of micropore group spacing, the 30 cm micropore group spacing increased the yield and WUE of tomato by 16.00% and 13.01% and 20.85% and 14.25%, respectively. With the increase of irrigation amount, the soil respiration rate in the root and nonroot zones and the yield of spring tomato and autumn tomato increased at first and then decreased, and the WUE showed a decreasing trend. As for the soil respiration rate, it is higher in the root zone than it is in the nonroot zone, and there exists a logarithmic function relationship between the soil respiration rate in the nonroot zone and the soil respiration rate in the root zone, and there is a quadratic relationship between the soil respiration rate and the yield in the nonroot zone. Considering our results comprehensively, the combination of the 30 cm micropore group spacing and 1.0 Epan is recommended for greenhouse tomato under MSPF in northwest China. This study can provide theoretical data and experimental support for water-saving irrigation and yield increase of facility agricultural tomato and analysis of the soil carbon cycle mechanism of MSPF. Although the results of this study explain the effects of different micropore group spacings and irrigation amounts on the soil respiration rate and yield of tomato, it is still an open question how to change the soil enzyme activity and microbial community structure to further drive the change of the soil respiration rate and yield under different micropore spacings and irrigation amounts.



---

*Source: 6658059-2021-01-21.xml*