Vol.:(0123456789) Potato Research (2024) 67:1603–1624 https://doi.org/10.1007/s11540-024-09695-3 1 3 Effects of Abiotic Stress Associated with Climate Change on Potato Yield and Tuber Quality Under a Multi‑environment Trial in New Zealand Allan Banasihan Siano1,2  · Nick Roskruge1 · Huub Kerckhoffs1 · Svetla Sofkova‑Bobcheva1 Received: 1 August 2023 / Accepted: 7 January 2024 / Published online: 6 March 2024 © The Author(s) 2024 Abstract In the 2018/19 growing season, a multi-environment trial in Opiki, Hastings, and Ohakune located in three different regions of the North Island of New Zealand was conducted to evaluate responses of selected potato cultivars to abiotic stress associ- ated with climate change. Heat and drought stresses were evident with supra-optimal temperatures (> 25 °C) in Opiki and Hastings and sub-optimal rainfall (< 500 mm) in Opiki, which influenced the different morpho-physiological characteristics of the potato crop, ultimately affecting yield and tuber quality. These abiotic stresses also increased the incidence of malformation, growth cracks, and second growth in tubers reducing the total and marketable tuber yields by 43% and 45%, respectively. In addition, the genotype × environment analysis showed that Ohakune had the most favourable environmental conditions for potato production since all cultivars in this site had superior marketable tuber yields. ‘Taurus’ was the most stable and adapta- ble cultivar across trial sites (wide adaptation), whilst ‘Hermes’ and ‘Snowden’ were more adapted under Opiki and Hastings conditions (specific adaptation), respec- tively. As established in this study, heat and drought stresses have significant effects on the morpho-physiology, yield, and tuber quality of commercial potato cultivars in New Zealand. Keywords Climate adaptation · Drought stress · Genotype · Heat stress · Marketable tuber yield · Physiological disorders * Allan Banasihan Siano allanbsiano@gmail.com 1 School of Agriculture and Environment, Massey University, Palmerston North, New Zealand 2 Philippine Council for Agriculture, Aquatic, and Natural Resources Research and Development, Los Baños, Laguna, Philippines http://orcid.org/0009-0002-4235-9516 http://crossmark.crossref.org/dialog/?doi=10.1007/s11540-024-09695-3&domain=pdf 1604 Potato Research (2024) 67:1603–1624 1 3 Introduction Potato (Solanum tuberosum L.) is an economically important vegetable crop grown in over 150 countries. It ranks fifth after sugarcane, maize, wheat, and rice as main crops in terms of production volume (FAO 2022). In New Zealand, potato is the top vegetable crop based on production area (10,901 ha), production volume (554,052.53 MT), and export value ($110 million) (FAO 2021). However, the New Zealand potato industry is faced with many challenges such as stringent government policies, high production costs due to the economic impacts of the global pandemic, and extreme weather events (Potatoes New Zealand 2022). Amongst these chal- lenges, the extreme weather events, e.g., droughts, heatwaves, and flooding, often associated with climate change are gradually impacting potato production in New Zealand. In general, potato production under climate change is expected to decline by 18–32% (without adaptations) and by 9–18% (with adaptations) (Hijmans 2003). The influence of environment, in the genotype × environment interaction, plays a significant role in the yield variation of potato crops. Agro-meteorological factors significantly contribute to the variability of potato crop yields under experimental (20–60%) (Kalbarczyk 2004) and field (40–80%) conditions (Kalbarczyk 2004; Trawczyński 2009). In New Zealand, the western and southern regions are anticipated to receive more rain and experience warmer temperatures, whilst eastern and northern regions will experience double to triple the frequency of drought episodes by 2040 (New Zea- land Agricultural Greenhouse Gas Research Centre 2012). In a 60-year (1958–2018) climate data analysis in selected potato-growing areas in New Zealand (Siano et.al. 2018), atmospheric temperature and rainfall showed an increasing and decreasing trend, respectively, which are consistent with reports from the Ministry of Environ- ment (2016) and the IPCC (2014). Before releasing a cultivar, genotypes of high potential are typically evaluated at different locations and over several years to determine their degree of adaptation (Acquaah 2012). However, in most cases, a thorough assessment of adaptation and yield stability characteristics of these cultivars is lacking before their recommenda- tion (FAO 2002). The multi-environment trial (MET) is a practice used to evalu- ate crops over multiple locations to identify superior cultivars (Acquaah 2012). The main effects of the environment (E), genotype (G), and associated G × E interaction are determined using appropriate statistical systems. Through MET, yield can be estimated, and stable cultivars in different environments can be identified (Crossa 1990). MET is often used as a tool for selecting superior genotypes for planting in new sites. In potatoes, MET has been applied in testing and screening for improved tuber yield and disease resistance (Paget et al. 2014). This study aimed to evaluate the performance of selected commercial potato cultivars, through MET, as affected by abiotic stresses in different regions in the North Island of New Zealand. 1605 1 3 Potato Research (2024) 67:1603–1624 Materials and Methods Site Description and Local Climate The multi-environment trial (MET) was established in the summer of 2018/19 in three commercial potato farms located in three regions (different ecozones) of the North Island of New Zealand i.e., Ohakune (39.500 ° S, 175.459 ° E, 593 m above sea level (masl), central North Island), Opiki (− 40.465 ° S, 175.481 ° E; 3 masl; southwest North Island), and Hastings (− 39.624 ° S, 176.738 °E; 17 masl; eastern North Island). The Ohakune site is a conventional-irrigated farm, ex-pasture, orthic allophanic-volcanic ash soil (Landcare Research 2020), with an average temperature of 10.2 °C and average annual rainfall of 1500 mm. The Opiki site is a conventional- rainfed farm, ex-pasture, humic organic soils (Landcare Research 2020), with an average temperature of 13.2 °C, and average annual rainfall of 990 mm. The Hast- ings site is an organic-irrigated farm, ex-cover crop, mottled fluvial recent soil-allu- vial sand silt or gravel soil (Landcare Research 2020), with an average temperature of 14.0 °C and average annual rainfall of 770 mm. Crop Management and Trial Layout Five offshore-bred potato cultivars—‘Agria’ (reference cultivar), ‘Hermes’, ‘Fianna’, ‘Taurus’, and ‘Snowden’—and two New Zealand-bred potato cultivars—‘Rua’ and ‘Ilam Hardy’—from certified seed suppliers were used in the trial. The experimental setup was a Randomized Complete Block Design with four blocks and cultivars were fully randomized in each block. Pre-sprouted seed potatoes 40–60 mm in diameter were manually planted in ridged rows 12 m long with between and within row spac- ing of 100 and 30 cm, respectively. The grower’s cultural management practices in each site, including land preparation, nutrient and water management, and pest and disease management, were adopted to simulate local growing practices. In the Hast- ings trial site, being an organic farm, a nylon mesh cover (0.6 mm pore size) was used as a physical barrier to exclude the tomato–potato psyllid (TPP) (Bactericera cockerelli), which is the insect vector of Candidatus Liberibacter solanacearum that causes zebra chip in potatoes. Weather Data The atmospheric temperature was recorded using a temperature logger DS1921G- F5# Thermochron (www. maxim integ rated. com) in Opiki and Harvest Air Tempera- ture Sensor (Harvest Electronics NZ Ltd.) in Hastings. The atmospheric temperature in Ohakune and the soil temperature and daily rainfall (within a 15-km radius of the trial sites) in all sites were obtained from CliFlo—New Zealand’s web-based http://www.maximintegrated.com 1606 Potato Research (2024) 67:1603–1624 1 3 national climate database. The accumulated growing degree days (GDD) at 4.4 °C base temperature was calculated using the CIP (2013) formula. Agronomic and Morpho‑physiological Data At 65  days after planting (DAP), the plant height (mm), stem diameter (mm), and total plant leaf area (cm2) were recorded. For the total plant leaf area, 10 fully developed leaf sub-samples per plant were randomly selected and weighed, and the leaf area was measured using LI-3100C Leaf Area Meter (LI-COR Bio- sciences, USA, www. licor. com). The total plant leaf weight was divided by the weight of the sub-samples and multiplied by the measured leaf area to get the total plant leaf area. The net photosynthesis (µmol CO2 m–2  s–1), transpiration rate (mmol H2O m–2 s–1), and stomatal conductance (mol H2O m–2 s–1) were measured using LI- COR portable photosynthesis system IRGA LI-6400 model (www. licor. com). At the same time, the chlorophyll content (SPAD value) was recorded using a port- able chlorophyll metre CCM-200plus (www. optis ci. com). All measurements were done on the 3rd or 4th fully expanded leaf from the shoot tip between 1000 and 1400 h, with eight samples measured per cultivar. Yield and Tuber Quality Data At 130 DAP, manual harvesting was done at the different sites. The total yield (t ha−1), marketable tuber yield (t ha−1), number of tubers per plant, number of market- able tubers per plant, dry matter content (DMC) (%), harvest index or HI (%), spe- cific gravity (SG), and percentage of tuber physiological disorders were recorded. Tubers were classified as marketable if within the marketable size standards and with minimal or no defects (physiological disorders) or unmarketable if undersized or had severe to very severe defects. Percent tuber physiological disorders (e.g., mal- formation, second growth, growth cracks, enlarged lenticels, feathering, and vascu- lar browning) were determined by dividing the number of tubers with physiological disorders over the total number of tubers per plant and multiplying by 100. Ten randomly selected marketable tubers from each cultivar were used to measure the DMC and SG. The DMC was determined by oven-drying the chopped potato tubers at 70  °C for five days. The fresh and dry weights were measured, and the DMC was calculated by dividing the dry weight by the fresh weight and multiplying by 100. The HI was estimated using the ratio of tuber dry weight: total plant dry weight and multiplying by 100 (MacKerron and Heilbronn 1985). The SG was determined by first weighing the tubers using a top-loading balance (weight in air) and then weighing them whilst suspended in water (weight in water) (Haase 2003). The weight in air was then divided by the difference between weight in air and weight in water. http://www.licor.com http://www.licor.com http://www.optisci.com 1607 1 3 Potato Research (2024) 67:1603–1624 Statistical Analysis All the statistical analyses were carried out in R statistical package version 3.4.0 (R Core Team 2017). The significant difference was analysed using Analysis of Variance (ANOVA), and mean values were compared amongst cultivars by using Tukey’s HSD test if the overall test was significant. Genotype × Environment Analysis The additive main effect and multiplicative interaction (AMMI) method was used, where an integrated ANOVA and principal component analysis (PCA) were applied to analyse the MET (Crossa et  al. 1991; Zobel et  al. 1988). ANOVA was used to determine the main effects of genotypes and environments and the PCA for the residual multiplicative interaction amongst genotypes and environments. After detecting a significant G × E interaction, the data were graphically analysed by biplot technique to simultaneously classify genotypes and environments (Hongyu et  al. 2014; Kempton 1984; Zobel et al. 1988). Through this method, cultivars with high productivity and wide adaptability (wide adaptation), as well as cultivars with spe- cific adaptability (specific adaptation), were identified. Results Weather Data In general, Ohakune had the most favourable weather condition amongst the sites, whilst Opiki and Hastings had more extreme weather conditions during the 2018/19 potato growing season (Table 1). The maximum average atmospheric temperature (AT) in Opiki (25.1  °C) was much higher than the historical average (20.4  °C). Similarly, the number of days with supra-optimal AT i.e., > 25 °C and > 30 °C was highest in Opiki at 103 and 24 days, respectively. Throughout the season from plant- ing (October) to harvesting (February), the average daily maximum AT increased by about 6 °C relative to the historical average. The use of mesh covering in Hast- ings created a microclimate resulting in about a 2 °C increase in AT. Correspond- ingly, the number of days with > 25  °C increased by 80% (from 49 to 61  days) and > 30 °C increased by 38% (from 8 to 21 days). In contrast, Ohakune had only 20 days of temperatures > 25 °C and one day of temperature > 30 °C. The increase in AT in Opiki and Hastings resulted in a faster accumulation of growing degree days or GDD (Fig. 1). In the same way, the seasonal average soil temperature (ST) was much higher in Opiki (19.5 °C) and Hastings (19.0 °C) than in Ohakune (15.7 °C). On the other hand, total rainfall was lower by 40% (from 642.1 to 388.0  mm) in Ohakune, by 10% (from 465.8 to 423.0 mm) in Opiki, and by 11% (from 455.4 to 406.0 mm) in Hastings. 1608 Potato Research (2024) 67:1603–1624 1 3 Ta bl e 1 W ea th er d at a in th e O ha ku ne , O pi ki , a nd H as tin gs tr ia l s ite s d ur in g th e 20 18 /1 9 po ta to g ro w in g se as on Tr ia l s ite s/ de ta ils M on th ly A tm . T em p. (° C ) A ve . A tm . Te m p. (° C ) N o. o f d ay s w ith su pr a- op tim al A tm . T em p So il Te m p. (° C ) To ta l r ai nf al l ( m m ) O ct N ov D ec Ja n Fe b M ar Se as on al (2 01 8/ 19 ) H ist or ic al A ve . (1 98 1– 20 10 ) O ha ku ne 38 8. 0 64 2. 1   S ea so n a v 9. 8 11 .3 14 .6 17 .2 15 .6 15 .1 13 .9 15 .7   H ist av 9. 7 11 .6 14 .0 15 .4 15 .4 13 .5 13 .3   M ax av 15 .1 16 .9 20 .4 23 .4 22 .4 21 .8 20 .0 19 .8   H ist av 14 .2 16 .7 19 .0 21 .2 21 .2 18 .9 18 .5   M in av 4. 5 5. 8 8. 8 10 .9 8. 8 8. 3 7. 9 12 .6   H ist av 5. 3 6. 4 8. 9 9. 7 9. 7 8. 1 8. 0   N o. o f d ay s > 25  °C 0 0 2 10 7 1 20   N o. o f d ay s > 30  °C 0 0 0 1 0 0 1 O pi ki 42 3. 0 46 5. 8   S ea so n a v 13 .8 16 .0 18 .3 19 .9 18 .4 18 .3 17 .5 19 .5   H ist av 12 .5 14 .2 16 .4 17 .8 18 .2 16 .6 16 .0   M ax av 18 .1 24 .5 26 .4 27 .3 25 .8 24 .7 25 .1 23 .2   H ist av 16 .6 18 .3 20 .6 22 .5 22 .9 21 .2 20 .4   M in av 8. 1 9. 4 12 .3 14 .2 11 .0 13 .4 11 .3 16 .5   H ist av 8. 4 10 .0 12 .1 13 .1 13 .5 11 .9 11 .5   N o. o f d ay s > 25  °C 1 18 24 25 19 16 10 3   N o. o f d ay s > 30  °C 0 3 7 9 4 1 24 H as tin gs 40 6. 0 45 5. 4   S ea so n a v 12 .1 15 .0 17 .5 20 .1 17 .9 17 .1 16 .6 19 .0   H ist av 13 .6 15 .3 17 .7 18 .8 18 .7 17 .1 17 .5   M ax av 18 .7 20 .1 22 .7 25 .9 24 .3 23 .2 22 .5 23 .1 1609 1 3 Potato Research (2024) 67:1603–1624 Ta bl e 1 (c on tin ue d) Tr ia l s ite s/ de ta ils M on th ly A tm . T em p. (° C ) A ve . A tm . Te m p. (° C ) N o. o f d ay s w ith su pr a- op tim al A tm . T em p So il Te m p. (° C ) To ta l r ai nf al l ( m m ) O ct N ov D ec Ja n Fe b M ar Se as on al (2 01 8/ 19 ) H ist or ic al A ve . (1 98 1– 20 10 )   H ist av 18 .6 20 .4 22 .7 23 .9 23 .5 22 .1 22 .5   M in av 5. 2 9. 4 12 .6 13 .8 11 .2 11 .4 10 .6 17 .0   H ist av 8. 5 10 .2 12 .8 13 .7 13 .9 12 .2 12 .6   N o. o f d ay s > 25  °C 0 1 8 19 12 9 49   N o. o f d ay s > 30  °C 0 0 0 4 4 0 8 H as tin gs (u nd er m es h co ve r)   S ea so n a v 13 .6 16 .2 19 .0 21 .2 19 .3 18 .6 18 .9 19 .0   H ist av 13 .6 15 .3 17 .7 18 .8 18 .7 17 .1 17 .5   M ax av 19 .7 20 .9 22 .4 28 .2 26 .6 22 .6 24 .1 22 .8   H ist av 18 .6 20 .4 22 .7 23 .9 23 .5 22 .1 22 .5   M in av 7. 5 11 .4 14 .2 15 .8 13 .6 13 .4 13 .7 15 .5   H ist av 8. 5 10 .2 12 .8 13 .7 13 .9 12 .2 12 .6   N o. o f d ay s > 25  °C 1 2 9 27 19 3 61   N o. o f d ay s > 30  °C 0 0 1 10 10 0 21 A ll hi sto ric al av er ag es w er e ba se d on a 3 0- ye ar p er io d (1 98 1– 20 10 ) Se as on av se as on al av er ag e, M ax av m ax im um a ve ra ge , M in av m in im um av er ag e, H is t av h ist or ic al av er ag e 1610 Potato Research (2024) 67:1603–1624 1 3 Abiotic Stress Effects on Morpho‑physiological Characteristics of Potato Crops Weather extremes in the 2018/19 season accompanied by drought and heat stress affected the morpho-physiology of the potato crops. The influence of the envi- ronment (trial sites) was significant on plant height (p < 0.001), stem diameter (p < 0.001), total plant leaf area (p < 0.01), net photosynthesis (PN) (p < 0.05), tran- spiration rate (E) (p < 0.001), stomatal conductance (gs) (p < 0.001), and chlorophyll content (ChlC) (p < 0.01) (Fig. 2, Supplementary Table 1). The influence of culti- vars was also significant for plant height (p < 0.001), stem diameter (p < 0.001), total plant leaf area (p < 0.01), and ChlC (p < 0.01) but not significant for PN, E, and gs (Fig. 3, Supplementary Table 1). Relative to Ohakune, plant height and total plant leaf area were higher in Hast- ings by 43% (112.1 cm) and by 59% (6892 cm2), respectively, but the stem diam- eter was lower by 13% (13.9 mm). In contrast, plant height and total plant leaf area were lower in Opiki by 4% (75.2  cm) and by 66% (1489 cm2), respectively, but higher stem diameter by 6% (16.9 mm) relative to Ohakune. Amongst the cultivars, ‘Fianna’ had the most vigour in terms of plant height (100.5 mm) and stem diameter a a a b c b c b c ab a a bb aa a c bb c Fig. 2 Morpho-physiological parameters of commercial potato cultivars as affected by environments in the Ohakune, Opiki, and Hastings trial sites during the 2018/19 potato growing season at 65 DAP (days after planting). Means ± standard error is presented and letters indicate significant differences calculated on Tukey’s HSD test (p ≤ 0.05) y = 11.243x - 80.193 R² = 0.9962 y = 14.464x - 140.04 R² = 0.9948 y = 14.806x - 87.029 R² = 0.9978 1 10 20 30 40 50 60 70 80 90 100 110 120 130 Fig. 1 Accumulated growing degree days (GDD) in the Ohakune, Opiki, and Hastings trial sites during the 2018/19 potato growing season at 4.4 °C base temperature 1611 1 3 Potato Research (2024) 67:1603–1624 (19.1 mm). On the other hand, ‘Hermes’ had the highest total plant leaf area (5605 cm2) but was not significantly different from ‘Ilam Hardy’ (5506 cm2). On physiological parameters, PN was lower by 16.5% (19.2 µmol CO2m–2 s–1) in Hastings but was not significantly different from Opiki (21.0 µmol CO2 m–2   s–1), and Opiki compared with Ohakune (23.0 µmol CO2 m–2 s–1). The E was higher by 61.8% (5.5 mmol H2O m–2 s–1) in Hastings and by 55.9% (5.3 mmol H2O m–2 s–1) in Opiki, compared to Ohakune. Correspondingly, gs was higher by 129.5% (0.342 mol H2O m–2 s–1) in Hastings and by 52.3% (0.227 mol H2O m–2 s–1) in Opiki relative to that in Ohakune. ChlC was reduced by 22.4% (25.8 SPAD) in Hastings, whilst increased by 12.2% (36.7 SPAD) in Opiki. On cultivar effects, ‘Agria’ had the high- est ChlC (37.1), followed by ‘Ilam Hardy’ (33.8) and ‘Taurus’ (33.8). Abiotic Stress Effects on Final Yield and Tuber Quality of Potato Crops The influence of the environment (trial sites) was very highly significant (p < 0.001) for total yield (TY), marketable tuber yield (MY), dry matter content (DMC), har- vest index (HI), and specific gravity (SG); highly significant for percent yield loss (YL); and significant (p < 0.05) for the number of tubers plant−1 (NTP) and number of marketable tubers plant−1 (NMTP) (Fig. 4, Supplementary Table 2). In the same way, the influence of cultivar was very highly significant (p < 0.001) for all param- eters, except for the YL and NTP, which is highly significant (p < 0.01) (Fig. 5, Sup- plementary Table 2). Ohakune had the highest TY (72.4 t ha−1), MY (57.8 t ha−1), NTP (24.7), NMTP (10.7), DMC (19.8%), HI (70.9%), and SG (1.078). On the other hand, Opiki had lower TY (41.0 t ha−1) and MY (34.8 t ha−1) but was not significantly different from Hast- ings’ TY (41.9 t ha−1) and MY (31.7 t ha−1). Hastings had the most YL (24%), followed by Ohakune (20%) and Opiki (15%). Opiki had 35.6% less NTP (15.9) and 14.9% NMTP (9.3), whilst Hastings had 18.6% less NTP (20.1) and 6.5% NMTP (10.0). In the same way, Opiki and Hastings had significantly lower DMC (16.7%, 18.5%), HI (66.5%, 62.6%), and SG (1.069, 1.071), respectively, compared to Ohakune. ababcbc a c d e a aa bc bcc c c b b ab b b aaa a a a a aa a aa ab a a a a cd c bc a a cd da a a a Fig. 3 Morpho-physiological parameters of commercial potato cultivars as affected by cultivars in the Ohakune, Opiki, and Hastings trial sites during the 2018/19 potato growing season at 65 DAP (days after planting). Means ± standard error is presented and letters indicate significant differences calculated on Tukey’s HSD test (p ≤ 0.05) 1612 Potato Research (2024) 67:1603–1624 1 3 Amongst the cultivars ‘Rua’ had the least TY (37.9 t ha−1), whilst ‘Hermes’ (56.2 t ha−1) had the highest TY, which was not significantly different from ‘Taurus’ (56.0 t ha−1), ‘Ilam Hardy’ (53.6 t ha−1), ‘Agria’ (53.3 t ha−1), ‘Fianna’ (53.1 t ha−1), and ‘Snowden’ (51.9 t ha−1). ‘Taurus’ had the highest MY (49.1 t ha−1) but was not significantly different from ‘Snowden’ (47.5 t ha−1), which also had the least YL (8.5%). The reference cultivar, ‘Agria’, had relatively low MY (37.9 t ha−1) and was amongst the cultivars with high YL (28.9%). ‘Hermes’ had the greatest NTP (27.3) but was not significantly different from ‘Ilam Hardy’ (24.0), and ‘Ilam Hardy’ from ‘Taurus’ (23.4). ‘Hermes’ also had the greatest NMTP (14.4), which was not signifi- cantly different from ‘Taurus’ (13.2) and ‘Snowden’ (12.3), whilst ‘Agria’ had the least NMTP (5.9). On the other hand, ‘Snowden’ had the highest DMC (20.4%) but was not sig- nificantly different from ‘Hermes’ (19.9%), whilst ‘Ilam Hardy’ had the lowest (16.5%). ‘Snowden’ also had the highest SG (76.3%) but was not significantly differ- ent from ‘Taurus’ (75.4%), whilst ‘Agria’ had the lowest (61.9%). On the other hand, ‘Hermes’ and ‘Snowden’ had the highest SG (1.080), but it was not significantly dif- ferent from ‘Taurus’ (1.077). a a a b b b a b cc b b a a b ca b ab c ba b b Fig. 4 Yield and tuber quality parameters of commercial potato cultivars as affected by environments in the Ohakune, Opiki, and Hastings trial sites during the 2018/19 potato growing season at 130 DAP (days after planting). Means ± standard error is presented and letters indicate significant differences calculated on Tukey’s HSD test (p ≤ 0.05) a a aa a a ab a ab a b bcbc c a a ab b c a b d c cd d f e f ab de de a ab abc bc cd dee f cd cac bc ab c a e a bc a a bc bc Fig. 5 Yield and tuber quality parameters of commercial potato cultivars as affected by cultivars in the Ohakune, Opiki, and Hastings trial sites during the 2018/19 potato growing season at 130 DAP (days after planting). Means ± standard error is presented and letters indicate significant differences calculated on Tukey’s HSD test (p ≤ 0.05) 1613 1 3 Potato Research (2024) 67:1603–1624 Abiotic Stress Effects on Tuber Physiological Disorders of Potato Crops The influence of the environment (trial sites) was very highly significant (p < 0.001) on the formation of tuber physiological disorders such as malformation, growth cracks, enlarged lenticels, feathering, and vascular browning, and highly significant (p < 0.01) for second growth (Fig. 6, Supplementary Table 3). In the same way, the influence of cultivars was also very highly significant (p < 0.001) for all parameters but not significant for vascular browning (Fig. 7, Supplementary Table 3). Feathering (28.1%), malformation (19.1%), and growth cracks (19.1%) in Ohakune; growth cracks (51.6%) and second growth (21.1%) in Opiki; and enlarged lenticels (33.1%) and second-growth (28.1%) in Hastings were the most common tuber physiological disorders observed (Fig. 8). Vascular browning was observed in Opiki and Hastings, whereas hollow heart though in a small number of tubers only was recorded in Hastings. a a a ba b b b c a a b a c b b b c Fig. 6 Tuber physiological disorders of commercial potato cultivars as affected by environments in the Ohakune, Opiki, and Hastings trial sites during the 2018/19 potato growing season at 130 DAP (days after planting). Means ± standard error is presented and letters indicate significant differences calculated on Tukey’s HSD test (p ≤ 0.05) a a b a b b ab bc c bc cd bcd cd c c ce ded e a a b c c d b c c cd e e d a a aaaaa Fig. 7 Tuber physiological disorders of commercial potato cultivars as affected by cultivars in the Ohakune, Opiki, and Hastings trial sites during the 2018/19 potato growing season at 130 DAP (days after planting). Means ± standard error is presented and letters indicate significant differences calculated on Tukey’s HSD test (p ≤ 0.05) 1614 Potato Research (2024) 67:1603–1624 1 3 ‘Agria’, ‘Rua’, ‘Fianna’, and ‘Ilam Hardy’ were more susceptible to tuber mal- formations. On the other hand, ‘Hermes’ and ‘Snowden’ were more susceptible to growth cracks, whilst ‘Taurus’ was less susceptible to growth cracks. ‘Ilam Hardy’ appeared to be susceptible to second growth i.e., heat sprouts and secondary tubers. Additionally, the following is the order of tolerance of the different cultivars to a cb d e f g h i j k l m n o Fig. 8 Tuber physiological disorders at final harvest (130  days after planting): in Ohakune, a feather- ing in ‘Agria’, b tuber cracking (superficial cracks) in ‘Hermes’, c second-growth (secondary tubers), and d tuber malformation in ‘Ilam Hardy’; in Opiki, e secondary tubers and f vascular browning in ‘Ilam Hardy’, g second-growth (heat sprouts) in ‘Agria’, h superficial cracks in ‘Hermes’, i second- growth (chain tubers) in ‘Ilam Hardy’; in Hastings j heat sprouts in ‘Ilam Hardy’, k enlarged lenticels in ‘Hermes’, l chain tubers in ‘Fianna’, m secondary tubers and n hollow heart in ‘Agria’, and o heat sprouts in ‘Hermes’ 1615 1 3 Potato Research (2024) 67:1603–1624 physiological disorders: ‘Taurus’ > ‘Snowden’ > ‘Hermes’ > ‘Fianna’ > Rua’ > ‘Ilam Hardy’ > ‘Agria’, during the 2018/19 season. Genotype × Environment Interaction Under the Multi‑environment Trial of Potato Crops The individual effects of genotype (G) and environment (E), and genotype × envi- ronment (G × E) interaction were very highly significant (p < 0.001) for marketable tuber yield based on ANOVA (Table 2). Further, the ANOVA of average marketable tuber yields across the seven cultivars and the three environments showed that 77.2% of the total SS (sum of squares) was attributed to environmental effects, 13.8% to genotypic effects, and 9.0% to G × E interaction. Moreover, the following is the order of the attained average marketable tuber yield of potato cultivars: ‘Taurus’ > ‘Snowden’ > ‘Hermes’ > ‘Fianna’ > ‘Agria’ > ‘Ilam Hardy’ > ‘Rua’, based on the generated biplot (Fig. 9). The average marketable tuber yields of ‘Taurus’ (49.1 t ha−1) and ‘Snowden’ (47.5 t ha−1) were well-above the mean (41.4 t ha−1) (intersect in the biplot); ‘Hermes’ (43.4 t ha−1) and ‘Fianna’ (42.9 t ha−1) were close to the mean, whilst ‘Agria’ (37.9 t ha−1), ‘Ilam Hardy’ (35.0 t ha−1), and ‘Rua’ (34.1 t ha−1) were below the mean. Discussion Weather Data The 2018/19 season had distinct variations and extreme weather patterns with Opiki and Hastings more negatively affected than Ohakune. It was the third hottest sum- mer on record, after 2017/18 and 1934/35 (NIWA 2020). The national average tem- perature reached 17.9  °C, which is 1.2  °C above the 1981–2010 average (NIWA 2020). The abnormally high temperature was due to marine heatwaves around New Zealand coastlines, especially in the Hawke’s Bay region (Salinger et al. 2020). The supra-optimal average atmospheric temperature (AT) (> 25 °C) in Opiki and Hast- ings for extended periods confirmed the presence of heat stress amongst the potato crops. Moreover, mesh coverings in Hastings aggravated this condition by creating Table 2 Analysis of variance for the genotype × environment interaction for average marketable tuber yield of seven commercial potato cultivars in the Ohakune, Opiki, and Hastings trial sites during the 2018/19 potato growing season Significance codes: ***, P < 0.001; ** P < 0.01 df Sum Sq Mean Sq F value Pr(> F) Environment (Env) 2 91,166 45,583 59.0470 6.692e-06 *** Rep Env 9 6948 772 2.9399 0.0021451 ** Genotype (Gen) 6 16,242 2707 10.3090 1.165e-10 *** Env × Gen 12 10,608 884 3.3666 0.0001043 *** Residuals 419 110,022 263 1616 Potato Research (2024) 67:1603–1624 1 3 a microclimate raising the temperature by 2 °C. In an earlier study, a 1 °C increase was observed under mesh-covered potato crops (Merfield et al. 2015). However, this was conducted on the South Island, where the AT is much lower than on the North Island. Soil temperature (ST) had the same trend as AT in all sites. Though, irriga- tion in Ohakune and Hastings could have had a cooling effect on ST, whilst the sub- optimal rainfall and non-irrigated Opiki could have aggravated the increase in ST. Total rainfall in all sites was below the optimum range (500–700 mm) for potato production. Rainfall was recorded below normal in the northern and western parts of Fig. 9 Additive Main Effect and Multiplicative Interaction (AMMI1) biplot showing the a Interaction Principal Component Analysis (IPCA1) (Factor1) vs. estimated average marketable tuber yield and b AMMI2 biplot showing the first two principal axes of interaction (IPCA2 vs. IPCA1) of seven commercial potato cultivars evaluated in three environments in the North Island of New Zealand during the 2018/19 growing season Ohakune Opiki Hastings Opiki Ohakune Hastings a b 1617 1 3 Potato Research (2024) 67:1603–1624 the North Island, including the Manawatu–Wanganui region (NIWA 2020). Across sites, particularly low rainfall was experienced in mid-spring (October) which was 40–80% lower than the historical average. This coincided with crop emergence and stolon initiation in Opiki. By early summer (December), heavy downpours occurred in Opiki and Hastings, caused by a storm and a wet December, respectively (NIWA 2020). However, this was followed by a generally drier mid-summer (January) that coincided with tuber bulking, up to early autumn (March). This resulted in a rapid decrease in soil moisture levels and progressed to severely dry soils in the northern and western regions of the North Island (NIWA 2020). Abiotic Stress Effects on Morpho‑physiological Characteristics of Potato Crops Weather extremes accompanied by heat and drought stress in the 2018/19 season have affected the morpho-physiology of the potato crops with Opiki and Hastings more negatively impacted than Ohakune. This was aggravated by mesh covering in Hastings that increased the plant height (but reduced stem diameter) and leaf area. This can be attributed to the relatively rapid vegetative growth and the increase in the number of smaller but more leaves under high temperatures (Gawronska et al. 1992; Hastilestari et al. 2018; Lafta and Lorenzen 1995; Rykaczewska 2015; Tang et al. 2018). In contrast, the drier conditions in Opiki reduced the plant height (but increased stem diameter), and leaf area. Drought stress can reduce the maximum leaf area index of potatoes by 29% in the field and by 63% in the greenhouse (Lahlou et al. 2003). Either or both the lowered leaf water potential that acts as a hydraulic signal and the hormonal signals (abscisic acid) produced in the root tips of potato crops may reduce leaf area expansion under drought conditions (Shahnazari et al. 2007), whereas the increase in stem diameter is a typical drought stress response where the accumulation of assimilates and conserving water into the stem are enhanced (Ohashi et al. 2008). On the other hand, the reduced PN in Hastings could be caused by heat stress that impaired the crop’s photosynthetic activity and increased the E and gs (Hammes and De Jager 1990; Hastilestari et al. 2018; Reynolds et al. 1990; Wolf et al. 1990). It was suggested that for every 5 °C increase above the optimum (20.0 °C) in leaf temperature of potatoes, the photosynthetic performance is reduced by 20% (Bur- ton 1981). The increased gs at high temperatures may improve photosynthesis either through increasing CO2 transport to the leaves or through the cooling effect on the leaves due to higher E (Wolf et al. 1990). However, if temperatures surpass a certain threshold, it may impair photosynthesis and reduce water use efficiency (Guoju et al. 2013). On the other hand, photosynthesis is hampered by the closure of stomata to slow down the transpiration rate under drought stress, which prevents the entry of CO2 into the leaves (Li et al. 2017) which may reduce growth rate and final yield. Meanwhile, the levels of photosynthetic pigments including chlorophylls a and b and carotenoids were reduced by up to 20% under heat stress (Hancock et al. 2014). In contrast, the increase in the ChlC under drought stress in Opiki could be due to the loss in leaf turgor pressure or a reduction in leaf growth (Ramírez et al. 2014). 1618 Potato Research (2024) 67:1603–1624 1 3 Abiotic Stress Effects on Final Yield and Tuber Quality of Potato Crops Ohakune had the highest TY, MY, NTP, HI, and SG, having temperatures close to optimum levels for potato growth and development, whilst the reduced rainfall was compensated by irrigation. In contrast, Hastings and Opiki had lower MYand tuber quality due to the presence of heat and drought stress. It is striking that the TY and MY in both sites were way below the national average of 51 t ha−1 (Potatoes New Zealand 2022). Heat stress alone impairs water use efficiency (Guoju et al. 2013), photosynthesis, and photosynthate production and enhances dry matter partitioning to the haulm, thereby promoting vegetative growth over tuber development (Gaw- ronska et al. 1992; Ghosh et al. 2000; Hancock et al. 2014; Hastilestari et al. 2018). As a consequence, the tuber fresh weight, dry weight, dry matter, and HI decreased (Ghosh et al. 2000; Hancock et al. 2014; Hastilestari et al. 2018). Heat stress at the early stage of potato crop growth can completely inhibit tuber initiation (Ghosh et al. 2000; Levy et al. 1991) and reduce the dry matter partitioning to the tubers, whilst at the latter stage, reducing SG (Ghosh et al. 2000), an important quality parameter for processing-type potatoes. On the other hand, drought stress at the early stage could reduce plant height (Shock et  al. 1992), delay canopy growth (Aliche et  al. 2018; Jefferies 1993; Martin et al. 1992), and decrease leaf area index (Lahlou et al. 2003); arrest stolon differentiation and tuber formation (Aliche et al. 2020); and at the latter stage, hasten crop maturity (Robins and Domingo 1956) and haulm senescence (van Loon 1981). The high MY observed for ‘Hermes’, ‘Taurus’, and ‘Snowden’ can be attributed to their vigorous plant growth e.g., greater leaf area that could enhance sunlight interception, higher DMC, NTP, and NMTP and lower non-marketable tubers com- pared to most tested cultivars. ‘Rua’ had the lowest TY and MY that agree with the recorded yield of the cultivar (McLeod 1973; Anderson et al. 1997). ‘Ilam Hardy’, on the other hand, had a relatively high TY but also had the highest YL. A common observation amongst varieties with higher MY, such as ‘Hermes’ and ‘Taurus’, is their inherent characteristics of producing a high NTP. HI and SG were reduced in Hastings and Opiki due to the elevated temperatures that resulted in a shift in biomass distribution away from the tubers and toward the aboveground parts (Hancock et al. 2014). High SG seemed to be a common char- acteristic of cultivars with relatively smaller tubers such as ‘Hermes’ (Amjad et al. 2016) and ‘Snowden’ (Baritelle and Hyde 2003; Canadian Food Inspection Agency 2020; Potato Association of America 2020). In contrast, SG seems to be lower amongst cultivars with larger tubers, such as ‘Agria’, ‘Fianna’, and ‘Ilam Hardy’. Abiotic Stress Effects on Tuber Physiological Disorders of Potato Crops The literature strongly suggests that many tuber physiological disorders are influenced by unfavourable environmental factors, including extreme atmospheric and soil tem- peratures, rainfall, and humidity (Hiller et  al. 1985; Lugt et  al. 1964; Rykaczewska 2015, 2017). In this study, the development of a specific tuber physiological disorder 1619 1 3 Potato Research (2024) 67:1603–1624 depended on the environment and the specific or combination of abiotic stress present during the growing season. In the cooler environment of Ohakune, feathering, a com- mon characteristic of immature tubers having thin skin that can easily be torn at har- vest and during handling (Hillerm and Thornton 2008) was the most common physio- logical disorder. Feathering could be due to the slower GDD accumulation throughout the growing season in Ohakune, as compared with Opiki and Hastings. In Opiki, the arid spring followed by a very wet early summer may have encour- aged growth crack formation. Whereas the moisture stress combined with supra- optimal AT and ST may have caused the premature breaking of tuber dormancy, which encouraged second-growth formation, i.e., heat sprouts, chain tubers, and sec- ondary tuber formation. A past study suggested that this could be an indirect effect of incomplete canopy coverage due to drought stress and high temperature (Shock et al. 1992). The common practice of ground storing in New Zealand could expose the ground stored-tubers to high ST which may enhance physiological disorders. The very highly significant influence of cultivars on the incidence of physi- ological disorders is indicative of strong differential susceptibility or tolerance to abiotic stress. Moreover, cultivars with naturally large and elongated tubers such as ‘Agria’, ‘Rua’, ‘Fianna’, and ‘Ilam Hardy’ were more susceptible to malforma- tion. In contrast, cultivars that have medium-sized circular tubers such as ‘Hermes’ and ‘Snowden’ were less susceptible to malformation but more susceptible to growth cracks. The rough and slightly netted skin characteristics of ‘Hermes’ and ‘Snowden’, appeared to be enhanced by abiotic stress, which developed into growth cracks. ‘Taurus’, on the other hand, with relatively smooth skin was less susceptible to growth crack formation. The reference cultivar ‘Agria’ was also susceptible to growth cracks and second growth. This study validates the preliminary findings in the 2017/18 season (Siano et al. 2018) where heat and moisture stresses were determined as likely causes of an array of tuber physiological disorders, which resulted in a reduction in marketable tuber yield. Genotype × Environment Interaction Under Multi‑environment Trial of Potato Crops The environments (trial sites) were diverse causing most of the variation in the average marketable tuber yield (MY). This can be indicated by the environment’s large sum of squares. Similarly, the very significant effect of the G × E interac- tion also showed that the genotypes (cultivars) had differential performance in the environments. According to Silveira et al. 2013, stable genotypes have lower influence on the G × E interaction and can be indicated by a lower IPCA1 value. In this MET, the cultivar ‘Taurus’ could be highlighted and identified as an ideal genotype because it had the highest average MY and stability (wide adaptation) compared with the reference cultivar ‘Agria’ and the other tested cultivars. In contrast, ‘Ilam Hardy’ was identified as an undesirable genotype owing to its low productivity and low sta- bility (IPCA1 values farthest from zero) (Gauch and Zobel 1988; Kempton 1984). 1620 Potato Research (2024) 67:1603–1624 1 3 On the other hand, ‘Agria’ showed intermediate stability only and an MY below the overall average. ‘Ilam Hardy’, ‘Snowden’, ‘Hermes’, and ‘Fianna’ were the most unstable cultivars. The unstable cultivars meant that they perform well under spe- cific environmental conditions only (specific adaptation). Genotypes and environ- ments located close to each other in the biplot have positive relations; thus, these enable the creation of a group of genotypes within an agronomic zone. For example, the cultivar ‘Hermes’ performed better in Opiki, whilst ‘Snowden’ in Hastings and ‘Fianna’ in Ohakune. All the environments had some cultivars with MY above the overall average; but only Ohakune had all cultivars above the overall average, indicating that this was the most favourable environment to obtain high yield. The possible reason for the high marketable tuber yields is the weaker presence of abiotic stress, the application of supplemental irrigation, and relatively higher natural fertility and favourable soil properties. Conclusions The multi-environment trial (MET) established that in the 2018/19 season, abiotic stresses associated with climate change i.e., heat and drought stress, were present and affected the morpho-physiology, yield, and tuber quality of commercial potato cultivars. Moreover, the cultivars exhibited differential susceptibility to abiotic stresses. Heat stress increases the plant height (but reduces stem diameter), leaf area, E, and gS, whilst drought stress suppresses plant height (but increases stem diame- ter), leaf area, and gS leading to a decrease in PN. Heat stress alone or in combination with drought stress reduces the TY, MY, NTP, DMC, and HI. The findings confirmed that abiotic stress increases the incidence of tuber physiological disorders, resulting in reduced MY. The occurrence of a specific physiological disorder depends on the environment and the type and timing of abiotic stress present during the growing season. Cultivars that have naturally large and elongated tubers are more susceptible to tuber malformation. In con- trast, cultivars that have medium-sized circular tubers are less susceptible to tuber malformation but more susceptible to growth cracks. The G × E analysis showed location-specific adaptations (agronomic zoning) amongst the tested commercial potato cultivars. The cultivar ‘Hermes’ is toler- ant to a combination of heat and drought stress, whilst ‘Snowden’ is tolerant to heat stress. ‘Taurus’, on the other hand, is the most stable and adaptable cultivar across the environments. Cultural practices that are unique in New Zealand potato production systems, such as the use of mesh covering as a physical barrier to protect crops from the tomato-potato psyllid, aggravate the increase in summer temperatures. Also, the practice of ground storing could expose the ground-stored potato tubers to abi- otic stresses (e.g., heat and drought stress), which dramatically affect the mar- ketable yield and tuber quality. Moreover, rainfed farming could be less feasible in the coming years with the changing rainfall patterns and frequent incidences of drought. Besides the application of appropriate cultural management (e.g., 1621 1 3 Potato Research (2024) 67:1603–1624 irrigation at critical stages, planting, and harvesting timing), identifying the right commercial potato cultivars with tolerance to abiotic stress (e.g., heat and drought stress) that would perform well either in a wide array of environments (wide adaptations) or specific locations (specific adaptations) should be conducted to ensure productivity and increase the resilience of the New Zealand potato indus- try to climate change. Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s11540- 024- 09695-3. Acknowledgements The authors would like to thank grower partners Mike Moleta of Freshpik Farms, Scott Lawson of True Earth Farm, and Keith Watson and Bryan Hart of Wilcox for providing farm inputs and land for the trials. Alex McDonald Ltd. and Eurogrow Ltd. are also acknowledged for providing the seed potatoes. The authors thank Professor Isaac Adeneyika for the statistical analysis. Author Contribution The research was part of A. Siano’s PhD dissertation. AS conceived and designed the experiments. N. Roskruge, H. Kerckhoffs, and S. Sofkova-Bobcheva, as panel of supervisors provided important inputs for the conceptualization and implementation of the different research activities. All of the authors in this paper read and approved the manuscript. Funding Open Access funding enabled and organized by CAUL and its Member Institutions. This work was carried out with research funding from the School of Agriculture and Environment, Massey Univer- sity, New Zealand, and Ph.D. scholarship funding from the New Zealand Aid Scholarship. Data Availability Data of agronomic, morpho-physiological, yield, and tuber quality of potato are included in this manuscript. Declarations Ethics Approval Not applicable. Consent to Participate Not applicable. Consent for Publication All authors give consent for publication. Conflict of Interest The authors declare no competing interests. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permis- sion directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/ licenses/by/4.0/. References Acquaah G (2012) Principles of Plant Genetics and Breeding. Wiley-Blackwell, USA Aliche EB, Oortwijn M, Theeuwen TPJM, Bachem CWB, Visser RGF, Linden CGVD (2018) Drought response in field grown potatoes and the interactions between canopy growth and yield. Agric Water Manag 206:20–30. https:// doi. org/ 10. 1016/j. agwat. 2018. 04. 013 https://doi.org/10.1007/s11540-024-09695-3 https://doi.org/10.1007/s11540-024-09695-3 http://creativecommons.org/licenses/by/4.0/ http://creativecommons.org/licenses/by/4.0/ https://doi.org/10.1016/j.agwat.2018.04.013 1622 Potato Research (2024) 67:1603–1624 1 3 Aliche EB, Theeuwen TP, Oortwijn M, Visser RG, van der Linden CG (2020) Carbon partitioning mechanisms in potato under drought stress. Plant Physiol Biochem 146:211–219. https:// doi. org/ 10. 1016/j. plaphy. 2019. 11. 019 Amjad A, Randhawa MA, Butt MS, Asghar M (2016) Assessing the processing quality of different potato cultivars during storage at various temperatures.  J Chem Soc Pak  38(5). https:// jcsp. org. pk/ Publi shedV ersion/ 2f8b4 2b6- 9c25- 4162- ad2c- 60196 6debd d8Man uscri pt% 20no% 2027,% 20Fin al% 20Gal ly% 20Pro of% 20of% 20110 99% 20(Adnan% 20Amj ad). pdf Anderson JAD, Lewthwaite SL, Genet RA, Braam WF, Gallagher DTP (1997) ‘Pacific’: a new fresh market and crisping potato cultivar. N Z J Crop Hortic Sci 25(3):303–305. https:// doi. org/ 10. 1080/ 01140 671. 1997. 95140 20 Baritelle AL, Hyde GM (2003) Specific gravity and cultivar effects on potato tuber impact sensitivity. Postharvest Biol Technol 29(3):279–286. https:// doi. org/ 10. 1016/ S0925- 5214(03) 00003-6 Burton WG (1981) Challenges for stress physiology in potato. Am Potato J 58:53–14. https:// doi. org/ 10. 1007/ BF028 55376 Canadian Food Inspection Agency (2020) Potato Varieties - ‘Snowden’. https:// www. inspe ction. gc. ca/ plant- health/ potat oes/ potat ovari eties/ ’Snowd en’/ eng/ 13125 87385 900/ 13125 87385 901. Accessed 23 May 2022 CIP (2013) Protocol for designing and conducting potato field experiments for modeling purposes. CIP, Peru Crossa J (1990) Statistical analysis of multilocation trials. Adv Agron 44:55–85. https:// doi. org/ 10. 1016/ S0065- 2113(08) 60818-4 Crossa J, Fox PN, Pfeiffer WH, Rajaram S, Gauch HG (1991) AMMI adjustment for statistical analysis of an international wheat yield trial. Theor Appl Genet 81(1):27–37. https:// doi. org/ 10. 1007/ BF002 26108 FAO (2021) FAOSTAT statistical database. https:// www. fao. org/ faost at/ en/# data. Accessed 30 Jan 2022 FAO (2022) FAOSTAT statistical database. https:// www. fao. org/ faost at/ en/# data. Accessed 23 May 2022 Gauch HG, Zobel RW (1988) Predictive and postdictive success of statistical analyses of yield trials. Theor Appl Genet 76(1):1–10. https:// doi. org/ 10. 1007/ BF002 88824 Gawronska H, Thornton MK, Dwelle RB (1992) Influence of heat stress on dry matter production and photo-assimilate partitioning by four potato clones. Am Potato J 69(10):653–665. https:// doi. org/ 10. 1007/ BF028 52678 Ghosh SC, Asanuma K, Kusutani A, Toyota M (2000) Effects of temperature at different growth stages on nonstructural carbohydrate, nitrate reductase activity, and yield of potato. Environ Control in Biol 38(4):197–206 Guoju X, Fengju Z, Zhengji Q, Yubi Y, Runyuan W, Juying H (2013) Response to climate change for potato water use efficiency in semi-arid areas of China. Agric Water Manag 127:119–123. https:// doi. org/ 10. 1016/j. agwat. 2013. 06. 004 Haase NU (2003) Estimation of dry matter and starch concentration in potatoes by determination of under-water weight and near infrared spectroscopy. Potato Res 46:117–127 Hammes PS, De Jager JA (1990) Net photosynthetic rate of potato at high temperatures. Potato Res 33(4):515–520. https:// doi. org/ 10. 1007/ BF023 58030 Hancock RD, Morris WL, Ducreux LJM, Morris JA, Usman M, Verrall SR, Fuller J, Simpson CG, Zhang R, Hedley PE, Taylor MA (2014) Physiological, biochemical, and molecular responses of the potato (Solanum tuberosum L.) plant to moderately elevated temperature. Plant, Cell Environ 37:439–450. https:// doi. org/ 10. 1111/ pce. 12168 Hastilestari BR, Lorenz J, Reid S, Hofmann J, Pscheidt D, Sonnewald U, Sonnewald S (2018) Decipher- ing source and sink responses of potato plants (Solanum tuberosum L.) to elevated temperatures. Plant, Cell Environ 41(11):2600–2616. https:// doi. org/ 10. 1111/ pce. 13366 Hijmans RJ (2003) The effect of climate change on global potato production. American Journal of Potato Res 80(4):271–279. https:// doi. org/ 10. 1007/ BF028 55363 Hiller LK, Koller DC, Thornton RE (1985) Physiological disorders of potato tubers. Potato Physiology 389–455. https:// doi. org/ 10. 1016/ B978-0- 12- 447660- 8. 50017-X Hillerm LK, Thornton RE (2008) Managing physiological disorders. In: Johnson DA (ed) Potato health management, 2nd edn. APS Press, USA Hongyu K, García-Peña M, Araújo LB, Santos D, Carlos T (2014) Statistical analysis of yield trials by AMMI analysis of genotype × environment interaction. Biometrical Letters 51(2):89–102. https:// doi. org/ 10. 2478/ bile- 2014- 0007 IPCC (2014) Climate change 2014: synthesis report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. In: Pachauri RK and Meyer LA (eds) Geneva, Switzerland https://doi.org/10.1016/j.plaphy.2019.11.019 https://doi.org/10.1016/j.plaphy.2019.11.019 https://jcsp.org.pk/PublishedVersion/2f8b42b6-9c25-4162-ad2c-601966debdd8Manuscript%20no%2027,%20Final%20Gally%20Proof%20of%2011099%20(Adnan%20Amjad).pdf https://jcsp.org.pk/PublishedVersion/2f8b42b6-9c25-4162-ad2c-601966debdd8Manuscript%20no%2027,%20Final%20Gally%20Proof%20of%2011099%20(Adnan%20Amjad).pdf https://jcsp.org.pk/PublishedVersion/2f8b42b6-9c25-4162-ad2c-601966debdd8Manuscript%20no%2027,%20Final%20Gally%20Proof%20of%2011099%20(Adnan%20Amjad).pdf https://doi.org/10.1080/01140671.1997.9514020 https://doi.org/10.1080/01140671.1997.9514020 https://doi.org/10.1016/S0925-5214(03)00003-6 https://doi.org/10.1007/BF02855376 https://doi.org/10.1007/BF02855376 https://www.inspection.gc.ca/plant-health/potatoes/potatovarieties/’Snowden’/eng/1312587385900/1312587385901 https://www.inspection.gc.ca/plant-health/potatoes/potatovarieties/’Snowden’/eng/1312587385900/1312587385901 https://doi.org/10.1016/S0065-2113(08)60818-4 https://doi.org/10.1016/S0065-2113(08)60818-4 https://doi.org/10.1007/BF00226108 https://www.fao.org/faostat/en/#data https://www.fao.org/faostat/en/#data https://doi.org/10.1007/BF00288824 https://doi.org/10.1007/BF02852678 https://doi.org/10.1007/BF02852678 https://doi.org/10.1016/j.agwat.2013.06.004 https://doi.org/10.1016/j.agwat.2013.06.004 https://doi.org/10.1007/BF02358030 https://doi.org/10.1111/pce.12168 https://doi.org/10.1111/pce.13366 https://doi.org/10.1007/BF02855363 https://doi.org/10.1016/B978-0-12-447660-8.50017-X https://doi.org/10.2478/bile-2014-0007 https://doi.org/10.2478/bile-2014-0007 1623 1 3 Potato Research (2024) 67:1603–1624 Jefferies RA (1993) Cultivar responses to water stress in potato: effects of shoot and roots. New Phytol 123(3):491–498. https:// doi. org/ 10. 1111/j. 1469- 8137. 1993. tb037 61.x Kalbarczyk R (2004) The relation between agrometeorological factors and the potato crop yields in dif- ferent regions of Poland. Acta Agrophys 4(2):339–350 Kempton RA (1984) The use of biplots in interpreting variety by environment interactions. The J Agric Sci 103(1):123–135. https:// doi. org/ 10. 1017/ S0021 85960 00433 92 Lafta AM, Lorenzen JH (1995) Effect of high temperature on plant growth and carbohydrate metabolism in potato. Plant Physiol 109:637–643. https:// doi. org/ 10. 1104/ pp. 109.2. 637 Lahlou O, Ouattar S, Ledent JF (2003) The effect of drought and cultivar on growth parameters, yield, and yield components of potato. Agronomie 23(3):257–268. https:// doi. org/ 10. 1051/ agro: 20020 89 Landcare Research (2020) Soils Portal. https:// soils. landc arere search. co. nz/. Accessed 13 Feb 2019 Levy D, Kastenbaum E, Itzhak Y (1991) Evaluation of parents and selection for heat tolerance in the early generations of a potato (Solanum tuberosum L.) breeding program. Theor. Appl. Genet 82(2):130–136. https:// doi. org/ 10. 1007/ BF002 26203 Li J, Cang Z, Jiao F, Bai X, Zhang D, Zhai R (2017) Influence of drought stress on photosynthetic charac- teristics and protective enzymes of potato at seedling stage. J Saudi Soc Agric 16(1):82–88. https:// doi. org/ 10. 1016/j. jssas. 2015. 03. 001 Lugt C, Bodlaender KBA, Goodijk G (1964) Observations on the induction of secondgrowth in potato tubers. Eur Potato J 7(4):219–227. https:// doi. org/ 10. 1007/ BF023 68253 MacKerron DKL, Heilbronn TD (1985) A method for estimating harvest indices for use in surveys of potato crops. Potato Res 28:279–282. https:// doi. org/ 10. 1007/ BF023 57582 Martin RJ, Jamieson PD, Wilson DR, Francis GS (1992) Effects of soil moisture deficits on yield and quality of ‘Russet Burbank’ potatoes. N Z J Crop Hortic Sci 20(1):1–9. https:// doi. org/ 10. 1080/ 01140 671. 1992. 10422 319 McLeod CC (1973) Potato Spacing Trials South Canterbury. Proceedings Agronomy Society of NZ 3:69–72. https:// www. agron omyso ciety. org. nz/ uploa ds/ 94803/ files/ 1973_ 11._ South_ Cante rbury_ potato_ spaci ng_ trials. pdf Merfield CN, Geary IJ, Hale RJ, Hodge S (2015) Field evaluation of the effectiveness of mesh crop cov- ers for the protection of potatoes from tomato potato psyllid. N Z J Crop Hortic Sci 43(2):123–133. https:// doi. org/ 10. 1080/ 01140 671. 2015. 10155 76 New Zealand Agricultural Greenhouse Gas Research Centre (2012) Impacts of Global Climate Change on New Zealand Agriculture. https:// www. nzagrc. org. nz/. Accessed 19 Oct 2019 NIWA (2020) Summer 2018–19. https:// niwa. co. nz/ clima te/ summa ries/ seaso nal/ summer- 2018- 19. Accessed 19 Oct 2019 Ohashi Y, Nakayama N, Saneoka H, Mohapatra PK, Fujita K (2008) Differences in the responses of stem diameter and pod thickness to drought stress during the grain filling stage in soybean plants. Acta Physiol Plant 31(2):271. https:// doi. org/ 10. 1007/ s11738- 008- 0229-4 Paget MF, Alspach PA, Genet RA, Apiolaza LA (2014) Genetic variance models for the evaluation of resistance to powdery scab (Spongospora subterranea f. sp. subterranea) from long-term potato breeding trials. Euphytica 197:369–385. https:// doi. org/ 10. 1007/ s10681- 014- 1073-9 Potato Association of America (2020) Varieties - Snowden. http:// potat oasso ciati on. org/ indus try/ varie ties/ white- varie ties/ snowd en- solan um- tuber osum. Accessed 25 May 2020 Potatoes New Zealand (2022) Potatoes New Zealand Inc. AGM Papers & Annual Report 2022 R Core Team (2017) R: A language and environment for statistical computing. R foundation for statistical computing. http:// www.R- proje ct. org/. Accessed 19 Oct 2019 Ramírez DA, Yactayo W, Gutiérrez R, Mares V, De Mendiburu F, Posadas A, Quiroz R (2014) Chloro- phyll concentration in leaves is an indicator of potato tuber yield in watershortage conditions. Sci Hortic 168:202–209. https:// doi. org/ 10. 1016/j. scien ta. 2014. 01. 036 Reynolds MP, Ewing EE, Owens TG (1990) Photosynthesis at high temperature in tuber-bearing solanum species: a comparison between accessions of contrasting heat tolerance. Plant Physiol 93(2):791– 797. https:// doi. org/ 10. 1104/ pp. 93.2. 791 Rykaczewska K (2015) The effect of high temperature occurring in subsequent stages of plant develop- ment on potato yield and tuber physiological defects. Am J Potato Res 92(3):339–349. https:// doi. org/ 10. 1007/ s12230- 015- 9436-x Rykaczewska K (2017) Impact of heat and drought stresses on size and quality of the potato yield. Plant Soil Environ 63(1):40–46. https:// doi. org/ 10. 17221/ 691/ 2016- PSE Salinger MJ, Diamond HJ, Behrens E, Fernandez D, Fitzharris BB, Herold N, Johnstone P, Kerckhoffs LHJ, Mullan AB, Parker AK, Renwick J, Scofield C, Siano AB, Smith RO, South PM, Sutton PJ, https://doi.org/10.1111/j.1469-8137.1993.tb03761.x https://doi.org/10.1017/S0021859600043392 https://doi.org/10.1104/pp.109.2.637 https://doi.org/10.1051/agro:2002089 https://soils.landcareresearch.co.nz/ https://doi.org/10.1007/BF00226203 https://doi.org/10.1016/j.jssas.2015.03.001 https://doi.org/10.1016/j.jssas.2015.03.001 https://doi.org/10.1007/BF02368253 https://doi.org/10.1007/BF02357582 https://doi.org/10.1080/01140671.1992.10422319 https://doi.org/10.1080/01140671.1992.10422319 https://www.agronomysociety.org.nz/uploads/94803/files/1973_11._South_Canterbury_potato_spacing_trials.pdf https://www.agronomysociety.org.nz/uploads/94803/files/1973_11._South_Canterbury_potato_spacing_trials.pdf https://doi.org/10.1080/01140671.2015.1015576 https://www.nzagrc.org.nz/ https://niwa.co.nz/climate/summaries/seasonal/summer-2018-19 https://doi.org/10.1007/s11738-008-0229-4 https://doi.org/10.1007/s10681-014-1073-9 http://potatoassociation.org/industry/varieties/white-varieties/snowden-solanum-tuberosum http://potatoassociation.org/industry/varieties/white-varieties/snowden-solanum-tuberosum http://www.R-project.org/ https://doi.org/10.1016/j.scienta.2014.01.036 https://doi.org/10.1104/pp.93.2.791 https://doi.org/10.1007/s12230-015-9436-x https://doi.org/10.1007/s12230-015-9436-x https://doi.org/10.17221/691/2016-PSE 1624 Potato Research (2024) 67:1603–1624 1 3 Teixeira E, Thomsen MS, Trought MCT (2020) Unparalleled coupled ocean-atmosphere summer heatwaves in the New Zealand region: drivers, mechanisms and impacts. Clim Change. https:// doi. org/ 10. 1007/ s10584- 020- 02730-5 Shahnazari A, Liu F, Andersen MN, Jacobsen SE, Jensen CR (2007) Effects of partial rootzone dry- ing on yield, tuber size and water use efficiency in potato under field conditions. Field Crops Res 100(1):117–124. https:// doi. org/ 10. 1016/j. fcr. 2006. 05. 010 Shock CC, Zalewski JC, Stieber TD, Burnett DS (1992) Impact of early-season water deficits on Russet Burbank plant development, tuber yield and quality. Am Potato J 69(12):793–803. https:// doi. org/ 10. 1007/ BF028 54186 Siano A, Roskruge N, Kerckhoffs L, Sofkova-Bobcheva S (2018) Yield and tuber quality variability in commercial potato cultivars under abiotic stress in New Zealand. Agronomy New Zealand 48:149– 163. https:// www. agron omyso ciety. org. nz/ files/ ASNZ_ 2018_ 14._ Potato_ yield_ and_ tuber_ quali ty. pdf Silveira LCID, Kist V, Paula TOMD, Barbosa MHP, Peternelli LA, Daros E (2013) AMMI analysis to evaluate the adaptability and phenotypic stability of sugarcane genotypes. Scientia Agricola 70:27– 32. https:// doi. org/ 10. 1590/ S0103- 90162 01300 01000 05 Tang R, Niu S, Zhang G, Chen G, Haroon M, Yang Q, Om PR, Li X-Q (2018) Physiological and growth responses of potato cultivars to heat stress. Botany 96(12):897–912. https:// doi. org/ 10. 1139/ cjb- 2018- 0125 Trawczyński C (2009) Wpływ nawadniania kroplowego i fertygacji na plon i wybrane elementy jakości bulw ziemniaka (Effect of drip irrigation and nitrogen fertilization on yield of root vegetables). Infrastructure and Ecology of Rural Areas 3:55–67 ((in Polish)) van Loon CD (1981) The effect of water stress on potato growth, development, and yield. Am Potato J 58(1):51–69. https:// doi. org/ 10. 1007/ BF028 55380 Wolf S, Olesinski AA, Rudich J, Marani A (1990) Effect of high temperature on photosynthesis in pota- toes. Ann Botany 65(2):179–185. https:// doi. org/ 10. 1093/ oxfor djour nals. aob. a0879 22 Zobel RW, Wright MJ, Gauch HG Jr (1988) Statistical analysis of a yield trial. Agron J 80:388–393. https:// doi. org/ 10. 2134/ agron j1988. 00021 96200 80000 30002x Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. https://doi.org/10.1007/s10584-020-02730-5 https://doi.org/10.1007/s10584-020-02730-5 https://doi.org/10.1016/j.fcr.2006.05.010 https://doi.org/10.1007/BF02854186 https://doi.org/10.1007/BF02854186 https://www.agronomysociety.org.nz/files/ASNZ_2018_14._Potato_yield_and_tuber_quality.pdf https://www.agronomysociety.org.nz/files/ASNZ_2018_14._Potato_yield_and_tuber_quality.pdf https://doi.org/10.1590/S0103-90162013000100005 https://doi.org/10.1139/cjb-2018-0125 https://doi.org/10.1139/cjb-2018-0125 https://doi.org/10.1007/BF02855380 https://doi.org/10.1093/oxfordjournals.aob.a087922 https://doi.org/10.2134/agronj1988.00021962008000030002x Effects of Abiotic Stress Associated with Climate Change on Potato Yield and Tuber Quality Under a Multi-environment Trial in New Zealand Abstract Introduction Materials and Methods Site Description and Local Climate Crop Management and Trial Layout Weather Data Agronomic and Morpho-physiological Data Yield and Tuber Quality Data Statistical Analysis Genotype × Environment Analysis Results Weather Data Abiotic Stress Effects on Morpho-physiological Characteristics of Potato Crops Abiotic Stress Effects on Final Yield and Tuber Quality of Potato Crops Abiotic Stress Effects on Tuber Physiological Disorders of Potato Crops Genotype × Environment Interaction Under the Multi-environment Trial of Potato Crops Discussion Weather Data Abiotic Stress Effects on Morpho-physiological Characteristics of Potato Crops Abiotic Stress Effects on Final Yield and Tuber Quality of Potato Crops Abiotic Stress Effects on Tuber Physiological Disorders of Potato Crops Genotype × Environment Interaction Under Multi-environment Trial of Potato Crops Conclusions Acknowledgements References