The Effects of Climate Change on Variability of the Growing Seasons in the Elbe River Lowland , Czech Republic

This research aimed to identify an approach for adaptation of agriculture to increased climate variability and projected changes, taking into account regional specificity of climate change. Changes in the timing of growing season (GS) parameters for both observation and models data were computed using daily mean temperatures for three thresholds that correspond to the physiological requirements of the vegetable types. This research included a new assessment of the potential impacts of climate change on the GS of vegetables grown in the Elbe River lowland, one of the largest farmed vegetable regions in Central Europe. To accomplish this, a comprehensive analysis was conducted of the spatiotemporal variability of the date of the beginning of the growing season (BGS), the date of the end of the growing season (EGS), and the length of the growing season (GSL) for the period 1961–2011. In addition, an assessment was made of the potential changes in the dates of the BGS, EGS, and GSL for the Elbe River lowland, simulated using the regional climate models. Prospective areas for growing thermophilic vegetables in the study region were also determined.


Introduction
The agricultural sector is increasingly exposed to risks, both environmental and economic, due to the phenomenon of climate change and climate variability.In this context, the identification of the best adaptation options is one of the major challenges to improve the risk management tools in agricultural sector at regional and local levels.The Summary for Policymakers of the Working Group I contribution to the IPCC Fifth Assessment Report shows that the risk of extreme events is growing due to global warming [1].Climate change is likely to affect agricultural systems very differently in various parts of Europe [2].In northern areas, climate change may have primarily positive effects through increases in productivity and in the range of varieties grown, while in southern areas the disadvantages of climate change will predominate, with lower harvestable yields, higher yield variability, and a reduction in suitable areas for traditional crops.The results of studies performed in Europe in recent years have shown consistent increases in temperature and various patterns of precipitation, with widespread increases in northern Europe and rather small decreases in southern Europe [3,4].These changes in climate patterns are expected to affect all components of Europe's agricultural ecosystems (e.g., crop suitability, yields, and crop protection).Thus, adaptation strategies should be introduced to reduce the negative effects and exploit the possible positive effects of climate change (e.g., changes in crop species, cultivars, and sowing dates).
Future climate change is projected to increase the length of the growing season [3,5].An increase in the length of the growing season, together with warmer climate during the growing season, may increase the potential for growing thermophilic vegetables in open fields in lowland areas in Central Europe and increase the potential number of harvests.For spring crops, climate warming will allow earlier planting or sowing than at present.Earlier planting in spring increases the length of the growing season; thus, earlier planting 2 Advances in Meteorology using long-season cultivars will increase yield potential if soil moisture is adequate and the risk of heat stress is low [6].All the reported studies concern Europe, where recent warming has clearly advanced a significant part of the agricultural calendar.Increased temperatures, associated with earlier late spring frost dates and delayed autumn frost dates, are clearly apparent in temperate regions of Europe [3,[7][8][9].Advances in production due to earlier springs have been documented for several types of crops [5,10], and the benefits to viticulture of longer growing seasons have been established [11,12].For example, according to [5], the growth stages of wild plants and fruit trees as perennial crops have advanced significantly since the beginning of the 1980s and were, together with true agricultural phases, occurring 1 to 2 weeks earlier by the 2000s.European viticulture has experienced direct effects of climate changes, mainly in terms of shifts in grapevine phenological events and stages.Reference [12], exploring viticulture and climate change relationships in Greece, noted a systematic shift in grape harvest dates (earlier harvests).Similar trends have been observed in Italy.Grape maturity dates have gradually shifted 19 d earlier over the period 1961-2009 for numerous varieties [11].
According to [13], the combination of increased air temperature and changes in the amount and annual cycle of precipitation will lead to further shifts in the area and location of individual agroclimatic zones in Central Europe.It is worth noting that, due to climate change, the breeding of new and improved vegetable crop varieties can lead to an extension of the border areas suitable for the profitable cultivation of vegetables.Some thermophilic vegetables that currently grow mostly in southern Europe (e.g., melons, eggplants, tomatoes, and peppers) can become more suitable for cultivation in lowland areas in Central Europe.The two important field vegetable regions in the Czech Republic are South Moravia and the Elbe River lowlands, in which partially different assortments of vegetables are grown.The Elbe River lowland has traditionally been a region of cultivation of brassica vegetables (kohlrabi, Savoy cabbage, white-headed cabbage, and cauliflower), while South Moravia is a profitable region for thermophilic vegetables (e.g., tomatoes and peppers).These differences in the assortments of cultivated vegetables grown in the two regions are mainly due to differences in the temperature conditions of these regions.In the warmest parts of the Elbe River lowland, growing thermophilic vegetables such as tomatoes (Solanum lycopersicum L.) and cucumbers (Cucumis sativus L.) are profitable only in warmer years.Based on our previous studies [9,14], we note that the combination of changes in European agricultural commodities and ongoing climate changes (increases in temperatures) can lead to higher costs for vegetable imports and stricter requirements for the maturity and quality of yields.In addition, favourable national agricultural policies could extend the areas suitable for thermophilic field vegetables from the hottest regions of South Moravia to the Elbe region.In addition to the current assortment of vegetables grown, nontraditional vegetables such as melons (Cucumis melo L.) and eggplants (Solanum melongena L.) could also be grown.
Vegetable crops are most sensitive to the timing of cold events at the beginning and end of the growing season.
Damage caused by late frosts in the spring or early frosts in the fall is a limiting factor, particularly for vegetables grown in Central Europe.In previous study [9], we examined the variations in the last spring frost, the first fall frost, and the length of the frost-free period for the Elbe River lowland in the Czech Republic as indicators of climate variations in this region.Our results demonstrated shifts toward earlier last spring frosts, later first fall frosts, and longer frost-free periods in the Elbe River lowland.Longer frost-free periods can be particularly beneficial for thermophilic vegetables in lowlands.In terms of the growth of field vegetables, however, a late spring frost remains a risk factor, although the degree of risk has decreased [9].
Despite the observed regional trends in agroclimatic characteristics [9,14] and the importance of the horticultural sector to the economy of the Czech Republic, there is an utter lack of research on the impacts of climate variability on the assortment of field vegetables grown in the Elbe River lowland.This subject continues to be of interest to agroclimatologists, although the range and yield of vegetable crops grown in Central Europe are primarily controlled by temperature.Therefore, the primary objective of this study was to evaluate the long-term changes in the timing of the growing season (GS) parameters over the Elbe River lowland for 1961-2011 at a high horizontal resolution using daily mean temperatures for three thresholds ( mean ≥ 5, 10, and 15 ∘ C) that correspond to the physiological requirements of the vegetable types.To provide more information on the potential for extension of the range of vegetables grown, the specific objectives of the present study were (i) to conduct a comprehensive analysis of the spatiotemporal variability of the dates of the beginning (BGS) and end (EGS) of the growing season and the length of the growing season (GSL) for three threshold temperatures for the period 1961-2011; (ii) to assess potential changes in the BGS, EGS, and GSL for three threshold temperatures in the Elbe River lowland as simulated using the regional climate models (ALADIN-Climate/CZ and RegCM) under the A1B SRES scenario (for the periods 2021-2050 and 2071-2100); and (iii) to determine prospective areas for growth of thermophilic vegetables in the study region based on projected climatic data provided by regional climate models.

Gridded Datasets and Quality Control.
The Elbe River lowland is one of the largest and most productive suppliers of vegetables in the Czech Republic (Figure 1).The study was based on gridded 116 grid points' daily series of mean air temperature ( mean ) data at a 10 km horizontal resolution for observed  and future (2021-2050 and 2071-2100) climate conditions.A regular gridded network (CZGRIDS, ALADIN-Climate/CZ) established by the Czech Hydrometeorological Institute (CHMI) was applied (Figure 1).High-density gridded datasets allow very precise and detailed delimitation of areas with advanced and/or delayed growing seasons compared with station network datasets.The gridded network was created using the technical series from 268 climatological stations and 774 rain gauge stations in the Advances in Meteorology 3 CHMI network (series after quality control, homogenization, and filling of all missing values in the period 1961-2011).New values for the regular 10 × 10 km grid network were interpolated from these technical series [15].The gridding and all data processing, including the analyses described in this paper, were performed using the ProClimDB database software for climatological dataset processing.
A remarkable aspect of vegetable production is that plants function within quite narrow temperature limits.The extreme temperature range is between the killing frost temperature of 0 ∘ C and death by heat and desiccation at 40 ∘ C [16].The wide range of vegetables grown in the studied region can be divided into three basic types according to their sensitivity to low temperatures at the onset of the growing season: thermophilic vegetables (e.g., tomato, pepper, pumpkins, and cucumber), cold-resistant vegetables (e.g., early kohlrabi, summer savoy cabbage, late cauliflower, late cabbage, late carrots, and celeriac), and frost-resistant vegetables (e.g., onion and root parsley).From an agronomic point of view, thermophilic vegetables should be planted after 15 May to minimise the risk of frost damage [9].The critical temperature for thermophilic vegetables after May 15 is min ≤ −0.1 ∘ C [17][18][19].A severe last spring frost after April 15 has occurred every second year during the sowing/planting of cold-resistant vegetables at higher altitudes.Typically, the date of planting/sowing of frost-resistant vegetables in the Czech Republic is at the end of March, and a gradual shift occurs from the hottest regions of South Moravia towards the Elbe region.

Regional Climate Models and Data
Postprocessing.Additionally, the ALADIN-Climate/CZ and RegCM regional climate models (RCMs) were adopted to calculate possible shifts in the start, end, and length of the climatological growing season under the SRES A1B scenario for two future periods (2021-2050 and 2071-2100) over the Elbe lowland.SRES scenario A1B is a baseline scenario referred to in [19].We used two different RCMs, which were driven by two different global circulation models (GCM), namely, ARP ÉGE-Climate [20] and ECHAM5 [21].
The RegCM model was originally developed [22] and then augmented and used in various reference and scenario simulations [23,24].RegCM family using RegCM transient ENSEMBLES run for whole Europe currently in 25 km resolution driven by transient run of ECHAM5, but within CECILIA-FP6 Project (Central and Eastern Europe Climate Change Impact and Vulnerability Assessment) was developed RegCM in 10 km resolution, because is important for impact studies in regional scale [25].
ALADIN-Climate/CZ performs quite well in comparasion to the other RCMs in impact studies [26,27].The model is able to capture the main features of the present climate of the Czech Republic and works well over smaller areas with rather complex orography (tablelands, valley, and hills), such as the Elbe River lowland.However, it should be always kept in mind that model simulations of future climate conditions involve many uncertainties, and it is necessary to apply some type of postprocessing of model outputs from the RCM simulation before using the results in other applications.Therefore, the observed daily temperature data were first transferred onto a regular grid of the RegCM and ALADIN-Climate/CZ models.Prior to the calculation of the start, end, and length of the climatological growing season for the RegCM simulations forced with the ECHAM5 GCM and ALADIN-Climate/CZ run forced by ARP ÉGE-Climate GCM, the daily temperature for each grid point has been bias corrected against the systematic errors induced by the GCM.The bias correction was applied to the scenario runs.The bias correction method is based on variable correction using individual percentiles whose relationships are derived from observations and a control RCMs simulation [28].After the correction, the model outputs (for 2021-2050 and 2071-2100) were fully compatible with the gridded observation dataset .Figure 2 shows the annual variation of bias corrected daily temperature means calculated for 30 years, corresponding to the RegCM (forced by ECHAM5) and ALADIN-Climate/CZ (forced by ARP ÉGE-Climate) in the A1B scenario runs (2021-2050 and 2071-2100), respectively.The results show that the projected temperature means for all months in the A1B scenario runs will increase compared to the current climate .
For the Czech Republic as a whole, ALADIN-Climate/CZ and RegCM simulated increases in temperature, especially in the summer and winter.Over the Elbe lowland region, the temperatures are projected to increase by the end of the 21st century compared to the mid-21st century and the reference period 1971-2000.The 30-year annual temperature means for the Elbe lowland are 8.4 ∘ C for the current climate (1971-2000) and 9.6 ∘ C for 2021-2050 and 11.4 ∘ C for 2071-2100.The mean summer temperature increases over all models are projected to be 1.7 ∘ C and 4.0 ∘ C in 2021-2050 and 2071-2100, respectively, relative to 1971-2000.The projected increases in the spring and autumn temperature means under the A1B scenario are 0.5 ∘ C and 2.3 ∘ C, respectively, in 2021-2050 and 1.1 ∘ C and 2.8 ∘ C, respectively, in 2071-2100.Overall, these projections suggest increase in the GSL in 2071-2100 relative to the present climate.Both RCMs simulations for the A1B SRES scenario reveal that the summer temperature, which has an important influence on the beginning and end of the GS ≥ 15 ∘ C in hilly areas, will increase by the end of the 21st century.Increasing temperatures can lead to earlier harvests of thermophilic vegetables.However, responses to temperature changes may differ among vegetable varieties, mainly in terms of their temperature thresholds for optimum ripening [16].

Determination of BGS, EGS, and GSL.
The effect of climate change on the GS of agricultural crops can be determined through phenological observations, the normalised difference vegetation index from satellite data, and surface air temperatures [7,[29][30][31][32][33][34][35][36].Definitions of the growing season using surface air temperatures may be considered valid in areas where the GS is largely limited by temperature.In our study, we used threshold temperatures for a predefined number of days to start and end the climatological GS, which has good spatial coverage over the regular grid network.However, there is no universal definition of the climatological GS depending on the geographic position and associated climate regimes, and there may be large disagreement in GSL changes depending on the definition used [34].The climatological GS can be defined as the entire period in which growth can theoretically take place [31].In an attempt to provide a common benchmark for investigators, [32] selected the 5 ∘ C threshold to standardise the definition of the GS as the period when daily temperatures remain above 5 ∘ C for more than 5 days.Reference [31] defined the start of the GS as >5 ∘ C for ≥5 days and the end as a 10-day period of mean daily temperatures below 5 ∘ C. In the Greater Baltic Area, [34] defined the start of the growing season as the last day of the first six-day period with daily mean temperatures above 5 ∘ C after the last winter/spring frost.In our study, the GS parameters were defined with respect to the physiological requirements of the vegetable types and the classification system of the Climate Atlas of Czechia [37] (the large growing season ( mean ≥ 5 ∘ C), the main growing season ( mean ≥ 10 ∘ C), and the summer growing season ( mean ≥ 15 ∘ C)).
The sowing/planting period differs among field vegetable types; therefore, we defined the BGS as the first five consecutive days during which the daily mean temperature equals or exceeds 5, 10, or 15 ∘ C and, subsequently, does not drop below  mean 5, 10, or 15 ∘ C for five days or more.The EGS is defined as the first day of the first 5-day period with a  mean below 5, 10, or 15 ∘ C. The GSL is the number of days between the BGS and EGS for the three threshold temperatures.Using the three thresholds of daily  mean is justified because they correspond to the vegetation period of the majority of field vegetables.BGS ≥ 5 ∘ C represents the average date of the sowing period of the majority of field frost-resistant vegetables, whereas the start of BGS ≥ 10 and 15 ∘ C corresponds to the sowing/planting periods of field cold-resistant and thermophilic vegetables, respectively, in the studied region.The EGS corresponds to the harvest period of different types of vegetable crops.The GSL is the time from sowing to harvest or the time from planting to harvest.
Using these definitions, we produced time series of the start, end, and length of the growing season for 116 grid points for both the observations and the model datasets.Furthermore, we created and statistically processed two datasets.First, to characterise the spatial patterns of the GS parameters, the annual BGS, EGS, and GSL series for the three thresholds were averaged at each grid point for the current climate (1961-2011) and two future climates (2021-2050 and 2071-2100).Second, to examine temporal variations in the GS parameters over the Elbe River lowland as a whole, the BGS, EGS, and GSL for the three thresholds were arithmetically averaged over all 116 grid points to obtain a time series for the current climate (1961-2011) and the two future climates (2021-2050 and 2071-2100).Then, projected changes in the dates of the start and end and the length of the growing season are analyzed by comparing the model simulations for the periods 2021-2050 and 2071-2100 with the reference period 1971-2000.
To investigate changes in the timing of the observed GS parameters, regional anomalies in the start, end, and length of the growing season for the three thresholds over the Elbe River lowland from 1961 to 2011 were calculated.Positive anomalies indicate dates later in the year (fall and spring occurrences) or a greater number of days (length).We then classified anomalies in the growing season onset, end, and length into two categories: (1) years with advances and/or increases in the growing seasons and (2) years with delays and/or decreases in the growing seasons.We also calculated linear trends and analysed the regression slopes of the GS parameters and assessed their significance using Student's one-tailed -test.Finally, to identify prospective areas for growing thermophilic vegetables, kriging interpolation tools in ArcGIS software were used to create maps of the spatial distributions of the GS parameters for the current and future climate conditions.Interpolation, where needed, is based on local linear regression (dependence of given meteorological element on altitude) and universal kriging interpolation method.To generate the maps at regional scale we used CLIDATA GIS environment and CLIDATA DEM (digital elevation model) method with a horizontal resolution of 500 m and 40 km radius of regression.  1  and 2).The onset of spring was very early in 1990 due to temperature anomalies from February to the second half of March that resulted in  mean deviations from +5.0 to +11.0 ∘ C. The GSL in 1990 also ranked among the four longest growing seasons between 1961 and 2011 (Table 2).In this year, the anomalies in the BGS in Central Europe (up to −27 days) were due to temperature anomalies of up to +4.5 ∘ C that occurred from February to April.These findings are similar to those reported by [7].The majority of these years appear to be associated with temperature anomalies in Central Europe [7,8,29].Our results are in line with these findings and offer additional insights into the complexity of these phenomena on a regional scale.Although the growing seasons have ended later, the BGS has made a greater contribution to the increase in GSL (1961, 1981, 1998, and 2009) (Table 2).A later start, early end, and shorter GSL occurred in 1970 and 1972.A shortening of the growing season occurred in 1971 and is attributed to an early end of EGS ≥ 10 ∘ C, whereas the shortening of the GS that occurred in 1965 and 1987 is attributed to a later onset of the growing season in those years.

Results and Discussion
The variability in strong anomalies in the start, end, and length of the growing season for  mean ≥ 15 ∘ C is similar to that for GS ≥ 5 and 10 ∘ C (Table 2).Positive anomalies for GSL ≥ 15 , which led to advanced dates of planting of field thermophilic vegetables in the studied region.These growing seasons appear to be associated with high positive temperature anomalies and high deficits in the water balance (i.e., exceptionally warm and dry spells) throughout the country.The year with the latest start of BGS ≥ 15 ∘ C, the earliest end of EGS ≥ 15 ∘ C, and the shortest GSL was 1965.This year has also been identified in previous study of the Czech Republic [14] as having been extremely wet and cold.Late-onset and shortened growing seasons also occurred in 1972, 1974, and 1980, while an early end and shorter GSL ≥ 15 ∘ C occurred in 1970, 1971, and 1978.
We observed that the majority of years with delays and/or decreases in the growing seasons occurred in the cool and wet decade of 1971 to 1980, which was characterised with persistently higher than normal precipitation and the lowest negative deviations in daily mean temperatures since 1961.During this decade, the regional average of the ending of mild spring frosts occurred in the second half of May, and the years had shorter frost-free periods [9].
Figure 3 displays the interannual variations of the areaaveraged BGS, EGS, and GSL for average daily temperatures of at least 5 ∘ C, 10 ∘ C, and 15 ∘ C in the Elbe River lowland during the period 1961-2011.As the left panel of Figure 3 shows, the BGS exhibited high interannual variation.Positive anomalies predominated during 1963-1973 for BGS ≥ 5 and 10 ∘ C and during 1970-1980 for BGS ≥ 15 ∘ C, whereas negative anomalies predominated during 1989-1997 and 2000-2011.Although the annual BGS ≥ 15 ∘ C significantly advanced during the 1990s, the signal was reduced in the 2000s.Overall, the long-term variations in BGS exhibited earlier spring features, which are consistent with the findings of other climatological research [3,34,36].Moreover, the majority of phenological studies show similar trends: earlier springs and longer periods of plant growth due to greater changes in the timing of the BGS rather than in the timing of the EGS [3,34,36].The annual changes in EGS also exhibit strong interannual variations (middle panel of Figure 3).In contrast to BGS ≥ 5 and 10 ∘ C, negative anomalies were observed during 1970-1980, whereas positive anomalies were observed during the 1990s and 2000s, which indicates a delayed EGS in recent decades.Asymmetric long-term variations in BGS and EGS led to increases in the total GSL (right panel of Figure 3).With the weakening of the earlier BGS trends, greater delays in EGS clearly led to longer GSLs from 1995 through the 2000s.Conversely, in the first decade of the 21st century, the GSL ≥ 15 ∘ C slightly shortened as a result of the delayed start of the growing season; however, the end of EGS ≥ 15 ∘ C occurred later, which may have had a positive effect on the ripening and quality of thermophilic vegetables.A later onset of GSL ≥ 15 ∘ C can be mitigated by agrotechnical measures.Using remote measurements of the vegetation index and daily mean temperatures, [35] found that the lengthening of the GS in Europe in the last decade could be attributed mainly to the delay in EGS rather than to the advance of BGS.Moreover, the temporal evolution of GS in different climatic regimes is characterised by large spatial differences.Climatic indicators (e.g., the frost-free period and the thermal growing season), remote sensing-based indicators, and phenological studies all indicate shifts in the timing and length of the growing season [3,5,7,8,10,31,32,34,36].

Trends of the Growing Season.
To obtain additional insights into the changes in the timing of the GS parameters in the area-averaged time series, the linear trends and regression slopes of the growing season onset, end, and length were calculated.The significance of these trends was tested using Student's one-tailed -test at the 95% significance level.The linear trends of the beginning, end, and length of the growing season are tabulated in Table 3. Increasingly negative BGS trends and positive EGS trends led to increasing GSL trends.Over a 51-year period, the GS ≥ 5 ∘ C started 10.7 days earlier and ended 4.6 days later, yielding an increase in the GSL of 15.3 days.The BGS ≥ 10 ∘ C advanced by 7.6 days on average, and the EGS ≥ 10 ∘ C was delayed by 3.6 days on average, resulting in an average lengthening of the GSL ≥ 10 ∘ C by 11.2 days.The GSL ≥ 15 ∘ C advanced by only 5.1 days, with smaller changes during the spring (1.0 day) than in the fall (4.1 days).
To understand the evolution of significant changes in the growing season series, the analysis period was divided into five decades (Table 3).The results suggest that the Elbe lowland has experienced the following changes in its growing season parameters over the last five decades.(1) Although advances in the BGS have been reported as the main factor in longer growing seasons, in this study, delays in the EGS were found to be a more important factor regulating GSL changes during the 1990s and 2000s at thresholds of  mean ≥ 15 ∘ C. (2) We found that the GSL ≥ 5 and 10 ∘ C increased considerably during 1991-2000 (>0.77 days yr −1 ), but the lengthening was less pronounced during 1961-1970 and 1971-1980. (3) The largest changes in the start of BGS ≥ 5 and 10 ∘ C (earlier growing season starts of approximately one week) in the 2000s occurred because of warmer temperatures and earlier ends of the last frosts.(4) General trends towards later growing season ends were observed, except for 1961-1970, for which shifts of −0.7 and −6.3 days decade −1 towards an earlier growing season end were observed.(5) In the Elbe River lowland, BGS and EGS exhibited shifts towards earlier  4).Overall, this is a spatial difference of 15 days.Thus, the BGS ≥ 5 ∘ C started 4 days later per 100 m of altitude.The map of BGS ≥ 5 ∘ C defines two main areas with a start of the growing season before March 21.These areas can be considered climatically close to the optimum for field vegetable production.However, the occurrence of severe frosts ( min ≤ −2.2 ∘ C) during the planting of field frost-resistant vegetables in the growing area of the Elbe lowland is, on average, 56.9% [9].It follows that, despite the considerable resistance of these vegetables to low temperatures, it is necessary to choose planting areas with southern exposures (or cover plants with nonwoven textiles).
Utilising resistant varieties and hardening seedlings before planting are advisable.Over the studied region, the GS ≥ 10 ∘ C starts nearly 4 weeks after BGS ≥ 5 ∘ C. Geostatistical analysis indicates that the BGS ≥ 10 ∘ C starts roughly 2 days later for every 100 m of elevation.Consequently, the BGS ≥ 10 ∘ C starts between April 15 and 30, and a start before April 15 can be observed in the traditional vegetable-growing regions (up to 300 m).In the warmest areas of the middle Elbe lowland, the middle Poohří, and the northeast of the Prague platform, the mean date of the BGS ≥ 10 ∘ C occurs between April 16 and 20.Conversely, in the hilly lands of the northern and eastern parts of the Elbe River lowland, the start of BGS ≥ 10 ∘ C occurs in the end of April.In addition, most vegetables grown in the warmest areas fall within the low-risk category for severe spring frosts.It will therefore be possible to extend the areas suitable for growing cold-resistant vegetables (i.e., mostly Brassicas) towards the northeast of the middle Elbe lowland.
The BGS ≥ 15 ∘ C starts between May 8 and 18, approximately 1.8 days later per 100 m of altitude.The latest beginning dates of BGS ≥ 15 ∘ C after May 15 are related to high altitudes in the region (higher than 300 m) and frost hollows, whereas the earliest starts before May 8 are observed in the lower areas (lower than 250 m).The risk of frost after May 15 in the traditional vegetable-growing regions (up to 250 m) is low.At higher altitudes, areas with zero incidences of negative minimum temperatures were found, which may allow for the possible expansion of the area of thermophilic vegetable cultivation [9].Over the Elbe River lowland, compared with the start of the growing season, the end of the growing season pattern follows the altitudinal gradient to a lesser degree.There are no large differences in the spatial distribution of the growing season end between  mean ≥ 5, 10, and 15 ∘ C (middle panel of Figure 4).The EGS ≥ 5 ∘ C dates range from 298 days (October 25) to 314 days (November 10), and their spatial distributions are inversely related to the climatology of the BGS dates.In 15% of the study area (hilly lands and frost hollows), the earliest end of the EGS ≥ 5, 10, and 15 ∘ C occurs before October 31, October 5, and September 5, respectively.The latest ending of EGS ≥ 5, 10, and 15 ∘ C (after November 6, October 11, and September 11) occurs in the northeastern Prague plateau and the middle Elbe River valley, which together cover 20% of the study area.The later end of the growing season in those regions is consistent with an observed earlier start of the growing season in the same region.
The pattern in GSL is consistent with the pattern in the onset of the growing season; GSL pattern also shows some local differences due to the heterogeneous end of the growing season pattern, especially for  mean ≥ 15 ∘ C (right panel of Figure 4).A total of 50% of the area of the Elbe River lowland has a GSL ≥ 5 ∘ C longer than 225, which satisfies the requirements of root vegetables (celeriac), which have the longest growing seasons.The GSL ≥ 10 ∘ C is between 150 and 180 days.65% of the area has a growing season shorter than 170 days.

Projected Changes in the Dates of the Start and End and the Length of the Growing Season
An assessment was made of the potential changes in the dates of the BGS, EGS, and GSL for the three threshold It is worth noting that, in addition to the types of vegetables traditionally grown in the Elbe River valley, extension of the cultivation of thermophilic varieties, especially Cucumis melo L., has already been observed [9].This extension is closely linked to warming of the climate system.Consequently, the cultivation of thermophilic vegetables will complement the cultivation of traditional assortment of vegetables.Thus, a detailed evaluation of the past long-term changes in the GS parameters is essential for predicting the effect of future climate variability on the range of vegetables that can be cultivated in this region.In general for both RCMs, the end of the growing season is projected to be significantly delayed (i.e., occur later) under the climate change scenario compared with the current climate.The medians of the EGS ≥ 5, 10, and 15 ∘ C by ALADIN-Climate/CZ are projected to occur on days 315 (11 Nov.), 284 (11 Oct.), and 255 (12 Sep.), respectively, during 2021-2050.The delay will be greatest during the period 2071-2100, when the mean dates of EGS ≥ 5, 10, and 15 ∘ C are projected to occur on days 318 (14 Nov.), 293 (20 Oct.), and 267 (24 Sep.), respectively.The latest dates of the end of EGS ≥ 5, 10, and 15 ∘ C for the period 2021-2050 are projected to occur on days 345 (11 Dec.), 302 (29 Oct.), and 273 (30 Sep.), respectively.The later end of the growing season will extend the vegetation period of field vegetables.The ALADIN-Climate/CZ and RegCM simulations project an increase in the average growing season length; however, high interannual variability is projected, as in the current climate (Table 4).Under the A1B scenario for the period 2021-2050, the shortest and longest GSL ≥ 5 ∘ C (GSL ≥ 10 and 15 ∘ C) vary between 182 (154 and 86 days) and 272 (207 and 140 days) days.Our regional projections for the end of the 21st century indicate that GSL ≥ 5 ∘ C (GSL ≥ 10 and 15 ∘ C) could be significantly extended by an average of 267 days (212 and 134 days).
Table 5 displays the predicted changes in the mean start, end, and length of the GS between the reference period  and the two future climate periods (2021-2050 and 2071-2100) for the Elbe lowland.As expected, the magnitudes of the changes for the two future periods are different.The projections based on two different RCMs by A1B scenario suggest that, during the mid-21st century, the dates of the BGS ≥ 5, 10, and 15 ∘ C will be earlier (by ≥5 days) and those of the EGS ≥ 5, 10, and 15 ∘ C will be later (by ≥5 days) across the studied region.Consequently, the GSL ≥ 5 ∘ C (GSL ≥ 10 and 15 ∘ C) will be increased under the A1B scenario by ≥11 days (≥13 and ≥11 days) compared with the current climate.
The advance of BGS ≥ 5, 10, and 15 ∘ C will be greater at the end of the 21st century (by ≥15 days).The ALADIN-Climate/CZ simulation for the period 2070-2100 projects delays in the EGS ≥ 5, 10, and 15 ∘ C by 12, 13, and 14 days, respectively, whilst RegCM simulation for the period 2070-2100 also projects delays in the EGS ≥ 5, 10, and 15 ∘ C by 12, 18, and 19 days, respectively.The ALADIN-Climate/CZ simulation for the BGS ≥ 5 ∘ C start will be advanced by 22 days, and the EGS ≥ 5 ∘ C end will be delayed by 12 days for the period 2071-2100.Therefore, the A1B scenario results in a projected lengthening of the GSL by 44 days on average compared with the current climate.

Spatial Distribution of BGS, EGS, and GSL Projections.
For the Czech Republic as a whole, the daily mean, maximum, and minimum temperatures simulated by ALADIN-Climate/CZ and RegCM for the A1B scenario are projected to increase more in Bohemia than in Moravia and at higher altitudes [27].These temperature changes affect the growing season parameters, and thus, ALADIN-Climate/CZ and RegCM project greater changes in those parameters for the hilly lands of the Elbe River lowland.According to the ALADIN-Climate/CZ simulation for the period 2021-2050, at higher elevations in the Elbe lowland the GSL will correspond to that of the warmest areas in the current climate (Figures 6(a) and 6(b)).The lengthening of the growing season is most apparent for GSL ≥ 5 and 10 ∘ C. At the end of the 21st century, the GSL in the hottest areas will be significantly lengthened.Additionally, the areas with suitable conditions for growing vegetables (i.e., significant advance and delay of GS) will shift from the lowlands of the Elbe River valley to higher elevations (Figures 6(a) and 6(b)).However, profitable cultivation of vegetables in these areas is complicated by the complex sloping terrain.Therefore, the shift of areas with longer GSL from lowlands to higher elevations could lead to higher potential agroclimatic productivity, which will, however, be difficult to utilise due to the lack of water resources for irrigation and the inaccessibility of these areas to machinery.Significant changes in the BGS can be expected mostly in hilly areas (Figures 6(a) and 6(b)).At the end of the 21st century, two main areas with different early growing season onsets will stand out: the area with the earliest onset of BGS, which will be located in the northwestern Elbe River lowland (i.e., in areas of frost hollows and hills), and a newly created region, which corresponds to currently peripheral areas of vegetable farmland (eastern part of Elbe River valley).The ALADIN-Climate/CZ and RegCM simulations reveal that, by the end of this century, Czech farmers could be growing their crops for up to two months longer.The Elbe River lowland can become a major producer of vegetables in the Czech Republic and significantly increase its competitiveness in the production of market vegetables.Our results show that, under the A1B scenario for the two future climate periods considered (2021-2050 and 2071-2100), the projected increases in the growing season are due to significant projected advances of the starts of BGS ≥ 5, 10, and 15 ∘ C. The BGS simulation data also show clear tendencies towards an earlier start, especially for GSL ≥ 5 and 10 ∘ C. The end of the growing season is projected to occur later on average, but these projected shifts are smaller than those for the beginning of the growing season.Based on these results, it is evident that regional climate changespecifically, increases in temperatures during the growing season-will affect Czech vegetable crops.The changes in the Elbe lowland climate are expected to significantly advance the sowing/planting dates of vegetables under current and future climate conditions.Similar findings have been reported for agricultural areas worldwide [4,6,11,12,30,33,35,36,38].

Summary and Concluding Remarks
In this study, a comprehensive analysis of the current climatic conditions (1961-2011) and possible changes in the climate in the near future (2021-2050) and at the end of the century (2071-2100) was undertaken using both observed gridded data and two RCMs driven by two different GCMs of the changes in the timing of GS parameters for the Elbe River lowland for the first time.The ALADIN-Climate/CZ and RegCM models, with a high horizontal resolution of 10 km, provided the daily step temperature series to project changes in the GS parameters.The date of sowing/planting and the harvest period depend on the vegetable variety, agroclimatic conditions, commercial targets, market constraints, and horticultural practices.This is the probable reason for the absence of data on the sowing/planting and harvest dates of vegetables in the Czech Republic and in other countries.Systematically recorded information on phenology is utterly lacking for most types of vegetables at both the local and the global scales [16].
This research aimed to identify an approach to select suitable sites for extension of new thermophilic assortment of vegetables in Elbe lowland, taking into account regional specificity of climate change, and also to determine prospective areas for growing thermophilic vegetables in the study region using regional climate models.Proposed study can be crucial in development of strategies on climate change adaptation for different varieties of thermophile crops for future climate change in different regions in order to increase productivity while reducing both the cost of farmers and the water footprint of agriculture per unit product.Moreover, our study is in line with the two aims defined in the EU strategy for adaptation of agricultural production system through scientific research and other actions to advance sustainable agriculture development that benefits local, regional, and European farmers and to develop climate-informed crop.The results can be also relevant to the European policy, especially to the Common Agricultural Policy by contributing to good agricultural practices and sustaining the rural communities.The high degree of uncertainty in expected temperature and precipitation changes implies investment risks for deciding on long-term changes in farming technology such as adopting irrigation or switching to new crops, requiring new machinery, infrastructure, and agronomic and market knowledge [39].The salient results are summarised below.
(i) The temporal evolution of anomalies in the GS onset, end, and length displays two distinct periods: a delay in the onset in spring and a shortening of the GSL in the 1960s and an intensified lengthening of the GSL since the 1990s.The majority of years with strong delays in the onset and/or decrease in the length of the growing season occurred in the decade of 1971-1980, which was cooler and wetter than the other decades of the reference period analysed.In the last two decades of the period (1991-2011), BGS and EGS exhibited significant shifts to earlier and later dates, respectively.Consequently, the end of the 20th century and the beginning of the 21st century will be a suitable period for the extension of the cultivation of varieties of vegetables with longer growing seasons and higher temperature requirements.
(ii) From 1961 to 2011, the Elbe River lowland experienced increasingly negative BGS and positive EGS trends that led to an increasing GSL trend.However, the roles of the changes in BGS and EGS differed depending on the time period and on the chosen temperature thresholds.The regional GS ≥ 5 ∘ C (≥10 and 15 According to the RCMs simulation, the most significant shifts in the dates of the beginning and end of the GS are projected to occur in hilly areas; in lowland areas, these changes are not projected to be as pronounced.
The lengthening of the GSL and the flat topography plateaus will create favourable conditions for the expansion of vegetables areas, mainly towards the eastern part of the Elbe River basin.The results also suggest potential for a northerly expansion of vegetables cultivation, although most of the lands to the north of the current frontier will remain only marginally suitable for growing field vegetables due to their complex terrain.
This study can be considered an initial step towards assessing the potential impacts of climate change on the types of vegetable crops grown in the Elbe River lowland.In our future work, we plan to use mechanistic crop simulation models (e.g., CROPGRO-Tomato model) and multimodel multiscenario ensemble (for the new emission scenarios RCP4.5 and RCP8.5) to project potential consequences of climate change on fresh-market tomato (Solanum lycopersicum L.) grown under open field conditions at farm scale in the Elbe lowland.

Figure 1 :
Figure 1: The boundaries of Elbe River lowland in the Czech Republic and location of the 116 grid points and their elevation (m a.s.l.) in the Elbe River lowland (a) and the relief map of Elbe River lowland (b).

Figure 3 :
Figure 3: Interannual variations of area-averaged BGS (beginning of the growing season), EGS (end of the growing season), and LGS (length of the growing season) for the 5, 10, and 15 ∘ C temperature thresholds from 1961 to 2011 in the Elbe River lowland.The black line indicates the climatology of each event (average), and the grey line indicates the 5-year moving average of each parameter.

Figure 4 :
Figure 4: Spatial variability of the mean values of BGS, EGS, and GSL for the 5, 10, and 15 ∘ C temperature thresholds from 1961 to 2011 in the Elbe River lowland.

4. 1 .
Temporal Variability of the Projection of BGS, EGS, and GSL.The projected dates of the start, end, and length of the growing season across the Elbe lowland as a whole for the three threshold temperatures and the two future periods considered (2021-2050 and 2071-2100) are presented in Table 4.For both RCMs the dominant dates of BGS ≥ 5, 10 and 15 ∘ C are projected to advance significantly (i.e., occur earlier) than under the current climate conditions.According to ALADIN-Climate/CZ, in the mid-21st century period, the medians BGS ≥ 5, 10, and 15 ∘ C are projected to advance to days 81 (22 Mar), 107 (17 Apr), and 141 (21 May), respectively.At the end of the 21st century, the medians are projected to advance significantly, to days 59 (28 Feb), 81 (22 Mar), and 133 (13 May), respectively.A similar pattern was shown by RegCM (Table4).These shifts reflect the higher projected temperatures in the A1B scenario.The earliest and latest starts of BGS ≥ 5, 10, and 15 ∘ C for the period 2021-2050 are projected to occur on days 51(20 Feb.), 80 (21 Mar.), and 120 (30 Apr.), respectively, and on days 107 (17 Apr.), 123 (3 May), and 166 (15 Jun.), respectively.The earliest onset of BGS ≥ 5 is projected to shift significantly towards the winter season by the end of the 21st century (2071-2100).

Figure 5 :
Figure 5: Hovmoller-type diagram of the spatiotemporal evolution of large growing season ( mean ≥ 5 ∘ C), the main growing season ( mean ≥ 10 ∘ C), and the summer growing season ( mean ≥ 15 ∘ C) projected from 116 grid points × 30 years × 365 days over the Elbe River lowland based on ALADIN-Climate/CZ and RegCM simulation data under the A1B SRES scenario runs (2021-2050 and 2071-2100).

Table 1 :
Statistical characteristics of the area-averaged growing season parameters of the studied area for 1961-2011. mean ≥ 5 ∘ C  mean ≥ 10 ∘ C  mean ≥ 15 ∘ C

Table 2 :
[7]a-averaged anomalies (sorted by the highest number of anomalies in days) in the start, end, and length of the growing season for the three thresholds of  mean ≥ 5, 10, and 15 ∘ C over the Elbe River lowland for 1961-2011.The beginning of active growth and development of the main field vegetables is considered to correspond to the stable transition of the average daily air temperature through 10 ∘ C. The mean dates of the start and end of the growing season throughout the Elbe lowland for  mean ≥ 10 ∘ C are April 21 and October 7, respectively, and the mean growing season length is 169 days (Table1).The earliest date of BGS ≥ 10 ∘ C was observed in 1998(30 March), an advance of 22 days caused by strong positive temperature anomalies (up to +6.0 ∘ C).The onset of spring was also very early this year throughout Central Europe[7].The reverse was observed in 1980, with the latest date of BGS ≥ 10 ∘ C observed (May 22, +31 days compared with the long-term average).Moreover, an early end of the growing season in the fall of 1980 resulted in extreme reductions in the GSL ≥ 5 ∘ C (−13 days), GSL ≥ 10 ∘ C (−31 days), and GSL ≥ 15 The period of  mean ≥ 15 ∘ C corresponds to the beginning of the transplanting of thermophilic vegetables.The mean dates of the start and end of the growing season for  mean ≥ 15 ∘ C are May 21 to September 7, and the mean length of the growing season is 109 days.The earliest date for BGS ≥ 15 ∘ C occurred in 1962 (21 April), and the latest date occurred in 1974 (26 June).The earliest end of the growing season occurred in 1978 and 1986 (19 August), earlier than the average by −19 days (Tables1 and 2).The latest end occurred in 1966(30  September), later than the average by +23 days.The greatest lengthening of the growing season, +42 days compared with the long-term mean, occurred in 2000, and the greatest shortening of the growing season, −42 days compared with the long-term mean, occurred in the cool year of 1965.In 2000, strong positive temperature anomalies of 3.5 ∘ C to 5.0 .2. Temporal Variability of the Beginning, End, and Length of Climatological Growing Seasons 3.2.1.Anomalies in the Start, End, and Length of the Growing Season.Regional anomalies in the start, end, and length of the growing season for the three threshold temperatures in the Elbe River lowland for 1961-2011 are shown in Table2.Growing seasons with stronger negative and positive anomalies were classified into two categories: (1) years with advances and/or increases in the growing season (i.e., an early start and/or late end) and (2) years with delays and/or decreases in the growing season (i.e., a later start and/or early end).During the second half of the 20th century and the first decade of the 21st century, the years with the largest anomalies in terms of an early start of the growing season (BGS ≥ 5 ∘ C and a later end of EGS ≥ 5 ∘ C, together with a late end of the growing season, occurred in 1961, 1990, and 2011 (Table2).However, the years 2010 and 2011 had early fall frosts that occurred during the harvest period for root vegetables (e.g., Apium graveolens).
[7,29]ater compared with the long-term mean, as a result of the strong positive  mean anomalies that occurred from September to October.The year 2006 had also the longest GSL ≥ 5 ∘ C.3.1.2.The Main Growing Season: mean ≥ 10 ∘ C. ∘ C (−30 days) compared with the long-term average (Table2).The length of the growing season significantly affects the vegetable sector.Its extension in the spring allows for early sowing/planting of crops, and its extension in the fall has a positive effect on the production quality, allowing a gradual harvest.The earliest date of EGS ≥ 10 ∘ C occurred in 1972 (September 16), and the latest date occurred in 2000 (October 28).Consequently, the longest and shortest growing seasons occurred in 2000 (195 days) and 1972 (136 days), respectively.3.1.3.The Summer Growing Season: mean ≥ 15 ∘ C. ∘ Cin the daily mean temperatures in April led to an advance in BGS ≥ 15 ∘ C of 28 days.In the same year, a significantly longer growing season was recorded, not only in terms of GSL ≥ 15 ∘ C but also in terms of GSL ≥ 5 ∘ C (+22 d) and GSL ≥ 10 ∘ C (+26 d).The growing season was extremely long that year throughout Europe[7,29].3

Table 3 :
[35]rved long-term changes (days per year) in growing season parameters over the Elbe River lowland.Marginally significant (0.01 <  ≤ 0.05); * * significant (0.001 <  ≤ 0.01); * * * highly significant ( ≤ 0.001).andlaterdates, respectively.However, from 1981 to 2000, the earlier BGS ≥ 5 and 10 ∘ C was more dominant than the delayed EGS ≥ 5 and 10 ∘ C, whereas from 2001 to 2010, for the growing season defined in terms of  mean ≥ 15 ∘ C, the later EGS (5.6 days decade −1 ) was more significant than the earlier BGS (0.1 days decade −1 ).Similar results at the European level were reported by[35], with the delay in EGS by 8.2 days being more significant than the advance in BGS by 3.2 days during 2000-2008.
3.3.Observed Spatial Variability of Start, End, and Length ofGrowing Seasons.The GS parameters exhibited clear spatial variability.For clarity of presentation, the locations of specific fields and the Elbe River across the lowlands were added to the maps.Over the past 51 years, the GS ≥ 5 ∘ C started between March 11 and 26, and a start before March 11 was calculated for the northeastern Prague plateau (Figure

Table 4 :
Projected dates of the start, end, and length of the regional growing season for two future periods (2021-2050 and 2071-2100) based on two different RCMs by A1B scenario run over the Elbe River lowland. mean ≥ 5 ∘ C  mean ≥ 10 ∘ C  mean ≥ 15 ∘ C

Table 5 :
Projected changes in the dates of the start, end, and length of the area-averaged growing season for two future periods (2021-2050 and 2071-2100) over the Elbe River lowland. mean ≥ 5 ∘ C  mean ≥ 10 ∘ C  mean ≥ 15 ∘ C ∘ C) started 10.7 days (7.6 d and 3.1 d) earlier and ended 4.6 days (3.6 d and 4.1 d) later, yielding an increase in GSL of 15.3 days (11.2 d and 7.2 d) from 1961 to 2011.The delay in EGS was found to be a more important factor than the advance in BGS in regulating the GSL changes during the 1990s and 2000s at the threshold of  mean ≥ 15 ∘ C. (iii) Under projected future climate conditions and temperature thresholds, the dominant dates of BGS and EGS for the entire study region are projected to be significantly advanced and delayed, respectively, compared with the current climate.A climate warming scenario suggests lengthening of the GSL in the coldest areas of the study region to the level of the warmest areas in the current climate.