Estimating Evapotranspiration of Rainfed Winegrapes Combining Remote Sensing and the SIMDualKc Soil Water Balance Model

1. Introduction

Estimating the evapotranspiration of orchards and vineyards using the soil water balance approach is often difficult due to the soil extent explored by the root system, especially under rainfed conditions [1,2,3]. Under such conditions, the root system varies between 0.5 and 6 m, although it is suggested that a small fraction of roots can grow to greater depths [4,5,6,7], at which soil water content measurements are not normally surveyed. Additionally, the root system of these crops tends to explore the soil further into the inter-row spacing. Grapevine root systems were studied across a wide range of climates and soil textures and showed that most active roots (about 80%) are located in the first 1.0 m of soil [7,8]. Vineyards have long been rainfed [6,9,10] as the species is inherently adapted to tolerate dry periods, typical of many wine-growing regions [6]. However, nowadays irrigation is commonly adopted because it allows improving the wine quantity and quality [11,12,13,14,15,16,17]. Vineyards have mechanisms to control transpiration, the most relevant of which is the pronounced stomatal control [5,6,11,18,19]. These plants also have a strategy to further explore the water reservoir in the soil by using a deep root system [5,6,13,20,21,22]. Both mechanisms allow for efficient use of available water [5,6,11,19,20,21,22], particularly in drier periods. Under rainfed conditions, the root system can be expected to develop further than in a comparable irrigated vineyard. There are few studies reporting a maximum rooting depths of vines greater than 6 m [23], however, it has been suggested that a small fraction of roots can grow to a depth of 20 m in the absence of impermeable barriers [4,5,6,7]. The root system of the vine depends on the characteristics of the soils, the planting density, and the crop management used. Soil characteristics that have the greatest impact on root distribution include the soil bulk density, resistance to soil penetration, and texture [19].

Soil water content is difficult to quantify due to the heterogeneity of the water distribution in the soil [3]. However, it is possible to assess the vineyard water status by analyzing the relationship between the actual crop evapotranspiration and its partitioning and the predawn leaf water potential (ψb) of grapevines. For instance, [24] reported a threshold of -0.3 MPa as limit for irrigating vineyards cv. Moscatel Galego-Branco. For water stress conditions, [25] also reported ψb ranging from -0.1 MPa to -0.36 MPa as indicative of weak to moderate water stress for rainfed vineyards. In addition, [18] quantified ψb varying of -0.3 MPa to -0.47 MPa for cv. Semillon and from -0.30 MPa to -1.07 MPa for cv. Muscat, both vineyards in rainfed conditions. [26] defined the water stress limits as follows: no water deficit (0 MPa ≥ ψb ≥ -0.2 MPa), mild to moderate water deficit (-0.2 MPa ≥ ψb ≥ -0.4 MPa), moderate to severe water deficit (-0.4 MPa ≥ ψb ≥ -0.6 MPa), severe to high water deficit (-0.6 MPa ≥ ψb ≥ -0.8 MPa), and high-water deficit (< -0.8 MPa).

A common approach to estimate crop water requirements is the FAO two-step approach which combines the climate conditions (reference crop evapotranspiration, ET_o) with the crop characteristics (crop coefficient, K_c) also named K_c-ET_o approach. This approach was proposed in FAO56 [27] and has been widely adopted [28]. The FAO56 K_c-ET_o approach can be applied with a single or a dual crop coefficient. In the first case, soil evaporation and crop transpiration are combined into a single K_c value for each crop stage; for the latter, daily plant transpiration is based on the basal crop coefficient (K_cb), while daily soil evaporation is estimated using an evaporation coefficient (K_e). Thus, ET_c is divided into crop transpiration (T_c = K_cb ET_o) and soil evaporation (E_s = K_e ET_o).

The standard tabulated values of K_c and K_cb allow the assessment of ETc under potential and well-watered conditions [27]. The standard tabulated K_c and K_cb values for trees and vines were recently reviewed and updated by [29]. However, under natural field conditions, the crop is often subject to biotic and abiotic stress due to water deficits caused by inadequate irrigation, improper management practices, soil quality and salinity, unsuitable crop varieties or rainfed conditions. Therefore, a water stress coefficient, K_s, in the range between 0 and 1.0 is introduced as a multiplicative factor to estimate actual values of K_c or K_cb, i.e., K_{c act} or K_{cb act} [27,28,30,31]. The actual crop evapotranspiration (ET_{c act}) is generally smaller than the potential value (ET_c) and can be defined as:

ET_{c act =} (K_s K_cb + K_e) ET_o = K_{c act} ET_o

(1)

The actual crop transpiration may be measured using the eddy covariance technique e.g. [32,33,34], the Bowen ratio energy balance method, soil water balance or a lysimeter as reviewed by [35]. The partitioning of evapotranspiration may be performed combining soil evaporation measured in micro-lysimeters or mini-lysimeters [36] with actual crop transpiration estimations with sap flow methods [37,38,39]. In addition, models properly calibrated allows assessing crop water requirements, to auxiliar adequately in the management of irrigation, and analysing the impact of crops and management, including under water deficit conditions. These models employ different functions and frequently do not use the FAO K_c-ET_o approach, e.g., the transient model Hydrus-2D [40]. Other examples have recently been reviewed [35]. Models applying the dual-K_c approach, include the SIMDualKc [41,42], which has been shown to be suitable for olive orchards [31,34,38,43] and other woody perennial crops, such as vineyards [44,45], peach trees [46] and almond, citrus, and pomegranate [43,47].

Alternatively to the models’ usage, [48] suggested an approach to obtain the basal crop and crop coefficients (K_cb, K_c) using a density coefficient (K_d), which calculated is estimated with information of fraction of the ground covered by the plants’ canopy (f_c) and crop height (h) [31,48]. In this approach, the estimation of K_cb considers two other parameters, the multiplier on f_c, which describes the influence of canopy density (M_L) and the resistance correction factor (F_r). M_L characterizes the transparency of the canopy to solar radiation while F_r is an empirical downward adjustment when vegetation exhibits tight stomatal control to transpiration. In the case of f_c, the value may be measured using ground or remote sensing data [35,49,50]. Application examples are provided by [38,50] and [35]. A full review of the parameter values for M_L and F_r is provided by [28].

In the last decades, remote sensing (RS) data have been used to estimate crop evapotranspiration considering two main approaches: (i) satellite-based surface energy balance models (SEB) [32,51,52,53,54] and (ii), spectral vegetation indices for estimation of basal crop coefficient (K_cb) based on FAO56 method, K_cb-VI [49,50,55,56,57]. The K_cb-VI approach requires information for a smaller number of variables than the SEB models, being based on elementary principles. Nevertheless, in contrast to the SEB models, the K_cb-VI approach ignores the effect of stomata closure related to the occurrence of water stress and assumes that this effect does not significantly affect the reduction in evapotranspiration compared to the effect of crop size [57].

The estimation of K_cb based on spectral vegetation indices has already been performed with different approximations. An overview on this topic is given by [49]. Two of the world’s most widely used spectral indices for this approximation are the Normalized Difference Vegetation Index, NDVI [58] and the Soil Adjusted Vegetation Index, SAVI [59]. The formulation of NDVI and SAVI combines red and near infrared (NIR) reflectance to provide an indirect measure of red-light absorption by chlorophyll (a and b) and reflectance of NIR by the mesophyll structure in leaves [49,60,61]. Compared to NDVI, SAVI is less sensitive to the soil backscatter effect, which is an advantage in crops with discontinuous ground cover, such as vineyards [59,62]. These spectral vegetation indices have also been used to determine the fraction of soil covered by the crop, f_c.

Significant gaps exist in the determination of the soil water available for the evapotranspiration process due to the uncertainty in the volume of soil explored by the root system of vines. Consequently, this uncertainty also occurs when modelling the soil water balance (SWB), particularly in the case of one-dimensional models. To improve the soil water balance and actual crop ET estimates in a rainfed managed vineyard, the current study uses a combined approach of ground and remote sensing data ingested in the SWB model SIMDualKc. To support the everyday decision making, the A&P approach was also tested using the estimates of SIMDualKc model. Thus, a comparison of the K_{cb SIMDualKc} with the K_cb derived from the observed values of f_c and h (K_{cb A&P}) [48] was included. The A&P approach is considered a useful approach to improve irrigation management and parameterize water balance models. The novelty of the current study lies in the precise determination of the soil water balance in a deep soil profile that can be explored by the roots, combined with a set of ground-truth data including predawn water potential measurements. The soil water content data allowed to perform an adequate parameterization, and calibration of a modelling tool. To the authors knowledge, this is the first time that the A&P approach has been used to estimate rainfed vineyard evapotranspiration.

2. Materials and Methods

In the present study, the ground data set was collected in a vineyard experimental site located in the Ribatejo region (Center Portugal). The remote sensing (RS) data characterizing the vineyard development was acquired by the Landsat 5 satellite (downloaded from the Earth Explorer portal; https://earthexplorer.usgs.gov/). The following subsections provide a detailed description of all data sets and modeling tool.

2.1. Study area and crop characterization

The experiments were developed in a 2.13 ha vineyard plot of Escola Superior Agrária (School of Agriculture), located in Santarém, Ribatejo region (39º 15' 38.15'' N, 8º 42' 37.7'' W) (Figure 1).

Data were collected in 1987, on a 20-year-old vineyard, cv. black Trincadeira, grafted on rootstock R110, planted at 3.0 m x 1.5 m, trained in a vertical shoot positioned trellis, pruned to the bead of 2 to 3 buds and a total load ranging between 6 to 10 buds, and with a high plant vigor. The vineyard was aligned, trellised, and oriented in the SSW-NNE direction. The timing of the phenological stages was recorded according to the Bagiolini scale as follows [19]: (i) cotton bud (bud swelling) in the first week of March, (ii) green tip (bud opening) occurred ten days later, (iii) onset flowering occurred between April 10 and 15, (iv) end of flowering between May 20 and 25, (v) pea-berry stage normally occurred in the first week of July and, (vi) complete ripening took place in the last week of September. These stages were then converted to those proposed by [27] as follows: initial stage from 05/03–09/04, development from 10/04–19/05, mid-season stage from 20/05–30/07, and late-season from 31/07–25/09 (Figure 2). Weed control in the inter-row was carried out with herbicides in April.

The vineyard predawn leaf water potential (ψb) and available soil water content relationship was analyzed as an indicator of the vineyard water stress level, and to better understand the soil-plant water relationships. Predawn leaf water potential measurements were carried out using a pressure chamber [63] on a total of nine dates, from June to September, on 7 to 10 leaves of different plants in each date. These measurements were considered to be representative of the vineyard. To minimize measurement errors, all ψb measurements were carried out close to the vines, with a time interval of around one minute between cutting the leaf and reaching the balance between the pressure in the chamber and the pressure in the xylem [19].

2.2. Climate Characterization

According to the Köppen-Geiger climate classification [64], the study area has a warm temperate climate (Csa), with dry and hot summers (T_air>22 ºC) and mild rainy winters, thus a typical Mediterranean climate. Average annual evapotranspiration ranges from 712 mm to 815 mm while average annual precipitation ranges from 674 mm to 844 mm, mostly concentrated in the winter months considering the years 1982-1987 [19].

Weather data were obtained from the Santarém weather station of (39º15'N, 8º42'W, 73 m a. s. l.), located 1 km away from the experimental site.

Daily data collected included maximum and minimum temperature (T_max and T_min, ºC), solar radiation (R_s, W m^-2), and average relative humidity (RH, %). The station is under reference site conditions [27], which allows the estimation of the reference evapotranspiration (PM-ET_o). Daily precipitation (mm) was also collected in the station. Based on [65], weather data were assessed in terms of quality and integrity. A brief characterization of the prevailing weather conditions in the study year is shown in Figure 3. Although the annual precipitation was 784 mm, it is clear that most of the precipitation was concentrated in the months of November to February, when the vineyards were dormant.

2.3. Soil characterization and soil water content monitoring

Soil properties were assessed by collecting disturbed soil samples at different depths. The vineyard soil in the studied area is an Anthrosol [66] with a sandy loam texture, down to a depth of 1.90 m, most of the sand being fine. Collection of undisturbed soil samples for determination of volumetric water content at field capacity was not possible due to the unstable soil surface layer and significant structural discontinuities. In addition, some soil layers exhibited very low saturated hydraulic conductivity (K_s<0.5 cm d⁻¹), which also made it difficult to determine soil moisture at field capacity. Thus, the volumetric soil water content at field capacity and at permanent wilting point were obtained using the pedotransfer functions determined for the Portuguese conditions [67] using the particle size and the soil bulk density.

A neutron probe (Humiterra, Laboratório Nacional de Engenharia e Tecnologia Industrial, Portugal) was used to monitor soil water content at each 0.10 m, up to 1.85 m, in two locations. Before use, the probe was calibrated following the procedures of [68] and [69]. As suggested by [68], special calibration of the probe was performed for the surface layers up to 0.15 m. Several transparent access tubes were placed surrounding the plants as shown in Figure 4. The measurements were carried out on seven dates every 16 days from 23 April to 8 September 1987.

2.4. Modelling approach

It was assessed the accuracy of a SWB model fed with ground truth data after proper calibration and a simplified approach that combines ground truth data with remote sensing data.

2.4.1. SIMDualKc modelling tool

The SIMDualKc modelling tool was developed with the aim of simplifying the implementation and computation of water use and irrigation requirements using the dual crop coefficient approach. The tool can be applied to a wide variety of situations including crops that do not fully cover the soil (e.g., orchards), intercropping systems, mulching impacts (plastic and organic), active ground cover, and soil and water salinity.

SIMDualKc is a one-dimension soil water balance modeling tool, that daily estimates the crop water requirements as follows [27,41], as:

$${\mathrm{D}}_{\mathrm{r},\mathrm{i}}={\mathrm{D}}_{\mathrm{r},\mathrm{i}-1}-{\left(\mathrm{P}-\mathrm{R}\mathrm{O}\right)}_{\mathrm{i}}-{\mathrm{I}}_{\mathrm{i}}-{\mathrm{C}\mathrm{R}}_{\mathrm{i}}+{\mathrm{E}\mathrm{T}}_{\mathrm{c}\mathrm{a}\mathrm{c}\mathrm{t},\mathrm{i}}+{\mathrm{D}\mathrm{P}}_{\mathrm{i}}$$

(2)

where D_r represents depletion in the root zone (mm), P represents the rainfall amount (mm), RO refers to the surface runoff (mm), I represents the net irrigation (mm), CR describes the capillary rise (mm) obtained of the groundwater table (mm), ET_{c act} is the actual crop evapotranspiration (mm), and DP is the is the root zone deep percolation (mm), all referring to days i or i-1.

The actual ET_c is partitioned into actual crop transpiration (T_{c act}) and soil evaporation (E_s) as:

$${\mathrm{T}}_{\mathrm{c}\mathrm{a}\mathrm{c}\mathrm{t}}=\left({\mathrm{K}}_{\mathrm{s}}{\mathrm{K}}_{\mathrm{c}\mathrm{b}}\right){\mathrm{E}\mathrm{T}}_{\mathrm{o}}={\mathrm{K}}_{\mathrm{c}\mathrm{b}\mathrm{a}\mathrm{c}\mathrm{t}}{\mathrm{E}\mathrm{T}}_{\mathrm{o}}$$

(3)

$${\mathrm{E}}_{\mathrm{s}}={\mathrm{K}}_{\mathrm{e}}{\mathrm{E}\mathrm{T}}_{\mathrm{o}}$$

(4)

The stress coefficient (K_s) is calculated daily as a linear function of the root zone depletion D_r (mm):

$${\mathrm{K}}_{\mathrm{s}}={\displaystyle \frac{\mathrm{T}\mathrm{A}\mathrm{W}-{\mathrm{D}}_{\mathrm{r}}}{\mathrm{T}\mathrm{A}\mathrm{W}-\mathrm{R}\mathrm{A}\mathrm{W}}}={\displaystyle \frac{\mathrm{T}\mathrm{A}\mathrm{W}-{\mathrm{D}}_{\mathrm{r}}}{\left(1-\mathrm{p}\right)\mathrm{T}\mathrm{A}\mathrm{W}}}\mathrm{for}{\mathrm{D}}_{\mathrm{r}}>\mathrm{R}\mathrm{A}\mathrm{W}$$

(5a)

$${\mathrm{K}}_{\mathrm{s}}=1\mathrm{for}{\mathrm{D}}_{\mathrm{r}}\le \mathrm{R}\mathrm{A}\mathrm{W}$$

(5b)

wherever TAW represents the total available soil water (mm) and RAW is the readily available soil water (mm), with RAW= p TAW, and p is the soil water depletion fraction for conditions of no-stress [27,48]. When the root zone depletion is higher than RAW, i.e., in situations of the water depleted fraction is larger than p, the stress coefficient is lower than one (K_s < 1.0).

The evaporation coefficient (K_e) will be the maximum immediately after rainfall or irrigation events occur and the wetted soil surface is being partially covered by plants canopy. This higher K_e value is dependent of K_c value be maximum (i.e., K_{c max}), which signifies the upper threshold of evaporation and transpiration from any cropped surface and reflects the natural restrictions imposed by the available energy on the transpiration coefficient K_cb on the same day, and thus on K_{c max} – K_cb [48,70]. Nevertheless, it changes day-to-day as a function of the total water available for evaporation in the soil surface layer; thus, the estimative of K_e is made by a dimensionless reduction factor, K_r, as:

$${\mathrm{K}}_{\mathrm{e}}={\mathrm{K}}_{\mathrm{r}}\left({\mathrm{K}}_{\mathrm{c}\mathrm{m}\mathrm{a}\mathrm{x}}-{\mathrm{K}}_{\mathrm{c}\mathrm{b}}\right)\le {\mathrm{f}}_{\mathrm{e}\mathrm{w}}{\mathrm{K}}_{\mathrm{c}\mathrm{m}\mathrm{a}\mathrm{x}}$$

(6)

being f_ew the wetted soil fraction exposed to direct solar radiation, which is calculate as f_ew= min (1 – f_c, f_w) according to the fractions of ground covered by plants’ canopy (f_c) and wetted by irrigation (f_w).

The value of K_r depends on the cumulative depth of water depleted from the topsoil and is daily calculated through the soil water balance considering the soil evaporation layer’s thickness (Z_e, m). Calculations are based on the assumption that the drying cycle happens in two stages, with the water in the topsoil and the energy available at the soil's surface limiting the first stage [36,48], thus with K_r reducing once the evaporated depleted water exceeds the readily evaporable water (REW). Thus, K_r is computed as:

$${\mathrm{K}}_{\mathrm{r}}=1\mathrm{for}{\mathrm{D}}_{\mathrm{e},\mathrm{i}-1}\le \mathrm{R}\mathrm{A}\mathrm{W}$$

(7a)

$${\mathrm{K}}_{\mathrm{r}}={\displaystyle \frac{\left(\mathrm{T}\mathrm{E}\mathrm{W}-{\mathrm{D}}_{\mathrm{e},\mathrm{i}-1}\right)}{\mathrm{T}\mathrm{E}\mathrm{W}-\mathrm{R}\mathrm{E}\mathrm{W}}}\mathrm{for}{\mathrm{D}}_{\mathrm{e},\mathrm{i}-1}>\mathrm{R}\mathrm{E}\mathrm{W}$$

(7b)

where D_e,i-1 is the amount of water evaporated from the evaporable soil layer at the end of day i-1.

The adjusted K_cb calculation already includes a water stress coefficient (K_s), that represents situations where water is a limiting factor. The soil evaporation coefficient (K_e) in turn is calculated by separating this coefficient into two parts: the K_ei relative to the fraction of soil wetted by irrigation and precipitation (f_ewi); and the K_ep, which is based on the fraction of soil wetted only by precipitation (f_ewp).

Mandatory SIMDualKc input data include: i) weather data required for calculation of the PM-ETo or computed elsewhere, precipitation and, in case of adjustment of K_c/K_cb to weather conditions, u₂ and RH_min; ii) crop characteristics data, which include the dates of the diverse crop growth stages, the crop height, rooting depth, the depletion fraction for no stress, the basal crop coefficients and the fraction of ground cover or the LAI for each defined stage; iii) soil characteristics, such as the soil water holding characteristics (soil water at field capacity and wilting point), or the TAW, or, as in the case of the current study, the use of pedotransfer functions (PTF) used to derive the TAW for diverse soil layers; iv) the surface layer characteristics: depth (m), total evaporable water (TEW, mm) and readily evaporable water (REW, mm); v) irrigation schedules that include different options, such as rainfed.

Other model functions were used to improve the SWB estimates. Therefore, other input data used with SIMDualKc in the current study included: i) the deep percolation or drainage fluxes (mm d⁻¹) is calculated with a decay function over time describing the soil water storage (W) close to saturation in the period following events of intense rainfall or irrigation. The SIMDualKc model uses a decay function [71] for the estimation of the deep percolation as follows:

$$\mathrm{W}=a_D{t}^{b_D}$$

(8)

where parameter a_D is the soil water storage value comprised between the field capacity and the saturation, and b_D < -0.0173 for soils draining quickly and b > -0.0173 for soils with slow drainage.

(ii) capillary rise computed using the parametric equations developed by Liu et al. (2006) which parameters (a₁, b₁, a₂, b₂, a₃, b₃, a₄ and b₄) depend on the soil’s physical and hydraulic characteristics, were subjected to the model’s calibration process [41,71] (Table 1); and (iii) surface runoff which is estimated with the curve number (CN, dimensionless) method [72]. CN ranges from 0 to 100, meaning there is no runoff, and all rainfall becomes runoff, respectively. Runoff (RO) was estimated using the curve number (CN) approach [51,72], in which CN value changes depending on the soil and vegetation type, and the antecedent moisture in soil. SIMDualKc model adjusts CN daily to represent the influences of adding or subtracting soil water content on soil infiltration properties, by linking the CN value to soil water depletion in the surface layer.

A complete description of the calibration and testing of the SIMDualKc model is provided by [41,42,73,74].

2.5. Estimating K_cb values using the A&P approach and remote sensing data

In the A&P approach [48], K_{cb A&P} is computed with the equation below [41,48] where impacts of plant density and/or leaf area are taken into consideration as a function of a density coefficient (K_d):

$${\mathrm{K}}_{\mathrm{c}\mathrm{b}\mathrm{A}\&\mathrm{P}}={\mathrm{K}}_{\mathrm{c}\mathrm{m}\mathrm{i}\mathrm{n}}+{\mathrm{K}}_{\mathrm{d}}\left({\mathrm{K}}_{\mathrm{c}\mathrm{b}\mathrm{f}\mathrm{u}\mathrm{l}\mathrm{l}}-{\mathrm{K}}_{\mathrm{c}\mathrm{m}\mathrm{i}\mathrm{n}}\right)$$

(9)

To consider the influence of crop height and density on K_cb, [48] proposed a density coefficient (K_d) (equation 9) that allows adjusting the K_cb of vegetation and ground cover conditions [48,51]. This K_d is calculated through the factors: f_{c eff} that corresponds to the fraction of ground effectively covered or shaded by plants at midday, M_L, which is a multiplicative factor of f_{c eff}, representative of the fraction of ET per unit area of horizontal vegetation to the ET_o on the same area, and h that is the vegetation height.

$${\mathrm{K}}_{\mathrm{d}}=min\left(1,\text{ ; }{M}_L{f}_{ceff},f_{ceff}^{\left({\displaystyle \frac{1}{1+h}}\right)}\right)$$

(10)

K_d has also been used in other studies to adjust K_cb to actual density conditions for full cover crops [35] and in partial cover crops using both ground [75] or remote sensing data [49,50,76].

The K_{cb full} represents a ceiling limit on K_{cb mid} for vegetation under no-water stress and with a full ground cover and LAI value higher than 3. Following [27], the K_{cb full} is estimated as a function of mean plant height (h) and corrected for climate conditions as:

$${\mathrm{K}}_{\mathrm{c}\mathrm{b}\mathrm{f}\mathrm{u}\mathrm{l}\mathrm{l}}={\mathrm{F}}_{\mathrm{r}}\left(\mathrm{m}\mathrm{i}\mathrm{n}\left(1.0+{\mathrm{k}}_{\mathrm{h}},1.20\right)+\left[0.04\left({\mathrm{u}}_2-2\right)-0.004\left({\mathrm{R}\mathrm{H}}_{\mathrm{m}\mathrm{i}\mathrm{n}}-45\right){\left({\displaystyle \frac{\mathrm{h}}{3}}\right)}^{0.3}\right]\right)$$

(11)

wherever u₂ is the daily mean wind speed (m s^-1) considered of 2 m above ground level through the growth period, RH_min (%) is the daily mean of minimum relative humidity over the growth cycle, and h is the average plant height (m) for the mid-season stage. Previously to the climatic correction, a value of 1.20 is considered as a superior threshold for K_{cb full}. The crop height influence is considered in the sum 1 + k_h h, with k_h = 0.1 for tree and vine crops, as recommended by [28,35].

Taller crops and situations where the local climate is windier or drier than typical (standard) climatic conditions (RH_min = 45% and u₂ = 2 m s^-1) are expected to have higher total K_{cb full} values. When vegetation shows greater stomatal adjustment transpiration, as is common in most annual crops, the Fr parameter applies an empirical correction (F_r ≤ 1.0). When crops show strong vegetative vigor and decline consequently by training and pruning, as well as when water availability is restricted, F_r is high for trees and vines. As proposed by [27], F_r, can be estimated as:

$${\mathrm{F}}_{\mathrm{r}}={\displaystyle \frac{∆+\gamma \left(1+0.34u_2\right)}{∆+\gamma \left(1+0.34u_2\frac{r_l}{r_{typ}}\right)}}$$

(12)

where r_l and r_typ are, respectively, the average of leaf resistance and the typical leaf resistance (s m^-1) for the considered plant, Δ is the slope of the saturation vapor pressure vs. air temperature curve, kPa, and γ is the psychrometric constant, kPa °C^-1, both for the same period when K_{cb full} is calculated. The earliest established version of that equation considered a fixed r_typ = 100 s m^-1, a normal value for yearly crops, however default F_r values for trees and vines and several annual crops were recently revised and updated by [28].

[35] provided examples of the K_{cb A&P} approach's application to a variety of annual and perennial crops, as well as active soil cover in the interrow for wine vineyards. The K_{cb A&P} technique does not require any calibration/test process when using the tabulated parameters [28,31].

Alternatively, in the A&P approach, the K_d can be estimated under actual conditions using remote sensing data of the fraction of soil covered by vegetation (f_{c VI}). In the present study, we have assumed that f_{c eff} is approximately equal to the f_{c VI} derived from spectral vegetation indices computed with RS data. This data was derived from spectral vegetation indices obtained through satellite data, namely Landsat 5 images for the period under analysis. Five Landsat-5 images covering the study area, relative to the path/row 204-033, were available and downloaded from Earth Explorer portal (https://earthexplorer.usgs.gov/), for the study period (dates 26/04/1987, 28/05/1987, 29/06/1987, 31/07/1987, and 17/09/1987). From each of these images, 12 pixels covering the vineyard under study were selected.

The Soil Adjusted Vegetation Index (SAVI) was considered to calculate the f_{c VI}, as it presents less sensitive to soil background reflectance variability than other VI [49,50].

$$\mathrm{S}\mathrm{A}\mathrm{V}\mathrm{I}={\displaystyle \frac{\left({\rho }_{NIR}-{\rho }_{red}\right)\left(1+L\right)}{\left({\rho }_{NIR}-{\rho }_{red}+L\right)}}$$

(13)

where the ρ_red and ρ_NIR are the reflectance at red and NIR spectral domains, and L is the soil conditioning factor varying between 0 and 1. A value of L equal to 0.5 is considered for the most common environmental conditions and was found to minimize the effects of soil brightness variation and eliminate the need for calibration under different soil conditions [49,50,77].

The f_{c VI} was computed as:

$${\mathrm{f}}_{\mathrm{c}\mathrm{V}\mathrm{I}}={\displaystyle \frac{\mathrm{S}\mathrm{A}\mathrm{V}\mathrm{I}-\mathrm{S}\mathrm{A}{\mathrm{V}\mathrm{I}}_{\mathrm{m}\mathrm{i}\mathrm{n}}}{{\mathrm{S}\mathrm{A}\mathrm{V}\mathrm{I}}_{\mathrm{m}\mathrm{a}\mathrm{x}}-{\mathrm{S}\mathrm{A}\mathrm{V}\mathrm{I}}_{\mathrm{m}\mathrm{i}\mathrm{n}}}}$$

(14)

where f_c is the fraction of ground cover by vegetation, SAVI, SAVI_max and SAVI_min correspond to the vegetation index for a specific date and pixel, vegetation index at maximum vegetation cover, and vegetation index at minimum vegetation cover (bare soil), respectively. The SAVI_max and SAVI_min values were obtained from the literature considering woody crops [78] and [49,50]. In the present study, the SAVI_max and SAVI_min values considered were 0.75 and 0.09, respectively.

When calculating K_d, it was assumed that the fraction of soil effectively covered at solar noon (f_{c eff}) was equal to the fraction of soil covered by vegetation (f_c). The M_L factor, used to calculate K_d was not changed and remained equal to 1.5 [48].

2.4.3. Parameterization and calibration procedures of SIMDualKc model

The first step in the calibration procedure was to define the upper and lower boundary conditions. A value of 1.20 is assumed as an upper limit for K_{cb full}, before climate adjustment. As above mentioned, for different conditions from those standards, i.e., considering taller crops and local climate windier/drier or windier than RH_min = 45% and u₂ = 2 m s^-1, higher K^{cb full} values can be obtained. The lower boundary conditions are K_{cb ini} equal to 0.15, which corresponds to bare soil conditions.

The observed parameters of weather data, dates of the crop growth stage, crop height (h), root depth (Z_r), and soil characteristics were input into the modelling tool. The weather data used were the daily maximum and minimum temperatures (°C), T_max and T_min, respectively, the daily average wind speed, measured by an anemometer at a height of 6 meters (km h^-1), the daily average relative humidity (%), RH_med, the daily number of sunshine hours (h) and the daily precipitation (mm), corresponding to the year of 1987.

The initial values for the non-observed parameters (K_cb, depletion factor p, a_D, b_D) were set using diverse sources of information. Based on the vineyard density, the training system and the f_c, the K_cb value was taken from [28]; for the depletion factor, the soil evaporative layer (Total Evaporable Water - TEW, Readily Evaporable Water - REW, evaporable layer - Z_e), and a_D, the values were calculated according to the soil textural characteristics (sandy loam soil); and for b_D, values were retrieved from [71].

The calibration procedure followed that described by [73]. Thus, the initial K_cb values used in the simulation were those tabulated by [28] for vineyards. The p factor was considered as a constant value throughout the cycle, considering that the study was performed in a 20-year rainfed vineyard and that, during these stages, the activation of deeper roots is to be expected, increasing the plant's ability to cope with water scarcity, as discussed by [79].

Regarding the evaporable soil layer (Z_e), the initial values for TEW and REW were derived from the characteristics of the soil, i.e., field capacity and the permanent wilting point values, and the textural characteristics, for the selected thickness of the evaporative layer.

In the following procedures, the DP, CN and CR parameter values were adjusted. The initial values for the deep percolation parameters, a_D and b_D were respectively 275 and -0.0173. CN value was assumed equal to 68 [72], according to soil characteristics. The water table depth (WTD) varied from 4.2 m to 4.61 m during the crop season. The measured values of WTD were taken from a piezometric station nearby the experimental plot which belongs to Portugal's National Water Resources Information System (https://snirh.apambiente.pt/). The values of the parameters a₁, a₂, b₁, b₂, a₃, b₃, a₄ and b₄ that vary according to soil texture were based on those proposed by [71] for loamy sandy soils. The maximum root depth was based upon measurements taken in open soil profiles in the studied vineyard, thus Z_{r max} = 1.85 m.

Since the weeds in the vineyard under study were controlled by tillage and herbicide, respectively in the inter-row and row, neither mulching, nor active ground cover were considered in the simulations.

The model was run until no further improvements were obtained and the best values of the variables were retained.

As already mentioned, two sets of measured SWC were used for the calibration and test of the SIMDualKc model. The data set called 1D (Figure 1) on SWC was used for calibration as it provided 10 observations, while the data set called 2A was chosen for testing (n=7). Due to restrictions in observed data sets it was therefore not possible to validate the model, but the calibrated parameters were tested using measured SWC from a different location within the same field.

2.4.4. Modeling tool accuracy assessment

The accuracy of the calibration and test was assessed using a set of goodness-of-fit indicators described below. These indicators have been used in other modeling performance evaluations, e.g. [45,44], and [47] for vineyard, [31] to evaluate olive grove, or by [30,80].

The calibration/test procedure was considered acceptable once the test results indicators’ accuracy was judged comparable to those obtained for the calibration dataset. The selected goodness-of-fit indicators are detailed in [73], including the targeted values, and include:

i) the regression coefficient (b₀, dimensionless) of the linear regression pushed to the origin between the measured and simulated variables; the goal value is b₀ = 1.0, which implies that the projected values are statistically equivalent to the observed ones;

ii) the coefficient of determination (R², dimensionless) of the ordinary least-squares regression between the simulated and observed values should be close to 1.0, indicating that the model explains the majority of the variance in the observed values;

iii) root mean square error (RMSE, mm), the desired value of which is zero, demonstrating a perfect match between simulated and observed variables, the value of which must be much smaller than the mean of the observed values;

iv) the normalized root mean square error (NRMSE, %), whose target value is zero, signifying a perfect match between simulated and observed variables, must be much smaller than the mean of the observed values;

v) the average absolute error (AAE, mm day^-1), represents the average error related to the estimations and ought to be significantly less than the average of the observed values;

vi) the target value of average relative error (ARE), which represents the magnitude of the error in relative terms, is zero, signifying a perfect match between the values of the simulated and observed variables;

vii) the ratio of the mean square error (MSE) to the variance of the observed variable is the Nash and Sutcliffe efficiency of modelling (EF, dimensionless) [81,82], whose target value of 1.0 denotes that MSE is negligible for observed variables.

3. Results and Discussion

3.1. Performance of the SIMDualKc Model in Calculating Soil Water Content

The crop characterization in the SIMDualKc model was done by entering the dates of its vegetative cycle with an indication of root depth, plant height (h), and soil water depletion fraction (p). The characterization of the crop relative to the K_cb values for the different phases of the cycle (initial, intermediate, and final), and soil characteristics (soil evaporation layer, runoff and deep percolation, capillary rise) are shown on Table 2. These initial values of K_cb, h and p were from [28]. The initial p value, proposed by FAO [27] was adjusted during calibration since this is a rainfed vineyard adapted to water stress.

Results in Table 2 show that, for most cases, the values of the parameters used for initializing the model are relatively close to those obtained after proper model calibration (the variation was 12.2% on average, with K_{cb end} and p factors responding for the higher variation, respectively increasing 23% and 25%). The depletion fraction p's calibrated values (Table 2) are greater than those suggested by [27,28] which relates to the fact that the vineyard was rainfed and, therefore, well adapted to low water availability conditions. The calibrated p value is in the range of those reported for rainfed vineyard [44,45]. A slightly lower p value was reported by [83] for irrigated cv. Loureiro vineyards. Differently, [55] reported a p value of 0.65 for a drip irrigated vineyard.

The textural characteristics of the soil and the water-retention capacity of the soil evaporation layer were used to estimate the initial values for Z_e, REW, and TEW calculations, as proposed in [27]. In this way, the values of these calibrated soil evaporation parameters are the same as those indicated by [27] for sandy loam textured soils. Similarly, the calibrated runoff parameter CN was kept unchanged from the initial value proposed by [51] and [72], considering the soil texture. One of the parameters associated with deep percolation, a_D, was estimated from the water available in the soil, at saturation and field capacity, and it was slightly changed from the initial values, decreasing 3.5%. The b_D parameter, estimated from the drainage characteristics of the soil [71], was left unchanged.

Figure 5 shows the simulated and observed available soil water (ASW) dynamics throughout the season. Results in Figure 5 show, that the model is able to adequately simulate the dynamics of the ASW observations. The differences between observations and simulations in the calibration set (Figure 5a), were lower in the intermediate and higher in the final stage of the crop. In other words, the ASW dynamics simulated with SIMDualKc agree well with the ASW observations. In addition, the results agree with the predawn leaf water potential (ψb) observations. The model simulates water stress from DOY 136 onwards while the ψb observations range from -0.2 to -0.4 MPa, thus presenting a slight water stress (Figure 6). The decreasing in the ψb, increase in the plant water stress, well matches the simulations of ASW from then onwards. The ψb reached values close to -1.0 MPa in the second week of September, and the daily minimum leaf water potential varied between -1.4 MPa and -1.6 MPa [19], which agrees well with the intense soil water stress simulated with SIMDualKc.

Rainfed vineyards show similar ψb dynamics during the June-August period in south Portugal (typical mediterranean climatic conditions) i.e., the occurrence of slight to moderate water stress, reaching severe water stress (-1.56 MPa and -1.96 MPa) by September as reported by [84]. Differently, in a rainfed vineyard in Douro Wine Region, located in the North of Portugal, where climatic demand is lower during July-September, [24] reported ψb ranging from -0.3 to -0.5 MPa, indicating slight to moderate water stress.

In a study developed in rainfed vineyards cv. Chardonnay in Catalonia, Spain [85,86] low ψb values up to -1.36 MPa. Similar finds were reported by [87], and [88], for mediterranean climatic conditions.

Figure 5 shows a slight mismatch for the extremes of the period considered, suggesting a more abrupt progression over time of the measured values. However, the statistical indicators showed a good performance of the model (Table 3). Both, the calibration and the test comparing available soil water simulated by the SIMDualKc model and available soil water measured using a neutron probe have a good correlation, with a regression coefficient (b₀) of 0.97 (Table 3). The coefficient of determination (R²) is equal to 1 in both cases, which indicates that the model can explain the variance of the data. Regarding the modelling efficiency (EF) the values also show a good result with values equal to or higher than 0.97, which indicates that the residual variance due to modelling is comparable to the variance of the measured data, as depicted in [81]. The values of root mean square error were quite similar in both the calibration and test, with RMSE of 11.1 and 11.9 mm respectively, while the mean absolute error (AAE, mm) was 9.5 to 10.2, respectively. Additionally, the NRMSE presents a good result with values below 13%. Similar to the other indicators, ARE showed a good agreement between simulated and observed SWC values, with an error of less than 0.60.

For vineyards, there are only a few results in the literature regarding the calibration of the SIMDualKc model using soil water content data. The RMSE and NRMSE (11.5 mm and 12%, respectively) obtained in the present study are higher than those verified by [83] for a rainfed vineyard cv. Loureiro in Ponte de Lima, northwest Portugal, by [44], working in a vineyard cv. Godello and cv. Mencía in Galicia, Spain and also higher than those estimated by [45], for a vineyard cv. Albariño in Pontevedra, Spain. Likewise in the present study, a small error (RMSE) was reported by [89] simulating water flow within the Hydrus (2D/3D) model towards a two-year set of soil water content measurements from a rainfed vineyard cv. Aglianico. Similarly, [90] found small values of RMSE comparing soil water content measured with TDR probes and estimated using the SWAP model for irrigated crops (passion fruit and pastures). On the contrary, [91] found RMSE values higher than those reported in the present study assessing soil moisture profiles with HYDRUS-1D model for several Mediterranean rainfed vineyards. Although the values of RMSE and NRMSE reported in the present study indicate errors of estimation, they represent less than 4.17% and 4.35% of the TAW for RMSE and NRMSE, respectively. [47] also obtained RMSE and NRMSE lower than the values verified in the present study for a rainfed vineyard in Monferrato, northern Italy. In the reported studies [44,45,47,83], the soil depth was 1.20 m, 1.0 m, 0.60 m, and 0.80 m, respectively, therefore lower than that considered in the present study (1.85 m). Nonetheless, given the very different edaphoclimatic characteristics, in places where annual rainfall is higher, the root system is expected to explore more superficial layers and the lower limit of the control volume in the modelling process can be located at a depth closer to the surface than in the current situation, unlike in places with lower rainfall and deeper soils. Moreover, in situations where the plants are under greater water stress, either due to more demanding environmental conditions or to the absence of irrigation (as in the present study), there will be a tendency for the root system to grow deeper, easily surpassing the common threshold (approx. 1.0 m) as depicted by several authors [4,5,6,7,23], sometimes deeper than 2 m (for an irrigated vineyard under extremely arid conditions) [92], or even higher than 6 m, as reported by [7].

3.2. Crop Coefficients Dynamics over the Season

The daily values of precipitation, the potential basal crop coefficient (K_cb), the actual basal crop coefficient (K_{cb act}), the soil evaporation coefficient (K_e), and the actual crop coefficient (K_{c act} = K_{cb act} + K_e) calculated by SIMDualKc model are shown in Figure 7. For both, calibration and test, the K_{cb act} and the K_cb values are coincident during most of the crop cycle, except for part of the mid-season (DOY 196), till harvesting, when rainfall events are not sufficient to avoid water stress (as seen in Figure 6), i.e., few rainfall events occurred during spring, winter, and autumn. The value of K_{cb ini} calibrated is equal to the standard value proposed in [28] for vineyards (0.15). However, during mid-season and late-season, the calibrated values of K_cb are quite different from standard values (K_{cb mid} = 0.65 and K_{cb end} = 0.40), being K_{cb mid} = 0.60 and K_{cb end} = 0.52 (as shown in Table 2). These differences are related to the rainfed conditions of the studied vineyard, which reflect its tolerance to water stress. The calibrated K_cb values (K_{cb ini}, K_{cb mid}, and K_{cb end}, respectively) in the present study are lower than those reported by [44] for vineyard cv. Godello (0.30, 0.80, and 0.60) and cv. Mencía (0.20, 0.75, and 0.60) in Galicia, Spain, probably related to the higher annual rainfall amount in that location. In the same way, the K_cb calibrated here are lower than those obtained by [45] for a vineyard cv. Albariño, under irrigation conditions, which can explain the differences in comparison to the values verified in the present study. The results of K_cb obtained in [45] for the K_{cb mid} (mean 0.57) with moderate stress, are quite similar to the K_{cb mid} values calibrated in the present study. [83] also reported calibrated K_cb values ranging from 0.09-0.13, 0.21-0.25, and 0.27-0.22 (for development, mid-season and late-season) for a vineyard cv. Loureiro under rainfed treatment, which are lower than those calibrated in the present study. The difference can be related to higher amounts of precipitation in the location reported in [83].

As verified for K_cb and K_{cb act}, the K_{c act} and K_e have a similar pattern, especially in the initial and final stages. The K_{c act} peaks occur mainly when the soil is wetted, after rainfall events, and the evapotranspiration process is dominated by soil evaporation, i.e., the transpiration amount is small in this period. Furthermore, in Figure 7, it is possible to find times (especially in the early growth phases) when the K_e values are very similar to the K_{c act} values, close to 1.2, in the end season with K_{c act} close to 1.0, and close to harvest, when the K_{c act} values are close to and higher than 1.2. This pattern is mainly related to the rainfall during these periods (Figure 7), as well as the absence of ground cover vegetation. Similar trends in K_e and K_{c act} values were observed by [44], for the cultivars cv. Godello and Mencía in Galicia, NW Spain and [83] to cv. Loureiro in Ponte de Lima, Portugal.

Numerous K_e peaks on the soil evaporation curves (Figure 7) show how the vineyard responded to rainfall events at its distinctive growth stages. These K_e peaks occur mostly at the initial, start of rapid growth and at the beginning of mid-season periods. K_e peaks respond to the total soil wetting that occurs from rain events (f_w = 1). During almost all mid-season, maturity and harvesting stages, K_e peaks are small and not very frequent due to the dry soil conditions with the absence of rain. This pattern K_e response to rainfall was also reported by [44,45] and [83] for vineyards cv. Godello and cv Mencía, cv. Albariño, and cv. Loureiro, respectively, under a rainfed treatment in Spain and Portugal.

3.3. Estimation of the Fraction of Ground Cover from Vegetation Indices

The fc retrieved from SAVI (f_{c VI}; Equation (14)) was estimated for five dates, those with cloud-free images for the period under study (Table 4). The lowest f_{c VI} value was estimated for the beginning of the growing season (DOY 116, f_c 0.174) and the highest f_{c VI} was verified for the mid-season (DOY 148, fc 0.286). The temporal evolution of SAVI is similar to the f_{c VI}, with the lowest SAVI value estimated for the beginning of the growing season (start rapid growth) and the highest value observed in the mid-season. This tendency was observed for table grapes using NDVI [93], and for other crops and VIs, like maize, barley, and olive, using SAVI and NDVI [50].

The f_{c VI} obtained in the present study (ranging from 0.17 to 0.28, Table 4) are lower than those reported by [93] for a f_{c VI} relationship of table grapes vineyards cv. Perlette and Superior in the semi-arid region of Northwest Mexico. However, [93] used NDVI for the computation of f_{c VI} instead of NDVI. Additionally, the f_c values of the present study are closer to the limit values presented in [27] for vines and other perennial crops. The authors reported f_c values ranging from 0.25 to 0.45, with K_{cb mid} values between 0.46 and 0.80 and K_{cb end} values varying from 0.20 to 0.60. On the other hand, [10] reported f_c values ranging from 0.04 to 0.28 and from 0.06 to 0.33, for 2017 and 2018 seasons, respectively, in a rainfed vineyard cv. Monastrell, in Albacete, Spain, with maximum values similar to those observed in the present study. [47] reported f_c values ranging from 0.15 to 0.28 for an Italian rainfed vineyard, which is similar compared to the f_c values verified in the present study. It important to take into account the spatial resolution of the Landsat-5 images available at the time of the study, 30 m x 30 m, which means that no more than 12 pixels per image can be used to estimate the f_c for the vineyard studied.

3.4. Comparison between K_cb obtained with SIMDualKc model and predicted with the A&P approach

The K_{cb A&P} values obtained through f_c, derived from remote sensing (f_{c VI}), and applying equation 9 [48], were compared with the K_cb obtained through the SIMDualKc model, to verify the accuracy of using the A&P approach combined with RS data. The five K_cb values obtained by SIMDualKc, for each day where K_{cb A&P} was estimated using the f_{c VI}, were compared.

Table 5 shows that the K_{cb A&P} values are close to the values obtained by the calibrated SIMDualKc (K_{cb SIMDualKc_1D}), especially for 26/05 and 29/06/1987. For these dates, the plants were not under water stress (shown in Figure 6) and both, SIMDualKc model and A&P approach, were able to precisely estimate the K_cb values. On the contrary, the K_{cb A&P} values for the final stages (after July 31) are greater than those for K_{cb SIMDualKc}. It is noteworthy that vegetation indices (VI) detect crop stress when a reduction in chlorophyll green cover or changes in canopy geometry occur, and thus do not immediately detect the slight water stress conditions in plants, when the effect is only noticeable in the transpiration rates and predawn water potential (Figure 6). During the two first weeks of August, the water stress increase was moderate and after the third week onwards it was severe (Figure 6), which can explain the higher deviation between K_{cb SIMDualKc} and K_{cb A&P} values occurred on September 17 (Table 5). For these stress conditions, it is necessary to add ancillary modelling approaches based on soil water balance, as also reported by [94] and [50], or based on a thermal signal, as demonstrated by [95], or even using a combined K_cb-VI approach, with soil water balance for assessing actual K_c, as demonstrated in [49,50].

Regarding the test process, the K_{cb A&P} and K_{cb SIMDualKc} values are very similar and the deviation between them is low, except for the value corresponding to the late-season stage, following the same pattern observed in the calibration (Table 5). As shown in Table 5, the deviation between K_{cb A&P} and K_{cb SIMDualKc} values starts when plants suffer slight water stress (confirmed by Figure 6), and the highest difference between K_{cb A&P} and K_{cb SIMDualKc_2A} occurs for the last crop stage, when the plant water stress is moderate, as also observed for calibration (K_{cb SIMDualKc_1D}).

The modelling efficiency was low for test (poor results with values between 0.03 and 0.38) and it presented a lower coefficient of determination (R²) than model calibration. However, the low number of Landsat-5 images available for the study period somehow affects the results of the quality of fit indicators. Nevertheless, the difference between the K_{cb A&P} and the K_{cb SIMDualKc} is small for most of dates, mainly when the plants are not under water stress. Differently, when moderate (or severe) water stress starts, (after 31/07/1987), the difference between the K_{cb A&P} and the K_{cb SIMDualKc} values is greater (specially for 31/07/1987 and 17/09/1987).

It is also worth mentioning that, while the results obtained with SIMDualKc refer to specific positions in the plot, namely the measurement position considered for the model calibration (1D) and the position considered for the test (2A), the results for the A&P approach, use the f_{c VI} values from SAVI derived from satellite images. The Landsat-5 satellite images, available in the study period, have a pixel size of 30 m x 30 m, thus representing an area on the ground where the variability of soil conditions is higher than the variability of the soil water content measurement positions.

The SIMDualKc model and the A&P approach were not able to explain the entire variance of the data, due to the small number of observations of the f_{c VI} data. This small number of available satellite images was pointed as a limitation in a study by [55]. Moreover, reduced efficacy in the A&P approach as a result of insufficient fc and h field observations was also reported by [31]. The current availability of a large number of satellite missions with improved spatial and temporal resolutions (e.g., Landsat 8, Landsat 9, Sentinel-2) can help overcome this data availability limitation.

3.5. Water Balance and Respective Components

Figure 8 summarizes the various soil water balance components for the studied year that were estimated using SIMDualKc model. For calibration (1D), regarding the runoff process (RO), it occurred only for two periods, in the initial crop stage (3.60 mm), and at the late season (with 7.70 mm). The capillary rise (CR) contribution to soil water balance, for the calibration process was estimated as a total of 25.16 mm from development till late-season stage, with 52% of the total CR occurring on mid-season. For test (2A), the CR contribution was 15.80 mm, from the mid-season up to the late-season stage. The difference in CR between the two positions (1D and 2A) can be explained by the variability of the soil in the vineyard area, which is composed, as described by [19], of calcareous and non-calcareous materials in the form of complexes, such as fine sandstones, clays, and argillites with a texture variation from loamy to clayey. Besides that, in this area, there is a variation in the slope of the terrain (approximately 2.5%), and position 1D (used for to calibrate the model) is located in a higher landscape position than position 2A (used for testing the model). Naturally, it is expected that greater accumulation of water in the soil will occur in lower topographical positions, associated with the soil's water retention capacity, as demonstrated in Figures 5 (available soil water) and 6 (predawn leaf water potential). Furthermore, [19] mentioned that this is an area intensely marked by crops and with frequent annual plowing up to 0.20-0.30 m, sometimes reaching 0.65-0.80 m. Thus, the natural soil horizons are mixed to a greater or lesser extent. Therefore, taking this into account and as the soil is an Anthrosol, it is understandable to have this variation in CR, as verified in the present study. Even so, it is worth mentioning that there was no significant difference in vine yield in the year under study between the two positions evaluated, with 5.22 kg plant^-1 in 1D and 3.80 kg plant^-1 in 2A, as described in [19].

Variability in precipitation and its distribution across the different stages of crop growth were verified for the study year (1987), which had less precipitation and was drier than the previous five years' average (1982–1987) (as depicted in Figure 3). The runoff (RO, mm) represents 5% of precipitation, and capillary rise (CR, mm) represents 7.02% (position 1D) and 11.18% (position 2A), respectively. These percentages show that rainfall was primarily utilized for crop evapotranspiration. The highest quantities of RO and CR were verified at the maturity and end-season crop phases, when crop evapotranspiration was lowest (Table 6).

The soil water balance values (average values for calibration and test) estimated with the SIMDualKc model are summarized in Table 6. The actual evapotranspiration (ET_{c act}) ranges from 68 mm to 55 mm, in the initiation stage to late-season, respectively. The higher value is observed in the mid-season (143 mm), which represents 40.4% of the total in the crop cycle. Most of the ET_{c act} (84.7%) occurs between the beginning and mid-season stages, with higher precipitation in March and April having a significant effect on this. A similar pattern occurs for the actual transpiration (T_{c act}), with 82% occurring from initiation to mid-season stages. Since the vineyard was not irrigated, the contribution from precipitation and the capacity of the soil to store water (Δ SWC, Figure 8), and to contribute to ET_{c act} (including capillary rise), are higher during mid-season, indicating that the water stored in the soil (especially in the deeper layers) played an important role in the actual crop evapotranspiration for the studied area. The daily actual transpiration (T_{c act}) and soil evaporation (E_s) for the year of study estimated with SIMDualKc model are presented on supplementary Figure S1 (Annex I).

By analyzing the consumptive use of water (Table 6), it can be observed that in the ET_{c act} partitioning, the soil evaporation (E_s) was the main component of the evapotranspiration process, particularly during the initial and development crop phases. E_s values were 83.8% and 51.6% of the total ET_c act on the aforementioned stages. This variation is related to higher precipitation on these periods, which kept the soil wet. On the other hand, E_s contribution to ET_{c act} was smaller than T_{c act} at mid-season and late-season stages, with values ranging from 3.5% to 23.6% of ET_{c act}. During these stages, the transpiration becomes more important and is the main consumptive water use, representing 97.2% and 74.54% of ET_{c act}. Once again, these values are explained by the influence of precipitation, soil water storage and capillary rise, therefore the rainfed condition of the vineyard. In general, it was found that E_s and T_{c act} values represents 36% and 64% of ET_{c act} during the entire crop season. Similar results for evapotranspiration partition occurring during the mid-season and, however less important, during the late season, were observed by [10] and by [47] for rainfed vineyards. [10] studied vineyards under two watering regimes that are frequently applied in the Mediterranean region for this crop (rainfed and deficit irrigation) and [47] evaluated a rainfed Italian vineyard, also in a Mediterranean condition. The differences between these studies and the present study can be explained by different amounts of rain and soil water storage capacity of the respective locations.

Results regarding the ET_{c act}/ET_c ratio (Table 6) show that ET_{c act} was smaller by 17% relative to the potential ET_c in the studied year. During the mid-season, that ratio ET_{c act}/ET_c was small, when ET_{c act} was smaller by 13% comparatively to the potential ET_c. These reductions are probably low in non-irrigated scenarios with low rainfall, as analyzed above, which did not affect yields (as reported by [19]), comparing calibration and test positions.

4. Conclusions

In the present study, soil water balance modelling using SIMDualKc was evaluated in a rainfed vineyard, in Santarém, Portugal, using as reference soil water content field data. These data were collected at a depth (1.85 m) greater than that commonly used for modelling and therefore, allowing to better represent the soil volume explored by the roots. Additionally, the actual crop coefficient (K_{cb act}) obtained in SIMDualKc model (K_{cb SIMDualKc}) was compared with K_cb estimated using the A&P approach [48] (K_{cb A&P}), where fc was estimated from spectral vegetation indices derived from Landsat 5 images to further estimate K_d.

For the specific characteristics of this study area, namely, the soil, the climate, and the studied vineyard cultural practices (spacing, crop height, inter-row management), the calibration and test of the SIMDualKc model were successfully performed, proving that it is possible to calibrate the model considering a soil profile depth greater than those that are normally used for vineyard ET modelling with soil water content.

Previous studies that adopted a similar calibration procedure of SIMDualKc with soil water content values were carried out in areas with higher annual rainfall and milder temperatures than those recorded in the region of Santarém. Although these studies provided good results for the application of SIMDualKc model, this was performed for less deeper soil profiles (up to 1.0 m).

For the study year (1987), a good fit was achieved between the soil water content values measured in the field with the neutron probe and the values simulated with the SIMDualKc model, both for the calibration and the test performed.

Furthermore, the comparison between K_{cb SIMDualKc} and the K_{cb A&P} showed a small difference for most values throughout the different crop stages. However, the reduced number of satellite images to estimate f_{c VI} and derive K_d in the A&P approach was a limitation in obtaining a good fit between the K_{cb act} simulated by the SIMDualKc model and the K_{cb act} estimated by the A&P approach. Nevertheless, nowadays increasingly available satellite missions can potentially contribute to better approximating the soil water balance modeling to the real conditions of vegetation development.

We can reckon that the present study allowed improving the knowledge about the application of the SIMDualKc model in vineyards, exploring the application of the model for conditions of increased root depth, as expected in rainfed conditions under Mediterranean influence, and thus potentially contributing to better water management.

Moreover, model fitting results (namely, modelling efficiency and variance explained by the model) were very good, indicating an enhanced accuracy of the model when using increased depth for calibration field data of vineyards.