Three-dimensional numerical simulation of flow in trapezoidal cutthroat flumes based on FLOW-3D

Danjie RAN, Wene WANG, Xiaotao HU

Front. Agr. Sci. Eng. ›› 2018, Vol. 5 ›› Issue (2) : 168-176.

PDF(1213 KB)
Front. Agr. Sci. Eng. All Journals
PDF(1213 KB)
Front. Agr. Sci. Eng. ›› 2018, Vol. 5 ›› Issue (2) : 168-176. DOI: 10.15302/J-FASE-2018217
RESEARCH ARTICLE
RESEARCH ARTICLE

Three-dimensional numerical simulation of flow in trapezoidal cutthroat flumes based on FLOW-3D

Author information +
History +

Abstract

To solve the common problem of flumes flow-measurement accuracy without sacrificing water head, a new type of trapezoidal cutthroat flume to measure the discharge in terminal trapezoidal channels is presented. Using the computational fluid dynamic method, three-dimensional flow fields in trapezoidal cutthroat flumes were simulated using the RNG k-ε three-dimensional turbulence model along with the TruVOF technique. Simulations were performed for 12 working conditions, with discharges up to 0.075 m3·s1 to determine hydraulic performance. Experimental data for the trapezoidal cutthroat flume in terminal trapezoidal channel were also obtained to validate the simulation results. Velocity distribution of the flume obtained from simulation analyses were compared with observed results based on time-averaged flow field and comparison yielded a solid agreement between results from the two methods, with relative error below 10%. The results indicated that the Froude number and the longitudinal average velocity increased along the convergence section and decreased in the divergent section. In the upper throat, the Froude number was less than 0.5, which meets the water measurement requirement, and the critical flow appeared near the throat section. The maximum water head loss of the trapezoidal cutthroat flume was less than 9% of the total head, compared to the rectangular cutthroat flume, and head loss of trapezoidal cutthroat flume was significantly less. Regression models developed for upstream depth versus discharge under different working conditions were satisfactory, with a relative error of less than 2.06%, which meets the common requirements of flow measurement in irrigation areas. It was concluded that trapezoidal cutthroat flumes can improve flow-measurement accuracy without sacrificing water head.

Keywords

distribution of velocity / Froude number / head loss / stage-discharge relationship / trapezoidal cutthroat flume

Cite this article

Download citation ▾
Danjie RAN, Wene WANG, Xiaotao HU. Three-dimensional numerical simulation of flow in trapezoidal cutthroat flumes based on FLOW-3D. Front. Agr. Sci. Eng., 2018, 5(2): 168‒176 https://doi.org/10.15302/J-FASE-2018217

1 Introduction

The arid region of north-west China receives abundant solar radiation and is suitable for viticulture. However, there is little plant-available water in this region and agricultural water use has been poorly managed. Moreover, most vineyards are furrow-irrigated and this makes grape production unsustainable[13]. Use of drip irrigation has increased in recent years. Population growth, increased food demand and reduced water resources have all increased the importance of effective water management in arid and semiarid areas. Thus, sustainable agriculture in the region requires water-saving management and efficient irrigation methods[46].
Knowledge of evapotranspiration (ET), evaporation (E) and transpiration (T), and of partitioning ET into E and T is of fundamental importance for optimizing water management in water-limited environments[7,8]. Moreover, T is beneficial water consumption or productive water use because T is associated with plant productivity whereas E is lost water (unproductive) because it does not directly contribute to plant development[7,9]. The soil surface in vineyards with widely-spaced rows forms a large proportion of total vineyard surface area and this makes E/ET relatively high[10,11]. Thus, in sparse vegetation, irrigation methods are a significant factor to consider in E. Lascano et al.[12] found T/ET to be only 0.23 in a flood-irrigated vineyard. Zhao et al.[13] reported E/ET values of 46%–62% and T/ET of 40%–43% in a furrow-irrigated vineyard. However, E/ET in a drip-irrigated vineyard was reported to be as low as 0.09, whereas T/ET reached nearly 0.91[11]. A high E/ET is unfavorable for the conservation of agricultural water, and decreasing E is also important in the optimization of crop water use efficiency[14]. Analysis of ET partitioning under different irrigation methods, especially in sparsely vegetated production systems such as vineyards, can indicate how much higher T/ET is under drip irrigation than under furrow irrigation. Adoption of the method will lead to increased plant accumulation of dry matter and crop water use efficiency.
The crop coefficient (Kc) can be partitioned into the evaporation coefficient (Ke) and the basal crop coefficient (Kcb) (Kc=Ke+Kcb)[15,16], both of which also reflect the ET partitioning. Numerous previous studies have demonstrated the effects of crop growth and canopy development on both Kc and Kcb[17,18]. However, few studies have examined the factors that influence Ke, especially in the case of sparse vegetation[13] and the influence of different irrigation methods on the response of Ke to the primary controllers. Thus, there is a need to identify the behavior of the crop coefficients using different irrigation methods in sparsely vegetated vineyards.
Energy partitioning is a key component in energy balance processes within the soil-canopy-atmosphere continuum and studying it will benefit the optimization of agricultural water use[19,20]. Various factors influence the energy partitioning of field production systems. Previous studies on vineyards show that leaf area index (LAI) and soil water content (SWC) were the two key controlling factors[20,21]. Air temperature, wind speed, vapor pressure deficit and other environmental factors, including irrigation, have also been reported to be correlated with field production system energy partitioning[10,20], as well as in other systems[2224]. Only a few studies have investigated the responses of energy partitioning to the LAI and SWC in production systems with sparse canopies subjected to different irrigation methods. Furthermore, in a sparsely vegetated production system with relatively high E/ET, investigation of LEs/Rn [ratio of latent heat flux (LE) by soil E to net radiation (Rn)], seems pivotal because it recognizes the transformation of radiant energy to soil latent heat as part of overall energy flux. However, few studies have investigated the factors affecting LEs/Rn in vineyards and under drip irrigation.
Here, we have compared ET and energy partitioning and their responses to the main controllers in two vineyards, one with furrow and the other with drip irrigation. The aims were to (1) investigate the differences in ET partitioning between the vineyards from the perspective of E/ET and T/ET, Kcb and Ke in quantity; (2) compare the responses of the ratios concerning ET partitioning to the primary controllers of the two vineyards; and (3) examine whether energy partitioning ratios respond differently to the biotic (LAI) and abiotic (SWC) factors under the two irrigation methods in order to determine a scientific basis for an irrigation program that optimizes crop water use efficiency.

2 Materials and methods

2.1 Study sites

The study sites are in the temperate zone in an arid region of north-west China (37°51′ N, 102°53′ E, 1585 m asl). The mean annual precipitation is 164.4 mm, annual mean temperature 8°C and total annual sunshine ˃ 3000 h. The region is very short of water, with the groundwater table below 25 m under the ground surface and the annual mean pan evaporation 1926 mm. The experiment was conducted at the Shiyanghe Experimental Station of China Agricultural University in the Shiyang River Basin. The early-ripening Merlot (Vitis vinifera) vineyard at Huangtai Station (which was studied in 2015) was 1650 m long by 1400 m wide. The vines were planted in 1999 in east–west rows. Row spacing was 2.7 m and plant spacing was 1 m. A vertical trellis system of three wires 0.5, 1.0 and 1.5 m high was used. The vines were trained in a bilateral cordon. The soil is a sandy loam with an average bulk density of 1.49 g·cm3 in the top 1 m of the soil profile. The saturated water content of the soil was 0.45 cm3·cm3. The field capacity and wilting point were 0.31 and 0.11 cm3·cm3 in the top 1 m of the profile. The vineyard was furrow-irrigated five times during the growing season on April 22, May 27, June 27, July 27, and August 30 with 70 mm applied on each occasion. The early-ripening Pinot Noir (V. vinifera) vineyard at Mogao Station was studied in 2017. The vines were planted in 2014 in north–south rows. Row spacing was 3.3 m and plant spacing was 0.5 m. A vertical trellis system of three wires 0.4, 0.7 and 1.0 m high was used. The vines were trained in a single cordon. The soil is a loamy clay with an average bulk density of 1.67 g·cm3 and the saturated SWC was 0.38 cm3·cm3. Field capacity was 0.22 cm3·cm3 in the top 1 m of the soil profile. The vineyard was drip-irrigated from 40 cm above the ground. The drip emitters had a discharge rate of 3 L·h1 with 40 cm between emitters. The vineyard was irrigated six times during the growing season on May 29, June 23, July 3, August 2 and 30, and September 16 with a total irrigation of 226 mm. The two vineyards are about 2 km apart.

2.2 Eddy covariance measurement and data correction

ET was measured in both vineyards using the same eddy covariance (EC) system in the growing seasons of 2015 and 2017. The instruments were mounted 4 m above ground level in the south-east of the vineyard and met the requirement for adequate fetch so that the ratio of the fetch to the height of installation was ˃ 100:1. The EC equipment consisted of a CSAT3 three-dimensional sonic anemometer (Campbell Scientific, Logan, UT), an LI-7500 open-path CO2 and H2O gas analyzer (Li-Cor Inc., Lincoln, NE), an NR-Lite net radiometer (Kipp and Zonen, Delft, Netherlands), an HMP45C temperature and humidity probe (Vaisala Inc., Helsinki, Finland), HFP01 soil heat flux (G) plates (Hukseflux, Delft, Netherlands), and a CR5000 data logger (Campbell Scientific). Wind speed, ultrasonic virtual temperature, and air and water vapor densities were measured every 0.1 s to obtain both latent and sensible heat observations. In 2015, two G plates were installed 5 cm below the soil surface in the ditch and on the ridge. In 2017, five G plates were installed 8 cm below in the soil surface at five points across the interrow. The sensors were connected to the CR5000 data logger and a value averaged over 30 min was recorded.
The correction of the raw EC data was conducted using EddyPro 4.2.1 (Li-Cor Inc.) in express mode. The correction comprised 30-min block averaging, time-lag compensation, statistical tests[25], WPL density correction[26], the planar fit coordination rotation[27,28], correction for the ultrasonic virtual temperature[29], spectral correction[30,31], data quality control flags[32] and footprint analysis[33]. EC data, after correction by EddyPro, were flagged with 0 (highest quality), 1 and 2 (lowest quality). LE and sensible heat flux (H) data with flag 2 were discarded. Data gaps caused by instrument problems, power failure, data filtration and correction were filled. Short gaps (≤2 h) were filled via linear interpolation and longer gaps (>2 h) using multiple linear regressions for the main controlling factors, Rn, LAI and SWC[34].
The closure of the measured energy budget component for the daytime EC data was forced using the β forced closure method in accordance with the assumption that the β was correctly measured by the EC system (the daytime and nighttime were identified through T/F (daytime/nighttime) flags obtained from EC data after correction via EddyPro). Assuming that H was accurately measured, the residual-LE closure method was used for the nighttime EC data, and residual LE was calculated by the energy balance equation LE=RnGH[35]. Numerous studies indicate the value of this method of data correction[13,36,37]. ET was calculated from LE after energy balance closure as LE=λET, where l is the latent heat of vaporization.
Energy balance closure has customarily been used to evaluate the accuracy of EC measurements through linear regression between turbulent fluxes (LE+H) and available energy (RnG)[35,38]. Linear regression using 0.5 h EC data filtered and corrected by EddyPro gave gradients of 1.10 and 0.84, respectively, for the 2015 and 2017 regression equations and the corresponding coefficients of determination (R2) were 0.91 and 0.90. The values were within the common results from previous studies[38,39].

2.3 Soil evaporation and transpiration

Daily E was measured using PVC tube microlysimeters (MLs) with an internal diameter of 10 cm and height of 20 cm. In 2015, three replicates of six MLs were used. They were pressed into the soil at distances of 0.05, 0.55, 1.05, 1.55, 2.05 and 2.55 m from the row center along a cross-section of the interrow from south to north. In 2017, six replicates of five MLs were installed at distances of 0.2, 0.5, 1.65, 2.8 and 3.1 m from the row center along a cross-section of the interrow from west to east. The tops of the MLs were deployed level with the soil surface and the MLs were weighed daily at 19:00 local time with an electronic scale with a precision of 0.1 g[13]. Daily vineyard E was calculated for both vineyards using the mean value of the daily E measured from the MLs.
A sap flow system (Flow32-1K, Dynamax, Houston, TX) was used to determine vine T. Three vines (2015) and eight vines (2017) were randomly selected in a circle with a radius of 15 m centered on the EC system. The types of sensors used were SGB 25 in 2015 and SGA 10 and 13 and SGB 16 and 25 in 2017. In 2015, the sensors were disconnected before irrigation and were reinstalled 2–3 d after furrow irrigation to prevent damage to the sensors from flooding[3,40]. Sap flow was measured every 60 s and the average over 15 min was recorded with a CR1000 data logger (Campbell Scientific, Logan, UT). Using the average ground area of each vine the sap flow (L·d1) was scaled to vine T (mm·d1) and the mean value for the monitored vines was used to calculate vineyard T[41].

2.4 Leaf area index, topsoil soil water content and meteorological parameters

LAI in 2015 was estimated using the relationship between shoot length and the total leaf area per shoot described by Ortega-Farias et al.[39] for all shoots on the vine. Branches were randomly selected and the total leaf area per shoot was measured with an AM300 portable leaf area meter (ADC BioScientific Ltd., Hoddesdon, UK) at the early and intermediate growth stages. The total shoot length per vine of the selected vines was measured every 7–10 d during the experimental period. The LAI of the vineyard was calculated from the average leaf area of the selected vines and plant spacing. Details of the calculation of LAI in the furrow-irrigated vineyard are given by Zhao et al.[40]. The LAI of the drip-irrigated vineyard was measured with an LAI-2200C plant canopy analyzer (Li-Cor Bioscience Inc., Lincoln, NE) on clear days between 18:30 and 19:30 local time every 5–7 d in 2017 following the instructions for row crops[42]. Six sets of measurements were undertaken during each LAI measurement. Each set of measurements consisted of seven readings with one reading above the canopy and the other six readings below the canopy taken at consecutive positions with the line sensor parallel to the row direction, running from one row center to another adjacent row center. An estimate of seasonal canopy development including unmeasured LAI days was obtained by linear regression between successive observations.
The topsoil volumetric SWC (10 cm depth) was continuously monitored using ECH2O soil moisture sensors (Decagon Devices, Inc., Pullman, WA) connected to EM50 data loggers (Decagon Devices, Inc.). Six (2015) and 15 (2017) sensors were used and SWC was calculated every 10 min. The mean value of the measured data was used to calculate topsoil SWC in each vineyard.
The meteorological parameters (air temperature, relative humidity, wind speed, precipitation, and solar radiation) were measured with a Hobo automatic weather station (Onset Computer Corporation, Bourne, MA, USA) located near the vineyards 2 m above the soil surface. Daily reference evapotranspiration (ETo) was calculated using the FAO Penman-Monteith equation[15]. The monthly average values of the main meteorological parameters during the 2015 and 2017 growing seasons are shown in Table 1, and Fig. 1 shows their seasonal dynamics. According to the average seasonal standardized precipitation index which was recommended for tracking meteorological drought, both of the growing seasons were classified as non-drought relative to the local climatic conditions[43].
Tab.1 Monthly mean air temperature (Ta), solar radiation (Rs), vapor pressure deficit (VPD), monthly total precipitation (P), soil water content (SWC, 10 cm depth), and leaf area index (LAI) in the furrow-irrigated vineyard studied in 2015 and the drip-irrigated vineyard studied in 2017
Year Month Ta/ºC Rs/(MJ·m–2·d–1) VPD/kPa P/mm SWC/(m3·m–3) LAI/(m2·m–2)
2015 Late Apr.–May 16.58 20.42 1.38 5.10 0.13 0.31
Jun. 19.27 19.84 1.17 43.00 0.15 1.53
Jul. 21.02 21.18 1.22 53.40 0.15 2.05
Aug. 20.25 19.66 1.30 8.60 0.11 1.55
Sep. 14.52 14.07 0.71 36.40 0.14 1.43
Late Apr.–Sep. 18.24 19.14 1.17 161.80 0.14 1.31
2017 Late Apr.–May 17.05 21.11 1.38 3.06 0.13 0.40
Jun. 20.20 21.24 1.38 21.80 0.13 1.06
Jul. 22.87 20.50 1.40 39.60 0.14 1.07
Aug. 20.62 17.39 0.99 44.40 0.17 1.19
Sep. 17.84 15.49 1.05 6.20 0.12 0.97
  Late Apr.–Sep. 19.67 19.19 1.24 133.40 0.14 0.93

2.5 Data analysis

To estimate the seasonal variation of E/ET, T/ET, Ke, Kcb and their magnitudes during each growing stage, gaps of E and T were filled. When the daily E was available and daily T was unavailable, T=ETE (the regression equation between the sums of observed E and T and the observed ET using EC was E+T=0.92ETfor the two growing seasons, R2 = 0.50, P<0.001). When the daily T was available and the daily E was unavailable, E=ETT. When both the daily E and T were unavailable, E was filled using the regression between the observed daily E and Rn and the SWC (10 cm), and T=ETE[13,44]. Missing E/ET values were calculated using interpolated E and T/ET=1E/ET. The interpolated values of E and T were used only to calculate the seasonal totals and the averages of different growing stages but were not used to analyze any correlation with the controlling factors.
The daily Kc, Ke and Kcb were calculated using values of daily ET, E, T and ETo (Kc=ET/ETo, Ke=E/ETo,Kcb=T/ETo)[15]. The analyses of the relationships between daily E, E/ET, Ke and the governing factors were undertaken using only observed data. The daily ET, Rn, LE, H and G were calculated as the sums of the 30-min values. Daily b was calculated as the ratio of average daily H to daily LE.
LAI was normalized to RLAI (RLAI=LAI/LAImax) to better reflect the growing stage and facilitate the analysis of the relationships between ratios concerning energy partitioning and LAI. SWC was normalized to RSWC (RSWC=SWC/SWCsat,SWCsat, i.e., the volumetric saturated SWC) to better reflect the SWC status of the two vineyards. Two-tailed significance tests were conducted using the SPSS 25.0 software package (SPSS Inc., Chicago, IL, USA).
Fig.1 Seasonal variation in daily (a, b) mean air temperature (Ta), (c, d) mean solar radiation (Rs), (e, f) mean vapor pressure deficit (VPD), and (g, h) total precipitation (P) in the furrow-irrigated vineyard studied in 2015 (left) and the drip-irrigated vineyard studied in 2017 (right).

Full size|PPT slide

3 Results and discussion

3.1 Evapotranspiration partitioning

Figure 2 shows the seasonal variation in daily ETo, ET and LAI for the 2015 and 2017 growing seasons together with daily precipitation and irrigation. The monthly averages of these parameters are listed in Table 2. The sums of precipitation and irrigation were 512 mm (2015) and 359 mm (2017). Cumulative ETo for the growing season was 611 mm (furrow irrigation) and 609 mm (drip irrigation) and cumulative ET was 404 mm (furrow irrigation) and 339 mm (drip irrigation). The seasonal trends of daily ET were similar for the two vineyards, maintaining a relatively high value in the middle of the season and low values early and late in the season.
Fig.2 Seasonal variation in daily evapotranspiration (ET) measured by eddy covariance and reference evapotranspiration (ETo) in the furrow-irrigated vineyard studied in 2015 (left) and the drip-irrigated vineyard studied in 2017 (right). Leaf area index (LAI), irrigation (I) and precipitation (P) are shown.

Full size|PPT slide

Tab.2 Monthly mean values of reference evapotranspiration (ETo), evapotranspiration (ET), evaporation (E), transpiration (T), E/ET, T/ET, crop coefficient (Kc), soil evaporation coefficient (Ke), and basal crop coefficient (Kcb) during the growing season in the furrow-irrigated vineyard studied in 2015 and the drip-irrigated vineyard studied in 2017
Year Month ETo/(mm·d–1) ET/(mm·d–1) E/(mm·d–1) T/(mm·d–1) E/ET T/ET Kc Ke Kcb
2015 Late Apr.–May 4.23 1.42 0.96 0.53 0.56 0.44 0.34 0.21 0.15
Jun. 3.97 2.94 1.62 1.19 0.50 0.50 0.74 0.38 0.31
Jul. 4.16 3.96 2.07 1.59 0.49 0.51 0.92 0.47 0.38
Aug. 3.73 2.33 0.71 1.30 0.31 0.69 0.62 0.20 0.34
Sep. 2.46 2.08 1.26 0.74 0.62 0.38 0.82 0.51 0.28
Late Apr.–Sep. 3.75 2.48 1.30 1.04 0.50 0.50 0.67 0.35 0.29
2017 Late Apr.–May 4.58 1.30 0.99 0.55 0.49 0.51 0.30 0.22 0.28
Jun. 4.31 2.47 1.07 1.58 0.41 0.59 0.57 0.25 0.36
Jul. 4.34 2.77 0.98 1.52 0.31 0.69 0.57 0.20 0.35
Aug. 3.42 3.14 1.66 1.79 0.50 0.50 0.89 0.48 0.55
Sep. 2.80 1.60 0.61 1.43 0.45 0.55 0.55 0.26 0.54
  Late Apr.–Sep. 3.91 2.23 1.06 1.36 0.42 0.58 0.57 0.29 0.43
Figure 3 shows the seasonal change in the proportional ET partitions, E/ET and T/ET for both vineyards. Table 2 shows the monthly averages of these ratios. For both vineyards E/ET was higher in the early and late growth stages, and lower in the middle stage with T/ET showing the opposite trend. T was low in the early growth stage because of the small canopy, with a large proportion of bare soil resulting in E/ET being higher than T/ET. Canopy cover (CC), increased as the vines grew and T increased, leading to an increase in T/ET. Leaf senescence in the later part of the growing season reduced daily T and increased the proportion of bare soil, thus E/ET increased and T/ET decreased. During the first half of the growing season in the furrow-irrigated vineyard (from May 20 to July 22, 2015) E/ET was higher because of greater precipitation[13]. E/ET was higher in the drip-irrigated vineyard from late July to mid-August because of intense rain events which resulted in a greater SWC, and hence led to a higher E. The average E/ET values of the vineyards were 0.50 (furrow irrigation) and 0.42 (drip irrigation), with corresponding T/ET values of 0.50 and 0.58. Notably, the seasonal average E/ET decreased and average T/ET increased under drip irrigation. López-Urrea et al.[45] also reported that the irrigation method affected E/ET. Yunusa et al.[46] found that seasonal average E/ET was higher (0.46–0.51) than average T/ET (0.18–0.19) (T/ET of cover crops ranged between 0.30 and 0.36) for a furrow-irrigated vineyard. However, in a drip-irrigated vineyard in the same region, Yunusa et al.[47] found lower seasonal average E/ET (0.41) and higher T/ET (0.55). Moreover, for drip-irrigated vineyards E/ET was found to range from 0.09 to 0.31, with T/ET ranging from 0.69 to 0.91 as reported by Kool et al.[11], Teixeira et al. [16], Poblete-Echeverría et al[48]., Poblete-Echeverría and Ortega-Farias[49]. Kool et al.[7] reported seasonal average E/ET of 0.41±0.21 and T/ET of 0.57±0.21 for vineyards, and seasonal average E/ET of 0.30±0.12 and T/ET of 0.69±0.13 for drip-irrigated vineyards. In sparsely vegetated vineyards, E is a relatively large component of ET because of the large proportion of bare soil in the interrow. It is therefore essential to choose the irrigation method that best reduces unproductive water consumption in sparse vegetation.
Fig.3 Seasonal variation in daily evapotranspiration (ET) partitioning (E/ET, T/ET) in the furrow-irrigated vineyard studied in 2015 (upper) and the drip-irrigated vineyard studied in 2017 (lower). E/ET, the ratio of soil evaporation (E) to ET; and T/ET, the ratio of transpiration (T) to ET.

Full size|PPT slide

Figure 4(a) shows that E/ET was linearly related to RSWC in production systems under a sparse canopy and, moreover, E/ET was higher in the furrow-irrigated vineyard than in the drip-irrigated vineyard with the same RSWC which can be explained by the differential responses of E to RSWC under the two irrigation methods as shown in Fig. 4(b). For any given RSWC, the proportion of the soil surface that was wet was greater in the furrow-irrigated vineyard and thus E increased due to the effect of solar radiation.
Fig.4 Relationships between (a) E/ET and (b) E and relative soil water content (RSWC, 10 cm depth) in the furrow-irrigated vineyard studied in 2015 and the drip-irrigated vineyard studied in 2017. E, evaporation; and ET, evapotranspiration. ***, P<0.001.

Full size|PPT slide

It is concluded that when considering total seasonal water use, more water is used productively with T, which indicates increased plant accumulation of dry matter and increased water use efficiency in a drip-irrigated vineyard than in a furrow-irrigated vineyard. Drip irrigation also decreased the proportion of unproductive water use (i.e., E) which also reduced the non-beneficial water consumption.

3.2 Crop coefficient and dual crop coefficient

Figure 5 shows seasonal variation in the daily Kc (Fig. 5(a) and 5(b)), the daily Ke and SWC (10 cm deep) (Fig. 5(c) and 5(d)), and the daily Kcb (Fig. 5(e) and 5(f)) during the two growing seasons in the furrow-irrigated vineyard in 2015 and the drip-irrigated vineyard in 2017. The monthly averages of the coefficients are shown in Table 2. Seasonal variation in Kc and Kcb coincided with canopy development. Kc increased abruptly in early September 2015 and in the period from late-July to mid-August 2017. The increases may be attributed to the favorable conditions for E and T during these periods as a result of more irrigation and precipitation events. The seasonal variation in Ke corresponded to the fluctuation in SWC. The average Kc for the two vineyards was 0.67 (furrow-irrigated) and 0.57 (drip-irrigated). The average Ke was 0.35 (furrow-irrigated) and 0.29 (drip-irrigated) and the average Kcb was 0.29 (furrow-irrigated) and 0.43 (drip-irrigated). Seasonal average Ke was larger than average Kcb in the furrow-irrigated vineyard but showed the opposite trend in the drip-irrigated vineyard. Lascano et al.[12] found that, due to high E, average Ke was up to 0.46 in a flood-irrigated vineyard. Teixeira et al.[16] reported an average range of Kcb of 0.66–0.69 but a much lower average Ke of 0.08–0.09 in a drip-irrigated vineyard in Brazil. Zhao et al.[13] reported that maximum Ke was reached after precipitation but not after irrigation as the wetted fraction of the soil was about 50% after irrigation but nearly 100% after precipitation in a furrow-irrigated vineyard. Yunusa et al.[47] found a Ke of 0.03 during a dry period (without rainfall or irrigation) and 0.26 during a wet period (irrigated days), indicating that the increased surface area of the wetted fraction of the soil, which can vary according to irrigation method, significantly affected the proportion of E to total water consumption.
Fig.5 Seasonal variation in daily (a, b) crop coefficient (Kc), (c, d) evaporation coefficient (Ke) and topsoil soil water content (SWC, 10 cm depth), (e, f) basal crop coefficient (Kcb) in the furrow-irrigated vineyard studied in 2015 (left) and the drip-irrigated vineyard studied in 2017 (right).

Full size|PPT slide

Figure 6 shows the relationship between daily vineyard Kc, Ke and RSWC in the furrow-irrigated vineyard in 2015 and the drip-irrigated vineyard in 2017. Both Kc and Ke increased linearly as RSWC increased and Ke was more closely correlated with RSWC than Kc. Notably, Kc was greater in the furrow-irrigated vineyard than in the drip-irrigated vineyard at any given RSWC, largely due to the greater Ke in the furrow-irrigated vineyard than in the drip-irrigated vineyard at the same RSWC. This is explained by the larger soil wetted fraction in the furrow-irrigated vineyard. The differential responses of Kc and Ke to RSWC under the two irrigation methods also indicates the effectiveness of a suitable irrigation method in saving water in irrigated vineyards[50].
Fig.6 Relationships between (a) daily crop coefficient (Kc) and (b) evaporation coefficient (Ke) and relative soil water content (RSWC, 10 cm depth) in the furrow-irrigated vineyard studied in 2015 and the drip-irrigated vineyard studied in 2017. ***, P<0.001.

Full size|PPT slide

3.3 Energy partitioning

Figure 7 shows the seasonal variation in daily vineyard proportional energy partitioning terms LE/Rn, H/Rn, G/Rn, evaporative fraction (EF), and b in the furrow-irrigated vineyard in 2015 and the drip-irrigated vineyard in 2017. The monthly averages of these ratios are given in Table 3. Seasonal LE/Rn and EF maintained relatively high values in the middle of the season and low values early and late in the season whereas H/Rn and b showed the opposite trend. Seasonal G/Rn exhibited no obvious variation. Seasonal average LE/Rn was 0.57 (furrow-irrigated) and 0.46 (drip-irrigated) and average H/Rn was 0.26 (furrow-irrigated) and 0.36 (drip-irrigated). Seasonal average EF (EF= LE / (Rn – G)) was 0.69 (furrow-irrigated) and 0.57 (drip-irrigated) and average b was 0.63 (furrow-irrigated) and 0.97 (drip-irrigated). Seasonal average G/Rn was 0.17 in both vineyards. The minimum monthly average H/Rn and b and the maximum monthly average LE/Rn and EF occurred in July (furrow-irrigated) and August (drip-irrigated) when the vines were at the high water-consuming growth stages of berry development and veraison. Greater seasonal average LE/Rn and EF values with lower H/Rn and b were observed in the furrow-irrigated vineyard than in the drip-irrigated vineyard. Teixeira et al.[16] reported an LE/Rn of 78%–83% and an H/Rn of 18% in two growing seasons in a drip-irrigated vineyard in Brazil, and a higher LE/Rn of 87%–91% and lower H/Rn of 12% in a microsprinkler-irrigated vineyard due to a more moist microclimate and soil conditions which were suitable for ET. Yunusa et al.[47] reported that H/Rn was 55% during a six-day warm dry period and 49% during a six-day cool humid period in a drip-irrigated vineyard in Austria and that EF was higher during cool humid days because of the greater availability of soil water. Kool et al.[11] reported an H/Rn of 42%, LE/Rn of 44%, and G/Rn of 3% in a drip-irrigated vineyard. The results further indicate that energy partitioning differed under different moisture conditions which could be induced by different irrigation methods, and drip irrigation resulted in a lower proportion of LE and higher proportion of H to the available energy than in the furrow-irrigated vineyard.
Fig.7 Seasonal variation in the proportional daily energy partitioning H/Rn, LE/Rn, G/Rn, EF, and b in the furrow-irrigated vineyard studied in 2015 (left) and the drip-irrigated vineyard studied in 2017 (right). LE/Rn, H/Rn, and G/Rn are the partitioning of net radiation (Rn) into latent heat flux (LE), sensible heat flux (H) and soil heat flux (G). EF, evaporative coefficient; and b, Bowen ratio.

Full size|PPT slide

Tab.3 Monthly mean energy fluxes of net radiation (Rn), latent heat flux (LE), sensible heat flux (H), and soil heat flux (G), proportional energy partitioning H/Rn, LE/Rn, G/Rn, evaporative fraction (EF) and Bowen ratio (β) in the furrow-irrigated vineyard studied in 2015 and the drip-irrigated vineyard studied in 2017
Year Month Rn/(MJ·m–2·d–1) LE/(MJ·m–2·d–1) H/(MJ·m–2·d–1) G/(MJ·m–2·d–1) H/Rn LE/Rn G/Rn EF β
2015 Late Apr.–May 8.75 3.19 3.55 2.02 0.41 0.35 0.23 0.45 1.49
Jun. 10.25 6.65 2.10 1.51 0.20 0.65 0.15 0.76 0.34
Jul. 11.50 8.96 1.28 1.26 0.13 0.75 0.12 0.84 0.21
Aug. 9.23 5.28 2.40 1.55 0.27 0.55 0.18 0.69 0.52
Sep. 6.34 4.69 1.26 0.39 0.20 0.68 0.12 0.79 0.30
Late Apr.–Sep. 9.20 5.60 2.21 1.39 0.26 0.57 0.17 0.69 0.63
2017 Late Apr.–May 8.74 3.03 5.68 0.04 0.59 0.29 0.11 0.35 2.23
Jun. 10.12 5.54 3.40 1.19 0.32 0.51 0.17 0.64 0.62
Jul. 10.10 5.61 3.00 1.49 0.31 0.45 0.25 0.57 0.90
Aug. 9.14 7.11 1.25 0.78 0.16 0.64 0.20 0.81 0.26
Sep. 6.88 3.42 2.79 0.67 0.39 0.46 0.15 0.57 0.84
  Late Apr.–Sep. 9.00 4.91 3.27 0.81 0.36 0.46 0.17 0.57 0.97

3.4 Biotic and abiotic control of energy partitioning

Figure 8 shows the relationship between energy partitioning terms LE/Rn, H/Rn, EF, β and RLAI in the furrow-irrigated vineyard in 2015 and the drip-irrigated vineyard in 2017. There were strong linear correlations between LE/Rn, H/Rn, EF, β and RLAI. LE/Rn and EF increased as RLAI increased whereas H/Rn and β decreased with rising RLAI. At a given RLAI, LE/Rn and EF were lower whereas H/Rn and β were higher in the drip-irrigated vineyard than the corresponding values in the furrow-irrigated vineyard, indicating that a greater proportion of Rn was partitioned into LE in the furrow-irrigated vineyard. Numerous previous studies find a robust relationship between energy partitioning terms and LAI. Zhao et al.[20] showed that CC was a major factor determining energy partitioning and that LAI explained 41% of the seasonal variation in energy partitioning. Shen et al.[51] found a positive correlation between EF and LAI in a winter wheat-summer maize rotation in Luancheng, Hebei Province, China. Burba and Verma[22] reported a linear correlation between ET/Rn and LAI when there was no water stress and that ET/Rn increased as LAI increased in a tallgrass prairie and cultivated wheat production system in Oklahoma in the United States. A linear correlation between EF and LAI was reported by Li et al.[23] above a grazing steppe in central Mongolia. Hammerle et al.[24] reported that in a mountain grassland in Austria, both LE/Rn and H/Rn were closely correlated with LAI, and before LAI reached around 4–5 m2·m2 LE/Rn increased with increasing LAI and this was attributed to increasing topsoil conductance comprising contributions from both the soil surface and the canopy conductance, and H/Rn decreased with increasing LAI, possibly induced by the decreased surface relative to the air temperature because of the vigorously transpiring canopy. Qiu et al.[19] obtained similar results with greenhouse-grown hot peppers under drip irrigation in north-west China, with lower LE/Rn and higher H/Rn than the corresponding ratios under furrow irrigation and both LE/Rn and H/Rn were correlated with LAI. In the present study, the difference in LE/Rn between vineyards may be ascribed largely to the LE values. The lower soil wetted fraction in the vineyard under drip irrigation might give lower LEs/Rn and thus lower LE/Rn. Higher H/Rn under drip irrigation might result from the greater temperature difference between soil and air than in the vineyard receiving furrow irrigation.
Fig.8 Relationships between proportional daily energy partitioning LE/Rn, H/Rn, EF, β and relative leaf area index (RLAI) in the furrow-irrigated vineyard studied in 2015 and the drip-irrigated vineyard studied in 2017. LE/Rn and H/Rn are the partitioning of net radiation (Rn) into latent heat flux (LE) and sensible heat flux (H). EF, evaporative coefficient; and β, Bowen ratio. ***, P<0.001.

Full size|PPT slide

Figure 9 shows the relationships between daily vineyard proportional energy partitioning terms LE/Rn, LEs/Rn, H/Rn, EF, b and RSWC (10 cm deep) in the furrow-irrigated vineyard in 2015 and the drip-irrigated vineyard in 2017. LE/Rn, LEs/Rn and EF increased with increasing RSWC. At any given RSWC, LE/Rn, LEs/Rn and EF were lower under drip irrigation than under furrow irrigation, indicating that a smaller fraction of Rn was partitioned into LE and LEs in the drip-irrigated vineyard. H/Rn and b decreased as RSWC increased, and the main factor might be more moist soil conditions, leading to a smaller difference between the soil surface temperature and the air temperature. Furthermore, at any given value of RSWC, H/Rn and b were higher under drip irrigation possibly because the SWC of the bare soil in the interrow was lower and thus there was a greater discrepancy between the soil surface temperature and the air temperature (data not shown). Kool et al.[11] reported that LE/Rn and H/Rn responded strongly to irrigation in a drip-irrigated vineyard, with highest LE/Rn and lowest H/Rn values occurring on irrigated days followed by a decrease in LE/Rn and an increase in H/Rn, showing that SWC influenced energy partitioning. Baldocchi et al.[21] found that water status influenced stomatal behavior and thus regulated energy partitioning and that SWC explained 32% of the seasonal variation in energy partitioning. Zhao et al.[20] reported that SWC explained 32% of the seasonal variation in energy partitioning in a furrow-irrigated vineyard. Burba and Verma[22] reported a positive correlation between ET/Rn and SWC under water stress in a native tallgrass prairie and cultivated wheat production system in Oklahoma in central USA. They further suggested that the influence of SWC on ET/Rn decreased as CC developed. Li et al.[23] found that EF increased with increasing SWC and further demonstrated that the soil received most rainfall and returned to the atmosphere as E and T, indicating that energy partitioning was strongly correlated with SWC, an effect also observed in the present study. Hammerle et al.[24] reported that energy partitioning was influenced by SWC; low SWC induced stomatal closure, leading to a decrease in latent heat which was often compensated for by an increase in sensible heat. The present study also confirms the strong correlation between energy partitioning and SWC and further indicates that furrow and drip irrigation can produce differences in energy partitioning in sparsely vegetated vineyards.
Fig.9 Relationships between proportional daily energy partitioning LEs/Rn, LE/Rn, H/Rn, EF, b and relative soil water content (RSWC, 10 cm depth) in the furrow-irrigated vineyard studied in 2015 and the drip-irrigated vineyard studied in 2017. LEs/Rn, LE/Rn, and H/Rn are the partitioning of net radiation (Rn) into the latent heat flux by soil evaporation (LEs), latent heat flux (LE) and sensible heat flux (H). EF, evaporative coefficient; and b, Bowen ratio. ***, P<0.001.

Full size|PPT slide

4 Conclusions

Here, we have made continuous measurements of ET, E and T to compare ET and energy partitioning, crop coefficients, and the main factors that control them in a furrow-irrigated vineyard and a drip-irrigated vineyard. The results show that drip irrigation reduced seasonal average E/ET and increased average T/ET. Daily Kc and Ke were linearly correlated with RSWC. Kc and Ke were lower in the drip-irrigated vineyard than in the furrow-irrigated vineyard at the same RSWC. Seasonal average Kcb was larger than Ke under drip irrigation but showed the opposite trend under furrow irrigation. The findings indicate that the irrigation method influenced ET partitioning and drip irrigation was favorable in increasing water use efficiency in a sparsely vegetated vineyard.
LE/Rn, H/Rn, EF, and b were determined by RLAI and RSWC. Daily LE/Rn and EF were lower whereas H/Rn and b were higher in the drip-irrigated vineyard than in the furrow-irrigated vineyard at the same value of RLAI. The positive linear correlations between LE/Rn, LEs/Rn, EF and RSWC and the lower values of these ratios at the same RSWC in the drip-irrigated vineyard together indicate that drip irrigation resulted in less unproductive water use compared to furrow irrigation.
This study highlights the importance of using drip irrigation, especially in production systems with sparse canopies in arid regions. Longer observation periods in the same field using different irrigation methods in future studies would provide further insights into this important topic in arid agricultural environments.

References

[1]
Galán-Martín A, Vaskan P, Vallejo A, Esteller L J, Guillén-Gosálbez G. Multi-objective optimization of rainfed and irrigated agricultural areas considering production and environmental criteria: a case study of wheat production in Spain. Journal of Cleaner Production, 2017, 140(2): 816–830
[2]
Wang Y B, Liu D, Cao X C. Agricultural water rights trading and virtual water export compensation coupling model: a case study of an irrigation district in China. Agricultural Water Management, 2017, 180(Part A): 99–106
[3]
Valipour M. Increasing irrigation efficiency by management strategies: cutback and surge irrigation. Journal of Agricultural and Biological Science, 2013, 8(1): 35–43
[4]
Wang C D. Water measurement technique and measure.Beijing: Water and Power Press, 2005 (in Chinese)
[5]
Samani Z, Magallanez H. Closure to “Simple flume for flow measurement in open channel” by Zohrab Samani and Henry Magallanez. Journal of Irrigation and Drainage Engineering, 2002, 128(2): 129–131 doi:10.1061/(ASCE)0733-9437(2002)128:2(132)
[6]
Cone V M. The Venturi flume. Journal of Agricultural Research, 1917, 6(4): 115–129
[7]
Parshall R L. The improved Venturi flume. Transactions of the American Society of Civil Engineers, 1926, 89(1): 841–851
[8]
Skogerboe G V, Hyatt M L. Rectangular cutthroat flow measuring flumes. Proceedings of the American Society of Civil Engineers, 1967, 93(IR4): 1–13
[9]
Clemmens A J, Bos M G, Replogle J A. RBC broad-crested weirs for circular sewers and pipes. Journal of Hydrology, 1984, 68(1–4): 349–368
CrossRef Google scholar
[10]
Lv H X, Pei G X, Yang L X. Hydraulics.Beijing: Agriculture Press, 2011(in Chinese)
[11]
Parshall R L. Parshall measuring flume. Colorado Experiment Station Bulletin, 1936, 423
[12]
Hager W H. Modified Venturi channel. Journal of Irrigation and Drainage Engineering, 1985, 111(1): 19–35
CrossRef Google scholar
[13]
Jesson M, Sterling M, Baker D. Application of ISO4359 for discharge calculation in a narrow flume. Flow Measurement and Instrumentation, 2017, 54: 283–287
CrossRef Google scholar
[14]
Das R, Nayek M, Das S, Dutta P, Mazumdar A. Design and analysis of 0.127 m (5”) cutthroat flume. Ain Shams Engineering Journal, 2017, 8(3): 295–303
CrossRef Google scholar
[15]
Mazumdar A, Dutta P, Nayek M, . Calibration and Discharge Measurement Using 0.127 Meter (5”) Parshall Flume. Iahr World Congress, 2017
[16]
Jing S Y, Wang L, Du H, Wei G. Applications of FLOW-3D in numerical simulation of fluid-structure interaction. The 13th National Hydrodynamic Academic Conference and the Twenty-Sixth National Hydrodynamics Seminar, 2017, 423–428
[17]
FLOW-3D® User Manual.FLOW-3D User Manual. Flow Science, 2016
[18]
Yakhot V, Orszag S A. Renormalization group analysis of turbulence I. Basic theory. Plenum Press, 1986
[19]
Xiao Y, Zhang J B, Yao B, Guan Y. Assembly and simulation analysis of shear-sheet machine based on Pro/E. Procedia Engineering, 2011, 16(1): 535–539
CrossRef Google scholar
[20]
Bayon A, Toro J P, Bombardelli F A, Matos J, López-Jiménez P A. Influence of VOF technique, turbulence model and discretization scheme on the numerical simulation of the non-aerated, skimming flow in stepped spillways. Journal of Hydro-environment Research, 2017. doi: org/10.1016/j.jher.2017.10.002
[21]
Wan B L, He W X, Chen C Y, . Assembly process bill of material construction technology for dpacecraft assembly and integration based on proengineer. Joint International Information Technology, Mechanical and Electronic Engineering Conference,2017
[22]
Pan Z B, Lü H X, Zhang X F. Experiment on airfoil-shaped measuring flume in trapezoidal canal. Transactions of The Chinese Society of Agricultural Machinery, 2009, 40(12): 97–100 (in Chinese)
[23]
Najafi-Jilani A, Niri M Z, Naderi N. Simulating three-dimensional wave run-up over breakwaters covered by antifer units. International Journal of Naval Architecture and Ocean Engineering, 2014, 6(2): 297–306
CrossRef Google scholar
[24]
Duguay J M, Lacey R W J, Gaucher J. A case study of a pool and weir fishway modeled with open foam and FLOW-3D. Ecological Engineering, 2017, 103: 31–42
CrossRef Google scholar
[25]
Cui W, Song H F. CFD simulation of fresh self-compacting concrete flow and casting process. Concrete, 2017, 1(8): 111–115
[26]
Xiao Y Z, Wang W N, Hu X, Zhou Y. Experimental and numerical research on portable short-throat flume in the field. Flow Measurement and Instrumentation, 2016, 47: 54–61
CrossRef Google scholar
[27]
Tekade S A, Vasudeo A D, Ghare A D, Ingle R N. Measurement of flow in supercritical flow regime using cutthroat flumes. Sādhanā, 2016, 41(2):265–272.
[28]
Samani Z. Three simple flumes for flow measurement in open channels. Journal of Irrigation & Drainage Engineering, 2017, doi:org/10.1061/(ASCE)IR.1943-4774.0001168
[29]
Hu H, Huang J, Qian Z, Yu G. Hydraulic analysis of parabolic flume for flow measurement. Flow Measurement and Instrumentation, 2014, 37: 54–64
CrossRef Google scholar
[30]
Zhang L, Wu P T, Zhu D L, Zheng C. Flow regime and head loss in a drip emitter equipped with a labyrinth channel. Journal of Hydrodynamics, 2016, 28(4): 610–616
CrossRef Google scholar

Acknowledgements

The authors would like to acknowledge the financial support given by the Special Fund for Agro-scientific Research in the Public Interest of China (201503125) and the National Key Research and Development Program of China (2016YFC0400200).

Compliance with ethics guidelines

Danjie Ran, Wene Wang, and Xiaotao Hu declare that they have no conflicts of interest or financial conflicts to disclose.
This article does not contain any studies with human or animal subjects performed by any of the authors.

RIGHTS & PERMISSIONS

The Author(s) 2018. Published by Higher Education Press. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0)
AI Summary AI Mindmap
PDF(1213 KB)

Accesses

Citations

Detail

Sections
Recommended

/