Next Article in Journal
Temporal Dynamics of the Hydropower Water Reservoirs of the Tocantins–Araguaia Basin, Brazil, Based on Remote Sensing and Hydrometeorological Station Datasets
Next Article in Special Issue
Simulation Study of CH4 and N2O Emission Fluxes from Rice Fields in Northeast China under Different Straw-Returning and Irrigation Methods Based on the DNDC Model
Previous Article in Journal
Characteristics of Sedimentary Organic Matter in Tidal Estuaries: A Case Study from the Minjiang River Estuary
Previous Article in Special Issue
WINDS Model Demonstration with Field Data from a Furrow-Irrigated Cotton Experiment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of the Midseason Crop Coefficient for the Evaluation of the Water Demand of Young, Grafted Hazelnut Trees in High-Density Orchards

1
Department of Agricultural, Food and Environmental Sciences, University of Perugia, Borgo XX Giugno, 74, 06121 Perugia, Italy
2
Institute of Research on Terrestrial Ecosystems (IRET), National Research Council (CNR), Via Marconi 2, 05010 Porano, Italy
3
National Biodiversity Future Center, 90133 Palermo, Italy
*
Authors to whom correspondence should be addressed.
Water 2023, 15(9), 1683; https://doi.org/10.3390/w15091683
Submission received: 15 February 2023 / Revised: 16 April 2023 / Accepted: 25 April 2023 / Published: 26 April 2023
(This article belongs to the Special Issue Model-Based Irrigation Management)

Abstract

:
Knowledge of crop water requirements is important in supporting irrigation management. Evapotranspiration (ET) is commonly measured with a variety of instruments and field procedures, but it is also typically computed or modeled using the FAO56 or FAO66 methods. The adoption of this approach requires the assessment of the crop coefficients. Some data are available for own-rooted hazelnut trees, but no data have been reported for young and grafted hazelnut trees. There is a need to update nut–tree crop coefficients, especially considering modern cultivars and production systems, such as those with a high tree density per ha−1. In this paper, the FAO66 crop transpiration coefficient Kc,Tr and the FAO56 dual crop coefficients Kcb were assessed for the mid-growing season of a young grafted hazelnut orchard. The field data were acquired manually and using UAV. The coefficients were determined for three tree densities and for two growing seasons. The crop coefficients, obtained using the FAO66 method, agreed with the literature data referring to low densities, while the FAO56 method could allow us to better define the crop coefficients for high-density hazelnut orchards.

1. Introduction

Water scarcity, in relation to climate change, has become the main barrier to maintaining crop production, especially in Mediterranean areas [1,2,3,4]. Therefore, it is important to develop sensible irrigation strategies to reduce water use while maintaining adequate yields [5]. Many studies on the crop water requirements of vegetables have been produced since the publication of the FAO56 guidelines in 1998 [6] and FAO66 guidelines in 2012 [7] for the yield response to the water of fruit trees and vines. A recent review has shown that the determination of crop evapotranspiration (ETc), using the crop coefficient Kc—the reference evapotranspiration ET0 approach—is the most widely used method for irrigation water management [8]. Another recent study reviewed [9] the research on the FAO56 crop coefficients of fruit trees and vines performed over the past twenty years. It showed that many studies have been published for table and wine grapes, oranges, and olives. In addition, a few recent studies have been conducted on almonds (Prunus dulcis Mill. D. A. Webb) [9].
To date, there is a lack of studies on the crop coefficients of hazelnuts (Corylus avellana L.), pecans (Carya illinoinensis L.), pistachios (Pistacia vera L.), and walnuts (Juglans regia L.) and their application in irrigation management. For hazelnut trees, only two recent studies have been identified [10,11]. The first one [10] was conducted on young trees from the “Itarski Dugi” cultivar in a micro-sprinkler orchard (833 plants ha−1) in Serbia. In this work, the soil water balance (SWB) was used to measure the reference evapotranspiration (ET0). The second [11] was conducted on hazelnut trees from the “Tonda Giffoni” cultivar in a drip-irrigated orchard (333 plants ha−1) in the Maule region in Chile. In this work, the eddy covariance system was used to determine the Kc. Past works used lysimeters to determine the hazelnut coefficients Kc of the cultivar “Negret”, grown in a Spanish orchard with 500 plants ha−1, and of the cultivars “Barcellona” and “Ennis”, which were grown in a French orchard with 416 plants ha−1 [12,13]. However, in the FAO56 [6] tabulated single basal crop coefficients for use with the FAO Penman–Monteith ET0, Kc values are not reported for hazelnuts, but only for walnuts and pistachios. In addition, all the aforementioned studies were carried out on own-rooted bush-trained hazelnut trees, produced by layering, with a low tree-planting density. Based on our knowledge, until now, no studies have been conducted on hazelnut trees grafted on not-suckering rootstock, such as “Corylus colurna” L., with a single trunk and a vase training system. A recent study showed that grafted plants may access deeper soil water layers in comparison with own-rooted ones and have the potential to more effectively tolerate thermal variations and drought conditions during the vegetative and productive seasons [14]. In light of the above factors, as already stressed by [9], there is a need to update nut trees’ crop coefficients, in particular for hazelnuts, while considering modern cultivars and production systems, especially those with increasing tree density.
The aim of this work was to assess and compare the FAO66 crop transpiration coefficient Kc,Tr, and the FAO56 dual crop coefficients, Kcb, for the mid-growing season of young, grafted hazelnuts from a new, recently patented cultivar (n. EU52157). We make reference to three tree densities, ranging from one that is more common in traditional growing areas (lowest density) to those with many more trees, which are common in modern hazelnut orchards, highly mechanized, and always irrigated [15].

2. Materials and Methods

2.1. Study Site Description and Tree Sampling

The research was carried out in an experimental orchard of the Department of Agricultural, Food, and Environmental Sciences of the University of Perugia, located in Central Italy (42°58′22.82″ N, 12°24′13.02″ E), 161 m a.s.l., during the growing season of 2021 and 2022. The orchard has an area of about 2300 m2, and it is composed of 6 rows 4 m away from each other, each containing 43 trees divided into 3 plantation densities as follows:
  • 625 trees ha−1, spaced 4 m between rows and 4 m on the row, used as a control treatment (the common density used by farmers);
  • 1250 trees ha−1, spaced 4 m × 2 m;
  • 2500 trees ha−1, spaced 4 m × 1 m.
The tree densities were chosen according to the past literature [12,16]. The low tree density is representative of traditional hazelnut orchards, while the medium and high densities are currently used in the new cultivation areas. Specifically, on the same row, the first 25 plants are spaced 1 m from each other, the next 12 are spaced 2 m apart, and the last 6 are placed 4 m apart. Each tree density is replicated three times, once per row.
All trees, grafted on no-suckering rootstock, were trained as single trunks with four main branches each.
The orchard is made up of trees belonging to two of the main Italian hazelnut cultivars, chosen for their higher productivity and early production: Tonda di Giffoni in the first three rows and Tonda Francescana® in the next three rows.
In this paper, the trees studied belonged only to the hazelnut cultivar Tonda Francescana®. The meteorological data were monitored using a Spectrum (Thayer Court, Aurora) WatchDog 2000 Series Weather Station [17] located within the experimental site and were recorded at a time resolution of 5 min.
The DJI (Shenzhen, China) Phantom 4 (P4) multispectral unmanned aerial vehicle UAV was used to assess the tree canopy and ground cover characteristics according to [18,19]. Moreover, the FAO66 procedure described in Section 2.3, calibrated for olives, was used to parametrize the Kc,Tr coefficients for the hazelnut orchard.

2.2. Manual Measurements

The following data were acquired to fill the equations reported in Section 2.3.
In other tree crops, corn (Zea mays L.), cotton (Gossypium arboretum L.), peanuts [20], soybean (Glycine max (L.) Merr.), and wheat (Triticum aestivum L.) LAI was remotely estimated using photogrammetry and UAVs; studies have been performed to extract the LAI [21]. Very few studies have been carried out on tree crops, mainly on vineyards using higher spatial resolution remote sensing. This species is trained with the vertical shoot position trellis system or pergola system with horizontal shoot position [21,22], while hazelnut plant has a unique plant architecture with bush growth habit and dense foliage, three or four main branches, and canopy shape similar to the cylinder; all of this makes it more difficult to implement LAI models using UAV. For these reasons, in the present study, data determined using UAV referred only to canopy characteristics such as tree height.
LAIplant (the leaf area index, that is, the ratio of the plant leaf area to the canopy cross-sectional area), LAIorchard (the ratio of the plant leaf area to the tree distance area), and LAD or DAF (the leaf area density) were measured and calculated following the procedure reported by [23] for 9 trees per growing season.
The LAI values were manually measured using an iron bar (2.5 m long, 0.1 m wide, 0.1 m thick) with holes, for a total of 25 measurements, positioned 10 cm apart when a metal rod is inserted at a perpendicular angle; the iron bar is equipped with a level to keep the rods perpendicular to the ground. The iron bar is moved from the trunk to the inter-row by an operator, who inserts the metal rods and counts the leaves at the intersection points on a grid set up in the canopy. Measurements with this instrument were performed in two canopy directions, to the east side and the west side. The mean leaf area per tree density plantation was calculated with SigmaScan Pro5 (Richmond, CA, USA) program [24] by measuring the area of 30 leaves (5 for each tree) collected from the whole canopy.
Moreover, ground cover was obtained as the ratio between the average canopy cross-sectional area and the ground area for each plant, equal to 16 m2, 8 m2, and 4 m2, respectively, from the highest to the lowest tree density plantation. The average canopy cross-sectional areas were calculated as reported by [19].
The radiation intercepted from the tree Qd was measured manually using a luxmeter; this is a kind of ceptometer but is larger in size, as it consists of metal bars: a fixed one that is planted on the ground and horizontal to it, which is 150 cm long, and another piece of furniture perpendicular to the first, which is 100 cm long. The latter has a support on top to hold a special bar equipped with a photosensitive sensor, which generally consists of a transducer (a photovoltaic cell). When reacting to light intensity, it generates an electric current that is detected by a specific galvanometer equipped with a screen and connected with cables calibrated in lux. The movable bar can therefore slide in different positions along the horizontal bar, from 0 cm (ground) to 50 cm or 100 cm, and it can also be increased in height by means of a screw, from 0 cm (ground) to 100 cm. In this way, the intensity of light is detected at different heights of the plant and is compared with the light measured in full sunlight, revealing the light interception by the tree canopy. This parameter is important because there is a direct relationship between cumulative ET and PAR (photosynthetically active radiation) light interception by the tree canopy; a higher PAR light interception value corresponds to higher ET [25]. To obtain an estimation of the light interception by the whole tree canopy, the canopy light interception is expressed as the mean value of all data recorded at different distances from the trunk and tree heights. A completely randomized design (CRD) is the model used so that each experimental unit has the same chance of receiving any one treatment. For the CRD, any difference among experimental units receiving the same treatment is considered an experimental error.

2.3. Evaluation of the Hazelnut Orchard’s Crop Coefficients Using the FAO66 Method

As assessed by FAO66 [7], the canopy size and stomatal conductance are the two main parameters that determine tree transpiration Tr. Therefore, the models needed to calculate Tr in different orchards, as a function of the size and conductance of the tree canopy, must be specific for each tree species. Orchard ETc is considered the net irrigation requirement and, in well-managed irrigation systems, is the major component of the irrigation requirements. Orchard ETc may be calculated using crop coefficients and ET0, or by estimating its individual components, namely, transpiration and evaporation from the soil.
Orchard or vineyard ETc may be calculated as the sum of four components [7]:
E T c = T r + T r C C + E w z + E d z
where Tr is tree transpiration, Trcc is cover crop (or actively growing weed) transpiration, Ewz is the surface evaporation from the soil wetted by the emitters, and Edz is the surface evaporation from the rest of the soil surface outside the emitter wetting pattern.

Evaluation of Tree Transpiration Tr

Steduto et al. [7] proposed a method to derive the tree transpiration Tr for an orchard, as follows:
T r = E T 0 · K c , T r
where Tr is the tree transpiration (mm); ET0 is the reference evapotranspiration (mm); and Kc,Tr is the transpiration coefficient, as follows:
K c , T r = ( Q d · F 1 ) · F 2
F1 depends on tree density, i.e., F1 = 0.66 for tree densities > 250 trees/ha. F2 is the monthly tabulated coefficient; in this case, F2 = 1.25 for June and July, and F2 = 1.20 for August.
The parameter Qd is the radiation intercepted by the tree (fraction), calculated as follows:
Q d = 1 e K e x t V u
with
K e x t = 0.52 + 0.00079 · d p 0.76 · e 1.25 · D A F
Kext is the radiation extinction coefficient; dp is the tree density (trees/ha); and DAF is the leaf area density:
D A F = 2 0.53 · ( V u 0.5 ) ,   ( note : D A F   must   be < 2 )
with
V u = V 0 · ( d p 10000 )
where Vu is the canopy volume per unit ground surface (m3/m2), and V0 is the canopy volume (m3/tree):
V 0 = 1 6 · π · D 2 · H
where D is the canopy’s average diameter (m), H is the canopy height (m), and e is the exponent equal to 2.718.
The coefficients of Equations (5) and (6) were calibrated and tested for olives in a Mediterranean environment [26,27].
For a preliminary evaluation of the method, Equations (5) and (6) were applied to the hazelnut orchard.
In Equation (7), the canopy volume V0 was substituted with the volume Vc, evaluated using the procedure described by [18,19]. Briefly, volume Vc was calculated considering the shape of the canopy obtained from the 3D point cloud reconstruction using a UAV survey. From the raster file of the canopy, the volume was obtained in a GIS environment by evaluating the volume between the DSM and a horizontal plane passing through the lowest point of the same. Thus, the evaluation of the canopy volume per unit ground surface V u (m3/m2) was made by substituting Equation (7) with the following formula:
V u = V c · d p 10000
The tree density, dp, was set to be equal to 625 trees/ha, 1250 trees/ha, and 2500 trees/ha.
The leaf area density, DAFFAO, was evaluated using Equation (6) with V u :
D A F F A O = 2 0.53 · ( V u 0.5 )
The DAFFAO was used for Equation (5), and by applying Equations (3) and (4), the K c , T r was obtained.
As described above, the coefficients of Equations (5) and (6) were derived for the olive orchards. Using manual measurements of the leaf area density for each tree density and for the whole growing season (reported later), the coefficients of Equation (6) were derived. Therefore, manual measurements of the radiation intercepted by the tree, Qd,m, were available. Thus, these were used in two ways:
  • To calculate the K c , T r using Equation (3);
  • To derive the values of the product K e x t · V u from Equation (4). Substituting the V u values, obtained from Equation (9), the radiation extinction coefficients K e x t were evaluated, and the coefficients of Equation (5) were determined for the hazelnut orchard.
Student’s t-test was conducted to evaluate the significant differences between results obtained for the different plant densities.

2.4. Estimation of the Hazelnut Orchard’s Crop Coefficients from the Fraction of Ground Cover and Height (Allen and Pereira [16])

The FAO56 procedure for the estimation of the dual crop coefficient Kcb shows that it is largely affected by the fraction of the soil surface covered by vegetation. Recently, Allen and Pereira [6,15] suggested using a density coefficient Kd as function of tree heights. Focusing on orchards, Allen et al. [28] and Allen and Pereira [15] combined the influence of crop density with standardized values for Kcb full, Kcmin, and Kcb cover to compute a Kcb value for orchards for both bare soil and active ground cover. This modification provided flexibility in adjusting the value for Kc and Kcb according to the fraction of ground shaded by canopy, which is highly variable between orchards [29].
The single Kc approach is not the most appropriate procedure for estimating crop ET for conditions with variable wetting of the soil surface, particularly when the ground cover is minimal. Changes in soil surface wetness are influenced by the frequency and duration of rainfall and irrigation events, as well as by the type of irrigation system or strategy employed, and these cannot be entirely accounted for using canopy cover estimates. This is particularly noticeable under the conditions of partial canopy cover, when soil evaporation (Es) is expected to be a considerable component of ET because increased energy reaches the soil’s surface. Due to this limitation, a dual Kc approach was developed to estimate crop ETc, which splits the Kc into a T (Kcb) component and an Es (Ke) component [6]:
E T c = ( K c b + K e ) · E T 0
The basal crop coefficient Kcb is defined as the ratio of the crop evapotranspiration over the reference evapotranspiration (ETc/ET0) when the soil surface is dry, but transpiration is occurring at a potential rate, i.e., water is not limiting transpiration. Therefore, ‘Kcb · ET0’ primarily represents the transpiration component of ETc.
The basal Kcb is correlated with the amount of vegetation because it mostly represents transpiration. It can be expressed in terms of a density coefficient Kd, where [15]
K c b = K c m i n + K d · ( K c b f u l l K c m i n )
where Kcb is the approximation for conditions represented by the density coefficient Kd, Kcb full is the estimated basal Kc during peak plant growth under the condition of nearly full ground cover (or LAI > 3), and Kc min is the minimum basal Kc for bare soil (Kc min = 0.05 for orchards [30]). The coefficient Kcb full represents a general upper limit for tall vegetation with full ground cover and LAI > 3 under a full water supply; Kcb full, for use with ET0, can be approximated as a function of the mean plant height and adjusted for climate following [6]:
K c b f u l l = F r · ( min 1.0 + k h · h , 1.20 + [ 0.04 · u 2 2 0.004 · R H m i n 45 ] · h 3 0.3 )
where h is the mean maximum plant height (m), u2 is the mean value for wind speed at the height of 2 m during the midseason (m/s), RHmin is the mean value for the minimum daily relative humidity during the midseason (%), and Fr is an adjustment factor relative to crop stomatal control, described below. As assessed by Pereira et al. [31], the effects of crop height are considered through the sum ( 1.0 + k h · h ) , with kh = 0.1 for tree crops.
For trees and vines, the parameter Fr is high when crops exhibit significant vegetative vigor and decreases with pruning and training, as well as with limited water supplies [29]. Adopting the definition in [6] for Fr, considering the variability of leaf resistance between annual and perennial crops, it is assumed that
F r = + γ ( 1 + 0.34 u 2 ) + γ ( 1 + 0.34 u 2 r l r t y p )
where rl and rtyp are, respectively, the mean leaf resistance and the typical leaf resistance [s m−1] for the vegetation in question, Δ is the slope of the saturation vapor pressure vs. the air temperature curve [kPa °C−1], and γ is the psychrometric constant, kPa °C−1, both relative to the period when Kcb full is computed. The original Fr was developed to empirically consider the effects of stomatal adjustment on Kcb full and, therefore, on Kcb, in response to water stress [6], since stomatal closure causes a decrease in transpiration and thus decreases Kcb and Kc. The stomatal adjustment process is well known for annual crops, while, for fruit trees, as for hazelnut orchards, it was less well known or unknown. Thus, the effect of the stomatal adjustment processes on transpiration could be characterized by the empirical adjustment of Fr. Some default Fr values were indicated by [15] for tree and vine crops.
The density coefficient Kd can be estimated as a function of the measured leaf area index, LAI [m2 m−2], defined as the area of leaves per area of ground surface averaged over a large area:
K d = ( 1 e 0.7 · L A I )
As described by [13], where estimates of the fraction of the ground surface covered by vegetation, fc, are available, the density coefficient Kd can be estimated as follows:
K d = m i n ( 1 , M L f c e f f , f c e f f 1 1 + h )
where fc eff is the effective fraction of ground covered or shaded by vegetation [0.01–1] near solar noon, ML is a multiplier of fc eff that describes the effect of the canopy density on shading and on the maximum relative ET per fraction of ground shaded [1.0–2.0], and h is the mean vegetation height (m). Thus, using the Kd values obtained from Equation (15) and the measurement data of fc eff and h, the ML coefficient was estimated.

3. Results

3.1. Application of the FAO66 Method

Using the measurements of the canopy’s average diameter D, the canopy height H, and the canopy volume obtained from the UAV survey [18,19], V u , the equations from (3) to (6) were applied to the hazelnut orchard for each plant density (Figure S1). The mean value for the midseason (June, July, and August) of the coefficients K c , T r , evaluated using the coefficients tested for olive orchards [7] are reported in the first row of Table 1. The resulting coefficients were significantly different for all the densities; in fact, Student’s t-test showed that the means were significantly different (α = 0.05).
For a tree density of 625 trees/ha, the values of DAFFAO were always higher than 2, and Equation (6) was set for DAF values lower than 2; as such, an adjustment of the coefficients of Equation (6) seems to be necessary. Using manual measurements of the leaf area density DAFm for each tree density, as reported in Table 2, and for the whole mid-growing season, the coefficients of Equation (6) were derived:
D A F m = 1.79 + 2.83 · V u
therefore, as described above, the manual measurements of the radiation intercepted by the tree, Qd,m, were available (Figure S2). Therefore, Equation (5) could be modified for the evaluation of the radiation extinction coefficient, Kext, as follows:
K e x t = 0.385 + 0.0047 · d p + 4.67 · e 0.909 · D A F
where dp is the tree density (trees/ha), DAF is the leaf area density, and e is equal to 2.718.
From the application of Equations (3), (4), (17) and (18), the coefficients K c , T r were evaluated and are reported in the second row of Table 1.
The coefficients were significantly different only between the densities of 625 trees/ha and 2500 trees/ha and between the densities of 1250 trees/ha and 2500 trees/ha, since Student’s t-test showed that the mean was significantly different from 0 (α = 0.05).

3.2. Application of the Allen and Pereira Method [15]

The procedure detailed in [15] was used to derive the Kcb. Equation (14) was not applicable because the value of the mean leaf resistance (rl) is not available in the literature for hazelnut orchards. Due to these missing data, in this study, assuming a similarity in the water sensibility behaviors of hazelnut and walnut orchards [32], the Fr value of walnuts reported in Table 2 of [15] was used (Fr = 0.90). The values of the mean tree height measured using the UAV surveys [18,19] were used for the evaluation of Kcb full using Equation (13). The density coefficient Kd was calculated using the measurements of LAI conducted during the midseason and applied to Equation (15) (Figure S3). Thus, Equation (12) was used to estimate the values of Kcb, with Kc min = 0.05 [31]. The values obtained for the three tree densities (625 trees/ha, 1250 trees/ha, and 2500 trees/ha) are reported in Table 3.
From the values of the density coefficients Kd, calculated using Equation (15), Equation (16) was inverted to estimate the ML coefficient as described by [8]. In particular, the fc eff values were obtained from the following equation:
f c e f f = f c s i n ( β )
where β is the mean angle of the sun above the horizon during the period of maximum ET, evaluated using the procedure described by [6]; fc is the fraction of the surface covered by the vegetation, measured from the UAV surveys, described in Vinci et al. [18,19]. The coefficients were significantly different only between the densities 625 trees/ha and 2500 trees/ha and between the densities 1250 trees/ha and 2500 trees/ha, since Student’s t-test showed that the mean was significantly different from 0 (α = 0.05).
The values of ML obtained were ML = 1.17 for 625 trees/ha, ML = 1.09 for 1250 trees/ha, and ML = 1.01 for 2500 trees/ha.

4. Discussion

Pereira et al. [8] suggest that the selection of Fr and ML values is performed through a numerical search for the values of the parameters that make Kcb match the tabulated Kcb TAB values. However, as noted by Rallo et al. [9], studies for hazelnut trees are lacking, and Kcb tabulated values are not available for hazelnut orchards. The values of fc measured and used for the analysis were in accordance with the values of fc reported for walnut by [8], and for these values, the Kc mid was around 1 (Figure 1).
Figure 1 shows the crop coefficients reported for hazelnuts and for walnuts in the literature and those determined in this research, for several tree densities. The highest tree density of 2500 trees ha−1, used in this study, was not reported in Figure 1 since it has never been used by farmers, and no data are available.
The Kc,Tr determined in this research for 625 trees ha−1 density was similar to that reported by Silvestri et al. [12] for 500 trees ha−1 for own-rooted young hazelnuts. On the other hand, it was higher than that suggested by FAO66 [7]; this may be because the Kc used was calibrated and tested for olives. In fact, olive trees are an evergreen and resilient species that are well adapted to the Mediterranean environment [26,27,33]; meanwhile, hazelnut is sensitive to water stress and prefers cooler and more humid environmental conditions [34,35]. However, for the medium density (1250 trees ha−1), the Kc,Tr determined for hazelnut in this research was very similar to that of FAO66 [7]. This could be related to the species-specific physiological responses under changing light absorption conditions and nutrient availability. The values of the coefficients Kcb determined in this study by using the Allen and Pereira method [6] resulted in the highest values among those reported in the literature for hazelnuts (Figure 1). We suppose that, in applying this method, we overestimated the Kcb values for hazelnuts, with values very similar to those reported for walnuts (Figure 1).
Allen and Pereira [6] suggested a relationship between Kc and plant height in connection with the net radiation increment. However, we did not find any differences between the Kcb values of the two plant densities shown in Figure 1. This could be due to the specific training shape of hazelnut trees, which have four branches and no single axis. Therefore, the hazelnut tree’s height does not reflect the actual tree shape and canopy volume, especially in young trees. On the other hand, the Kc,Tr determined in this research for both crop densities was lower than the values reported in the literature for hazelnuts (Figure 1). The higher Kc values reported by Mačkić et al. [10] were probably due to the higher planting density and the specific orchard irrigation method. The higher Kc values reported by Mingeau and Rousseau [13] and by Ortega-Farias et al. [11] could be related to the specific orchard characteristics and management practices.
However, since the Kc,Tr determined in this research for grafted plants turned out to be lower (0.61–0.65), albeit with a higher planting density, than those reported in the literature for own-rooted plants (Figure 1), we cannot overlook the effect of grafting on hazelnut crop coefficients. In fact, as reported by Portarena et al. [14], the deeper root systems of grafted plants compared with own-rooted plants lead to differences in the depth of the root water uptake and in physiological functioning at the whole-plant level. It is likely that grafted hazelnut trees also have better radiation interception with a lower LAI (Figure S3) compared to those recorded for own-rooted trees of the same age by other authors [36]. All these factors reflect the stomatal control capacity of trees and the entire orchard’s water requirements, improving orchard sustainability.

5. Conclusions

For the first time for grafted hazelnut trees on non-suckering rootstock, Kcmid has been assessed for the evaluation of the water demand of young hazelnut trees in a high-density orchard.
The values of the density coefficient Kd did not differ among the three tested densities, a result that conflicted with the method of Allen and Pereira (2009), who suggested a relationship between Kd and the tree height. Moreover, the Kcb estimated in our grafted orchard was lower than the values reported in the literature for own-rooted hazelnuts. These lower values could be primarily related to the specific local environmental conditions, the cultivar, and the orchard management system, including the grafting practices employed.
Further studies are needed to determine the Kcb for grafted adult hazelnut trees. Against this background, further research will be conducted to evaluate the Kc ini and Kc end values of grafted hazelnut trees in order to schedule a well-planned and efficient irrigation practice based on the orchard’s water requirements.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/w15091683/s1. Figure S1: Canopy volume (Vc) in the three plant densities (means ± s.e.); Figure S2: Intercepted radiation by the tree, expressed as percentage of the external radiation, in the three plant densities (means ± s.e.); Figure S3: Leaf area index (LAI) in the three plant densities (means ± s.e.).

Author Contributions

Conceptualization, A.V. and D.F.; methodology, A.V.; validation, A.V. and D.F.; formal analysis, A.V. and D.F.; investigation, A.V., C.T. and D.F.; resources, A.V. and D.F.; data curation, A.V., C.T. and D.F.; writing—original draft preparation, A.V. and D.F.; writing—review and editing, A.V., D.F. and S.P.; project administration, D.F.; funding acquisition, A.V. and D.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “Ricerca di Base 2020, grant code VINRICBASE2020”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Noto, L.V.; Cipolla, G.; Pumo, D.; Francipane, A. Climate Change in the Mediterranean Basin (Part II): A Review of Challenges and Uncertainties in Climate Change Modeling and Impact Analyses. Water Resour. Manag. 2023. [Google Scholar] [CrossRef]
  2. Vergni, L.; Vinci, A.; Todisco, F. Effectiveness of the new standardized deficit distance index and other meteorological indices in the assessment of agricultural drought impacts in central Italy. J. Hydrol. 2021, 603, 126986. [Google Scholar] [CrossRef]
  3. Vergni, L.; Vinci, A.; Todisco, F.; Santaga, F.S.; Vizzari, M. Comparing Sentinel-1, Sentinel-2, and Landsat-8 data in the early recognition of irrigated areas in central Italy. J. Agric. Eng. 2021, 52, 1265. [Google Scholar] [CrossRef]
  4. Vinci, A.; Vergni, L.; Todisco, F.; Mannocchi, F. Analysis of rainfed alfalfa evapotranspiration measured by an eddy covariance system. J. Agric. Eng. 2013, 44, 779–784. [Google Scholar] [CrossRef]
  5. Levidow, L.; Zaccaria, D.; Maia, R.; Vivas, R.; Todorovic, M.; Scardigno, A. Improving water-efficient irrigation: Prospects and difficulties of innovative practices. Agric. Water Manag. 2014, 146, 84–94. [Google Scholar] [CrossRef]
  6. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop evapotranspiration. guidelines for computing crop water requirements. In FAO Irrigation and Drainage Paper; FAO: Rome, Italy, 1998; p. 300. Available online: https://www.fao.org/3/X0490E/x0490e00.htm#Contents (accessed on 10 January 2023).
  7. Steduto, P.; Hsiao, T.C.; Fereres, E.; Raes, D. Crop yield response to water. In FAO Irrigation and Drainage Paper; FAO: Rome, Italy, 2012; Available online: https://www.fao.org/3/i2800e/i2800e00.htm (accessed on 10 January 2023).
  8. Pereira, L.S.; Paredes, P.; Lopez-Urrea, R.; Hunsaker, D.J.; Mota, M.; Mohammadi Shad, Z. Standard single and basal crop coefficients for vegetable crops, an update of FAO56 crop water requirements approach. Agric. Water Manag. 2021, 241, 106197. [Google Scholar] [CrossRef]
  9. Rallo, G.; Paço, T.A.; Puig, A.; Paredes, P.; Massai, R.; Provenzano, G.; Pereira, L.S. Updated single and dual crop coefficients for tree and vine fruit crops. Agric. Water Manag. 2021, 250, 106645. [Google Scholar] [CrossRef]
  10. Mačkić, K.; Pejić, B.; Belić, M.; Janković, D.; Pavlović, L. Hazelnut (Corylus avellana L.) response to microsprinkler irrigation in climatic conditions of Vojvodina province. Res. J. Agric. Sci. 2016, 48, 1–7. [Google Scholar]
  11. Ortega-Farias, S.; Villalobos-Soublett, E.; Riveros-Burgos, C.; Zúñiga, M.; Ahumada-Orellana, L.E. Effect of irrigation cut-off strategies on yield, water productivity and gas exchange in a drip-irrigated hazelnut (Corylus avellana L. cv. Tonda di Giffoni) orchard under semiarid conditions. Agric. Water Manag. 2020, 238, 106173. [Google Scholar] [CrossRef]
  12. Silvestri, C.; Bacchetta, L.; Bellincontro, A.; Cristofori, V. Advances in cultivar choice, hazelnut orchard management, and nut storage to enhance product quality and safety: An overview. J. Sci. Food Agric. 2021, 101, 27–43. [Google Scholar] [CrossRef]
  13. Mingeau, M.; Rousseau, P. Water use of hazelnut trees as measured with lysimeters. Acta Hotic. 1992, 351, 315–322. [Google Scholar] [CrossRef]
  14. Portarena, S.; Gavrichkova, O.; Brugnoli, E.; Battistelli, A.; Proietti, S.; Moscatello, S.; Famiani, F.; Tombesi, S.; Zadra, C.; Farinelli, D. Carbon allocation strategies and water uptake in young grafted and own-rooted hazelnut (Corylus avellana L.) cultivars. Tree Physiol. 2022, 42, 939–957. [Google Scholar] [CrossRef] [PubMed]
  15. Allen, R.G.; Pereira, L.S. Estimating crop coefficients from fraction of ground cover and height. Irrig. Sci. 2009, 28, 17–34. [Google Scholar] [CrossRef]
  16. Tombesi, A. Effect of Light Penetration on High-Density Filbert Planting; Institute of Fruit Culture, University of Perugia: Perugia, Italy, 1977; Volume 32, p. 33. [Google Scholar]
  17. Spectrum Technologies Inc. WatchDog 2000 Series Mini Stations product Manual. Available online: https://www.specmeters.com/assets/1/22/2000_Series_Mini_Loggers_2021.pdf (accessed on 13 February 2023).
  18. Vinci, A.; Traini, C.; Farinelli, D.; Brigante, R. Assessment of the geometrical characteristics of hazelnut intensive orchard by an Unmanned Aerial Vehicle (UAV). In Proceedings of the 2022 IEEE Workshop on Metrology for Agriculture and Forestry (MetroAgriFor), Perugia, Italy, 3–5 November 2022; pp. 218–222. [Google Scholar] [CrossRef]
  19. Vinci, A.; Brigante, R.; Traini, C.; Farinelli, D. Geometrical Characterization of Hazelnut Trees in an Intensive Orchard by an Unmanned Aerial Vehicle (UAV) for Precision Agriculture Applications. Remote Sens. 2023, 15, 541. [Google Scholar] [CrossRef]
  20. Sarkar, S.; Cazenave, A.B.; Oakes, J.; McCall, D.; Thomason, W.; Abbott, L.; Balota, M. Aerial high-throughput phenotyping of peanut leaf area index and lateral growth. Sci. Rep. 2021, 11, 21661. [Google Scholar] [CrossRef] [PubMed]
  21. Mathews, A.J.; Jensen, J.L. Visualizing and quantifying vineyard canopy LAI using an unmanned aerial vehicle (UAV) collected high density structure from motion point cloud. Remote Sens. 2013, 5, 2164–2183. [Google Scholar] [CrossRef]
  22. Ilniyaz, O.; Du, Q.; Shen, H.; He, W.; Feng, L.; Azadi, H.; Kurban, A.; Chen, X. Leaf area index estimation of pergola-trained vineyards in arid regions using classical and deep learning methods based on UAV-based RGB images. Comput. Electron. Agric. 2023, 207, 107723. [Google Scholar] [CrossRef]
  23. Farinelli, D.; Boco, M.; Tombesi, A. Influence of canopy density on fruit growth and flower formation. Acta Hortic. 2005, 686, 247–252. [Google Scholar] [CrossRef]
  24. SigmaScan Pro5, Inpixon Richmond, California (USA). Available online: https://systat.co.uk/ (accessed on 14 February 2023).
  25. Souto, C.; Lagos, O.; Holzapfel, E.; Ruybal, C.; Bryla, D.R.; Vidal, G. Evaluating a Surface Energy Balance Model for Partially Wetted Surfaces: Drip and Micro-Sprinkler Systems in Hazelnut Orchards (Corylus avellana L.). Water 2022, 14, 4011. [Google Scholar] [CrossRef]
  26. Testi, L.; Goldhamer, D.A.; Iniesta, F.; Salinas, M. Crop Water Stress Index Is a Sensitive Water Stress Indicator in Pistachio Trees. Irrig. Sci. 2008, 26, 395–405. [Google Scholar] [CrossRef]
  27. Orgaz, F.; Testi, L.; Villalobos, F.J.; Fereres, E. Water requirements of olive orchards II: Determination of crop coefficients for irrigation scheduling. Irrig. Sci. 2006, 24, 77–84. [Google Scholar] [CrossRef]
  28. Allen, R.G.; Pereira, L.S.; Smith, M.; Raes, D.; Wright, J.L. FAO-56 dual crop coefficient method for estimating evaporation from soil and application extensions. J. Irrig. Drain. Eng. 2005, 131, 2–13. [Google Scholar] [CrossRef]
  29. Pereira, L.S.; Allen, R.G.; Smith, M.; Raes, D. Crop evapotranspiration estimation with FAO56: Past and future. Agric. Water Manag. 2015, 147, 4–20. [Google Scholar] [CrossRef]
  30. Pereira, L.S.; Paredes, P.; Melton, F.; Johnson, L.; Wang, T.; López-Urrea, R.; Cancela, J.J.; Allen, R.G. Prediction of crop coefficients from fraction of ground cover and height. Background and validation using ground and remote sensing data. Agric. Water Manag. 2021, 241, 106197. [Google Scholar] [CrossRef]
  31. Pereira, L.S.; Paredes, P.; Hunsaker, D.J.; Lopez-Urrea, R.; Mohammadi Shad, Z. Standard single and basal crop coefficients for field crops. Updates and advances to the FAO56 crop water requirements method. Agric. Water Manag. 2021, 243, 106466. [Google Scholar] [CrossRef]
  32. Calvo, F.E.; Trentacoste, E.R.; Silvente, S.T. Vegetative growth, yield, and crop water productivity response to different irrigation regimes in high density walnut orchards (Juglans regia L.) in a semi-arid environment in Argentina. Agric. Water Manag. 2022, 274, 107969. [Google Scholar] [CrossRef]
  33. Portarena, S.; Leonardi, L.; Scartazza, A.; Lauteri, M.; Baldacchini, C.; Farinelli, D.; Famiani, F.; Ciolfi, M.; Brugnoli, E. Combining analysis of fatty acid composition and δ13C in extra-virgin olive oils as affected by harvest period and cultivar: Possible use in traceability studies. Food Control 2019, 105, 151–158. [Google Scholar] [CrossRef]
  34. Cincera, I.; Frioni, T.; Ughini, V.; Poni, S.; Farinelli, D.; Tombesi, S. Intra-specific variability of stomatal sensitivity to vapour pressure deficit in Corylus avellana L.: A candidate factor influencing different adaptability to different climates? J. Plant Physiol. 2019, 232, 241–247. [Google Scholar] [CrossRef]
  35. Di Lena, B.; Curci, G.; Vergni, L.; Farinelli, D. Climatic Suitability of Different Areas in Abruzzo, Central Italy, for the Cultivation of Hazelnut. Horticulturae 2022, 8, 580. [Google Scholar] [CrossRef]
  36. Altieri, G.; Maffia, A.; Pastore, V.; Amato, M.; Celano, G. Use of high-resolution multispectral UAVs to calculate projected ground area in Corylus avellana L. tree orchard. Sensors 2022, 22, 7103. [Google Scholar] [CrossRef]
Figure 1. Crop coefficients reported for hazelnut and for walnut in the literature [7,9,10,11,12,13] and those determined in this research, for several tree densities.
Figure 1. Crop coefficients reported for hazelnut and for walnut in the literature [7,9,10,11,12,13] and those determined in this research, for several tree densities.
Water 15 01683 g001
Table 1. Values of the K c , T r coefficients derived from the FAO66 method [7] and using the coefficients derived from the hazelnut orchard.
Table 1. Values of the K c , T r coefficients derived from the FAO66 method [7] and using the coefficients derived from the hazelnut orchard.
625 Trees/ha1250 Trees/ha2500 Trees/ha
Equations (3)–(8)0.41 c0.64 a0.61 b
Equations (3), (4), (17) and (18)0.61 b0.65 ab0.74 a
Note: In each row, means followed by different letters are significantly different, with p < 0.05.
Table 2. Values of the manual measurements of leaf area density DAFm, the intercepted radiation Qd, and the ground fraction cover fc.
Table 2. Values of the manual measurements of leaf area density DAFm, the intercepted radiation Qd, and the ground fraction cover fc.
DAFmQdfc
625 trees/ha1.560.5350.74
2.540.635
2.780.837
2.660.835
2.490.856
2.380.799
1.600.882
2.470.856
5.630.768
3.990.644
4.560.535
1250 trees/ha3.130.8630.795
1.470.828
1.480.714
1.960.856
2.350.737
1.800.844
2.400.934
2.410.806
5.890.532
3.260.889
7.010.748
2500 trees/ha4.500.8290.895
6.410.858
2.580.895
4.480.902
5.060.885
3.510.912
2.810.920
4.260.943
3.850.936
4.090.935
4.530.897
3.620.824
Table 3. Values of the coefficients K c b derived from the application of the Allen and Pereira [13] method ( K c b A&P).
Table 3. Values of the coefficients K c b derived from the application of the Allen and Pereira [13] method ( K c b A&P).
625 trees/ha1250 trees/ha2500 trees/ha
Equation (13)0.98 b0.98 b1.04 a
Note: Means followed by different letters are significantly different for p < 0.05.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vinci, A.; Traini, C.; Portarena, S.; Farinelli, D. Assessment of the Midseason Crop Coefficient for the Evaluation of the Water Demand of Young, Grafted Hazelnut Trees in High-Density Orchards. Water 2023, 15, 1683. https://doi.org/10.3390/w15091683

AMA Style

Vinci A, Traini C, Portarena S, Farinelli D. Assessment of the Midseason Crop Coefficient for the Evaluation of the Water Demand of Young, Grafted Hazelnut Trees in High-Density Orchards. Water. 2023; 15(9):1683. https://doi.org/10.3390/w15091683

Chicago/Turabian Style

Vinci, Alessandra, Chiara Traini, Silvia Portarena, and Daniela Farinelli. 2023. "Assessment of the Midseason Crop Coefficient for the Evaluation of the Water Demand of Young, Grafted Hazelnut Trees in High-Density Orchards" Water 15, no. 9: 1683. https://doi.org/10.3390/w15091683

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