Next Article in Journal
Estimating the Reducing Power of Carbon Nanotubes and Granular Activated Carbon Using Various Compounds
Previous Article in Journal
Case Study of Effects of Mineral N Fertilization Amounts on Water Productivity in Rainfed Winter Rapeseed Cultivation on a Sandy Soil in Brandenburg (Germany) over Three Years
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relationships between Leaf Area Index and Evapotranspiration and Crop Coefficient of Hilly Apple Orchard in the Loess Plateau

1
College of Natural Resources and Environment, Northwest A&F University, No. 3, Taicheng Road, Yangling, Xianyang 712100, China
2
Key Laboratory of Plant Nutrition and the Agri-Environment in Northwest China, Ministry of Agriculture, No. 3, Taicheng Road, Yangling, Xianyang 712100, China
*
Author to whom correspondence should be addressed.
Water 2021, 13(14), 1957; https://doi.org/10.3390/w13141957
Submission received: 16 June 2021 / Revised: 10 July 2021 / Accepted: 12 July 2021 / Published: 16 July 2021
(This article belongs to the Section Hydrology)

Abstract

:
Drought and water shortage are the key factors that restrict the sustainable development of the apple industry in the Chinese Loess Plateau. The accurate prediction of ET can provide a scientific basis for water management of apple orchards. A study on the relationship between LAI, ET and crop coefficient Kc under water deficit is particularly necessary for the accurate prediction of ET in apple orchards. In this work, the crop coefficient Kc under water deficit was defined as the product of the crop coefficient KcI under no water stress and the water stress coefficient Ks, namely Kc = KcI × Ks. LAI and ET of the hilly apple orchard were measured from April to September in 2019 and 2020. The results showed: (1) The LAI of the apple orchard showed a trend of rapid increase—moderate increase—declined during the growth period, with 0.26–2.16 [m2 m−2] variation range; (2) The ET of the orchard was greater than the rainfall, the maximum ET was in July or August. The maximum components of ET in the apple orchard was E, with 47.8–49.1% of ET; T accounted for 42.5–43.9% of ET; Ic accounted for only 9.1–9.6% of ET; (3) There was a significant exponential relationship between the LAI and T or ET. The crop coefficient KcI under no water stress changed with the development of the apple tree canopy. The variation of water stress Ks was basically consistent with the variation of rainfall; (4) There is a significant exponential relationship between LAI and crop coefficient Kc under water deficit (Kc = 0.1141e1.0665LAI, R2 = 0.7055, p < 0.01). This study demonstrates that LAI could be used to estimate the crop coefficient Kc of apple orchards under water deficit in the Loess Plateau, and the actual evapotranspiration of apple orchards in this region could be predicted.

1. Introduction

The Loess hilly region of northern Shaanxi is a high-quality apple (Malus domestica) production region due to its vast land resources, abundant sunlight, large temperature difference, good ventilation and low pollution, and the apple area has exceeded 1.3 million hm2 by 2020, which increased farmers’ income and reduced soil loss. However, less rainfall (usually only 450–550 mm), large soil evaporation and strong transpiration in this area result in serious soil moisture deficits and widespread deep soil drying in the hilly apple orchard. Soil water deficit has a great negative effect on apple production.
Evapotranspiration (ET) is an important component of water balance in apple orchards, which consists of soil evaporation (E), transpiration (T) and canopy interception (Ic) [1]. The main factors affecting ET are the canopy area and architecture, evaporative power of the atmosphere (ET0) and stomatal conductance and soil type [2]. ET0 is defined as the evapotranspiration rate at which water would be removed from a reference surface where water is not limited or a limited factor [2]; it is necessary to accurately measure/quantify it under varied data-availability conditions. Despite the necessity of taking decisions on sectoral water distributions, irrigation scheduling, groundwater recharge, and reservoir management at both the spatial and temporal scales, in situ measurement of ET is still a difficult task [3,4]. The widely used equations, which have been developed to estimate ET0, include the California Irrigation Management Information System Penman equation, the Penman–Monteith equation standardized by the Food and Agriculture Organization (FAO), the Reference Evapotranspiration equation standardized by the American Society of Civil Engineers, the Hargreaves equation and Priestley-Taylor model [5,6,7,8,9,10], besides, indirect ET estimation methods, such as moderate resolution imaging spectroradiometer (MODIS), satellite-based remote-sensing techniques and the water-budget approach built into the semidistributed variable infiltration capacity (VIC-3L) land-surface model [3]. The FAO’s suggested approach of the FAO-56 Penman-Monteith model (FAO-ET0) is recommended as the standard model world-wide for estimating ET0 [2,5]. The FAO-ET0 require meteorological data from a well-watered, 12-cm-high grassy surface that fully covers the ground to estimate ET0 [11], as the use of inappropriate data for ET0 estimation from non-ideal surfaces, leads to significant and systematic cumulative errors introducing uncertainties when determining the crop water requirements in a region [12]. The major limitation of the FAO-ET0 is the high requirement of several meteorological variables that are often incomplete and/or unavailable [13,14]. The application of the FAO-ET0 equation has been reported to lead to errors, especially when estimated variables as an alternative to measured values were used as inputs in ET0 computations [15,16]. The most important meteorological variables for ET0 estimation are air temperature and solar radiation [17,18], and daily temperature range could be related to relative humidity and cloudiness [19]. In spite of advection, which depends on the interaction between temperature, relative humidity, vapor pressure, and wind speed [20], the air temperature is the most widespread monitored meteorological variable among the required inputs for ET0 estimation [21].
Measuring crop evapotranspiration (ETc) and relating it to reference evapotranspiration (ET0) is the standard procedure for the determination of crop coefficient (Kc) used for skilled irrigation management [11,22]. Kc is defined as the ratio between the actual crop evapotranspiration (ETc) and reference evapotranspiration (ET0), i.e., Kc = ETc/ET0. Since biotic and abiotic stress on the crop may affect its water consumption [2,23], standard Kc needs to be determined on plants that are disease-free, well-fertilized and achieving full production and grown in a large field under optimum soil water conditions [11]. Kc varies along the growing season as a function of leaf area index (LAI) dynamics, the solar radiation intercepted by the canopy and the phenological stage of the crop [11,22]. The Kc values may vary with agricultural practices [24]. Kc has two components–Ke, soil evaporation, and Kt, plant transpiration, i.e., Kc = Ke + Kt [11].
At present, ETc in a cultivated horticultural field has been measured and estimated with different techniques, such as the microclimatological method [25], soil water balance method [26] and the method of combining measurements of soil evaporation with plant transpiration [27,28,29]. In addition, remote sensing [30,31] and lysimeter [32,33,34,35] techniques have also been applied to the determination of ETc. Lysimetry is considered the standard technique for measuring ETc. The soil water balance is a simple method for estimating ETc; however, it requires accurate estimates of its components, such as deep drainage and change in soil water content. There are difficulties in using the micrometeorological method because of canopy heterogeneity of orchards [36]. As well, remote sensing and lysimeter techniques have high costs of equipment. The combination of soil evaporation and plant transpiration measurements was proven to be more accurate than the microclimatological method and soil water balance method [37].
LAI is closely related to photosynthesis, transpiration, evapotranspiration and productivity [38]. Vegetation has a major impact on the different hydrological components of the water cycle which alters the surface energy balance. Heterogeneity in vegetation affects the physical characteristics, such as albedo, atmospheric transmissivity, root zone soil moisture, crop growth, evapotranspiration, leaf area index (LAI), stomatal conductance and surface runoff [39,40,41]. Crop coefficient Kc reflects the comprehensive influence of crop growth conditions, biological characteristics, yield level and other conditions on water demand. Crop water demand can be calculated according to reference evapotranspiration and crop coefficient Kc [42]. Crop coefficient Kc is mainly affected by plant growth and development with the change of the growing season, and is related to leaf area index (LAI) dynamics, canopy intercepted solar radiation and crop phenology [43].
Since water stress can affect crop growth and yield, it is important to estimate evapotranspiration accurately, especially in semi-arid regions [44]. Munitz, etc. [43] reported that the canopy area of wine grapevines is a reliable approach for estimating their Kc in semi-arid and arid regions with a limited water supply. Wang, etc. [45] developed appropriate region-specific crop coefficients (Kc, Ks), and both Kc and Ks to LAI is a linear regression, but to canopy conductance is an exponential one. LAI is a better indicator than canopy conductance when it is used to predict Kc and Ks of vineyards in arid Northwest China. Du, etc. [46] reported that a significant correlation was detected between LAI and transpiration at the stages of leaf expansion and fruit expanding; however, there was no correlation between LAI and transpiration at the stages of bud development and flowering, and fruit maturing. Average LAI across the whole apple growth season was significantly higher in the treatments with irrigation amounts of 500 mm than 400 mm in arid Northwest China. Gong, etc. [47] reported that the crop coefficient Kc showed a strong linear dependence on the leaf area index of apple trees. The water stress coefficient Ks was approximately 1.0 when soil moisture was above 23% and started to decrease linearly after that. The prediction of evapotranspiration in apple orchards can be made using the Food and Agriculture Organization’s crop coefficient method from commonly available meteorological data in arid Northwest China.
Crop coefficient Kc was greatly affected by climate, soil, water stress degree and other environmental conditions in different regions [48]. Significant uncertainties in estimating ET and its components have been reported in previous studies when directly using the FAO-proposed coefficient values in different regions [49]. Thus, it is necessary to identify the specific values of Kc and Ks across different agricultural ecosystems and environmental conditions, as these values provide important guidance for local irrigation practices and can be used to improve water-use efficiency [50]. The relationships of both Kc and Ks to various ecological and environmental factors (e.g., canopy conductance, LAI, and water deficit) need to be investigated in detail for the accurate prediction of evapotranspiration of apple orchards on Chinese Loess Plateau. However, as far as we know, such studies are relatively few.
In this study, the crop coefficient Kc was defined as the product of the crop coefficient KcI under no water stress and the water stress coefficient Ks under the condition of water stress, namely Kc = KcI × Ks. We selected nine apple trees 5–9-years-old in Mizhi County to measure the LAI and evapotranspiration of apple orchards in hilly apple orchards for two consecutive years. The objectives of this study were to: (1) Identify the dynamic variation characteristics of LAI in hilly apple orchards; (2) Analyze the change of ET of apple orchards and quantify the components of ET; (3) Reveal the relationship between LAI and ET; and (4) Determine the crop coefficient KcI under no water stress and the water stress coefficient Ks at different growth stages of apple trees, and establish the relationship between LAI and Kc under water deficit. It provides scientific reference for accurate prediction of evapotranspiration and guidance for water balance management of apple orchards in Chinese Loess Plateau.

2. Materials and Methods

2.1. Study Area

The field study was conducted from April to September 2019 and 2020 at the Mizhi Dangta Mountain Apple Technology Demonstration Base (38°08′32″ N, 109°57′20″ E), located in the Yulin Mountain Apple Test Station, Shanxi Province of northwest China. Mizhi County is located in the Loess Plateau with an altitude of 847.2–1255.2 m. The site is a semi-arid climate in the temperate zone, where the annual average temperature and total radiation are 8.9 °C and 587.1 kJ cm−2, respectively, and the annual average rainfall is 40–500 mm, where most of the rainfall from April to June was an ineffective rainfall of less than 5 mm, and more than 60% of the rainfall occurred from July to September. Moreover, the rainfall is mostly in the form of heavy rain with high intensity. The total frost-free period was 162 days. The soil type of the experimental site is yellow cotton with an average bulk density of 1.20–1.35 g cm−3. The field water capacity and soil organic matter content were 22.4% and 4.2 g∙kg−1, respectively.
The area of the experimental apple orchard is 13.5 hm2; the main varieties are Red Fuji, pollinating varieties are Gala and Qinguan and so on. The plant row spacing was 4 m × 5 m. Nine Red Fuji apple trees of different sizes were selected as sample plants; their canopy structure is shown in Table 1, of which 3 apple trees with the same age and growth status were irrigated (the number is 7, 8 and 9) in order to determine the crop coefficient KcI under no water stress. In order to avoid water stress in the apple tree growing season, timely irrigation should be carried out when the soil water content is less than or close to 60% field capacity. Considering that the uneven leaf coverage of apple trees leads to the heterogeneity of solar radiation on the soil surface, 5 micro-lysimeters were installed around each fruit tree, and the location is shown in Figure 1. A field weather station was installed in the center of the orchard to monitor meteorological data. Management measures, such as fertilization, weed and insect control were the same for all sample plants.

2.2. Leaf Area Index Measurement

The leaf area index (LAI) was measured using a standard digital SLR camera equipped with an optical fish-eye lens and was used to take 8 pictures of each apple tree [51]. LAI was measured every 15 days during the apple growth period from April to September in 2019 and 2020. In order to make sure the camera was level and shoot up or down depending on the height of the plant, and to ensure the accuracy of the measurements, canopy photographs were taken before sunrise or after sunset on cloudy and sunny days. The captured images were processed using CAN_EYE Software V 6.47 and the average LAI was obtained using CAN_EYE V 5.1.

2.3. Monitoring Soil Water Contents

Soil samples were collected at intervals of 20 cm in a 0–200 cm soil profile using a soil auger. Soil moisture was measured at 5–7 day intervals by oven at 105 °C drying method during the experiment, increased the number of times before and after irrigation and delayed by a day or two when rainfall events occurred. Soil water storage was obtained by integrating soil water content, soil bulk density and respect to depth.

2.4. Measurement of Evapotranspiration

The total evapotranspiration of apple trees was calculated by adding the transpiration (T), soil evaporation (E) and canopy interception (Ic). Evapotranspiration can be expressed as:
E T = E + T + I c

2.4.1. Soil Evaporation

Soil evaporation of 0–20 cm was measured with a micro-lysimeter [52] made from PVC tubes, with 11 cm in diameter and 20 cm in height during non-rainfall periods. Five micro-lysimeters were installed and numbered around each apple tree, respectively (Figure 1), and there were 45 micro-lysimeters in total. The micro-lysimeters were weighed every day at 08:00 using an electric balance with a precision of 0.1 g. The soil in the micro-lysimeters was replaced every 3–5 days. Soil evaporation (E, mm d−1) was calculated as follows:
E = 10 × Δ W / ρ π ( D / 2 ) 2
where the value of 10 is the conversion factor; ΔW(g) is the weight change of the micro-lysimeter in 24 h; ρ is the density of water (g cm−3); D is the internal diameter of the micro-lysimeters (cm).

2.4.2. Transpiration

The SF velocities (L h−1) of selected experimental apple trees were measured routinely by using thermal-dissipation probes [53,54]. A pair of probes (model: TDP-10) with 10 mm in length and approximately 1.2 mm in diameter were inserted into the trunk sapwood of each tree (one heated and one unheated). The monitoring data were recorded every half an hour. The sap-flux velocity of individual trees was calculated as [55,56]:
SF = A s × 0.0119 ( Δ T max Δ T Δ T ) 1.231 × 3600
where As is the sapwood area (cm2); ∆Tmax is the maximum temperature difference (°C) of the two probes (the heated and unheated) when there is no liquid flow before dawn; ∆T is the difference in instantaneous temperature difference (°C) between the heated and unheated probes. The results of ∆Tmax and ∆T are obtained automatically by the instrument.
In order to avoid damaging the sampled trees and ensure the continuity of the experiment, the sapwood area was measured outside the experimental plot and the relationship was determined between the sapwood area and diameter at breast height (DBH) for calculating daily cumulative SF [57] (Figure 2).
The daily transpiration of apple (Tday, mm) can be calculated as:
T day = n × SF day S = n × 24 × SF ¯ S
where n is the number of selected apple trees; SFday is the daily liquid flux of a single tree (L day−1); S is the experiment plot area (m2); S F ¯ is the the average sapwood liquid flux for 24 h (L h−1).

2.4.3. Canopy Interception

The canopy interception was estimated by subtracting stemflow and throughfall from precipitation; the calculation formula is as follows:
I c = P S t F T s
where P is gross precipitation (mm); Ts is throughfall (mm); StF is stemflow (mm).
Throughfall (Ts): During the study period, a homemade rain gauge was used to measure penetrating rainfall under the canopy of selected trees. Ts was recorded within 30 min of each gross precipitation event. The recorded values were converted to water depths (mm) per unit area [58].
Stemflow (StF): Corrugated plastic tubes (3.0 cm diameter) were longitudinally halved, the contact between the lower edge of the plastic tubes and the trunk of the tree was sealed with silicone, and the end of the plastic tubes was connected to a 20 L plastic bucket. The data were recorded 30 min after the end of the rainfall event. The measured water quantity divided by the average canopy projected area of a single apple tree was converted to stemflow depth (mm) [59].

2.5. Meteorological

Precipitation, net radiation, air temperature, relative humidity, wind speed and soil heat flux were continuously monitored by a small farmland weather station produced by Campbell Scientific, located in the center of the experimental orchard, and the data were recorded every 30 min. In order to meet the requirements of ET0 estimated by the PM equation, irrigation was carried out irregularly around the weather station. Since late April, grass with a height of more than 10 cm was always fully covered on the surface.

2.6. Reference Evapotranspiration

Daily reference crop evapotranspiration (ET0, mm day−1) was calculated using the widely used FAO56-Penman-Monteith equation [2]. The calculation is as follows:
ET 0 = 0.408 Δ ( R n G ) + γ 900 T + 273 U 2 ( e s e a ) Δ + γ ( 1 + 0.34 U 2 )
where Rn is the net radiation (mJ m−2 d−1); G is soil heat flux density (MJ m−2 d−1); γ is the hygrometer constant (kPa °C−1); T is the mean air temperature at 2 m height (°C); U2 is the wind speed at 2 m height (m s−1); es − ea is the difference of saturated vapor pressure (kPa); Δ is the slope of saturated vapor pressure with temperature (kPa °C−1).

2.7. Crop Coefficient Calculations

Crop coefficient KcI under no water stress can be determined by the ratio of evapotranspiration to reference evapotranspiration under irrigation. The calculation is as follows:
Kc 1 = ET I ET 0
where ETI is the evapotranspiration under irrigation conditions (mm d−1); ET0 is the reference evapotranspiration calculated by the Penman-Monteith equation (mm d−1).
The water stress coefficient Ks can be calculated according to the measured evapotranspiration and reference evapotranspiration under the condition of water deficit and the crop coefficient KcI under no water stress [2]. The calculation is as follows:
Ks = ET D ET 0 × Kc I
where ETD is the evapotranspiration under water deficit (mm d−1); ET0 is the reference evapotranspiration calculated by Penman-Monteith equation (mm d−1); KcI is the crop coefficient without water stress.
Crop coefficient Kc can be calculated from the crop coefficient KcI and water stress coefficient Ks under the condition of no water stress, and the calculation is as follows:
Kc = KcI × Ks
where KcI is the crop coefficient without water stress; Ks is the water stress coefficient.

2.8. Statistical Analysis

The measured test data were statistically analyzed by Excel in WPS, plotted by Excel and R language, and analyzed by SPSS19.0 for correlation and regression.

3. Results

3.1. Leaf Area Index

The change of the LAI in the apple orchard showed a rapid increase—moderate increase—declined trend from April to September in the Loess hilly region of northern Shaanxi in China (Figure 3). The change of LAI was small at the bud development period, then from mid-April to June, the LAI increased rapidly from 0.38 [m2 m−2] to 1.67 [m2 m−2] in 2019 and from 0.26 [m2 m−2] to 1.70 [m2 m−2] in 2020, due to new branch growth and leave enlargement, and at this stage apple trees need enough water to maintain normal physiological activities [24]. From July to mid-August, the LAI slowly increased from 1.67 [m2 m−2] to 1.83 [m2 m−2] in 2019 and from 1.70 [m2 m−2] to 1.86 [m2 m−2] in 2020; this is because the new shoots stopped growing and the leaf size was stable. After mid-August, the autumn treetop of apple trees grew to produce new leaves, which gradually expanded and LAI increased slightly. By mid-September, the LAI reached its maximum value of 2.06 [m2 m−2] in 2019 and 2.16 [m2 m−2] in 2020 for the entire growing season; after that, LAI slightly went down, which may be attributed to rainstorm, strong wind and other climatic conditions.

3.2. Precipitation and Reference Evapotranspiration

The rainfall distribution of the experimental orchard from April to September in 2019 and 2020 is shown in Figure 4, and 49 and 57 rainfall events were observed from April to September in 2019 and 2020, respectively, with accumulative rainfall totals of 399.5 mm and 411.9 mm, respectively. The maximum rainfall occurred on 22 July 2019 and 5 August 2020, with 53.5 mm and 95 mm, respectively. The range of monthly rainfall in 2019 and 2020 was 8.9–101.1 mm and 10.7–245.2 mm, respectively (Table 2), and the rainfall distribution was extremely uneven. The rainfall from July to September accounted for 70.8% and 90.0% of the total rainfall in the growth period (from April to September), respectively. Rainfall was the least in May and the most in July and August, with high intensity and long duration, which provided sufficient water for the soil and effectively alleviated the water stress.
Daily ET0 in 2019 and 2020 ranges from 1.98–6.88 mm d−1 and 1.68–7.76 mm d−1, respectively, with an average of 4.48 mm d−1 and 4.52 mm d−1, respectively. The maximum daily ET0 occurred on 23 August 2019 and 22 August respectively. Total ET0 in 2019 and 2020 was 829.9 mm and 842.6 mm, respectively (Table 2). The range of monthly ET0 in 2019 and 2020 was 125.5–151.9 mm and 114.2–192.7 mm, respectively (Table 2). From April to September, the monthly ET0 showed a trend of first increase and then decrease, with the maximum monthly ET0 in August and the minimum monthly ET0 in April.

3.3. Evapotranspiration

Monthly soil evaporation (E), transpiration (T), canopy interception (Ic) and evapotranspiration (ET) of the apple orchard during the study period in 2019 and 2020 are shown in Table 3. Annual ET was 424.8, 463.2 mm for 2019 and 2020, respectively. The ET in both 2019 and 2020 was greater than the respective rainfall, indicating that the orchard presents a negative water balance. Monthly ET ranged from 29.3 to 122.3 and 26.7 to 175.6 mm in 2019 and 2020, respectively. The trends of monthly ET were similar in each year. Monthly ET was highest in August for both years. The maximum components of ET in the apple orchard was soil evaporation in both 2019 and 2020, which were 203.3 mm and 227.5 mm, respectively, accounting for 47.8% and 49.1% of the total ET, respectively. The annual T were 180.6 mm and 203.2 mm for 2019 and 2020, which accounted for 42.5% and 43.9% of the annual ET, respectively. The minimum components were Ic, which were 41.0 mm and 42.3 mm for 2019 and 2020, respectively, accounting for 9.6% and 9.1% of the ET, respectively.
In April 2019, the apple trees entered the germination and flowering period, the branches and leaves were gradually developing, due to heavy rainfall, E was relatively intense and the ET was also larger, accounting for 8.4% of the total ET. The ET decreased in May, only accounting for 6.9%. Since June to August entered the rainy season, the rainfall increased and the temperature rose rapidly. The water consumption of the apple trees was intense, the ET gradually increased and reached the maximum value of 122.3 mm in August, accounting for 28.8% of the total ET. In September, the ET decreased to 89.0 mm. In the two years, the ET in April 2020 was the smallest, only 26.7 mm, accounting for 5.8% of the total annual ET. The ET increased continuously from May to July in 2020, accounting for 6.3%, 11.3% and 20.6% of the total ET, respectively. In August, the rainfall was larger and the ET increased, which accounted for 37.9% of the total ET. The ET in September decreased, and only accounted for 18.2% of the total ET.

3.4. Leaf Area Index and Evapotranspiration Relationship

The exponential relationship between LAI and evapotranspiration (ET) or transpiration (T) was strong and significant (Figure 5, p < 0.01). This indicates that LAI is an important factor affecting T and ET. In 2019, the fitting equation of LAI and T was as follows: Y = 1.0615e1.5494X (R2 = 0.8206, n = 72); the fitting equation of LAI and ET was y = 9.2836e0.7983x (R2 = 0.5334, n = 72). In 2020, the fitting equation of LAI and T was as follows: Y = 0.6911e1.6856x (R2 = 0.8011, n = 72); the fitting equation of LAI and ET was as follows: Y = 7.7238e0.9009x (R2 = 0.5417, n = 72).

3.5. Leaf Area Index and Crop Coefficient Relationship

Rainfalls, ET0, ET under different water conditions, crop coefficient KcI without water stress and water stress coefficient Ks and Kc under water deficit from April to September in 2019 and 2020 are shown in Table 4. As can be seen from Table 4, ET0 > ETI > ETD were shown in each year of the experiment, indicating that water conditions had a great influence on evapotranspiration, and irrigation increased the evapotranspiration of orchards in the study area.
The crop coefficient KcI under no water stress changed with the development of the apple tree canopy. Because of canopy development, the wet evaporation surface area increases; thus, ETI also increases under the same atmosphere demand, resulting in an increase in the crop coefficient KcI. In 2019, the KCI rose from 0.43 in April to 0.74 in July, and decreased to 0.60 in September. In 2020, the KCI rose from 0.40 in April to 0.98 in August, and decreased to 0.68 in September. The reason for the decrease of KC after August may be related to the decrease of LAI. In both 2019 and 2020, the variation of Ks was basically consistent with the variation of rainfall, and its variation range is 0.53 to 1.15. From April to June, there was less rainfall and Ks was low, and apple trees suffered serious water stress. After July, the rainfall increased sharply and Ks was approximately 1.0, indicating that the growth and development of apple trees would not be affected by water stress. The variation trend of crop coefficient Kc under water deficit was basically the same as that of KcI and Ks, which were small at the early growth stage and were 0.28 and 0.23, respectively. With the growth of apple trees, the crop coefficient Kc increased continuously, and reached the maximum value in August of 2019 and 2020, which was 0.81 and 0.91, respectively, and decreased to 0.64 and 0.67 in September.
LAI can reflect the growth status of crops, so the relationship between LAI and crop coefficient Kc can realize the dynamic simulation of evapotranspiration under the condition of water deficit. There is a significant exponential relationship between LAI and crop coefficient Kc (Figure 6), and the fitting equation was as follows: Kc = 0.1141e1.0665LAI, R2 = 0.7055, p < 0.01. These results indicated that LAI could be used to estimate the crop coefficient Kc of apple orchards under the condition of water deficit in the Loess Plateau, and the actual evapotranspiration of apple orchards in this region could be predicted.

4. Discussion

4.1. Leaf Area Index

LAI is an important parameter that characterizes the interface between a vegetation canopy and the atmosphere. It has been commonly used in plant growth, climatic and canopy reflectance models, as well as in much other research [60]. In the Loess Plateau region of China, the LAI of apple orchards is closely related to topography, soil, meteorological factors (such as rainfall and rainfall distribution, illumination, temperature, etc.) and management measures (such as density, canopy structure, soil moisture conservation measures, irrigation, fertilization and pest control). In this study, the variation trend of LAI in experiment orchards during the whole growing season was consistent with other results from previous studies [47,61,62,63]; LAI was smaller than the observation of others [47,61,62,63], which is determined by the growth of the apple tree under the local natural ecological conditions. The rainfall during the apple tree growing season was 399.5 mm and 411.9 mm in 2019 and 2020, respectively, the rainfalls were normal in experiment orchards and can represent the general rainfall conditions in the region. During the survey, apple tree pruning, fertilization, flower management, pest control, weeding and so on were normal. Therefore, we believe that the size and variation trend of the leaf area index investigated by us can represent the leaf growth status of local apple trees.

4.2. Precipitation and Reference Evapotranspiration

In the growth period of apples in 2019 and 2020, the rainfall distribution was not uniform. Most of the rainfall was concentrated in July to September, with the minimum rainfall in May and the maximum rainfall in July to August. In both 2019 and 2020, the evapotranspiration during the apple growth period was greater than the annual rainfall, which indicated that the rainfall during the growth period could not meet the demand of the water consumption of the apple orchard. The daily ET0 of the apple orchard in 2019 and 2020 ranges from 1.98–6.88 mm d−1 and 1.68–7.76 mm d−1, respectively, and the mean values were 4.48 mm d−1 and 4.52 mm d−1, respectively. As can be seen from Table 2, rainfall and ET0 are not closely related. The rainfall in April 2019 was 65.5 mm, 39.2 mm more than that in April 2020, but ET0 was only 10.8 mm more. One reason may be that the low humidity soil absorbed and stored part of the rainfall, and the other reason was basically no grass covered on the ground at this stage, which affected the accuracy of ET0 simulation by the PM equation [13], Alexandris et al. [14] pointed out that the FAO56-PM method, when applied above bare, non-irrigated soil, present overestimation. In addition, this period was the germination and leaf spreading period of apple trees, and the water consumption of transpiration was small, indicating that the large ET0 was mainly caused by soil evaporation. Rainfall is low in May and June, and the low soil moisture (or completely dry soil surface) and the persistence of dry conditions in the overlying to the soil surface air layer, results in energy allocation from the ground heat flux, attributing to the overlying to the soil’s surface air layer. But grass had already covered the ground, so the estimated values of ET0 in the experimental orchard can be considered as very close to ET0 (grass reference height 8–12 cm). Strong winds are frequent in the experimental orchard, so the possible reason for the high ET0 is wind speed. There was more rainfall and higher temperatures in July and August, the ground was full of weeds in the experimental orchard, ET0 was high due to high transpiration consumption and high soil surface evaporation. After September, the temperature drops, the sunshine hours are shortened, and ET0 becomes smaller.
The variation of daily and monthly ET0 for the two years had no obvious regularity, partly because ET0 is the result of a combination of meteorological factors, and the sensitivity of ET0 to climatic factors varies on daily and seasonal scales [57]. Annual ET0 was higher for 2019–2020 than that reported by Wang et al. [57] for the Loess Plateau. We speculate that the higher ET0 was mostly caused by higher wind speed, low-relative humidity and more hours of sunshine. This was not in agreement with the views of Li et al. [64] and Wang et al. [56]. It should be noted here that the main influencing factors of ET0 on the Loess Plateau vary with time and the regional climate environment. In addition, the PM equation was sensitive to errors induced by non-representative weather data; experimental weather stations situated above a dry soil will measure higher temperature and correspondingly, lower relative humidity values, compared to a station on well-irrigated soil [65,66]. This may result in an overestimation of ET0.

4.3. Evapotranspiration

Evapotranspiration (ET) consists of soil evaporation (E), transpiration (T) and canopy interception (Ic). Soil evaporation could not be precisely measured on rainy days, so we present the results only for monthly ET. Quantification of the components of evapotranspiration is the key to designing management strategies for improving water use efficiency [47]. The factors that influence E, T and Ic would thus influence ET. In this study, E and T were the largest proportions of ET. E/ET was greater than T/ET, which is consistent with the research results of Li et al. [67] and Wang et al. [57]. The main reason is that the LAI was low in the orchard of our study, and the surface soil water after rainfall was the main source of ET loss [68]. However, both Wang [69] and Yu [70] believed that the largest component of ET was T, which may be related to tree age, density, rainfall and other natural ecological conditions. Ic was always small in the whole growth period, and Ic only accounted for 1.7–13.6% of ET. Ic showed different characteristics in different months, which was consistent with the research results of Li et al. [71]; however, Ic should not be ignored [61], because the physiological and ecological characteristics of plants (leaf area index, plant height) and various meteorological factors (rainfall, rainfall intensity and wind speed) will affect the plant canopy interception. If an evapotranspiration assessment does not consider canopy interception volume in apple orchards, it may be misleading regarding the relationship between water supply and demand. During the two years of the experiment, the maximum monthly ET in the apple orchard occurred in August. One reason is that this stage is during the rapid expansion of fruit and fruit development consumed a lot of soil water; another reason is more evaporation due to high temperatures and more rainfall. Because the apple trees in August were in the period of fruit expansion, during the fast growth rate, the fruit consumed a lot of water [46]. The minimum monthly ET appeared in May and April, This was because transpiration consumption was lower due to smaller leaf area and LAI, although less rainfall and strong winds caused more evaporation.
In this study, T and E were measured separately by thermal-dissipation probes and micro-lysimeters at localized scale. Each method for estimating T or E of whole apple orchards faced its own problems due to sampling constraints [57,72]. The limitation of estimation of actual evapotranspiration in apple orchards was mainly related to the monitoring of soil evaporation by microlysimeters. Firstly, replacing the soil in the micro-lysimeters would be essential every 3–5 days and after a heavy rain (above 5 mm d−1), which can become time consuming. Secondly, the time-resolution of the method was constrained by the weighing frequency of the micro-lysimeters [57].

4.4. Leaf Area Index and Evapotranspiration Relationship

Evapotranspiration intensity is jointly affected by soil moisture, solar radiation and crop LAI [73]. Our results showed that there was a positive correlation between LAI and evapotranspiration in apple orchards in the Loess hilly region, compared with the relationship between LAI and evapotranspiration, LAI is more closely related to transpiration (Figure 5), however, there was no significant correlation between LAI and soil evaporation. Our results were consistent with most previous research results. For example, Du et al. [46] reported that there was a positive linear relationship between apple transpiration and LAI at the stages of leaf expansion and fruit expanding in apple orchards in Northwest China under irrigation conditions, however, there was no correlation between LAI and transpiration at the stages of bud development and flowering, and fruit maturing because the leaves were not developed during the stage of bud development and flowering, and the leaves withered during the stage of fruit maturing. Liu et al. [63] reported that there was a positive linear correlation between LAI and sapflow in apple orchards in arid areas of Northwest China. Testi et al. [29] reported that there was a significant linear growth relationship between evapotranspiration and LAI of olive trees in southern Spain when there was no moisture restriction in summer. Wang et al. [74] reported evapotranspiration, canopy interception, and transpiration were all significantly positively correlated with the canopy LAI of the plantation, while soil evaporation was significantly negatively correlated with canopy LAI. Juhász and Hrotkó [75] reported that canopy transpiration was positively correlated with LAI of sweet cherry in Hungarian. Benyon et al. [76] found that there was a positive correlation between annual transpiration and LAI of Radiata Pine and Eucalyptus globulus in south-eastern Australia. Almeida et al. [77] reported that forest evapotranspiration was positively correlated with LAI of Eucalyptus grandis on the Atlantic coast of Brazil. Tian et al. [78] studied the response of water yield to key stand structure and site factors in the Liupanshan larch plantation, established and fitted canopy LAI and stand evapotranspiration models in the growing season, and calculated water yield in the growing season using water budget method. Zhang et al. [79] reported that LAI of winter wheat was positively correlated with evapotranspiration in Changwu tableland, Shaanxi.
However, there are a few studies that differ from our results. Juul et al. [80] found that there was a significant negative correlation between evapotranspiration and LAI of Betula latticola in the Netherlands. Li et al. [81] reported evapotranspiration and transpiration showed a curve growth with LAI of apple trees under sufficient water supply.
In this study, the rainfall during the apple tree growing season was 399.5 mm and 411.9 mm in 2019 and 2020, respectively and the rainfalls were the general rainfall conditions in the region. Rainfall varies widely from year to year, usually by approximately 20% from year to year. This difference in rainfall will inevitably lead to differences in the growth of apple branches and leaves. Therefore, the relationship between LAI and ET in rain-deficient years and rain-rich years needs to be further studied.

4.5. Leaf Area Index and Crop Coefficient Relationship

Crop coefficient Kc reflects the comprehensive influence of crop growth, biological characteristics, yield level and other conditions on water demand. The crop coefficient approach is based on the assumption that the reference evapotranspiration ET0 accounts for the atmospheric effect on the evapotranspiration process, whereas the crop coefficient Kc describes the integrated effects of aerodynamic and canopy resistances to water vapor flux, and soil evaporation [2]. Crop water demand can be calculated according to ET0 and crop coefficient Kc [42]. In this study, the crop coefficient Kc was defined as the product of the crop coefficient KcI under no water stress and the water stress coefficient Ks under water stress, namely Kc = KcI × Ks.
The crop coefficient KcI under no water stress changes with the development of the apple tree canopy. Because of canopy development, the wet evaporation surface area increases; thus, ET also increases under the same atmosphere demand, resulting in an increase in the crop coefficient KcI. The variation trend of KcI during the whole growing season was consistent with LAI, and it agreed with other results from previous studies [42,47]. There was an excellent linear relationship between crop coefficient KcI and LAI, indicating that crop coefficient KcI can be estimated using LAI. In this study, the crop coefficient KcI under no water stress was small, as compared to the results of Gong et al. [47]. LAI is also smaller than the results of Gong et al. [47], indicating the smaller KcI may be related to smaller LAI.
If there is water stress under natural conditions, it is necessary to consider the water stress coefficient Ks to estimate the actual crop evapotranspiration. In this study, the variation of Ks is basically consistent with the variation of rainfall. From April to June, there was less rainfall and Ks was low, and apple trees suffered serious water stress. After July, the rainfall increased sharply and Ks was approximately 1.0, indicating that the growth and development of apple trees would not be affected by water stress. The reason of Ks from April to June is that a low soil water availability increases the difficulty of root uptake water, which stimulates the production of abscisic acid (ABA) in roots. An increase in root-source ABA content in leaves should decrease the opening of leaf stomata, which limits the transpiration of plants [47]. Our results showed that there was a significant exponential relationship between LAI with crop coefficient Kc (p < 0.01) under water deficit. It indicated that the index model established by LAI can accurately estimate the crop coefficient Kc, which further provides a theoretical basis for the prediction of evapotranspiration under the condition of water deficit in Chinese Loess Plateau.
It should be noted that this study was carried out under specific site conditions. However, in the Loess Plateau region of China, the site conditions of apple orchards, such as elevation, topography, slope and slope position vary greatly, which inevitably leads to great differences in meteorological data. Therefore, the relationship between LAI and KC under different site conditions needs to be further studied.

5. Conclusions

The change of LAI in apple orchards showed a rapid increase—moderate increase—declined trend from April to September in the Loess hilly region of northern Shaanxi in China. The variation range of LAI was 0.26–2.16 [m2 m−2].
The ET of the orchard was greater than the rainfall; the maximum ET was in July or August. The maximum components of ET in apple orchards was E, and E accounted for 47.8–49.1% of ET; T accounted for 42.5–43.9% of ET and Ic accounted for only 9.1–9.6% of ET. There was a significant exponential relationship between LAI and T or ET. The crop coefficient KcI under no water stress changed with the development of the apple tree canopy. The variation of water stress Ks was basically consistent with the variation of rainfall. There is a significant exponential relationship between LAI and crop coefficient Kc under water deficit (Kc = 0.1141e1.0665LAI, R2 = 0.7055, p < 0.01). LAI could be used to estimate the crop coefficient Kc of apple orchards under water deficit in the Loess Plateau, and the actual evapotranspiration of apple orchards in this region could be predicted. This study was carried out under normal rainfall conditions and specific site conditions. Therefore, the relationship between LAI and KC should be further studied under the conditions of low rainfall years and rich rainfall years, as well as different site conditions.

Author Contributions

Conceptualization, Q.J. methodology, Q.J. and Y.-P.W.; formal analysis, Q.J.; data curation, Q.J.; writing-original draft preparation, Q.J.; writing-review and editing, Q.J. and Y.-P.W.; supervision, Y.-P.W.; project administration, Y.-P.W.; funding acquisition, Y.-P.W. Both authors have read and agreed to the published version of the manuscript.

Funding

Contract/grant sponsor: National Natural Science Foundation of China; contract/grant number: 41571218, 41401613. Contract/grant sponsor: Science-Technology Innovation Project of Shaanxi; contract/grant number: 2021NY-204.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

We are grateful for the support from the staff of the Mizhi Experimental Station of Northwest A&F University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Raz-Yaseef, N.; Yakir, D.; Schiller, G.; Cohen, S. Dynamics of evapotranspiration partitioning in a semi-arid forest as affected by temporal rainfall patterns. Agric. For. Meteorol. 2012, 157, 77–85. [Google Scholar] [CrossRef]
  2. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements. FAO Irrig. Drain. Pap. 1998, 56, 300. [Google Scholar]
  3. Srivastava, A.; Sahoo, B.; Raghuwanshi, N.S.; Singh, R. Evaluation of variable-infiltration capacity model and MODIS-terra satellite-derived grid-scale evapotranspiration estimates in a River Basin with Tropical Monsoon-Type climatology. J. Irrig. Drain. Eng. 2017, 143, 04017028. [Google Scholar] [CrossRef] [Green Version]
  4. Sheffield, J.; Wood, E.F. Projected changes in drought occurrence under future global warming from multi-model, multiscenario, IPCC AR4 simulations. Clim. Dyn. 2008, 31, 79–105. [Google Scholar] [CrossRef]
  5. Temesgen, B.; Eching, S.; Davidoff, B.; Frame, K. Comparison of Some Reference Evapotranspiration Equations for California. J. Irrig. Drain Eng. 2005, 131, 73–84. [Google Scholar] [CrossRef]
  6. Walter, I.A.; Allen, R.G.; Elliott, R.; Jensen, M.E.; Itenfisu, D.; Mecham, B.; Howell, T.A.; Snyder, R.; Brown, P.; Echings, S.; et al. ASCE’s Standardized Reference Evapotranspiration Equation. In Proceedings of the Watershed Management and Operations Management 2000, Fort Collins, CO, USA, 20–24 June 2000; pp. 1–11. [Google Scholar]
  7. Howes, D.J.; Fox, P.; Hutton, P. Evapotranspiration from Natural Vegetation in the Central Valley of California: Monthly Grass Reference-Based Vegetation Coefficients and the Dual Crop Coefficient Approach. J. Hydrol. Eng. 2015, 20, 04015004. [Google Scholar] [CrossRef]
  8. Gurski, B.C.; Jerszurki, D.; De Souza, J.L.M. Alternative reference evapotranspiration methods for the main climate types of the state of Paraná, Brazil. Pesqui. Agropecu. Bras. 2018, 53, 1003–1010. [Google Scholar] [CrossRef] [Green Version]
  9. Akumaga, U.; Tarhule, A.; Yusuf, A.A. Validation and testing of the FAO AquaCrop model under different levels of nitrogen fertilizer on rainfed maize in Nigeria, West Africa. Agric. For. Meteorol. 2017, 232, 225–234. [Google Scholar] [CrossRef]
  10. Gundalia, M.; Dholakia, M. Modelling daily reference evapotranspiration in Middle South Saurashtra Region of India for monsoon season using dominant meteorological variables and the FAO-56 penman-monteith method. Int. J. Sustain. Water Environ. Syst. 2016, 8, 101–108. [Google Scholar]
  11. Allen, R.G.; Pereira, L.; Raes, D. Crop Evapotranspiration. Guidlines for Computing Crop Water Requirments; FAO: Rome, Italy, 2006. [Google Scholar]
  12. Alexandris, S.; Proutsos, N. How signifificant is the effect of the surface characteristics on the Reference Evapotranspiration estimates? Agric. Water Manag. 2020, 237, 106181. [Google Scholar] [CrossRef]
  13. Jabloun, M.D.; Sahli, A. Evaluation of FAO-56 methodology for estimating reference evapotranspiration using limited climatic data: Application to Tunisia. Agric. Water Manag. 2008, 95, 707–715. [Google Scholar] [CrossRef]
  14. Almorox, J.; Quej, V.H.; Martí, P. Global performance ranking of temperature-based approaches for evapotranspiration estimation considering Köppen climate classes. J. Hydrol. 2015, 528, 514–522. [Google Scholar] [CrossRef]
  15. Awal, R.; Habibi, H.; Fares, A.; Deb, S. Estimating reference crop evapotranspiration under limited climate data in West Texas. J. Hydrol. Reg. Stud. 2020, 28, 100677. [Google Scholar] [CrossRef]
  16. Kwon, H.; Choi, M. Error assessment of climate variables for FAO-56 reference evapotranspiration. Theor. Appl. Clim. 2011, 112, 81–90. [Google Scholar] [CrossRef]
  17. Samani, Z. Estimating Solar Radiation and Evapotranspiration Using Minimum Climatological Data. J. Irrig. Drain. Eng. 2000, 126, 265–267. [Google Scholar] [CrossRef]
  18. Jensen, M.E.; Allen, R.G. Evaporation, Evapotranspiration, and Irrigation Water Requirements; American Society of Civil Engineers: Reston, VA, USA, 2016. [Google Scholar]
  19. Di Stefano, C.; Ferro, V. Estimation of Evapotranspiration by Hargreaves Formula and Remotely Sensed Data in Semi-arid Mediterranean Areas. J. Agric. Eng. Res. 1997, 68, 189–199. [Google Scholar] [CrossRef]
  20. Vanderlinden, K.; Giraldez, J.; Van Meirvenne, M. Assessing Reference Evapotranspiration by the Hargreaves Method in Southern Spain. J. Irrig. Drain. Eng. 2004, 130, 184–191. [Google Scholar] [CrossRef]
  21. Mendicino, G.; Senatore, A. Regionalization of the Hargreaves coefficient for the assessment of distributed reference evapotranspiration in Southern Italy. J. Irrig. Drain. Eng. 2013, 139, 349–362. [Google Scholar] [CrossRef]
  22. 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]
  23. Netzer, Y.; Shenker, M.; Schwartz, A. Effects of irrigation using treated wastewater on table grape vineyards: Dynamics of sodium accumulation in soil and plant. Irrig. Sci. 2014, 32, 283–294. [Google Scholar] [CrossRef]
  24. Williams, L.E.; Fidelibus, M.W. Measured and estimated water use and crop coefficients of grapevines trained to overhead trellis systems in California’s San Joaquin Valley. Irrig. Sci. 2016, 34, 431–441. [Google Scholar] [CrossRef]
  25. Carrasco-Benavides, M.; Ortega-Farías, S.; Lagos, L.O.; Kleissl, J.; Morales, L.; PobleteEcheverría, C.; Allen, R.G. Crop coeffiffifficients and actual evapotranspiration of a drip-irrigated Merlot vineyard using multispectral satellite images. Irrig. Sci. 2012, 30, 485–497. [Google Scholar] [CrossRef]
  26. Palomo, M.J.; Díaz-Espejo, A.; Férnandez, J.E.; Moreno, F.; Girón, I.F. Water balance in an olive orchard. Acta Hortic. 2000, 537, 573–580. [Google Scholar] [CrossRef] [Green Version]
  27. Dragoni, D.; Lakso, A.N.; Piccioni, R.M.; Tarara, J.M. Transpiration of grapevines in the humid northeastern United States. Am. J. Enol. Vitic. 2006, 57, 460–467. [Google Scholar]
  28. Intrigliolo, D.S.; Lakso, A.N.; Piccioni, R.M. Grapevine cv. ‘riesling’ water use in the northeastern United States. Irrig. Sci. 2009, 27, 253–262. [Google Scholar] [CrossRef]
  29. Testi, L.; Villalobos, F.J.; Orgaz, F. Evapotranspiration of a young irrigated olive orchard in southern Spain. Agric. For. Meteorol. 2004, 121, 1–18. [Google Scholar] [CrossRef]
  30. Rozenstein, O.; Haymann, N.; Kaplan, G.; Tanny, J. Estimating cotton water consumption using a time series of Sentinel-2 imagery. Agric. Water Manag. 2018, 207, 44–52. [Google Scholar] [CrossRef]
  31. Vanino, S.; Pulighe, G.; Nino, P.; De Michele, C.; Bolognesi, S.; D’Urso, G. Estimation of evapotranspiration and crop coeffiffifficients of tendone vineyards using multi-sensor remote sensing data in a mediterranean environment. Remote Sens. 2015, 7, 14708–14730. [Google Scholar] [CrossRef] [Green Version]
  32. Girona, J.; del Campo, J.; Mata, M.; Lopez, G.; Marsal, J. A comparative study of apple and pear tree water consumption measured with two weighing lysimeters. Irrig. Sci. 2010, 29, 55–63. [Google Scholar] [CrossRef]
  33. Lorite, I.; Santos, C.; Testi, L.; Fereres, E. Design and construction of a large weighing lysimeter in an almond orchard. Span. J. Agric. Res. 2012, 10, 238. [Google Scholar] [CrossRef] [Green Version]
  34. García-Tejero, I.F.; Hernandez, A.; Rodriguez, V.; Ponce, J.; Ramos, V.; Muriel, J.; DuranZuazo, V.H. Estimating almond crop coeffiffifficients and physiological response to water stress in semiarid environments (SW Spain). J. Agric. Sci. Technol. 2015, 17, 1255–1266. [Google Scholar]
  35. Ben-Gal, A.; Kool, D.; Agam, N.; van Halsema, G.E.; Yermiyahu, U.; Yafe, A.; Presnov, E.; Erel, R.; Majdop, A.; Zipori, I.; et al. Whole-tree water balance and indicators for short-term drought stress in non-bearing ‘Barnea’ olives. Agric. Water Manag. 2010, 98, 124–133. [Google Scholar] [CrossRef]
  36. 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]
  37. Trambouze, W.; Bertuzzi, P.; Voltz, M. Comparison of methods for estimating actual evapotranspiration in a row-cropped vineyard. Agric. For. Meteorol. 1998, 91, 193–208. [Google Scholar] [CrossRef]
  38. Sonnentag, O.; Talbot, J.; Chen, J.M.; Roulet, N.T. Using direct or indirect measurements of leaf area index to chaeacterize the shrub canopy in an ombrotrophic peat land. Agric. For. Meteorol. 2007, 144, 200–212. [Google Scholar] [CrossRef]
  39. Srivastava, A.; Kumari, N.; Maza, M. Hydrological Response to Agricultural Land Use Heterogeneity Using Variable Infiltration Capacity Model. Water Resour. Manag. 2020, 34, 3779–3794. [Google Scholar] [CrossRef]
  40. Aghsaei, H.; Dinan, N.M.; Moridi, A.; Asadolahi, Z.; Delavar, M.; Fohrer, N.; Wagner, P.D. Effects of dynamic land use/land cover change on water resources and sediment yield in the Anzali wetland catchment, Gilan, Iran. Sci. Total Environ. 2020, 712, 136449. [Google Scholar] [CrossRef] [PubMed]
  41. Liu, Y.B.; Hu, Z.H.; LI, Q.; Zhang, X.S.; Zhang, Q. Effect of leaf area index change on evapotranspiration and water yield in northern China. Chin. J. Eco-Agric. 2017, 25, 1206–1215. [Google Scholar]
  42. Meng, Q.Q. Soil Moisture Consumption Pattern and Growth Response of Apple Orchard in the Loess Plateau. Ph.D. Thesis, Northwest A&F University, Yangling, China, 2011. [Google Scholar]
  43. Munitz, S.; Schwartz, A.; Netzer, Y. Water consumption, crop coeffiffifficient and leaf area relations of a Vitis vinifera cv. ‘Cabernet Sauvignon’ vineyard. Agric. Water Manag. 2019, 219, 86–94. [Google Scholar] [CrossRef]
  44. Mansour, T.; Daniel, M.; Prasanna, G.; Jean, L.; Steiner, H.K.; Tadesse, A.M.; Nelson, P.S. Simultaneous calibration of evapotranspiration and crop yield in agronomic system modeling using the APEX model. Agric. Water Manag. 2018, 208, 299–306. [Google Scholar]
  45. Wang, S.; Xia, D.; Ma, J.; Han, T.; Zhu, G. The characteristics of evapotranspiration and crop coefficients of an irrigated vineyard in arid Northwest China. Agric. Water Manag. 2019, 212, 388–398. [Google Scholar] [CrossRef]
  46. Du, S.Q.; Kang, S.Z.; Li, F.S.; Du, T.S. Water use efficiency is improved by alternate partial root-zone irrigation of apple in arid northwest China. Agric. Water Manag. 2017, 179, 184–192. [Google Scholar] [CrossRef]
  47. Gong, D.Z.; Kang, S.Z.; Yao, L.M.; Zhang, L. Estimation of evapotranspiration and its components from an apple orchard in northwest China using sap flow and water balance methods. Hydrol. Process. 2006, 21, 931–938. [Google Scholar] [CrossRef]
  48. Wang, Y.N.; Hang, X.M.; Lu, L.; Gu, N.; Wang, Z.L.; Liu, M.; Wang, G.Q. Estimation of crop coefficient and evapotranspiration of summer maize by path analysis combined with BP neural network. Trans. Chin. Soc. Agric. Eng. 2020, 36, 109–116. [Google Scholar]
  49. Sánchez, J.M.; López-Urrea, R.; Rubio, E.; González-Piqueras, J.; Caselles, V. Assessing crop coefficients of sunflflower and canola using two-source energy balance and thermal radiometry. Agric. Water Manag. 2014, 137, 23–29. [Google Scholar] [CrossRef]
  50. Yang, P.; Hu, H.; Tian, F.; Zhang, Z.; Dai, C. Crop coeffiffifficient for cotton under plastic mulch and drip irrigation based on eddy covariance observation in an arid area of northwestern China. Agric. Water Manag. 2016, 171, 21–30. [Google Scholar] [CrossRef] [Green Version]
  51. Zhang, R.; Wang, D.; Yang, Z.; Seki, K.; Singh, M.; Wang, L. Changes in rainfall partitioning and its effect on soil water replenishment after the conversion of croplands into apple orchards on the Loess Plateau. Agric. Ecosyst. Environ. 2021, 312, 107342. [Google Scholar] [CrossRef]
  52. Ritchie, J.T. Model for predicting evaporation from a row crop with incomplete cover. Water Resour. Res. 1972, 8, 1204–1213. [Google Scholar] [CrossRef] [Green Version]
  53. Granier, A. A new method of sap flow measurement in tree stems. Ann. Sci. For. 1985, 42, 193–200. [Google Scholar] [CrossRef]
  54. Granier, A. Evaluation of transpiration in a Douglas-fir stand by means of sapflow measurements. Tree Physiol. 1987, 3, 309–320. [Google Scholar] [CrossRef]
  55. Oren, R.; Pataki, D.E. Transpiration in response to variation in microclimate and soil moisture in southeastern deciduous forests. Oecologia 2001, 127, 549–559. [Google Scholar] [CrossRef]
  56. Santiago, L.S.; Goldstein, G.; Meinzer, F.C.; Fownes, J.H.; Mueller-Dombois, D. Transpiration and forest structure in relation to soil waterlogging in a Hawaiian montane cloud forest. Tree Physiol. 2000, 20, 673–681. [Google Scholar] [CrossRef] [Green Version]
  57. Wang, D.; Wang, L. Dynamics of evapotranspiration partitioning for apple trees of different ages in a semiarid region of northwest China. Agric. Water Manag. 2017, 191, 1–15. [Google Scholar] [CrossRef]
  58. Wei, S.P.; Wang, L.; Wu, F.Q. Hydrological propeties of canopy of acacia in loess hilly and gully region. J. Nanjing For. Univ. 2008, 32, 43–48. [Google Scholar]
  59. Lu, J.F.; Ma, Q.Y.; Liu, S.H.; Tian, P.; Chen, Y. Rainfall interception capacity of a Pinus tabularformis plantation in Miyun, Beijing. J. Beijing For. Univ. 2005, 27, 129–132. [Google Scholar]
  60. Liu, C.; Sun, G.; McNulty, S.G.; Kang, S. An Improved Evapotranspiration Model for an Apple Orchard in Northwestern China. Trans. ASABE 2015, 58, 1253–1264. [Google Scholar]
  61. Guo, X.-M.; Yang, X.-T.; Chen, M.-X.; Li, M.; Wang, Y.-A. A model with leaf area index and apple size parameters for 2.4 GHz radio propagation in apple orchards. Precis. Agric. 2014, 16, 180–200. [Google Scholar] [CrossRef]
  62. Wang, D.; Wang, L. Canopy interception of apple orchards should not be ignored when assessing evapotranspiration partitioning on the Loess Plateau in China. Hydrol. Process. 2018, 33, 372–382. [Google Scholar] [CrossRef]
  63. Liu, C.; Du, T.; Li, F.; Kang, S.; Li, S.; Tong, L. Trunk sap flow characteristics during two growth stages of apple tree and its relationships with affecting factors in an arid region of northwest China. Agric. Water Manag. 2012, 104, 193–202. [Google Scholar] [CrossRef]
  64. Li, Z.; Zheng, F.-L.; Liu, W. Spatiotemporal characteristics of reference evapotranspiration during 1961–2009 and its projected changes during 2011–2099 on the Loess Plateau of China. Agric. For. Meteorol. 2012, 154–155, 147–155. [Google Scholar] [CrossRef]
  65. Yoder, R.E.; Ley, T.W.; Elliott, R.L. Measurement and reporting practices for automatic agricultural weather stations. In Proceedings of the 4th National Irrigation Symposium, ASAE, St. Joseph, MI, USA, 14–16 November 2000. [Google Scholar]
  66. Jensen, D.T.; Hargreaves, G.H.; Temesgen, B.; Allen, R.G. Computation of ET0 under nonideal conditions. J. Irrig. Drain. Eng. 1997, 123, 394–400. [Google Scholar] [CrossRef]
  67. Li, J.Y. Soil Water Dynamics and Evapotranspiration of Hilly Apple Orchard in the Loess Plateau. Ph.D. Thesis, Northwest A&F University, Yangling, China, 2017. [Google Scholar]
  68. Cavanaugh, M.L.; Kurc, S.A.; Scott, R.L. Evapotranspiration partitioning in semiarid shrubland ecosystems: A two-site evaluation of soil moisture control on transpiration. Ecohydrology 2011, 4, 671–681. [Google Scholar] [CrossRef]
  69. Wang, Y.P. Water Balance and Evapotranspiration of Apple Orchards on the Gully Region of Loess Plateau; Northwest A&F University: Yangling, China, 2014. [Google Scholar]
  70. Yu, J.F. Study on the Characteristics of Apple Tree Transpiration and Orchard Evapotranspiration on the Loess Plateau. Ph.D. Thesis, Northwest A&F University, Yangling, China, 2010. [Google Scholar]
  71. Li, M.X.; Du, S.N.; Bai, G.S.; Geng, G.J. Effects of renewal pruning on soil moisture and growth of apple tree. J. Zhejiang Univ. 2012, 38, 467–476. [Google Scholar]
  72. Zhang, Y.Q.; Kang, S.Z.; Ward, E.J.; Ding, R.S.; Zhang, X.; Zheng, R. Evapotranspiration Compon. Determ. by sapflow and microlysimetry techniques of a vineyard in northwest China: Dynamics and inflfluential factors. Agric. Water Manag. 2011, 98, 1207–1214. [Google Scholar] [CrossRef]
  73. Yang, X.G.; Liu, H.L.; Wang, Y.L.; Yu, H.N. Research on evapotranspiration of field ecological system of summer maize in Northern Plain. Chin. J. Eco-Agric. 2003, 11, 71–73. [Google Scholar]
  74. Wang, Y.R.; Wang, Y.H.; Yu, P.T.; Xiong, W.; Du, A.P.; Li, Z.H.; Liu, Z.B.; Ren, L.; Xu, L.H.; Zuo, H.J. Simulated responses of evapotranspiration and runoff to changes in the leaf area index of a Larix principis-rupprechtii plantation. Acta Ecol. Sin. 2016, 36, 6928–6938. [Google Scholar]
  75. Juhász, Á.; Hrotkó, K. Comparison of the transpiration part of two sources evapotranspiration model and the measurements of sap flow in the estimation of the transpiration of sweet cherry orchards. Agric. Water Manag. 2014, 143, 142–150. [Google Scholar] [CrossRef]
  76. Benyon, R.G.; Theiveyanathan, S.; Doody, T.M. Impacts of tree plantations on groundwater in south-eastern Australia. Aust. J. Bot. 2006, 54, 181–192. [Google Scholar] [CrossRef]
  77. Almeida, A.C.; Soares, J.V.; Landsberg, J.J.; Rezende, G.D. Growth and water balance of Eucalyptus grandis hybrid plantations in Brazil during a rotation for pulp production. For. Ecol. Manag. 2007, 251, 10–21. [Google Scholar] [CrossRef]
  78. Tian, A.; Wang, Y.H.; Webb, A.A.; Liu, Z.; Ma, J.; Yu, P.; Wang, X. Water yield variation with elevation, tree age and density of larch plantation in the Liupan Mountains of the Loess Plateau and its forest management implications. Sci. Total Environ. 2021, 752, 141752. [Google Scholar] [CrossRef] [PubMed]
  79. Zhang, J.J.; Liu, W.Z.; Han, X.Y. Characteristics of evapotranspiration of winter wheat in Changwu Plateau, Shaanxi. Jiangsu Agric. Sci. 2013, 41, 72–75. [Google Scholar]
  80. Juul, L.; Milena, H.; Jacobs, C.M.J.; Sjoerd, E.A.T.M.; Karofeld, E.; Berendse, F. How does tree density affect water loss of peatlands? A mesocosm experiment. PLoS ONE 2014, 9, e91748. [Google Scholar]
  81. Li, X.C.; Jiang, T.; Wen, S.S.; Wang, Y.J.; Qiu, X.F. Analysis of Temporal and Spatial Changes and Influencing Factors of Actual Evapotranspiration in the Pearl River Basin. J. Trop. Meteorol. 2014, 30, 483–494. [Google Scholar]
Figure 1. Location of micro-lysimeters.
Figure 1. Location of micro-lysimeters.
Water 13 01957 g001
Figure 2. Relationship between diameter and sapwood area at the beginning of study in 2019.
Figure 2. Relationship between diameter and sapwood area at the beginning of study in 2019.
Water 13 01957 g002
Figure 3. Dynamics of leaf area index during apple tree growth period in 2019 and 2020.
Figure 3. Dynamics of leaf area index during apple tree growth period in 2019 and 2020.
Water 13 01957 g003
Figure 4. Dynamics of daily precipitation and reference evapotranspiration of apple orchards during the growth period in 2019 and 2020.
Figure 4. Dynamics of daily precipitation and reference evapotranspiration of apple orchards during the growth period in 2019 and 2020.
Water 13 01957 g004
Figure 5. Relationship between LAI and transpiration (T) or evapotranspiration (ET) in 2019 and 2020.
Figure 5. Relationship between LAI and transpiration (T) or evapotranspiration (ET) in 2019 and 2020.
Water 13 01957 g005
Figure 6. The relationship between leaf area index of apple trees and crop coefficient Kc during the growth period in 2019 and 2020.
Figure 6. The relationship between leaf area index of apple trees and crop coefficient Kc during the growth period in 2019 and 2020.
Water 13 01957 g006
Table 1. The canopy structure of apple trees.
Table 1. The canopy structure of apple trees.
Tree NumberTree Age
(Years)
Canopy Diameter
(cm)
Canopy Height
(cm)
Stem Height
(cm)
Diameter
(cm)
Canopy Volume
(m3)
Land Area
(m2)
17401.0340.075.010.828.612.6
25317.5320.056.06.716.97.9
37442.5310.073.012.431.815.4
46315.0220.079.07.311.47.8
59458.5216.090.011.223.816.5
67399.5231.081.013.019.312.5
77401.5270.077.012.226.614.3
87413.5282.083.014.128.716.2
97420.1295.082.013.829.417.1
Table 2. Variation of monthly precipitation (P) and monthly reference evapotranspiration (ET0) of apples from April to September in 2019 and 2020.
Table 2. Variation of monthly precipitation (P) and monthly reference evapotranspiration (ET0) of apples from April to September in 2019 and 2020.
YearIndiceAprilMayJuneJulyAugustSeptemberSum
2019P (mm)65.58.942.490.5101.191.1399.5
ET0 (mm)125.5131.9133.8148.1151.9138.7829.9
2020P (mm)16.310.714.258.3245.267.2411.9
ET0 (mm)114.2125.0135.6148.8192.7126.3842.6
Table 3. Variation of monthly evapotranspiration of apples from April to September in 2019 and 2020.
Table 3. Variation of monthly evapotranspiration of apples from April to September in 2019 and 2020.
YearMonthE (mm)T (mm)Ic (mm)ET (mm)Percentage (%)
2019429.14.02.535.68.4
518.210.30.929.36.9
623.520.53.747.711.2
740.249.411.2100.823.7
852.357.212.8122.328.8
940.039.29.989.021.0
Total203.3180.641.0424.8100.0
2020421.93.61.226.75.8
516.411.61.029.06.3
626.724.90.952.511.3
740.647.67.095.220.6
886.375.523.8175.637.9
935.740.08.484.218.2
Total227.5203.242.3463.2100.0
Note: The percentage (%) indicates monthly ET/the total ET from April to September.
Table 4. Variation of crop coefficient Kc under different soil moisture.
Table 4. Variation of crop coefficient Kc under different soil moisture.
YearMonthRainfall (mm)ETD (mm)ETI (mm)ET0 (mm)KcIKsKc
2019465.535.654.0125.50.430.66 0.28
58.929.355.4131.90.420.53 0.22
642.447.766.3133.80.500.72 0.36
790.5100.8109.6148.10.740.92 0.68
8101.1122.3106.3151.90.701.15 0.81
991.189.183.2138.70.601.07 0.64
Total399.5424.8474.8829.9///
2020416.326.746.0114.20.400.58 0.23
510.729.051.7125.00.410.56 0.23
614.252.579.5135.60.590.66 0.39
758.395.2117.6148.80.790.81 0.64
8245.2175.6188.8192.70.980.93 0.91
967.284.285.9126.30.680.98 0.67
Total411.9463.2569.6842.6// /
Note: ETD and ETI represent evapotranspiration under water stress and irrigation, respectively.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jia, Q.; Wang, Y.-P. Relationships between Leaf Area Index and Evapotranspiration and Crop Coefficient of Hilly Apple Orchard in the Loess Plateau. Water 2021, 13, 1957. https://doi.org/10.3390/w13141957

AMA Style

Jia Q, Wang Y-P. Relationships between Leaf Area Index and Evapotranspiration and Crop Coefficient of Hilly Apple Orchard in the Loess Plateau. Water. 2021; 13(14):1957. https://doi.org/10.3390/w13141957

Chicago/Turabian Style

Jia, Qiong, and Yan-Ping Wang. 2021. "Relationships between Leaf Area Index and Evapotranspiration and Crop Coefficient of Hilly Apple Orchard in the Loess Plateau" Water 13, no. 14: 1957. https://doi.org/10.3390/w13141957

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