Next Article in Journal
Single- versus Double-Species Cover Crop Effects on Soil Health and Yield in Mississippi Soybean Fields
Next Article in Special Issue
The Impact of Drought Stress on Soil Microbial Community, Enzyme Activities and Plants
Previous Article in Journal
Modelling the Effect and Variability of Integrated Weed Management of Phalaris minor in Rice-Wheat Cropping Systems in Northern India
Previous Article in Special Issue
Adoption of Conservation Agriculture in Rwanda: A Case Study of Gicumbi District Region
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Evapotranspiration and Crop Coefficient of Rain-Fed Tea Plants under a Subtropical Climate

1
Key Laboratory of Tea Biology and Resource Utilization of Tea (Ministry of Agriculture), Tea Research Institute, Chinese Academy of Agriculture Sciences, Hangzhou 310008, China
2
Graduate School of Chinese Academy of Agriculture Sciences, Beijing 100081, China
3
Lishui Academy of Agricultural and Forestry Sciences, Lishui 323000, China
4
Lishui Vocational and Technical College, Lishui 323000, China
5
Xihu National Agricultural Experimental Station for Soil Quality, Hangzhou 310008, China
*
Author to whom correspondence should be addressed.
Submission received: 12 October 2021 / Revised: 12 November 2021 / Accepted: 16 November 2021 / Published: 18 November 2021
(This article belongs to the Special Issue Agroecology and Organic Horticulture)

Abstract

:
Crop coefficient (Kc) is one of the most significant parameters for crop water demand prediction and irrigation scheduling. However, there is a lack of knowledge about water vapor and Kc in tea plantation ecosystems. This study explored and determined the actual evapotranspiration (ETc) and Kc of two cultivars of tea (Camellia sinensis), clone variety Baiye1 (BY1) and Longjing43 (LJ43), based on lysimeter data. An estimation was made for both for ETc and adjusted ETc, and the corresponding Kc and adjusted Kc. The results showed that the adjusted ETc and Kc values revealed a minor fluctuation when compared to the ETc and Kc values during the experimental period, which indicated that the adjusted ETc and Kc values were more precise and practical to field conditions. The average adjusted Kc values were 0.71 (range of 0.43–1.02) for BY1 and 0.84 (range of 0.48–1.22) for LJ43. Additionally, heavy pruning can decrease ETc and Kc values, possibly due to the lower level of LAI after pruning. Moreover, it is clearly manifested that BY1 consistently had lower ETc and Kc values than those of LJ43 because of plant growth status differences between BY1 and LJ43. Overall, our study proposed a reliable reference of Kc in tea plantation, and illuminated the effects of pruning and plant growth differences on Kc, which could provide a strong basis for precise irrigation in tea plantations in a subtropical climate.

1. Introduction

As a perennial evergreen plant originated from Southwestern China, tea plants (Camellia sinensis) favor a warm and humid climate [1]. However, to meet the water demand of tea tree growth, the temporal distribution of rainfall over the year is more vital [2,3,4]. Droughts usually cause irreparable yield losses, as well as prematurely aged tea plants [5,6,7,8,9], while heavy rainfall, on the other hand, can cause water logging, runoff, topsoil erosion, and deep percolation of critical soil nutrients [10]. In subtropical regions of China, where most tea plants are cultivated, the monthly rainfall always unevenly distributes, with frequent droughts in autumn and excessive rainfall in summer [11]. Therefore, irrigation, as a countermeasure to the droughts, is of great importance to tea plant growth.
It is necessary to schedule a proper irrigation regime according to the water demand of tea plants. Crop water demand is related to its real evapotranspiration (ETc), which combines water loss from the ground surface (evaporation) and plants (transpiration) with the atmosphere in the shape of vapor [12]. The accurate determination of ETc and crop coefficient (Kc) is doubtlessly a benefit for assessing the water requirement of tea. It could help tea producers to decide on the optimal time and water amount of irrigation. This is of great significance for attaining the growers’ purposes. In addition, the accurate determination of ETc and Kc depends on the local conditions, and can increase tea yields, water use efficiency, and profits, or reduce energy use, costs, and negative environmental influences [13]. Considering the fact that the majority of Chinese tea plantations are situated in mountain areas, where tea plants can only receive either rainwater or irrigated rainwater, much attention needs to be paid to the determination of ETc and Kc in the rain-fed tea plants, so as to improve the water use efficiency of the limited or uneven rainwater resource.
In practice, the crop ETc could be estimated as the proportion of the reference evapotranspiration (ET0), which is always calculated from the Penman–Monteith equation. Kc is defined as the ratio of ETc to ET0 [14,15]. The Kc of tea plants, recommended by the FAO paper, is from 0.9 to 1.0 [14]. However, the Kc is prone to being affected by climate, crop growth, and the reference method [16,17]. In situ and fast frequency measurement techniques are always used in the determination of water flux, e.g., Bowen ratio–energy balance method [18], static chamber/infrared gas analyzer [19], and eddy covariance technique [20]. With the advent of new satellite technology, remote sensing approaches, together with various surface energy balance algorithms, are also employed to estimate the ETc and Kc from landscapes and regions [21,22,23,24]. However, the above-mentioned methods lack accuracy, although they exhibit high temporal resolution and are suitable for annual cereal crops. Alternatively, the direct water balance method, i.e., lysimeter method, could obtain a more accurate result [4,25,26]. However, there was no Kc of the tea obtained through the lysimetric data reported in China, even though China is the largest tea production country in the world. Furthermore, the effect of tea pruning practice on Kc is still an open question, since periodical pruning has been reported to alter the model of energy sharing in tea plantations [27].
In this study, an attempt was made to estimate the ETc and Kc of tea plantations in a subtropical region, based on in situ annual lysimeter monitoring. The objectives of this study were as follows: (1) to estimate the crop coefficient of tea (Camellia sinensis) through the water balance equation, based on a field lysimeter study; (2) to check the effect of tea cultivars and tea pruning practice on the crop coefficient of tea.

2. Materials and Methods

2.1. Study Site

The field experiment was conducted at Shengzhou Experimental Station of Tea Research Institute (29.74° N, 120.82° E, 23 m above sea level), affiliated with Chinese Academy of Agricultural Sciences (TRI-CAAS), Shaoxing, China. The experimental site was in a typical subtropical region, with average annual rainfall of 1200 mm and an average annual temperature of 12.6 °C, which is favorable for tea plant growth. Soil properties before the experiment are shown in Table 1.
Seedlings (1 year) were transplanted in the studied field in November 2015. The tea cultivars were clone variety Baiye1 (BY1) and Longjing43 (LJ43) (Figure S1a), respectively. The tea trees were planted in rows with 1.5 m between rows and 0.33 m between plants at a density of approximately 20,000 plants per hectare. Each variety was replicated in 9 rectangular plots (36 m2 in area) (Figure S1b,c). After spring tea harvest, tea pruning was uniformly conducted in May 2019 and 2020 at a height that was 45 cm above the ground.

2.2. Sampling and Measurement

The 18 lysimeters (2 varieties and 9 replicated plots) were installed in July 2015 before the tea transplantation. Each lysimeter was 1.5 m in length, 1.0 m in width and 1.0 m in depth, respectively, and installed in the middle of the row with 3 tea plants in the covered area. The sketches of the installation of the lysimeter and collection of the leachate are shown in Figure S2a,b.
The leachate below 1.0 m was collected in the 50 L nylon tank below the lysimeter. In case of side seepage of soil solution, each lysimeter was surrounded by a piece of 5-mm-thick plastic film (high-density polyethylene, HDPE). Each lysimeter had two 1.5 m pipes, one of which was for air flow and the other one was to pump out the leachate in the 50 L nylon tank. The leachate was pumped out after the rainfall events by a mechanical vacuum pump at 25–30 kPa (PHDmate Co., Ltd, Suzhou, China). The leachate volume was determined gravimetrically and subsamples were collected for further analysis.
An automatic weather station (SP200, LSI Lastem, Milano, Italy) located ~100 m away from the plots provided the weather parameters, e.g., temperature, humidity, rainfall, wind and solar radiation. In addition, three tube access probes (JZSX Corp., Beijing, China) based on time domain reflectometry [28] were elaborately installed in the field to measure soil volumetric water content at 20 cm depth intervals until 1 m depth (Figure S2c).
Soil particle size was tested in pastes of 1:5 (w/v) in 0.2 M sodium hydroxide solution with Mastersizer 3000 Laser particle size analyzer (Malvern Ltd., Worcestershire, UK). Soil pH was measured in pastes of 1:2.5 (w/v) in deionized water with an ORION 3 STAR pH meter (ThermoFisher Ltd., Swedesboro, NJ, USA). Soil bulk density was determined by excavation method by dividing the mass of soil excavated by the volume of water needed to fill the hole [29]. Soil organic carbon determination was conducted via spectrophotometric absorption measurement after organic matter oxidation with potassium dichromate in concentrated sulfuric acid [30].

2.3. Calculations and Statistical Analysis

With the obtained daily weather parameters, reference evapotranspiration (ET0) was calculated according to FAO-56 Penman–Monteith equation. Using the daily weather, the calculation of ET0 can be simplified as follows [14]:
ET 0 = 0.408 Δ ( R n G ) + r 900 T + 273 u 2 ( e s e a ) Δ + r ( 1 + 0.34 u 2 )
where ET0 is the reference evapotranspiration (mm day−1), Rn is the net radiation at the crop surface (MJ m−2 day−1), G is the soil heat flux density (MJ m−2 day−1), T is the mean daily air temperature at 2 m height (°C), u2 is the wind speed at 2 m height (m s−1), es is the saturation vapor pressure (kPa), ea is the actual vapor pressure (kPa), (es − ea) is the saturation vapor pressure deficit (kPa), ∆ is the slope vapor pressure curve (kPa °C−1), γ psychrometric constant (kPa °C−1), and 900 is the conversion factor.
The water balance method considered rainfall and irrigation as the primary water inputs, and subtracted the stored water content at the root zone and deep percolated water below the root zone to calculate ETc [31,32]. As rainfall was the only water input in this study, the daily ETc was estimated by the following equation [33]:
ET c = ( P D ± i = 1 n ( θ 1 θ 2 ) Δ S i ) 1 Δ t
where P is the rainfall (mm), D is the deep percolation (mm), n is the number of soil layers, ∆S is the thickness of each soil layer (mm), θ1 and θ2 are the volumetric soil water content at times one and two, respectively (m3 m−3), and ∆t is the time interval between two consecutive sampling times (days).
Kc was computed using the estimated ETc derived from water balance equation divided by ET0 calculated via Penman–Monteith method through the following equation:
K c = ET c ET 0
where ETc and ET0 are crop evapotranspiration and reference crop evapotranspiration at different sampling intervals in this study, respectively.
F test using one-way ANOVA tested the effect of the tea variety on the tea bush structure, in terms of plant height, plant width, leaf thickness and leaf area index (LAI), and the significance of the means were calculated using Duncan’s multiple range tests at significance level of 0.05 [34]. To ascertain the statistical relationship between rainfall and drainage, we used the results of cumulative rainfall and drainage at every single sampling interval. The linear plus platform model was applied to analyze the relationship between rainfall and drainage. The model can be expressed by the following equations:
y = bx + a1, x < C
y = a2, x ≥ C
where a1, a2 and b are the parameters, and C is the critical value for rainfall.
Data processing and plotting were performed using SPSS22 (SPSS Inc., New York, NY, USA) and Sigmaplot12.5 (Systat Software Inc., Milpitas, CA, USA).

3. Results

3.1. Reference Evapotranspiration (ET0)

The daily ET0, together with the primary meteorological variables, including relative humidity, daily rainfall, wind speed, and air temperature during the measurement period from June 2018 to October 2020, are presented in Figure 1. The air temperature appeared to be relatively low in winter, particularly in December and January, with a mean value of 7.16 °C; whereas, the air temperature in summer was high, with the highest mean value of 28.53 °C during July to August. Similarly, the daily ET0 showed low values in winter, with a mean value of 0.80~0.96 mm day−1, while it reached the peak during summer, with a mean value of 4.84 mm day−1 in August, during the three years.

3.2. Precipitation and Drainage

The time series of the cumulative rainfall and drainage distribution during the study period are illustrated in Figure 2. From June to December 2018, the total rainfall amount was 604 mm, 87% of which was in July to September. The total annual rainfall amount in 2019 was 1565 mm, 42% of which was evenly distributed during February to June, but more frequent and extreme rainfall events happened during July to December. During January to September 2020, the total rainfall was 1985 mm, 73% of which was in the period of May to September.
During the 28 months of the field experiment, there were 39 leachate collection events in total. In the period of 2018, the total leachates were 144 mm and 128 mm for BY1 and LJ43, respectively. Due to the longer measurement period in 2019, the annual total leachate was 252 mm for BY1 and 261 mm for LJ43. However, it decreased to 241 mm and 250 mm in 2020, for BY1 and LJ43, respectively.
In this study, it was found that there was an inflection point in the relationship between the cumulative rainfall and the collected leachate amount (Figure 3). It showed that the collected leachate amount remained at 20.63 mm and 21.45 mm when the cumulative rainfall increased to 78.02 mm and 90.98 mm, for BY1 and LJ43, respectively, despite the increasing cumulative rainfall.

3.3. Changes of ETc and Kc during the Experiment Period

The trends of the average ETc and adjusted ETc (calculated using critical rainfall as the water input item when the cumulative rainfall increased to and above the critical value) during the period from March 2019 to September 2020 are illustrated in Figure 4. The adjusted ETc values of BY1 and LJ43 revealed a minor fluctuation when compared to the ETc during the period of this experiment. Moreover, for both BY1 and LJ43, the summer season always exhibited a higher daily ETc and adjusted ETc, while the winter season revealed the reverse result. For BY1, the maximum ETc value was 17.63 mm day−1, while the maximum adjusted ETc value was 4.83 mm day−1, and the minimum ETc value was 0.56 mm day−1, whereas the minimum adjusted ETc value was 0.52 mm day−1. For LJ43, the maximum ETc value was 17.50 mm day−1, whereas the maximum adjusted ETc value was 4.71 mm day−1, and the minimum ETc value was 0.76 mm day−1, while the minimum adjusted ETc value was 0.75 mm day−1.
Regarding the Kc value and adjusted Kc value (ratio of adjusted ETc to ET0), for both BY1 and LJ43, the adjusted Kc value had tiny fluctuations, while the Kc value waved sharply during the whole experiment period (Figure 5). The Kc values of BY1 and LJ43 ranged from 0.52 to 4.41 (mean = 1.61) and from 0.55 to 4.39 (mean = 1.62), respectively, whereas the adjusted Kc values of BY1 and LJ43 ranged from 0.43 to 1.02 (mean = 0.71) and from 0.48 to 1.22 (mean = 0.84), respectively. Unlike the ETc value and adjusted ETc value, the maximum Kc values of BY1 and LJ43 were present in the autumn season (September 2019), while the maximum adjusted Kc values of BY1 and LJ43 appeared in about summer time (June 2020).

3.4. Monthly Variation in ET0, ETc and Adjusted Kc

According to the sampling date and interval, each sampling date was separated to calculate the monthly ET0, ETc, and adjusted Kc values (Table 2). The values of ET0 increased from 2.35 mm day−1 in March 2019 to its maximum value of 4.72 mm day−1 in August 2019, and then decreased to its minimum value of 0.99 mm day−1 in January 2020. A similar change in monthly ET0 was found in 2020.
For both BY1 and LJ43, the high ETc values often occurred in summer time, especially in July and August, both for BY1 and LJ43, while the low ETc values often appeared in winter time. The maximum and minimum values of monthly ETc were 3.67 mm day−1 and 0.52 mm day−1 for BY1, and 4.55 mm day−1 and 0.75 mm day−1 for LJ43, respectively.
The ranges of monthly adjusted Kc values were 0.43–0.92 mm day−1 for BY1 and 0.55–1.15 mm day−1 for LJ43. The maximum adjusted Kc values for both BY1 and LJ43 occurred in June 2020, while the minimum values appeared in June 2019 for BY1 and in June 2019 for LJ43. Overall, the monthly total ET0 during the experiment was 1034 mm, and the ETc values for BY1 and LJ43 were 711 mm and 830 mm, respectively. The average monthly adjusted Kc values were 0.68 for BY1 and 0.82 for LJ43.

4. Discussion

In our study, the pattern of the cumulative rainfall–drainage curve fitted well with the linear plus platform model (Figure 3), indicating a critical amount rainfall that led to maximum drainage. This result may be attributed to the strong surface runoff in big rainfall events, because rainwater firstly arrives at the soil surface, and then runoff occurs if the soil surface has a slope. In addition, strong rainfall usually results in surface runoff happening faster than in rainwater infiltrate soils [35]. Previous research has proven that rainfall intensity and kinetic energy are key factors for determining the formation of runoff, and the rate of crust formation and groundwater infiltration [36,37,38]. Heavy rainfall intensity could shorten the formation time of runoff and sharply decrease the infiltrate curve [39,40]. However, previous studies also revealed that rainfall water prefers to infiltrate into the soil when the rainfall event is initiated, and finally forms the runoff when the soil water-adsorption capacity is saturated [41,42,43]. Thus, determining the saturation of soil water-adsorption capacity is vital for us to monitor the subsequent water movement in the field. In addition, an estimation of the runoff amount should be taken into account when applying the water balance equation to estimate drainage or evapotranspiration, which could improve the accuracy of estimating the drainage. Given that the runoff amount could not be estimated in our study, and considering that runoff only occurs in heavy rainfall situations, a linear plus platform model was applied to determine the relationship between cumulative rainfall and drainage, and the critical amount of rainfall was successfully estimated in our study. Consequently, we suggested that using the linear plus platform model to estimate the critical cumulative rainfall could be regarded as an alternative for this estimation, especially because runoff monitoring is absent in big rainfall events.
In the present study, the mean Kc values of BY1 (0.71) and LJ43 (0.84) were slightly lower than the recommended value (0.9–1.0) for tea plants, according to the FAO Report [14], but were comparable to the mean value of 0.72 reported by Dabral for tea in Assam at different irrigation levels [44]. This result may be due to the young age of the tea plant, since plant growth stage is one of the most important factors for Kc value, and lower Kc values represent slower plant growth and lower canopy cover [45]. Similar results were found in citrus and grapevine orchard; for example, Alves et al. reported that the Kc value of citrus orchard increased from 0.57 in a one-year-old plant to 1.12 in a four-year-old plant [46]. Similarly, the Kc value of grapevine orchard increased from 0.4 (young shoot germination stage) to 1.0 (flourishing period) under freshwater irrigation [47].
As for the seasonal variation in Kc values (Table 2), it seemed that a higher Kc value often appeared during July to September, while lower values occurred during October to January. This phenomenon was in line with the trends of precipitation distribution and air temperature during the study period (Figure 1 and Figure 2), suggesting that rainfall and air temperature could be responsible for the results. A previous study showed that the variation in Kc values between seasons could be attributed to several reasons, such as the trend of canopy development, rainfall amounts in different seasons, and differences in energy absorption characteristics [4]. When plants received adequate water, together with increased temperature, this resulted in a rapid growth rate and peak of ETc and Kc [4]. These findings were consistent with our results that the maximum Kc value appears in July to November. In addition, a relatively lower Kc value (mean = 0.71) was found in November to April; this may be due to the limited rainfall. This is because the limited rainfall caused soil moisture stress, and then resulted in declining growth rate, which, in turn, decreased the ETc and Kc values. Moreover, Squire et al. reported that tea plant growth rate became very slow when the temperature was below 12–13 °C [48]. This is a good explanation of the higher Kc in July to September and the lower Kc values in October to January, because the mean temperature was 28.53 °C during July to August and 7.16 °C during December to January.
During the experimental period, we found that heavy tea tree pruning decreased the ETc value and adjusted Kc value of BY1 and LJ43 (Figure 4 and Figure 5). This result was similar to the previous findings reported by Lv et al. [19], which were that the mean daily ET rates could decrease by 36.7–48.3% after pruning in a subtropical tea plantation, by using the static chamber/infrared gas analyzer (IRGA). Similarly, Geng et al. also found that strong pruning interference resulted in the highest differences between ETa and ET0 in the pruning period, compared to other seasons with relatively low mean monthly ETa/ET0 values (average 0.76) [20]. The reason for this phenomenon may be related to the heavy pruning of tea resulting in a lower leaf area index (LAI), which, in turn, caused a decrease in the means of water transport, causing a decline in the ETc and Kc values [20]. Moreover, from the dual crop coefficient perspective, Kato and Kamichika [49] found that the highly LAI-dependent basal crop coefficient (Kcb), which represents crop transpiration, exponentially increased with the LAI, while the soil water evaporation coefficient (Ke) showed a positive relationship with soil water content and a negative linear relationship with the LAI. Thus, intense tea pruning, on the one hand, lessened the LAI of tea bushes, which subsequently caused a reduction in the Kcb value. On the other hand, a large quantity of pruned litter was discarded and covered in the tea garden’s soil, which substantially curbed the soil water evaporation and the Ke value. Therefore, we suggested that heavy tea pruning in May weakened the evapotranspiration capacity of the tea plantations.
Interestingly, our results found that the ETc value and Kc value of BY1 were consistently lower than those of LJ43 (Table 2). This could be attributed to the differences in the plant growth status between BY1 and LJ43; the plant growth status can be easily reflected by the canopy structure, such as plant height, canopy width, leaf thickness, and LAI. The tea plant height, width, leaf thickness, and LAI of LJ43 were significantly greater than BY1 (Table 3). Therefore, LJ43 has a relatively well-structured canopy, which indicated that LJ43 had higher canopy cover and a greater LAI than BY1, resulting in higher values of ETc and Kc. These findings were in good agreement with the former study on pomegranate crops [14]. In addition, previous studies also revealed that the Kc value has a strong relationship with LAI and soil coverage percentage in many fruit plantations, such as maize, bean, apple, and pecan [50,51,52,53]. Therefore, we suggested that the significantly higher LAI of LJ43, compared to BY1, could be another reason for the lower values of ETc and Kc of BY1.
In our study, the experiment was conducted in a rain-fed subtropical region with uneven rainfall distribution during the study period. Therefore, the interval of drainage sampling from the lysimeter could not be scheduled uniformly, due to a lack of water input for certain periods, which occurred more often from November to January. Thus, the absence of data for the ETc and Kc calculations during this period made it difficult to precisely reveal the monthly variation in water use status in the tea plantation, which needs to be further studied by alternative approaches for the monitoring and estimation of ETc and Kc. In addition, as RS- and GIS-based remote sensing approaches and algorithms have become important tools for ETc and Kc estimation at a regional scale, we suggest that further studies should enhance the remote sensing measurements of the energy and water fluxes in tea plantations; for example, improving the monitoring accuracy via the integration of innovated methods to obtain the remote sensing data and meteorological variables, and subsequently using a readily available, open-source model (such as SEBALIGEE) to calibrate it with the ground-based lysimeter could produce larger spatiotemporal assessments.

5. Conclusions

In the present study, the linear plus platform model was applied to analyze the relationship between cumulative rainfall and lysimeter drainage, and critical cumulative rainfall values (78.02 mm for BY1 and 90.98 mm for LJ43) were obtained to estimate ETc and Kc via the water balance equation when the cumulative rainfall was more than the critical value. Finally, ETc and adjusted ETc, and the corresponding Kc and adjusted Kc, were estimated by using the original cumulative rainfall. The Kc values for BY1 and LJ43 ranged from 0.52 to 4.41 and from 0.55 to 4.39, respectively, whereas the adjusted Kc values of BY1 and LJ43 ranged from 0.43 to 1.02 and from 0.48 to 1.22, respectively. The adjusted ETc and Kc values revealed a minor fluctuation when compared to the ETc and Kc values during the experimental period. In addition, periodic pruning significantly decreased the LAI of the tea plant, and resulted in a reduction in ETc and Kc values. Furthermore, plant growth difference is a key factor contributing to the change in ETc and Kc values between different tea cultivars. Collectively, our study provided a reliable reference of Kc in a tea plantation, and illuminated the effects of pruning and plant growth difference on Kc.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/agronomy11112332/s1, Figure S1: (a) Original images of BY1 and LJ43, (b) Illustration of experiment layout, lysimeter installation and soil volumetric water content probe site, and (c) Tea bushes on the study site, Figure S2: (a) Illustrations of lysimeter and leachate collecting from lysimeter, (b) Top view of lysimeter in-stalled, (c) Illustrations of the weather station and soil volumetric water content probe.

Author Contributions

Conceptualization, S.Z. and J.R.; writing—original draft preparation, S.Z.; writing—review and editing, K.N., L.J. and J.R.; formal analysis, C.Z. and H.C.; resources, X.Y. and W.H.; funding acquisition, J.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Earmarked Fund for China Agriculture Research System (CARS 19); the Agricultural Science and Technology Innovation Program of the Chinese Academy of Agricultural Sciences (CAAS-ASTIP-2018-TRICAAS, CAAS-XTCX2016015).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hasimoto, M. The origin of the tea plant. In Proceedings of the 2001 International Conference on O-CHA (Tea) Culture and Science, Shizuoka, Japan, 5–8 October 2001; Session I. pp. J5–J7. [Google Scholar]
  2. Piyashee, M.; Tuhin, G. Impact of Climate on Tea Production: A Study of the Dooars Region in India. Res. Square 2021, 1–21. [Google Scholar] [CrossRef]
  3. Chen, X.H.; Zhuang, C.G.; He, Y.F.; Wang, L.; Han, G.Q.; Chen, C. Photosynthesis, yield, and chemical composition of Tieguanyin tea plants (Camellia sinensis (L.) O. Kuntze) in response to irrigation treatments. Agric. Water Manag. 2010, 97, 419–425. [Google Scholar] [CrossRef]
  4. Sikka, A.; Madegowda, M.; Sahoo, D.C. Determination of crop evapotranspiration of tea (Thea sinensis) using weighing lysimeter for the Nilgiris. J. Agrometeorol. 2009, 11, 144–147. [Google Scholar]
  5. Ruiming, F.; Shijie, S. Daily reference evapotranspiration prediction of Tieguanyin tea plants based on mathematical morphology clustering and improved generalized regression neural network. Agric. Water Manag. 2020, 236, 106177. [Google Scholar] [CrossRef]
  6. Ma, J.; Jia, X.; Zha, T.; Bourque, C.P.A.; Tian, Y.; Bai, Y.; Zhou, C. Ecosystem water use efficiency in a young plantation in Northern China and its relationship to drought. Agric. For. Meteorol. 2019, 275, 1–10. [Google Scholar] [CrossRef]
  7. Ahmed, S.; Griffin, T.; Cash, S.B.; Han, W.Y.; Matyas, C.; Long, C.; Xue, D. Global Climate Change, Ecological Stress, and Tea Production. In Stress Physiology of Tea in the Face of Climate Change; Springer: Singapore, 2018; pp. 1–23. [Google Scholar] [CrossRef]
  8. Ahmed, S.; Orians, C.M.; Griffin, T.S.; Buckley, S.; Unachukwu, U.; Stratton, A.E.; Kennelly, E.J. Effects of water availability and pest pressures on tea (Camellia sinensis) growth and functional quality. AoB Plants 2014, 6, plt054. [Google Scholar] [CrossRef] [Green Version]
  9. Ahmed, S.; Stepp, J.R.; Orians, C.; Griffin, T.; Matyas, C.; Robbat, A.; Kennelly, E. Effects of Extreme Climate Events on Tea (Camellia sinensis) Functional Quality Validate Indigenous Farmer Knowledge and Sensory Preferences in Tropical China. PLoS ONE 2014, 9, e109126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Wijeratne, M.A. Vulnerability of Sri Lanka tea production to global climate change. Water Air Soil Pollut. 1996, 92, 87–94. [Google Scholar] [CrossRef]
  11. Chen, B.; Yang, J.; Pu, J. Statistical Characteristics of Raindrop Size Distribution in the Meiyu Season Observed in Eastern China. J. Meteorol. Soc. Jpn. Ser. II 2013, 91, 215–227. [Google Scholar] [CrossRef] [Green Version]
  12. Samuel, A.; Girma, A.; Zenebe, A.; Ghebreyohannes, T. Spatio-temporal variability of evapotranspiration and crop water requirement from space. J. Hydrol. 2018, 567, 732–742. [Google Scholar] [CrossRef]
  13. Payero, J.O.; Irmak, S. Daily energy fluxes, evapotranspiration and crop coefficient of soybean. Agric. Water Manag. 2013, 129, 31–43. [Google Scholar] [CrossRef]
  14. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; pp. 2–15. [Google Scholar]
  15. Bhantana, P.; Lazarovitch, N. Evapotranspiration, crop coefficient and growth of two young pomegranate (Punica granatum L.) varieties under salt stress. Agric. Water Manag. 2010, 97, 715–722. [Google Scholar] [CrossRef]
  16. Jiang, X.; Kang, S.; Tong, L.; Li, F.; Li, D.; Ding, R.; Qiu, R. Crop coefficient and evapotranspiration of grain maize modified by planting density in an arid region of northwest China. Agric. Water Manag. 2014, 142, 135–143. [Google Scholar] [CrossRef]
  17. Kang, S.; Gu, B.; Du, T.; Zhang, J. 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]
  18. Buttar, N.; Hu, Y.; Shabbir, A.; Lakhiar, I.; Ullah, I.; Ali, A. Estimation of Evapotranspiration Using Bowen Ratio Method. IFAC-PapersOnLine 2018, 51, 807–810. [Google Scholar] [CrossRef]
  19. Lv, W.; Yang, G.; Wan, R. Daily variations of evapotranspiration rates of different ages tea plantation in Yixing city based on static chamber/igra. Res. Environ. Yangtze Basin 2012, 21, 1370. (In Chinese) [Google Scholar]
  20. Geng, J.; Li, H.; Pang, J.; Zhang, W.; Chen, D. Dynamics and environmental controls of energy exchange and evapotranspiration in a hilly tea plantation, China. Agric. Water Manag. 2020, 241, 106364. [Google Scholar] [CrossRef]
  21. Allen, R.G.; Masahiro, T.; Ricardo, T. Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC)—Model. J. Irrig. Drain. Eng. 2007, 133, 380–394. [Google Scholar] [CrossRef]
  22. Anderson, M.C.; Kustas, W.P.; Norman, J.M.; Hain, C.R.; Mecikalski, J.R.; Schultz, L.; González-Dugo, M.P.; Cammalleri, C.; d’Urso, G.; Pimstein, A.; et al. Mapping daily evapotranspiration at field to continental scales using geostationary and polar orbiting satellite imagery. Hydrol. Earth Syst. Sci. 2011, 15, 223–239. [Google Scholar] [CrossRef] [Green Version]
  23. Allam, M.; Mhawej, M.; Meng, Q.; Faour, G.; Abunnasr, Y.; Fadel, A.; Xinli, H. Monthly 10-m evapotranspiration rates re trieved by SEBALI with Sentinel-2 and MODIS LST data. Agric. Water Manag. 2021, 243, 106432. [Google Scholar] [CrossRef]
  24. Mhawej, M.; Nasrallah, A.; Abunnasr, Y.; Fadel, A.; Faour, G. Better irrigation management using the satellite-based adjusted single crop coefficient (aKc) for over sixty crop types in California, USA. Agric. Water Manag. 2021, 256, 107059. [Google Scholar] [CrossRef]
  25. Dagg, M. A study of the water use of tea in East Africa using a hydraulic lysimeter. Agric. Meteorol. 1970, 7, 303–320. [Google Scholar] [CrossRef]
  26. Sikka, K.; Sahoo, D.C.; Madhu, M.; Selvi, V. Determination of Crop Coefficient of Tea. J. Agric. Eng. 2009, 46, 41–45. [Google Scholar]
  27. Lv, W.; Yang, G.; Wan, R.R.; Li, W.L. Comparative analysis of evapotranspiration rates of tea crops before and after pruning in the Western Hills of the Tai Lake Basin. Chin. J. Eco-Agric. 2013, 21, 184–191. [Google Scholar] [CrossRef]
  28. Villoro, A.; Latorre, B.; Tormo, J.; José Jiménez, J.; Victoria López, M.; Manuel Nicolau, J.; Moret-Fernández, D. A TDR wireless device for volumetric water content sensing. Comput. Electron. Agric. 2021, 181, 105939. [Google Scholar] [CrossRef]
  29. Page-Dumroese, D.S.; Brown, R.E.; Jurgensen, M.F.; Mroz, G.D. Comparison of Methods for Determining Bulk Densities of Rocky Forest Soils. Soil Sci. Soc. Am. J. 1999, 63, 379–383. [Google Scholar] [CrossRef]
  30. Walinga, I.; Kithome, M.; Novozamsky, I.; Houba, V.J.G.; van der Lee, J.J. Spectrophotometric determination of organic carbon in soil. Commun. Soil Sci. Plant Anal. 1992, 23, 1935–1944. [Google Scholar] [CrossRef]
  31. Sepaskhah, A.R.; Andam, M. Crop coefficient of sesame in a semi-arid region of I.R. Iran. Agric. Water Manag. 2001, 49, 51–63. [Google Scholar] [CrossRef]
  32. Flumignan, D.L.; de Faria, R.T.; Prete, C.E.C. Evapotranspiration components and dual crop coefficients of coffee trees during crop production. Agric. Water Manag. 2011, 98, 791–800. [Google Scholar] [CrossRef]
  33. Jensen, M.E.; Burman, R.D.; Allen, R.G. Evapotranspiration and Irrigation Water Requirements; ASCE Manuals and Reports on Engineering Practice No. 70; American Society of Civil Engineers: New York, NY, USA, 1990; p. 332. [Google Scholar]
  34. Zotarelli, L.; Scholberg, J.M.; Dukes, M.D.; Muñoz-Carpena, R. Monitoring of Nitrate Leaching in Sandy Soils. J. Environ. Qual. 2007, 36, 953–962. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Dexter, A.R. Advances in characterization of soil structure. Soil Tillage Res. 1988, 11, 199–238. [Google Scholar] [CrossRef]
  36. Wang, H.; Gao, J.E.; Zhang, M.; Li, X.; Zhang, S.; Jia, L. Effects of rainfall intensity on groundwater recharge based on simulated rainfall experiments and a groundwater flow model. CATENA 2015, 127, 80–91. [Google Scholar] [CrossRef]
  37. Liu, G.; Zheng, F.; Jia, L.; Jia, Y.; Zhang, X.; Hu, F.; Zhang, J. Interactive effects of raindrop impact and groundwater seepage on soil erosion. J. Hydrol. 2019, 578, 124066. [Google Scholar] [CrossRef]
  38. Alizadehtazi, B.; Gurian, P.L.; Montalto, F.A. Impact of successive rainfall events on the dynamic relationship between vegetation canopies, infiltration, and recharge in engineered urban green infrastructure systems. Ecohydrology 2020, 13, e2185. [Google Scholar] [CrossRef]
  39. Cuomo, S.; Della Sala, M. Rainfall-induced infiltration, runoff and failure in steep unsaturated shallow soil deposits. Eng. Geol. 2013, 162, 118–127. [Google Scholar] [CrossRef]
  40. Assouline, S.; Ben-Hur, M. Effects of rainfall intensity and slope gradient on the dynamics of interrill erosion during soil surface sealing. CATENA 2006, 66, 211–220. [Google Scholar] [CrossRef]
  41. Rao, K.P.; Steenhuis, T.; Cogle, A.; Srinivasan, S.; Yule, D.; Smith, G. Rainfall infiltration and runoff from an Alfisol in semi-arid tropical India. II. No-till systems. Soil Tillage Res. 1998, 48, 61–69. [Google Scholar] [CrossRef] [Green Version]
  42. Wang, J.; Chen, L.; Yu, Z. Modeling rainfall infiltration on hillslopes using Flux-concentration relation and time compression approximation. J. Hydrol. 2018, 557, 243–253. [Google Scholar] [CrossRef]
  43. Hueso-González, P.; Martínez-Murillo, J.F.; Ruiz-Sinoga, J.D. Benefits of adding forestry clearance residues for the soil and vegetation of a Mediterranean mountain forest. Sci. Total Environ. 2018, 615, 796–804. [Google Scholar] [CrossRef]
  44. Dabral, P.P.; Rao, K.A. Estimation of crop coefficient and irrigation requirement under various irrigation levels for tea during dry period. Indian J. Soil Conserv. 1997, 25, 233–235. [Google Scholar]
  45. Wang, T.; Melton, F.S.; Pôças, I.; Johnson, L.F.; Thao, T.; Post, K.; Cassel-Sharma, F. Evaluation of crop coefficient and evapotranspiration data for sugar beets from landsat surface reflectances using micrometeorological measurements and weighing lysimetry. Agric. Water Manag. 2021, 244, 106533. [Google Scholar] [CrossRef]
  46. Alves, J., Jr.; Folegatti, M.V.; Parsons, L.R.; Bandaranayake, W.; da Silva, C.R.; da Silva, T.J.A.; Campeche, L.F.S.M. Determination of the crop coefficient for grafted ‘Tahiti’ lime trees and soil evaporation coefficient of rhodic kandiudalf clay soil in Sao Paulo, Brazil. Irrig. Sci. 2007, 25, 419–428. [Google Scholar] [CrossRef]
  47. Netzer, Y.; Yao, C.; Shenker, M.; Bravdo, B.A.; Schwartz, A. Water use and the development of seasonal crop coefficients for Superior Seedless grapevines trained to an open-gable trellis system. Irrig. Sci. 2009, 27, 109–120. [Google Scholar] [CrossRef]
  48. Squire, G.R. Weather, physiology and seasonality of tea (Camellia sinensis L.) yields of Malawi. Exp. Agric. 1979, 15, 321–330. [Google Scholar] [CrossRef]
  49. Kato, T.; Kamichika, M. Determination of a crop coefficient for evapotranspiration in a sparse sorghum field. Irrig. Drain. 2006, 55, 165–175. [Google Scholar] [CrossRef]
  50. Shao, G.; Han, W.; Zhang, H.; Liu, S.; Wang, Y.; Zhang, L.; Cui, X. Mapping maize crop coefficient Kc using random forest algorithm based on leaf area index and UAV-based multispectral vegetation indices. Agric. Water Manag. 2021, 252, 106906. [Google Scholar] [CrossRef]
  51. De Medeiros, G.A.; Arruda, F.B.; Sakai, E.; Fujiwara, M. The influence of crop canopy on evapotranspiration and crop coefficient of bean (Phaseolus vulgaris L.). Agric. Water Manag. 2001, 49, 211–224. [Google Scholar] [CrossRef]
  52. Gong, D.; Kang, S.; Yao, L.; Zhang, L. Estimation of evapotranspiration and its components from an apple orchard in northwest China using sap flow and water balance methods. Hydrol. Process. 2007, 21, 931–938. [Google Scholar] [CrossRef]
  53. Wang, J.; Sammis, T.W.; Andales, A.A.; Simmons, L.J.; Gutschick, V.P.; Miller, D.R. Crop coefficient of open-canopy pecan orchards. Agric. Water Manag. 2007, 88, 253–262. [Google Scholar] [CrossRef]
Figure 1. Daily relative humidity, rainfall, wind speed, mean air temperature and reference evapotranspiration (ET0) during the experiment period from June 2018 to October 2020 on the research site.
Figure 1. Daily relative humidity, rainfall, wind speed, mean air temperature and reference evapotranspiration (ET0) during the experiment period from June 2018 to October 2020 on the research site.
Agronomy 11 02332 g001
Figure 2. Cumulative rainfall and drainage distribution during the study period from June 2018 to September 2020. Drainage entries represent the mean values in different N treatments. These values are mean for all 9 barrels. Inserted error bars denote standard error of the mean.
Figure 2. Cumulative rainfall and drainage distribution during the study period from June 2018 to September 2020. Drainage entries represent the mean values in different N treatments. These values are mean for all 9 barrels. Inserted error bars denote standard error of the mean.
Agronomy 11 02332 g002
Figure 3. Correlation and change point between rainfall and drainage.
Figure 3. Correlation and change point between rainfall and drainage.
Agronomy 11 02332 g003
Figure 4. Variations in crop evapotranspiration (ETc) and adjusted ETc for BY1 and LJ43 during the period from March 2019 to September 2020. These values are mean of 9 lysimeters. Inserted error bars denote standard error of the mean. Vertical gray arrows denote tea pruning practice.
Figure 4. Variations in crop evapotranspiration (ETc) and adjusted ETc for BY1 and LJ43 during the period from March 2019 to September 2020. These values are mean of 9 lysimeters. Inserted error bars denote standard error of the mean. Vertical gray arrows denote tea pruning practice.
Agronomy 11 02332 g004
Figure 5. Variations in crop coefficient (Kc) and adjusted Kc for BY1 and LJ43 during the period from March 2019 to September 2020. These values are mean of 9 lysimeter barrels. Inserted error bars denote standard error of the mean. Vertical gray arrows denote tea pruning practice.
Figure 5. Variations in crop coefficient (Kc) and adjusted Kc for BY1 and LJ43 during the period from March 2019 to September 2020. These values are mean of 9 lysimeter barrels. Inserted error bars denote standard error of the mean. Vertical gray arrows denote tea pruning practice.
Agronomy 11 02332 g005
Table 1. Soil properties before the experiment.
Table 1. Soil properties before the experiment.
Depth (cm)Particle Size Distribution (%)pHOrganic C
(mg g−1)
Bulk Density
(g cm−3)
SandSiltClay
0–201.0728.9469.994.475.711.17
20–401.9131.0667.034.665.471.12
40–602.5830.9966.434.404.091.13
60–803.1434.3262.534.375.021.19
80–1001.4632.8265.734.344.771.16
Table 2. Monthly reference evapotranspiration ET0, crop evapotranspiration ETc and adjusted crop coefficient Kc both for BY1 and LJ43 during the period from March 2019 to September 2020. Data are means and standard errors for 9 lysimeter barrels.
Table 2. Monthly reference evapotranspiration ET0, crop evapotranspiration ETc and adjusted crop coefficient Kc both for BY1 and LJ43 during the period from March 2019 to September 2020. Data are means and standard errors for 9 lysimeter barrels.
MonthsET0 (mm Day−1)ETc (mm Day−1)Adjusted Kc
BY1LJ43BY1LJ43
March 20192.35 ± 0.021.59 ± 0.071.71 ± 0.080.68 ± 0.030.73 ± 0.03
April 20193.47 ± 0.052.23 ± 0.142.30 ± 0.130.65 ± 0.050.68 ± 0.05
May 20194.35 ± 0.092.36 ± 0.042.77 ± 0.050.54 ± 0.010.64 ± 0.02
June 20193.66 ± 0.001.58 ± 0.022.01 ± 0.010.43 ± 0.010.55 ± 0.00
July 20194.46 ± 0.233.59 ± 0.264.06 ± 0.160.79 ± 0.020.92 ± 0.01
August 20194.72 ± 0.033.66 ± 0.313.89 ± 0.140.77 ± 0.060.82 ± 0.03
September 20193.23 ± 0.062.42 ± 0.212.95 ± 0.090.73 ± 0.050.91 ± 0.01
October 20192.78 ± 0.001.20 ± 0.052.05 ± 0.060.43 ± 0.020.74 ± 0.02
November/December 20191.19 ± 0.000.52 ± 0.010.75 ± 0.000.44 ± 0.010.63 ± 0.00
January 20200.99 ± 0.000.56 ± 0.050.76 ± 0.070.56 ± 0.050.76 ± 0.07
Feburary/March 20201.34 ± 0.000.97 ± 0.021.17 ± 0.000.72 ± 0.010.87 ± 0.00
April 20202.96 ± 0.152.30 ± 0.042.59 ± 0.060.82 ± 0.050.93 ± 0.07
May 20204.16 ± 0.053.42 ± 0.103.57 ± 0.120.83 ± 0.040.87 ± 0.04
June 20203.98 ± 0.033.67 ± 0.074.55 ± 0.040.92 ± 0.021.15 ± 0.02
July 20203.46 ± 0.063.02 ± 0.053.43 ± 0.130.87 ± 0.031.00 ± 0.05
August 20204.64 ± 0.063.42 ± 0.124.11 ± 0.090.73 ± 0.020.88 ± 0.01
September 20202.84 ± 0.001.92 ± 0.032.31 ± 0.010.68 ± 0.010.81 ± 0.00
Monthly total1033.63 ± 18.94701.60 ± 29.67830.24 ± 26.760.68 ± 0.0370.82 ± 0.036
Monthly mean3.21 ± 0.292.26 ± 0.262.65 ± 0.290.68 ± 0.040.82 ± 0.04
Table 3. Comparison of tea stands’ properties between three different tea cultivars. Data are means and standard errors of 9 plots of tea stands. Different letters indicate significance (p < 0.05).
Table 3. Comparison of tea stands’ properties between three different tea cultivars. Data are means and standard errors of 9 plots of tea stands. Different letters indicate significance (p < 0.05).
Tea CultivarsPlant Height (cm)Plant Width (cm)Leaf Thickness (cm)LAI (m2 m−2)
DG94.53 ± 1.22 a111.58 ± 1.43 a41.22 ± 0.63 a8.65 ± 0.13 a
LJ4388.83 ± 1.01 b110.14 ± 1.23 a39.94 ± 0.48 a8.72 ± 0.17 a
BY177.22 ± 1.03 c95.08 ± 1.47 b33.89 ± 0.70 b7.64 ± 0.12 b
Note: DG stands for Dangui, another tea cultivar besides BY1 and LJ43.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zheng, S.; Ni, K.; Ji, L.; Zhao, C.; Chai, H.; Yi, X.; He, W.; Ruan, J. Estimation of Evapotranspiration and Crop Coefficient of Rain-Fed Tea Plants under a Subtropical Climate. Agronomy 2021, 11, 2332. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11112332

AMA Style

Zheng S, Ni K, Ji L, Zhao C, Chai H, Yi X, He W, Ruan J. Estimation of Evapotranspiration and Crop Coefficient of Rain-Fed Tea Plants under a Subtropical Climate. Agronomy. 2021; 11(11):2332. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11112332

Chicago/Turabian Style

Zheng, Shenghong, Kang Ni, Lingfei Ji, Chenguang Zhao, Hongling Chai, Xiaoyun Yi, Weizhong He, and Jianyun Ruan. 2021. "Estimation of Evapotranspiration and Crop Coefficient of Rain-Fed Tea Plants under a Subtropical Climate" Agronomy 11, no. 11: 2332. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11112332

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