Variations in the Peczely Macro-Synoptic Types (1881–2020) with Attention to Weather Extremes in the Pannonian Basin

: Daily Peczely circulation types are investigated over 140 years (1881–2020). After presenting monthly mean frequencies and durations of the 13 circulation types, two further questions are investigated: (i) How do the circulation types inﬂuence local weather extremes?; (ii) Are there signiﬁcant trends in the frequency of the original and the grouped circulation types in the recent monotonically warming 50 year period (1971–2020)? The answers are as follows: (i) Four local weather extremes were investigated in nine grid-points of the Pannonian Basin and analyzed in the central months of the seasons. It was established that high precipitation and wind maxima occur in almost all circulation types and months, whereas for both high temperature maxima and low temperature minima, there are six circulation types, where no extremity occurred in one, two, or three investigated months. (ii) In the last 50 years, 37% of the linear seasonal frequency trends have been signiﬁcant. However, these trends are rarely signiﬁcant in the shorter monotonously warming (1911–1940) and cooling (1941–1970) 30-year periods. Therefore, the signiﬁcant trends of the last 50 years are unlikely to be the direct consequences of the monotonous hemispherical warming. Since these hemispherical temperature trends are most likely caused by different sets of physical reasons, the reality of the presented circulation frequency trends needs to be validated by climate models.


Introduction
Weather extremes are connected with large-scale circulation patterns, including the mesoscale objects embedded in them. The present paper deals with two main aspects of atmospheric circulation, as follows: • the classification of weather maps into pre-defined discrete types of circulation; • the impact of global warming on large-scale circulation.
Macro-synoptic (or macro-circulation) classification is based on spatial fields of sea level pressure or mid-tropospheric geopotential (e.g., [1][2][3] in case of the Peczely types). Since its first appearance in the literature, probably connected with the names of van Bebber and Köppen [4], macro-synoptic types often meet application in diagnostic and sometimes in prognostic tasks of descriptive or applied climatology, weather forecasting (e.g., [5][6][7]), and in air-pollution meteorology [8].
The description and catalogue of Hungary's macro-synoptic fields appeared in the work of Péczely [1]. This catalogue contains a daily classification of the sea-level pressure fields from 1877 until 1956. The next publication of the updated code series was issued 26 years later [2].
Peczely's classification has been widely used for diagnostic purposes (e.g., [9,10]) to understand connections between the circulation and local weather, sometimes including its variability, too. In air-pollution meteorology, where relatively short series of measurements Both approaches make it possible to provide analyses of temporal changes [56,57]. In [56], long-term variations are established at multi-decadal time scales, but they are also different in the different seasons. At the same time, [57] finds in several global climate model experiments that no unequivocal changes can be established in Central Europe either in summer or in winter.
This result based on global climate models contradicts several other investigations that explain a more persistent flow over Europe, explaining this tendency with weakening of the temperature gradient, i.e., weakening of the tropospheric zonal flow, as a result of human-induced climate change [58][59][60][61].
The present paper is further structured as follows: Section 2 comprehends the analyzed circulation and local weather data including those that characterize the hemispherical warming during the last 50 years. Section 3 presents the statistical results in three sub-Sections starting with the primary results of synoptic climatology. The second part of the results represent the distribution of selected daily extreme values among the circulation types. Finally, in the third part of the results linear trends of the various circulation types are analyzed for the 1971-2020 period of hemispherical warming, compared to three other 30-year periods of the total 140 years of existence of the Peczely circulation types. The results are discussed in Section 4, followed by a short conclusion in Section 5 and the References.

Materials and Methods
The 13 macro-circulation types defined by Peczely [1,2] are listed in Table 1, in the original grouping, based on the typical wind direction. The coding takes place subjectively, but the 140-year archive of diurnal types (1881-2020) is the product of two experts, only. Until 1983, the author of the classification performed the coding, continued by his colleague, co-author of the present paper, C. Karossy. Table 1. The 13 circulation types of the Peczely (1957) macro-synoptic classification [1], originally defined for Hungary (Hu). The original table by Peczely was edited in German, to characterize the 13 circulation types. The brief descriptions of each type were translated to English by one of the authors (JM).

Meridional Types Zonal and Central Types
Types Connected with Western Current: Types Connected with Northern Current: zC-zonal flow, slightly cyclonic influence mCc-Hu is in the rear of a West-European cyclone Aw-anticyclone extending from the west AB-anticyclone over the British Isles As-anticyclone to the south from Hungary CMc-Hu is in the rear of a Mediterranean cyclone Types Connected with Eastern Current: An-anticyclone to the north from Hungary Types Connected with Southern Current: AF-anticyclone over Fenno-Scandinavia mCw-Hu is in the fore of a West-European cyclone Types of Pressure Centers: Ae-anticyclone to the east from Hungary A-anticyclone center over Hungary CMw-Hu is in the fore of a Mediterranean cyclone C-cyclone center over Hungary All methods and results involve this classification. Figure 1 presents examples of each circulation type, grouped according to the typical direction of the driving current, as presented above in Table 1.
The majority of the calculations apply to the Excel software package, except for the significance of the linear trends in Section 3.3, which applied Statistical Package for the Social Sciences (SPSS ver. 20).
In the first set of results (see in Section 3.1) univariate monthly statistics of average frequency and average duration is calculated for the most recent 1991-2020 period. These calculations follow the definition of relative frequency and duration. The majority of the calculations apply to the Excel software package, except for the significance of the linear trends in Section 3.3, which applied Statistical Package for the Social Sciences (SPSS ver. 20).
In the first set of results (see in Section 3.1) univariate monthly statistics of average frequency and average duration is calculated for the most recent 1991-2020 period. These calculations follow the definition of relative frequency and duration.
The results in Section 3.2 will present the frequency distribution of selected weather extremes among the 13 circulation types. Grid-points A-I, selected to characterize extreme values in the Pannonian Basin, are presented in Figure 2. The coordinates of the central grid-point are 46.9 °N 20.4 °E with 1.5 deg. longitude or 1.0 deg. latitude difference between the neighboring grid-points. These differences correspond to ca. 100 × 100 km spatial resolution. Five grid-points fall into Hungary, two of them into Romania and one each into Serbia and Croatia. The results in Section 3.2 will present the frequency distribution of selected weather extremes among the 13 circulation types. Grid-points A-I, selected to characterize extreme values in the Pannonian Basin, are presented in Figure 2. The coordinates of the central grid-point are 46.9 • N 20.4 • E with 1.5 deg. longitude or 1.0 deg. latitude difference between the neighboring grid-points. These differences correspond to ca. 100 × 100 km spatial resolution. Five grid-points fall into Hungary, two of them into Romania and one each into Serbia and Croatia.
The nine grid-point data are taken from the CarpatClim daily database [62,63]. Carpat-Clim is a freely available, high resolution gridded database. This database contains homogenized, interpolated, daily-scale meteorological parameters between 1961 and 2010. Climatological grids cover the area between latitudes 44 • N and 50 • N and longitudes 17 • E and 27 • E. Spatial resolution of the database is 0.1 • × 0.1 • geographical latitudes and longitudes, i.e., ca. 11 × 7 km in the region of our 9 grid-points. The CarpatClim database contains meteorological data for various parts of seven countries including Croatia, Hungary, Poland, Romania, Serbia, Slovakia and Ukraine. The nine grid-point data are taken from the CarpatClim daily database [62,63]. Car-patClim is a freely available, high resolution gridded database. This database contains homogenized, interpolated, daily-scale meteorological parameters between 1961 and 2010. Climatological grids cover the area between latitudes 44 °N and 50 °N and longitudes 17 °E and 27 °E. Spatial resolution of the database is 0.1° × 0.1° geographical latitudes and longitudes, i.e., ca. 11 × 7 km in the region of our 9 grid-points. The CarpatClim database contains meteorological data for various parts of seven countries including Croatia, Hungary, Poland, Romania, Serbia, Slovakia and Ukraine.
The gridded database has been created from the observed data at the meteorological stations by homogenization (MASH) and interpolation procedure (MISH) [64].
Technical details (e.g., types of observing devices and number of stations) are detailed in [65], grouped according to the participating countries and the meteorological variables.
The central months of each season, i.e. January, April, July and October are presented, only. Note that in Section 3.3 the seasons mean the conventional three-monthly periods, i.e., December-February for winter, March-May for spring, June-August for summer and September-November for autumn.
By selection of the grid-points, we intended to characterize mostly the plain areas of the Pannonian Basin. Therefore, seven grid-points, from the selected nine, correspond to 73-136 m altitude above sea level. Only the points A and B have 549 and 233 m altitude, respectively. Since even point A does not really represent a mountainous area, we did not specifically analyze the differences among the nine grid-points.
Extremity is defined by the 5% threshold values comprehended in Table 2. Higher than the upper thresholds are investigated for daily precipitation, maximum temperature and wind speed values, whereas for the minimum temperature, the values cooler than the lower 5% thresholds are analyzed. The gridded database has been created from the observed data at the meteorological stations by homogenization (MASH) and interpolation procedure (MISH) [64].
Technical details (e.g., types of observing devices and number of stations) are detailed in [65], grouped according to the participating countries and the meteorological variables.
The central months of each season, i.e., January, April, July and October are presented, only. Note that in Section 3.3 the seasons mean the conventional three-monthly periods, i.e., December-February for winter, March-May for spring, June-August for summer and September-November for autumn.
By selection of the grid-points, we intended to characterize mostly the plain areas of the Pannonian Basin. Therefore, seven grid-points, from the selected nine, correspond to 73-136 m altitude above sea level. Only the points A and B have 549 and 233 m altitude, respectively. Since even point A does not really represent a mountainous area, we did not specifically analyze the differences among the nine grid-points.
Extremity is defined by the 5% threshold values comprehended in Table 2. Higher than the upper thresholds are investigated for daily precipitation, maximum temperature and wind speed values, whereas for the minimum temperature, the values cooler than the lower 5% thresholds are analyzed. The applied 5% threshold values are presented in Table 2 for the nine grid-points displayed in Figure 2. In addition to it, the absolute maxima or minima for the same nine grid-points are included in Table 3. The values of both Tables are derived independently from the circulation codes, and they can be considered important additions to the climate of the Pannonian Basin. The results presented in Section 3.2 indicate the distributions of the 5% extremes among the 13 circulation types calculated as the average of the nine grid-points in the central months of each season.
The last part of the investigations, detailed in Section 3.3, is the linear trend analysis of the original circulation types in the preliminarily selected 50 years (1971-2020) of monotonic increase of annual mean hemispherical averages of near-surface air temperature, together with three 30-year periods, characterized by neutral (1881-1910), warming , and cooling  periods, as presented in Figure 3. The most important Northern Hemispherical values, characterizing the above 30-and 50-year periods are included in Table 4. The regression coefficients of this Table are expressed in • C/10 years.  -1910, 1911-1940, 1941-1970, 1971-2020. As concerns the longest monotonous period of 50 years, its warming tendency is "very likely" influenced by the increasing concentration of greenhouse gases [16]. Hence, this 50-year period can be treated as a natural experiment, simulating regional features of global warming. The external factors of the shorter warming, cooling, and neutral periods are most likely influenced by different global reasons.  1881-1910, 1911-1940, 1941-1970, 1971-2020. The source of information is https://crudata.uea.ac.uk/cru/data/temperature/Had-CRUT5.0Analysis.pdf [66]. The link for the licensing conditions is seen in the description of [66], in the list of References.
Finally, the significance of the linear trends for the local circulation frequency in the 50 years were established if both the correlation coefficients and the linear regression coefficient were significant according to the Z-test (Table A9 in [67]) and t-test [68] applied by SPSS, respectively. Since these two tests were statistically significant, for the same circulation types and seasons in the 50-year period, further below significance according to the correlation coefficients were considered in case of the 30-year periods, too. The 95% significance threshold for the correlation coefficients is 0.27 for the 50-year samples and 0.35 for the 30-year samples. Finally, the significance of the linear trends for the local circulation frequency in the 50 years were established if both the correlation coefficients and the linear regression coefficient were significant according to the Z-test (Table A9 in [67]) and t-test [68] applied by SPSS, respectively. Since these two tests were statistically significant, for the same circulation types and seasons in the 50-year period, further below significance according to the correlation coefficients were considered in case of the 30-year periods, too. The 95% significance threshold for the correlation coefficients is 0.27 for the 50-year samples and 0.35 for the 30-year samples.

Monthly Mean Frequency and Duration of the Circulation Types
The frequency distribution among the circulation types and the average duration are often applied characteristics of synoptic climatology. Table 5 and Figure 4 present these values for the 13 types in monthly resolutions. In the case of the frequency distribution, the annual mean values are also presented.   The average duration of the circulation types ( Figure 4) is grouped according to Peczely [1,2], considering the main direction of near-surface wind, as displayed above in Table 1.
The smallest average duration can be 1.0 days, which occurs when all appearances of the given type remain valid for one day only. This situation is valid in several months for the CMc type with 1-2 months of occurrence for CMw, zC, and C types. The anticyclones mostly represent the longest duration of the given wind direction group with a characteristic duration of 1.5-2.5 days on average. The only exception to this is the type zC in September, which had over 2.5 days' average duration.

Distribution of Weather Extremes among the Circulation Types
Distribution of selected daily extreme values have been analyzed in nine grid-points presented above in Figure 2. High precipitation, high diurnal temperature maxima, low diurnal minima, and high diurnal wind speed values are considered. Threshold values of the 5% frequency and the absolute maxima are presented above in Tables 2 and 3, respectively.
In this sub-Section, the frequency distributions that occurred in all nine grid-points are presented in Figure 5. the 5% frequency and the absolute maxima are presented above in Tables 2 and 3, respectively.
In this sub-Section, the frequency distributions that occurred in all nine grid-points are presented in Figure 5   The high diurnal precipitation sums occur most frequently in mCc, CMw, and Aw, whereas rather rarely in the zC, As, and A macro-circulation types. The high frequency of much rain in the Aw situation is not surprising, as in the praxis of storm warning for the Lake Balaton this is a dangerous situation if a hidden cold front is contained by the so-called "nose of the anticyclone" spreading ahead from the main body of the high pressure to the west from the Pannonian Basin [69]. Otherwise, one can establish that almost all macro-circulation types may be accompanied by heavy rains in all central months of the seasons.
For the daily wind speed (averaged for the whole day from 10 min periods) similar behavior can be established: the diurnal wind speed maxima occur in almost each original circulation type in each investigated month. The highest proportion of extremely strong winds occur in the Aw types of circulation. In July this frequency reached 40%, which means that this proportion of the extremes occur in days with this type of circulation.
As concerns the temperature extremes, there are several macro-circulation types and seasons when no extreme values occur in any of the nine grid-points. For both the high maximum temperatures and for the low minimum temperatures, there are six types for both temperature extremes when no extreme occurs in one, two, or three investigated months. The highest maximum temperature extremes occur in the days with Ae types, whereas for the lowest temperature minima are different in the different investigated months. In July and in October over 30% of these extremes belong to Aw and Ae types, respectively.

Trends in Monotonous Hemispherical Warming and Cooling Periods
From the period 1971-2020, only those circulation types are involved in the following Figures and Tables that are significant according to both the correlation and the regression coefficients. As can be seen in Table 6, these two criteria coincide for both statistical parameters. In other words, both parameters found identical seasons and circulation types to be significant. More exactly, the Table indicates identical numbers of significance of both coefficients in each season. The pair-wise equivalence of significance could be established from the detailed matrices of significance according to the two statistical indicators. As can be seen from Table 6, 37% of the original types are significant. This means 19 significant trends from the 13 types × four seasons = 52 possible investigated coefficients. The 19 coefficients are the sum of the corresponding lines in the table, indicated in the various seasons.
The 37% of the significant linear trends is much higher than the random occurrence (5%). Hence, in Figures 6-9, the graphs of significant linear trends are presented with the best fitted linear regression lines. The four figures contain the seasonal trends of frequency from winter to autumn. The numerical values of these linear trends are presented later in this Section. There are only two types that have significant linear trends in all four seasons. These types are the mCc (increasing) and CMc (decreasing), which exhibit significant trends in all seasons. The zC is represented by a significant trend in three seasons. Three types, mCw, Ae, and CMw, exhibit significant trends in two seasons, whereas An is characterized by a significant trend only in one season. As mentioned above, all the trends from zC to An are negative (decreasing).
It is worth mentioning that for a given type, the signs of the significant trends are the same in all seasons. From the total of 19 significant coefficients in the four seasons, only mCc exhibits an increasing linear trend, whereas the other seven types are characterized by decreasing trends parallel to the hemispherical warming in 1971-2020.        Concerning the above statistically established increasing and decreasing trends in the recent 50 years, it would be ideal if one could establish whether or not they are direct consequences of hemispherical warming. This question can finally be answered by climate modelling, which is not available for the authors of the present paper. Therefore, in the last two Tables of this study the circulation frequency trends are also established for the other three 30-year periods specified in Table 4 of Section 2. The three periods, taken from the latter Table, are the neutral 1881-2010, the warming 1911-1940, and the cooling 1941-1970 years, based on hemispherical mean temperature, in the hope that the warming 30 years produce trends similar to the 50 years of warming, and the cooling 30 years yields opposite results, with no significant circulation trends in the neutral 30 years. Table 7 presents the proportions of significant linear trends of circulation types in these 30-year periods for comparison to the trends in the recent 1971-2020 period. One can establish that the proportions of significant linear trends are much smaller: 10 and 13% for the shorter warming and cooling periods, respectively. In addition to this, the neutral period exhibits similarly low frequency, 12%. The significant regression coefficients are presented for the 13 types in Table 8. For the 30-year periods, only those circulation types and seasons are included that are significant in the 50-year period.
From the latter 19 regression coefficients, found significant for the 13 circulation types in 1971-2020, there are only seven that are significant in any other 30-year period. From these coefficients, five occurred in warming or cooling periods, signs of which are appropriate in four cases, only (Table 8).  Figures 6-9 display the seasonal linear trends found equally significant according to both tests concerning correlation (Z-test) and regression (t-test) described in Section 2. Only eight macro-synoptic types exhibit significant frequency trends in any of the seasons (see below), and in five types there are no significant trends at all. They are the AB, As, AF, A, and C types.
There are only two types that have significant linear trends in all four seasons. These types are the mCc (increasing) and CMc (decreasing), which exhibit significant trends in all seasons. The zC is represented by a significant trend in three seasons. Three types, mCw, Ae, and CMw, exhibit significant trends in two seasons, whereas An is characterized by a significant trend only in one season. As mentioned above, all the trends from zC to An are negative (decreasing).
It is worth mentioning that for a given type, the signs of the significant trends are the same in all seasons. From the total of 19 significant coefficients in the four seasons, only mCc exhibits an increasing linear trend, whereas the other seven types are characterized by decreasing trends parallel to the hemispherical warming in 1971-2020.
Concerning the above statistically established increasing and decreasing trends in the recent 50 years, it would be ideal if one could establish whether or not they are direct consequences of hemispherical warming. This question can finally be answered by climate modelling, which is not available for the authors of the present paper. Therefore, in the last two Tables of this study the circulation frequency trends are also established for the other three 30-year periods specified in Table 4 of Section 2. The three periods, taken from the latter Table, are the neutral 1881-2010, the warming 1911-1940, and the cooling 1941-1970 years, based on hemispherical mean temperature, in the hope that the warming 30 years produce trends similar to the 50 years of warming, and the cooling 30 years yields opposite results, with no significant circulation trends in the neutral 30 years. Table 7 presents the proportions of significant linear trends of circulation types in these 30-year periods for comparison to the trends in the recent 1971-2020 period. One can establish that the proportions of significant linear trends are much smaller: 10 and 13% for the shorter warming and cooling periods, respectively. In addition to this, the neutral period exhibits similarly low frequency, 12%. The significant regression coefficients are presented for the 13 types in Table 8. For the 30-year periods, only those circulation types and seasons are included that are significant in the 50-year period. From the latter 19 regression coefficients, found significant for the 13 circulation types in 1971-2020, there are only seven that are significant in any other 30-year period. From these coefficients, five occurred in warming or cooling periods, signs of which are appropriate in four cases, only (Table 8).
In summary, we can conclude that the relatively high 37% proportion of significant tendencies in the circulation frequency are much higher than their random occurrence. However, the question of whether or not these trends are universal consequences of the hemispherical warming is discussed below.

Discussion
As was mentioned in the Introduction, both empirical and objective methods are applied to establish circulation classifications and their daily coding. In the case of empirical coding of the Peczely types, it is advantageous that the coding is performed by two experts, only, who even worked together as a professor and his assistant for some years [3].
Monthly mean relative frequencies and durations are presented for the most recent 30-year period (1991-2020). Although these values are often published [8,9,13,15] based on various time sequences, the values presented in this paper are the most recent ones for a standard 30-year period.
The interesting difference between the distribution of precipitation and wind extremes, on one side, and the temperature maxima and minima, on the other side are that the extremes of the first group may occur practically in all circulation types, whereas those on the second side occur in all seasons only in seven types. In other words, there are six types from the 13, where in one, two or three season-central months where no 5% extremes occur.
The most likely explanation of this behavior is that the lifetime of the mesoscale objects leading to heavy rain and strong wind is much shorter than the 24 h for which the macro-synoptic codes are defined [5,69,70]. This means that these mesoscale objects may develop practically under any macro-synoptic code. On the other hand, formation of the temperature extremes needs a longer time at the temperate latitudes, hence in case of several macro-synoptic situations the conditions are not favorable for the temperature extremes.
The applied 50-year (1971-2020) period exhibits a steep hemispherical warming tendency (+0.28 • C/10 yr) also characterized by a very high (0.954) correlation coefficient. The previous warming period  also exhibits steep (+0.20 • C/10 yr) regressions and high (0.843) correlation coefficients. Both periods are recommended to use as first approximations of regional features for global climate change, although the set of external forcing factors are probably different in the two warming periods [16].
Concerning the frequency trends analyzed in the recent 50 years (1971-2020) there are three questions to discuss: (i) How are the significant trends related to the former analyses of circulation changes? (ii) Can we validate these trends by performing a similar investigation for shorter time periods? (iii) How could one involve the climate models in such a validation process?
In relation to (i), one can establish that 37% of the 13 macro-types and seasons provided significant trends at the 95% significance level, which is much more than their random occurrence, i.e., 5%. However, the former expectations for the increase in persistent circulation [58][59][60][61] parallel to global warming are not supported by our results, which rather support the conclusions of [56,57], stating seasonal differences in the statistical investigations and a lack of unequivocal changes in the circulation projected by global climate models.
Concerning (ii), we should establish that the circulation frequency trends of the shorter periods occur much more rarely than in the 50-year period. This means that by using only this statistical technique of linear trend analysis, we cannot decide whether or not they are universal consequences of global warming. Even if the proportion of significant trends is much smaller, the external forcing factors can also be sources of these big differences between the results of the 50-year period and the 30-year periods.
For question (iii), objective classification of the daily weather maps [52][53][54][55] would be very useful, especially in terms of those meteorological fields that are available in the climate models. Here, both the General Circulation Models (GCM) and the Regional Climate Models (RCM) can be considered. Validation of climate models in terms of circulation and the physical reasoning of the derived frequency trends could be performed based on this objective classification.
The significant trends in frequency of the macro-synoptic types (Figures 6-9) can be combined with the occurrence of weather extremes in the various macro-types ( Figure 5) to assess how the frequency of the given weather extremes are expected to change parallel to global warming. For example, the observed increase in type mCc in all seasons during the recent 50 years may increase the frequency of high precipitation events in the future. However, as mentioned above in this Section, one cannot be sure that the observed circulation trends are universal consequences of global warming.

Conclusions
Finally, let us briefly list the conclusions of the presented investigations: As the result of subjective coding by only two experts, one of the longest diurnal series of macro-circulation series exist for the Pannonian Basin.
The monthly frequency and duration statistics are available for each of the 13 macrosynoptic types for the last 30 years' standard period .
Nine grid-points are selected from the CarpatClim database to obtain a ca. 100 × 100 km resolution network to identify four diurnal weather extremes in the central months of each season and to establish their occurrence parallel to the macro-circulation types. The 5% threshold values and the absolute extremity are presented for each grid-point.
The distribution of 5% extreme values among the macro-circulation types is presented as the average of the nine grid-points. An unexpected difference was found between distribution of precipitation and wind extremes, on one hand, and the temperature maxima and minima, on the other hand. More specifically, the extremes of the first group may occur practically in all circulation types, whereas those of the second group are concentrated in seven circulation types with no appearance in one-three season-center months in the other six types. For the likely reason for this difference, see the Discussion.
One 50-year period (1971-2020), characterized by monotonous hemispherical warming, has been defined, as well as three 30-year sequences of warming , cooling , and neutral (1881-1910) periods. They are tackled as natural experiments for expected further global warming. This approach, however, has not been supported by our trend estimations, as is discussed in Section 4. Therefore, the established trends detected in the circulation types can be considered facts of the given 50-year period but not necessarily general consequences of global warming, as discussed in Section 4.
Author Contributions: J.M. carried out the conceptualization and the majority of the statistical computations. C.K. provided the daily codes of the macro-types, selected and created Figures 1 and 2, and collected most of the Hungarian references. L.L. contributed with research ideas, provided or improved parts of figures, and significantly contributed to the collection and evaluation of the international references. All authors have read and agreed to the published version of the manuscript.