Evaluation of the Effect of Average Annual Temperatures in Slovakia between 1971 and 2020 on Stresses in Rigid Pavements

: The scientiﬁc community in Central Europe often discusses the extent to which temperature changes over the last two decades have contributed to changing the stresses induced in structures. In the ﬁeld of road and environmental engineering, this question is especially pertinent for pavements. The pavement structure must ﬁrst be deﬁned methodologically by identifying and deﬁning the types of parameters that change with time. Additionally, it is important to identify the areas of Central Europe that are most affected by climate change. The most important parameters must be described statistically for these areas. Slovakia is one of the countries that may be able to contribute to the solution of this issue due to its location in the middle of Europe. This paper provides a statistical analysis for the period from 1971 to 2020 in Slovakia. A concrete pavement, which is the most commonly used type of pavement, must be used as an example to numerical assess the situation. The conclusions and discussion in this scientiﬁc ﬁeld are directed towards the evaluation of the measurement results in the context of the designed pavement composition and the calculations using the different methods speciﬁed in the standards.


Introduction
Increased average annual temperatures are linked to global warming, a worldwide problem. Central Europe (CE) is a very specific area. It is characterized throughout by a wide range of temperature zones, as well as different traffic loads and infrastructure. There has been a significant urbanization process in these countries that influences temperature changes. This process involves covering the Earth's surface with different types of roads. As a result, grassland areas are being reduced. This effect is called the urban heat island effect (UHI) [1,2]. It is related to the increasing surface temperature of concrete pavements and their accumulation capacity. A rigid pavement and its surroundings form a single, interdependent system. Climate change (temperature, precipitation, groundwater level, freezing and thawing cycles) influences the structural properties and interactions of this system [3][4][5][6].
The most significant factor is temperature change, which can be observed in the evolution of the annual mean temperature in Slovakia between 1971 and 2020. Temperature trends in these years are comparable to those of the Visegrad Four. Currently, this development can hardly be slowed [7,8]; therefore, it is necessary to analyze the phenomena that arise as a result of these changes. Rigid pavements are increasingly loaded with different types of traffic. Vehicle loads induce stresses and deformations in the cement-concrete slab (CCS). The CCS is influenced by both the traffic load and the temperatures of the individual According to European Environment Agency (EEA), an agency of the EU, data from 1976 to 2006 [6], the following average warming values occurred in CE: annual mean 0.4 to 0.6 °C, summer 0.2 to 0.4 °C, winter 0.6 to 0.8 °C (Figure 2). The EEA provides sound, independent information on Europe's environment through the provision of timely, targeted, relevant, and reliable information to policymaking agents [17].  [6]. Copyright 2021 by the European Environment Agency.
We note with grave concern the significant gap between the aggregate effect of parties' mitigation pledges in terms of global annual emissions of greenhouse gases by 2020 and aggregate emission pathways consistent with having a likely chance of holding the increase in global average temperature below 2 degrees Celsius or 1.5 degrees Celsius above pre-industrial levels [27]. High concrete temperatures increase the rate of hydration, thermal stresses, permeability, and the tendency towards dry shrinkage cracking and decrease long-term concrete strength and durability as a result of said cracking. Data analysis from the Texas Rigid Pavement database showed that these disorders occur especially in the case of using aggregates: limestone and siliceous river gravel [28].
The results of the analysis emphasize the importance of concrete temperature control during concrete pavement construction in hot weather conditions. Most states specify a  According to European Environment Agency (EEA), an agency of the EU, data from 1976 to 2006 [6], the following average warming values occurred in CE: annual mean 0.4 to 0.6 °C, summer 0.2 to 0.4 °C, winter 0.6 to 0.8 °C (Figure 2). The EEA provides sound, independent information on Europe's environment through the provision of timely, targeted, relevant, and reliable information to policymaking agents [17].  [6]. Copyright 2021 by the European Environment Agency.
We note with grave concern the significant gap between the aggregate effect of parties' mitigation pledges in terms of global annual emissions of greenhouse gases by 2020 and aggregate emission pathways consistent with having a likely chance of holding the increase in global average temperature below 2 degrees Celsius or 1.5 degrees Celsius above pre-industrial levels [27]. High concrete temperatures increase the rate of hydration, thermal stresses, permeability, and the tendency towards dry shrinkage cracking and decrease long-term concrete strength and durability as a result of said cracking. Data analysis from the Texas Rigid Pavement database showed that these disorders occur especially in the case of using aggregates: limestone and siliceous river gravel [28].
The results of the analysis emphasize the importance of concrete temperature control during concrete pavement construction in hot weather conditions. Most states specify a  [6]. Copyright 2021 by the European Environment Agency.
We note with grave concern the significant gap between the aggregate effect of parties' mitigation pledges in terms of global annual emissions of greenhouse gases by 2020 and aggregate emission pathways consistent with having a likely chance of holding the increase in global average temperature below 2 degrees Celsius or 1.5 degrees Celsius above preindustrial levels [27]. High concrete temperatures increase the rate of hydration, thermal stresses, permeability, and the tendency towards dry shrinkage cracking and decrease long-term concrete strength and durability as a result of said cracking. Data analysis from the Texas Rigid Pavement database showed that these disorders occur especially in the case of using aggregates: limestone and siliceous river gravel [28].
The results of the analysis emphasize the importance of concrete temperature control during concrete pavement construction in hot weather conditions. Most states specify a maximum concrete temperature at placement to mitigate the detrimental effects of placement during hot weather. Changes in ambient air temperature during concrete placement are associated with the risk of early cracking. This risk is much greater when temperatures drop from 21 • C to 7 • C than when temperatures drop from 38 • C to 24 • C. Further information can be found in [29]. Concern has been expressed in Florida that, because of a nonlinear temperature gradient in a Portland cement concrete (PCC) pavement, internal stresses could be developed such that the life of the pavement would be seriously reduced. These results were compared with those obtained from the AASHO Test Road and with Bergstrom's prediction method. The results indicated that the nonlinearity of the temperature gradient in a PCC pavement did not have a significant impact on its performance [30]. The territory of Slovakia is in CE at the interface of the parts of the continent that are influenced by continental and oceanic climates. Apart from traffic load, a change in climatic conditions is one of the constant external factors with adverse effects on the physical and mechanical properties of different structural multi-layers in pavements. The 20th century CE temperature increase (+1.2 • C) evolved stepwise with the first peak near 1950 and the second increase (1.3) starting in the 1970s [31].
Since the 19th century, Slovakia has experienced a growth in its annual average temperature of about 1.5 • C, precipitation changes, decreased relative humidity, and changes in solar radiation [7].  [12,15] presents the development and the present understanding of the territory of Central Europe (for clarity, this part is highlighted in Figure 3) for which the authors apply objective research results. Many geographical terms either lack or come to lack a precise meaning and the term Central Europe is a typical example. The development of the perception of Central Europe (Middle Europe, Mitteleuropa, Europe Centrale, Central Europe, etc.) according to various authors is shown in Figure 3a.

Delimitation of the Territory of Central Europe
maximum concrete temperature at placement to mitigate the detrimental effects of placement during hot weather. Changes in ambient air temperature during concrete placement are associated with the risk of early cracking. This risk is much greater when temperatures drop from 21 °C to 7 °C than when temperatures drop from 38 °C to 24 °C. Further information can be found in [29]. Concern has been expressed in Florida that, because of a nonlinear temperature gradient in a Portland cement concrete (PCC) pavement, internal stresses could be developed such that the life of the pavement would be seriously reduced. These results were compared with those obtained from the AASHO Test Road and with Bergstrom's prediction method. The results indicated that the nonlinearity of the temperature gradient in a PCC pavement did not have a significant impact on its performance [30]. The territory of Slovakia is in CE at the interface of the parts of the continent that are influenced by continental and oceanic climates. Apart from traffic load, a change in climatic conditions is one of the constant external factors with adverse effects on the physical and mechanical properties of different structural multi-layers in pavements. The 20th century CE temperature increase (+1.2 °C) evolved stepwise with the first peak near 1950 and the second increase (1.3) starting in the 1970s [31].
Since the 19th century, Slovakia has experienced a growth in its annual average temperature of about 1.5 °C, precipitation changes, decreased relative humidity, and changes in solar radiation [7].

Delimitation of the Territory of Central Europe
Figure 3 [12,15] presents the development and the present understanding of the territory of Central Europe (for clarity, this part is highlighted in Figure 3) for which the authors apply objective research results. Many geographical terms either lack or come to lack a precise meaning and the term Central Europe is a typical example. The development of the perception of Central Europe (Middle Europe, Mitteleuropa, Europe Centrale, Central Europe, etc.) according to various authors is shown in Figure 3a. The graphic borders of Central Europe, shown in Figure 3a, were available only in some cases. In other cases, they are drawn based on written descriptions in relevant texts. The location of the geographical center of Europe depends on the definition of the borders of Europe, mainly whether remote islands are included to define the extreme points of The graphic borders of Central Europe, shown in Figure 3a, were available only in some cases. In other cases, they are drawn based on written descriptions in relevant texts. The location of the geographical center of Europe depends on the definition of the borders of Europe, mainly whether remote islands are included to define the extreme points of Europe, and on the method of calculating the final result. Thus, several places claim to house this hypothetical center: the town of Kremnické Bane or the neighboring town Krahule, near Kremnica, in Central Slovakia; the small town of Rakhiv or the town of Dilove, near Rakhiv, in Western Ukraine; the town of Girija, near Vilnius, in Lithuania; a point on the island of Saaremaa in Estonia; a point near Polotsk, or in Vitebsk, or near Babruysk, or the nearby Land 2022, 11, 764 5 of 24 lake Sho in Belarus; and a point near the town of Tallya, in Northeastern Hungary ( Figure 4). The towns of Krahule and Kremnicke Bane are generally understood by the Slovak authors as the geometric center of Continental Europe, as well as by foreign authors [13,14]. house this hypothetical center: the town of Kremnické Bane or the neighboring town Krahule, near Kremnica, in Central Slovakia; the small town of Rakhiv or the town of Dilove, near Rakhiv, in Western Ukraine; the town of Girija, near Vilnius, in Lithuania; a point on the island of Saaremaa in Estonia; a point near Polotsk, or in Vitebsk, or near Babruysk, or the nearby lake Sho in Belarus; and a point near the town of Tallya, in Northeastern Hungary ( Figure 4). The towns of Krahule and Kremnicke Bane are generally understood by the Slovak authors as the geometric center of Continental Europe, as well as by foreign authors [13,14].
In the following sections, the development of the average annual air temperature at 10 climatological stations in Slovakia as a Central European country is presented for this reason.  Regarding research on correlation dependencies, we briefly describe the basic statistics-correlation coefficient and coefficient of determination. In mathematics, the Spearman correlation coefficient is usually used to quantify how much two columns of data linearly depend on each other. Charles Edward Spearman (1863-1945) was an English In the following sections, the development of the average annual air temperature at 10 climatological stations in Slovakia as a Central European country is presented for this reason. Regarding research on correlation dependencies, we briefly describe the basic statisticscorrelation coefficient and coefficient of determination. In mathematics, the Spearman correlation coefficient is usually used to quantify how much two columns of data linearly depend on each other. Charles Edward Spearman (1863-1945) was an English psychologist known for his work on statistics as a pioneer of factor analysis and Spearman's rank correlation coefficient. The values of Spearman's rank correlation coefficient according to Table 1 are most often mentioned in the technical, as well as behavioral, literature [32][33][34]. Karl Pearson built on Francis Galton's research, and Auguste Bravais developed and published the mathematical formula in 1844 [35][36][37].  [38,39].

Spearman Correlation Coefficient
where cov(X, Y) is the covariance; σ X and σ Y are the standard deviations of X and Y, respectively; n is the sample size; X i , Y i are the individual sample points; and X, Y are the sample means. The correlation coefficient ranges from −1 to 1. The value R(X, Y) = 0 implies that there is no linear dependency between variables (X, Y). A relatively high correlation coefficient means that there is a high linear dependency between the variables. This does not necessarily mean that there is also a high causal dependency. The degree of causal dependence is expressed by the coefficient of determination R 2 , a key output in regression analysis. It is the square of the correlation coefficient between variables based on the sample values. It gives valid results when the observations are evaluated correctly without measurement errors. Based on the coefficient of determination, we can assess to what extent the regression model fits the observed data.

Correlation Dependencies of Average Annual Temperatures T a for the Period 1971-2020
In this section, the authors present the latest research results in the field of the objective correlation dependences of the average annual temperature from the altitude of the assessed pavement. For practical purposes, the flow of daily temperature is expressed by the average daily temperature, calculated as shown in Equation (2): where indexes 7, 14, and 21 are the times at which air temperature T is measured. Considering the standard in Equation (2), the average daily air temperature is calculated as the average of the four measured values T 7 , T 14 , and twice the value of T 21 . Twice the value of T 21 is included because there are no night measurements available [40]. The temperature of earth structures changes along with the change in air temperature. Average annual air temperature T a is expressed by the formula: Air temperature has a cyclical character repeated in daily and annual cycles with an approximately sinusoidal shape. Figure 4 presents the mean annual temperatures T a of  Table 2.
of earth structures changes along with the change in air temperature. Average annual air temperature Ta is expressed by the formula: Air temperature has a cyclical character repeated in daily and annual cycles with an approximately sinusoidal shape. Figure 4 presents the mean annual temperatures Ta of 11 climatological stations with significantly different altitudes; from Hurbanovo (115 m asl) and Bratislava (131 m asl) through Žilina (365 m asl) to Oravska Lesna (858 m asl) for the period from 1971 to 2020 and from 1971 to 2017. Figure 5 presents the correlation dependences of the development of average annual air temperatures Ta from year Y for the period 1971-2017 and 1971-2020. According to the Spearman correlation coefficient, all identified dependencies, with the exception of Tatranska Lomnica, show a strong degree of correlation. When evaluated according to Pearson's correlation coefficient for altitudes above 350 m, correlation dependencies show high positive correlation and, below 350 m, moderate positive correlation. An overview of the correlation coefficients of the development dependence of the considered 10 climatic stations in Slovakia for the period from 1971 to 2020 (2017) is given in Table 2.    Table 2 present the results of the author's 20-year research in the field of prediction of the development of the average annual air temperature as one of the limiting factors in the assessment of cement concrete covers of rigid pavements. Table 2 Periods 1971Periods -2000Periods , 1971Periods -2010, and 1971-2020 Figure 5 and Table 2 present the results of the author's 20-year research in the field of prediction of the development of the average annual air temperature as one of the limiting factors in the assessment of cement concrete covers of rigid pavements. Table 2

Design Principles of CCP Design According to Slovak Road Act
There are three conventional types of cement-concrete pavements (CCP): jointed plain concrete pavements (JPCPs; Figure 7a,b), jointed reinforced concrete pavements

Design Principles of CCP Design According to Slovak Road Act
There are three conventional types of cement-concrete pavements (CCP): jointed plain concrete pavements (JPCPs; Figure 7a,b), jointed reinforced concrete pavements (JRCPs; Figure 7c), and continuously reinforced concrete pavements (CRCPs; Figure 7d). All three common rigid pavements carry traffic loading through the flexural strength of the concrete. JPCPs are the most common type of rigid pavement due to their cost and simplicity in the United States of America [29,41], as well as in tunnel pavements in Slovakia [10]. Examples of the employment of JPCP in other countries-Moldova and Romania-are presented in [42,43].
Typically, a flexible (or asphalt) pavement is produced by continuously reducing the stiffness parameters of the individual layers to minimize stresses and deformations within the road body. The stiffest part of rigid pavements is the slab, under which layers with lower stiffness parameters can be applied. Furthermore, even when the stress in the rigid slab is increased, the overall deformation of the structure tends to decrease. The materials and thicknesses are chosen based on site-specific weather conditions.  Figure 7c), and continuously reinforced concrete pavements (CRCPs; Figure 7d). All three common rigid pavements carry traffic loading through the flexural strength of the concrete. JPCPs are the most common type of rigid pavement due to their cost and simplicity in the United States of America [29,41], as well as in tunnel pavements in Slovakia [10]. Examples of the employment of JPCP in other countries-Moldova and Romania-are presented in [42,43]. Typically, a flexible (or asphalt) pavement is produced by continuously reducing the stiffness parameters of the individual layers to minimize stresses and deformations within the road body. The stiffest part of rigid pavements is the slab, under which layers with lower stiffness parameters can be applied. Furthermore, even when the stress in the rigid slab is increased, the overall deformation of the structure tends to decrease. The materials and thicknesses are chosen based on site-specific weather conditions. Under CE conditions, the pavement is perceived as a paved part of the road intended for vehicles and consisting of a surface, a road foundation, and a capping layer. Cementconcrete pavement (CCP) represents a rigid pavement, i.e., a pavement substantially constructed of high-strength cement concrete [44]. The basis of the design of the CCP structure is the selection of materials for the layers of the base system, the design of the thickness of these layers, and the design of the dimensions of the CCS, their thickness, width, and length. The currently relevant regulations for the design of CCP is technical standard TP 098 [45], which specifies the procedure for designing, calculating, and assessing the structures of pavements with CC cover (wearing course) in Slovakia. The technical conditions for the design of CCP in the Slovak Republic, TP 098, are intended for investment organizations, designers, and construction organizations, as well as for the administration and maintenance of roads. The design of the structure, the choice of materials, as well as the criteria by which the dimensions of the pavements are assessed are differentiated according to the size of the traffic load and the type of cement-concrete cover.
The principles of designing roads with CCP in the Czech Republic are very similar to those in Slovakia. For the design of cement concrete coverings, the average annual air temperature is determined according to Annex A of ČSN 73 6114: 1995 [46]. This can be precisely determined from the data of meteorological stations in the given area for a mean return period of 10 years. Under CE conditions, the pavement is perceived as a paved part of the road intended for vehicles and consisting of a surface, a road foundation, and a capping layer. Cementconcrete pavement (CCP) represents a rigid pavement, i.e., a pavement substantially constructed of high-strength cement concrete [44]. The basis of the design of the CCP structure is the selection of materials for the layers of the base system, the design of the thickness of these layers, and the design of the dimensions of the CCS, their thickness, width, and length. The currently relevant regulations for the design of CCP is technical standard TP 098 [45], which specifies the procedure for designing, calculating, and assessing the structures of pavements with CC cover (wearing course) in Slovakia. The technical conditions for the design of CCP in the Slovak Republic, TP 098, are intended for investment organizations, designers, and construction organizations, as well as for the administration and maintenance of roads. The design of the structure, the choice of materials, as well as the criteria by which the dimensions of the pavements are assessed are differentiated according to the size of the traffic load and the type of cement-concrete cover.
The principles of designing roads with CCP in the Czech Republic are very similar to those in Slovakia. For the design of cement concrete coverings, the average annual air temperature is determined according to Annex A ofČSN 73 6114: 1995 [46]. This can be precisely determined from the data of meteorological stations in the given area for a mean return period of 10 years.

Assessment of Mechanical Efficiency of Cement Concrete Slab
According to Act No. 135/1961 Coll. on Roads (Road Act), the design of the pavement is carried out in accordance with valid Slovak technical standards, technical regulations, objectively determined results of research and development for road infrastructure, or similar technical specifications [44]. According to this act, the structural safety design of a CCS must satisfy condition (4). The limiting criterion is the assessment for a single load [45].
Here, R i,n (MPa) is the flexural tensile strength of the cement-concrete, σ max,OT (MPa) is the maximum flexural tensile stress from a single load on the longitudinal or transverse joint of the CCS, σ DA,115kN is the stress due to pure force load (design axle), and σ T represents stress due to thermal load [45]. The safety ratio between its strength and its actual stress response is ξ i . This value is necessary to establish when assessing a cement-concrete pavement (ξ i > 1.10 for highways and ξ i > 1.05 for Class 1 to III roads). The maximum flexural tensile stress σ max,OT must be greater than or equal to the flexural tensile strength of the cement concrete R i,n at the specified utilization factor. Usually, the highest value of the maximum stress is achieved at the transverse or longitudinal joint of the CCS and therefore the correct design for these two load positions is assessed in Figure 8.
CCS must satisfy condition (4). The limiting criterion is the assessment for a single load [45].
Here, Ri,n (MPa) is the flexural tensile strength of the cement-concrete, σmax,OT (MPa) is the maximum flexural tensile stress from a single load on the longitudinal or transverse joint of the CCS, σDA,115kN is the stress due to pure force load (design axle), and σT represents stress due to thermal load [45]. The safety ratio between its strength and its actual stress response is ξi. This value is necessary to establish when assessing a cement-concrete pavement (ξi > 1.10 for highways and ξi > 1.05 for Class 1 to III roads). The maximum flexural tensile stress σmax,OT must be greater than or equal to the flexural tensile strength of the cement concrete Ri,n at the specified utilization factor. Usually, the highest value of the maximum stress is achieved at the transverse or longitudinal joint of the CCS and therefore the correct design for these two load positions is assessed in Figure 8. In the contraction joints and at the edge of the cement concrete slab, extreme normal stresses are expected. Therefore, the load is positioned at two critical points by the rear axle of the load vehicle. We consider one position in which the rear axle is most closely located to the intersection of the transverse and longitudinal contraction joints. It is always  In the contraction joints and at the edge of the cement concrete slab, extreme normal stresses are expected. Therefore, the load is positioned at two critical points by the rear axle of the load vehicle. We consider one position in which the rear axle is most closely located to the intersection of the transverse and longitudinal contraction joints. It is always possible to consider the second position to be either the intersection of the transverse contraction joint or closest to the edge of the cement concrete slab. The current legislation requires the calculation to be conducted using the principles of FEM.
The stresses from temperature change effects (σ T,x(y) also TS temperature value ) can be solved using Equation (5). In this equation, E CC,T is the modulus of elasticity for permanent load, α is the coefficient of thermal expansion, C x ( y ) (−) are the values of the coefficients (which are identified from Figure 9, according to TP 098 [45]), and h c (m) is the critical design parameter-proposed slab thickness. Empirical value 12.44 is in • C and multiplication constant 28 has a physical dimension • C/m. The reduced modulus of elasticity of cement concrete is considered in the calculation of stresses induced by temperature loading. For longer time actions, this modulus of elasticity corresponds to its curing at E CC,T = 27,500 (MPa) level. The coefficient of thermal expansion is α = 0.000012 ( • C −1 ).
crete is considered in the calculation of stresses induced by temperature loading. For longer time actions, this modulus of elasticity corresponds to its curing at ECC,T = 27,500 (MPa) level. The coefficient of thermal expansion is α = 0.000012 (°C −1 ).
= √ , ⋅ ℎ 3 12 ⋅ (1 − 2 ) ⋅ 4 (6) Equation (6) determines the radius of relative stiffness under thermal stress lT (m). ν is the Poisson's constant, hc (m) is the proposed slab thickness, ECC,T is the modulus of elasticity for permanent load, and k is the modulus of reaction of the pavement system. The radius of relative stiffness under thermal stress lT is an important value in the graph for the determination of the Cx (Cy) coefficient (Figure 9), the dimensions of the CCS Lx and Ly (m) and the Poisson's constant ν = 0.2. The modulus of reaction of the pavement system k (MN·m −3 ) is recommended to be between 100 and 300 MN·m −3 . Figure 9. Graph for the determination of the Cx (Cy) coefficient according to TP 098. Reprinted with permission from [45]. Copyright 2022 by the Slovak road administration.

Determination of the Average Annual Air Temperature for the Calculation of Thermal Stresses
It is stated in TP 098 that the design of pavement with a cement-concrete cover must comply with the following: road protection against adverse effects of subgrade freezing, the ratio of concrete strength and tensile stress at bending in CCS due to one-time stress, the ratio of concrete flexural strength and tensile stress at bending in CCS due to repeated stress. Due to the great influence of temperature on deformation and stresses in CCSs (thermal stresses), it is always necessary to consider the combination of traffic load and the influence of the temperature regime of CCP in the calculations. Figure 9. Graph for the determination of the C x (C y ) coefficient according to TP 098. Reprinted with permission from [45]. Copyright 2022 by the Slovak road administration. Equation (6) determines the radius of relative stiffness under thermal stress l T (m). ν is the Poisson's constant, h c (m) is the proposed slab thickness, E CC,T is the modulus of elasticity for permanent load, and k is the modulus of reaction of the pavement system. The radius of relative stiffness under thermal stress l T is an important value in the graph for the determination of the C x (C y ) coefficient (Figure 9), the dimensions of the CCS L x and L y (m) and the Poisson's constant ν = 0.2. The modulus of reaction of the pavement system k (MN·m −3 ) is recommended to be between 100 and 300 MN·m −3 .

Determination of the Average Annual Air Temperature for the Calculation of Thermal Stresses
It is stated in TP 098 that the design of pavement with a cement-concrete cover must comply with the following: road protection against adverse effects of subgrade freezing, the ratio of concrete strength and tensile stress at bending in CCS due to one-time stress, the ratio of concrete flexural strength and tensile stress at bending in CCS due to repeated stress. Due to the great influence of temperature on deformation and stresses in CCSs (thermal stresses), it is always necessary to consider the combination of traffic load and the influence of the temperature regime of CCP in the calculations.
At the time of processing the submitted paper, to determine the average annual air temperature, together with the traffic loads of the decisive characteristics of the CCP assessment, STN 73 6114:1997 [40] was use, which contains the map presented in Figure 10.
For assessment, according to TP 098 [45], we must know the average annual air temperature T a . Figure 11 shows the correlation dependences of the air temperatures of interest on the altitude SL of CCP found from 1971 to 2017 (2011). Specifically, the dependences include T max -the maximum measured air temperature at 14 h for the period 1971-2011, T d, max -the maximum average daily air temperature for the period 1971-2011, T a,2021 -the predicted average annual temperature in 2021 determined according to the equations in Figure 5, and T a, LLTZ -the lower limit of the temperature zone determined from the map STN 736114: 1997 on the SL of the assessed pavement.
At the time of processing the submitted paper, to determine the average an temperature, together with the traffic loads of the decisive characteristics of the C sessment, STN 73 6114:1997 [40] was use, which contains the map presented in Fig   Figure 10. For assessment, according to TP 098 [45], we must know the average annual perature Ta. Figure 11 shows the correlation dependences of the air temperatures o est on the altitude SL of CCP found from 1971 to 2017 (2011). Specifically, the depen include Tmax-the maximum measured air temperature at 14 h for the period 197 Td, max-the maximum average daily air temperature for the period 1971-2011, Ta,2 predicted average annual temperature in 2021 determined according to the equa Figure 5, and Ta, LLTZ-the lower limit of the temperature zone determined from t STN 736114: 1997 on the SL of the assessed pavement.     Table 2. Reprinted with permission from [40]. Copyright 2019 by the Slovak Office of Standards, Metrology and Testing.
For assessment, according to TP 098 [45], we must know the average annual air temperature Ta. Figure 11 shows the correlation dependences of the air temperatures of interest on the altitude SL of CCP found from 1971 to 2017 (2011). Specifically, the dependences include Tmax-the maximum measured air temperature at 14 h for the period 1971-2011, Td, max-the maximum average daily air temperature for the period 1971-2011, Ta,2021-the predicted average annual temperature in 2021 determined according to the equations in Figure 5, and Ta, LLTZ-the lower limit of the temperature zone determined from the map STN 736114: 1997 on the SL of the assessed pavement.   Figure 12 presents detailed results of sensitivity analyses of calculated thermal stress values for CCSs thicknesses of 12 to 32 cm and the modulus of the subgrade reaction values k = 100, 200, and 300 MN·m −3 determined by the procedure given in TP 098. The modulus of subgrade reaction k is used as a primary input for rigid pavement design. It estimates the support of the layers below a rigid pavement surface course (the PCCS). The k-value can be determined by field tests or by correlation with other tests. There is no direct laboratory procedure for determining the k-value. The modulus of subgrade reaction came about because work conducted by Westergaard during the 1920s led to the development of the k-value as a spring constant to model the support beneath the slab [24]. value can be determined by field tests or by correlation with other tests. There is no direct laboratory procedure for determining the k-value. The modulus of subgrade reaction came about because work conducted by Westergaard during the 1920s led to the development of the k-value as a spring constant to model the support beneath the slab [24].

Calculation of the Response of a Real Cement Concrete Pavement in the Conditions of Slovakia
For the design and analysis of cement concrete pavements, there are a number of methods and processes currently available. In many developed countries around the world, numerical modelling is the most appropriate method. The finite element method (FEM) is one of the most effective general-purpose, numerical, variation methods for solving continuum mechanics boundary value problems. It is universal to the type, geometry, boundary conditions, and loading of the structure. The deformation variant is the most commonly used solution. Slovak legislation even requires pavements to be analyzed using FEM. Different types of FEM models are suitable for different application purposes, including the calculation of static and dynamic responses by means of their shapes, boundary conditions, types of finite elements, types of contacts, etc. By examining the analytical and experimental results that were presented, the FEM model is supposed to investigate the validity of its use in practice [47,48].
An example of a real pavement numerical analysis using FEM is presented in this section. Based on the above results of the sensitivity analysis of the calculated thermal stress values for CCSs with thicknesses from 12 to 32 cm, the thickness of the CCS was determined. This is approximately the average value in the range (19 cm). In this example, the temperature is not directly included in the stress calculation. By comparing both sections' analyses, we can determine the influence of temperature. The CCS is the most important part of the model. In the model, only the contact stresses in the wheel tire on the slab surface are considered.
It has been determined which section of the roadway corresponds to the boundary conditions for the CC pavement thickness set out above. It has been selected as a plausible illustration of the use of numerical modelling in practice. This is the existing real section of the motorway near the town of Žilina, before the Ovčiarsko tunnel. In both sections, before and inside the tunnel, the proposed pavement composition is the same.

Solution Method and Numerical System
For this example, the IDA NEXIS system was used. It was chosen for its simplicity and user-friendliness. As a unique program, it overcomes barriers for ordinary users. For common purposes, its use is simple, natural, and intuitive and it provides diverse and rich features. The program uses the deformation variant of FEM to deal with 1D and 2D elements, plane strain, material and geometric nonlinearities, dynamics, stability analysis, and geotechnics. It is used to numerically analyze problems in civil engineering [49]. In this FEM system, the following steps are necessary in order to create a credible model: • Stage 1 = Preprocessing-Preprocessing is performed for the preparation of the modelcontrol using FEM. It includes the generation of the mesh, either surface or volume, and the following operations: automatic meshing, mesh crossing, operations (shifting, copying, mirroring, etc.), and mesh clipping. It defines the geometry of the finite element mesh, the individual boundary conditions, and the load. • Stage 2 = Processing-The FEM equations are processed and solved. The system creates the equations and stiffness matrices using the data from the first stage automatically. This numerical procedure is automatically processed by a computer and does not require a user interface. • Stage 3 = Postprocessing-Visualization and graphical display of the results of the entire analysis. These results can be presented numerically or as a video simulation.
The form of the feature function generally depends on the shape and type of the individual analysis [50]. Element convergence in FEM is also related to the element form characteristics of the function. This is very useful information for understanding heterogeneous shape functions and their properties. When developing a model, it is important to know when the elements of the structure interact with each other and when they behave as separate parts. This interaction can be identified for different variants of the model: For numerical pavement models, the interaction with the subgrade is most important. The soil-in module in NEXIS allows us to solve this type of problem. It includes the conversion of the geotechnical characteristics of the subgrade to stiffness C i compressibility parameters. The calculation is performed on a nonlinear basis using iterations converging on the final response solution.
The module calculates the C i parameters for the interaction between the foundation slab and the subgrade, considering the load distribution and intensity, the contact stress at the structure-soil interface, the foundation geometry, and the local geological conditions. The use of this module also offers other advantages such as multi-parametric interaction between the foundation slab and the subgrade, subgrade input using data obtained from exploratory boreholes, and others [51].

Numerical Calculation of Real CCP Structure in High Altitude
The location of this section ideally corresponds to the topic of the presented article. The Ovčiarsko tunnel, located on the D1 motorway in the section Hričovské Podhradie-Lietavská Lúčka, is a two-tube highway tunnel belonging to the medium-length tunnel category (2360 m). It forms part of the Žilina bypass. On motorways, due to their longer service life and lower maintenance costs, cement concrete pavements are designed and mechanically considered as thin slabs (Kirchhoff) with the elastic properties of reinforced or unreinforced concrete, which must satisfy the basic criteria for assessing the effects of cyclic loading by truck passes. As part of the design process, it is always advisable to assess the structural response of the various design options in terms of the structural layers [52].

Mechanical Characteristics of the Model
The main, stiffest part of the pavement structure of the model consists of a 190 mmthick CB I cement concrete slab. This is followed by additional layers-40 mm thick AC 16 B asphalt concrete for the roadbed layer, under which is 180 mm thick CBGM C5/6 cementitious stabilization. The last 200 mm thick layer is ND gravel, under which the sub-base is placed. The total thickness of the modelled pavement is 610 mm. The composition of the structural layers and the summarized mechanical properties used for the modelled elements can be seen in Table 3. The IDA NEXIS computational system allows the individual layers of the computational model to be modelled using the soil-in module, which allows the individual mechanical properties to be entered separately for each layer. The layer thicknesses are entered relative to the bottom edge of the cement concrete slab [53].

Considered Loads and Model Geometry
The numerical analysis considers the maximum realistic load, which is assessed using the basic criteria given in TP 098 [45]. The modelled pavement was loaded at the wheel pressure points of the front and rear axle of a TATRA T815-S3 26 208 truck [54]. This is a triple-folding truck, designed primarily for the transport of bulk materials up to 15,700 kg ( Figure 13). The effect of the rear double axle is modelled as the wheel pressure of one rear axle in accordance with the regulation.

Results of Numerical Analysis of the CCP Highway Section Using FEM
The results of the numerical calculations are represented by the response of the pavement structure of the section in front of the tunnel to the above-specified static loads. This response can be divided into terms of stresses and strains. Both are then related to the assessment of the first and second limit states, as is generally the case for building structures. The results show two variants of the solution-Variant 1, loading at the outer edge of the modelled region and Variant 2, loading at the inner expansion ( Figure 14).
Since the CCP cover slab is the most important part of the model, only the contact stresses of the wheel tire on the slab surface are considered. Because they are assumed to be linear after the thickness of the slab (see elasticity theory). These stresses are represented along the median surface of the slab in the x-axis direction along the length-σx and the y-axis direction across-σy. Another surface stress is the stress τxy, which represents the shear stress in the horizontal plane. Figure 15 shows the waveforms of all three of these stresses σx, σy, and τxy for the variants considered, as well as the contact stresses between the CCP slab and the asphalt concrete for the surface layers σz. Another important parameter in the modelling task is the contact area between the wheel and the road, which is in fact not an ideal rectangle but an idealized one. This area depends on several parameters. It is influenced by the weight exerted by the vehicle on the wheel, the speed of the moving vehicle, the pressure and type of tire used, the road roughness, the coefficient of friction, etc. Its determination is dealt with as concisely as possible in the literature [54]. This problem can also be solved by numerical simulations based on FEM. As an example, a computational model with the contact stress distribution of the tire-pavement in ABACUS is presented (Figure 13). For a standard vehicle tire width of 300 mm, a con-contact area length of 100 mm was considered for the front axle and 150 mm for the rear axle. The modelled highway section has a length of 5 × 6 m and a width of 2 × 4 m (contraction line continuous joints bordered). From the axle loads of each axle, the load on the contact patch in the vertical direction was expressed and then the load obtained was distributed evenly over the contact patch of the individual wheels. This is the wheel pressure in kN/m 2 .

Results of Numerical Analysis of the CCP Highway Section Using FEM
The results of the numerical calculations are represented by the response of the pavement structure of the section in front of the tunnel to the above-specified static loads. This response can be divided into terms of stresses and strains. Both are then related to the assessment of the first and second limit states, as is generally the case for building structures. The results show two variants of the solution-Variant 1, loading at the outer edge of the modelled region and Variant 2, loading at the inner expansion ( Figure 14). Since the CCP cover slab is the most important part of the model, only the contact stresses of the wheel tire on the slab surface are considered. Because they are assumed to be linear after the thickness of the slab (see elasticity theory). These stresses are represented along the median surface of the slab in the x-axis direction along the length-σ x and the y-axis direction across-σ y . Another surface stress is the stress τ xy , which represents the shear stress in the horizontal plane. Figure 15 shows the waveforms of all three of these stresses σ x , σ y , and τ xy for the variants considered, as well as the contact stresses between the CCP slab and the asphalt concrete for the surface layers σ z .
In terms of deformations, Figure 16 shows the vertical displacement or deflections of the slab and the rotations of the tangent surfaces. All outputs are in the form of isosurfaces. The color scale on the sides of the figures shows the maxima and minima of the individual values. The final table (Table 4) summarizes the results of the static calculation using FEM.
The following partial conclusions are evident from both modelled options:

1.
Surface normal stresses σ x -take their highest values on the underside of the slab in the axle area and gradually disappear at a distance of about 1 m from the axle. Values at greater distances are negligible. Their progression in the transverse direction is approximately constant between the wheels.

2.
Surface normal stresses σ y -take on significant peaks directly at the contact surface of the wheels. Their values are the most pronounced and most influenced by the principal normal stresses. They decrease at very small distances and take negligible values.

3.
Surface shear stresses τ xy -change sign between axles and also between wheels. They increase towards the force and decrease to zero when moving away from the force. Their values are the least significant.

4.
Centre-plane deflections-it is evident from the deflections that the largest deflections are directly under the wheels of the vehicle. The values gradually decrease in the transverse direction between the wheels and also in the longitudinal direction between the axles. They disappear approximately 1 m in front and 1 m behind the vehicle.

5.
Rotations of the tangent surfaces-about the x, y axes, these values are related to the deflections and the deflection surface. They correspond to deflections because they are their first derivative.

Discussion
The authors present the results of 25 years of research in the field of adaptability of transport structures, especially road pavements, to climate change in Slovakia, located in CE. Based on the analysis of the development and current state of perception of Central Europe, for the purposes of this article, they defined its schematic territory according to

Discussion
The authors present the results of 25 years of research in the field of adaptability of transport structures, especially road pavements, to climate change in Slovakia, located in CE. Based on the analysis of the development and current state of perception of Central Europe, for the purposes of this article, they defined its schematic territory according to Figure 4. The article presents a linear and polynomic correlation dependence of average annual temperature T a ( • C) on the SL altitude above sea level (m) of the designed pavements of Slovakia evaluated for the periods 1971-2000, 1971-2010, and 1971-2020. For technical practice purposes, we recommend using the following equation determined from the period 1971 to 2020 (Figure 11).
T a,1971−2020 = −0.0067 × SL + 11.24 (7) The average annual temperatures of 10.6 to 3.8 • C correspond to the determined range of altitudes inhabited by the Central European territories of 100 to 1100 m. For the purposes of judging Central European road pavements, the following dependencies were specifically used ( Table 5)  It is clear from Table 5 that climatic stations (CS) at airports show a significantly higher increase in T a (2.7 • ) than other CS (1.9 • ), which is probably caused by large areas of asphalt or cement-concrete pavements. Between 1971 and 2020, an average increase in T a of 2.1 • C was found in Slovakia, which is about 0.7 • C higher (Figure 1) than the average increase reported by the Intergovernmental Panel on Climate Change (IPCC). By comparing the research results presented in this paper and [7,46] with the findings of the European Environment Agency from 1976 to 2006 (Figure 2), full compatibility of the outputs was identified. For the period 1976 to 2006, the EEA presents an average increase T a in Central Europe of 1.2 to 1.5 • C and the authors found a value of 1.3 • C in Slovakia ( Table 5).
The accuracy of the CCP design is currently limited in Slovakia due to the accuracy of the determination of T a , since the effects of traffic load can be calculated with high accuracy by means of FEM (Section 5). The important conclusion is Figure 11, where the authors present the dependences of thermal stresses from a single load of 115 kN designed for thicknesses from 12 to 32 cm of cement-concrete cover h c on the annual average air temperature T a from 4.5 to 11.5 • C for the modulus of subgrade reaction k = 100, 200, and 300 MN·m −3 [55].
To illustrate the significant effect of T a on thermal stresses (TS), we give corresponding numerical values for h c = 20 cm and dimensions of CCSs 4 × 4 m. The summarized results are shown in Table 6 (TS Ta=4.5 = TS 4.5 ). The above results indicate that the analytical solution for the stress profile of the standard load and the actual temperature load from the focus research part of the paper is valid. Surface maximum stress values of normal stresses obtained by FEM calculation represent 61% to 82% of all maximum stress values obtained by analytical calculation. Due to this, it is evident that reinforced concrete slabs are the most effective structure regardless of the stiffness of the subgrade. This type of road body structure is less vulnerable to the effect of temperature variations.
Based on the above research results, it can be stated that the current climatic characteristics used for pavement design in Slovakia do not take into account current climate change. The authors are fully aware of the serious negative effects of climate change on pavement structures, but these changes have a positive effect on the structural design. Global warming causes a significant increase in T a , which allows for the structural design of road pavements with lower overall thicknesses, as well as h c in the case of cement-concrete pavements (CCP) [56].
In Central Europe, more attention has been paid in the last 10 years to the development of innovative climate-adaptable asphalt mixture, [57][58][59][60], as well as cement-bound granular mixtures [60,61]. In all presented innovative materials [57][58][59][60][61] temperatures play a significant role in both the mix design and structure design of pavements. For this reason, we actually focus on the implementation of the presented knowledge to the development of composite foam concrete for the base layers of civil engineering structures [62].

Conclusions
The results of temperature change monitoring presented in the substantive sections of this paper are relevant to studies dealing with the effect of temperature change on stress increase in CCS [63][64][65], which are presented in world-renowned and highly cited publications. In terms of conclusions, the results need to be divided into the time-lapse original data, the mathematical analysis, the correlation and approximation of dependencies that led to the predictive models, the numerical-analytical solution of the response to force, and temperature effects of the load. In these results and their detailed analysis, it is possible to find a solution for the specific participation of stress effects in the computational models of this part of the pavement structure. Since this is a structural part that is exposed to weather effects and is also directly impacted by traffic loads, one of the main conclusions is the directly specified proportion of temperature and force loads. We can see a few differences between the standard solution of the stress response presented in the literature [66][67][68] and the findings presented in Section 3, Section 4, and Section 5 of the paper. Taking these comparisons into account, it can be concluded that CE should be considered a specific area, especially when it comes to some of the regions mentioned in Section 2 (in Slovakia). In this context, it is a more sensitive approach to assessing stresses that arise in CCP. In addition, it is also important to define the necessity to modify standards and regulations to take into account location-related phenomena. Another interesting conclusion relates to the potential for creating a predictive mathematical model based on correlations of results between the long-term monitoring of temperature changes (years 1971 to 2020). Further implementation of the model can be undertaken to assess the suitability of CCP use or to define more precisely the impact of temperature on design over a longer time horizon.
The contributions of the paper must be restricted to the scientific community in Central Europe engaged in road construction research. The long-term monitoring of temperature changes in Slovakia has produced significant original results. Throughout the territory of CE, the variability of mountain ranges and plains is affected in different ways by global climate change. This is a characteristic of Slovakia. The results reported in the experimental sections of the paper will be used to conduct research on advanced materials and additives for cement-concrete pavements. One of the main conclusions is that their use is most appropriate in terms of mitigating the effects of temperature on CCPs.
The results of the measurements can also influence the CCP design methodology to introduce mandatory monitoring and diagnostics of pavement structures for specified territories most significantly impacted by temperature changes.