Clim. Past 12, 1421-1434, 2016 www.clim-past.net/12/1421/2016/ doi:10.5194/cp-12-1421 -2016 © Author(s) 2016. CC Attribution 3.0 License. o Climate 1 of the Past I (KU April-August temperatures in the Czech Lands, 1499-2015, reconstructed from grape-harvest dates Martin Možný1, Rudolf Brázdil2'3, Petr Dobrovolný23, and Miroslav Trnka3'4 ^zech Hydrometeorological Institute, Doksany Observatory, 41182 Doksany, Czech Republic institute of Geography, Masaryk University, Kotlářská 2, 61137 Brno, Czech Republic 3Global Change Research Institute AS CR, Bělidla 986/4a, 603 00 Brno, Czech Republic 4Faculty of Agronomy, Mendel University in Brno, Zemědělská 1, 603 00 Brno, Czech Republic Correspondence to: Martin Možný (martin.mozny@chmi.cz) and Rudolf Brázdil (brazdil@sci.muni.cz) Received: 9 February 2016 - Published in Clim. Past Discuss.: 10 February 2016 Revised: 1 June 2016 - Accepted: 2 June 2016 - Published: 5 July 2016 Abstract. Viticulture has long been essential to the commercial and social well-being of parts of the Czech Lands (now the Czech Republic), and detailed records have been kept for centuries of the timing and relative success of the grape crop. Using such documentary data from the Bohemian winegrowing region (mainly northwest of the capital, Prague), series of grape-harvest dates (GHDs) were created for the 1499-2015 period. Because the link between harvest dates and temperatures is strong, GHD series, together with instrumental mean temperature series starting in 1801, were used to reconstruct mean April-August temperatures for the region from 1499 to 2015. Linear regression (LR) and variance scaling (VS) methods were used for calibration and compared in terms of explained variance and their ability to capture extreme values. It emerged that LR does not significantly underestimate temperature variability. However, VS shows far greater capacity to capture extremes. GHDs explain 64 % of temperature variability over the full calibration period. The 1986-2015 period was identified as the warmest 30-year period of the past 514 years, an observation consistent with recent global warming. The highest April-August temperatures appeared in a reconstruction for the year 1540, which was warmer than the next two very warm, and far more recent, seasons in 2003 and 2015. The coldest period occurred at the beginning of the 20th century (1900-1929). The series reconstructed for the Czech Lands is in close agreement with other (central) European reconstructions based on other proxies. The series created here makes an important contribution to a better understanding of long-term spatiotemporal temperature variability in central Europe. 1 Introduction A specific set of climatic circumstances is crucial to the achievement of complete vine grape maturation, and climate variability determines year-to-year differences in the yield and quality of the wine produced (Jones and Hell-man, 2003). Temperature and solar irradiance are particularly critical variables, since they have a direct effect on the length of the growing season, phenological stages, grape yields, and the synthesis and accumulation of sugars, organic acids, polyphenols, vitamins, and aromatic compounds in the berries (Jones et al., 2005; Keller, 2010). The warming trend in recent climate change has accelerated ripening and induced changes in grape composition (Jones et al., 2012). Documentary data may provide important proxies for the study of past climate evolution (Brázdil et al., 2005, 2010; Jones et al., 2009). Grape-harvest dates (GHDs) are well suited to the reconstruction of interannual variations in spring and summer temperatures over large parts of Europe. The prevailing temperatures during the stages before flowering and veraison (color and profound composition changes) have been established as the most influential upon GHDs (Chuine et al, 2004; Garcia de Cortázar-Atauri et al, 2010). GDHs have been used for the reconstruction of temperature series in France (Le Roy Ladurie and Baulant, 1980; Chuine et al, 2004; Le Roy Ladurie, 2005; Menzel, 2005; Etien et al, 2008; Gamier et al, 2010), Switzerland (Pfister, 1981, 1984; Burkhardt and Hense, 1985; Meier et al, 2007; Wetter and Pfister, 2013), Germany (Glaser and Hagedorn, 1991), Austria (Strommer, 2003; Maurer et al, 2009), Italy (Mariáni Published by Copernicus Publications on behalf of the European Geosciences Union. 1422 M. Možný et al.: April-August temperatures in the Czech Lands et al., 2009), Hungary (Kiss et al., 2011), and Luxembourg (Molitor et al, 2016). Daux et al. (2012) created an open-access dataset of GHD data from different parts of Europe. However, when addressing GHD series, a number of uncertainties should be taken into account, largely associated with changes in varieties, viticulture techniques, climatic impacts on vines on interan-nual and decadal scales, and missing data (Meier et al., 2007; Garcia de Cortázar-Atauri et al., 2010; Krieger et al., 2011). Wine production in the Czech Lands (since 1993 the Czech Republic) is limited to two regions - southern Moravia and northwestern Bohemia. In an analysis of the long-term viticultural data available in the Czech Lands, Brázdil et al. (2008) concluded that temperature reconstruction prior to AD 1800 was rendered difficult by insufficient data coverage. However, newly acquired records of grape-harvest dates from the Bohemian wine-growing region have now allowed a longer-term chronology of GHDs to be created, and it makes extended temperature reconstruction possible. In the light of this, the aim of this study is a reconstruction of April-August temperatures for the Czech Lands over the 1499-2015 period. This contribution describes the area studied, and outlines basic data and methods used for reconstruction (Sects. 2 and 3), presents a new Czech temperature reconstruction along with its fluctuations and some comparisons with other temperature reconstructions in Central Europe (Sect. 4), and discusses the results achieved (Sect. 5). The final section contains some concluding remarks. 2 Area and data The Czech Republic is located in central Europe and is characterized by a moderate, humid climate, with the year divided into four distinct annual seasons (Tolasz et al., 2007). The Bohemian wine-growing region is comparatively small, situated near the capital Prague in the northwest of the country, near the Vltava, Ohře and Elbe rivers and is protected by the Krušné hory, the Doupovské hory, and the Czech central mountains to the northwest, all of which combine to create a rain-shadow effect. It features protected slopes facing south, southwest and southeast at elevations between 170 and 260 m. The vineyards are located around the city of Prague and the towns of Litoměřice, Roudnice nad Labem, Mělník, Most, Slaný, Kutná Hora, Louny and Karštejn (Fig. la). In the environmental conditions prevailing in recent times in the Czech Republic, the significant phenophases of vine (Vi-tis vinifera) occur on average on the following days: leaf-bud burst on 20 April-8 May, beginning of flowering on 5-15 June, end of flowering on 15-26 June, and softening of berries on 1-29 August (Hájková et al., 2012). The mean annual temperature reaches 8.7-8.9 °C (15-15.3 °C for April-August season) and the mean annual precipitation total is 480-540 mm (255-280 mm for April-August season) for the 1961-2000 period in the area studied. The entire Czech wine 1500 1550 1600 1650 1700 1750 1800 1850 1900 1950 2000 Year Figure 1. (a) Map of the Bohemian wine-growing region; (b) the number of sites of grape-harvest observations per year for the 1499-2015 period in individual years. region falls into Region I according to the standard growing degree days (GDDs), representing the climate type "Cool" according to the Huglin Index (HI) and cool climate maturity according to the average growing season temperature index (GST). These indices are used to characterize the climatic conditions suitable for specific wine-grape cultivars and the general types of wine that can be produced (Jones et al., 2010). Legend has it that the first Czech vineyard, near the town of Melnik, was established by St. Wenceslas (907-935), King of Bohemia, considered the patron saint of the Bohemian wine-growing region (Fig. 2a). In the year 1358, Emperor Charles IV (1316-1378) contributed an essential stimulus when he decreed that vineyards should be established wherever suitable. He brought Burgundy grape varieties to Bohemia and gave a 12-year exemption from taxes to those who started vineyards (Kraus et al., 2005). For 1845 to the present, the relevant GHDs were taken from the PHENODATA database of the Czech Hydrom-eteorological Institute (CHMI), which includes systematic phenological observations from several hundred sites across the territory of the recent Czech Republic (Svitakova et al., 2005). For the period prior to this, GHDs were obtained from chronicles, documentary records kept by town authorities, Clim. Past 12, 1421-1434, 2016 www.clim-past.net/12/1421 /2016/ M. Možný et al.: April-August temperatures in the Czech Lands 1423 Figure 2. (a) Oil painting of St. Wenceslas in a vineyard, 18th century; (b) flag of the guild of wine cultivators, 1762; (c) insignia of the guild of wine cultivators, 1676 (from the collections of the Mělník Regional Museum; photo M. Možný). documentary records from aristocratic sources (in particular from the Lobkowicz and Nostitz families), ecclesiastical documentary records (Archdiocese of Prague, Order of Cistercians), farming calendars, personal diaries, farming records, and reports. Among the earliest records are those from the towns of Litoměřice (Smetana, 1978) and Mělník (Böhm, 1892) and financial records kept for Bohemian viticulture (e.g. Schams, 1835; Šimáček, 1891; Honzík, 1897). For example, documentation for the year 1602 records the following: "Beginning of grape harvest on 14th October, harvest lasted for 9 days, good wine" (Lobkowicz records); "Grape harvest on 15th October, grape-picking lasted 8 days, wages paid for 6 men and 42 women" (Mělník records); and "Harvest of sacramental wine on 16th October" (Archdiocese of Prague records). Each town had a wine growers' guild that organized the beginning of the grape harvest. For example, the Mělník town guild started in 1542; numerous relics of its work survive, including written records (Fig. 2b-c). Records of GHDs for the Bohemian wine-growing region were thus ascertained for the 1499-2015 period (Fig. lb). In total, 5431 harvest dates from 61 sites were collected. The new Czech GHD series is made up of grape-harvest dates recorded in contemporary manuscripts (1499-1844) and those observed by institutional phenological networks (1845-2015). Documentary data earlier than AD 1584 have been redated to the Gregorian calendar by adding 10 days to the previous Julian calendar. The majority of observations go back to 1921, whereas relatively fewer data are available before 1700. There are no temporal gaps; at least one entry was found for each year. While up to 1844, timing differences in GHD between individual places averaged 1-3 days, a further increase in the number of sites in the phenological network (see Fig. lb) led to differences of 1-8 days. The mean GHD time series for the 1499-2015 series was taken as a median value for all site data in the given year. Mean Czech temperatures for 1800-2010 (Brázdil et al., 2012a, b), further updated to 2015, were used as a target climatological series for temperature reconstruction during the pre-instrumental records for GDHs. To put the Czech April-August temperature reconstruction into a wider European context, the following temperature series were selected for comparison (only those starting in the 16th century): 1. AMJJA (April-August), central Europe, documentary data (Dobrovolný et al., 2010); 2. JJA, central Europe, various proxies (Luterbacher et al., 2004) ; 3. JJA, Low Countries, documentary data (van Engelen et al, 2001); 4. AMJJA, Swiss Plateau, grape-harvest dates (Meier et al, 2007); 5. AMJJA, Burgundy, France, grape-harvest dates (Chuine et al., 2004); 6. MJJ, Vienna and Klosterneuburg region, Austria, grape-harvest dates (Maurer et al., 2009); 7. MJJ, Közseg, Hungary, grape-harvest dates (Kiss et al., 2011); 8. AMJJ, Switzerland, grape-harvest dates (Wetter and Poster, 2013); 9. MAMJJ, northern Switzerland and southwestern Germany, grain-harvest dates (Wetter and Poster, 2011); 10. MAMJJ, Czech Lands, winter wheat-harvest dates (Možný etal, 2012); 11. JJA, Alps, multiproxies (Trachsel et al., 2012); 12. AMJJAS, western Europe, multiproxies (Guiot et al., 2005) ; 13. JJA, Europe, multiproxies (Luterbacher et al., 2016); 14. JJA, western and central Europe, tree-ring widths (Bünt-gen etal, 2011); 15. JJAS, Alps, tree-ring widths (Büntgen et al, 2006); 16. JJ, central-eastern Alps, tree-ring widths (Wilson and Topham, 2004). www.clim-past.net/12/1421/2016/ Clim. Past, 12, 1421-1434, 2016 1424 M. Možný et al.: April-August temperatures in the Czech Lands 3 Methods The mean GHD series for the Bohemian wine-growing region was further used as a predictor for temperature reconstruction for the Czech Lands in the 1499-2015 period. Two approaches were taken to calibration: simple linear regression (LR) and variance scaling (VS). LR uses least squares estimation to minimize mean square error (MSE) between target data and reconstructed temperature. However, regression-based reconstructions usually underestimate the variability of reconstructed temperatures to a significant extent (Esper et al., 2005). Moreover, their capacity for capturing extreme values is limited. In contrast to regression, variance scaling (matching) rescales proxy data to fit target measured temperatures by providing them with the same mean and variance over the calibration period (McCarroll et al., 2015). Thus, VS reconstruction usually captures the full range of measured temperatures and also better reproduces extreme values. However, a higher MSE is a drawback of VS in comparison with regression. Two statistical measures were employed herein to allow an objective comparison of the two approaches above (McCarroll et al., 2015). In linear regression, the squared correlation coefficient (r2) compares MSE based on regression estimates with MSE based on the mean (climatology). Similarly, the equivalent variance explained (r2s) compares MSE based on variance scaling with MSE based on the mean (climatology). The r2s has a straightforward interpretation - if its value falls below zero, variance-scaled reconstruction has less skill than simply using the mean values of data over the calibration period. Moreover, differences between r2 and r2s measure the degree of loss in variance accounted for by VS-based reconstruction compared with LR-based reconstruction. Extreme-value capture employs a similar logic, comparing how many extreme years (defined, for example, as 10 % of the highest and 10 % of the lowest measured temperatures in the full calibration period) are captured by either LR-based or VS-based reconstructions. The generally lower variability of regression-based reconstruction also means that its ability to capture the same extreme years is lower, while variance scaling may improve the ability to identify the same extreme years that occur in measured data. The extreme-value capture approach tests whether the number of identified extremes in a reconstruction is significantly higher than it would be by chance alone. For more details of the calculation of equivalent variance and extreme-value capture, see McCarroll et al. (2015). Final calibration was performed for the full period of overlap (1811-2010) - for verification purposes divided into early (1811-1910) and late (1911-2010) sub-periods - and the two calibration-verification schemes for LR and VS deployed. Calibration was done on one of the sub-periods and verification undertaken on the other. The whole process was then repeated with calibration-verification sub-periods switched. The quality of calibration was evaluated by means of the square of Pearson's correlation coefficient (r2), which measures the degree of common variance explained, and with Durbin-Watson test statistics (DW), which diagnose the degree of first-order autocorrelation within the regression residuals; values between 1.5 and 2.5 are generally acceptable, while those outside this range indicate problems with reconstructing multi-decadal variations (Montgomery et al., 2012). DW is calculated as DW = T(et-ei-i)2 1=2_ n ' 1=1 (1) where et is the residual associated with the observation at time t. Verification was based on r2, reduction of error (RE), and the coefficient of efficiency (CE). RE indicates whether a reconstruction provides a better estimate of temperature variability than simply using the mean value of target temperatures in the calibration period. CE is similar to RE, but it tests reconstruction skill against the mean value of target temperatures in the verification period. RE is calculated as RE = 1 "žli = l(Xi ~ Xc)2 (2) where and are target and reconstructed temperatures for the verification period, i.e. the mean value of target temperatures in the calibration period CE is similar to RE, but it tests reconstruction skill against the mean value of target temperatures in the verification period (Xu): CE= 1 X!i=i(Xi— x«)2 (3) Both RE and CE provide values between 1 and negative infinity. Positive RE and CE values indicate that the method examined has some potential for reconstruction, and the result is better than simply using the mean of a given calibration-verification period as a "reconstruction". Final LR-based reconstruction was completed with uncertainty estimates defined as 2 standard errors (SE) of estimates from the corresponding regression model, approximating to 95 % confidence limits. The series reconstructed were compared with other European reconstructions using Pearson correlations for the whole common period as well as with running 31-year correlations to reveal periods with significant agreement in reconstructed temperatures. The significance level was estimated by means of the Monte Carlo simulation to oversee the effect of multiple comparisons (Livezey, 2010). Running correlations were calculated for 1000 permutations of the time series compared, and the highest correlation for each trial was established. Finally, the 95th quantile Clim. Past, 12, 1421-1434, 2016 www.clim-past.net/12/1421 /2016/ M. Možný et al.: April-August temperatures in the Czech Lands 1425 1500 1550 1600 1650 1700 1750 Year Figure 3. Fluctuations in median grape-harvest day (in days from 1 September) for the Bohemian wine-growing region in the 1499-2015 period. Smoothed by Gaussian filter over 30 years. í 06 ..JIHL ■ * Month Figure 4. Pearson correlation coefficients between grape-harvest dates and Czech temperatures in the 1811-2010 period for months and seasons; broken horizontal line indicates significant correlations for one-tail test (p = 0.05, N = 200). from the 1000 values of the highest correlation coefficients determined the critical value at p — 0.05 for the one-tailed test. 4 Results 4.1 Czech temperature reconstruction based on grape-harvest dates Fluctuations of GHDs in the Bohemian wine-growing region during the 1499-2015 period indicate that earlier harvest dates (9 days earlier than the 1961-1990 mean) were found in 1986-2015 (Fig. 3). On a decadal scale, the earliest harvests occurred in 2000-2009 (mean, 28 September) and the latest in 1814-1823 (29 October). The earliest GHD of all was recorded in 1540 (11 September) and the absolute latest in 1919 (2 November). GHDs as early as 30 September occurred only very rarely: 10 times in the 16th century; twice in the 17th century; 5 times in the 18th century; twice in the 19th century; and 13 times in the 20th and 21st centuries. GHDs later than the second half of October occurred more often: 29 times in the 16th century; 49 times in the 17th century; 52 times in the 18th century; 45 times in the 19th century; and 43 times in the 20th and 21st centuries. The beginning of the grape harvest depends on the weather patterns of the preceding months and the closest relationship with GHDs in the 1811-2010 period appeared for mean April-August temperatures (r2 = 0.64). Based on the results of the correlation analysis (Fig. 4), GHDs were calibrated to mean AMJJA Czech temperature series (Brázdil et al., 2012a, b). Summaries of the calibration-verification processes for early (1811-1910) and late (1911-2010) sub-periods and for the full calibration period (1811-2010) and for LR and VS methods of reconstruction appear in Table 1 and Figs. 5-6. The LR statistical measures applied confirm high reconstruction skill for GHDs in both sub-periods. The AMJJA temperatures explain 51 and 76 % of GHD variability for early and late calibration, respectively. Results of the Durbin-Watson test indicate no problems with autocorrelation in residuals, and standard error does not differ signifi- 1850 1900 1950 2000 "(b) Verilicalion Calibration i 1850 1900 1950 2000 "(C) \t/^^~^^^ 1850 1900 1950 2000 Figure 5. Comparison of measured (1) and reconstructed (2) mean Czech April-August temperatures: anomalies with respect to the 1961-1990 period. Linear regression method of reconstruction was performed for (a) early (1811-1910) and (b) late (1911-2010) calibration periods; (c) running correlations (window 21 years) between mean Czech April-August temperatures and grape-harvest dates during the 1811-2010 period. cantly for either of the sub-periods. Validation of these calibration results shows that RE and CE are highly positive both for early and late sub-periods and indicates that GHDs have a high potential for temperature reconstruction. Full calibration with the LR model explains 64 % of common variability between GHDs and AMJJA temperatures. The VS method, verified by the same RE and CE metrics, also confirms high reconstruction skill for GHDs. The two metrics give lower results for late verification than those from the LR method (Table 1), but they remain highly positive. Equivalent variance explained (r2s) for the VS method and the full calibration period is 0.61, which, compared with r2 of the LR model, indicates that the scaling (variance matching) applied to GHD data involves only 6 % of lost variance explained compared with LR. www.clim-past.net/12/1421/2016/ Clim. Past, 12, 1421-1434, 2016 1426 M. Možný et al.: April-August temperatures in the Czech Lands Table 1. Summary of calibration and verification statistics used for April-August temperature reconstruction based on grape-harvest dates for the Czech Lands: r2 - explained variance; SE - standard error of estimate; DW - Durbin-Watson test; RE - reduction of error; CE -coefficient of efficiency. RE and CE were calculated for linear regression (LR) and variance scaling (VS). Sub-periods Calibration statistics | Verification statistics r2 SE DW | r2 RELR CELR REVs CEVs Early calibration (1811-1910)/ 0.51 0.52 1.9 0.76 0.66 0.53 0.71 0.60 late verification (1911-2010) Early verification (1811-1910)/ 0.76 0.47 1.8 0.51 0.51 0.21 0.50 0.19 late calibration (1911-2010) Full calibration (1811-2010) 0.64 0.54 1.9 0.64 - 1850 1900 1950 2000 1850 1900 1950 2000 Year Figure 6. Comparison of measured (1) and reconstructed (2) mean Czech April-August temperatures: anomalies with respect to the 1961-1990 period. Variance scaling method of reconstruction was performed for (a) early (1811-1910) and (b) late (1911-2010) calibration periods. p ■ 1 -4 3 • 2 - (a) Regression 9 of 21 ■Já- 1 • 0 A fuK A a A.A ,AA IffV* 1 -2 - [-4—8of20 - Measured r | ¥ " Reconstructed 0 50 100 150 200 3 -2 ' (b) Scaling 17of21 1 • 0 ' 1 -2 - ! 1«—11 of 20 - Measured - Reconstructed 0 50 100 150 200 Rank order Figure 7. Comparison of measured (continuous blue line) and reconstructed temperatures using regression (a - green line) and scaling (b - red line) presented in rank order; dashed boxes indicate the upper and lower 10 % portions of measured temperatures and numbers of these extremes reproduced in reconstructions. Much greater differences between the two reconstruction methods become evident in a comparison of the overall variability of reconstructed values with those of target data in the full calibration period. AMJJA temperature anomalies varied from -2.2 °C (1919) to 2.7 °C (2003) in the 1811-2010 period. The corresponding values for the LR model are — 1.9 °C (1919) and 2.5 °C (2000), respectively, and they showed that the regression reconstruction is biased towards the mean and underestimates the true variability of the target data. According to VS reconstruction, AMJJA temperature anomalies fluctuate from -2.4°C (1919) to 3.2°C (2000). Thus, the range of the VS-reconstructed temperatures exceeds those measured in the common period of overlap (1811-2010). Finally, both reconstructions were compared in terms of their ability to capture extremely low or high temperatures (Fig. 7). The same period (1811-2010) was examined and the 10% of years with the lowest (20 years) and the highest (21 years) temperatures were defined as extreme. The maxima contain one more year than the minima because 2 years in the former reached equally high temperatures. The measured and reconstructed temperatures were arranged in or- der of rank, and a comparison was made to establish how many reconstructed extremes corresponded to their measured equivalents (Fig. 7). While the LR reconstruction captured 8 extremely cold years and 9 extremely warm years, the VS reconstruction proved able to reproduce 11 cold and 17 warm extremes. All the figures for the two reconstructions exceed the critical thresholds for the extreme-value capture test (p<0.05), after McCarroll et al. (2015). However, reconstruction based on scaling always performs far better than regression. 4.2 Fluctuations in Czech April-August temperatures Annual and multi-decadal variations in the Czech April-August temperatures since AD 1499, reconstructed by the two methods, appear in Fig. 8. Whereas certain differences on the high-frequency (annual) scale follow from the methods used and are described in the previous paragraphs, both approaches provide almost identical results on decadal and multi-decadal scales. This further confirms the high potential of GHD series for temperature reconstruction and also confirms the degree of robustness of reconstruction results. Clim. Past, 12, 1421-1434, 2016 www.clim-past.net/12/1421 /2016/ M. Možný et al.: April-August temperatures in the Czech Lands 1427 1500 1550 1600 1650 1700 1750 1800 1850 1900 1950 2000 1500 1550 1600 1650 1700 1750 1800 1850 1900 1950 2000 Year Figure 8. Measured (1) and reconstructed (2) mean Czech April-August temperatures in the 1499-2015 period using linear regression (a) and variance scaling (b); temperatures indicate anomalies with respect to the 1961-1990 period. Thick blue and black lines represent measured and reconstructed temperatures smoothed by Gaussian filter for 30 years; grey area approximates 95 % confidence interval constructed as two units of standard error (SE) of the estimate from the linear regression. Certain irregular variations in warmer and colder periods occurred during the 1499-2015 period. Significantly warmer periods (based on unsmoothed data) were found for the years 1516-1559 and 1971-2015 (Fig. 8). The first half of the 16th century was exceptionally warm, as indicated by early grape harvests. Four cycles of particularly chilly periods, beginning around 1649, 1740, 1805, and 1913, each separated by slightly warmer intervals, are detectable in reconstruction. In contrast, the coldest periods were recorded around 1820, 1840, and 1920. The warmest 30-year periods occurred in 1986-2015 (0.7 °C above the mean for the 1961-1990 reference period) and the coldest in 1900-1929 (-0.5 °C; both values based on LR reconstruction). Temperature fluctuations show great interannual and interdecadal variability and a positive trend for the instrumental part of the series from the 19th century onwards, particularly pronounced since the 1970s. 4.3 Czech temperature reconstruction in the European context The mean Czech April-August temperature reconstruction may be further verified by comparison with other Czech and European temperature reconstructions, as listed in Sect. 2. Such a comparison may be influenced by the distances between regions or areas for which the reconstructions are representative; in general, correlation coefficients decrease with increasing distance (see Fig. 5 in Dobrovolný et al., 2010). Direct comparison is also complicated by variations in sea- sonality, i.e. the selection of months used for averaging in individual reconstructions. Any common variability in reconstructed series was investigated by running correlations (31-year window) with 15 different series (Fig. 9). A rough categorization into five groups may be made, according to the dominant types of proxy: documentary-based, grape-harvest dates, grain-harvest dates, multiproxy, and tree rings. For Czech territory, GHD reconstruction may be correlated with mean Czech March-June temperatures reconstructed from winter wheat-harvest dates (WWHDs) (Možný et al., 2012), giving an overall correlation of 0.52 (p <0.05). The relationship between these two Czech reconstructions becomes weaker in the light of several considerable drops in running correlation coefficients around 1600, 1675, 1730, 1770, 1900, and 1980. The earliest WWHDs occurred in 2007 (4 July), 1516 (6 July), and 1521,1535,1539, and 1540 (7 July). Grain-harvest days in 1516 were affected by an extremely warm winter (Pejml, 1966; Brázdil et al., 2013). It is remarkable that the grain harvest in 1540 was, like the grape harvest, extremely early. The reconstructed Czech April-August temperatures correlated best (r =0.71) with the central European temperature series reconstructed from documentary data in Germany, Switzerland, and the Czech Lands compiled by Dobrovolný et al. (2010). These relationships become weaker in the light of several considerable drops in 31-year running correlations occurring around 1750 and 1900. Similar periods of lower coherency are also evident for the summer reconstruction based on temperature indices derived from documentary data for the more distant Low Countries (van Engelen et al., 2001) and for a gridded summer reconstruction for central Europe (Luterbacher et al., 2004) which, however, employed more types of proxy than documentary evidence alone. Correlation coefficients with tree-ring-based temperature reconstructions are distinctly lower; however, two of them (Wilson and Topham, 2004; Buntgen et al., 2006) represent mountain positions in the Alps. 5 Discussion The Bohemian wine-growing area lies among the northernmost parts of the European viticulture region, between the latitudes 49.9 and 50.6° N, and the quality of its product is thus significantly vulnerable to temperature variability between seasons and years. A cold growing season results in a below-mean vintage, with unripe and inferior grapes of low quality, while warmer growing seasons improve the quality of grapes considerably. In terms of minimum temperatures, spring frosts that damage developing shoots and frosts after bud burst, reducing the season's overall yield, are among the extreme events that have detrimental effects. Winter frosts below —12 °C may severely damage or even kill fully dormant vines (Kraus et al., 2005); the greater frequency of www.clim-past.net/12/1421/2016/ Clim. Past, 12, 1421-1434, 2016 1428 M. Možný et al.: April-August temperatures in the Czech Lands hard winters in the 18th and 19th centuries adversely affected both the revenues and profitability of viticulture in the Bohemian wine-growing region (Pejml, 1966). At the beginning of the 16th century, the Bohemian vineyards covered around 5500 ha. With deteriorating weather conditions, this area steadily declined to 3300 ha in the mid-18th century, falling to 700 ha at the end of the 19th century, a figure that holds to the present day (Kraus et al., 2005). In the context of recent climate change (global warming), conditions for growing wine in the Bohemian area are improving, with positive effects on the quality of the wine. The main grapevine varieties (Pinot, Traminer, Gouais, Riesling) have been cultivated in the Bohemian winegrowing region for the last 6 centuries (Sotolář, 2009). The Pinot varieties (Pinot noir, Pinot blanc, and Pinot gris) have always been dominant in the region. In the past 30 years, the Müller Thurgau and St. Laurent varieties have also expanded. The long-term prevalence of Pinot varieties in the area more or less excludes the possible influence of changing varieties on temperature reconstruction. At first sight, the LR and VS methods used in this analysis appear to have provided a very high level of reconstruction skill, as summarized in Table 1. However, closer examination of Figs. 5 and 6, where reconstructed and target temperatures for early and late calibration periods are compared, reveals certain inconsistencies. When the earlier period (1811-1910) is used for calibration, reconstructed April-August temperatures clearly underestimate those measured in the most recent decades. This is more pronounced for the LR method (Fig. 5) than for the VS method (Fig. 6). This probably relates to the fact that there has been a highly significant positive temperature trend since the 1970s, whereas the calibration period (1811-1910) shows lower temperatures, lacking any long-term positive trend. This anomaly is a common feature of regression-based methods, since they are able to reproduce well only that which they have "experienced". Moreover, differences between the measured and reconstructed temperatures in recent decades are also partly related to the fact that regression-based methods suppress the variability of reconstructed values - a feature discussed in Sect. 3. The effect described here is partly eliminated by the final reconstruction being based on calibration employing the full overlap period (1811-2010). In spite of this, lower reconstructed temperatures compared to those measured are evident in Fig. 8, especially in smoothed series for recent decades. This feature demonstrates further the exceptionality of recent global warming on a central European scale over the last 500 years. Comparison of Czech April-August temperatures with other proxy-based temperature reconstructions reveals possible inconsistencies in the first part of the 16th century, as well as around 1650,1750, and 1900 (Fig. 9). The first half of the 16th century, in response to warmer climate, saw particularly consistent early starts for the grain harvest in the Czech Lands (Možný et al., 2012). This warming allowed a northward extension of the Bohemian wine- and hop-growing re- 1500 1550 1600 1650 1700 1750 1800 1850 1900 1950 2000 - 4 (.60) -; 5 (-51) : 6(59^-7 (.56) - 8 (.58)1 1500 1550 1600 1650 1700 1750 1800 1850 1900 1950 2000 ](C) -9(46) - 10(52) *Y O 1500 1550 1600 1650 1700 1750 1800 1850 1900 1950 2000 V- 11 (.40)-12 (.46) 13(60)'^ 1500 1550 1600 1650 1700 1750 1850 1900 1950 2000 0 '■ Q-2 0 0 0 2 =) V V- 14(31) -15(26) V* 16(34) 1500 1550 1600 1650 1700 1750 1800 1850 1900 1950 2000 Year Figure 9. Running 31-year correlation coefficients between GHD-based reconstructed Czech AMJJA temperatures and selected European proxy reconstructions, (a) Documentary: (1) AMJJA Central Europe (Dobrovolný et al., 2010); (2) JJA central Europe (Luter-bacher et al., 2004); (3) summer Low Countries (van Engelen et al., 2001). (b) Grape-harvest dates: (4) AMJJA Swiss Plateau (Meier et al., 2007); (5) AMJJA Burgundy (Chuine et al., 2004); (6) MJJ Vienna and Klosterneuburg region (Maurer et al., 2009); (7) MJJ Közseg, Hungary (Kiss et al., 2011); (8) AMJJ Switzerland (Wetter and Poster, 2013). (c) Grain-harvest dates: (9) MAMJJ northern Switzerland and southwestern Germany (Wetter and Poster, 2011); (10) MAMJJ Czech Lands (Možný et al., 2012). (d) Multiproxies: (11) JJA European Alps (Trachsel et al., 2012); (12) AMJJAS western Europe (Guiot et al., 2005); (13) JJA Europe (Luterbacher et al., 2016). (e) Tree rings: (14) JJA western and central Europe (Büntgen et al., 2011); (15) JJAS Alps (Büntgen et al., 2006); (16) JJ central-eastern Alps (Wilson and Topham, 2004). Numbers in brackets are overall correlations for the common period in question. Broken horizontal lines represent critical values of correlation coefficient estimated using Monte Carlo simulation to oversee the effect of multiple comparisons and approximate 95 % confidence level for one-tailed test. Vertical columns A-D identify periods in which GHDs might have been influenced by the social and political processes discussed in the text: A - series of bad grape yields; B - the Thirty Years' War; C - the War of the Austrian Succession and the Seven Years' War; D - series of bad grape yields. gions and extensive cultivation of quite thermophilic plants, such as melons and almonds (Pejml, 1966). High summer temperatures in the Czech Lands during the first half of the 16th century are further confirmed by the results of temperature reconstructions for this century using series of temperature indices (Brázdil et al., 2013). Clim. Past, 12, 1421-1434, 2016 www.clim-past.net/12/1421 /2016/ M. Možný et al.: April-August temperatures in the Czech Lands 1429 The Czech April-August temperature reconstruction was based on a series of continuous GHDs for 1499-2015, using the innovative statistical tool of variance scaling to optimize the capture of extreme temperatures. This reconstruction indicates extremely warm years more readily than cold ones; the percentage of warm extremes confirmed from GHDs (61 %) is significantly higher compared to cold extremes (39 %). This finding is interesting in the light of comparison with a number of hydroclimate reconstructions (drought, precipitation) based on dendro-climatological data (Büntgen et al, 2011; Bronisz et al, 2012; Dobrovolný et al, 2015). As follows from Sect. 4.3, overall correlations between the new GHD-based temperature series and 16 different European proxy reconstructions are relatively high and statistically significant. However, those relations are not stable over time, and periods with evident drops in running correlations occur quite frequently. The most evident periods of lower shared variance refer to the same time intervals in most of the series compared. This indicates that there must be some common causes, examined below for the four most obvious periods (Fig. 9). The recent Czech reconstruction shows a temperature decrease in a critical period for wine production in the late 16th century (Fig. 9a) (Brázdil et al., 2013). Starting in the mid-15808, a number of years produced small amounts of wine, of very inferior quality. This trend also appears in Land-steiner (1999), who documents a sudden drop in wine production in central Europe. A series of bad grape yields started in 1585 in Switzerland, in 1586 in Württemberg (southwestern Germany), and in 1587 in Lower Austria and western Hungary. This decline was reflected in a significant drop in revenue derived from viticulture for the Austrian state. Moreover, wine cost more at the time and tasted less sweet, leading to a shift in consumer choice from wine to beer. GHDs in the Czech Lands around 1650 (see Fig. 3) were enormously influenced by the Thirty Years' War (1618— 1648) and its subsequent hardships. The war totally devastated the Czech Lands. The population decreased by a fifth (Wilson, 2009) and some historical reports estimate population declines as catastrophic as 80 % in certain regions (Bal-bfn, 1986). Moreover, GDHs were negatively affected by an increased frequency of late spring frosts, as has been documented for the 1640s and 1660s by Brázdil et al. (2008), for example (Fig. 9b). Temperatures decreased over the second half of the 17th century and the growing season shortened all over Europe (Hudson and Keatley, 2010). This period, corresponding to the Maunder minimum of solar activity (Eddy, 1976), is recognized as the coldest phase of the Little Ice Age (Grove, 2004; Matthews and Briffa, 2005). The period following 1740 was also characterized by intense sociopolitical conflicts associated with the War of the Austrian Succession (1740-1748) and the Seven Years' War (1756-1763), with many of the key field campaigns of the period taking place on Czech territory (Fig. 9c). During the war years, grapes were harvested at the earliest possible onset of maturity to preclude damage from enemy troops (Kilián, 2009). These periods of strife render the year-to-year variability of GDHs around 1650 and 1750 far less visible (Fig. 3). Although the socioeconomic consequences and loss of human life were not as acute as they had been ~ 100 years earlier, the countryside took the brunt of the damage. Peasant farms were burdened with increasing taxes, direct war damage, and army recruitment, while many owners simply left their farms. The non-sustainability of the situation was expressed in an unusually high number of peasant uprisings demanding the reduction of service-based tenure (socage) and consequent violent suppression of all rebellions. This constantly deteriorating situation culminated in what became known as the "hunger years" of 1770-1772, the result of adverse weather patterns, very bad harvests, and socioeconomic problems, when the population of Bohemia fell by about 10 % (Brázdil et al., 2001; Pfister and Brázdil, 2006). Decreases in running correlations between Czech April-August temperatures and other proxy-based central European temperature reconstructions in the mid-19th century and at the end of the 19th century (slightly earlier for the tree-ring group) appear to a greater or lesser extent in all five groups of the reconstructions compared (Fig. 9d). This indicates that there may be some problems in the quality of the GHD series for these periods. Warmer and wetter weather in August-September in the mid-18th century led to a higher incidence of fungal diseases on hops and vines, which influenced timing (early harvest) and yields (Možný et al., 2016). Effective synthetic pesticides and fungicides were not available until the early 1920s. The approach to vine growing and fermenting grapes changed in the late 19th century, affecting later harvest dates (Kilián, 2012). Vine growers experimented for several years with harvesting a few days later, since sunny weather in October may have helped augment sugar content. Comparison of Czech April-August temperatures reconstructed from GHDs and Czech March-June temperatures reconstructed from WWHDs is also worthwhile (Možný et al., 2012). Possible inconsistencies in the two reconstructions based on biophysical series are disclosed for the years around 1675, 1730, 1770, 1900, and 1980 (Fig. 9). Explanations may be related to recorded winter and spring frosts, which affected the harvest of winter wheat and grapes in different ways (Pejml, 1966, 1974; Možný et al., 2016). The earliest grape harvest recorded for the Bohemian wine-growing region during the 1499-2015 period occurred in 1540. The harvest may have been brought forward by very warm weather, but this could also have been a response to drought. Pfister and Brázdil (1999) identified the 1530s and 1550s as the warmest decades for summers in the 16th century in central Europe, and the 1530s to 1550s as the driest. In fact, 1540, with its very warm, dry weather, which has been documented and described as probably the most outstanding of its kind in the past 500 years in central Europe, was documented and described based on many documentary www.clim-past.net/12/1421/2016/ Clim. Past, 12, 1421-1434, 2016 1430 M. Možný et al.: April-August temperatures in the Czech Lands data (Glaser et al., 1999; Pfister, 1999; Brázdil et al., 2013). The temperature estimates for Switzerland from a great number of coherent qualitative documentary evidence confirm exceptional heat; Wetter and Pfister (2013) considered the spring-summer temperatures for 1540 probably more extreme in neighboring regions of western and central Europe than those of 2003 when between 22000 and 35 000 heat-related deaths were recorded across Europe (e.g. Beniston, 2004; Schär and Jendritzky, 2004; Stott et al, 2004; Fischer et al, 2007). More recently, Wetter et al. (2014) employed the term "megadrought" to the 1540 weather patterns in Europe, a coinage that led to some discussion in the tree-ring (Büntgen et al., 2015) and documentary data (Pfister et al., 2015) communities. Periods with statistically significant correlations were addressed by means of the Monte Carlo simulation to take into account the degree of freedom from autocorrelations. This approach confines the probability that significant correlations might be found only "by chance" (i.e. type I errors) to an acceptable level. The significant running correlations between our new GHD series and 16 existing reconstructions may be further affected by other factors, such as significant local trends in compared series. This may be the case for the second half of the 20th century, which shows high and largely significant running correlations for all the series compared. More obvious delimitation of the time intervals for which the new GHD reconstruction shows the best agreement with other proxy reconstructions (Fig. 9) may be expressed by the mean of running correlations (r-bar). In general, the best performance was found for the 1530s-1630s, the 1780s-1830s, and from the 1950s onwards. Higher correlations with documentary-based reconstructions compared with those based on tree rings underlines the importance of a multi-proxy approach. In this sense, direct comparison of the new GHD Czech reconstruction with the European multi-proxy JJA temperature reconstruction of Luterbacher et al. (2016) provides promising results. Running correlations are mostly significant, with the exception of only two relatively short periods (the 1650s and mid-18th century), which are discussed above (Fig. 9b and c). This result also demonstrates the distinctly spatially representative nature of the new Czech GHD reconstruction, which is not confined to the territory of the Czech Lands. Moreover, the GHD series also shows a high potential for the study of the long-term variability of droughts (Cook and Wolkovich, 2016; Možný etal, 2016). 6 Conclusions Central Europe is a region very rich in the types of documentary evidence that facilitate climate reconstructions for the pre-instrumental period (Brázdil et al., 2005, 2010). Such reconstructions, applying standard paleoclimatological approaches (e.g. Cook et al., 1994; Dobrovolný et al., 2009), are based either on series of temperature indices (Glaser and Riemann, 2009; Dobrovolný et al, 2010; Molitor et al, 2016) or biophysical series (Meier et al., 2007; Maurer et al., 2009; Kiss et al, 2011; Wetter and Pfister, 2011; Možný et al, 2012). In this context, April-August temperature reconstruction for the Czech Lands based on GHDs in the 1499-2015 period constitutes a further important contribution to the better understanding of long-term spatiotemporal temperature variability in central Europe, as well as in the broader European context. Advantages of the Czech temperature reconstruction include the very long overlap period (1811-2010) used for calibration-verification, a very high explained variability of 64 % in the instrumental period of overlap, the consistent dominance of Pinot varieties through time, and the stability of vineyard management throughout the period reconstructed. On the other hand, variable density of records over the past 500 years and some social bias related to long-term development of viticulture, together with the sociopolitical situation in the Czech Lands over time, may partly influence the reconstruction. In the context of discussions related to the exceedingly warm and dry year of 1540 (Wetter and Pfister, 2013; Wetter et al., 2014; Büntgen et al, 2015; Pfister et al., 2015), it is important that this new reconstruction confirms the year as the warmest April-August in 1499-2015 in the Czech Lands. Moreover, recent global warming (Stocker et al., 2013) finds a local reflection in the new Czech series, with 1986-2015 being the warmest period in the past 517 years. The Supplement related to this article is available online at doi:10.5194/cp-12-1421-2016-supplement. Acknowledgements. We acknowledge the CzechAdapt -System for Exchange of Information on Climate Change Impacts, Vulnerability and Adaptation Measures on the Territory of the Czech Republic project (EHP-CZ02-OV-1-014-2014), which was supported by a grant from Iceland, Liechtenstein, and Norway. R. Brázdil and R Dobrovolný acknowledge the support of the Grant Agency of the Czech Republic for the project no. 13-19831S. We extend our thanks to Richard Telford (University of Bergen) for his valuable comments and for the R-code estimation of significant running correlations. Tony Long (Svinošice) helped improve the English. Edited by: J. Luterbacher References Balbfn, B.: Wonders and wealth of Bohemia, Excerpt from Miscellanea historica regni Bohemiae by B. Balbfn, 1679, Panorama, Prague, 1986. Clim. Past, 12, 1421-1434, 2016 www.clim-past.net/12/1421 /2016/ M. Možný et al.: April-August temperatures in the Czech Lands 1431 Beniston, M.: The 2003 heat wave in Europe: A shape of things to come? An analysis based on Swiss climatological data and model simulations, Geophys. Res. Lett, 31, L02202, doi:10.1029/2003GL018857, 2004. Böhm, L.: Královské věnné město Mělník a okres Mělnický, The Royal Endowment Town of Mělníky and the Mělník District, Jelen, Mělník, 1892. Brázdil, R., Valášek, H., Luterbacher, J., and Macková, J.: Die Hungerjahre 1770-1772 in den böhmischen Ländern. Verlauf, meteorologische Ursachen und Auswirkungen, Österr. Z. Geschwiss., 12, 44-78, 2001. Brázdil, R., Pfister, C, Wanner, H., von Storch, H., and Luterbacher, J.: Historical climatology in Europe - the state of the art, Climatic Change, 70, 363^30, doi:10.1007/sl0584-005-5924-l, 2005. Brázdil, R., Zahradníček, P., Dobrovolný, P., Kotyza, O., and Valášek, H.: Historical and recent viticulture as a source of climatological knowledge in the Czech Republic, Geografie, 113, 351-371, 2008. Brázdil, R., Dobrovolný, P., Luterbacher, J., Moberg, A., Pfister, C, Wheeler, D., and Zorita, E.: European climate of the past 500 years: new challenges for historical climatology, Climatic Change, 101, 7^0, doi:10.1007/sl0584-009-9783-z, 2010. Brázdil, R., Bělínová, M., Dobrovolný, P., Mikšovský, J., Pišoft, P., Řezníčková, L., Štěpánek, P., Valášek, H., and Zahradníček, P.: History of Weather and Climate in the Czech Lands IX: Temperature and Precipitation Fluctuations in the Czech Lands During the Instrumental Period, Masaryk University, Brno, 2012a. Brázdil, R., Zahradníček, P., Pišoft, P., Štěpánek, P., Bělínová, M., and Dobrovolný, P.: Temperature and precipitation fluctuations in the Czech Lands during the period of instrumental measurements, Theor. Appl. Climatol., 110, 17-34, doi:10.1007/s00704-012-0604-3, 2012b. Brázdil, R., Kotyza, O., Dobrovolný, P., Řezníčková, L., and Valášek, H.: History of Weather and Climate in the Czech Lands X: Climate of the Sixteenth Century in the Czech Lands, Masaryk University, Brno, 2013. Bronisz, A., Biak, S., Bronisz, K, and Zásada, M.: Climate influence on radial increment of oak (Quercus SP.) in central Poland, Geochronometria, 39, 276-284, doi:10.2478/sl3386-012-0011-7, 2012. Büntgen, U, Frank, D. C, Nievergelt, D., and Esper, J.: Summer temperature variations in the European Alps, A.D. 755-2004, J. Climate, 19, 5605-5623, doi:10.1175/JCLI3917.1, 2006. Büntgen, U, Tegel, W., Nicolussi, K, McCormick, M., Frank, D. , Trouet, V, Kaplan, J., Herzig, F., Heussner, U., Wanner, H., Luterbacher, J., and Esper, J.: 2500 years of European climate variability and human susceptibility, Science, 331, 578-582, doi:10.1126/science.H97175, 2011. Büntgen, U, Tegel, W., Career, M., Krusic, P. J., Hayes, M., and Esper, J.: Commentary to Wetter et al. (2014): Limited tree-ring evidence for a 1540 European "Megadrought", Climatic Change, 131, 183-190, doi:10.1007/sl0584-015-1423-l, 2015. Burkhardt, T. and Hense, A.: On the reconstruction of temperature records from proxy data in Mid Europe, Arch. Meteor. Geophy. B., 35, 341-359, 1985. Chuine, I., Yiou, P., Viovy, N., Seguin, B., Daux, V, and Le Roy Ladurie, E.: Historical phenology: Grape ripening as a past climate indicator, Nature, 432, 289-290, doi:10.1038/432289a, 2004. Cook, B. I. and Wolkovich, E. M.: Climate change decouples drought from early wine grape harvests in France, Nature Climate Change, doi:10.1038/nclimate2960, online first, 2016. Cook, E. R., Briffa, K. R., and Jones, P. D.: Spatial regression methods in dendroclimatology: a review and comparison of two techniques, Int. J. Climatol., 14, 379-402, doi:10.1002/joc.3370140404, 1994. Daux, V, Garcia de Cortazar-Atauri, I., Yiou, P., Chuine, I., Gamier, E., Le Roy Ladurie, E., Mestre, O., and Tardaguila, J.: An open-access database of grape harvest dates for climate research: data description and quality assessment, Clim. Past, 8, 1403-1418, doi:10.5194/cp-8-1403-2012, 2012. Dobrovolný, P., Brázdil, R., Valášek, H., Kotyza, O., Macková, J., and Halíčková, M.: A standard paleoclimatological approach to temperature reconstruction in historical climatology: an example from the Czech Republic, A.D. 1718-2007, Int. J. Climatol., 29, 1478-1492, doi:10.1002/joc.l789, 2009. Dobrovolný, P., Moberg, A., Brázdil, R., Pfister, C, Glaser, R., Wilson, R., van Engelen, A., Limanówka, D., Kiss, A., Halíčková, M., Macková, J., Riemann, D., Luterbacher, J., and Böhm, R.: Monthly and seasonal temperature reconstructions for Central Europe derived from documentary evidence and instrumental records since AD 1500, Climatic Change, 101, 69-107, doi: 10.1007/s 10584-009-9724-x, 2010. Dobrovolný, P., Rybníček, M., Kolář, T, Brázdil, R., Trnka, M., and Büntgen, U: A tree-ring perspective on temporal changes in the frequency and intensity of hydroclimatic extremes in the territory of the Czech Republic since 761 AD, Clim. Past, 11, 1453-1466, doi:10.5194/cp-ll-1453-2015, 2015. Eddy, J. A.: The Maunder minimum, Science, 192, 1189-1202, 1976. Esper, J., Wilson, R. J. S., Frank, D. C, Moberg, A., Wanner, H., and Luterbacher, J.: Climates: past ranges and future changes, Quaternary Sei. Rev., 24, 2164-2166, doi: 10.1016/j .quascirev.2005.07.001, 2005. Etien, N, Daux, V, Masson-Delmotte, V, Stievenard, M., Bernard, V, Durost, S., Guillemin, M. T, Mestre, O., and Pierre, M.: A bi-proxy reconstruction of Fontainebleau (France) growing season temperature from A.D. 1596 to 2000, Clim. Past, 4, 91-106, doi:10.5194/cp-4-91-2008, 2008. Fischer, E. M., Seneviratne, S. I., Vidale, P. L., Lüthi, D., and Schär, C: Soil moisture-atmosphere interactions during the 2003 European summer heat wave, J. Climate, 20, 5081-5099, doi:10.1175/ICLI4288.1, 2007. Garcia de Cortázar-Atauri, I., Daux, V, Gamier, E., Yiou, P., Viovy, N, Seguin, B., Boursiquot, J. M., Parker, A., van Leeuwen, C, and Chuine, I.: An assessment of error sources when using grape harvest date for past climate reconstruction, Holocene, 20, 599-608, doi:10.1177/0959683609356585, 2010. Gamier, E., Daux, V, Yiou, P., and Garcia de Cortázar-Atauri, I.: Grapevine harvest dates in Besancon (France) between 1525 and 1847: Social outcomes or climatic evidence?, Climatic Change, 104, 703-727, doi:10.1007/sl0584-010-9810-0, 2010. Glaser, R. and Hagedorn, H.: The climate of lower Franconia since 1500, Theor. Appl. Climatol., 43, 101-104, 1991. Glaser, R. and Riemann, D.: A thousand year records of climate variation for Central Europe at a monthly resolution, J. Quaternary Sei., 24, 437^49, doi:10.1002/jqs.l302, 2009. www.clim-past.net/12/1421/2016/ Clim. Past, 12, 1421-1434, 2016 1432 M. Možný et al.: April-August temperatures in the Czech Lands Glaser, R., Brázdil, R., Pfister, C, Dobrovolný, P., Barriendos Valivé, M., Bokwa, A., Camuffo, D., Kotyza, O., Limanówka, D., Rácz, L., and Rodrigo, F. S.: Seasonal temperature and precipitation fluctuations in selected parts of Europe during the sixteenth century, Climatic Change, 43, 169-200, 1999. Grove, J. M.: Little Ice Age: Ancient and Modern, Routledge, London, 2004. Guiot, J., Nicault, A., Rathegeber, C, Edouard, J. L., Guibal, F., and Pichard, G: Last-millennium summer-temperature variations in western Europe based on proxy data, Holocene, 15, 489-500, doi: 10.1191/0959683605hl819rp, 2005. Hájková, L., Voženílek, V., Tolasz, R., Kohut, M., Možný, M., Nekovář, J., Novák, M., Reitschläger, J. D., Stříž, M., Vávra, A., and Vandrákova, A.: Atlas fenologických poměrů Česka (Atlas of Phenological Conditions in Czechia), Český hydrometeorologický ústav, Univerzita Palackého v Olomouci, Praha, Olomouc, 2012. Honzík, K.: Weinbau der Gräflich Sylva-Tarouca-Nostitz'schen Do-maine Czernosek, Prag, Karl Honzík, 1897. Hudson, I. L. and Keatley, M. R. (Eds.): Phenological research. Methods for environmental and climate change analysis, Springer, Dordrecht, 2010. Jones, G. V. and Hellman, E.: Site assessment and Oregon winegrowing regions, in: Oregon Viticulture, edited by: Hellman, E., Oregon State University Press, Corvallis, Oregon, 44-50, 51-69, 2003. Jones, G. V, Duff, A., Hall, A., and Myers, J.: Spatial analysis of climate in winegrape growing regions in the Western United States, Am. J. Enol. Viticult, 61, 313-326, 2010. Jones, G. V, White M. A., Cooper, O. R., and Storchmann, K.: Climate change and global wine quality, Climatic Change, 73, 319-343, doi:10.1007/sl0584-005-4704-2, 2005. Jones, G. V, Reid, R., and Vilks, A.: Climate, grapes, and wine: structure and suitability in a variable and changing climate, in: The Geography of Wine: Regions, Terroir and Techniques, edited by: Dougherty, P. H., Springer, New York, Berlin, 109-133, 2012. Jones, P. D., Briffa, K. R., Osborn, T. J., Lough, J. M., van Om-men, T. D., Vinther, B. M., Luterbacher, J., Wahl, E. R., Zwiers, F. W., Mann, M. E., Schmidt, G. A., Ammann, C. M., Buckley, B. M., Cobb, K. M., Esper, J., Goosse, H., Graham, N, Jansen, E., Kiefer, T., Kuli, C, Küttel, M., Mosley-Thompson, E., Overpeck, J. T., Riedwyl, N, Schulz, M., Tudhope, A. W., Villalba, R., Wanner, H., Wolff, E., and Xoplaki, E.: High-resolution paleoclimatology of the last millennium: a review of current status and future prospects, Holocene, 19, 3-49, doi: 10.1177/0959683608098952, 2009. Keller, M.: Managing grapevines to optimise fruit development in a challenging environment: a climate change primer for viticultur-ists, Aust. J. Grape Wine R., 16, 56-69, 2010. Kilián, J.: Vinařství mělnických měštanů za třicetileté války (Viticulture of the inhabitants of Mělník during the Thirty Years War), in: Víno a vinařství v českých zemích ve středověku a v raném novověku, edited by: Kilián, J., Regionální museum, Mělník, 123-131, 2009. Kilián, J.: Nástin vývoje mělnického vinařství v 19. a 20. století (Outline of the development of Mělník viticulture during the 19th and 20th centuries), in: Víno a vinařství v českých zemích po roce 1800, edited by: Babíková, P, Regionální museum, Mělník, 44-53, 2012. Kiss, A, Wilson, R., and Bariska, I.: An experimental 392-year documentary-based multi-proxy (vine and grain) reconstruction of May-July temperatures for Koszeg, West-Hungary, Int. J. Biometeorol., 55, 595-611, doi:10.1007/s00484-010-0367-4, 2011. Kraus, V, Foffová, Z., and Wurm, B.: Nová encyklopedie českého a moravského vína, (New Encyclopaedia of Bohemian and Moravian Wine, Praga Mystica), Praha, 2005. Krieger, M., Lohmann, G, and Laepple, T.: Seasonal climate impacts on the grape harvest date in Burgundy (France), Clim. Past, 7, 425^135, doi:10.5194/cp-7-425-2011, 2011. Landsteiner, E.: The crisis of wine production in late sixteenth-century Central Europe: Climatic causes and economic consequences, Climatic Change, 43, 323-334, doi:10.1007/978-94-015-9259-8_12, 1999. Le Roy Ladurie, E.: Canicules, fraicheurs, vendanges (France 15ěme-19ěme siěcles), C. R. Biol., 328, 213-218, 2005. Le Roy Ladurie, E. and Baulant, M.: Grape harvests from the fifteenth through the nineteenth centuries, J. Interdiscipl. Hist., 10, 839-849, 1980. Livezey, R. E.: Field intercomparison, in: Analysis of Climate Variability, edited by: von Storch, H. and Navarra, A., Applications of Statistical Techniques, Springer, Berlin, Heidelberg, 161-178, 2010. Luterbacher, J., Dietrich, D., Xoplaki, E., Grosjean, M., and Wanner, H.: European seasonal and annual temperature variability, trends and extremes since 1500, Science, 303, 1499-1503, doi: 10.1126/science. 1093877, 2004. Luterbacher, J., Werner, J., Smerdon, J., Fernandez-Dondo, L., Gonzalez-Rouco, J. F, Barriopedro, D., Ljungqvist, F. C, Bunt-gen, U, Zorita, E., Wagner, S., Esper, J., McCarroll, D., Toreti, A., Frank, D. C, Jungclaus, J. H., Barriendos, M., Bertolin, C, Bothe, O., Brázdil, R., Camuffo, D., Dobrovolný, P., Gagen, M., García-Bustamante, E., Ge, Q., Gomez-Navarro, J. J., Guiot, J., Hao, Z., Hegeři, G. C, Holmgren, K, Klimenko, V. V, Martin-Chivelet, J., Pfister, C, Roberts, N, Schindler, A., Schurer, A., Solomina, O., von Gunten, L., Wahl, E. R., Wanner, H., Wetter, O., Xoplaki, E., Yuan, N, Zanchettin, D., Zhang, H., and Zerefos, C: European summer temperatures since Roman times, Environ. Res. Lett, 11, 1-12, 2016. Mariáni, L., Parisi, S., Cola, F. O., Zoia, G, and Bonardi, L.: Tirani (1624-1930): a long time series of harvest dates for grapevine, Ital. J. Agrometeorol., 1, 7-16, 2009. Matthews, J. A. and Briffa, K. R.: The "Little Ice Age": re-evaluation of an evolving concept, Geogr. Ann., 87, 17-36, doi: 10.1111/j.0435-3676.2005.00242.X, 2005. Maurer, C, Koch, E., Hammerl, C, Hammerl, T., and Pokorný, E.: BACCHUS temperature reconstruction for the period 16th to 18th centuries from Viennese and Klosterneuburg grape harvest dates, J. Geophys. Res., 114, D22106, doi:10.1029/2009JD011730, 2009. McCarroll, D., Young, G, and Loader, N.: Measuring the skill of variance-scaled climate reconstructions and a test for the capture of extremes, Holocene, 25, 618-626, doi:10.1177/0959683614565956, 2015. Meier, N, Rutishauser, T., Pfister, C, Wanner, H., and Luterbacher, J.: Grape harvest dates as a proxy for Swiss April to August tem- Clim. Past, 12, 1421-1434, 2016 www.clim-past.net/12/1421 /2016/ M. Možný et al.: April-August temperatures in the Czech Lands 1433 perature reconstructions back to AD 1480, Geophys. Res. Lett., 34, L20705, doi:10.1029/2007GL031381, 2007. Menzel, A.: A 500 year pheno-climatological view on the 2003 heatwave in Europe assessed by grape harvest dates, Meteorol. Z., 14, 75-77, doi:10.1127/0941-2948/2005/0014-0075, 2005. Molitor, D., Udelhoven, T., Ney, S., Hoffmann, L., and Poster, L.: Historical vintage descriptions from Luxembourg - an indicator for the climatic conditions in the past?, Vitis, 55, 23-30, doi:10.5073/vitis.2016.55.23-30, 2016. Montgomery, D. C, Peck, E. A., and Vinning, G. G: Introduction to Linear Regression Analysis, John Wiley & Sons, New York, 2012. Možný, M., Brázdil, R., Dobrovolný, P., and Trnka, M.: Cereal harvest dates in the Czech Republic between 1501 and 2008 as a proxy for March-June temperature reconstruction, Climatic Change, 110, 801-821, doi:10.1007/sl0584-011-0075-z, 2012. Možný, M., Brázdil, R., Dobrovolný, P., Trnka, M., Potopová, V., Hlavinka, P., Bartošová, L., Zahradníček, P., Štěpánek, P., and Žalud, Z.: Drought reconstruction based on grape harvest dates for the Czech Lands in AD 1499-2012, Clim. Res., in review, 2016. Možný, M., Hájková, L., Vráblík, T, Novák, J., Reitschläger, J., and Bareš, D.: Počasí a výnosy zemědělských plodin v České republice (Weather and Yields of Agricultural Crops in the Czech Republic), Sborník prací ČHMÚ, Praha, 2016. Pfister, C: Die Fluktuationen der Weinmosterträge im schweizerischen Weinland vom 16. bis ins frühe 19. Jahrhundert: klimatische Ursachen und sozioökonomische Bedeutung, Schweizerische Z. Gesch., 31, 445^91, 1981. Pfister, C: Klimageschichte der Schweiz 1525-1860, Das Klima der Schweiz von 1525-1860 und seine Bedeutung in der Geschichte von Bevölkerung und Landwirtschaft, Bd. 1, Paul Haupt, Bern, Stuttgart, 1984. Pfister, C: Wetternachhersage: 500 Jahre Klimavariationen und Naturkatastrophen, Paul Haupt, Bern, 1999. Pfister, C. and Brázdil, R.: Climatic variability in sixteenth-century Europe and its social dimension: a synthesis, Climatic Change, 43, 5-53, doi:10.1023/A:1005585931899, 1999. Pfister, C. and Brázdil, R.: Social vulnerability to climate in the "Little Ice Age": an example from Central Europe in the early 1770s, Clim. Past, 2,115-129, doi:10.5194/cp-2-115-2006, 2006. Pfister, C, Wetter, O., Brázdil, R., Dobrovolný, P., Glaser, R., Luter-bacher, J., Seneviratne, S. I., Zorita, E., Alcoforado, M. J., Bar-riendos, M., Bieber, U., Burmeister, K. H., Camenisch, C, Con-tino, A., Grünewald, U., Herget, J., Himmelsbach, I., Labbé, T, Limanówka, D., Litzenburger, L., Kiss, A., Kotyza, O., Nordli, 0., Přibyl, K, Retsö, D., Riemann, D., Rohr, C, Siegfried, W., Spring, J.-L., Söderberg, J., Wagner, S., and Werner, J. P.: Tree-rings and people - different views on the 1540 Megadrought. Reply to Büntgen et al. 2015, Climatic Change, 131, 191-198, doi:10.1007/sl0584-015-1429-8, 2015. Pejml, K: Příspěvek ke kolísání klimatu v severočeské vinařské a chmelařské oblasti od r. 1500-1900 (A Contribution to Climate Fluctuations in the North Bohemian Wine- and Hop-growing Region in 1500-1900), Český hydrometeorologický ústav, Praha, 1966. Pejml, K: Příspěvek ke znalosti kolísání klimatu v Čechách v 16. až 18. stol (A contribution to knowledge of climate fluctuations in Bohemia from the 16th to the 18th centuries), Meteorol. Zprav., 27, 90-95, 1974. Schams, F.: Vollständige Beschreibung sämtlicher berühmten Weingebirge in Oesterreich, Mähren und Böhmen in statistisch-topographisch-natur-historischer und ökonomischer Hinsicht, 1835. Schär, C. and Jendritzky, G: Hot news from summer 2003, Nature, 432, 559-560, doi:10.1038/432559a, 2004. Smetana, J.: Nejstarší kronikárske záznamy litoměřických radních písařů (The earliest chronicle records of scribes to the Litoměřice councillors), Litoměřicko, 14, 119-142, 1978. Sotolář, R.: Historické odrůdy révy na našem území (The historical varieties of grapevine in the Czech Lands), in: Víno a vinařství v českých zemích ve středověku a v raném novověku, edited by: Kilián, J., Regionální muzeum, Mělník, 16-23, 2009. Stocker, T. F., Qin, D., Plattner, G-K, Tignor, M. M. B., Allen, S. K, Böschung, J., Nauels, A., Xia, Y, Bex, V, and Midgley, P. M. (Eds.): Climate Change 2013: The Physical Science Basis, Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, UK, 2013. Stott, P. A., Stone, D. A., and Allen, M. R.: Human contribution to the European heatwave of 2003, Nature, 432, 610-614, doi:10.1038/nature03089, 2004. Strömmer, E.: Klima-Geschichte. Methoden der Rekonstruktion und historische Perspektive. Ostösterreich 1700 bis 1830, Franz Deuticke, Wien, 2003. Svitáková, Z., Kott, L, and Nekovář, J.: Fenologická data za posledních 150 let (Phenological data for the recent 150-year period), in: Bioclimatology present and future, Krtiny, 12-14 September 2005. Simáček, J.: Vinařství v království českém: spis jednající o počátcích, rozvoji a nynějším stavu vinařství/ Der Weinbau im Königreiche Böhmen: Schilderung der Anfänge, der Entwicklung und des gegenwärtigen Standes des Weinbaues, A. Reinwart, Praha, 1891. Tolasz, R., Míková, T, Valeriánová, A., and Voženílek, V. (Eds.): Climate Atlas of Czechia, Czech Hydrometeorological Institute and Palacký University, Prague, Olomouc, 2007. Trachsel, M., Kameník, C, Grosjean, M., McCarroll, D., Moberg, A., Brázdil, R., Büntgen, U., Dobrovolný, P., Esper, J., Frank, D. C, Friedrich, M., Glaser, R., Larocque-Tobler, L, Nicolussi, K., and Riemann, D.: Multi-archive summer temperature reconstruction for the European Alps, AD 1053-1996, Quaternary Sei. Rev., 46, 66-79, doi:10.1016/j.quascirev.2012.04.021, 2012. van Engelen, A. F. V, Buisman, J., and Ijnsen, F: A millennium of weather, winds and water in the Low Countries, in: History and Climate: Memories of the Future?, edited by: Jones, P. D., Ogilvie, A. E. J, Davies, T. D., and Briffa, K. R., Kluwer Academic/Plenum Publishers, New York, Boston, Dordrecht, London, Moscow, 101-124, 2001. Wetter, O. and Pfister, C: Spring-summer temperatures reconstructed for northern Switzerland and southwestern Germany from winter rye harvest dates, 1454-1970, Clim. Past, 7, 1307-1326, doi:10.5194/cp-7-1307-2011, 2011. Wetter, O. and Pfister, C: An underestimated record breaking event - why summer 1540 was likely warmer than 2003, Clim. Past, 9, 41-56, doi:10.5194/cp-9-41-2013, 2013. www.clim-past.net/12/1421/2016/ Clim. Past, 12, 1421-1434, 2016 1434 M. Možný et al.: April-August temperatures in the Czech Lands Wetter, O., Pfister, C, Werner, J. P., Zorita, E., Wagner, S., Senevi-ratne, S. I., Herget, J., Grünewald, U., Luterbacher, J., Alco-forado, M. J., Barriendos, M., Bieber, U., Brázdil, R., Burmeis-ter, K. H., Camenisch, C, Contino, A., Dobrovolný, P., Glaser, R., Himmelsbach, I., Kiss, A., Kotyza, O., Labbé, T., Li-manówka, D., Litzenburger, L., Nordli, 0., Přibyl, K., Retsö, D., Riemann, D., Rohr, C, Siegfried, W, Söderberg, J., and Spring, J.-L.: The year-long unprecedented European heat and drought of 1540 - a worst case, Climatic Change, 125, 349-363, doi:10.1007/sl0584-014-1184-2, 2014. Wilson, P. H.: The Thirty Years War: Europe's Tragedy, Harvard University Press, Harvard, 2009. Wilson, R. and Topham, J.: Violins and climate, Theor. Appl. Chmatol., 77, 9-24, doi:10.1007/s00704-003-0025-4, 2004. Clim. Past, 12, 1421-1434, 2016 www.clim-past.net/12/1421 /2016/