Projections of Climate Change Impacts on Flowering-Veraison Water Deﬁcits for Riesling and Müller-Thurgau in Germany

: With global warming, grapevine is expected to be increasingly exposed to water deﬁcits occurring at various development stages. In this study, we aimed to investigate the potential impacts of projected climate change on water deﬁcits from the ﬂowering to veraison period for two main white wine cultivars (Riesling and Müller-Thurgau) in Germany. A process-based soil-crop model adapted for grapevine was utilized to simulate the ﬂowering-veraison crop water stress indicator (CWSI) of these two varieties between 1976–2005 (baseline) and 2041–2070 (future period) based on a suite of bias-adjusted regional climate model (RCM) simulations under RCP4.5 and RCP8.5. Our evaluation indicates that the model can capture the early-ripening (Müller-Thurgau) and late-ripening (Riesling) traits, with a mean bias of prediction of ≤ 2 days and a well-reproduced inter-annual variability for more than 60 years. Under climate projections, the ﬂowering stage is advanced by 10–20 days (higher in RCP8.5) between the two varieties, whereas a slightly stronger advancement is found for Müller-Thurgau than for Riesling for the veraison stage. As a result, the ﬂowering-veraison phenophase is mostly shortened for Müller-Thurgau, whereas it is extended by up to two weeks for Riesling in cool and high-elevation areas. The length of phenophase plays an important role in projected changes of ﬂowering-veraison mean temperature and precipitation. The late-ripening trait of Riesling makes it more exposed to increased summer temperature (mainly in August), resulting in a higher mean temperature increase for Riesling (1.5–2.5 ◦ C) than for Müller-Thurgau (1–2 ◦ C). As a result, an overall increased CWSI by up to 15% (ensemble median) is obtained for both varieties, whereas the upper (95th) percentile of simulations shows a strong signal of increased water deﬁcit by up to 30%, mostly in the current winegrowing regions. Intensiﬁed water deﬁcit stress can represent a major threat for high-quality white wine production, as only mild water deﬁcits are


Introduction
Grapevine is a major fruit crop of economic and cultural importance worldwide.Sustainable development of the wine grape industry is of particular importance for Europe, where annual wine production accounted for about 63% of the world's total wine production in 2020 [1].However, the observed trend in climate warming [2] has already altered complex climate-soil-variety interactions.Projected climate change is expected to differently affect the grapevine agroecosystems in distinct agro-environmental contexts [3][4][5][6][7].In the Mediterranean basin, indicated as a potential hotspot, projected increases in drought and heat stress events in the future will likely result in lower wine grape production and may thus threaten the suitability for viticulture [8][9][10].Conversely, wine regions in central Europe with cooler growing seasons may become progressively viable for grapevine cultivation due to increased temperature [8,11].In the meantime, the increasing temperature can lead to an increase in the frequency and duration of drought episodes [2].Even in cool-climate wine regions, vine growth is already subject to some degree of water deficit [7,[12][13][14].
Water deficit stress is generally the focus of studies carried out in Mediterranean regions with warm and dry climates [15][16][17] but is less frequently considered for regions with temperate climates, such as in Germany.Conventionally, water availability does not represent a major challenge for traditional German winegrowing regions because the summer temperature is moderate, with a nearly uniform distribution of monthly precipitation throughout the year.However, in recent years, drought stress in summer has increasingly gained attention in many German winegrowing regions [18].Increased drought stress mainly results from the effects of (i) increased temperatures, especially in the post-flowering phenophase, leading to a higher potential evapotranspiration rate [19], and (ii) increased chance of occurrence of prolonged heat waves [20,21].The number of days per year with intense drought stress for grapevine is projected to increase in the future for shallow-soil, steep-slope vineyards in the Rheingau region [13].However, this issue has not yet been investigated for more humid conditions and/or different soil types with varying soil water holding capacities in the other German regions.Therefore, there is a clear need to extend the studied wine regions in Germany for drought-risk assessment in order to support the development of appropriate adaptation strategies under climate change to become more climate-resilient and remain economically viable.
Grapevine water status is one of the key factors determining grapevine yield and wine quality [7, 14,[22][23][24].There is a consensus in the literature that water deficit impacts vine yields because it reduces berry size [22,25,26] and bud fertility [27].Lower yields, when not compensated by higher selling prices for wine, reduce the profitability of wine production.However, mild water deficits are required for the production of high-quality wines, especially for red wine production, as it promotes the accumulation of skin phenolics (e.g., anthocyanins) in grape berries [7, 23,24].Nonetheless, severe water deficits are detrimental for both red and white wine quality [7,14,[22][23][24].On the other hand, the effects of water deficit vary according to timing, intensity and duration during successive grapevine development stages.It is widely acknowledged that pre-veraison water deficits are prone to affect berry size and composition [15,16,23,24,28,29].One reason could be the changes in hydraulic connections between vine and berry after veraison; pre-veraison berries are connected to the vine through xylem and are thus highly sensitive to the vine's water status and prone to drought-induced shriveling, whereas post-veraison berries are mainly connected to the phloem and become less sensitive to drought [16,24].The floweringveraison phenophase represents an important pre-veraison period for berry response to water deficits [15,17,30,31].Yang et al. [30] showed that flowering-veraison water availabil-ity can be critical to potential cluster weight at harvest for both temperate European and southern Mediterranean wine regions, for which the sensitivity of such relationship tends to vary annually, depending on the site and genotypes.
There have been many studies devoted to field experiments and trials to evaluate the effects of seasonal water deficits at different timings and intensities on grapevine growth performance, yield and fruit composition [15][16][17]24].Overall, these studies suggest the effects of water deficits vary according to varieties, soil characteristics, weather conditions and in-situ vineyard management (e.g., irrigation timing and amount).Soil-crop system models play a prominent role in agri-environmental research, as they can capture the complex interactions among genotypes, management and environment [32][33][34].In particular, these models can consider the magnitude, timing and duration of weather effects at subsequent development stages for different cultivars, integrating the influence of soil characteristics and field practices [32,34,35].These features make them particularly suitable for assessing phase-dependent drought risk at the field and regional scale [13,30,36].However, these models are not always well adapted to vineyards due to the inherent complexity of the viticultural system compared to that of annual crops [4].One soil-crop system model that adapts well to grapevine systems is STICS, which has demonstrated its ability to predict field-measured soil water content, phenology, biomass production and yield components for grapevine [37][38][39][40].Furthermore, STICS was shown to be able to accurately reproduce observed leaf water potential at pre-dawn (Ψpd) and midday (Ψmid) [38,39].However, Ψpd and Ψmid measurements for many years are often not available at a large geographic scale, which is a scale compatible for agricultural policy decision making.Therefore, to better support agricultural drought monitoring and assessment, the appropriate selection of drought stress index to model is crucial.Typically, the drought stress index is composed of single or multiple meteorological variables and/or hydrologic variables [41].In particular, a soil-moisture-based drought index integrated with evapotranspiration allows for an approximation of water deficits and linking to crop growth potential [41].By using the crop model, the crop growth response to soil moisture deficits can be further integrated into the drought index, considering the weather-soil-crop interactions [33,34].One such index that already applies to different crop models is the crop water stress index (CWSI), taking into account variety-specific sensitivity in crop drought response [30,36].
To assess potential climate change impacts, a previous study implemented STICS at a regional level by coupling to downscaled outputs of climate models to estimate likely effects of climate change on grapevine water deficits during the reproductive phase (fruit settingmaturity) across different European countries [42].However, water deficits during the flowering-veraison period have not been specifically investigated under projected climates.In this study, we aim to assess the impacts of projected changes in climate, as simulated by multiple state-of-the-art regional climate models (RCMs), on vine water deficits during the flowering-veraison period based on CWSI for one early-ripening (Müller-Thurgau) and one late-ripening (Riesling) white variety in Germany, as modelled by STICS.The spatial variability of projected changes in the intensity of water deficits is of particular focus.Based on the framework developed by Fraga et al. [42], we further enlarged the ensemble size of RCMs with a novel bias-adjustment method [43] and integrated a recent soil hydraulic dataset [44] for increasingly realistic soil-crop-water balance simulations.A long time series (>60 years) of observed phenology data was used to calibrate the model for local varieties.

Viticulture in Germany
In Germany, major viticulture regions are currently located in the southwest of the country [11,45].However, our analysis was carried out throughout Germany, assuming that the suitable viticulture regions can be extended as a result of climate change [8,9].For example, Neumann and Matzarakis [11] projected expanded areas suitable for growing both Müller-Thurgau and Riesling based on the Huglin index.The topography of Germany shows distinct elevations between the north (flat lowlands) and south (plateaus and mountains) (Figure S1).The predominant training system in Germany is vertical shoot positioning (VSP), usually with one or two canes per plant and planting densities of approximately 4000 plants per ha.However, in recent years, training systems such as the semi-minimal pruned hedge are increasingly gaining importance [46,47].The most common rootstocks are V. berlandieri × V. riparia siblings (e.g., SO4, 5BB, 125 AA).The majority of vineyards are non-irrigated, as annual precipitation is sufficient in most years/regions.Yields are restricted by legal regulations in marketable hl/ha of must/wine, depending on the region, cultivar and desired wine quality.

Characterization of Müller-Thurgau and Riesling
Müller-Thurgau is characterized as a high-yielding and early-ripening variety, with low berry sugar concentration, even at full ripeness [45].Riesling is regarded as a moderately high-yielding and late-ripening variety [48], with moderately low sugar concentration at full ripeness (yet higher than Müller-Thurgau) [45].When comparing the phenology the two varieties, the precocity for budbreak and flowering is similar (Müller-Thurgau is, on average, 1-3 days earlier), but Müller-Thurgau reaches veraison 10-15 days earlier than Riesling [45,49].Müller-Thurgau (syn.Rivaner in Luxembourg) and Riesling are presently the two most widely cultivated grape varieties in Germany.However, vineyard areas with Riesling have slightly increased worldwide and have been stable in Germany, whereas plantations of Müller-Thürgau sharply declined both worldwide and in Germany (Table 1).In 2010, Müller-Thurgau and Riesling accounted for 13% of and 22% of winegrowing areas in Germany, respectively [50].The impact assessment was performed between a 30-year historical baseline (1976-2005) and a 30-year future period (2041-2070).Several high-resolution (~12.5 km) regional climate models (RCMs) (Table 2) produced within the EURO-CORDEX initiative [51] were bias-adjusted (detailed in the subsequent sections) to provide the required weather variables at a daily timestep.We utilized 4 different RCMs driven by 6 global climate models (GCMs), resulting in a total of 13 unique GCM-RCM combinations (Table 2).The required meteorological variables (for impact models) were daily minimum and maximum air temperature ( • C), solar radiation (MJ m −2 day −1 ), precipitation (mm), wind speed (m s −1 ) and actual vapour pressure (mbar).These variables were obtained from RCM simulations for both the baseline and the future period.For future climate projections, a set of four representative concentration pathway (RCP) scenarios were commonly used, namely RCP2.6,RCP4.5, RCP6 and RCP8.5 [52,53].RCP2.6 denoted a mitigation scenario with an ambitious target to reduce greenhouse gas (GHG) emissions over the 21st century; RCP4.5 represented an intermediate-level emission scenario with emissions peaking around 2040 and declining afterward; RCP6 corresponded to a relatively high-emission scenario compared to RCP4.5, with emissions peaking around 2080 and declining afterward; RCP8.5 represented a business-as-usual high-emission scenario, assuming the continuous rise of GHG emissions throughout the 21st century [52,53].Since RCP2.6 is a very stringent pathway, we only selected RCP4.5 and RCP8.5 in this study to encapsulate the possible range of future evolutions.The atmospheric CO 2 concentration (ppm) was also a required crop model input, which was derived from the recommended global annual average by the Coupled Model Intercomparison Project 5 (CMIP5) for the selected historical period and the respective future period of each RCP [54].

Bias-Adjustment Method
Simulations from high-resolution RCMs can still have important biases [43,55,56].This precludes their direct usage in most climate impact studies.A common solution is to apply a suitable bias-adjustment method to reduce biases.Herein, we applied a stateof-the-art parametric quantile mapping method, namely the ISIMIP3BASD method [43].ISIMIP3BASD was specifically designed to reduce bias over the tail of the distribution, allowing for a robust adjustment of bias at extreme values [43] (Table 2).The method is also trend-preserving across all quantiles, i.e., not adding an artificial trend during bias adjustment [43,57].The reference dataset used for bias adjustment was provided by the German Meteorological Service.

Uncertainty Estimations
A multi-model ensemble of RCMs was utilized to estimate the range of possible future climate change impacts concerning different change signals projected by each RCM.The 5th to 95th percentile interval was chosen to represent the range of 'very likely' changes, according to the IPCC likelihood definition [2].However, care needed to be taken when calculating these ensemble statistics, since the RCM ensemble utilization was imbalanced, e.g., some RCMs appeared multiple times in the ensemble due to multiple driving GCMs.To account for this imbalance, statistics of simulated impacts under all GCMs of a particular RCM were averaged, ending up with one average for each RCM.Then, the median, 5th and 95th percentile were computed based on these averages.It should be noted that our uncertainty estimations did not take into account the temporal variability and the variability resulting from different driving GCMs.

Description of STICS-Grapevine Modules
STICS model followed a modular structure, with each module corresponding to a process: phenology development, LAI and root growth, radiation interception, biomass conversion, yield formation and water and nitrogen balance, among others.Herein, we only briefly described STICS (v9.1) modules concerning grapevine phenology and soil-crop water balance that were most relevant for this study.A more detailed presentation can be found in Brisson et al. [35,[58][59][60], Brisson and Perrier [61], García de Cortázar-Atauri [37] and García de Cortázar-Atauri et al. [62].

Related Phenology Modules
Accurate phenology simulation represents a critical step to quantify water deficits during a particular phenophase.The model simulated phenology development as a function of temperature accumulation.The phenophase between dormancy and budburst (BBCH09) was simulated using the BRIN model [62], taking into account the genotypedependent chilling-and forcing-requirement for dormancy-dormancy break phase and dormancy break-budburst phase, respectively.The subsequent phenology stages were all simulated using the growing degree day (GDD, • C degree.day −1 ) model with a base temperature of 10 • C [35].The required thermal threshold in GDD from budburst to fruit set, flowering (BBCH65) to fruit set, fruit setting to berry softening (BBCH85) and fruit setting to maturation were all defined as cultivar-specific parameters in the model [35].Although fruit maturity was mainly simulated based on the GDD concept, the simulation of harvest date was based on the choice of desired wine quality (e.g., berry water content) [63].

Related Soil-Crop Water Balance Modules
Soil-crop water balance corresponded to an essential component of soil-crop system modelling, providing the estimations of water deficits experienced by the plants.Water balance considers the initial soil water content, water input (precipitation + irrigation, if any), surface runoff, soil drainage and actual crop evapotranspiration (ET) [35].In particular, ET simulation was subject to the influence of all the other aforementioned factors and was of particular importance in this study.The model divided the ET computation into actual soil evaporation (E) and plant transpiration (T).E was computed by a semiempirical model: soil evaporation is maximum (E max ) until a cumulative threshold is reached; afterward, E max is reduced depending on the weather and soil type [61].T computations required first solving the term of maximum transpiration (T max ), which represented the plant water requirement.T max was computed using the Shuttleworth and Wallace (S-W) model, which characterizes the soil-plant-atmosphere system with a resistance network [64].The S-W model incorporates the influence of both physiological (e.g., canopy height and LAI) and meteorological variables (e.g., temperature, radiation, wind speed and vapour pressure) and proves to be effective in explaining the energy budget [65].Subsequently, T/T max ratio is simulated as a bilinear function of the available soil water content in the root zone [35].

Incorporation of Soil Data
Besides meteorological data, soil information is also critical for the model-based impact assessment.Soil data were grouped into two categories: soil physical characteristics and soil hydraulic properties.The former mainly consisted of soil texture, particle size distribution, surface dry albedo and bulk density, which were all derived from the Harmonized World Soil Database (HWSD, v1.2) at about 1 km resolution [66].Soil surface slope was also considered to account for the complex terrain conditions in Germany, and was estimated from the European Digital Elevation Model (EU-DEM, v1.1) at 25 m spatial resolution (source: Copernicus service).Soil hydraulic properties, including soil field capacity and wilting point, were obtained from the EU-SoilHydroGrids 1 km (v1.0), which provided recent EU-wide estimations for soil depth up to 2 m [44].To harmonize the input datasets with different spatial resolutions, the simulation grid was interpolated (nearest-neighbor) to a resolution consistent with that of RCM output (~12.5 km).Computation workload was reduced by considering only the dominant soil type within each grid, thereby neglecting potentially larger spatial variability attributed to soil variability within each grid box.Note that there were empty grids (in white color) in the results, which were either identified as non-land areas or areas without soil information according to HWSD data.More detailed information on how raw soil information was extracted and converted into model-ready inputs can be found in Yang et al. [67][68][69][70].

Model Applications 2.5.1. Model Calibrations for Müller-Thurgau and Riesling
For calibration, a 64-year-long time series (1956-2019) of observed data (BBCH09,65,81) was collected for Müller-Thurgau and Riesling at Eltviller Sonnenberg (lon 8.126 • , lat 50.038 • , elevation 165 m), Rheingau, Germany.The model was calibrated to simulate the budburst (BBCH09), flowering (BBCH65) and veraison (BBCH81) stage of the two studied white wine varieties.However, the model only simulated BBCH85 instead of BBCH81.Therefore, BBCH81 was estimated by linear interpolations between BBCH65 and BBCH85 (Section 2.3.1)using obtained GDD values at these two stages.The calibration attempted to minimize the mean squared error (MSE) [71][72][73] between observations and simulations for the three successive stages using the grid-search method.Note that although BBCH09 was not related to the water deficit assessment, the model was still calibrated for this stage in order to ensure that a good fit for BBCH65 and BBCH81 was not at the expense of poor fit for an earlier stage, such as BBCH09.Indeed, calibration followed a consecutive scheme, thus with more emphasis on the fit for BBCH09.More details on the statistical assumptions and implementation of the calibration are provided in Yang et al. [74].The resulting calibrated parameters are shown in Table S1, which successfully characterize Müller-Thurgau and Riesling as early-and late-ripening varieties, respectively.The 64-year time series of observed and simulated data is shown in Figure 1 and reveals a low bias in mean (1-2 days difference) with similar standard deviations between observations and simulations.Further evaluations for Müller-Thurgau and Riesling indicated that 65-86% and 72-86% of the observed variance, respectively, can be explained by predictions, whereas predicted stages were generally within one-week difference (MAE and RMSE) of observed stages for both varieties (Figure 1).

Simulations of Flowering-Veraison CWSI
In the field of remote sensing, the crop water stress indicator (CWSI) has been frequently applied to assess the vine water status by calculating the difference between canopy and air temperature [75,76].However, canopy temperature data for vineyards spanning many years and a large geographic scale remains scarce.Instead, several studies have reformulated CWSI to characterize the crop water deficit as the balance between crop water demand and water consumption, which can be simulated by many different biophysical models [30,36,77].CWSI was computed as follows: where ET and ET max (sum of E max and T max ) are daily model outputs (see Section 2.3.2).The CWSI ranged from 0 (no stress) to 1 (severe stress), aiming to quantify the extent to which the actual ET failed to meet the potential water consumption (i.e., ET max ) from the soil-crop system, thus estimating the magnitude of water deficits experienced by plants.Simulations of CWSI enabled the integration of the effects of soil, climate and crop characteristics.CWSI was computed as the mean over the flowering-veraison phenophase (hereafter CWSI) in order to preserve its physical meaning, as suggested by Zhu et al. [36].

Software Environment
A Python-based (v3.7) simulation interface was established at the High-Performance Cluster at PIK (HPC@PIK) for grid-based crop model simulations.The interface automated dataset preprocessing (including any arbitrary soil and climate datasets) and model setup and read and wrote the simulation outputs for each target grid point.Furthermore, it used the Python message-passing interface (MPI) to parallelize the gridded simulations, thus considerably facilitating large-scale crop model applications, especially when multiple climate model projections were utilized.For data processing and result analysis, the involved python libraries were mainly numpy (v1.21.2), panda (v1.3.4) and xarray (0.19.0).To plot the results, matplotlib (v3.4.3) and cartopy (v0.20.1) were used.

Projected Seasonal (April-October) Temperature and Precipitation Changes
The mean seasonal average temperature over the baseline period  shows a range between 13 • C and 16 • C, except for some high-elevation and mountainous regions and an area located in the extreme north of Germany, where these temperatures are generally below 13 • C (Figure 2a).The mean seasonal precipitation sum of the baseline period shows a range between 300 and 500 mm for most of the area, with only some high-elevation areas in the south and central regions exceeding 500 mm (Figure 2b).For climate change projections over 2041-2070, the mean seasonal temperature is projected (by ensemble median) to mainly increase by 1.5-2 • C under RCP4.5 and by 1.5-2.5 • C (a stronger increase for the southern parts of Germany) under RCP8.5 (Figure 2a).Seasonal precipitation under RCP4.5 is projected to have a marginal decrease of up to 6% for the southwestern part of the country (currently the major viticulture regions) and a slight increase of up to 5% for the remaining areas (Figure 2b).However, under RCP8.5, mean seasonal precipitation is projected to experience a widespread increase (mainly ≤10%)(Figure 2b).

Projected Phenology (Flowering & Veraison) Changes
The projected (temporal) mean changes in the flowering stage (BBCH65), veraison stage (BBCH81) and the phenophase (BBCH65-BBCH81) are shown in Figure 3.For Müller-Thurgau, a widespread spatial pattern of advancement of 10-15 days and 15-20 days is projected for BBCH65 and BBCH81, respectively, under RCP4.5 (Figure 3a).For the same variety under RCP8.5, a prevalent pattern of earlier occurrence by 10-20 days and by 20-25 days is projected for BBH65 and BBCH81, respectively (Figure 3a).As a result, a very similar pattern of shortened BBCH65-BBCH81 phenophase is found between the two scenarios for Müller-Thurgau; mean phenophase length is mostly reduced by about one week (≤7 days) (Figure 3a).For Riesling, nearly the same spatial pattern as that for Müller-Thurgau is projected for BBCH65, i.e., a widespread advancement by 10-15 days under RCP4.5 and by 10-20 days under RCP8.5 (Figure 3b).Meanwhile, BBCH81 is projected to have a more heterogeneous spatial pattern: for some areas in the north, central and southern high-elevation regions, the advancement is limited to ≤10 days for both scenarios, whereas the remaining areas tend to have earlier occurrence by 10-20 days in RCP4.5 and by 10-25 days in RCP8.5 (Figure 3b).The resulting mean BBCH65-BBCH81 phenophase tends to be extended by up to 14 days (>14 days for some areas in RCP8.5) in those identified areas with a small advanced veraison stage and shortened by within a week (<7 days) for the remaining area (Figure 3b).Although the ensemble median of projected change is shown here, a notable spread (5th and 95th percentile) among RCMs was discovered and is presented in Figures S2-S4.

Projected Flowering-Veraison Temperature and Precipitation Changes
The flowering-veraison mean temperature for Müller-Thurgau is projected to mainly increase by 1-2 • C in both scenarios, despite the fact that the projected higher temperature increases by 2-2.5 • C in some central and southern plateaus and mountainous regions are more pronounced in RCP8.5 than in RCP4.5 (Figure 4a).For Riesling, the mean floweringveraison temperature is projected to increase by 1.5-2.5 • C for the majority of areas under both scenarios, whereas, similarly to Müller-Thurgau, the projected temperature increase is higher (2-2.5 • C) in central and southern high-elevation regions (more pronounced in RCP8.5 than in RCP4.5) (Figure 4a).Concerning the flowering-veraison precipitation sum, there is a dominant spatial pattern of mean decrease by up to 15% for Müller-Thurgau under both scenarios (Figure 4b).In contrast, the distribution is more heterogeneous for Riesling.The north and central regions mainly show increases of up to 15%, whereas the rest of the study area tends to show a decrease of up to 15% under RCP4.5 (Figure 4b).In RCP8.5, a stronger (than RCP4.5) signal of increased precipitation was discovered; mean increases of 15% are predominant, whereas there are more areas with higher precipitation increases of 15-45% (Figure 4b).The ensemble median projections are shown here, while a marked spread is also found that is particularly pronounced for the projected precipitation change (Figures S5 and S6).These projected changes can be partially explained by how the shifted phenophase is exposed to projected monthly mean temperature and precipitation changes, as shown for both RCP4.5 (Figure 5a) and RCP8.5 (Figure 5b).With a higher magnitude of changes under RCP8.5, the projected monthly average temperature increase (by ensemble median) is relatively higher in summer months (July-September) than in spring months (April-June): 1.7-2.3• C, 1.9-2.5 • C and 2.3-2.7 • C as averaged over the north, central and south regions, respectively, for summer months, as compared to those with <1.8 • C in spring months (Figure 5b).The respective (90% uncertainty range) spread of monthly temperature increase among RCMs is generally larger in April, June, September and October (up to 1.4 • C) and much smaller in July and August (<0.5 • C) (Figure 5b).For the precipitation changes (by ensemble median) in RCP8.5, the mean increases of monthly precipitation occur in April-June and in October, i.e., 8-14% in the north, 5-12% in the central and southern regions, whereas the changes in July-September are more uncertain, as smaller mean changes with a larger spread are generally found in these months than in the other months (Figure 5b).The flowering-veraison phenophase of both varieties is shifted earlier, but the two varieties tend to have different exposures to projected monthly meteorological changes.For Müller-Thurgau, the shifted phenophase is accompanied by reduced length of the phase (Figure 3a), spanning from mid-June to early August for all three regions, thus being substantially exposed to summer months with increased mean temperature and uncertain precipitation changes (although this may partially capture the mean increased June precipitation) (Figure 5b).For Riesling, the shifted phenophase mainly spans from mid-June to late August but is accompanied by varying lengths of phenophase across all regions (Figure 3b).Thus, Riesling is more exposed to summer months as compared to Müller-Thurgau (Figure 5b).A similar pattern is also projected for RCP4.5 (Figure 5a) but with an overall smaller magnitude of changes as compared to RCP8.5 (Figure 5b).In the line plot, the lines represent regional mean values with the respective colour band indicating their 90% uncertainty range among RCMs.In the bar plot, bars indicate the regional mean values with the respective error bars indicating their 90% uncertainty range among RCMs.The regional mean flowering stage and veraison stage are also shown for Müller-Thurgau (triangular) and Riesling (diamond).Note that for the same variety in each subplot, the phenology symbols occurring earlier correspond to the flowering stage, and the next symbol is the veraison stage.The phenology symbols with black, yellow and red are for the baseline, RCP4.5 and RCP8.5, respectively.

Projected Impacts of Climate Change on Flowering-Veraison CWSI
The mean flowering-veraison CWSI is projected to have a similar pattern for the two varieties, in which a mean increase of up to 15% is projected for most of the area (including the major winegrowing regions) in both scenarios, except for the northeastern and central areas, which are projected to decrease by up to 15% (Figure 6).However, there is a more consistent increase (by ensemble median) of CWSI in RCP4.5 than in RCP8.5, i.e., RCP8.5 projections show more areas (mainly in the northeastern and central regions) with decreased CWSI (Figure 6).However, these results are subject to considerable uncertainties, as the lower (5th) percentile simulations of RCMs (both varieties) tend to show widespread mean reductions of CWSI of up to 15% in RCP8.5, with a more heterogeneous pattern for RCP4.5 (Figure 6).Conversely, the higher (95th) percentile simulations of RCMs (both varieties) largely show a consistent spatial pattern of mean increased CWSI of up to 30% in both scenarios (Figure 6).The individual model projections are provided in Appendix A. Although both varieties tend to have a similar spatial pattern for projected CWSI changes, they differ in their projections of ET max and ET (CWSI components).For Müller-Thurgau, both the sum of ET max and that of ET over the flowering-veraison phase tend to show widespread reductions by up to 15% (relative to baseline), as indicated by both the ensemble median and 5th/95th percentile (Figures S7 and S8).In contrast, for Riesling, both ET max and ET sums primarily show an increase by up to 15% (some areas >30%) over most of Germany, whereas some areas show a decrease by up to 15% (Figures S7 and S8).

Evaluation of Model Performance for Phenology Simulations
Accurate phenology simulations are critical to assess the water deficits calculated over the flowering-veraison phase.First, the model was calibrated to simulate the flowering (BBCH65) and veraison (BBCH81) stages using a long time series (1956-2019) of observed data for the two varieties.Satisfactory model performance was obtained for both varieties, with a mean bias of only 1-2 days and well-reproduced interannual variability (Figure 1).When the model was further driven by gridded historical RCM simulations, it is corroborated that the simulated regional mean DOY of flowering and of veraison (Figure 5) agree well with local observed mean timings (mean over 1968-2010) on 22th June and 11th August (mean 50 day interval) for Müller-Thurgau and on 24th June and 24th August (mean 61 day interval) for Riesling [78].

Phenology Projections
With increased temperatures under climate change projections, the results indicate generally advanced phenology stages, although a greater advancement of phenology was discovered for RCP8.5 than for RCP4.5 (Figure 3).A similar magnitude of advancement (10-20 days across scenarios) was found for the flowering stage between varieties, whereas a slightly greater advancement of veraison is projected for Müller-Thurgau (15-25 days across scenarios) than for Riesling (10-25 days across scenarios) (Figure 3).In comparison, Molitor and Junk [19] projected advanced flowering and veraison stages by 9 and 10 days for Müller-Thurgau and by 9 and 11 days for Riesling, respectively, over 2031-2060 under the A1B (comparable to RCP6.0) emission scenario in Luxembourg.The relatively higher magnitude of advancement in our study can be attributed to projected higher seasonal mean temperature increases by 1.5-2.5 • C (depending on the scenario) compared with <1.5 • C projected by Molitor and Junk [19].Many studies suggest that differences in the magnitude of earlier occurrence of simulated phenology stages under global warming can be explained by differences in the studied varieties, the phenology models applied and the extent of temperature increases [19,31,45,79].Moreover, smaller advanced veraison stages (≤ 10 days) of Riesling (under both scenarios) are identified in areas with either a cooler reference climate (Figure 2a) or mountainous terrain (>400 m) (Figure S1).A particularly higher temperature increase is projected in these areas (Figure 4a).Hence, there could have been a delay in the veraison date due to projected high temperatures (mainly in August for Riesling) that frequently exceed the defined upper threshold (37 • C) above which phenology development is halted [35].As a result of such asymmetric phenology advancements (a greater advancement for flowering than for veraison stage), the floweringveraison phenophase in those areas is extended by up to 14 days (some >14 days under RCP8.5)(Figure 3b).In Alsace (France), delayed veraison projections with increased temperature were also obtained for Riesling using a similar GDD model and with a 33 • C upper threshold as compared with projections without a threshold [79].Similarly, Ramos et al. [31] also report delayed phenology timing under the hottest conditions for vineyard plots located at a high elevation compared to plots at a low elevation.

Projected Temperature and Precipitation Changes
The projected higher magnitude of seasonal temperature increase in RCP8.5 than in RCP4.5 (Figure 2) is expected, since RCP8.5 assumed a much higher thermal forcing than that in RCP4.5 [52,53].For seasonal precipitation projections, there is a widespread increase in RCP8.5 compared with more uncertain changes (−6% to 5%) in RCP4.5 (Figure 2).Similar findings were also obtained in several other studies [9,51,80,81].For instance, Dietrich et al. [80] showed the projected mean annual precipitation changes for Germany over 2041-2050 (relative to 1971-2000) are of −2.2% and 6.7% in RCP4.5 and RCP8.5, respectively, for which a higher magnitude of monthly precipitation increase was also found in RCP8.5 than in RCP4.5, spanning most of the months.On the other hand, the projected changes in mean temperature and precipitation sum during the flowering-veraison phase (Figure 4) should be interpreted with care, as they are simultaneously affected by three factors: (1) advanced phenology with shifted phenophase towards the cooler spring season for both varieties; (2) shortened or prolonged phenophase, depending on the variety (Figure 3); (3) different exposures to projected monthly mean temperature and precipitation changes, depending on the variety (Figure 5).However, it is beyond the scope of this article to disentangle each of these factors and quantify the contribution of each of them to the projected changes across Germany.Results show Riesling will remain a late-ripening variety as compared to Müller-Thurgau under future climate conditions, although differences in the flowering stage remain small (Figure 5).The gap between the veraison stage for Müller-Thurgau and Riesling is enlarged in the future period under both scenarios (Figure 5).Therefore, Riesling is more exposed to the high-temperature increase in August, whereas the monthly precipitation change is more uncertain (Figure 5).The extended floweringveraison phenophase (Figure 3b) can be assumed as the major factor (compared to the other two factors) in the projected increases of flowering-veraison precipitation for Riesling (Figure 4b).The decreased precipitation for Müller-Thurgau (Figure 4b) can be mainly attributed to the shortened flowering-veraison phase (Figure 3a).Nevertheless, it should be emphasized that small absolute changes in precipitation amount may appear as large percentage changes, particularly given the relatively small cumulative precipitation amount of the flowering-veraison phase.As Molitor and Junk [19] demonstrated, temperature changes in calendar-based time frames cannot reveal the actual temperature changes in a specific development stage.In this study, the changes for a phase-dependent meteorological condition were assessed; however, this approach is complex, especially when analysing a large geographic area.

Projected Changes of Flowering-Veraison Water Deficits
The impacts of projected climate change on the total flowering-veraison ET max and ET tend to differ between the two varieties (Figures S7 and S8), which represent the water demand and consumption of the investigated soil-crop system, respectively [36,61,65].For Müller-Thurgau, the overall reduction in the sum of ET max in the flowering-verasion phase (Figure S7) can be mainly attributed to the reduced length of the phenophase (Figure 3a), reduced potential evapotranspiration due to shifted phenophase towards a cooler spring and projected flowering-veraison temperature increase (Figure 4a).For Riesling, the increased sum of ET max (Figure S7) is primarily due to the combined effect of an extended floweringveraison phase (Figure 3b), higher flowering-veraison (relative to Müller-Thurgau) temperature increase (Figure 4a) and shifted phenophase towards the spring when potential evapotranspiration is lower.For Riesling, we assume that the length of phenophase plays a major role in determining the changes in ET max , as the identified areas with decreased ET max (Figure S7) are generally consistent with areas with a shortened phenophase (Figure 3b).Besides, it should be noted that the effects of elevated CO 2 levels (in two scenarios) on reducing ET max (mainly via T max ) have already been explicitly considered by STICS for both varieties under the two scenarios investigated [35,65].The ET pattern generally follows that of ET max but is more difficult to interpret, as it additionally interacts with projected precipitation changes and plant growth response (e.g., heat stress limitations) [35].As a result, both varieties show a similar spatial pattern of projected CWSI changes (Figure 6), partly because of the offset of changes between ET max and ET for each variety (Figures S7 and S8).An overall increased CWSI by up to 15% (captured by the ensemble median) is obtained over most of the study area for both varieties (Figure 6).This result is consistent with results reported by Hofmann et al. [13] and Fraga et al. [42].At the vineyard level, Hofmann et al. [13] projected a significantly increased frequency of the occurrence of severe drought stress over a fixed calendar period (1st May-30th September) for Riesling by using a different vineyard model under multiple RCMs.At a large geographic scale, Fraga et al. [42] suggested that vine water stress conditions for Pinot noir from fruit set until maturity are projected to have a particularly high increase in central European wine regions.
Assessing potential climate change impacts can help support winegrowers to implement suitable adaptations for improved viticulture management.The effect of water deficits on red wine quality is well documented; except for situations where it is particularly severe, quality is increased [22,82].The effect of water deficits on white wine quality is less well documented.The threshold of water deficit where white wine quality decreases is probably reached earlier compared to that of red wine, although recent researches suggest that mild levels of water deficit increase the content of some key aromatic compounds in white wines [82,83].Therefore, in areas with projected intensifications of CWSI, shifting to red varieties might be a good option, e.g., replacing the declining variety Müller-Thurgau with Pinot noir, Cabernet-Sauvignon, Syrah or Tempranillo [45], as their quality can benefit from increased water deficits [22,29].This is especially relevant for some current winegrowing regions in the southwest of Germany, where an upper (95th) percentile simulation potentially shows a strong signal of intensified water deficit by around 30% (some areas even >45%) (Figure 5).On the other hand, the identified areas in the northeastern and central regions with alleviated CWSI (by up to 15%) should preferentially be planted with white varieties (Figure 5).Several other adaptation strategies can also be implemented to reduce the potentially detrimental effects of excessive water deficit stress.Grafting vines onto more drought-resistant rootstocks is a cost-effective and environmentally friendly adaptation to increased drought [84].Predominant rootstocks currently used in Germany are V. berlandieri × V. riparia siblings (e.g., SO4, 5BB, 125 AA).Rootstocks from the family V. berlandieri × V. rupestris (in particular 110R) are much more drought-resistant [85].Substantial differences also exist among grapevine varieties regarding their tolerance to drought [86,87].Reduced planting density also relieves water deficit in vineyards because it decreases ET [88].Vines are currently dry-farmed in Germany, and implementing irrigation can also be considered, in particular to limit potential yield losses [14,16].However, irrigation has a high environmental cost because it significantly increases the blue water footprint of wine production [89] and should be used as the last adaptation strategy [7].

Strengths and Limitations
Potential climate change effects on water deficit/scarcity of a drought-sensitive period (flowering-veraison) were assessed for two main German varieties, i.e., Müller-Thurgau and Riesling.The viticulture model employed in this study is one of few models that can capture the interplay between phenology timing and vine growth response to drought stress [35,37,59].Rising temperature alone impacts both phenology and plant water demand, whereas changes in the precipitation distribution during the phenophase additionally affect soil water availability [7].A holistic approach was implemented to address climate change impacts, considering vine development to projected changes not only in temperature and precipitation but also the covariance between other climatic variables (radiation, wind speed and humidity, i.e., these variables can additionally modify vine water demand).Besides, we specifically analysed drought risk during an important phenophase instead of the calendar season, which can produce more relevant information for viticulturalists and winemakers.Furthermore, our findings demonstrate that phenophase-dependent water deficits are affected by multiple and correlated factors, among which changes in the length of the phenophase may play an important role.The results obtained in this study can be applied to any variety with similar precocity for flowering and veraison as Müller-Thurgau (early-ripening) and Riesling (late-ripening).For example, the flowering and veraison dates of Riesling are comparable to those of Merlot [49,90], so the obtained result for Riesling might be applied to Merlot as well.
Nevertheless, when assessing climate change impacts on viticulture, some limitations should be considered [4].Firstly, although the projected magnitude of increase in summer temperature is quite consistent across RCMs, considerable uncertainties are found in projected summer-month precipitation changes (Figure 5).The spread among RCMs is still considerable when simulating the flowering-veraison CWSI (Figure 6, Appendix A).This again illustrates the complexity of climate change impact assessment; projections already diverge substantially for mean changes, and the assessment of interannual variability is even more difficult.Secondly, only one grapevine model was used in this study, and there is a clear need to consider potential variability in simulations through the use of multiple crop models in impact assessment [32,34,91].This point mainly concerns the question of how reliable the S-W module (Section 2.3.2) is for simulation of ET max and ET [65].Since the two varieties are only characterized by their differences in phenology (Table S1), additional lysimeter measurements for these water use variables could be useful to verify the model simulations.Lastly, quantitative information about natural variation and in situ vineyard management is lacking, which does not allow for consideration of the effects of fine-scale

Figure 2 .
Figure 2. Ensemble median of climate model projections for the mean seasonal (April-October) (a) average temperature ( • C) and (b) precipitation sum (mm) in the baseline period (1976-2005) and respective mean changes in the future period (2041-2070) under RCP4.5 and RCP8.5.

Figure 3 .
Figure 3. Ensemble median projections for the mean DOY changes of flowering (BBCH65) and veraison stages (BBCH81) and for the BBCH65-BBCH81 phenophase between the future period (2041-2070) under RCP4.5 and RCP8.5 and the baseline period (1976-2005) for grape varieties of (a) Müller-Thurgau and (b) Riesling.For the phenophase, positive values indicate the extended period, whereas negative values denote the shortened period.DOY: day of year.Negative DOY denotes the advanced days of the phenology stages.

Figure 5 .
Figure 5. Projected monthly average temperature changes ( • C) (in line plot) and precipitation sum changes (%) (in bar plot) over the growing season (April-October) for the future period (2041-2070) under (a) RCP4.5 (yellow) and (b) RCP8.5 (red) as compared to the baseline period (1976-2005) in northern (Latitude 52 • -55 • ), central (Latitude 50 • -52 • ) and southern (Latitude 47 • -50 • ) Germany.In the line plot, the lines represent regional mean values with the respective colour band indicating their 90% uncertainty range among RCMs.In the bar plot, bars indicate the regional mean values with the respective error bars indicating their 90% uncertainty range among RCMs.The regional mean flowering stage and veraison stage are also shown for Müller-Thurgau (triangular) and Riesling (diamond).Note that for the same variety in each subplot, the phenology symbols occurring earlier correspond to the flowering stage, and the next symbol is the veraison stage.The phenology symbols with black, yellow and red are for the baseline, RCP4.5 and RCP8.5, respectively.

Table 2 .
List of used bias-adjusted global climate model-regional climate model (GCM-RCM) chains from EURO-CORDEX.