UniPhen “PIWI”—high-resolution simulation of the phenological development of 13 fungus-tolerant cultivars based on a broad observation data set from Central Europe
Abstract
Due to their high potential for fungicide reduction, fungus-tolerant “PIWI” cultivars are increasingly gaining interest in European viticulture. This investigation aimed (i) to obtain phenological observation datasets of PIWI cultivars across a broad set of European locations, (ii) to apply the temperature sum based “UniPhen” model allowing for a precise simulation of the phenological development at all BCCH stages between bud swelling and berries ripe for harvest and (iii) to discuss the potential implementations in applied viticulture as well as in viticultural climate change impact research.
Four consecutive years of complete data sets of phenological observations originating from eight locations in Central Europe for 13 PIWI and 3 traditional cultivars were used to apply the UniPhen model using a cumulative degree day approach with three temperature thresholds (lower threshold: 10 °C; upper threshold: 20 °C; heat threshold: 30 °C). Muscaris showed the thermal-temporally earliest budburst, while Solaris was earliest in beginning of flowering and harvest ripeness. The latest budburst and the latest beginning of flowering in PIWI cultivars were observed in Pinotin, while harvest ripeness was reached latest by Calardis blanc. The average normalised standard deviation (SD15 °C) over all stages, locations and cultivars was 5.5, corresponding to 5.5 days at 15 °C, with the lowest SD15 °C values around budburst and flowering stages. The highest SD15 °C values were observed in the bunch closure and post-veraison stages.
UniPhen “PIWI” enables a precise simulation of all 31 BBCH stages between the beginning of the bud swell (01) and berries being ripe for harvest (89) for 13 fungus-tolerant cultivars and can be extended to additional cultivars. The model can be applied (i) as a bioclimatic indicator describing the suitability of a location/region for the cultivation of specific cultivars under present and future climate conditions, (ii) for the simulation of cultivar-specific phases of highest susceptibility for fungal diseases and, indirectly, the timing of fungicide treatments as well as (iii) for the classification of the relative late frost risk depending on the thermal-temporal precocity of budburst. This knowledge helps to lower the barrier to growing PIWI cultivars and helps to pave the way for a more sustainable, climate change-resilient viticulture.
Introduction
Grapevine (Vitis vinifera L.) phenology is mainly temperature-driven (Jones, 2013). Air temperature is the main factor influencing the timing of phenological phases (e.g., Gladstones, 2011; Jones, 2013) if water and radiation requirements of the plants are fulfilled (Nendel, 2010; Webb et al., 2007). Consequently, grape phenology models in viticulture are usually based on air temperature observation data as a single input parameter (e.g., Caffara & Eccel, 2010; Cola et al., 2014; Molitor et al., 2014b; Parker et al., 2011; Parker et al., 2020).
Many phenological models assume linear forcing effects of air temperature by accumulating temperatures above a defined threshold to temperature sums or cumulative degree days (e.g., Amerine & Winkler, 1944; Duchêne et al., 2010; Nendel, 2010; Parker et al., 2011). However, in such unlimited temperature sum approaches the physiological evidence that above plant-specific threshold temperatures the forcing effect of the temperature will not increase further or even tends to decrease is neglected (Bonhomme, 2000; Molitor et al., 2014b; Wang & Engel, 1998; Yan & Hunt, 1999). Consequently, Molitor et al. (2014b) incorporated (i) an upper threshold temperature and (ii) a heat threshold temperature in their temperature sum-based model approach and demonstrated a significantly improved model accuracy compared to the uncapped cumulative degree day approaches previously published. Meanwhile, this approach has been proven to be of high accuracy in a broad range of cultivars using a unified threshold triplet of 10, 20 and 30 °C as lower, upper and heat temperature thresholds (Molitor et al., 2020). This “UniPhen” approach is covering all phenological stages according to the BBCH scale and is open for the incorporation of any other grape cultivar based on high-resolution observation data (Molitor et al., 2020).
Climate change is threatening viticultural production (van Leeuwen et al., 2024). Global near-surface mean air temperature has increased by over 1 °C over the past century, and all major climate projections are projecting a further increase by up to 3 °C by 2100 (IPCC, 2021). Moreover, global warming is expected to lead to a general increase in the frequency and intensity of extreme weather events (IPCC, 2021) such as droughts and heat waves (Fraga et al., 2020), threatening grape yield as well as wine quality and typicity.
In viticulture, weather and climate control grapevine growth, physiology, yield, and berry composition (Santos et al., 2020). Hence, climate change is a direct threat to the socio-economic sustainability of many viticultural regions (e.g., Lereboullet et al., 2013; Mosedale et al., 2016). Furthermore, a large portion of the total active ingredient volume of pesticides used in agriculture is applied in vineyards, and the use of fungicides in viticulture alone accounts for more than 70 % of the total fungicide use in Europe (Wingerter et al., 2021). Via the application of fungicides, viticulture contributes to climate change by releasing greenhouse gases during pesticide synthesis, transport, and application in the field via viticultural machinery. Addressing climate change in viticulture requires significant efforts to (i) adapt to changing environmental conditions as well as (ii) mitigate greenhouse gas emissions by promoting more sustainable forms of vineyard management.
In this context, the cultivation of new fungus-tolerant or fungus-resistant grape cultivars (Töpfer & Trapp, 2022), the so-called “PIWIs - pilzwiderstandsfähige Rebsorten”, has received increasing attention in recent years as they bear the potential to address climate change adaptation and mitigation simultaneously. PIWIs are new, cross-bred cultivars bearing resistance genes stemming from American or Asian Vitis species (Pedneault & Provost, 2016) and display reduced susceptibility towards major fungal diseases in viticulture, namely downy and powdery mildew. Consequently, PIWI cultivation allows for a substantial reduction in the use of fungicides (minus 50–80 %) (Töpfer & Trapp, 2022), which represent 96 % of pesticides used in viticulture (Marinho et al., 2020).
By increasing PIWI cultivation, the wine industry could thus actively contribute to reducing greenhouse gas emissions and ensure the sustainability of the sector under changing climatic conditions, while reducing its overall environmental footprint. Accordingly, increased PIWI cultivation may not only contribute to achieving pesticide reduction targets, but also to the EU Biodiversity strategy for 2030. Furthermore, the cultivation of PIWI cultivars represents an economical risk reduction strategy under increasingly frequent extreme weather conditions such as long-lasting rain periods in summer, causing severe calamites of fungal diseases and hampering access to the vineyards with machines.
Despite their potential benefits, PIWIs are presently cultivated only in 0.5 % of the viticultural area of Luxembourg (Anonymous, 2023) while in countries like Switzerland or Germany, the PIWI cultivars exceed 3 % of the nationwide planted area (Destatis, 2023; OFAG, 2023). In the relatively new and rapidly growing grape-growing country, Belgium, meanwhile, PIWIs represent 29.6 % of the viticultural area (Anonymous, 2024). Anyhow, there is still a lack of detailed information about the viticultural traits of these new cultivars, including the response of grape phenology to heat consumption that is key to consider the suitability of PIWI cultivars across different climatic conditions. The growing interest in PIWI cultivars demonstrated the need to integrate them into existing models. Hence, there is a need for more detailed and transregional studies on the differences in phenological development of different PIWI cultivars, especially in the context of climate change and fungal disease pressure.
The aim of the present investigation was (i) to obtain datasets of high-resolution phenological observations of PIWI cultivars (and traditional reference cultivars) at a broad range of European locations, (ii) to apply the temperature sum based UniPhen model for all cultivars of observation allowing for a precise simulation of the phenological development at all BCCH stages between beginning of bud swelling and berries ripe for harvest and (iii) to discuss the potential implementations of UniPhen “PIWI” in applied viticulture as well as in viticultural climate change impact research.
Materials and methods
1. Experimental vineyards
The phenological monitoring was carried out at 15 experimental vineyards across Europe (Figure 1). All sites are located along a topographical gradient with Ath (Belgium) and Haidegg (Austria) as the two extremes with elevations of 60 and 600 m a.s.l., respectively. All sites provide a fair representation of the climatic conditions and soil types that characterise the vineyards cultivated in Europe outside the Mediterranean region. According to the Köppen-Geiger climate distribution in Europe (Beck et al., 2018), the sites Wädenswil and Haidegg characterised the continental cold climate with dry season and warm summer (Dfb), Laimburg characterised the temperate climate with hot summer and no dry season (Cfa), and the rest of the sites the temperate climate with warm summer and no dry season (Cfb). The main soil types are represented by a Cambisol (CM) soil unit that is present in eight of the experimental vineyards (Table 1). The Luvisol (LV) is found in Ath (Belgium) and Kreuznach (Germany), and the soil of the remaining four sites are Fluvisols (FV), Phaeozems (PH), Planosols (PL), and Regosols (RG).
Systematic phenological assessments took place in the years 2021 to 2024 in the following locations (Table 1):
Location | Country | Coordinates | Years of observation | Years with complete data sets | Soil type (EUSIS, 2023) |
Remich | Luxembourg | 49.55 N, 06.36 E | 2021–2024 | 4 | Cambisol (CMvr) |
Geilweilerhof | Germany | 49.22 N, 08.05 E | 2021–2024 | 4 | Phaeozems (PHlv) |
Wädenswil | Switzerland | 47.23 N, 08.68 E | 2021–2023 | 0 | – |
Ath | Belgium | 50.62 N, 03.77 E | 2021–2024 | 4 | Luvisol (LVha) |
Marcelin | Switzerland | 46.52 N, 06.48 E | 2021–2024 | 4 | Cambisol (CMeu) |
Kindel | Germany | 49.97 N, 07.06 E | 2021–2024 | 4 | Cambisol (CMdy) |
Neustadt | Germany | 49.37 N, 08.18 E | 2021–2024 | 4 | Regosol (RGca) |
Freiburg | Germany | 47.97 N, 07.83 E | 2021–2024 | 4 | Cambisol (CMdy) |
Ebringen | Germany | 47.96 N, 07.77 E | 2021–2024 | 4 | Cambisol (CMca) |
Geisenheim | Germany | 49.98 N, 07.93 E | 2022–2024 | 3 | Fluvisol (FLeu) |
Changins | Switzerland | 46.40 N, 06.23 E | 2022–2023 | 0 | Cambisol (CMeu) |
Bernkastel | Germany | 49.92 N, 07.04 E | 2023 | 0 | Cambisol (CMdy) |
Laimburg | Italy | 46.36 N, 11.28 E | 2023 | 1 | Cambisol (CMeu) |
Haidegg | Austria | 46.63 N, 15.50 E | 2023–2024 | 2 | Planosol (PLeu) |
Kreuznach | Germany | 49.86 N, 07.84 E | 2023 | 1 | Luvisol (LVha) |
Only those locations providing four years of complete data sets were selected for model calibration.
Phenology of the following cultivars was recorded: Cabernet blanc, Cabernet Cortis, Cabertin, Calardis blanc, Johanniter, Monarch, Muscaris, Pinotin, Regent, Sauvignac, Solaris, Souvignier gris, Villaris, Divico, Helios, Bronner, Cabaret noir, Cabernet Carbon, Divona, Rondo, Hibernal, Prior, Baron, Floreal, Satin noir, Vidoc, Artaban, Laurot, VB 1-22, Donauriesling, Müller-Thurgau, Pinot noir, Riesling, Chardonnay, Chasselas and Dornfelder. Only those location × cultivar combinations with complete data sets for all four years of observation, leading to at least 8 complete data sets (covering all BBCH stages between 01 and 89) per cultivar, were used as input data for model calibration (Table 2).
Cultivar | Observation data sets | Complete data sets |
Cabernet blanc | 25 | 20 |
Cabernet Cortis | 21 | 16 |
Cabertin | 13 | 12 |
Calardis blanc | 25 | 16 |
Johanniter | 17 | 12 |
Monarch | 12 | 8 |
Muscaris | 25 | 16 |
Pinotin | 14 | 12 |
Regent | 19 | 8 |
Sauvignac | 26 | 20 |
Solaris | 30 | 24 |
Souvignier gris | 24 | 20 |
Divico | 18 | 12 |
Müller-Thurgau | 30 | 24 |
Pinot noir | 42 | 32 |
Riesling | 24 | 20 |
In total, calibration data sets consisted of the following complete data sets (covering all BBCH stages between 01 and 89 in all four years of observation) originating from the following locations (Table 3).
Location | Cabernet blanc | Cabernet Cortis | Cabertin | Calardis blanc | Johanniter | Monarch | Muscaris | Pinotin | Regent | Sauvignac | Solaris | Souvignier gris | Divico | Müller-Thurgau | Pinot noir | Riesling |
Remich | X | X | X | X | X | X | X | X | X | X | X | X | X | |||
Geilweilerhof | X | X | X | X | X | X | X | X | X | X | X | X | X | X | ||
Ath | X | X | X | X | X | X | X | X | X | X | ||||||
Marcelin | X | X | ||||||||||||||
Kindel | X | X | X | X | X | X | ||||||||||
Neustadt | X | X | X | X | X | X | X | X | ||||||||
Freiburg | X | X | X | X | X | X | X | X | ||||||||
Ebringen | X | X | X | X | X | X | X |
All cultivars were trained in a cane-pruned vertical shoot positioning system (VSP). Six plants were monitored per cultivar; the plants of observation were the same in all four years.
2. Assessment of phenological stages
All phenological growth stages according to the BBCH (Biologische Bundesanstalt, Bundessortenamt und Chemische Industrie) scale as defined by Lorenz et al. (1995) (Table S1) between BBCH 01 (beginning of bud swelling) and BBCH 85 (softening of berries) were recorded when 50 % of the vines or shoots exhibited the respective stage. The assessment intervals ranged from two to three days. Usually, records were taken at the same location by the same person in all years and cultivars. The date of BBCH 89 was defined as the DOY (day of the year) when 60° Oechsle (= 14.17 °Brix) was reached or exceeded for the first time in the respective year following the definition of Molitor et al. (2020). Maturity control took place at weekly intervals between veraison (BBCH 81) and harvest by collecting approximately 50 randomly selected berries (clusters from different positions of the canopy; berries from different positions in the cluster) per cultivar (avoiding berries with visible bunch rot symptoms). The date of reaching a specific phenological stage was noted as day of the year (DOY).
3. Meteorological measurements
Daily average air temperatures at 2 m height were recorded during the period of examination by standard automatic weather stations in proximity (less than 200 m) of the experimental vineyards. Average growing season (April to October) temperatures per location and year are given in Table S2.
The highest growing season temperatures were observed in the year 2022 (average 17.3 °C), which was on average 2.5 °C warmer than in 2021 (average 14.8 °C). On average, Ebringen was the location with the highest average growing season temperatures (17.5 °C), while the lowest average growing season temperatures were recorded in Ath (15.4 °C).
4. Application of the unified high-resolution phenological model approach “UniPhen”
In the present study, the approach of the high-resolution cumulative degree day-based phenological model, as previously developed for the Rivaner cultivar (Molitor et al., 2014b), parameterised for Riesling (Molitor et al., 2016) and Pinot noir (Molitor & Junk, 2019), as well as unified by using a unified global optimised temperature threshold triplet (lower threshold: 10 °C, heat threshold: 20 °C, heat threshold: 30 °C) (Molitor et al., 2020) was applied to all selected cultivars.
In this “UniPhen” approach, degree days (DD) for a specific day are calculated based on the daily average air temperatures, applying the following conditions:
DD10,20,30 = 0, if t ≤ 10 °C or if t ≥ 40 °C
DD10,20,30 = (t – 10), if 10 °C ≤ t ≤ 20 °C
DD10,20,30 = 10, if 20 °C ≤ t ≤ 30 °C
DD10,20,30 = 10 – (t – 30), if 30 °C ≤ t ≤ 40 °C
where DD10,20,30 is the degree day for a specific day and t is the daily average temperature.
For every growth stage, degree days (DD) were summed up to cumulative degree days (CDDs) using the Equation (1):
(1)
where CDD10,20,30 are the cumulative degree days relative to the date of BBCH 09 in Riesling; DD10,20,30 the value of the degree day for the day i; m the DOY when BBCH 09 was observed in Riesling; and n the DOY when the target phenological stage was reached in the respective cultivar.
The reference date was the day of the year (DOY) at which BBCH 09 was reached in Riesling (following Molitor et al. (2020)).
If no observation data for BBCH 09 for Riesling were available (Marcelin, Freiburg, Ebringen), its virtual DOY was calculated based on the observed DOY of BBCH 09 for Pinot noir. According to the data presented by Molitor et al. (2020), Pinot noir reaches BBCH 09 6.1 CDD10,20,30 prior to Riesling. Hence, the date of BBCH 09 in Riesling was fixed at the DOY at which a temperature sum of 6.1 CDD10,20,30 after to the DOY of observation of BBCH 09 in Pinot noir was reached or exceeded.
5. Data analysis
To determine general temperature sum thresholds to simulate all phenological stages between BBCH 01 and BBCH 89 in all cultivars, average CDD10,20,30 were calculated for each cultivar, relative to the date of reaching BBCH 09 in Riesling. Standard deviations of the cumulative degree days were calculated for each phenological stage in each cultivar and each location. To normalise the relative value of the standard deviations caused by the daily DD values, standard deviations were divided by the theoretical DD at 15 °C (approximate growing season temperature in Remich in the period of observation, selected following Molitor et al. (2020)). These normalised standard deviations, SD15° are given for all phenological stages, all cultivars and all locations.
The R script used to calculate the CDD per cultivar is available in Jiménez-Rodríguez (2024) as UNIPHEN-Tool v1.0.0.
To assess the goodness of fit of the model, mean bias errors (MBE) and mean absolute errors (MAE) were calculated for all years, locations and cultivars as an average of all BBCH stages. A positive MBE value indicates that, on average, the model overestimates actual observations; a negative value indicates that the model underestimates actual observations under specific annual or local conditions (Caffara & Eccel, 2010).
Results
1. Phenological observation data
Days of the year (DOY) on which the phenological stages BBCH 01 to 89 were reached in the years 2021 to 2024 at the selected locations (complete data sets over four years) in the 16 cultivars of observation are shown in Table S3.
2. Cultivar specific temperature sum thresholds for the different phenological stages
Figure 2 shows the average CDD10,20,30 threshold values relative to BBCH 09 in Riesling for the different cultivars. On average, the lowest CDD10,20,30 value for BBCH 09 was recorded for Muscaris (–13) and the highest for Pinotin (3). The lowest CDD10,20,30 value for BBCH 61 was observed for Solaris (195). Divico showed the lowest CDD10,20,30 value for BBCH 81 (601), Solaris for BBCH 89 (797). For all those three stages highest values were observed for Riesling (BBCH 61: 251, BBCH 81: 859, BBCH 89: 1012) (Figure 2). Figure 3 displays graphically the differences between the different cultivars in CDD10,20,30 values for the respective stages.
Lowest SD15 °C values (except BBCH 09 for Riesling as reference) were observed in case of all stage × cultivar combinations for Muscaris at BBCH 11 (0.6 = 0.6 days at 15 °C) and highest for Riesling at BBCH 79 (19.2). On average of all cultivars, lowest SD15 °C values were recorded for BBCH 09 (1.2) and highest for BBCH 79 (13.4). On average of all stages, the lowest SD15 °C values were observed in Sauvignac (5.0) and the highest in Divico (6.7), with 5.5 being the average SD15 °C value over all 16 cultivars (Figure 4).
Highest MBE values (recording later than simulated) were recorded in Cabernet Cortis in Geilweilerhof (20.1 CDD10,20,30) and lowest (recording earlier than simulated) in Pinot noir in Marcelin (–21.4 CDD10,20,30). On average of all cultivars highest MBEs were observed in case of Geilweilerhof (7.6 CDD10,20,30) and lowest in case of Freiburg (–12.7 CDD10,20,30).
Highest MAE values were observed for Pinot noir in Marcelin (23.1 CDD10,20,30) and lowest for Sauvignac in Neustadt (7.0 CDD10,20,30). Likewise, on average of all cultivars highest MAE values were reached in Marcelin (17.6 CDD10,20,30) and lowest in Neustadt (10.2 CDD10,20,30) (Table 4).
MBE | MAE | |||||||||||||||
Cultivar | Re | Gf | At | Ma | Ki | Nw | Fr | Eb | Re | Gf | At | Ma | Ki | Nw | Fr | Eb |
Cabernet blanc | 1.3 | 1.1 | –4.8 | 9.1 | –6.7 | 11.8 | 8.8 | 9.2 | 10.6 | 15.7 | ||||||
Cabernet Cortis | 20.1 | –10.2 | –16.7 | 6.8 | 20.6 | 14.2 | 16.7 | 11.8 | ||||||||
Cabertin | –11.3 | 10.2 | 1.1 | 14.3 | 13.3 | 8.2 | ||||||||||
Calardis blanc | –6.3 | 10.9 | 10.9 | –15.4 | 14.3 | 14.5 | 13.0 | 15.4 | ||||||||
Johanniter | –3.9 | 3.5 | 0.4 | 7.9 | 10.0 | 9.5 | ||||||||||
Monarch | 6.8 | –6.8 | 8.9 | 8.9 | ||||||||||||
Muscaris | –0.3 | 4.8 | –6.3 | 1.9 | 7.9 | 12.8 | 11.7 | 10.0 | ||||||||
Pinotin | –8.4 | 5.1 | 3.3 | 14.3 | 9.4 | 10.8 | ||||||||||
Regent | 7.8 | –7.8 | 10.4 | 10.4 | ||||||||||||
Sauvignac | –9.2 | 4.7 | –6.6 | 2.9 | 8.2 | 15.2 | 7.7 | 9.4 | 7.0 | 15.8 | ||||||
Solaris | 2.0 | 10.7 | –10.8 | 2.0 | –14.9 | 11.1 | 9.3 | 13.2 | 11.3 | 10.6 | 15.1 | 12.2 | ||||
Souvignier gris | 6.4 | 9.5 | –3.4 | –18.4 | 5.9 | 13.3 | 14.9 | 9.3 | 18.6 | 13.6 | ||||||
Divico | 7.3 | –3.2 | –4.0 | 13.4 | 12.0 | 12.8 | ||||||||||
Müller-Thurgau | 5.3 | 5.0 | 4.7 | –11.0 | –10.9 | 6.8 | 12.4 | 12.5 | 13.5 | 12.9 | 12.4 | 11.1 | ||||
Pinot noir | –0.2 | 10.3 | 5.6 | –21.4 | –6.0 | 12.6 | –8.7 | 7.7 | 15.1 | 12.8 | 14.2 | 23.1 | 9.3 | 13.0 | 16.1 | 11.5 |
Riesling | 1.1 | 6.0 | 8.5 | –14.4 | –1.2 | 8.0 | 12.6 | 18.4 | 15.6 | 8.2 | ||||||
Average | –1.2 | 7.6 | –1.3 | –12.3 | –7.7 | 4.7 | –12.7 | 5.7 | 12.1 | 12.0 | 12.7 | 17.6 | 11.0 | 10.2 | 15.3 | 13.1 |
Concerning the vintage effect, highest MBE values were recorded in Muscaris in 2024 (12.9 CDD10,20,30) and lowest in Regent in 2021 (–23.0 CDD10,20,30). On average of all cultivars, highest MBE were observed in 2023 (7.3 CDD10,20,30) and lowest in 2021 (–11.9 CDD10,20,30). In 2021, the MBE values of all cultivars were negative.
Highest MAE values were observed for Regent in 2023 (25.1) and lowest for Calardis blanc in 2022 (7.0 CDD10,20,30). On average of all cultivars highest MAE values were recorded for 2021 (16.3 CDD10,20,30) and lowest for 2022 (9.7 CDD10,20,30) (Table 5).
MBE | MAE | |||||||
Cultivar | 2021 | 2022 | 2023 | 2024 | 2021 | 2022 | 2023 | 2024 |
Cabernet blanc | –14.2 | 0.0 | 5.1 | 9.1 | 18.6 | 7.9 | 9.3 | 12.6 |
Cabernet Cortis | –15.7 | 1.0 | 9.2 | 5.5 | 18.5 | 9.0 | 11.8 | 11.5 |
Cabertin | –10.8 | 3.4 | 3.7 | 3.7 | 19.7 | 11.6 | 7.8 | 11.8 |
Calardis blanc | –5.8 | –2.7 | 11.0 | –2.6 | 12.1 | 7.0 | 12.1 | 7.6 |
Johanniter | –11.8 | 1.8 | 12.5 | –2.5 | 15.2 | 9.0 | 12.6 | 11.7 |
Monarch | –14.4 | 5.1 | 7.2 | 2.2 | 18.3 | 11.3 | 8.1 | 12.6 |
Muscaris | –15.4 | –4.0 | 6.6 | 12.9 | 18.4 | 8.4 | 8.3 | 16.1 |
Pinotin | –13.3 | –1.0 | 13.5 | 0.8 | 19.8 | 8.2 | 14.4 | 12.4 |
Regent | –23.0 | 11.3 | 8.4 | 3.3 | 25.1 | 12.4 | 13.8 | 15.8 |
Sauvignac | –8.3 | –8.1 | 8.2 | 8.1 | 11.4 | 10.6 | 12.8 | 11.8 |
Solaris | –7.3 | –7.4 | 6.6 | 8.1 | 12.0 | 9.3 | 7.8 | 12.4 |
Souvignier gris | –5.3 | –2.6 | –4.6 | 12.5 | 8.7 | 8.5 | 12.3 | 18.6 |
Divico | –3.3 | –10.4 | 8.0 | 5.7 | 12.3 | 15.0 | 16.1 | 10.1 |
Müller-Thurgau | –14.3 | 2.8 | 7.8 | 3.7 | 16.8 | 8.7 | 9.6 | 8.9 |
Pinot noir | –12.6 | –4.0 | 8.4 | 8.2 | 15.4 | 8.6 | 11.6 | 11.5 |
Riesling | –15.7 | 2.2 | 4.8 | 8.8 | 17.9 | 9.1 | 8.0 | 12.4 |
Average | –11.9 | 0.8 | 7.3 | 5.5 | 16.3 | 9.7 | 11.0 | 12.4 |
Discussion
1. Model application for PIWI cultivars
The application of the UniPhen for 13 PIWI cultivars based on at least eight observation data sets revealed the average CDD10,20,30 thresholds relative to BBCH 09 in Riesling until the respective BBCH stage is reached per cultivar. It allows for a classification of the relative precocity of all 13 PIWI cultivars (plus three traditional reference cultivars) in each of the 31 stages of the BBCH scale between 01 and 89. As described before (Molitor et al., 2020), the relative precocity of the different cultivars is not stable over the different phenological stages, for example, Calardis blanc is one of the latest cultivars in leaf development, while one of the earliest in beginning of flowering and again the latest PIWI cultivar to reach harvest maturity (Figure 2). According to their precocity in budburst, Muscaris and Solaris could be classified as “early” while Pinotin and Cabertin are rather late, even slightly later than the traditional reference Riesling (consequences are discussed in the paragraph Practical applications).
The present investigations confirmed the general suitability of the cumulative degree day approach using three threshold temperatures as proposed first by Molitor et al. (2014b) and unified later by the introduction of the optimised temperature sum threshold triplet 10, 20 and 30 °C in the UniPhen approach (Molitor et al., 2020). Observed global average normalised standard deviations correspond to 5.5 days at 15 °C (or 2.75 days at 20 °C) demonstrating the strong robustness of the UniPhen model also for the PIWI cultivars.
Generally, for all cultivars, normalised standard deviations SD15 °C are lowest around budburst. In this period normalised standard deviations correspond to 1 to 2 days, indicating the high model precision. While in stages of inflorescence emergence standard deviations slightly increase, during flowering (BBCH 61-69) and fruit-set (BBCH 71) the SD15 °C values correspond to standard deviations of less than 6 days. As described before (Molitor et al., 2020) for traditional cultivars, higher average normalised standard deviations were observed in the bunch closure stages 77 and 79. Here, as well as in the stages of inflorescence emergence the perception of the observer is more subjective than in case of leaf emergence or flowering (Verdugo-Vasquez et al., 2017). Furthermore, the degree of fruit set, water and nutrient availability, as well as the yield potential have an influence on the cluster architecture and, in consequence, the moment of fruit set. Standard deviations corresponding to 13.3 days at 15 °C were observed for stage BBCH 89. These standard deviations can partly be explained by the observation intervals (7 days (weekly maturity control) versus 2–3 days for the other stages) as well as annual differences in sugar accumulation caused by differences in crop load (Santesteban & Royo, 2006).
The effects of the locations as well as the vintages were evaluated using the mean bias error and the mean average error. While in Geilweilerhof on average of all years and cultivars the stages were recorded 7.6 CDD10,20,30 later than simulated, for Freiburg this was the case 12.7 CDD10,20,30 earlier (corresponding to 2 to 3 days at 15 °C). These minor effects might be explained by the differences in the perception of the observers or potentially micro-climatic effects caused by the position of the weather station relative to the vineyard of observation. Observed average mean average errors at the different locations between 10.2 CDD10,20,30 in Neustadt and 17.6 CDD10,20,30 in Marcelin might be reasoned by differences in the observation interval or the consistency of the observations.
Interestingly, concerning the deviations caused by the vintage, all cultivars showed in 2021 a negative mean bias error, meaning that all stages were on average earlier observed than simulated, with deviations corresponding to an average of 2 days at 15 °C. While this deviation of 2 days might be acceptable in the application of the model, this systematicity might have its origin in the consistently lower temperatures in 2021 compared to the two following years and/or non-linear dependencies between the temperatures recorded by the weather stations and the apparent temperatures in the vineyard canopy.
2. Practical applications
Generally, the UniPhen approach is designed to allow the simulation of phenological phases, where environmental conditions or crop cultural measures, for example, have a distinct influence on factors such as yield formation (Molitor & Keller, 2016), wine typicity (Molitor & Junk, 2019) or susceptibility towards pests or diseases (Molitor et al., 2016; Molitor et al., 2020).
Usually, PIWI cultivars are not completely resistant against fungal infections, but of (i) lower susceptibility and/or (ii) their period of susceptibility of specific organs is shorter than in traditional cultivars since ontogenetic resistance appears phenologically earlier. Frequently, a limited number of fungicide applications in PIWI cultivars is recommended in the period of highest susceptibility for berry infections (Schumacher et al., 2024). Periods of highest susceptibility, however, differ in their length between different PIWI cultivars, for example, Cabernet Cortis and Solaris are susceptible towards Plasmopara viticola until BBCH 71-73, whereas Souvignier gris is still susceptible at BBCH 75 (Schumacher et al., 2024). Hence, simulating the phenological development supports a more targeted timing of fungicide treatments in PIWIs depending on their phenology driven status of susceptibility. The knowledge about the phenological status of different PIWI cultivars potentially contributes to a further reduction of pesticide use if applications outside the period of susceptibility are abandoned, for example, some PIWI cultivars might need to be treated at a specific point in time while other cultivars passed already their cultivar specific phenological growth stage of ontogenetic resistance. Here, the combination of the knowledge of cultivar-specific phases of susceptibility with UniPhen “PIWI” might deliver valuable decision support for users.
The obtained classification of relative cultivar precocity at different stages opens new chances for cultivar selection under practical conditions, especially in regions where late frost damage might constitute a serious threat for sustainable wine production (e.g., Kartschall et al., 2015; Leolini et al., 2018; Molitor et al., 2014a; Mosedale et al., 2015), such as in many cool-climate viticulture regions including the new winegrowing regions emerging in recent years. Especially in these regions, PIWIs are frequently the cultivars of choice due to their reduced susceptibility as well as the missing customer’s habit to prefer traditional cultivars. Under high late frost risk conditions cultivars with early bud development and early bud burst might be avoided based on the information derived from the present model.
Solaris and Muscaris were identified as thermal-temporally early in budburst and, hence, more likely to be damaged by late frost events at dates where other cultivars are not yet susceptible, since green shoots and leaves have not yet emerged. Taking into consideration a long-term average (1991–2020) April to May (period where late frost damage takes place) temperature in Remich of 12.4 °C, the observed difference of 10 (Solaris) or even 13 (Muscaris) CDD10,20,30 compared to Riesling could be translated to a precocity of four days for budburst and, hence, a four or five days longer risk period for late frost damage. This fact should be considered when thinking about cultivar selection in regions or locations where late frost damage is likely. Concerning the late frost damage risk, Pinotin or Cabertin might be interesting candidates since they are the only PIWI cultivars with a slightly later budburst than the traditional reference cultivar Riesling. Indeed, under cool climate conditions such as in the newly emerging Northern European grape growing region, the ideotype of a PIWI would be one that exhibits a late/delayed bud burst combined with an early full maturity. In fact, the observations as well as the methodology of the present paper could also be useful for breeders in selecting parental material for the breeding of new cultivars with those ideal characteristics.
Generally, PIWI cultivars could potentially pave the way for long-term adaptation strategies to achieve higher climate resilience (Terleth & Tavernar, 2022; Tscholl et al., 2024). Consequently, the percentage of PIWI cultivars in new plantations is increasing in recent years – especially, but not exclusively, in the new winegrowing regions and at higher altitudes, which have in past not been suitable for viticulture due to limited heat consumption. However, so far little was known about the cultivar specific heat demand of PIWI cultivars to produce grapes with adequate maturity for wines of high quality. The BBCH scale defined the stage BBCH 89 as “berries ripe for harvest”. Hence, the temperature sum necessary to reach BBCH 89 is linked to the annual heat demand for a specific cultivar to produce fully mature grapes. Using UniPhen “PIWI” as a bioclimatic indicator allows to describe the suitability of different locations or altitudes for the cultivation of specific cultivars under changing climatic conditions.
Especially the classification of relative precocity of the maturity period (BBCH 81-89) might be used for the selection of PIWI cultivars in climate change adaptation strategies. Particularly high temperatures in the maturation period have been demonstrated to negatively affect wine fruitiness and its aroma (Duchêne et al., 2010). Consequently, PIWI cultivars with a thermal-temporally early ripening period leading to hot temperatures during maturation such as Solaris should be avoided in the traditional Central European wine growing regions – especially since a twofold increase of the ripening period temperatures is expected with ongoing climate change (Molitor & Junk, 2019). On the other hand, early ripening PIWI cultivars are of highest interest for cool winegrowing regions outside the traditional regions as well for viticulture at higher altitudes. Here, it is recommended to check prior to the planting of a vineyard if the annual heat consumption of a region or location fits to the heat demand of the cultivar of interest. The information about the general suitability of cultivars for cultivation under specific climatic conditions can be derived from UniPhen “PIWI”. Based on long-term observation data, (i) virtual dates of reaching specific phenological stages as well as CDD10,20,30 values at the end of the season can be calculated and (ii) potentially fitting cultivars might be selected, for example, if the CDD10,20,30 value at the end of the season reaches 800, only Solaris might be a suitable candidate out of the 13 PIWI cultivars tested here. Since the cultivar selections are long-term decisions, additionally, future climate projections might be incorporated to simulate the future heat consumption for a specific region.
Generally, the observed systematic biases caused by the location were low and might be caused by differences in the perception of the observers. The observed mean bias errors between –12.7 and 7.6 CDD10,20,30 refer to less than 3 days earlier or less than 2 days later observations than simulated, respectively. Consequently, the model proved to be valid under present Central European conditions. Under completely different climatic conditions, the model may need to be validated first with local observation data before being applied. However, the chance that cultivars with heat demands corresponding to the conditions in Central Europe (as the present ones), might be cultivated under clearly hotter conditions appears unlikely. Generally, it must be taken into consideration that extreme weather events, like late frost damage or hailstorm damage, destroying the green plant tissues, might lead to a temporary stop or even a reset of the phenological development. In such cases the model might not be applicable.
Conclusions
Based on a broad set of observation data, the UniPhen model was applied and allows now for a precise simulation of all 31 BBCH stages between the beginning of bud swell (01) and berries ripe for harvest (89) for 13 PIWI cultivars. UniPhen “PIWI” could be used (i) to simulate periods of highest susceptibility towards fungal diseases and contribute to reduce fungicide use to a minimum as well as (ii) a bioclimatic indicator for the suitability of a location/region for the cultivation of specific PIWI cultivars, considering late frost risk and minimum heat demand to produce high quality wines. Present data and the model itself might help to alleviate lowering the existing barriers for grape growers (due to a lack of experience and knowledge) to grow these alternatives to traditional, susceptible cultivars.
UniPhen “PIWI” as well as the UniPhen approach in general is open to be extended for other grape cultivars and application under different or future climatic conditions.
Acknowledgements
The authors thank M. Schultz, D. Dam, C. Simon, P. Lopes (IVV), T. Kaltenbach, C. Mertes (WBI Freiburg), A.-K. Ertel (JKI Geilweilerhof), L. Kuenzler (Agroscope Wädenswil) for supporting the phenological assessments.
This work was partly funded by the Luxembourgish Ministry of Agriculture, Food and Viticulture in the framework of the research projects “VinoManAOP”, “VinoManAOP2” and “PIWI3”, as well as by funds of the Federal Ministry of Food and Agriculture (BMEL) based on a decision of the parliament of the Federal Republic of Germany via the Federal Office for Agriculture and Food (BLE) under the Federal Programme for Ecological Farming in the framework of the research project “VITIFIT”.
References
- Amerine, M. A., & Winkler, A. J. (1944). Composition and quality of musts and wines of California grapes. Hilgardia, 15(6), 493–675. https://doi.org/10.3733/hilg.v15n06p493
- Anonymous (2023). Das Weinjahr 2022 und seine Ernteergebnisse. Veröffentlichung der Abteilung Weinbau des Weinbauinstitutes in Remich.
- Anonymous (2024). La viticulture belge en chiffres. SPF Economie, P.M.E.
- Beck, H. E., Zimmermann, N. E., McVicar, T. R., Vergopolan, N., Berg, A., & Wood, E. F. (2018). Present and future Koppen-Geiger climate classification maps at 1-km resolution. SCIENTIFIC DATA, 5, 180214. https://doi.org/10.1038/sdata.2018.214
- Bonhomme, R. (2000). Bases and limits to using ‘degree.day’ units. European Journal of Agronomy, 13(1), 1–10. https://doi.org/10.1016/S1161-0301(00)00058-7
- Caffara, A., & Eccel, E. (2010). Increasing the robustness of phenological models for Vitis vinifera cv. Chardonnay. International Journal of Biometeorology, 54(3), 255–267. https://doi.org/10.1007/s00484-009-0277-5
- Cola, G., Mariani, L., Salinari, F., Civardi, S., Bernizzoni, F., Gatti, M., & Poni, S. (2014). Description and testing of a weather-based model for predicting phenology, canopy development and source-sink balance in Vitis vinifera L. cv. Barbera. Agricultural and Forest Meteorology, 184, 117–136. https://doi.org/10.1016/j.agrformet.2013.09.008
- Destatis (2023). Mit Keltertrauben bestockte Rebfläche: Deutschland, Jahre, Rebsorte [Dataset]. https://www-genesis.destatis.de/genesis/online?operation=table&code=41252-0001&bypass=true&levelindex=0&levelid=1728310746247#abreadcrumb
- Duchêne, E., Huard, F., Dumas, V., Schneider, C., & Merdinoglu, D. (2010). The challenge of adapting grapevine varieties to climate change. Climate Research, 41(3), 193–204. https://doi.org/10.3354/cr00850
- EEA (2016). Elevation map of Europe based on GTOPO30. Prod-ID: DAT-36-en [Map]. European Environment Agency (EEA). https://www.eea.europa.eu/en/datahub/datahubitem-view/357f686f-1939-4652-8e79-c9ac7a1c5da6
- EUSIS (2023). Soil Geographical Database of Eurasia at Scale 1:1,000,000. European Soil Database v2.0 (vector and attribute data) European Soil Data Centre (ESDAC). [Dataset].
- Fraga, H., Molitor, D., Leolini, L., Santos, J. A, (2020). What is the impact of heatwaves on European viticulture? A modelling assessment. Applied Sciences, 10(9), 3030. https://doi.org/10.3390/app10093030
- Gladstones, J. (2011). Wine, Terroir and Climate Change. Wakefield Press. https://doi.org/10.1063/pt.5.010145
- IPCC (2021). The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (Camebridge University Press, p. 2391). Cambridge. https://www.ipcc.ch/report/sixth-assessment-report-working-group-i/
- Jiménez-Rodríguez, C. D. (2024). Vine-Tools: a set of simple tools developed with R for the analysis of vineyard data (1.0.0). Zenodo. https://doi.org/10.5281/zenodo.14408480
- Jones, G. V. (2013). Winegrape phenology. In M. D. Schwartz (Ed.), Phenology: An integrative environmental science. Second edition, pp. 563–580. Springer. https://doi.org/10.1007/978-94-007-6925-0_30
- Kartschall, T., Wodinski, M., von Bloh, W., Oesterle, H., Rachimow, C., & Hoppmann, D. (2015). Changes in phenology and frost risks of Vitis vinifera (cv Riesling). Meteorologische Zeitschrift, 24(2), 189–200. https://doi.org/10.1127/metz/2015/0534
- Leolini, L., Moriondo, M., Fila, G., Costafreda-Aumedes, S., Ferrise, R., & Bindi, M. (2018). Late spring frost impacts on future grapevine distribution in Europe. Field Crops Research, 222, 197–208. https://doi.org/10.1016/j.fcr.2017.11.018
- Lereboullet, A.-L., Beltrando, G., & Bardsley, D. K. (2013). Socio-ecological adaptation to climate change: A comparative case study from the Mediterranean wine industry in France and Australia. Agriculture, Ecosystems & Environment, 164, 273–285. https://doi.org/10.1016/j.agee.2012.10.008
- Lorenz, D.H., Eichhorn, K.W., Bleiholder, H., Klose, R., Meier, U. and Weber, E. (1995). Phenological growth stages of the grapevine, Vitis vinifera L. ssp. vinifera. Codes and descriptions according to the extended BBCH scale. Australian Journal of Grape and Wine Research, 1, 100-103. https://doi.org/10.1111/j.1755-0238.1995.tb00085.x
- Marinho, M. da C., Salazar Diogo, B., Lage, O. M., & Antunes, S. C. (2020). Ecotoxicological evaluation of fungicides used in viticulture in non-target organisms. Environmental Science and Pollution Research, 27(35), 43958–43969. https://doi.org/10.1007/s11356-020-10245-w
- Molitor, D., Caffara, A., Sinigoj, P., Pertot, I., Hoffmann, L., & Junk, J. (2014a). Late frost damage risk for viticulture under future climate conditions: A case study for the Luxembourgish winegrowing region. Australian Journal of Grape and Wine Research, 20(1), 160–168. https://doi.org/10.1111/ajgw.12059
- Molitor, D., Junk, J., Evers, D., Hoffmann, L., & Beyer, M. (2014b). A high-resolution cumulative degree day-based model to simulate phenological development of grapevine. American Journal of Enology and Viticulture, 66(1), 72–80. https://doi.org/10.5344/ajev.2013.13066
- Molitor, D., & Keller, M. (2016). Yield of Müller-Thurgau and Riesling grapevines is altered by meteorological conditions in the current and previous growing seasons. OENO One, 50(4), 245–258. https://doi.org/10.20870/oeno-one.2016.50.4.1071
- Molitor, D., Baus, O., Hoffmann, L., & Beyer, M. (2016). Meteorological conditions determine the thermal-temporal position of the annual Botrytis bunch rot epidemic on Vitis vinifera L. cv. Riesling grapes. OENO One, 50(4), 231–244. https://doi.org/10.20870/oeno-one.2016.50.4.36
- Molitor, D., & Junk, J. (2019). Climate change is implicating a two-fold impact on air temperature increase in the ripening period under the conditions of the Luxembourgish grapegrowing region. OENO One, 53, 409–422. https://doi.org/10.20870/oeno-one.2019.53.3.2329
- Molitor, D., Fraga H., & Junk, J. (2020). UniPhen –a unified high resolution model approach to simulate the phenological development of a broad range of grape cultivars as well as a potential new bioclimatic indicator. Agricultural and Forest Meteorology, 291. https://doi.org/10.1016/j.agrformet.2020.108024
- Mosedale, J. R., Wilson, R. J., & Maclean, I. M. D. (2015). Climate change and crop exposure to adverse weather: Changes to frost risk and grapevine flowering conditions. PLoS ONE, 10(10), e0141218. https://doi.org/10.1371/journal.pone.0141218
- Mosedale, J. R., Abernethy, K. E., Smart, R. E., Wilson, R. J., & Maclean, I. M. D. (2016). Climate change impacts and adaptive strategies: Lessons from the grapevine. Global Change Biology, 22(11), 3814–3828. https://doi.org/10.1111/gcb.13406
- Nendel, C. (2010). Grapevine bud break prediction for cool winter climates. International Journal of Biometeorology¨, 54(3), 231–241. https://doi.org/10.1007/s00484-009-0274-8
- OFAG (2023). Rapport agricole 2023. Office fédéral de l`agriculture. https://www.agrarbericht.ch/fr/services/documentation/publications
- Parker, A. K., de Cortázar-Atauri, I. G., van Leeuwen, C., & Chuine, I. (2011). General phenological model to characterise the timing of flowering and veraison of Vitis vinifera L. Australian Journal of Grape and Wine Research, 17(2), 206–216. https://doi.org/10.1111/j.1755-0238.2011.00140.x
- Parker, A. K., de Cortázar-Atauri, I. G., Gény, L., Spring, L., Destac, A., Schultz, H. R., Molitor, D., Lacombe, T., Graça, A., Monamy, C., Stoll, M., Storchi, P., Trought, M. C. T., Hofmann, R. W., & van Leeuwen, C. (2020). Temperature-based grapevine sugar ripeness modelling for a wide range of Vitis vinifera L. cultivars. Agricultural and Forest Meteorology, 285. https://doi.org/10.1016/j.agrformet.2020.107902
- Pedneault, K., & Provost, C. (2016). Fungus resistant grape varieties as a suitable alternative for organic wine production: Benefits, limits, and challenges. Scientia Horticulturae, 208, 57–77. https://doi.org/10.1016/j.scienta.2016.03.016
- Santesteban, L. G., & Royo, J. B. (2006). Water status, leaf area and fruit load influence on berry weight and sugar accumulation of cv. ‘Tempranillo’ under semiarid conditions. Scientia Horticulturae, 109(1), 60–65. https://doi.org/10.1016/j.scienta.2006.03.003
- Santos, J. A., Fraga, H., Malheiro, A. C., Moutinho-Pereia, J., Dinis, L. T., Correia, C., Moriondo, M., Leolini, L., Dibari, C., Costafreda-Aumedes, S., Kartschall, T., Menz, C., Molitor, D., Junk, J., Beyer, M., & Schultz, H. R. (2020). A review of the potential climate change impacts and adaptation options for European viticulture. Applied Sciences, 10(9). https://doi.org/10.3390/app10093092
- Schumacher, S., Mertes, C., Kaltenbach, T., Bleyer, G., & Fuchs, R. (2024). A method for phenotypic evaluation of grapevine resistance in relation to phenological development. Scientific Reports, 14(1), 915. https://doi.org/10.1038/s41598-023-50666-4
- Terleth, J., & Tavernar, A. (2022). Prüfung von PIWI Sorten am Versuchszentrum Laimburg. Obstbau Weinbau - Fachmagazin Des Südtiroler Beratungsringes, 59(9), 44–48.
- Töpfer, R., & Trapp, O. (2022). A cool climate perspective on grapevine breeding: Climate change and sustainability are driving forces for changing varieties in a traditional market. Theoretical and Applied Genetics, 135(11), 3947–3960. https://doi.org/10.1007/s00122-022-04077-0
- Tscholl, S., Candiago, S., Marsoner, T., Fraga, H., Giupponi, C., & Egarter Vigl, L. (2024). Climate resilience of European wine regions. Nature Communications, 15 (1), 1-12 https://doi.org/10.1038/s41467-024-50549-w
- van Leeuwen, C., Sgubin, G., Bois, B., Ollat, N., Swingedouw, D., Zito, S., & Gambetta, G. A. (2024). Climate change impacts and adaptations of wine production. Nature Reviews Earth & Environment, 5(4), 258–275. https://doi.org/10.1038/s43017-024-00521-5
- Verdugo-Vasquez, N., Panitrur-De la Fuente, C., & Ortega-Farias, S. (2017). Model development to predict phenological scales of table grapes (cvs. Thompson, Crimson and Superior Seedless and Red Globe) using growing degree days. OENO One, 51(3), 277–288. https://doi.org/10.20870/oeno-one.2017.51.2.1833
- Wang, E. L., & Engel, T. (1998). Simulation of phenological development of wheat crops. Agricultural Systems, 58(1), 1–24. https://doi.org/10.1016/S0308-521X(98)00028-6
- Webb, L. B., Whetton, P. H., & Barlow, E. W. R. (2007). Modelled impact of future climate change on the phenology of winegrapes in Australia. Australian Journal of Grape and Wine Research, 13(3), 165–175. https://doi.org/10.1111/j.1755-0238.2007.tb00247.x
- Wingerter, C., Eisenmann, B., Weber, P., Dry, I., & Bogs, J. (2021). Grapevine Rpv3-, Rpv10- and Rpv12-mediated defense responses against Plasmopara viticola and the impact of their deployment on fungicide use in viticulture. BMC Plant Biology, 21(1). https://doi.org/10.1186/s12870-021-03228-7
- Yan, W. K., & Hunt, L. A. (1999). An equation for modelling the temperature response of plants using only the cardinal temperatures. Annals of Botany, 84(5), 607–614. https://doi.org/10.1006/anbo.1999.0955
Attachments
Supplementary data
Download