Irrigation Management Scale and Water Application Method to Improve Yield and Water Productivity of Field-Grown Strawberries

: Improvements in water productivity are of primary importance for maintaining agricultural productivity and sustainability. Water potential-based irrigation management has proven e ﬀ ective for this purpose with many di ﬀ erent crops, including strawberries. However, problems related to spatial variability of soil properties and irrigation e ﬃ ciency were reported when applying this management method to strawberries in soils with rock fragments. In this study, a ﬁeld-scale experiment was performed to evaluate the impacts of three irrigation management scales and a pulsed water application method on strawberry yield and water productivity. An analytical solution to Richards’ equation was also used to establish critical soil water potentials for this crop and evaluate the e ﬀ ects of the variability in the soil properties. Results showed that spatial variability of soil properties at the experimental site was important but not enough to inﬂuence crop response to irrigation practices. The studied properties did not present any spatial structure that could allow establishing speciﬁc management zones. A four-fold reduction in the size of the irrigation management zones had no e ﬀ ect on yield and increased the water applications. Pulsed application led to signiﬁcant yield (22%) and water productivity (36%) increases compared with the standard water application method used by the producer at the experimental site.


Introduction
In the face of a growing world population and a decreasing availability of fresh water, it is crucial to rationalize agricultural water use [1]. Crop irrigation is estimated to account for 45% of freshwater uptake within the OECD countries, with values as high as 90% for individual countries [2]. Hence, any improvement in irrigation efficiency could greatly benefit the global balance of freshwater.
Strawberries (Fragaria × ananassa Duch.) are generally considered to have high water requirements [3]. In North American commercial field productions, day-neutral strawberries are generally grown in raised beds covered with polyethylene mulch and irrigated with a subsurface drip irrigation (SDI) system [4]. This system has a potential water application efficiency of 85% to 90%. However, the efficiency can be limited by the effects of: (1) Spatio-temporal variability in soil properties and plant water requirements [5][6][7][8], (2) moisture distribution patterns of the irrigation system [9], and (3) the use of inadequate irrigation management criteria [10].
The measurement of the soil water potential (SWP) with tensiometers has proven to be an effective irrigation initiation criterion that improves crop yield and water productivity (WP) for different crops under a wide range of soils, climatic conditions, and agricultural practices [11]. Studies have shown

Materials and Methods
An experiment was conducted on a commercial site located near Québec City, QC, CAN (46 • 54 N, 70 • 56 W, 29 m AMSL), during the 2012 growing season. According to the Canadian Soil Classification System [40], soils in the area are silt loams to silty clay loams from the Orléans series. The field had a constant 4% slope in the north-south direction and was divided into three sections of equal width in the east-west direction by natural wind barriers. Bare root day-neutral strawberry plants (cv. seascape) were planted on 20 May 2012 in raised beds 140 cm apart (center to center) with two plant rows per bed for a density of 54,800 plants per hectare. The beds were 90 cm wide, 30 cm high, and covered with black polyethylene mulch. The area of the experimental site was 2.67 ha, and the area of the individual plots ranged between 0.11 and 0.45 ha. Except for irrigation management, all agricultural operations (fertilization, pest and disease control, pruning of flowers and runners, and general maintenance) were conducted in a similar manner on all plots per the commercial producer standard procedures.
The experimental area was divided into 12 irrigation zones of equal width (nine beds or 12 m) but whose length varied between 180 and 90 m ( Figure 1A). In all zones, irrigation water was delivered by a single drip lateral per bed (1.6 cm diameter, 20 cm emitter spacing, and 0.5 L h −1 flow rate (Aquatraxx PBX, Toro-Irrigation, Bloomington, MN, USA)), buried 3-4 cm under the soil surface at the apex of the bed. For each individual zone, all the drip laterals were connected to a single tertiary water pipe. Upstream from the drip laterals, the water flow in each tertiary pipe was controlled by a manual gate valve (model IV 0818-10, Dubois Agrinovation, St-Rémi, QC, Canada), and the pressure head was adjusted to 62 kPa by a pressure-regulating valve (model PRV.75-112, Dubois Agrinovation, St-Rémi, QC, Canada). A water totalizer (model LR-F-208, Les Compteurs Lecompte ltd, Ste-Hyacinthe, QC, Canada) was also installed on the tertiary pipe upstream from the pressure-regulating valves ( Figure 1B). This setup allowed to individually control and measure the water flow to each individual zone while maintaining the pressure head in the drip laterals within the acceptable operating range for the drip emitters. Due to elevation differences within the zones, the pressure head in the drip laterals was expected to vary considerably, especially for longer zones. However, pressure variations as a function of the drip lateral length were calculated with the Hazen-Williams equation [41]. It showed that while the pressure gain was more important for longer zones, it did not significantly increase the variability of the emitter flow rate. As shown in Figure A1, the average and the coefficient of variation of emitter flow rate and the total flow rate per unit drip lateral length did not vary significantly for the two zone lengths considered in this study. For the remainder of this paper, a zone will be defined as a part of the field in which the water pressure and water flow can be controlled independently (as previously described), while a plot (discussed further in the following section) will be defined as a part of the field to which a specific irrigation management treatment was applied.
Four treatments were applied in a complete randomized block design with 3 replications. The blocks consisted of field sections separated with natural wind barriers. To evaluate the effect of irrigation management scale on yield and WP, the first three experimental treatments differed only by the size of the plot managed based on SWP measurements from a single monitoring station. The first treatment (Global) reproduced a large-scale management typically practiced by producers in the area where the experiment took place. For this treatment, a plot was 360 m long (two 180 m zones). The second treatment (Inter) aimed to reproduce an intermediate management. For this treatment, a plot consisted of a 180 m long irrigation zone. The third treatment (Local) reproduced a small-scale management. For this treatment, a plot consisted of a 90 m long zone. A fourth treatment (Pulse) was applied to evaluate the effect of water application method on yield and WP at the smallest management scale. For this treatment, a plot also consisted of a 90 m long zone.
significantly increase the variability of the emitter flow rate. As shown in Figure A1, the average and the coefficient of variation of emitter flow rate and the total flow rate per unit drip lateral length did not vary significantly for the two zone lengths considered in this study. For the remainder of this paper, a zone will be defined as a part of the field in which the water pressure and water flow can be controlled independently (as previously described), while a plot (discussed further in the following section) will be defined as a part of the field to which a specific irrigation management treatment was applied.  The irrigation management of each plot was based on SWP measurements. A monitoring station consisting of two buried wireless tensiometers (Hortau, QC, CAN) was installed in 12 of 15 zones ( Figure 1A). SWP was measured and transmitted online at five-minute intervals. A shallow probe, used to initiate irrigation, was installed at a depth of 15 cm at the outer limit of the wetting bulb from a drip emitter, i.e., a horizontal distance of 5 to 10 cm from the emitter. An additional probe was installed at a depth of 30 cm directly under the emitter and was used to adjust the irrigation duration to avoid leaching under the root zone. For all treatments, irrigation was initiated when a predefined irrigation threshold of −18 kPa was reached at a depth of 15 cm. This threshold was selected based on the results of previous studies at nearby sites [10,16]. Detailed information about the irrigation initiation and adjustment of irrigation duration is provided in Létourneau et al. [10]. For Global, Inter and Local, water was applied when the threshold was reached through a single irrigation event lasting between 40 and 60 min depending on crop development. For Pulse, water application was divided into 2 events separated by 2 to 3 h. Each event lasted for half of the irrigation time applied to the other treatments. For Global, all three replicates were simultaneously irrigated when the average SWP from all three replicates reached the threshold. This reproduced standard management of a large field based on average measurements from randomly positioned monitoring stations. For all other treatments, each replicate was individually monitored and irrigated to reproduce site-specific management. The experimental setup is summarized in Figure 2. Henceforth, Global, Inter and Local will be referred to as the management scale treatments.
The management scale treatments were established in this manner because when the experiment was set up, in early May 2012, a detailed analysis of the spatial variability of soil properties was not yet available. However, soil survey maps of the area and preliminary experiments on nearby sites had shown that a gradient in arable soil depth, soil texture, or in the proportion of rock fragments would most likely be in the north-south direction. Additionally, to maximize the cultivated surface within the field, to facilitate agricultural operations, and to protect against dominant winds, the beds were almost exclusively oriented in the north-south direction in the area surrounding the experimental site. Hence, the most practical way to reduce the irrigation management scale in this field was to reduce the length of the irrigation management zones in the north-south direction.
to the other treatments. For Global, all three replicates were simultaneously irrigated when the average SWP from all three replicates reached the threshold. This reproduced standard management of a large field based on average measurements from randomly positioned monitoring stations. For all other treatments, each replicate was individually monitored and irrigated to reproduce sitespecific management. The experimental setup is summarized in Figure 2. Henceforth, Global, Inter and Local will be referred to as the management scale treatments.  Two to three times a week, depending on weather conditions and crop development, the marketable yield was harvested from six full-length beds per plot. Weekly, two samples of 50-100 fruits were randomly collected from each plot. These fruits were counted and weighed to calculate the average fruit size. Within these samples, two fruits of equal size and color were selected for a firmness measurement with a penetrometer (FT-2, QA supplies, Norfolk, VA, USA) and sugar content measurement (brix index) with a refractometer (PAL-1, ATAGO CO, Tokyo, Japan). The volume of applied irrigation water was obtained from weekly surveys of flow totalizers installed on the supply pipe of every plot. Early in the season, three healthy plants of average size were selected and marked at two randomly placed sampling sites per plot. From early June to mid-July, these plants were used for a weekly measurement of leaf coverage (the longest leaf-covered distance in two perpendicular directions) and crown diameter. These measurements were halted when leaves from adjacent plants began to intertwine. The linear part of the growth curves was then used to calculate the leaf coverage and crown diameter growth rates by linear regression.
Soil samples were collected from the beds between 25 May and 1 June 2012. Georeferenced samples were collected from four main transects in the east-west direction ( Figure 1A). Each main transect consisted of six sampling locations approximately 12 m apart (9 beds). Six additional northwest-southeast or northeast-southwest diagonal sub-transects, each consisting of six sampling locations 2.5 m apart, were positioned to cover most of the field ( Figure 1A), for a total of 60 sampling locations.
The soil texture was measured at every other point along the main transects at a depth of 15 cm and at every third point at a depth of 30 cm following the standard hydrometer method [42]. Soil cores (54 mm high, 83 mm i.d.) were collected at all points (15 cm depth) and every third point (30 cm depth) from the main transects to measure the bulk density (ρ bulk ) and the saturated hydraulic conductivity (Ks-core) using the constant head method [43]. Thirteen of the cores collected at a 15 cm depth were randomly selected to obtain the water retention curves with pressure cells [44], and the unsaturated hydraulic conductivity with the multistep outflow method [45]. The easily available water (EAW) was calculated from the water retention curves as the difference between the water content at bed capacity (−2.5 kPa) and the water content at the irrigation initiation threshold (−18 kPa). The SWP at bed capacity is the average SWP 2 to 3 hours after an irrigation event. Details about this concept are given in Létourneau et al. [10]. The in situ hydraulic conductivity (Ks-in situ) was also measured at all sampling points with a constant head infiltrometer [46]. The proportion of rock fragments, i.e., particles with a diameter greater than 2 mm in the 0-30 cm horizon, was measured from samples collected at all points on the main transects using the procedure described by Grossman and Reinsch [47]. The SWP at which root water uptake should become limiting for root uptake was calculated from the 13 unsaturated hydraulic conductivity curves at a depth of 15 cm, following the approach presented in Rekika et al. [15]. As a first step, the parameters from the Gardner (1958) hydraulic model, where K (LT-1) is the unsaturated hydraulic conductivity, K s G (LT-1) is the fitted saturated hydraulic conductivity, h (L) is the SWP, and α* (L-1) is a fitted shape parameter, were obtained by inverse modelling of the experimental data with the minpack.lm package in R [48]. The optimized parameters were then used as inputs in the equations proposed by Rekika et al. [15], which were obtained by adapting an analytical solution of a steady-state form of Richards' equation considering evaporation, root water uptake, and SWP distribution with depth [33] to the determination of the critical SWP for root uptake (h c ), where h 1 is the asymptotic soil pressure head to maintain throughout all the domain (cm), L is the depth of the root zone (cm), z obs is the depth for which h c is calculated, S0 is the root water uptake rate (cm h −1 ) and q 0 is the steady-state soil surface evaporation (cm h −1 ). The model parameters were selected to represent conditions observed in the field. Hence, L was set to 30 cm, z obs was set to 15 cm, S 0 was set to 8.33 × 10 −3 cm h −1 (transpirative demand of 4 mm over a period of 12 hours with a crop coefficient of 0.75 and a rooting depth of 30 cm), and q 0 was set to 1.38 × 10 −4 cm h −1 (evaporation rate of 0.5 mm per day). The same approach was also used to calculate h c for theoretical soils from 12 major classes. Twelve theoretical unsaturated hydraulic conductivity curves ( Figure 3) were obtained by direct calculations with the class-averaged parameters of the van Genuchten (1980) hydraulic model, as predicted by the ROSETTA model [49]. Then, the Gardner model parameters were obtained by inverse modelling with the theoretical curves for h c calculations.
Agronomy 2019, 9, x FOR PEER REVIEW 6 of 18 where h1 is the asymptotic soil pressure head to maintain throughout all the domain (cm), L is the depth of the root zone (cm), zobs is the depth for which hc is calculated, S0 is the root water uptake rate (cm h −1 ) and q0 is the steady-state soil surface evaporation (cm h −1 ). The model parameters were selected to represent conditions observed in the field. Hence, L was set to 30 cm, zobs was set to 15 cm, S0 was set to 8.33 × 10 −3 cm h −1 (transpirative demand of 4 mm over a period of 12 hours with a crop coefficient of 0.75 and a rooting depth of 30 cm), and q0 was set to 1.38 × 10 −4 cm h −1 (evaporation rate of 0.5 mm per day). The same approach was also used to calculate hc for theoretical soils from 12 major classes. Twelve theoretical unsaturated hydraulic conductivity curves ( Figure 3) were obtained by direct calculations with the class-averaged parameters of the van Genuchten (1980) hydraulic model, as predicted by the ROSETTA model [49]. Then, the Gardner model parameters were obtained by inverse modelling with the theoretical curves for hc calculations. Statistical analysis of the marketable yield results was performed by using a linear mixed-effects repeated measures model with the nlme package in R [50]. For water use, fruit quality, and growth indicator measurements, a standard generalized linear model ANOVA and least significant Statistical analysis of the marketable yield results was performed by using a linear mixed-effects repeated measures model with the nlme package in R [50]. For water use, fruit quality, and growth indicator measurements, a standard generalized linear model ANOVA and least significant difference (LSD) multiple comparison tests were conducted using R [51]. For soil properties (when enough data points were available), after verification that anisotropy did not have any effect on the spatial structure of the variance [52], an omni-directional experimental variogram (γ(h)) was computed using where γ(h) is the semi-variance, x i is the sampling location, h (meters) is the lag distance between observation pairs, N is the number of observation pairs and Z is the value of the soil property. Standard variogram models were fitted to the experimental data with a restricted maximum likelihood approach using the geoR package in R [53], and the best fit model was chosen based on the Akaike information criterion (AIC). When the semi-variogram showed a significant structure, the spatial variability of the soil properties was interpolated and mapped by ordinary kriging.

Soil Properties
As shown in Table 1, the soil texture at the experimental site, particularly the sand content, showed significant variability at both measurement depths. A substantial (up to 35%) and variable (coefficient of variation = 85%) fraction of rock fragments 2 mm was observed in the 0-30 cm horizon, which is consistent with the description of the Orléans series. An average of 82.7% of the rock fragments were within the 2.0-9.5 mm range, with most of the remaining fraction smaller than 25.2 mm. As shown in Figure 4, there was an important gradient in the proportion of rock fragments in the north-south direction. It must be noted that the interpolation of the rock fragment proportion over the experimental area presented in Figure 4 was not based on a good fit of the variogram model to the experimental variogram. Hence, the interpolation standard error is important, and Figure 4 only provides a rough visual description of the variability of this property. Still, this is consistent with studies reporting that the variability of the rock fragment fraction is greater along than across the slope of a rock formation [54,55]. This also confirms that reducing the length of the experimental plots in the north-south direction was a proper way to evaluate the impacts of various management scales.
Such variability in the rock fragment proportion was expected to affect hydraulic conductivity [54,56,57]. As shown in Table 2, the saturated hydraulic conductivities, whether measured in situ or by laboratory measurements on soil cores at two depths, were very high and varied as much as 200 cm h −1 between sampling sites. At a 15 cm depth, the coefficient of variation for hydraulic conductivity was similar to that of the rock fragment proportion (Table 1), and it was much higher at 30 cm depth. Such variability in saturated hydraulic conductivity could considerably influence the leaching of water and fertilizers under the root zone following irrigation events.  Significant variability was also observed in soil water retention and unsaturated hydraulic conductivity ( Figure 5). The variability in the water retention curves mostly affected the saturated and residual water contents, which varied by as much as 20% and 10%, respectively. However, all the individual curves presented a similar shape ( Figure 5A). A very sharp decrease in water content was observed between soil saturation and a SWP around −10 kPa, at which point the curves almost become flat, and the water content slowly decreases with decreasing SWP. Despite the variability in water retention curves, EAW was similar between the sampling sites, with an average value of 0.05 cm 3 cm −3 (Table 2). In all cases, water availability could rapidly become a limiting factor at SWP lower than −10 to −15 kPa, as from this point great differences in SWP are required to induce small changes  Significant variability was also observed in soil water retention and unsaturated hydraulic conductivity ( Figure 5). The variability in the water retention curves mostly affected the saturated and residual water contents, which varied by as much as 20% and 10%, respectively. However, all the individual curves presented a similar shape ( Figure 5A). A very sharp decrease in water content was observed between soil saturation and a SWP around −10 kPa, at which point the curves almost become flat, and the water content slowly decreases with decreasing SWP. Despite the variability in water retention curves, EAW was similar between the sampling sites, with an average value of 0.05 cm 3 cm −3 (Table 2). In all cases, water availability could rapidly become a limiting factor at SWP lower than −10 to −15 kPa, as from this point great differences in SWP are required to induce small changes in water content. The unsaturated hydraulic conductivity curves ( Figure 5B) exhibited a similar bi-modal pattern. The shape of the curves is similar to what would be expected for a coarse sand, which suggests that soil texture alone is not a good indicator of the soil hydraulic behavior in the presence of rock fragments. The shape of the water content and hydraulic conductivity curves suggests that gravitational and preferential water flows have a major influence on soil hydrodynamics at this site. Gravitational flows generally take place in macro-pores, whose proportions are affected by the particle size distribution and orientation of soil and rock particles.
The soil bulk density (Table 2) was also variable but to a lesser extent than the other measured properties. This is most likely because ρbulk measurement is a simple procedure less prone to experimental error than water retention or hydraulic conductivity properties. However, in these specific conditions, ρbulk is most likely influenced by the proportion and orientation of schist fragments with a diameter >2 mm, which in turn could influence water availability and movement through the soil matrix [54,[58][59][60]. The unsaturated hydraulic conductivity curves ( Figure 5B) exhibited a similar bi-modal pattern. The shape of the curves is similar to what would be expected for a coarse sand, which suggests that soil texture alone is not a good indicator of the soil hydraulic behavior in the presence of rock fragments. The shape of the water content and hydraulic conductivity curves suggests that gravitational and preferential water flows have a major influence on soil hydrodynamics at this site. Gravitational flows generally take place in macro-pores, whose proportions are affected by the particle size distribution and orientation of soil and rock particles.
The soil bulk density (Table 2) was also variable but to a lesser extent than the other measured properties. This is most likely because ρ bulk measurement is a simple procedure less prone to experimental error than water retention or hydraulic conductivity properties. However, in these specific conditions, ρ bulk is most likely influenced by the proportion and orientation of schist fragments with a diameter >2 mm, which in turn could influence water availability and movement through the soil matrix [54,[58][59][60].

Spatial Variability of Soil Properties
For all soil properties measured on this site (ρ bulk , Ks-core, K s -in situ, EAW, particle size distribution, rock fragments proportion), the geostatistical analysis showed no spatial structure (not shown). The experimental semi-variograms mainly exhibited an important nugget effect and no range, meaning that the spatial variability of the soil properties was of a random nature at the measurement scale used in this experiment. This could be explained by the fact that tillage and cultural bed preparation practices have an important influence on the extent of the within-bed variability of these properties [54,59,61]. Because these operations are typically repeated every season, a random temporal variability component is added, making it unrealistic and impractical to establish irrigation management zones based on soil properties mapping for this site. This result differs from previous studies on the use of geostatistics for precision irrigation purposes [5,7,62].

Root Uptake Limiting SWP
The root uptake limiting SWP (h c ) calculated with the Rekika et al. [15] equations from the unsaturated hydraulic conductivity curves obtained from measurements at the experimental site and for the theoretical class averaged soils are presented in Table 3. The h c values from the unsaturated conductivities measured at the experimental site varied between −8 and −17.9 kPa with an average of −12.1 ± 2.4 kPa. This is consistent with critical SWP between −8 and −12 kPa previously reported for strawberries [10]. The variographic analysis for hc showed no consistent spatial structure and no spatial correlation between hc and the proportion of rock fragments or any other measured property could be found. Still the average h c values strongly contrasted from what would be expected for a theoretical silty clay loam (−23.2 kPa) and resembled more the h c value expected for a theoretical clay or a silty clay (Table 3). This suggests that, as expected, the rock fragment proportion strongly modifies the hydraulic behavior of the soil, but that its influence is randomly distributed throughout the experimental field, most likely due to the random effects of tillage and bed preparation. Table 3. Fitted Gardner (1958) hydraulic model parameters and calculated root uptake limiting SWP calculated from the observed properties and the theoretical soils (K sG and α* are the fitted Gardner (1958) model parameters for soil classes: S = sand, Si = silt, C = clay and L = loam).  Table 4 presents the proportion of time during the growing season spent above the average h c value. No difference between treatments were observed even if h c was on average exceeded for a considerable (24%) proportion of the season for all treatments. This suggests that an irrigation set point higher than −18 kPa would most likely have been more appropriate for strawberry production in this soil. This also suggests that even with the extent of the variation observed in soil properties at the experimental site, their effects on hydraulic behavior were random and most likely affected all treatments equally, regardless of the management scale. In this specific case, no benefits regarding the avoidance of condition limiting root water uptake could be obtained by applying a site-specific management. This also suggests that any location of the monitoring stations used for irrigation management would be as adequate as the other. As shown by the h c values obtained for the theoretical soils (Table 3), greater within-field differences in soil properties would probably be required to observe changes in the response to site-specific management. In this case, any effect of the management scale on crop production and water use are most likely attributed to practical effects of the management scale, for instance differences in the hardship to maintain the soil in the intended humidity range. Table 4. Number of irrigations, average achieved irrigation threshold (IT), average soil water potential (SWP) and water use resulting from the application of the experimental treatments. Values followed by a same letter do not statistically differ according to a least significant difference (LSD) multiple comparison test. The last row presents the p value from a standard ANOVA.

Irrigation Management and Water Use
From planting until early June, precipitation and the weekly application of fertilizers through the drip irrigation system provided sufficient water to maintain the soil water potential above the predefined irrigation threshold of −18 kPa. Consequently, the treatment application began in mid-June and was continued until early October. During that 104-day period, total rainfall was 324 mm (including 111 mm from two isolated rainstorms), and there were 80 days without rain. Figure 6 shows the SWP with respect to time for one replicate of each treatment in the month of September, which is the peak strawberry production period in the experimental area. It shows that all the treatments were maintained in the intended SWP range for most of the period. However, because of some practical delays to initiate irrigation (displacement time of the irrigator on the farm, valve adjustments, pump maintenance and availability), the achieved seasonal threshold (Table 4) was 2-4 kPa lower than the intended threshold for all scale management treatments. The intended threshold was exceeded more frequently in the Global treatment ( Figure 6A). It appeared that irrigation management based on the mean value of three SWP measurements caused the intended threshold to be frequently exceeded within the individual plots. Conversely, for the treatments with individually managed plots (Inter, Local, Pulse), the intended threshold was exceeded much less frequently. However, as discussed previously, the SWP was above h c for the same proportion of time for all treatments. Thus, these threshold exceedances should not have caused differences in plant root water uptake between the management scale treatments. However, a more precise evaluation of the effects of successions of short periods where the SWP is below h c could require more attention in further studies and would probably be easier to identify in controlled environment experiments. managed plots (Inter, Local, Pulse), the intended threshold was exceeded much less frequently. However, as discussed previously, the SWP was above hc for the same proportion of time for all treatments. Thus, these threshold exceedances should not have caused differences in plant root water uptake between the management scale treatments. However, a more precise evaluation of the effects of successions of short periods where the SWP is below hc could require more attention in further studies and would probably be easier to identify in controlled environment experiments. The Pulse treatment ( Figure 6D) allowed a better agreement between intended and achieved thresholds. Throughout the experiment, SWP changed more gradually between −10 and −18 kPa for the Pulse treatment than for the other treatments, thus allowing a longer period for intervention (displacement to the field, valve adjustments, pump startup) when SWP was nearing the intended threshold, which facilitated the initiation of irrigation at the intended threshold. This slower change suggests that the Pulse treatment maintained more available water within the root zone, most likely by reducing leaching or allowing greater horizontal distribution of water within the beds. This would be consistent with the simulation results obtained in other studies [30,36]. A more specific investigation would be required to confirm it.
The number of irrigation events and the total seasonal irrigation time (Table 3) did not differ significantly among treatments. An average of 0.73 irrigation events were initiated per day without rain (for the Pulse treatment, 1 irrigation initiation event implies 2 pulses). However, the total volume of applied irrigation water was significantly lower for the Global treatment compared with the Local treatment and represented water savings of 30%. The Inter and Pulse treatments required intermediate volumes of water, but were not statistically different from other treatments. The lower volume of water application over a similar irrigation time in the Global treatment can be partially explained by the longer delay between irrigation initiation and the achievement of the optimal emitter pressure, and hence the nominal flow rate, throughout the drip laterals. The average SWP for each treatment is also presented in Table 3. The SWP was significantly lower in the Global treatment as a result of the lower volume of water applied. The average SWP in the scale management treatments was within the optimal range for strawberry production reported in various studies [10,16,18,19], so there were practical benefits to reducing the management scale as it was easier to achieve the intended threshold, and significantly higher SWP was maintained at smaller management scales ( Table 3). The Pulse treatment had the added benefit of an increased intervention time to initiate irrigation when the SWP values approached the threshold. The Pulse treatment ( Figure 6D) allowed a better agreement between intended and achieved thresholds. Throughout the experiment, SWP changed more gradually between −10 and −18 kPa for the Pulse treatment than for the other treatments, thus allowing a longer period for intervention (displacement to the field, valve adjustments, pump startup) when SWP was nearing the intended threshold, which facilitated the initiation of irrigation at the intended threshold. This slower change suggests that the Pulse treatment maintained more available water within the root zone, most likely by reducing leaching or allowing greater horizontal distribution of water within the beds. This would be consistent with the simulation results obtained in other studies [30,36]. A more specific investigation would be required to confirm it.
The number of irrigation events and the total seasonal irrigation time (Table 3) did not differ significantly among treatments. An average of 0.73 irrigation events were initiated per day without rain (for the Pulse treatment, 1 irrigation initiation event implies 2 pulses). However, the total volume of applied irrigation water was significantly lower for the Global treatment compared with the Local treatment and represented water savings of 30%. The Inter and Pulse treatments required intermediate volumes of water, but were not statistically different from other treatments. The lower volume of water application over a similar irrigation time in the Global treatment can be partially explained by the longer delay between irrigation initiation and the achievement of the optimal emitter pressure, and hence the nominal flow rate, throughout the drip laterals. The average SWP for each treatment is also presented in Table 3. The SWP was significantly lower in the Global treatment as a result of the lower volume of water applied. The average SWP in the scale management treatments was within the optimal range for strawberry production reported in various studies [10,16,18,19], so there were practical benefits to reducing the management scale as it was easier to achieve the intended threshold, and significantly higher SWP was maintained at smaller management scales ( Table 3). The Pulse treatment had the added benefit of an increased intervention time to initiate irrigation when the SWP values approached the threshold.

Yield, WP and Growth Indicators
As indicated in Table 5, the irrigation management scale had no effect on marketable yield. Conversely, the Pulse water application method significantly improved yield by 22.4% compared to the other three treatments. This increase in yield is within the range observed by Wang et al. [63] for potatoes under high-frequency irrigation. WP showed a tendency to decrease as the management scale was reduced, which differs from the results of King et al. [20] and Périard et al. [14]. At the same management scale, the Pulse treatment led to a significant (36.8%) increase in WP relative to the Local treatment, which is consistent with the results of previous studies on potatoes and bell peppers [39,63]. In summary, the pulsed irrigation method was shown to increase yield by 22% and water productivity by 36% relative to conventional irrigation practices. More experiments are required to understand if this improvement is due to a lower average seasonal SWP, a better distribution of water within the growth beds, or a reduction in water leaching under the root zone. Table 5. Marketable yield, water use, fruit quality and growth indicators from the production experiment. Values followed by a same letter do not statistically differ according to a LSD multiple comparison test. The last row presents the p value from a standard ANOVA. Standard statistical annotation one asterisk means significant, two mean very significant. In contrast with studies reporting an influence of water application on fruit size [17,19], no difference was observed among any of the treatments in this experiment. As noted by Giné Bordonaba and Terry [64], this result can be attributed to varietal specificities. In their study, the fruit size of two out of five cultivars was not influenced by the irrigation management. The increase in yield associated with the Pulse treatment is thus attributable to a greater number of fruits. No significant differences in fruit quality (soluble sugar content or firmness) were observed among any of the treatments, but the fruits from the Local treatments showed a nearly significant tendency (p = 0.054) for a lower brix index. These results are consistent with observations by Giné Bordonaba and Terry (2010) [64], who found an increase in fruit sugar content in response to deficit irrigation, but they differ with the findings of Hoppula and Salo [18], who reported an increase in firmness and sugar content with a higher SWP for one out of three harvest periods.
The leaf coverage and crown diameter growth rates were similar among all treatments. This is consistent with other studies in which no effect of irrigation management on vegetative plant growth was observed [18,65]. The effect of irrigation on these parameters is also likely to be specific to the cultivar under study and is difficult to evaluate under field conditions because of the random effects of rain, wind, and damage from agricultural operations. In studies conducted under tunnels or greenhouses, leaf area and leaf architecture were reported to be affected by irrigation [66,67].

Conclusions
This study demonstrated that in a field where soil properties exhibited a great level of randomly distributed variability, intensive soil sampling and geostatistical analysis could not delineate irrigation management zones. High-definition sensing methods like electrical conductivity mapping with mobile instruments and water content mapping from spectral imagery might be more adequate for this purpose. Even if calculated critical irrigation thresholds did not present a structured spatial variability, they provided realistic estimates when based on the measured hydraulic properties and could be used as irrigation management guidelines. However, when calculated based on soil texture and pedotransfer functions, critical irrigation thresholds were unrealistic for strawberry production. More work is required to derive hydraulic properties from simple properties or from high-density remote sensing data.
Contrary to previous reports with other crops, the irrigation management scale had no influence on yield, and there was a trend for a decrease in water productivity with a reduction in management scale. However, pulsed irrigation water application based on SWP increased yield and WP considerably without affecting fruit quality. This suggests that further studies aiming to improve yield and WP in field strawberry production in similar conditions should focus on water application methods and not on the size of irrigation management units. More field-and/or simulation-based experiments designed to evaluate the impacts of different combinations of irrigation thresholds and water application methods are still needed to optimize irrigation practices in similar soil conditions. Author Contributions: Conceptualization, writing-review and editing, G.L. and J.C. methodology, validation, formal analysis, investigation, writing-original draft preparation, G.L. supervision, project administration, resources, funding acquisition, J.C. Acknowledgments: The authors would like to express their gratitude to Ferme Onésime Pouliot for providing the experimental site as well as for their technical and financial contributions to the research project. We also thank Jacques Gallichand and Claire Depardieu for their helpful comments.

Conflicts of Interest:
No conflict of interests to report.

Appendix A
with mobile instruments and water content mapping from spectral imagery might be more adequate for this purpose. Even if calculated critical irrigation thresholds did not present a structured spatial variability, they provided realistic estimates when based on the measured hydraulic properties and could be used as irrigation management guidelines. However, when calculated based on soil texture and pedotransfer functions, critical irrigation thresholds were unrealistic for strawberry production. More work is required to derive hydraulic properties from simple properties or from high-density remote sensing data.
Contrary to previous reports with other crops, the irrigation management scale had no influence on yield, and there was a trend for a decrease in water productivity with a reduction in management scale. However, pulsed irrigation water application based on SWP increased yield and WP considerably without affecting fruit quality. This suggests that further studies aiming to improve yield and WP in field strawberry production in similar conditions should focus on water application methods and not on the size of irrigation management units. More field-and/or simulation-based experiments designed to evaluate the impacts of different combinations of irrigation thresholds and water application methods are still needed to optimize irrigation practices in similar soil conditions. Author Contributions: The authors contributed to the present manuscript as follows: conceptualization, writing-review and editing, G.L. and J.C. methodology, validation, formal analysis, investigation, writingoriginal draft preparation, G.L. supervision, project administration, resources, funding acquisition, J.C. Acknowledgments: The authors would like to express their gratitude to Ferme Onésime Pouliot for providing the experimental site as well as for their technical and financial contributions to the research project. We also thank Jacques Gallichand and Claire Depardieu for their helpful comments.  Discharge rate of the emitters was calculated with:

Conflicts of Interest
where: • Q i = Discharge rate at the ith emitter (l h −1 ) • k = Unit conversion constant • x = Manufacturer supplied emitter exponent (0.5)