Effects of Irrigation Management on Chipping Potato ( Solanum tuberosum L.) Production in the Upper Midwest of the U.S.

: Irrigation is required for proﬁtable commercial potato ( Solanum tuberosum L.) production. Excessive or deﬁcit soil water availability during the growing season can have adverse effects on tuber yield, quality, and storability. A ﬁeld study was conducted during the 2018 and 2019 ﬁeld and storage seasons in Central Wisconsin, a region in the U.S. with a high volume of potato production, to evaluate the impacts of different irrigation rates on three chipping potato varieties, Hodag, Lamoka, and Snowden. The treatments were implemented during the late-tuber bulking and tuber maturation growth stages, and consisted of irrigation at 125%, 100%, 75%, and 50% of crop evapotranspiration (ET). Irrigation before the treatment period was at 100%ET for all plots. With the industry standard irrigation practice being at 100%ET, other treatments were designated as over-irrigation or deﬁcit irrigation. The impact of these watering rates on tuber yield and quality was evaluated at harvest, and tuber storage quality was assessed by measuring chip fry color and sugar concentrations at 0, 4, and 8 months of storage. It was found that compared to the standard practice, the over-irrigation treatment at 125%ET when tubers reached late bulking resulted in no signiﬁcant increase in total yield, marketable yield, tuber quality at harvest and during storage, as well as reduced irrigation efﬁciency (IE) and water-use efﬁciency (WUE). This treatment also increased nitrate leaching potential in both years. In comparison, deﬁcit irrigation at 75%ET or even 50%ET during the late season had no impact on tuber growth, could increase IE and WUE in one of the two years, and showed reduced drainage. In both years, irrigation rate had no signiﬁcant effects on hollow heart incidence, tuber speciﬁc gravity at harvest, and fry quality during the 8-month storage period. This study suggested that over-irrigation was not beneﬁcial for potato production in Central Wisconsin of the U.S., and deﬁcit irrigation during late tuber bulking and tuber maturation stages could potentially result in more sustainable water use while not penalizing tuber yield, quality and storability of chipping potatoes.


Introduction
Water is the most important compound in an active plant and constitutes more than 80% of the growing tissue [1]. Water for crop irrigation is becoming both scarce and expensive, making it essential for the growers to utilize precious water resource more sustainably [2]. Wisconsin ranks as No. 3 in the U.S. for its potato production, and the Central Sands region produces the majority of the state's potatoes [3,4]. The Central Sands is a sandy outwash plain featured by extremely coarse texture (90% sand) with a low soil organic matter content (<10 g kg −1 ) and a relatively shallow groundwater table (0.9-11 m) [5]. Nearly all of the irrigation water in the Central Sands is derived from groundwater. The impact of irrigation on depth of groundwater has become increasingly scrutinized in recent years [6]. Groundwater is a critical resource that must be preserved for long-term viability of the potato and vegetable industry and the local economy, and efficient water management will be crucial for preservation of the invaluable groundwater resource in this region [7]. field map of the 2018 study can be found in the Supplemental figure. A weather station (Campbell Scientific, Logan, UT, USA) was located within 500 m of the experimental plots to record precipitation, maximum and minimum soil temperature, and other weather conditions on a daily basis. Table 1. Irrigation treatments that were applied in this study.

Irrigation Rate Treatments
1-over irrigation Irrigation at 125% of adjusted evapotranspiration (ET) 2-standard irrigation Irrigation at 100%ET, standard treatment 3-deficit irrigation 1 Irrigation at 75%ET 4-deficit irrigation 2 Irrigation at 50%ET 5 *-deficit irrigation 3 Irrigation at 75%ET starting on 25 July, and then at 50% ET starting on 20 August * Rate 5 was only applied in 2019.
Seed potato pieces with an average weight of 71 g were planted on 1 and 2 May in 2018, and on 3 and 6 May in 2019. Emergence occurred on or around 20 May, and tuber initiation took place in the first week of June in each year. Irrigation management was performed by a linear overhead nozzle-controlled variable rate irrigation (VRI) system. Watering at 100% evapotranspiration (ET) is a standard practice at HARS as well as in the Wisconsin potato industry. In both years, all plots received irrigation at 100%ET before 25 July, which was the beginning of the irrigation treatment period. 25 July was chosen as it represented the average date across the three varieties to reach late tuber bulking. 20 August was chosen as the date to change application rates in treatment 5 (Table 1) because it was during the tuber maturation stage. All irrigation treatments ended at vine kill, which occurred on 10 September in 2018 and 9 September in 2019. The vine kill date was 14 days before harvest. Daily potential ET was estimated by the Wisconsin Irrigation Scheduling Program (WISP), where air temperature, relative humidity, wind speed, and cloud cover are used as inputs in the Priestley-Taylor equation. Potential ET values are then adjusted using the current crop canopy cover percentage. This ET estimation method was originally developed for potatoes, although it is not crop specific. WISP uses the latitude and longitude of the field location to pull relevant weather information from the University of Wisconsin-Extension's system of weather databases, which is a compilation of data from weather stations across the region.
For the standard irrigation treatment, the target soil moisture was 80% of field capacity, which is necessary to meet crop water needs when irrigating at 100%ET. Using WISP to make decisions, irrigation applications occurred when the farm manager determined that the soil water content would soon drop below the allowable depletion range. At that time, the crop would be irrigated to 80% field capacity.
Daily soil water balance was calculated to estimate percolation below the rooting zone using the approach outlined by Bohman et al. [21]. In this approach, daily drainage (D i ) was estimated to be the difference between input of daily precipitation (P i ) and daily irrigation (IRR i ), output of daily crop evapotranspiration (ET i ), and daily change in soil moisture storage (∆S i ) (Equation (1)). In this framework, the sum of ∆S i and the soil moisture storage on the previous day (S i−1 ) cannot exceed field capacity (FC). This implies that D i equals to the quantity of soil moisture in excess of FC, and leaching of excess soil moisture takes place within a 24-h period: All other cultural practices, including fertilization and pest management, were based on the University of Wisconsin-Madison Extension's recommendation [22]. Total nitrogen applied was 340 kg N ha −1 , and application was split between planting, hilling, tuber initiation, canopy closure, early tuber bulking, and mid tuber bulking.
To monitor soil moisture content during the growing season, one water content reflectometer sensor (WCR) (model CS616, Campbell Scientific, Logan, UT, USA) was calibrated to the soil of the field and was installed at each of the two soil depths (30 cm as topsoil and 60 cm as subsoil) in representative plots for each irrigation treatment, so there were 8 sensors in 2018 and 10 sensors in 2019. Each WCR sensors recorded soil volumetric water content at 15-min intervals to a datalogger. Due to high operation cost, there were not multiple replications of each set of WCR sensors. As a result, these data were used for general management decisions and not for statistical analysis.

Harvest Practices
All sub-plots were harvested on 24 September 2018 and 23 September 2019 using a two-row digger (Gallenberg Technologies, Antigo, WI, USA). The tubers were washed and graded for their individual weight using an AgRay Vision X-ray grading machine (AgRay Vision, Acampo, CA, USA). Tubers that were misshapen, rotten, or green due to sunlight exposure were culled before being graded. Cull weights were recorded in their own category and counted towards total yield from each plot.
Grading was categorized following the size profile based on tuber diameter: 0-4.76 cm, 4.76-8.89 cm, and greater than 8.89 cm. Marketable yield was then calculated as the yield of non-culled tubers with diameter between 4.76 and 8.89 cm.
Ten marketable tubers from each sub-plot were randomly chosen to determine specific gravity as well as internal defects. Specific gravity was measured using a Weltech PW-2050 Dry Matter Assessment System. After taking the weight of a 2000-3000 g sample of tubers in the air and in water, the system used the following formula to calculate the specific gravity for the sample (Murphy and Goven, 1959): weight in air (weight in air − weight in water) To assess internal defects, a separate sample of 10 tubers (4.76-8.89 cm) from each sub-plot were cut in half from stem end to bud end to evaluate incidence of hollow heart. Irrigation efficiency (IE) was defined as the ratio of the yield produced by the crop to the amount of water supplied to the crop through irrigation. IE was calculated as yield (in megagrams per hectare) divided by the total amount of irrigated water (with a unit of cm) applied over the course of the entire growing season. Water-use efficiency (WUE) was defined as the ratio of the yield produced by the crop to the amount of water (irrigation and precipitation) supplied to the crop. WUE was calculated as yield divided by the total amount of precipitation and irrigation in the entire growing season. Both total yield and marketable yield were used for these calculations.

Storage Management
After harvest, about 20 kg of healthy potatoes from each sub-plot were placed into the University of Wisconsin Hancock Potato Storage Research Facility, which was maintained at 12.8 • C for the purpose of wound healing for 30 days. Temperatures were then decreased at a rate of 0.3 • C per day to a final holding temperature of 8.9 • C. The sprout inhibitor, isopropyl N-(3-chlorophenyl)carbamate (CIPC, Decco Chemical, at 78.6% active ingredient), was applied to all tubers as a thermal aerosol at a rate of 22 mg kg −1 right before ramping temperature down to the holding temperature. Relative humidity was maintained at 95 ± 3% and CO 2 level was kept below 2000 mg m −3 throughout the entire duration of storage. Fry quality data of the three varieties were collected at harvest, and 16 and 32 weeks after harvest in both years.

Fry Color Analysis
At each sampling point, 10 tubers from each sub-plot were cut in half lengthwise. Three 1-mm thick slices through the stolon attachment point were removed from one half-tuber using a customized cutter and fried in vegetable oil at 185 • C for 2 min 10 s or until bubbling ceased. Fried chips were crushed for fry color analysis indicated by Hunter L-values using a D25LT colorimeter (HunterLab, Reston, VA, USA). Hunter Lvalues from 0 to 50 indicate dark chip fry color, and those from 51 to 100 indicate light chip fry color. Samples with Hunter L-values of 55 or higher are acceptable to the U.S. potato chip processing industry. Fry color at harvest in 2019 was not tested due to equipment breakdown.

Sugar Analysis
Sucrose and glucose concentrations were determined on the same 10-tuber sample used for the fry color determination by the method of Sowokinosetal [23] with modifications. Two hundred grams of tuber samples were cut from the center of the tubers. Samples were then juiced using an Acme Supreme Juicerator (Acme Equipment, Spring Hill, FL, USA) mixed with a 50 mM phosphate buffer. The samples were placed in a refrigerator for a minimum of 30 min, then a 1 mL sample from the top of the container was extracted. A YSI 2900 Select Biochemistry Analyzer (Sigma-Aldrich, St. Louis, MO, USA) was used to analyze sucrose and glucose concentration.

Data Analysis
Data analysis for the data collected at harvest was performed following a split-plot design. Fixed effects included irrigation rate, variety, and their interaction. Random effects included replicated blocks within the whole plot. Data from 2018 and 2019 seasons were analyzed separately due to significant interaction between year and the irrigation treatment. Analysis of variance (ANOVA) was conducted in R (version 3.6, R Foundation for Statistical Computing, Vienna, Austria) using the 'sp.plot' function from the 'agricolae' package [24,25]. The sp.plot function automatically specifies fixed and random effects in a split-plot design based on the order that the factors are input into the script. Means were separated using a protected Tukey's honestly significant difference (HSD) test at the α = 0.05 level. This test was performed with the 'HSD.test' function from the 'agricolae' package [25]. To assess the effect of storage length on storage quality of tubers harvested from each replication of irrigation rate by variety combination, data were analyzed using the SAS (version 3.8, SAS Institute Inc., Cary, NC, USA) PROC MIXED procedure as a split-split plot in time, with repetition*irrigation specifying the whole plot and repetition*irrigation*variety specifying the subplot as random and 'storage length' time repeated with an Autoregressive(1) covariance structure [26].

In-Season Weather, Irrigation, and Soil Moisture
Accumulated growing degree days (GDD) from 1 May to 1 July from 2011 to 2019 are shown in Figure 1. For potatoes grown in Wisconsin, the period between 1 May and 1 July is the most critical for tuber sprouting, tuber initiation and early tuber bulking, all of which are directly correlated with final tuber yield and quality. The year of 2018 had the most GDD between planting and early tuber bulking during this 9-year period, whereas the year of 2019 showed the least. The result suggests that the two seasons when the study was conducted covered a range of extreme growing conditions in central Wisconsin over the last 9 years.
Total precipitation between 1 May and 30 September was 49.80 cm in 2018 and 50.72 cm in 2019. Total precipitation during the treatment period (25 July to vine kill) was 27.61 cm in 2018 and 9.96 cm in 2019. In 2018, there were three periods of heavy precipitation (>2.5 cm per rain event) that largely disrupted the irrigation (Figure 2a). These periods took place between 14 June and 22 June, between 19 July and 23 July, as well as between 21 August and 9 September. The intense rainfall between 21 August and 9 September caused the second irrigation treatment of treatment 5, which was expected to start on 20 August at the rate of 50%ET, to be completely disrupted and not to be imposed. In 2019, there were erratic rainfall events (<2.5 cm) over the course of the treatment period (Figure 2b).   Volumetric water content (VWC) at two soil depths are shown in Figures 3 and 4 for the 2018 and 2019 growing seasons, separately. In both years, the soil water content of the topsoil (0-30 cm) and subsoil (30-60 cm) had similar values among treatments before 25 July, which was the start of the irrigation treatment period. Thereafter treatment 1 (over-irrigation) had higher VWC values than rate 2 (standard irrigation), particularly at the subsoil depth. Rate 3 and rate 4 (the two deficit irrigation levels) had lower VWC values at the top soil depth compared to rate 2 (standard irrigation), and the difference was more pronounced in 2018. In 2019, rate 5 had similar soil VWC changes as rate 4 although the topsoil values were slightly higher under rate 5 during most of the treatment period ( Figure 4).
Throughout the growing season, the topsoil had more changes in VWC values within each day, which was indicated by the daily spike approximately at midnight, than the subsoil. The different irrigation treatment effects were more closely reflected by the VWC data in 2018 than in 2019, likely due to precipitation in 2019 throughout the treatment period that can disrupt the irrigation treatments (Figure 2b).
Weekly water balance values under each irrigation rate of this experiment are displayed in Figures 5 and 6. Compared to the weekly ET values of May and June in 2018, the values in 2019 were considerably lower. This is likely related to the lower growing degree days in the early portion of the 2019 growing season (Figure 1). During the treatment period, it is clear that compared to over-irrigation, deficit irrigation at 75%ET or 50%ET when tubers reach late bulking can substantially reduce drainage, which is directly related to nitrate leaching potential to the groundwater. The high drainage values in late August and early September of 2018 should be caused by the series of large rainfall events during that period (Figure 2a). Overall, total drainage under 75%ET and 50%ET was 74.9 mm and 112.4 mm less than that under over-irrigation in 2018. Total drainage under 75%ET, 50%ET, and 75-50%ET was 69.2 mm, 103.8 mm, 81.6 mm less than the value under over-irrigation in 2019.

Yield and Tuber Size
ANOVA results for traits measured at harvest were separately shown for 2018 and 2019 in Table 2. Different irrigation rates did not result in any significant difference of total yield or marketable yield in either year. However, there was significant varietal difference of total yield and marketable yield in both years (p < 0.001). In 2018, Snowden showed the highest total yield, and Hodag and Snowden showed similar marketable yield (Table 3). In 2019, Hodag had the highest numbers for both yield categories (Table 3). This indicated the premium yield potential of the new Hodag variety. In comparison, Lamoka performed the most poorly for both total and marketable yield in both years (Table 3). There were not any significant two-way interactions between irrigation treatment and variety for any of the measurements done at harvest in either year (Table 2).  All three varieties had lower total and marketable yield in 2019 than in 2018, possibly due to varying growing conditions (Figure 1). For Hodag, total yield was 18% lower in 2019, with marketable yield 50% lower. For Snowden, total yield was 25% lower in 2019, with marketable yield 77% lower. For Lamoka, average total yield was similar between the two years, but the average marketable yield was 74% lower in 2019 than in 2018 (Table 3).

Tuber Quality at Harvest
In both years, irrigation rate did not have any significant impact on incidence of hollow heart (Table 2). Overall, hollow heart incidence in all varieties was lower than 20%. In addition, no significant difference of specific gravity was observed under different irrigation rates (Table 2).
Comparatively, varietal difference led to significant difference in hollow heart incidence (Table 2), with Hodag being higher than Snowden and Lamoka in both years. In 2019, Hodag was the only variety that had any occurrence of hollow heart, due to its larger tuber size (Table 3). Varietal difference also had a significant effect on specific gravity ( Table 2). In both years, Hodag and Snowden had similar specific gravity values, while the values of Lamoka were significantly lower (Table 3).

Water Use
For both total yield and marketable yield across the varieties, irrigation efficiency (IE) was significantly lower under the over-irrigation rate than under all other irrigation rates in 2018 (Tables 4 and 5). In 2019, IE for total yield under over-irrigation and standard irrigation was lower than those under other treatments (Table 4), but irrigation rate had no effect on IE for marketable yield in this year ( Table 2). Compared to the other treatments, over-irrigation consistently showed the lowest IE in both years (Tables 4 and 5). Water-use efficiency (WUE), which considers use of both precipitation water and irrigation water, was also significantly affected by irrigation rate. WUE for total and marketable yield in 2018 was significantly lower under the over-irrigation rate than all the other treatments (Tables 4 and 5). In 2019, WUE for total yield was lower in the overirrigation treatment than under the rate 4 at 50%ET or rate 5 that had two tiers of irrigation reduction during the treatment period (Tables 4 and 5).
It is observed that in 2019, there was significant varietal difference of IE and WUE for producing total yield (Table 2), with Hodag being similar to Snowden, and Lamoka being the lowest (Table 6).  Table 7 shows the ANOVA results for the storage quality in 2018 and 2019, respectively. Irrigation rate did not influence sugar concentrations or chip fry color in either year. In 2018 the interaction between storage length and variety had a significant effect on sucrose and glucose ( Table 7). The sugar concentrations, particularly glucose, significantly increased in Snowden at the eight-month sampling point, whereas the numbers stayed low over the course of the storage season in Hodag and Lamoka (Tables 8 and 9). In 2019, the interaction between storage length and variety significantly impacted sucrose, glucose, and fry color (Table 7). Snowden again showed significant spikes of sucrose and glucose and a symptom of darkening fry color, indicated by lower Hunter-L values, at the end of the storage season, but Hodag and Lamoka maintained low sugar levels and light fry color (Table 10). Table 7. ANOVA results of sucrose, glucose and chip fry color for 2018 and 2019 separately. Values that were not significant at α = 0.05 are indicated as '-'.

Source
Sucrose Concentration (mg g −1 FW) It is interesting to note that in both years of the study, Lamoka showed the lowest levels of sucrose and glucose and lightest fry color after eight months of storage compared to other varieties at other sampling points (Tables 8-10), suggesting its premium longterm storability. Alternatively, Snowden showed a substantial increase of sucrose and glucose and worsening of fry color after eight months of storage (Tables 8-10), due to its susceptibility to senescent sweetening.

Deficit Irrigation Applied during Late Season Provides Potential for Water-Use Sustainability
A 2013 study [27] conducted in Wisconsin on 7 round white chipping varieties showed that the impact of reducing irrigation by 25% during crop maturation was less than expected as yield difference of individual varieties between non-reduced irrigation and reduced irrigation was not significant. Bussan [27] suggested that reducing irrigation while not penalizing yield would possibly provide an economic benefit to growers because of reduced pumping costs and retained productivity of the crops. Hassan et al. [28] examined the effects of deficit irrigation at different potato growth stages in Bangladesh and found that in both years of the study, deficit irrigation at the tuber maturation stage resulted in higher yield compared with the treatment with no deficit irrigation. Additionally, they found that treatments that applied deficit irrigation during both the late tuber bulking stage and the tuber maturation stage did not lead to a significant yield loss when compared to the non-deficit treatments.
In our study, deficit irrigation at 75%ET or 50%ET when tubers reached the late-bulking stage did not affect yield, at-harvest tuber quality, and tuber storage quality (Tables 4 and 9) when compared to the standard rate at 100%ET or over-irrigation at 125%ET. The increased irrigation efficiency and water-use efficiency at 50%ET during the late growing season showed that irrigating below the standard level at 100%ET when tubers mature results in more efficient use of water (Tables 4 and 5). In addition, reducing irrigation rate can reduce drainage, which is directly related to nitrate leaching to the groundwater (Figures 5 and 6). Comparatively, irrigating at 125%ET could lead to significantly lower irrigation efficiency and water-use efficiency (Tables 4 and 5), and higher nitrate leaching potential (Figures 5  and 6) compared to the standard practice. In both years, the highest irrigation efficiency and water-use efficiency was found at the lowest deficit irrigation level at 50%ET. This was true for water use for both total yield and marketable yield. These results are in line with findings from Badr et al. [29].
Aside from overall water savings, irrigation management can interact with nitrogen management and nitrogen losses to groundwater, because the combination of intense irrigation and nitrogen requirements, together with the shallow rooting depth of potato, can result in high nitrate leaching potential [30]. Waddell et al. [31] showed that allowing the soil to become drier between irrigation events was beneficial in reducing N leaching, especially when rains occurred right after irrigation. Sexton et al. [32] showed that if applying irrigation to bring the soils back to field capacity, significant nitrate leaching occurred when thunderstorms followed soon after irrigation or fertilization. Saffigna et al. [33] allowed the soil to become drier between irrigation events, thus providing additional soil-water storage to retain unforeseen rainfall. The smaller spikes of volumetric water content after large rainfall event (>2.5 cm) in both topsoil and subsoil at the 50%ET rate compared to other irrigation treatments (Figures 3 and 4) also indicated the lower nitrate leaching potential under the lowest deficit irrigation level in our study.
Overall our VWC data generally showed the trends that were expected under different irrigation rates, there were times when the trends were not reflected in the data. In particular, the standard irrigation rate at 100%ET in 2018 showed higher topsoil VWC than expected. This finding illustrated the spatial variability in soil hydrologic properties that exist throughout a field, especially in the topsoil. This might suggest the weakness of using soil sensors for closely monitoring water status throughout a field, and conveys the benefits of utilizing remote sensing that covers spatial variability.

Over-Irrigation Is Not Beneficial for Potato Yield, Quality, and Storage Quality
Our study showed that changes in late season irrigation amount did not impact total yield or marketable yield for all three chipping varieties in both years. Over-irrigation at 125%ET did not result in any improvement in total yield, marketable yield, internal quality, and specific gravity in either year of the study. Therefore, it is suggested that over-irrigation did not provide benefits to potato production in the context of this research. Over-irrigation has been considered by some farmers as an "insurance practice" to ensure crop yield gain. However, when the soil is maintained as excessively wet, root and tuber needs for oxygen are not met [13]. Tubers that develop in frequently saturated soil show swollen lenticels that provide openings for soft rot bacteria to enter [34], and are subject to secondary growth [35]. Tuber yields and marketable grade can decrease with over-irrigation [13,36]. Additionally, late-season excessive soil moisture could lead to tuber storage quality issues such as pink rot and Pythium leak 1997).
Further research is needed in the Central Sands region of Wisconsin to investigate potato performance under over-irrigation in dry years when precipitation does not substantially disrupt irrigation treatments. One strategy is to cover field plots with movable hoop houses during precipitation events, however, significant amount of labor and precise weather forecasting will be needed in order to make this strategy successful.

Extreme Weather Conditions Increase Year-to-Year Potato Crop Production Variation
When a precipitation event is larger than 1.3 cm in the central sands region of Wisconsin, irrigation is not needed [37], preventing the irrigation treatments of this study from being effectively applied. In 2018, the long period of precipitation from 21 August to 9 September substantially disrupted rate 5, the split deficit irrigation rate that was supposed to start with 50%ET on 20 August. A significant amount of commercial potato fields in that year were also affected by this long stretch of heavy rainfall that caused enlarged tuber lenticels associated with pathogen infection and storage quality issues, or crops not to be harvested due to saturated field conditions [38].
In 2019, the abnormally low amount of growing degree days from 1 May to 1 July contributed to lower yield and smaller tuber sizes compared to the previous year. It is important to note that hollow heart was almost non-existent in the 2019 tubers, most likely due to the reduced growing degree days and the consequently poor yield in this year.
A recent study has found that weather condition extremes, such as temperatures that are too hot and too cold, drought, and heavy precipitation, accounted for 18-43% of interannual variations in crop yield [39]. Extreme weather events are expected to affect the volatility of crop yields and are seen as the principle immediate threat to global crop production in the near future [40,41]. Among the various adaptation measures, those related to new crop variety development and modern irrigation technologies are considered more critical [42]. In this study, Hodag and Snowden showed higher IE and WUE than Lamoka (Table 8). This provides insights for future breeding efforts of developing new outstanding potato varieties that use water more efficiently.
Modern irrigation technologies, such as precision irrigation, are being widely promoted to resolve challenges regarding improving crop productivity under current variable weather patterns. Variable rate irrigation (VRI) technology used in this study is an innovative system that applies varying rates of irrigation amount along the length of the center pivots to address specific soil, crop, and/or other conditions [43]. VRI changes the spatial distribution of water applied to a given field without affecting the amount and timing of application. VRI is a potentially useful tool for enhancing farm profitability, irrigation water productivity (yield produced per unit of irrigation water applied), and water quality [44,45]. O'Shaughnessy et al. [46] had demonstrated the utility of VRI in deficit irrigation management to improve the efficiency of irrigation (e.g., reducing the fraction of irrigation that runs off, drains below the root zone or is otherwise not available for extraction by plant roots), because VRI may allow for a more spatially uniform managed level of water stress across a given field, by accounting for soil and other variability.

Increased Understanding of the Difference of Agronomic Performance between Chipping Varieties
Although irrigation treatment did not have a significant impact on traits measured at harvest and during storage, substantial varietal difference was observed in both total and marketable yield in both years (Table 2), which provided information for future agronomic research for new potato variety development. For example, Hodag showed higher total and marketable yield in both years, whereas Lamoka was the poorest yielder (Table 3). On the other hand, Hodag tubers had a higher incidence of hollow heart than Snowden and Lamoka tubers, supporting its features of greater yield potential and growing larger tubers.
Hodag and Lamoka are found to be more premium for retaining long-term storage quality, because Snowden showed senescent sweetening with increased sugar concentrations and worsened fry color during the late storage season (Tables 8-10). Senescent sweetening is a process that potato tubers in storage undergo a developmentally controlled, age-dependent shift in metabolism that results in an accumulation of sucrose, glucose, and fructose. Senescent sweetening negatively affects potato chip color and is manifested as a blush of dark color that begins in the center of the chip and progresses outward as sweeten-ing becomes more severe [47]. Increasing the wound-healing temperature at the beginning of the storage season can mitigate the senescent sweetening issue on Snowden [38].

Conclusions
Sustainable irrigation is required for profitable commercial potato production. Our study suggested that over-irrigation at 125%ET when tubers reached late bulking was not beneficial for potato production in Central Wisconsin of the U.S., and deficit irrigation at 75%ET or even 50%ET during the late season could potentially result in more efficient and sustainable water use while not penalizing tuber yield, quality, and storability of chipping potatoes. While reductions in irrigation did not significantly affect production, choice in variety had major impacts, with Hodag showing the greatest potential for yield and tuber size, and Snowden showing the greatest risk of late season senescent sweetening.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.