Next Article in Journal
Experimental Study on the Characteristics of Sediment Transport and Sorting in Pressurized Pipes
Next Article in Special Issue
Evapotranspiration and Quantitative Partitioning of Spring Maize with Drip Irrigation under Mulch in an Arid Region of Northwestern China
Previous Article in Journal
OpenHi.net: A Synergistically Built, National-Scale Infrastructure for Monitoring the Surface Waters of Greece
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Evapotranspiration and Crop Coefficient of Chinese Cabbage Using Eddy Covariance in Northwest China

1
Center for Agricultural Water Research in China, China Agricultural University, Beijing 100083, China
2
Shiyanghe Experimental Station for Improving Water Use Efficiency in Agriculture, Ministry of Agriculture and Rural Affairs, Wuwei 733000, China
3
College of Land Science and Technology, China Agricultural University, Beijing 100083, China
*
Authors to whom correspondence should be addressed.
Submission received: 12 August 2021 / Revised: 12 September 2021 / Accepted: 23 September 2021 / Published: 7 October 2021

Abstract

:
Chinese cabbage is a key vegetable crop in northwest China. It is of great significance to study the evapotranspiration (ET) and crop coefficient (Kc) for agricultural water-saving management in this area. Eddy covariance (EC) was used to measure the ET and Kc of Chinese cabbage in northwest China from 1 May to 16 October 2020, in order to analyze the characteristics of these variables under plastic mulch. The results showed that the average Kc of the first crop of cabbage was higher in the middle and late stages, with values of 1.08 and 1.09 during the heading and maturity stages, respectively. The average Kc of the second crop of cabbage was higher in the middle stage, with values of 1.10 and 1.13 during the rosette and heading stages, respectively. The average annual Kc of Chinese cabbage was 0.81. Although Kc was higher in the middle and late periods, the water use efficiency was still 28.96 kg·ha−1·mm−1. The annual ET of Chinese cabbage was 505.3 mm. The study revealed the variation pattern of ET and Kc of Chinese cabbage, which provides an important scientific basis for the irrigation management of Chinese cabbage and is of great significance to guide the practice of water-saving vegetable planting.

1. Introduction

In recent decades, with the continuous population growth and development of the social economy, the demand for water resources has increased rapidly in world. Water resource deficiency, particularly in northwest China, leads to environmental degradation and limits regional economic development and social progress [1,2]. In the arid region of northwest China, about 10% of the area receives agricultural irrigation, and irrigation accounts for the largest part of available water consumption. At present, the area mainly pumps groundwater for irrigation and crop growth. How to improve the agricultural water use efficiency (WUE) has become the focus of scientists around the world [3,4,5,6,7,8,9]. Evapotranspiration (ET) is a main unknown variable for understanding the ecohydrological system. In arid areas, ET can reach 95% of the water balance [10]. The main form of agricultural water consumption is crop evapotranspiration (ET). Therefore, the key to water saving in agriculture is to reduce crop evapotranspiration. The determination of agricultural water management practices, the design of irrigation systems and irrigation regimes, and the calculation of crop yields are all based on accurate ET estimation [11]. A better understanding of ET can help make decisions about irrigation and improve the efficiency of available water [12,13,14,15]. The ratio of crop potential evapotranspiration to reference evapotranspiration (ET0) is defined as the crop coefficient (Kc) [11]. Kc is a key parameter for managing agricultural irrigation systems. Thus, studying the evapotranspiration process and estimating the crop coefficient (Kc) are key aspects to improve water use efficiency and save agricultural water in this area.
Previous studies on ET and Kc have mainly focused on field food crops such as maize and wheat in northwest China. Kang [16] showed that the average seasonal ET and Kc of maize in northwest China during the last decade were 424.0 mm and 1.04, respectively, according to lysimeter measurements. Li [17] measured the ET and Kc of spring maize in northwest China in 2007 and studied their characteristics under plastic film mulching. Guo [18] used the eddy covariance (EC) method to study the dynamics of crop coefficient of spring maize with plastic mulch in this region over the past 12 years. Yang [19] studied the effects of drip irrigation technology on soil evaporation (E), crop transpiration (T), and wheat total ET in this area. Studies on vegetable crops are more concentrated in terms of crop coefficient prediction and plant physiology. Pereira [20] predicted crop coefficients from a fraction of ground cover and height, as well as used ground and remote sensing data for background and verification. Pereira [21] also pointed out their practical application to vegetable, field, and fruit crops with a focus on parameterization. Wu [22] investigated the comparative responses to silicon and selenium in relation to different physiological processes in flowering Chinese cabbage under cadmium stress. Sun [23] studied the effect of Piriformospora indica on the drought tolerance of Chinese cabbage leaves.
At present, some changes have taken place in the crop planting structure in north-west China. The planting structure in this region has gradually changed from food crops to vegetable crops. Chinese cabbage, a subspecies of Brassica pekinensis (Lour.) Rupr., has been widely planted in Northwest China due to its high economic value. A lack of water is still a barrier for growing crops in this area. However, there are few reports about Chinese cabbage ET and Kc based on field observation. Water balance methods, such as the weighing lysimeter, enable direct estimations of ET and Kc [24]. However, the two measurements need to be at least one day apart, and their representativeness is poor. The EC method is regarded as a standard method for measuring ET. It can accurately capture ET information in a large range within a short period of time (e.g., 10 min) [25]. Therefore, the EC method was adopted in our study to measure the ET Chinese cabbage in northwest China. Then, Kc and other water consumption indices were calculated.
The aims of this experiment were (1) to compare the results of the water balance (WB) method and eddy covariance (EC) measurement of ET on Chinese cabbage, (2) to study the daily ET and seasonal ET variation of Chinese cabbage and its influencing factors, and (3) to analyze the seasonal variation of Kc.

2. Materials and Methods

2.1. Experimental Site and Design

The research was implemented at the Shiyanghe Experimental Station for Water Saving in Agriculture and Ecology of the China Agricultural University, located in Wuwei City, Gansu Province, northwest China (N 37°52′, E 102°50′, altitude 1581 m) from 1 May to 16 October 2020. The study area is in a typical temperate continental climate zone. The annual average temperature is 8 °C, the annual accumulated temperature (>0 °C) is 3550 °C, the annual precipitation is 164.4 mm, the average annual pan evaporation is about 2000 mm, the annual number of average sunshine hours is 3000 h, and the frost-free period is 150 days. The average groundwater depth is 25 m, and the data of soil texture are shown in Table 1. The region, suitable for the growth of half-hardy vegetables, is characterized by high altitude, low temperature, long sunshine duration, and suitable light and heat resources. Chinese vegetables, with yellow inner leaves and a short growth cycle, are suitable for growth in the environment of 5–25 °C. The vegetable has the characteristics of high yield and excellent quality.
Our study area spanned 500 m from north to south and 250 m from east to west. The planting density of Chinese cabbage is about 108,000 plants·ha−1. The row spacing was 40 cm and the planting spacing was 6.7 cm. In this study, the growth period of the first crop of Chinese cabbage was from 17 April to 1 July, and that of the second crop was from 1 August to 16 October. The study area was irrigated 14 times in 2020 on 3 May, 11 May, 19 May, 29 May, 3 June, 13 June, 19 June, 1 August, 12 August, 18 August, 27 August, 5 September, 14 September, and 24 September by drip irrigation with plastic mulch. The total volume of precipitation was 144.4 mm, and the total volume of irrigation water was 293.9 mm throughout the whole study period. The precipitation (P) and irrigation (I) volumes throughout the observation period are shown in Figure 1.

2.2. Using Eddy Covariance to Measure and Correct the Measured Value

The eddy covariance system was installed at about 320 m from the northernmost side of the experimental site and 145 m from the westernmost side, adjusting the relative height between the Chinese cabbage canopy and sensors at a constant 2.5 m. Chinese cabbage is the main crop type on experimental site, and the experimental field is large enough to provide an adequate fetch length for EC measurements. The smallest fetch length was 100 m. The mean wind direction was northwest. The eddy covariance system consisted of a temperature and humidity sensor (model HMP45C, Vaisala, Vantaa, Finland), a krypton hygrometer (model KH20, Campbell Scientific, Inc., Logan, UT, USA), and a 3D sonic thermometer/anemometer (model CSAT3, Campbell Scientific, Inc., USA). The HMP45C sensor can measure average temperature and humidity at 10 min intervals. The KH20 sensor and CSAT3 sensor can measure vertical fluctuations in temperature, water vapor density, and wind speed at 0.1 s intervals. All data are transferred to the same data logger (model CR5000, Campbell Scientific Inc, 100 USA), generating 10 min of statistical data (average, variance, and covariance). Continuous observations were conducted from 1 May to 16 October 2020.
The latent heat flux and sensible heat flux were calculated using the eddy correlation method [25].
λ E T = ρ a λ w q ¯ ,
H = C p ρ a w T ¯ ,
where λET and H are the latent and sensible heat flux (W·m−2), w q ¯ is the covariance between fluctuations of vertical wind speed w′ (m·s−1) and humidity q′ (kg·kg−1), w T ¯ is the covariance between fluctuations of w′ and sonic temperature T′ (K), ρa is the air density (kg·m−3), Cp is the specific heat of dry air at constant pressure (J·kg−1·K−1), λ is the latent heat of water vaporization (J·kg−1), and ET is the crop evapotranspiration (kg·m−2·s−1).
In this study, the original data measured by EC at 10 Hz were recalculated by Eddpro software (version 4.0, LI-COR Inc., Lincoln, NE, USA) to obtain the average latent heat flux, sensible heat flux, and momentum flux at time intervals of 30 min. The original data observed using the EC method were modified according to the Guide for Observation of Eddy Correlation Flux [26,27,28,29]. The following two interpolation methods were used to obtain the missing data in this study: (1) the linear interpolation method when the missing period was within 0–2 h; (2) the diurnal average method when the missing period was between 2 and 24 h.

2.3. Other Measurements

A net radiometer (model NR-LITE, Kipp & Zonen, Delft, The Netherlands), which was set at a relative height of 1.5 m above the canopy, was applied to measure the net radiation (Rn). Before the experiment, a net pyrgeometer (model CG2, Kipp & Zonen, Delft, The Netherlands) and a high-precision albedometer (model CM7B, Kipp & Zonen, Delft, The Netherlands) were used to calibrate the radiometer.
A soil heat flux plate (model HFP0, Hukseflux, The Netherlands) was set at a soil depth of 50 mm below the surface. Surface soil heat flux was calculated by correcting the value at 50 mm for the heat storage above the sensors, determined by the change in soil temperature above the heat flux sensors. The temperature above the soil heat flux plate was measured by four pairs of thermocouples (model 105T, Campbell Scientific, USA) located directly above the soil heat flux plate at depths of 20 mm and 40 mm. All sensors were sampled every 5 s, and the average values for soil heat flux and net radiation every 30 min were collected and calculated in CR5000. Observations were conducted continuously from 1 May to 16 October 2020.
Soil volumetric water content was monitored by soil moisture sensors (model CS616, Campbell Scientific, USA), and soil temperature was monitored by soil temperature sensors (model 109L, Campbell Scientific, USA). Soil temperature was monitored from 1 May to 16 October 2020 at soil depths of 20 cm, 40 cm, 60 cm, and 80 cm. Soil volumetric water content was monitored from 1 May to 17 July 2020 at soil depths of 20 cm, 40 cm, 60 cm, and 100 cm. It was also monitored from 17 July to 16 October 2020 at soil depths of 20 cm, 40 cm, 60 cm, and 80 cm. At the same time of monitoring, soil samples were collected regularly in the observation area at depths of 10 cm, 20 cm, 40 cm, 60 cm, 80 cm, and 100 cm. The soil moisture content of each layer was measured by the drying method to calibrate the data measured by the instrument. The diurnal variation in soil water content (SWC) during the whole observation period is shown in Figure 1.
White plastic mulch with a width of 1.2 m was applied in the field, and the spacing of each application was 0.4 m. Six weighing micro-evaporation barrels were applied to measure soil evaporation (ES) under the film and bare soil [30]. During the growth period of Chinese cabbage, four plants were selected to measure the growth indicators (plant height, canopy height, leaf length, leaf width, and number of leaves) every week. Lastly, Kc was calculated according to FAO 56, i.e., the ratio of the ET(ETEC) measured by the EC method to the reference ET(ET0) calculated by the Penman–Monteith formula [11].

2.4. Applying the Water Balance Method to Estimate ET

The water balance (WB) method based on mass conservation is a simple method to estimate ET [24,31,32], which can be expressed by Equation (5).
P + I + Q g E T D = Δ W ,
where P is the precipitation, I is the irrigation, Qg is the contribution from the water table, ET is the actual ET, D is the deep drainage, and ΔW is the change in water storage in the 0–100 cm soil layer, which can be calculated as
Δ W = W t 2 W t 1 ,
where W t 1 and W t 2 are the average water storage in the 0–100 cm soil layer at time t 1 and t 2 , respectively.
We ignored deep drainage in drip irrigation. The buried depth of the groundwater level is more than 40–50 m; thus, the contribution of the water table (Qg) could also be ignored. Thus, Equation (3) can be simplified as
E T W B = P + I Δ W .

2.5. Water Use Efficiency Calculation

Water use efficiency based on cumulative evapotranspiration is calculated as
W U E = Y E T ,
where Y is the yield (kg·ha1), and ET is the cumulative evapotranspiration (mm).

2.6. Data Analysis Method

Excel data analysis software was used to carry out polynomial fitting analysis on ET, Es, and Kc in the observation area, and regression analysis was conducted on ET and its influencing factors. The statistical index used was the determination coefficient (R2).
R 2 = i = 1 N O i O ¯ P i P ¯ i = 1 N O i O ¯ 2 i = 1 N P i P ¯ 2 2 ,
where O i is the measured value, P i is the fitted value, O ¯ is the average of the measured value, and P ¯ is the average of the fitted value.

3. Results and Discussion

3.1. The Examination of ET Measured by Eddy Covariance

The water balance method is generally regarded as an accurate method for calculating ET over a long period of time [25]. In our study, the results from the water balance method were used to test whether the ET measured by the eddy covariance method was reliable. Figure 2 shows that the weekly ETEC and ETWB were consistent at different growth stages (ETEC = 0.94ETWB, R2 = 0.59, n = 19).
As can be seen from Table 2, the total ETEC and ETWB throughout the growth period of the first crop were 275.6 mm and 236.3 mm, respectively. The ratio of ETEC to ETWB at the seedling stage, rosette stage, heading stage, and maturity stage were 0.48, 3.66, 1.05, and 1.34, respectively, and the ratio throughout the growth period was 1.17. The total ETEC and ETWB of the second crop were 229.7 mm and 192.8 mm throughout the growth period, respectively. The ratio of ETEC to ETWB at the seedling stage, rosette stage, heading stage, and maturity stage were 2.84, 0.91, 1.38, and 0.74, respectively, and the ratio throughout the growth period was 1.19. These results showed that the ETEC was in good agreement with ETWB, especially throughout the growth period. This proves that the EC method could accurately measure the Chinese cabbage’s ET in the study area.

3.2. Diurnal and Seasonal Variations of Chinese Cabbage’s ET and Influencing Factors

Figure 3a,b show the daily evapotranspiration change of two crops at different growth stages, which presented a tendency to rise first and then to fall. The maximum daily ET at the seedling stage, rosette stage, heading stage, and maturity stage of the first crop was 0.12 mm·h−1, 0.17 mm·h−1, 0.22 mm·h−1, and 0.27 mm·h−1, respectively, and the corresponding times were 1:00, 12:00, 2:00, and 1:00 p.m., respectively. The maximum daily ET at the seedling stage, rosette stage, heading stage, and maturity stages of the second crop was 0.21 mm·h−1, 0.22 mm·h−1, 0.22 mm·h−1, and 0.22 mm·h−1, respectively, and the corresponding times were all 1:00 p.m. This may have been caused by the change in sunshine duration and solar radiation during the different periods. It was observed that ET at night was approximately zero at different stages, with the exception (below zero) of the first crop seedling period, which may have been due to the large temperature difference between day and night in early summer and the condensation of water vapor at night.
As can be seen from Figure 4a,b, the variation in ET throughout the growth period of the two crops showed a cubic parabolic trend with the relative growing day (first crop: ET = −12.53RGD3 + 9.66RGD2 + 7.56RGD + 0.57, R2 = 0.67; second crop: ET = −9.39RGD3 + 0.51RGD2 + 7.61RGD + 1.40, R2 = 0.63). The daily ET gradually increased from the seedling stage to the heading stage but decreased after the heading stage. The maximum value of ET in the first crop was 8.91 mm·day−1 at the heading stage (4 June), and the minimum value of ET in the first crop was 0.67 mm·day−1 at the seedling stage (3 May). In the second crop, the maximum value of ET was 5.45 mm·day−1 at the heading stage (11 September), and the minimum value of ET was 0.16 mm·day−1 at the maturity stage (17 October), which may have been mainly due to the late harvesting of the mature Chinese cabbage.
Figure 5a,b also showed a parabolic trend of daily Es with relative growing day throughout the growth period of the two crops (first crop: Es = 7.45RGD2 − 12.42RGD + 5.72, R2 = 0.34; second crop: Es = 75.72RGD2 − 82.12RGD + 22.56, R2 = 0.35). Es decreased from the seedling to late rosette stage and increased in the late growth stage. This change was mainly caused by the fact that the enlarged leaf area covered the ground. In the later stage, the increase in Es was due to the senescence of green leaves and the decomposition of plastic film.
As can be seen from Table 2, the total ETEC of the first crop at the seedling, rosette, heading, and maturity stages was 32.3 mm, 53.9 mm, 66.0 mm, and 123.5 mm, respectively, and the corresponding average daily ETEC was 1.8 mm·day−1, 4.5 mm·day−1, 6.6 mm·day−1, and 5.9 mm·day−1, respectively. Throughout the growth period, the amount of total ET was 275.6 mm, and the daily average ET was 4.5 mm·day−1. The total ETEC of the second crop at the seedling stage, rosette stage, heading stage, and maturity stage was 83.8 mm, 43.7 mm, 35.7 mm, and 66.6 mm, respectively, and the corresponding daily average ETEC was 2.8 mm·day−1, 4.4 mm·day−1, 4.5 mm·day−1, and 2.4 mm·day−1, respectively. Throughout the growth period, the total ET was 229.7 mm, and the daily average ET was 3.0 mm·day−1.
Guo [18] showed that the annual average ET of spring maize was 481.7 mm. The study of Yang [19] estimated that the annual average ET of wheat in this area was 430.0 mm. In this study, the annual ET of Chinese cabbage was 505.3 mm, only 23.6 mm higher than that of maize, and 75.3 mm higher than that of wheat. Furthermore, the water use efficiency (WUE) of Chinese cabbage was 28.96 kg·ha−1·mm−1, which is higher than the previous result of 25.2 kg·ha−1·mm−1 for spring maize [17]. The water consumption of field vegetable crops is not much higher than that of field food crops, although the water use efficiency of field vegetable crops is higher than that of food crops. Therefore, a change in the planting structure in this area would not cause a significant increase in agricultural irrigation water consumption.
Crop ET is related to many factors, such as crop growth, surface characteristics, and meteorological factors. To clarify the influencing factors of Chinese cabbage’s ET in the study area, Table 3 lists the relationship of net radiation (W·m−2)–ET and temperature (°C)–ET for two crops of Chinese cabbage. Daily ET had a linear relationship with Rn (first crop: R2 = 0.53; second crop: R2 = 0.47) and an exponential relationship with Ta (first crop: R2 = 0.34; second crop: R2 = 0.14). Table 3 also shows the results of multiple linear regression between ET and these factors (first crop: ET = 0.03Rn + 0.19Ta − 2.81, R2 = 0.6; second crop: ET = 0.03Rn − 0.04Ta + 0.82, R2 = 0.48). Law [33] pointed out that Rn was the main influencing factor on ET, which is consistent with our conclusion. Temperature had no significant effect on ET of Chinese cabbage in the study area. In addition, the multiple regression equation established in the study can be used as an experiential formula to calculate ET in this area.

3.3. Seasonal Variation of Chinese Cabbage’s Kc and Its Influencing Factors

Figure 6a,b show that the Kc of the two crops throughout the growing period exhibited a cubic parabolic trend with the change in relative growing days (first crop: Kc = −1.42RGD3 + 0.46RGD2 + 1.83RGD + 0.11, R2 = 0.87, n = 62; second crop: Kc = −4.65RGD3 + 3.71RGD2 + 0.76RGD + 0.35, R2 = 0.86, n = 77). The daily Kc of the first crop showed an increasing trend from the seedling to heading stage and a small decreasing trend after the heading stage. The maximum value of Kc was 1.19 at the heading stage (8 June), and the minimum value of Kc was 0.15 at the seedling stage (1 May). The maximum Kc on 8 June was probably caused by continuous rainfall on 7 and 8 June. The daily Kc was higher in the middle and late periods, ranging from 0.93 to 1.19. The daily Kc of the second crop showed an increasing trend from the seedling stage to heading stage and a gradual decreasing trend after maturity stage. The maximum value of Kc was 1.24 at the rosette stage (7 September) and the minimum value of Kc was 0.07 at the maturity stage (16 October). The daily Kc was higher in the rosette and heading periods, ranging from 0.98 to 1.24.
As can be seen from Table 2, the daily average Kc of the first crop of Chinese cabbage at the seedling stage, rosette stage, heading stage, and maturity stages was 0.35, 0.93, 1.08, and 1.09, respectively, and the daily average Kc throughout the growth period was 0.85. The average daily Kc of the second crop at the seedling stage, rosette stage, heading stage, and maturity stage was 0.62, 1.10, 1.13, and 0.84, respectively, and the daily average Kc throughout the growth period was 0.80. Kc of the first crop in the seedling stage was far smaller than that of the second crop. That situation may have come from the fact that the period just after sowing was missed, because the first crop was sowed on 17 April, but the observation started on 1 May. There was no significant difference in Kc value between the two crops at the rosette and the heading stages. At the maturity stage before harvest, the Kc value of the first crop of Chinese cabbage was still high, which might have been due to the antiaging effect of plastic film mulching on Chinese cabbage. The Kc of Chinese cabbage was a bit small in the second crop, which may have been due to the late harvest of the second crop of Chinese cabbage for cost saving, which led to Kc dropping rapidly. Guo [18] showed that the average Kc of spring maize was 0.88. The average Kc of Chinese cabbage was 0.81, which was 0.07 lower than that of maize. A previous study found that Kc mainly depends on crops growth, but is also affected by soil moisture [34]. All these influencing factors may underlie the complex changes in Kc.
Table 4 shows that irrigation (I) greatly improved Kc. Kc of the first crop of Chinese cabbage changed from 0.16 to 0.45, from 0.37 to 0.40, from 0.61 to 0.83, from 0.78 to 1.04, from 1.10 to 1.14, and from 1.12 to 1.10 on the days after the first, second, third, fourth, fifth, and sixth irrigations, respectively. Unlike the irrigation at the five other timepoints, the Kc value decreased after the sixth irrigation, which may have been related to the small irrigation amount at the maturity stage. For the second crop, the Kc value increased from 0.13 to 0.56, from 0.56 to 0.60, from 0.58 to 0.47, from 0.80 to 0.99, from 1.02 to 1.19, and from 1.11 to 1.17 on the days after the first, second, third, fourth, fifth, and sixth irrigations, respectively. The value remained basically unchanged after the seventh irrigation.

4. Conclusions

According to our research, the following conclusions can be drawn:
  • The weekly ETEC and ETWB at different growth periods had good consistency, and the consistency was better throughout the growth period. The ratio of ETEC to ETWB throughout the growth period of the first crop was 1.17 and that of the second crop was 1.19.
  • The diurnal variation of ET presented a tendency to rise first and then to fall. The total ET of the first crop throughout the growth period was 275.6 mm, the average daily ET was 4.52 mm·day−1, and the WUE was 31.58 kg·ha−1·mm−1. The total ET of the first crop throughout the growth period was 229.7 mm, the average daily ET was 3.02 mm·day−1, and the WUE was 25.91 kg·ha−1·mm−1. Daily ET had a linear relationship with Rn and an exponential relationship with Ta. Rn was the main influencing factor on ET.
  • Daily Kc showed a cubic curve trend throughout the growth period. The average daily Kc at the four growth stages of the first crop was 0.35, 0.93, 1.08, and 1.09, respectively, and the average Kc throughout the growth period was 0.85. The average daily Kc at the four growth stages of the second crop was 0.62, 1.10, 1.13, and 0.84, respectively, and the average Kc throughout the growth period was 0.80. The lower Kc in the early stage of the first crop was due to the late observation time, while the lower Kc in the late stage of the second crop was due to the late harvest. Kc was greatly and positively affected by irrigation.
Our study revealed the variation pattern of ET and Kc of Chinese cabbage, which pro-vides an important scientific basis for the irrigation management of Chinese cabbage. The eddy covariance method was innovatively applied to reveal the law of water consumption of Chinese cabbage, which provides a practical basis for the design of vegetable crop irrigation systems. This is of great significance to guide the practice of water-saving vegetable planting. Our experimental observations were relatively short, with only two crop growth periods. In the future, more observational studies will be conducted to improve the credibility of the study and further propose field management measures.

Author Contributions

Conceptualization, J.D. and S.L.; methodology, J.D.; validation, J.D., S.L. and H.W.; investigation, C.W., Y.Z. and D.Y.; resources, S.L.; data curation, J.D.; writing—original draft preparation, J.D.; writing—review and editing, S.L and H.W.; supervision, C.W.; project administration, S.L.; funding acquisition, S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Chinese National Natural Science Fund (51879262, 51725904, 51790530).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

We wish to thank the staff of the Shiyang River Basin Experimental Station for Agriculture and Ecological Water Saving of China Agriculture University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kang, S.Z.; Su, X.L.; Tong, L.; Shi, P.Z.; Yang, X.Y.; Abe, Y.K.; Du, T.S.; Shen, Q.L.; Zhang, J.H. The impacts of human activities on the water-land environment of the Shiyang River basin, an arid region in northwest China. Hydrol. Sci. J. 2004, 49, 413–427. [Google Scholar] [CrossRef]
  2. Deng, X.P.; Shan, L.; Zhang, H.P.; Turner, N.C. Improving agricultural water use efficiency in arid and semiarid areas of China. Agric. Water Manag. 2005, 80, 23–40. [Google Scholar] [CrossRef]
  3. Molden, D.; Oweis, T.; Steduto, P.; Bindraban, P.; Hanjra, M.A.; Kijne, J. Improving agricultural water productivity: Between optimism and caution. Agric. Water Manag. 2009, 97, 528–538. [Google Scholar] [CrossRef]
  4. Pereira, L.S.; Cordery, I.; Iacovides, I. Improved indicators of water use performance and productivity for sustainable water conservation and saving. Agric. Water Manag. 2012, 108, 39–51. [Google Scholar] [CrossRef]
  5. Levidow, L.; Zaccaria, D.; Maia, R.; Vivas, E.; Todorovic, M.; Scardigno, A. Improving water-efficient irrigation: Prospects and difficulties of innovative practices. Agric. Water Manag. 2014, 146, 84–94. [Google Scholar] [CrossRef] [Green Version]
  6. Du, T.S.; Kang, S.Z.; Zhang, J.H.; Davies, W.J. Deficit irrigation and sustainable water-resource strategies in agriculture for China’s food security. J. Exp. Bot. 2015, 66, 2253–3369. [Google Scholar] [CrossRef] [PubMed]
  7. Kang, S.Z.; Hao, X.M.; Du, T.S.; Tong, L.; Su, X.L.; Lu, H.N.; Li, X.L.; Huo, Z.L.; Li, S.E.; Ding, R.S. Improving agricultural water productivity to ensure food security in China under changing environment: From research to practice. Agric. Water Manag. 2017, 179, 5–17. [Google Scholar] [CrossRef]
  8. Ryosuke, M.; Fumitaka, A.; June, K.; Yuuri, T.; Keisuke, T.; Hisato, K.; Yoichi, S.; Kousuke, H.; Hisashi, T.; Jun, K.; et al. Tuning water-use efficiency and drought tolerance in wheat using abscisic acid receptors. Nat. Plants 2019, 5, 153–159. [Google Scholar]
  9. Kool, D.; Agam, N.; Lazarovitch, N.; Heitman, J.L.; Sauer, T.J.; Ben-Gal, A. A review of approaches for evapotranspiration partitioning. Agric. For. Meteorol. 2014, 184, 56–70. [Google Scholar] [CrossRef]
  10. Wilcox, B.P.; Breshears, D.D.; Seyfried, M.S. Encyclopedia of Water Science; Marcel Dekker, Inc.: New York, NY, USA, 2003; pp. 791–794. [Google Scholar]
  11. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for Computing Crop Requirements; Irrigation and Drainage Paper No. 56; FAO: Roma, Italy, 1998; 331p. [Google Scholar]
  12. Kite, G. Using a basin-scale hydrological model to estimate crop transpiration and soil evaporation. J. Hydrol. 2000, 229, 59–69. [Google Scholar] [CrossRef]
  13. Zhao, N.N.; Liu, Y.; Cai, J.B.; Paredes, P.; Rosa, R.D.; Pereira, L.S. Dual crop coefficient modelling applied to the winter wheat–summer maize crop sequence in North China Plain: Basal crop coefficients and soil evaporation component. Agric. Water Manag. 2013, 117, 93–105. [Google Scholar] [CrossRef]
  14. Pereira, L.S.; Paredes, P.; López-Urrea, R.; Hunsaker, D.J.; Mota, M.; Mohammadi Shad, Z. Standard single and basal crop coefficients for vegetable crops, an update of FAO56 crop water requirements approach. Agric. Water Manag. 2020, 243, 106196. [Google Scholar]
  15. Pereira, L.S.; Paredes, P.; Hunsaker, D.J.; López-Urrea, R.; Mohammadi Shad, Z. Standard single and basal crop coefficients for field crops. Updates and advances to the FAO56 crop water requirements method. Agric. Water Manag. 2020, 243, 106466. [Google Scholar]
  16. Kang, S.Z.; Gu, B.J.; Du, T.S.; Zhang, J.H. Crop coefficient and ratio of transpiration to evapotranspiration of winter wheat and maize in a semi-humid region. Agric. Water Manag. 2003, 59, 239–254. [Google Scholar] [CrossRef]
  17. Li, S.E.; Kang, S.Z.; Li, F.S.; Zhang, L. Evapotranspiration and crop coefficient of spring maize with plastic mulch using eddy covariance in northwest China. Agric. Water Manag. 2008, 95, 1214–1222. [Google Scholar] [CrossRef]
  18. Guo, H.; Li, S.E.; Kang, S.Z.; Du, T.S.; Hao, X.M.; Ding, R.S. Crop coefficient for spring maize under plastic mulch based on 12-year eddy covariance observation in the arid region of Northwest China. J. Hydrol. 2020, 588, 125108. [Google Scholar] [CrossRef]
  19. Yang, D.N.; Li, S.E.; Kang, S.Z.; Du, T.S.; Guo, P.; Mao, X.M.; Tong, L.; Hao, X.M.; Ding, R.S.; Niu, J. Effect of drip irrigation on wheat evapotranspiration, soil evaporation and transpiration in Northwest China. Agric. Water Manag. 2020, 232, 106001. [Google Scholar] [CrossRef]
  20. Pereira, L.S.; Paredes, P.; Melton, F.; Johnson, L.; Wang, T.; López-Urrea, R.; Cancela, J.J.; Allen, R. Prediction of crop coefficients from fraction of ground cover and height. Background and validation using ground and remote sensing data. Agric. Water Manag. 2020, 241, 106197. [Google Scholar] [CrossRef]
  21. Pereira, L.S.; Paredes, P.; Melton, F.; Johnson, L.; Mota, M.; Wang, T. Prediction of crop coefficients from fraction of ground cover and height. Practical application to vegetable, field and fruit crops with focus on parameterization. Agric. Water Manag. 2020, 241, 106663. [Google Scholar]
  22. Wu, Z.C.; Liu, S.; Zhao, J.; Wang, F.H.; Du, Y.Q.; Zou, S.M.; Li, H.M.; Wen, D.; Huang, Y.D. Comparative responses to silicon and selenium in relation to antioxidant enzyme system and the glutathione-ascorbate cycle in flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis) under cadmium stress. Environ. Exp. Bot. 2017, 133, 1–11. [Google Scholar] [CrossRef]
  23. Sun, C.; Johnson, J.M.; Cai, D.G.; Sherameti, I.; Oelmüller, R.; Lou, B. Piriformospora indica confers drought tolerance in Chinese cabbage leaves by stimulating antioxidant enzymes, the expression of drought-related genes and the plastid-localized CAS protein. J. Plant Physiol. 2010, 167, 1009–1017. [Google Scholar] [CrossRef]
  24. Baldocchi, D.D. Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems: Past, present and future. Glob. Chang. Biol. 2003, 9, 479–492. [Google Scholar] [CrossRef] [Green Version]
  25. Boast, C.W.; Robertson, T.M. A ‘micro-lysimeter’ method for determining evaporation from bare soil: Description and laboratory evaluation. Soil Sci. Soc. Am. J. 1982, 46, 689–696. [Google Scholar] [CrossRef]
  26. Yu, G.R.; Sun, X.M. Principles of Flux Measurement in Terrestrial Ecosystems; Higher Education Press: Beijing, China, 2006. [Google Scholar]
  27. Qin, S.J.; Li, S.E.; Kang, S.Z.; Du, T.S.; Tong, L.; Ding, R.S.; Wang, Y.H.; Guo, H. Transpiration of female and male parents of seed maize in northwest China. Agric. Water Manag. 2019, 213, 397–409. [Google Scholar] [CrossRef]
  28. Wang, Y.H.; Li, S.E.; Qin, S.J.; Guo, H.; Yang, D.N.; Lam, H.M. How can drip irrigation save water and reduce evapotranspiration compared to border irrigation in arid regions in northwest China. Agric. Water Manag. 2020, 239, 106256. [Google Scholar] [CrossRef]
  29. Wang, Y.H.; Li, S.E.; Liang, H.; Hu, K.L.; Qin, S.J.; Guo, H. Comparison of Water- and Nitrogen-Use Efficiency over Drip Irrigation with Border Irrigation Based on a Model Approach. Agronomy 2020, 10, 1890. [Google Scholar] [CrossRef]
  30. Rana, G.; Katerji, N.; Mastrorilli, M. Environmental and soil-plant parameters for modeling actual evapotranspiration under water stress conditions. Ecol. Modeling 1997, 101, 363–371. [Google Scholar] [CrossRef]
  31. Rana, G.; Katerji, N.; Mastrorilli, M.; Moujabber, M.E.; Brisson, N. Validation of model of actual evapotranspiration for water stressed soybeans. Agric. For. Meteorol. 1997, 86, 215–224. [Google Scholar] [CrossRef]
  32. Rana, G.; Katerji, N. Measurement and estimation of actual evapotranspiration in the field under Mediterranean climate: A review. Eur. J. Agron. 2000, 13, 125–153. [Google Scholar] [CrossRef]
  33. Law, B.E.; Falge, E.; Gu, L. Environmental controls over carbon dioxide and water vapor exchange of terrestrial vegetation. Agric. For. Meteorol. 2002, 113, 97–120. [Google Scholar] [CrossRef] [Green Version]
  34. Han, S.M.; Yang, Y.H.; Li, H.L.; Yang, Y.M.; Wang, J.S.; Cao, J.S. Determination of crop water use and coefficient in drip-irrigated cotton fields in arid regions. Field Crop. Res. 2019, 236, 85–95. [Google Scholar] [CrossRef]
Figure 1. The precipitation (P) and irrigation (I) volumes and the diurnal variation in soil water content (SWC) of Chinese cabbage throughout the observation period.
Figure 1. The precipitation (P) and irrigation (I) volumes and the diurnal variation in soil water content (SWC) of Chinese cabbage throughout the observation period.
Water 13 02781 g001
Figure 2. Comparison of weekly ET of Chinese cabbage from eddy covariance measurement and water balance method.
Figure 2. Comparison of weekly ET of Chinese cabbage from eddy covariance measurement and water balance method.
Water 13 02781 g002
Figure 3. (a) Diurnal ET variation of Chinese cabbage at different growth stages of the first crop; (b) diurnal ET variation of Chinese cabbage at different growth stages of the second crop.
Figure 3. (a) Diurnal ET variation of Chinese cabbage at different growth stages of the first crop; (b) diurnal ET variation of Chinese cabbage at different growth stages of the second crop.
Water 13 02781 g003
Figure 4. (a) Seasonal ET variation of Chinese cabbage throughout the growing period of the first crop. Relative growing days (RGDs) refer to the normalized growing days, e.g., “1” represents 61 days after sowing; (b) seasonal ET variation of Chinese cabbage throughout the growing period of the second crop.
Figure 4. (a) Seasonal ET variation of Chinese cabbage throughout the growing period of the first crop. Relative growing days (RGDs) refer to the normalized growing days, e.g., “1” represents 61 days after sowing; (b) seasonal ET variation of Chinese cabbage throughout the growing period of the second crop.
Water 13 02781 g004
Figure 5. (a) Seasonal Es variation throughout growing period of the first crop; (b) seasonal Es variation throughout growing period of the second crop.
Figure 5. (a) Seasonal Es variation throughout growing period of the first crop; (b) seasonal Es variation throughout growing period of the second crop.
Water 13 02781 g005
Figure 6. (a) Seasonal Kc variations of Chinese cabbage throughout the growing period of the first crop; (b) seasonal Kc variations of Chinese cabbage throughout the growing period of the second crop.
Figure 6. (a) Seasonal Kc variations of Chinese cabbage throughout the growing period of the first crop; (b) seasonal Kc variations of Chinese cabbage throughout the growing period of the second crop.
Water 13 02781 g006
Table 1. Soil texture of experimental sites.
Table 1. Soil texture of experimental sites.
Depth (cm)Soil TextureDry Bulk Density (g·cm−3)
0–10Sandy loam1.60
10–20Silty loam1.63
20–40Silty loam1.57
40–60Silty loam1.41
60–80Silt1.43
80–100Silty loam1.69
Table 2. ET of Chinese cabbage from water balance method, as well as ET and Kc from eddy covariance measurement, over the whole growth period.
Table 2. ET of Chinese cabbage from water balance method, as well as ET and Kc from eddy covariance measurement, over the whole growth period.
CropGrowth StagePeriodDaysCumulative ET (mm)Daily Average ET (mm·Day−1)Average Kc
ETWBETECETWBETECET0Kc
First cropSeedling1–19 May1866.832.33.71.891.50.35
Rosette20–31 May1214.753.91.24.557.80.93
Heading1–10 June1062.666.06.36.660.81.08
Maturity11 June–1 July2192.1123.54.45.9113.41.09
Whole growth stage1 May–1 July61236.3275.63.94.5323.50.85
Second cropSeedling1–31 August3029.583.81.02.8135.60.62
Rosette1–10 September1048.143.74.84.439.71.10
Heading11–18 September825.835.73.24.531.691.13
Maturity19 September–16 October2889.466.63.22.479.60.84
Whole growth stage1 August–16 October76192.8229.72.53.0286.50.80
Table 3. Regression analysis between daily ET and its influencing factors.
Table 3. Regression analysis between daily ET and its influencing factors.
First CropRegression EquationR2n
Net radiation (Rn)ET = 0.03Rn − 0.300.5361
Air temperature (Ta)ET = 0.57e0.10Ta0.3461
Multiple linear regressionET = 0.03Rn + 0.20Ta − 2.810.6061
Second CropRegression EquationR2n
Net radiation (Rn)ET = 0.0247Rn + 0.49880.4777
Air temperature (Ta)ET = 1.1975e0.0468Ta0.1477
Multiple linear regressionET = 0.03Rn − 0.04Ta + 0.820.4877
Whole CropRegression EquationR2n
Net radiation (Rn)ET = 0.03Rn − 0.010.47168
Air temperature (Ta)ET = 1.63e0.03Ta0.03168
Multiple linear regressionET = 0.03Rn − 0.06Ta + 0.740.48168
Table 4. Irrigation effect on Kc.
Table 4. Irrigation effect on Kc.
First CropIrrigation3rd DBI *2nd DBI1st DBIIrrigation Day1st DAI #2nd DAI3rd DAI
Kc1st irrigation 0.150.220.160.450.400.42
2nd irrigation0.410.290.240.370.400.460.45
3rd irrigation0.320.370.590.610.830.790.84
4th irrigation1.061.041.000.781.041.051.03
5th irrigation1.051.031.091.101.141.231.22
6th irrigation0.961.121.151.121.100.921.10
Second CropIrrigation3rd DBI2nd DBI1st DBIIrrigation Day1st DAI2nd DAI3rd DAI
Kc1st irrigation0.100.120.230.130.560.470.48
2nd irrigation0.380.400.290.560.600.550.43
3rd irrigation0.430.560.770.580.470.720.55
4th irrigation0.860.800.840.800.990.771.00
5th irrigation1.051.001.021.121.091.241.11
6th irrigation1.210.981.111.171.151.131.14
7th irrigation0.881.011.111.091.081.101.00
* DBI represents days before irrigation; # DAI represents days after irrigation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ding, J.; Li, S.; Wang, H.; Wang, C.; Zhang, Y.; Yang, D. Estimation of Evapotranspiration and Crop Coefficient of Chinese Cabbage Using Eddy Covariance in Northwest China. Water 2021, 13, 2781. https://0-doi-org.brum.beds.ac.uk/10.3390/w13192781

AMA Style

Ding J, Li S, Wang H, Wang C, Zhang Y, Yang D. Estimation of Evapotranspiration and Crop Coefficient of Chinese Cabbage Using Eddy Covariance in Northwest China. Water. 2021; 13(19):2781. https://0-doi-org.brum.beds.ac.uk/10.3390/w13192781

Chicago/Turabian Style

Ding, Jie, Sien Li, Hongshuo Wang, Chunyu Wang, Yunxuan Zhang, and Danni Yang. 2021. "Estimation of Evapotranspiration and Crop Coefficient of Chinese Cabbage Using Eddy Covariance in Northwest China" Water 13, no. 19: 2781. https://0-doi-org.brum.beds.ac.uk/10.3390/w13192781

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

Article Metrics

Back to TopTop