Water Use and Soil Water Balance of Mediterranean Vineyards under Rainfed and Drip Irrigation Management: Evapotranspiration Partition and Soil Management Modelling for Resource Conservation

: Vineyards represent complex Mediterranean agrosystems that deliver signiﬁcant ecosystem services to society. Yet, many vine-growers still need to assimilate the importance of crop and soil management to the conservation of soil and water resources. The main objective of this study was to evaluate water use and the water balance terms in rainfed and irrigated vineyards in Italy and Portugal, respectively, in both cases aiming at the sustainability of natural resources use. The SIMDualKc model is used for both sites after calibration and validation by ﬁtting soil water content measurements. The Italian case study focused on the impacts of inter-row conservation management in hillslope vineyards while the Portuguese case study analyzed irrigation water management under scarcity in ﬂat vineyards. For the Italian vineyards, the model results focused on the evapotranspiration ﬂuxes and their partition, control of surface runoff, and soil water recharge provided by the inter-row soil management using cover crops. Model results of the Portuguese case study showed the need for improving irrigation water use and the terms of water balance, namely referring to percolation and soil water evaporation. Both case studies further demonstrated the advantages of using computational tools to better cope with climate variability in the Mediterranean region and made evident the beneﬁts of improved crop and soil management practices in counteracting land degradation and valuing the use and conservation of natural resources.


Introduction
Viticulture is one of the most diffused cultivations in the world and has been practiced in the Mediterranean area for millennia [1]. In this region, the vineyard agricultural system is potentially well suited for delivering provisioning services such as grapes for table or wine production, the regulation of climate and the hydrologic cycle [2], and the preservation and enhancement of cultural heritage including landscape and aesthetic values [3,4]. Contrasting, vineyards are often associated with several environmental problems resulting from the intensification of production systems, which evidence the need for better management of soil and water resources.
it is crucial to know the response of vineyards to rainfall distribution in terms of water availability to plants when considering the presence of cover crops in the inter-row, both in rainfed and irrigated vineyards, since the dynamics of evapotranspiration and water use is insufficiently known in both cases. Moreover, related measures and practices are among those required for adaptation and resilience to climate change impacts.
Costa et al. [27] assumed a long-term perspective when focusing on genetics, selection of varieties and plant amelioration of both the grapevines and the rootstocks, aimed at improving grapevine responses to heat stress and drought, thus in addition to the measures and practices reported. These authors largely referred to policies required to make effective the development of farm responses of vines to heat and drought and focused on water issues. Moreover, they assumed that future strategies to optimize the environmental performance of the wine sector in the Mediterranean must be focused on water and irrigation. Hannah et al. [28] also assumed that attempting to maintain wine grape productivity and quality under climate change implies increased water use for irrigation and to cool grapes through misting or sprinkling, which creates potential impacts on freshwater conservation. Thus, freshwater habitats may be affected where climate change undermines growing conditions for already established vineyards. Hannah et al. [28] concluded that climate change adaptation strategies are required for creating a positive future for producers, wine makers, and vine ecosystems.
In Italy and Portugal, which represent the 1st and 5th largest wine producers in Europe [29], vineyards were traditionally rainfed since, as in many European countries, irrigation was historically not used, forbidden by regulations for many quality wines. Nowadays, in Italy, many IGT, DOC and DOCG regulations admit only emergency irrigation. However, since the soil water stress strongly affects the growth and production of vines and the quality of berries [30], irrigation has become an increasingly frequent practice especially in dry areas of Southern Europe [31,32]. Yet, poor irrigation practices, namely the inadequacy of irrigation depths to crop water needs, climate, and soil and irrigation system characteristics may cause water percolation and leaching of fertilizers and pesticides, thus not promoting environmental friendliness [33][34][35], nor greater and more stable yields and quality of wine [16]. Better knowledge of the crop water requirements relative to the various crop growth stages are critical to better preserve the vineyards natural resources, namely soil and water.
Accurate knowledge of the soil water balance is fundamental for improving the soil and water management of the vineyard system and to further cope with the challenges resulting from climate change. In Piedmont, Italy, long-term simulations performed over 60 years (1950-2009) confirmed that the climate change is already influencing local vineyards since 1980 [36,37]. Furthermore, Fraga et al. [38] observed that mean phenological timings were projected to undergo significant advancements in Portugal (e.g., budburst and harvest may be anticipated by 1 month or more), with implications also in the corresponding pheno-phase intervals. Impacts of climate change on viticulture likely will be significant [39] in terms of increasing temperature combined with extreme events such as droughts and/or short-term storms [40]. Hence, there is the need for accurate tools for estimating crop water and irrigation requirements, namely soil water balance (SWB) models using weather data and soil water observations and/or other data such as eddy covariance and sap flow data [41,42]. In addition to local ground observations, also remote sensing data may be used in models. Pôças et al. [43] used hyperspectral reflectance data derived from a handheld spectroradiometer to estimate the predawn leaf water potential to assess the water status of grapevine cultivars in the Port wine region to be used in irrigation scheduling. A two-source model was used by Ortega-Farias et al. [44] to estimate vineyard energy balance using thermal images acquired by an unmanned aerial vehicle (UAV). Romero et al. [45] also reported on using multispectral imagery from an UAV platform for water status estimation.
Through modelling, farmers and/or farm advisers use models to support irrigation scheduling under diverse climates and management scenarios. Models may also be used to assess the impact of management decisions on soil and water resources, e.g., SIMDualKc model, which is able to both the partition of crop evapotranspiration (ET) into transpiration and soil evaporation, and the partition of transpiration into the fractions relative to the crop and the understory vegetation [46,47]. Another example is given by Cellete et al. [48] that developed the WaLIS model to simulate water resources partitioning and to estimate ET and the water use for both grapevine and the inter-row cover crop. Phogat et al. [49] applied a Richards-based mechanistic model for computing soil water fluxes and improving agricultural practices in irrigated vineyards, and Kustas et al. [50] and Kool et al. [51] applied the thermal-based two-source energy balance model for monitoring daily ET in vines and the inter-rows. A biosensing IoT platform for water management in vineyards is referred by Loddo et al. [52]. Moreover, there is abundant bibliography for ET studies with grapevines as reviewed by Rallo et al. [53], but not with a cover crop in the inter-row irrigated or rainfed.
The FAO 56 dual-K c approach is a widely used method that refers to the determination of crop evapotranspiration (ET c ) as the product of a K c value for a specific crop stage and the grass reference evapotranspiration (ET o ) computed with the FAO Penman-Monteith equation [54,55]. The K c value is partitioned into the basal crop coefficient (K cb ) referring to transpiration and the soil evaporation coefficient (K e ), providing thus separate estimates of the ET c components: crop transpiration (T c ) and soil evaporation (E s ). Examples of applications to vine and fruit crop systems can be found in the recent literature review provided by Rallo et al. [53] and Pereira et al. [56]. The SIMDualKc model [57] has adopted the dual-K c approach for estimating daily ET fluxes of crops grown under different climatic regions and management practices, including grapevine [46,47,58], peach [59], and olive [60,61]. Yet, despite the growing use of the FAO56 dual-K c approach for computing crop evapotranspiration fluxes and improving soil water management in complex agricultural systems, there is still the need for extending research focusing on the search for standard K c and K cb values as well the impact of active ground cover, cover crops, and mulches on crop evapotranspiration of vine systems [53].
Considering the insufficient knowledge on vineyards evapotranspiration as related to cover crops, the current study aims to estimate evapotranspiration fluxes and water use in two Mediterranean vineyards, a northern Italy sloping rainfed crop planted in a soil with large soil water holding capacity, and an irrigated one in southern Portugal cropped in a flat area with a sandy soil. Differences between vineyards allow perceiving differences in management and in water responses to both the rainfall regime and the ET. Objectives consist of (i) calibration and validation of the water balance model SIMDualKc [57], already proved for vine and tree crops; (ii) evaluation of the influence of the inter-row cover crop on the evapotranspiration dynamics and the soil water balance of the rainfed vineyard (2016-2019); (iii) evaluation of the evapotranspiration dynamics, the soil water balance, and the irrigation management applied to the drip irrigated vineyard (2018-2020); (iv) with support of the referred model, demonstrating the importance of improved soil and water conservation measures and practices for both Italian and Portuguese cases. The goal of the study focusses on the sustainable use of soil and water resources, thus on the sustainability and resilience of the vineyards production systems taking in account the challenges of climate change.

Description of the Study Areas
The Italian case study was located at the "Tenuta Cannona" Experimental Vine and Wine Center of Agrion Foundation (44 • 40 N, 8 • 37 E, 296 m a.s.l.), in the municipality of Carpeneto (AL), in the southern part of the Monferrato hilly area, known as "Alto Monferrato", North-West Italy. Data was collected from January 2016 to December 2019. The climate is alpine sublitoranean. According to records from the nearest weather station over the period 1951-1990 (Ovada, 187 m a.s.l.), the average annual precipitation is 965 mm, mainly concentrated in Autumn (October and November) and Spring (March), while the driest month is July [62]. At the experimental site, over the period 2000-2019, the average annual precipitation was slightly lower (881 mm), ranging from a maximum of 1455 mm (2019) to a minimum of 493 mm (2017). The mean annual air temperature was 13 • C ( Figure 1). The Cannona vineyards lie on Pleistocenic fluvial terraces in the Tertiary Piedmont Basin, including highly altered gravel, sand, and silty clay deposits, with red alteration products [63]. The main physical characteristics of the studied soils are presented in Table 1. The soils had clay to clay-loam texture and were classified as Dystric Cambisols [64]. Particle size distribution, soil bulk density (ρ b ), and soil water contents at saturation (θ s ) and at field capacity (θ FC ) were obtained from undisturbed soil samples taken at different depths according to Blake and Hartge [65] and Cavazza [66], whereas soil water contents at the wilting point (θ WP ) were obtained using the Rosetta pedotransfer functions [67] and the particle size distribution and ρ b as input ( Table 1).
The Portuguese case study was located at Companhia das Lezírias, Samora Correia, southern Portugal (38.808 • N, 8.900 • W, 45 m a.s.l.). Data were collected from January 2018 to October 2020. The climate in the region is dry sub-humid, with mild winters and hot, dry summers. The mean annual precipitation is 669 mm, mainly concentrated between October and May, while the mean annual temperature is 16.8 • C. The weather data for the study area was obtained from the local weather station and is reported in Figure 2.  Weather data for Monferrato plots during the study period (P, precipitation, ETo, reference evapotranspiration; Tmax and Tmin, maximum and minimum air temperatures, RHavg, mean relative humidity; Rs, solar radiation, u2, wind speed at 2 m height).   Weather data for Monferrato plots during the study period (P, precipitation, ET o , reference evapotranspiration; T max and T min , maximum and minimum air temperatures, RH avg , mean relative humidity; R s , solar radiation, u 2 , wind speed at 2 m height).   The soil was classified as a Haplic Fluvisol [64], with main physical characteristics presented in Table 1. The particle size distribution was obtained using the pipette method for particles having diameters <2 µm (clay fraction) and between 2-20 µm (silt), and by sieving for particles between 20-200 µm (fine sand) and between 200-2000 µm (coarse sand). These textural classes follow the Portuguese classification system [68] and are based on international soil particle limits (Atterberg scale). The dry bulk density (ρ b ) was obtained by drying volumetric soil samples (100 cm 3 ) at 105 • C for 48 h. The soil water contents at saturation (θ s ), field capacity (θ FC ), and the wilting point (θ WP ) were obtained from pedotransfer functions using the particle size distribution as input [69,70].

Vineyards Management
As typical in the Monferrato area (Italy), vineyards were rainfed. The experiment was carried out in a vineyard planted in 1988 with Barbera vines, managed according to conventional farming for wine production. Two vineyard plots of 1221 m 2 (16.5 m wide and 74 m long) each were considered. Each plot was composed of 6 rows aligned along the slope (SE aspect, average slope 15%), spaced 2.75 m, where the vines were spaced 1.0 m along the row and grown on vertical shoot positioned trellis (VSP). Since 2000, the soil in the two plots has been managed with different techniques: (i) conventional tillage (CT, hereafter) cultivation with chisel (at a depth of about 0.25 m); and (ii) controlled grass cover (GC), i.e., mulching of the spontaneous grass cover. Both practices were usually carried out twice a year, in spring and autumn. Then, the grass cover in the inter-row was mowed additional times during the season, if necessary. Weeds under the rows were controlled with Glyphosate application in spring, on the surface, 0.6 m across the vine row. Wood pruning was carried out in winter, and residues were chipped in one inter-row out of two. Most of the farming operations in the vineyard were carried out using tracked or tyre tractors carrying or towing implements, with intensification of passages from spring to grape harvest time (from 14 to 27 passages per year). The dates of the main crop stages during the four growing seasons are reported in Table 2. Note: GDD is the cumulative growing degree-days.
In the 4-years study period, i.e., from January 2016 to December 2019, weather variables, runoff amounts, soil losses, and hourly soil water content related to 89 runoff events were recorded for the two plots. Rainfall, air temperature, and humidity were recorded at 10-min intervals by a weather station placed near the plots. Daily values of solar radiation and wind speed and direction were obtained from stations located 10 km from the vineyard (Acqui Terme, Basaluzzo), belonging to the Regional Environmental Agency [71]. Each vineyard plot was hydraulically "isolated". The runoff generated by rainfall was collected separately for each plot by a channel connected with a sedimentation trap. A tipping bucket device measured the hourly volumes of runoff (RO, mm) in both CT e GC. Runoff samples were collected to obtain sediment concentration for erosive events, and, if sedimentation occurred in the channels and sediment trap, then the sediment yield was collected and weighted (see Biddoccu et al. [72] for details). Soil water contents were recorded every hour from an average of 1-min by indirect method [73] measurements of capacitance/frequency domain sensors (ECH 2 O-5TM sensors, Decagon Devices Inc., Pullman, WA, USA), gravimetrically calibrated, placed at 0.1, 0.2, and 0.3 m depth, and stored by a Decagon EM50 datalogger. Soil water content measurements were carried out in the two plots, both in the track position (T), which is the portion of inter-row affected by the passage of tractor wheels or tracks, where the compressive effects tend to concentrate [74], and in the middle of the inter-row, identified as the no-track position (NT), which is not affected by direct contact with tractor wheels or tracks, and then averaged for comparison with model output.
The selected field in Samora Correia (Portugal), planted in 2008, was relatively flat, and part of a larger vineyard (130 ha). The field was a drip-irrigated plot, 5 ha in size, planted with different varieties of wine grapes but where Touriga Nacional was dominant. The plants were grown on VSP trellis, with wood pruning during the dormant period in winter. Plants were at a row distance of 1.0 m and a row spacing of 2.8 m, thus a plant density of approximately 3571 plants ha −1 , with an orientation in the east-west direction. The dates of the main crop stages during the three growing seasons can also be found in Table 2. Harvest of grapes in Samora Correia vineyards was performed during September. The inter-row was covered with spontaneous grass from autumn to spring. When necessary, weeds under the rows were controlled with Glyphosate application in spring. Irrigation was delivered through a drip system, with management practices performed according to the standard practices in the region and decided by the farmer. Drippers were spaced 1 m apart, and the drip line was placed on the trellis 0.5 m above the soil surface. The total water applied through irrigation summed 470, 625, and 465 mm in 2018, 2019, and 2020, respectively. The application depth during irrigation events varied from 1 to 12 mm. Soil water contents were continuously monitored in two locations at depths of 10, 20, 30, 40, 50, 60, 70, and 80 cm using EnviroPro MT (MAIT Industries, Australia) capacitance probes. Probes were installed in the crop row, distanced approximately 30 cm from emitters, and measurements were averaged for comparison with model output. The SIMDualKc model [57] computes the daily soil water balance at the field scale as follows: where D r is the root zone depletion (mm), P is the rainfall (mm), RO is the runoff (mm), I is the net irrigation depth (mm), CR is the capillary rise from the groundwater table (mm), DP is the deep percolation (mm), and ET c act is the actual crop evapotranspiration (mm), all referring to day i or i−1. In this study, CR was not considered as the groundwater table was too deep (>5 m) in both case studies and could not contribute to crop evapotranspiration. The SIMDualKc model adopts the FAO56 dual K c approach for computing crop evapotranspiration [42,54,55]. In this approach, the components relative to crop transpiration (T c , mm) and soil evaporation (E s , mm) are computed separately as follows: where K cb (-) is the standard basal crop coefficient that refers primarily to crop transpiration although some diffusive soil evaporation may also be included, particularly during the initial crop stage, K e is the evaporation coefficient (-) that describes direct evaporation from the surface soil layer of depth Z e (cm), and ET o is the reference evapotranspiration (mm) computed with the FAO Penman-Monteith equation [54]. T c values are reduced when water stress occurs through a multiplier stress coefficient (K s ) to K cb : where T c act is the actual crop transpiration (mm) and K cb act is the actual basal crop coefficient (-). When only the matric potential is constraining T c values, K s is computed as follows [42,54]: where TAW and RAW are, respectively, the total and readily available soil water relative to the rooting depth (mm). These are computed as: where θ FC and θ WP are the soil water contents at field capacity and the wilting point (m 3 m −3 ), respectively, Z r is the root depth (m), and p (-) is the soil water depletion fraction for no stress. When the depletion exceeds p, i.e., the soil water content drops below RAW, T c values are reduced due to water stress (K s < 1.0), otherwise K s = 1.0. Soil evaporation is limited by the amount of energy available at the soil surface in conjunction with the energy consumed by transpiration, and by water availability in the surface soil evaporation layer [54,55]. The evaporation coefficient (K e ) is then computed as: (8) where K r is the evaporation reduction coefficient (0-1), K c max is the maximum value of K c (-) (i.e., K cb + K e ) following a rain or an irrigation event (−), and f ew is the fraction of the soil that is both exposed to radiation and wetted by rain or irrigation. f ew depends upon the effective fraction of ground covered or shaded by vegetation near solar noon (f c eff ). K r is calculated using the two-stage drying cycle approach where the first stage is the energy limited stage, and the second is the water limited stage [42,54,55,75]: where TEW is the maximum depth of water that can be evaporated from the evaporation soil layer when it has been completely wetted (mm), REW is the depth of water that can be easily evaporated without water availability restrictions (mm), and D e is the evaporation layer depletion at the end of day i−1 (mm). D e is computed through a daily water balance of the evaporation soil layer, with the evaporation decreasing as the evaporable soil water decreases in the evaporation soil layer beyond REW. DP is calculated using a time decay function relating the soil water storage near saturation with the time after the occurrence of heavy rain or irrigation [76]: where W a is the actual soil water storage in the root zone (mm), a D is the soil water storage comprised between saturation (θ S ) and θ FC (mm), b D is an empirical dimensionless parameter (-), and t is the time after irrigation or rain that produces storage above field capacity (days). Further information on the estimation of a D and b D is provided by Liu et al. [76]. RO is estimated using the curve number (CN) approach [77], with the CN value depending on the soil type, vegetation type, and antecedent soil moisture. In SIMDualKc, CN is adjusted each day to reflect the impact of increasing or decreasing of the soil water content on soil infiltration properties by relating the CN value to the soil water depletion in the surface layer (D e ). A different CN value is thus computed according to the correspon-dence between D e and the antecedent soil water conditions AWC I, AWC II, and AWC III (respectively, representing dry, average, and wet soil conditions) [77,78]: with: (16) and: where CN I , CN II , and CN III are, respectively, the curve numbers associated with the antecedent soil water conditions AWC I (dry), AWC II (average), and AWC III (wet) (0-100), and D e-AWC I , D e-AWC II , and D e-AWC III are, respectively, the depletion of the evaporative layer at AWC I, AWC II, and AWC III conditions (mm). CN II corresponds to tabulated values available in USDA-SCS [77] and Allen et al. [78]. The K cb values are described for four growth stages: the initial stage or start of the crop season, the rapid growth or development stage, the mid-season stage, and the late-season stage. Yet, the K cb curve is defined by three K cb values corresponding to the initial (K cb ini ), mid-(K cb mid ), and end-season (K cb end ). For tree crops and vineyards, these are computed based on crop characteristics [42,53,79,80]. When the inter-row has an active ground cover, the K cb relative to the crop and the understory vegetation may be estimated with the following Equation [57,79]: where K cb gcover is the K cb of the ground cover vegetation in the absence of tree foliage, K cb full is the estimated basal K c during peak plant growth for conditions having nearly full ground cover, and K d is the crop density coefficient. The second term of the max function reduces the estimate for K cb during the mid-season stage by half the difference between K cb full and K cb cover when this difference is negative. This accounts for impacts of the shading of the surface cover by overstory vegetation having a K cb that is lower than that of the ground cover due to differences in stomatal conductance. When no ground cover exists or when the cover crop dries out becoming a less dense residuals mulch, the previous equation 19 is simplified by replacing K cb cover with the minimum K c for bare soil (K c min = 0.15). The K cb full is estimated primarily as a function of crop height and then adjusted for tree crops using a reduction factor (F r ; from Pereira et al. [56]) estimated from the mean leaf stomatal resistance, as follows [80]: where u 2 is the average daily wind speed (m s −1 ) at a height of 2 m above ground level during the crop growth period, RH min (%) is the average daily minimum relative humidity during the growth period, and h is the mean plant height (m) during the mid-season.
Before climatic adjustment, an upper limit for K cb full is assumed 1.20. The effect of the crop height is considered through the sum (1 + k h h), with k h = 0.1 for tree and vine crops [80]. Higher K cb full values are expected for taller crops and when the local climate is drier or windier than the standard climate conditions (RH min = 45% and u 2 = 2 m s −1 ). When the vegetation shows more stomatal adjustment upon transpiration, the parameter F r requires an empirical adjustment (F r < 1.0), otherwise F r = 1.0. For trees and vines, F r is closer to 1.0 when crops exhibit great vegetative vigor and decreases under a limited water supply and due to pruning and training [80]. The density coefficient (K d ) is estimated from the fraction of ground cover as follows [79]: (21) where f c eff is the effective fraction of ground covered or shaded by vegetation near solar noon (-), M L is a multiplier on f c eff (1.5-2.0) describing the effect of the canopy density on shading and on maximum relative evapotranspiration per fraction of ground shaded (to simulate the physical limits imposed on water flux through the plant root, stem, and leaf systems), and h is the mean height of vegetation (m).
The K cb gcover is computed through the following set of equations [57]: f c eff gcover = density. f c gcover (24) K d gcover = min 1, M L f c eff gcover , f ( 1 1+hgcover ) c eff gcover (25) K cb gcover = K c min + K d gcover K cb gcover full − K c min (26) where the variables shown are the same as defined earlier but refer to the active ground cover (subscript gcover). The equations 19 through 26 are not used when there is no need for partition of transpiration between the crop and the active ground cover.

Model Setup
The input data required by the SIMDualKc model were the following: • Climatic data: daily weather data from local meteorological stations, namely the maximum and minimum temperatures (T min and T max , • C), minimum and maximum relative humidity (RH max and RH min , %), global solar radiation (R s , MJ m −2 day −1 ), wind speed at 2 m height (u 2 , m s −1 ), and rainfall (P, mm) (Figures 1 and 2). • Soil data: θ S , θ FC , and θ WP as well as the particle size distribution of the different soil layers defined in the rootzone domain (Table 1). • Soil evaporation parameters: the evaporable layer depth (Z e ), with TEW and REW being then estimated using the textural and water holding characteristics of the Z e [54,55]. • Initial conditions: initial (observed) values of the soil water content in both the root zone (% of TAW) and the evaporation layer (% of TEW) ( Table 3).
• Crop data: the dates defining the different stages of vineyards during growing seasons ( Table 2) and respective values of K cb , p, h, the fraction of ground covered (f c ), M L , and Z r . Values of h and f c were observed in the field and are presented in Table 4. Z r was set to 1.0 m in both case studies. The M L value is unique for all crop stages and was set in both case studies to 1.5 following Pereira et al. [56]. • Active ground cover: the density of the active ground cover (28-30% and 15-20% for the GC and CT plots in Italy, respectively, and 15-20% for the Portuguese plot); the initial values of active ground cover fraction (f c gcover ) and height (h gcover ), which were then updated for each crop stage according to Table 4; and M L , which was also assumed to be 1.5 following Fandiño et al. [58]. • Deep percolation: the deep percolation parameters a D and b D relative to the parametric equation of Liu et al. [76] were defined according to the soil texture data, and θ S and θ FC values in the soil profiles (Table 1). • Runoff (more relevant for the Italian case studies): the CN value for each inter-row condition [77]. • Irrigation: the dates of irrigation events and irrigation depths, which were specified according to observations; and the fraction of the soil surface wetted by irrigation (f w ), which was measured in the field and set to 0.13. Note: TAW, total available water; TEW, total evaporable water.

Model Calibration and Validation
The SIMDualKc model was calibrated by adjusting model parameters one at a time within reasonable ranges of values, using a trial-and-error procedure, until deviations between measured soil water contents and model predictions were minimized. For Monferrato plots, data from the 2019 growing season were used for calibration while for the Samora Correia plot the 2020 growing season data were used. The calibrated parameters were (i) the CN value (more relevant for the Italian plots); (ii) the K cb and p values relative to the initial, mid-, and end-season; (iii) the depth of the evaporative soil layer (Z e ) and the TEW and REW values relative to soil evaporation; and (iv) the parameters relative to the deep percolation function, a D and b D .
The calibrated model parameters were validated using independent data sets corresponding to soil water contents measured in the Italian plots during the 2016-2018 growing seasons, and relative to the growing seasons 2018-2019 measured at Samora Correia. The resulting K cb values are adjusted for climate. Following Pereira et al. [81], model performance was considered acceptable when the goodness-of-fit indicators relative to the validation were within a range of "20% variation" relative to calibration.
The goodness-of-fit indicators adopted for comparing SIMDualKc model simulations with soil water content observations were the coefficient of determination (R 2 ) of the ordinary least squares regression, the coefficient of regression (b 0 ) of the linear regression forced through the origin, the root mean square error (RMSE), the ratio of the RMSE to the standard deviation of observed data (NRMSE), the percent bias of estimation (PBIAS); and the modeling efficiency (EF), respectively, given as: where O i and P i are, respectively, the observed and model predicted values at time i, O and P are the respective mean values, and n is the number of observations. R 2 values close to 1 indicate that the model explains well the variance of observations. b 0 target is 1.0 with values b 0 < 1.0 indicating under-estimation of the predicted values and b 0 < 1.0 indicating over-estimation. RMSE and NRMSE close to zero indicate small estimation errors and good model predictions [82]. PBIAS values close to zero indicate that model simulations are accurate, while positive or negative values indicate under-or over-estimation bias, respectively. EF values close to one indicate that the residuals' variance is much smaller than the observed data variance, hence the model predictions are good. On the contrary, when EF is close to zero or negative, there is no gain in using the model [83]. Table 4. Fraction of ground cover f c , crop height (h), fraction of the active ground cover (f c gcover ) and height of the active ground cover (h gcover ) in every plot and all growing seasons.

Monferrato (CT) Monferrato (GC) Samora Correia Vineyard
Inter-Row Vineyard  Table 5 presents the model parameters calibrated for the Italian and Portuguese plots with data relative to 2019 and 2020, respectively. These parameters were then validated for 2016-2018 relative to the Monferrato vineyards and for 2017 and 2019 relative to the Samora Correia vineyard. The initial K cb values were obtained from tabulated data provided by Pereira et al. [56]. For Monferrato vineyards in year 2019, the calibrated K cb ini , K cb mid , and K cb end values were, respectively, 0.20, 0.47, and 0.34 in the conventional tillage CT plot, and, respectively, 0.35, 0.47, and 0.40 in the ground cover GC plot. The active ground cover promoted an increase of the K cb values in the initial and end season, while no differences were noticed in the mid-season. During this period, the inter-row grass was at the end of its annual cycle or already dried out. In the Samora Correia plot in year 2020, the calibrated K cb ini , K cb mid , and K cb end were 0.17, 0.47, and 0.39, with the same K cb mid as in Monferrato but with different K cb ini and K cb end . Note: CT, conventional tillage; GC, grass cover; K cb , basal crop coefficient for the initial (K cb ini ), mid (K cb mid ), and end season (K cb end ); K cb nongrowing , basal crop coefficient during the non-growing period; p, depletion fraction for no stress during the initial (p ini ), mid (p mid ), and end season (p end ); TEW, total evaporable water; REW, readily evaporable water; Z e , depth of the soil evaporation layer; a D and b D , parameters of the deep percolation; CN, curve number.

Model Parametrization
Vineyards are heterogeneous, sparsely vegetated surfaces with complex canopies. For that, the K cb values may vary among locations due to differences in varieties, age, training, irrigation, soil cover, and crop management as demonstrated in the recent literature review by Rallo et al. [53]. These authors updated the K cb values proposed by Allen et al. [54] for vines and other perennial crops. Rallo et al. [53] provided standard K cb for diverse training systems used in vineyards combined with the fraction of the ground cover (f c ) but not including active ground cover calculations. After reviewing literature relative to wine grapes, it was concluded that for the training system adopted in the all the studied vineyards (vertical shoot positioned trellis), the f c values were reported to vary from 0.25 to 0.45, with K cb mid values ranging from 0.46 to 0.80 and K cb end values varying from 0.20 to 0.60. The K cb values in both Monferrato and Samora Correia vineyards were therefore close to the lower limit of the range of variation of VSP trained vineyards reported by Rallo et al. [53].  Table 4) were closer to the lower f c values reported by Rallo et al. [53], hence indicating that it is likely that the low K cb mid and K cb end values found in the current study are due to a small f c observed in the vineyards under study, so in agreement with reported data of those authors.
The calibrated K cb values listed in Table 5 correspond to the entire vineyard system (vine plants + active ground cover). Table 6 shows the partitions of those values based on formulations given in Section 2.3 and crop density data (f c and h) given in Table 4 for each crop stage and growing season. Contribution of the active ground cover to the K cb value was obviously higher in the Monferrato GC plot, but the K cb values of vine plants were smaller, especially during the mid and end season stages since, as explained earlier, the cover crop may influence plant's vigor and reduce canopy leaf area and T c act rates of the grapevine in the earlier development stages which will then also affect the mid-season period. When looking at the partition of the K cb values, it becomes also the higher K cb mid value at Samora Correia (Portugal), resulting from higher air temperature registered in that location. The calibrated p ini , p mid , and p end values (Table 5) matched those proposed by Allen et al. [54] in the Monferrato CT plot but were slightly smaller in the Monferrato GC and the Samora Correia vineyard. The GC plot showed the lowest p values, likely due to the competition for water by the active ground cover.
The Z e , TEW, and REW, as well as the a D and b D calibrated values reflect the hydraulic properties of the soils of the studied vineyards, a loamy sand to sandy textured soil in Samora Correia vs. a clay to clay-loam soils at Monferrato. Lastly, the CN values were set to 70 and 55 in Monferrato CT and GC plots, respectively, with the inter-row grass cover clearly impacting the calibrated values of CN. Those values are lower than the CN reported by Gaudin et al. [84] from an analysis of three runoff events (CN = 82) for a bare soil. For olive orchards, Romero et al [85] reported CN = 89 for bare soil, CN = 82 and 70 for well-established cover crop strips with, respectively, 1 and 3 m wide, and CN = 88 and 87 for degraded cover crop strips with the same widths. For the Samora Correia vineyard, where runoff is less relevant due to the flat topography, the CN = 65 value is, as expected, smaller. Romero et al. [85] pointed out the need for improvements in the SCS-CN method aiming at a better understanding of the effects of soil management practices on runoff formation. As shown in Biddoccu et al. [7], different soil management implementations for the same treatment can lead to high variability in the soil roughness and coverage (namely the RUSLE C-factor), for the same region or among different countries.  Figure 4 shows the measured and simulated SWC in the root zone of the Samora Correia vineyard during 2018-2020. Both figures further provide the dates and depths of rainfall events and, for the Portuguese vineyard, also of the irrigation dates and depths. Both figures show that the measured values of SWC were kept between θ FC and θ p during most of the year. In Monferrato, this fact resulted from precipitation, with large or successive events occurring during the autumn and spring, leading to SWC values temporarily above θ FC . However, in the summer dry season, the soil moisture dropped below θ p for some extended periods when rainfall was lacking. The year 2017 had the lowest SWC values as annual precipitation was uncharacteristically low during that year, amounting to only 493 mm. Indeed, among the last 20 years, 2017 was the least rainy one, with only 56% of the mean annual precipitation, while summer was very dry, with less than 25 mm of precipitation from the beginning of June to the end of August (the average for this period in 20 years is about 100 mm).      In Samora Correia (Figure 4), drip irrigation was commonly applied at small depths (up to 12 mm) aimed at maintaining SWC above θp during the summer dry season. This type of irrigation management is uncommon in Portuguese vineyards of southern Portugal, where deficit irrigation is usually practiced during the summer dry period to better control shoot vigor and ripening, as well as the quality of fruit attributes [86,87]. Thus, the case observed at Samora Correia is rare and likely results from the very low water holding capacity of the soil, which calls the growers to avoid SWC < θp, at least during somewhat large periods.

Model Performance
The statistical indicators used to evaluate the agreement between measured and simulated SWC values are presented in Table 7    Simulated SWC were higher in GC plot relative to CT for all years, in agreement with the average of measured values. Likely, this would result from higher infiltration of rainfall water as already observed [5,11,12]. Nevertheless, Novara et al. [14], relative to a review of studies about the effects of a cover crop on SWC, reported that such effects often result in a decrease of soil moisture at a soil depth from 0-0.5 to 0-1.0 m, i.e., below the roots of the grass crop, between bloom and veraison, which can threaten yields in semiarid conditions.
The statistical indicators used to evaluate the agreement between measured and simulated SWC values are presented in Table 7 for both case studies. The SIMDualKc model performed well when simulating SWC in the Monferrato plots for calibration and for validation. The regression coefficient b 0 varied from 0.95 to 1.06 in the CT plot and from 0.91 to 1.03 in the GC plot, i.e., close to the 1.0 target indicating that simulated values are close to the observed ones. The value of R 2 was relatively high in both plots (0.76 to 0. 90), showing that the model could explain most of the variability of the observed data. The errors of the estimates were small, resulting in a RMSE value ranging from 0.01 to 0.03 m 3 m −3 and a NRMSE value ranging between 0.32 and 0.81. In agreement with b 0 , the PBIAS were small, not expressing a particular over-or under-estimation trend in simulating the measured data. The EF values were higher for CT than for GC (respectively, 0.53 to 0.90 and 0.34 to 0.80), indicating that the variance of the residuals was smaller than the measured data variance particularly for the calibration set. The goodness-of-fit indicators show better results for the CT plot relative to the GC case. This indicate that a supplemental effort for characterizing the active ground cover of the GC plots is required. Nevertheless, results obtained are reasonable enough for the pretended analysis. Table 7. Goodness-of-fit indicators for the adjustment between measured and simulated values in all case studies.

Site
Year The performance of the SIMDualKc model during the calibration period (2018) for the Samora Correia vineyard was comparable to results in the Italian case studies, with the PBIAS values, close to zero, also indicating accuracy in predicting the measured SWC. The regression coefficient b 0 varied 1.00 to 1.02, thus leading to assume that there is no noticeable trends in the simulation, which is confirmed by the data. However, both the R 2 and EF values are not high, indicating that improvements in describing vineyards crop data are required, particularly relative to the ground cover that changes from active to dry residues. The errors of the estimate RMSE were small, but the NRMSE are high because soil water depletion by evapotranspiration are low.
Issues related to the quality of the measured dataset and the empirical one-dimensional modeling approach used here to describe the soil water balance in a three-dimensional drip irrigation system may also help explaining why the goodness-of-fit indicators are not as high as in former studies, e.g., Paço et al. [88]. Nevertheless, the reported goodness-of-fit indicators are within the range of values reported in the literature for soil water content simulations of field crops [89,90], vegetable crops [91,92], and perennial crops [58][59][60][61] using the same model.

Soil Water Balance in the Hillslope Vineyards of Monferrato, Northern Italy
In Figure 5 are presented the daily values of the potential non-stressed basal crop coefficients (K cb ), the actual basal crop coefficients (K cb act ), the soil evaporation coefficients (K e ), and the actual crop coefficients (K c act = K cb act + K e ) computed by SIMDualKc. The rainfall events are also depicted in Figure 5. Table 8 presents then the components of the  soil water balance in the Monferrato plots from 2016 to 2019.   2016  CT  0  779  −156  372  262  234  27  90  GC  0  779  −85  417  298  234  118  43  2017  CT  0  493  35  395  275  206  42  8  GC  0  493  −17  445  300  203  8  2  2018  CT  0  1190  −100  347  313  307  310  178  GC  0  1189  −86  386  326  306  315  92  2019  CT  0  1476  −68  353  313  274  467  355  GC  0  1476  −40  403  339  284  600  198 Note: I, irrigation; P, precipitation; ΔSW, soil water storage variation; Tc, potential transpiration; Tc act, actual transpiration; Es, evaporation; DP, deep percolation; RO, runoff. The K cb act matched the K cb values during most of crop season duration except for the dry summers ( Figure 5). The resulting seasonal T c act values were estimated to vary from 262 to 313 mm in the CT plot and 298 to 339 mm in the GC plot (Table 8). T c act values were up to 14% higher in the GC plot because comprising the transpiration by the grass cover in the inter-row. The T c act reduction due to water stress started during mid-season and extended to the late-season due to low soil water availability during these dry months. The T c act /T c ratios are always smaller during the late season (Table 9) and the lowest value was observed in the late season of 2017 (0.30) due to the exceptional lack of rainfall in this year. Model results even suggested the need for irrigation in such dry conditions despite northern Italy is a traditionally rainfed wine-growing region. The impact of water stress on transpiration rates as given by the T c act /T c ratio was quite similar in both CT and GC plots with a general ratio value slightly smaller in the GC plot over three out of the four seasons. Despite the computed difference is small, it is noticeable in both the mid-and late-seasons, which explains the reluctance of winegrowers in using active ground cover in the inter-row due to concerns over soil water competition between the vines and the grass cover.  Note: T c , potential transpiration; T c act , actual transpiration; T c act gcover , actual transpiration of the active ground cover; E s , evaporation. Table 9 further shows the T c act values relative to the active ground cover (T c act gcover ). In the GC plot, the weight of T c act gcover relative to T c act of the entire vineyard system was always high, amounting to a maximum of 63.6% to 77.3% of the T c act during the non-growing periods and to a minimum of 17.2% to 40.7% in the mid-season. In the CT plot, the ratios T c act gcover /T c act were naturally always smaller than those reported for the GC plot. It should be noticed that the f c values for the vineyard in the GC plot were always slightly lower than those in the CT plot (Table 4), which further exacerbates the importance of the cover crop to the soil water balance and preservation of soil and water resources in these hillslope areas.
The K e values reached their maximum with the increase of soil moisture in the soil surface layer after rainfall events, rapidly decreasing in the following days as the soil dried out ( Figure 5). For this reason, the K e values have an enormous variability throughout the time. E s values were relatively similar in the GC and CT plots, with the small differences resulting from the different Z e set during calibration. The number of rainfall events influencing soil evaporation in each growing season is highly varied. E s values in the CT plot varied from 206 to 307 mm while in the GC plot they varied from 203 to 306 mm. Soil evaporation occurred mostly during the non-growing periods and early crop stages (Table 9) and was slightly larger in the CT plot because a fraction of the potential E s was used for ground cover transpiration.
In the CT plot, the seasonal simulated surface runoff ranged from 8 mm in 2017 to 355 mm in 2019. These values contrast with those estimated for the GC plot, which varied from 2 to 198 mm, i.e., 44% to 75% lower than in the CT plot. As shown in Figure 6, which compares the simulated and measured seasonal surface runoff in both plots during 2016-2019, the model provided a narrower range of runoff values when compared to field measurements. However, differences between CT and GC were well noticed. Computing runoff using a single CN value led the model to underestimate measured values during rainy years and overestimate them during drier seasons. Similarly, Celette et al. [48], when modelling the water balance in vineyards, observed that the use of a constant CN value for the surface runoff could not consider changes in the soil surface and their effects on runoff and water infiltration. Nevertheless, the SIMDualKc was able to catch the seasonal trends as well as the effect of inter-row management on runoff since measured values in the GC plot were also 50% to 78% lower than in the CT plot, thus very similar to the trend in model predictions. The model results are also in line with results of previous studies about cover crop and water conservation in olive orchards and vineyards in the Mediterranean area as reviewed by Novara et al. [14]. These authors showed that, on average, the annual runoff coefficient was reduced by 27% with the application of cover crop management. However, its effect on water conservation was higher in rainier regions (rainfall higher than 700 mm), whereas in drier regions it was not relevant. In sloping vineyards, the grass cover results not only in a reduction of runoff, but it contributes to decreasing significantly soil losses, up to 74% during very intense and erosive events [10,93]. On the one hand, improving the description of the soil coverage, either with the CT or the GC conditions, should improve the input data to the model and, therefore, its simulation capability; on the other hand, it is likely that simulations capabilities of the model may be enhanced. Nevertheless, current results are satisfactory and support the interpretation of the differences in behavior of vineyards, namely aspects that relate with the conservation of soil and water resources.
Lastly, the estimated percolation summed 27 to 467 mm in the CT plot and 8 to 600 mm in the GC plot. Except for 2017, which results were uncommon due to the very low precipitation, percolation was always higher 1.6% to 337% in the GC plot, with the grassed inter-row promoting higher infiltration rates, less runoff and soil erosion, thus higher groundwater recharge. Through direct measurements in vineyards of Montpellier, also with Mediterranean climate, Gaudin et al. [84] also observed that the presence of a permanent cover crop had a positive effect on water infiltration, favoring the refilling of the soil water profile. In any case, they observed that the soil water refilling at budbreak time was subjected to the rainfall pattern during the fall season. ter 2022, 14, x FOR PEER REVIEW 23 of 31 Lastly, the estimated percolation summed 27 to 467 mm in the CT plot and 8 to 600 mm in the GC plot. Except for 2017, which results were uncommon due to the very low precipitation, percolation was always higher 1.6% to 337% in the GC plot, with the grassed inter-row promoting higher infiltration rates, less runoff and soil erosion, thus higher groundwater recharge. Through direct measurements in vineyards of Montpellier, also with Mediterranean climate, Gaudin et al. [84] also observed that the presence of a permanent cover crop had a positive effect on water infiltration, favoring the refilling of the soil water profile. In any case, they observed that the soil water refilling at budbreak time was subjected to the rainfall pattern during the fall season.  Table 10. The exception was only a small period during the rapid growth stage of the 2018 growing season where mild water stress was noticed. Irrigation was fundamental for such results, which otherwise would lead to a reduction of the Tc act values due to the high atmospheric demand observed particularly during the mid-and late season stages. The cover crop had only a small impact on Tc act values, with Tc act gcover being only expressive during the non-growing period (Table  11).

Soil Water Balance in the Irrigated Vineyard of Samora Correia, Southern Portugal
Like in the Italian case study, the Ke values reached their maximum after rainfall events, to then drop with the depletion of the evaporation soil layer. During the irrigation season, Ke values were maintained at a much lower range, close to 0.2 (Figure 7), resulting from the low fw value (0.13) describing evaporation from a small wetting bulb around emitters as opposite to rainfall events where the entire soil surface is wetted, and soil evaporation is maximized. The seasonal Es values then ranged from 235 to 281 mm (Table  10), corresponding to 40% to 45% of the ETc, and were distributed in relatively similar proportions throughout the different crop stages (Table 11).  Figure 7 presents the daily values of the K cb , K cb act , K e , and K c act computed with the SIMDualKc for the years of 2018-2020. Rainfall and irrigation events are also depicted. The K cb act matched the K cb values throughout the growing seasons, meaning that the T c act was always close to its potential values during the different stages of crop development as shown by the T c act /T c ratio always close to 1.0 in Table 10. The exception was only a small period during the rapid growth stage of the 2018 growing season where mild water stress was noticed. Irrigation was fundamental for such results, which otherwise would lead to a reduction of the T c act values due to the high atmospheric demand observed particularly during the mid-and late season stages. The cover crop had only a small impact on T c act values, with T c act gcover being only expressive during the non-growing period (Table 11).

Soil Water Balance in the Irrigated Vineyard of Samora Correia, Southern Portugal
Like in the Italian case study, the K e values reached their maximum after rainfall events, to then drop with the depletion of the evaporation soil layer. During the irrigation season, K e values were maintained at a much lower range, close to 0.2 (Figure 7), resulting from the low f w value (0.13) describing evaporation from a small wetting bulb around emitters as opposite to rainfall events where the entire soil surface is wetted, and soil evaporation is maximized. The seasonal E s values then ranged from 235 to 281 mm (Table 10), corresponding to 40% to 45% of the ET c , and were distributed in relatively similar proportions throughout the different crop stages (Table 11).
The computed percolation was high, with values varying from 265 to 342 mm in the considered years (Table 10) which relates with the high infiltration and hydraulic conductivity of the sandy soil. In 2020, 23% (65 mm) of the seasonal percolation resulted from irrigation excess, corresponding to 14% of the water applied. In contrast, in 2019, seasonal percolation resulted mostly from irrigation surplus (92%; 243 mm), mainly during the late season, which corresponded to 39% of the total water applied through irrigation. These findings evidence the need for improving irrigation management in the considered field. Therefore, an optimization of irrigation schedules was considered by setting irrigation depths to 5 mm and the management allowed depletion (MAD) equal to the depletion fraction for no stress (p). This application revealed a substantial reduction of percolation losses during the monitored seasons without affecting the remaining components of the soil water balance (Table 10). This is visible in Figure 8, which compares the daily values of the different components of the soil water balance monitored in the field with those computed by the model when optimizing the irrigation schedules. Nonetheless, irrigation water management in vineyards is relatively more complex than the analysis shown herein since impacts on fruit quality are not considered. Hence, in future studies, it would be of interest to relate the SIMDualKc simulated irrigation schedules with both grapes yield and main characteristics of fruit quality for a better understanding of the impacts of irrigation management on crop development and yields.

Conclusions
This paper presents two study cases, one at Monferrato, northern Italy, analyzing alternative soil management issues, conventional tillage (CT) and grass cover (GC) as active ground cover, the other at Samora Correia, southern Portugal, drip irrigated. The vineyards are cultivated in different climates, humid at Monferrato, where soils have a large water holding capacity, and dry subhumid at Samora Correia, where the soil is sandy and therefore has a small water holding capacity. The SIMDualKc model was successful used for both cases to simulate water use (crop evapotranspiration) and the soil water balance. The model was calibrated and validated using, respectively, 4 and 3 years of soil moisture measurements, thus making it able to simulate soil and water management conditions observed in northern Italy and southern Portugal case studies.
The SIMDualKc model, adopting the dual Kc approach, was able to partitioning crop evapotranspiration into transpiration and evaporation from the soil. Moreover, it was able to compute distinctively the transpiration from the vine plants and the active ground cover. Results are satisfactory but it was observed that a better description of both the active ground cover and the mulch residues resulting from drying the former in summer is required, as well as improving related model simulation.

Conclusions
This paper presents two study cases, one at Monferrato, northern Italy, analyzing alternative soil management issues, conventional tillage (CT) and grass cover (GC) as active ground cover, the other at Samora Correia, southern Portugal, drip irrigated. The vineyards are cultivated in different climates, humid at Monferrato, where soils have a large water holding capacity, and dry subhumid at Samora Correia, where the soil is sandy and therefore has a small water holding capacity. The SIMDualKc model was successful used for both cases to simulate water use (crop evapotranspiration) and the soil water balance. The model was calibrated and validated using, respectively, 4 and 3 years of soil moisture measurements, thus making it able to simulate soil and water management conditions observed in northern Italy and southern Portugal case studies.
The SIMDualKc model, adopting the dual Kc approach, was able to partitioning crop evapotranspiration into transpiration and evaporation from the soil. Moreover, it was able to compute distinctively the transpiration from the vine plants and the active ground cover. Results are satisfactory but it was observed that a better description of both the active ground cover and the mulch residues resulting from drying the former in summer is required, as well as improving related model simulation. The K cb values estimated over the different crop stages in the plot having the inter-row grass cover (GC) were higher than in the tilled CT plot, indicating higher transpiration of both the vines and the grass cover. The use of grass cover in the inter-row did not affect the ratio between the actual and potential root water uptake rates, i.e., the cover crop did not cause additional water stress to the grapevine except in conditions of great weather dryness of a peculiar year. On the contrary, the cover crop promoted infiltration of surface water and groundwater recharge, thus decreasing the risk of soil erosion and land degradation in addition to increase soil water availability, which confirmed previous studies reported in literature. Adopting the SIMDualKc model may support further analysis of problems and respective issues for improvement. Long term economic gains for adopting active ground cover are important relative to resource conservation, water and soil, due to controlling runoff and erosion with small increases in water stress. Nevertheless, financial gains in the short term are likely small, which limits active ground cover adoption by growers. Thus, it is adequate to search innovative approaches to adopt an active ground cover in the inter-rows, including searching the most adequate grass to be used, as well as with better consideration of the impacts of field slopes, always focusing both the soil and water conservation and the economic and financial return of farming.
Relative to the drip irrigation of sandy soils, the computed components of the soil water balance evidenced the need for improving irrigation management, mainly irrigation schedules, and the drip irrigation design, despite not dealt in the current study but known as limiting the performance of the irrigation system. In fact, crop transpiration needs to be increased relative to the total water use, which refers to increase the wetted bulbs to favor roots development and water extraction from the soil in combination with controlling percolation losses, so also favoring resource conservation. Using SIMDualKc as a decision support system should be appropriate.