Climate Regionalization of Asphalt Pavement Based on the K-Means Clustering Algorithm

)e climate regionalization of asphalt pavement plays an active role in ensuring the good performance and service life of asphalt pavement. In order to better adapt to the climate characteristics of a region, this study developed a multi-index method of climate regionalization of asphalt pavement. First, meteorological data from the research region were statistically analyzed and the major climate variables were identified. )en, a principal component analysis (PCA) was used to eliminate any correlation between the major climate variables.)ree principal components were extracted by the PCA as cluster factors, namely, the temperature factor, precipitation factor, and radiation factor. )e research region was divided into the following four asphalt pavement climate zones via the K-means clustering algorithm. )ose zones are affected by the climate comprehensively: an inland zone with high temperatures, little rainfall, and radiation, a coastal zone with high temperatures, and a rainy mountainous zone.)e results of the climate regionalization were compared with the results of on-site investigations. )e pavement degradation in each climatic zone was related to the climate characteristics of the region. Probabilistic neural network (PNN) and support vector machine (SVM) climate regionalization predictive models were established with MATLAB. )e clustering factors were used as the input data to identify the climate zones, and the identification accuracy rate was determined to be over 90%. )e climate regionalization of pavement can provide reference and guidance for the selection of reasonable technical measures, parameters, and building materials in highway projects with similar climatic conditions.


Introduction
e performance of asphalt pavement is affected by many environmental factors such as temperature, precipitation, and solar radiation, during its service period. In addition, the bearing capacity and service performance of asphalt pavements are attenuated continuously [1]. e asphalt pavement structure should be coordinated with the natural environment and adapted to the climate characteristics of the region. In this way, typical service performance and service life can be guaranteed [2]. erefore, many researchers try to analyze and study the climate conditions expected during the road construction process and service period. And then, local plans of the climate regionalization of asphalt pavements are proposed.
ese plans are used to guide the selection of asphalt materials, mixture design, and reasonable construction methods for the pavements. In recent years, many countries have focused on the temporal and spatial variability of climate elements and established the standards of the regionalization of asphalt pavement.
e Performance Grade (PG) method and standard for matrix asphalt were proposed by the Strategic Highway Research Program in America. Asphalt is divided into four grades based on the highest average pavement temperature for seven continuous days and the extreme minimum air temperature in a year [3]. Subsequently, other countries began to introduce and carried out native analysis of applicability. Based on the high temperature, low temperature, and precipitation indexes, the method of the climate regionalization for asphalt pavement performance was proposed in China [4]. e climatic regionalization of asphalt under temperature conditions in Poland was established based on the standard of PG [5,6]. Similarly, other countries such as Jordan [7], Italy [8], ailand [9], and Iran [10] have developed the plans of the climate regionalization according to the extreme temperature range that the asphalt will be subjected to. In these methods, a Geographic Information System (GIS) is used to draw intuitive maps of the climate regionalization [11]. Spatial interpolation methods are used to process complex meteorological data [12]. Most standard methods used in previous studies are applied according to temperature and precipitation directly. However, the climate regionalization established only based on few climate factors is obviously not reasonable enough. More climate factors should be considered in the regionalization plans. But there are often correlations between the complex climate factors considered. is will not make the plans of the climate regionalization also adequately reveal the different effects of climatic conditions on asphalt pavement. And the results of the climate regionalization are also relatively inaccurate and cannot completely meet practical needs. erefore, other accurate methods of the climate regionalization that comprehensively consider multiple climatic factors need to be developed. Principle clustering analyses (PCA) and clustering algorithms have been widely used to define different climate zones in processing meteorology and geography data [13].
e main impact factors can be selected from the complex climate impact factors determined by PCA and defined as principal components [14]. And the homogeneous regions can be identified using clustering techniques [15,16]. It is necessary to integrate the two methods described here into the climate regionalization of asphalt pavement. is comprehensive method can eliminate the correlation between evaluation indexes to obtain more accurate clustering results.
To this end, this study proposes a multi-index method for the climate regionalization of asphalt pavements based on the K-means clustering algorithm. is method takes into account various climatic factors and eliminates the correlation between those factors. e climate data of Liaoning Province from 2001 to 2014 are statistically analyzed. e effects of high temperatures, low temperatures, precipitation, and solar radiation are considered comprehensively. To verify the accuracy of the regionalization effect, we conducted a road degradation survey in the research region. e results of the climate regionalization were found to be consistent with the actual road conditions as established by a comparative analysis. In addition, a probabilistic neural network (PNN) prediction model and the support vector machine (SVM) prediction model were constructed and found to accurately predict the partition results. is method can be used to guide the selection of pavement materials and design of asphalt mixtures of construction.

Principal Component Analysis.
A PCA is a dimensionality reduction method based on a linear transformation. is method can convert a high-dimensional data set into low-dimensional space by some linear transformation. And then the problem to be solved is converted into one of determined covariance matrix eigenvalues and eigenvectors of the sample data. e goal of a PCA is to reflect the variation in the sample data with a data set characterized by lower dimensionality. It also completes the transformation of a set of potentially correlated variables into a set of linearly uncorrelated variables, thereby simplifying the calculations [17].
Each sample, X, has n-dimensional characteristics, X � (x 1 , x 2 , . . . , x n ) T in the sample space. At the same time, each feature meets the conditions x i ∈ R n and (1/n) n i�1 x i � 0. e covariance matrix of the sample is shown in where e eigenvalues of the covariance matrix, λ 1 ≥ λ 2 ≥ · · · λ n ≥ 0, and the corresponding eigenvectors, u 1 , u 2 , . . . , u n , are obtained according to the characteristic equation CU � λu.
e ith principal component is F i by projection. e calculation is shown in (2) is the cumulative contribution rate of F i . e greater the cumulative contribution rate, the smaller the body of information that is lost. When the cumulative contribution rate is more than 85%, the principal components selected contain the most information of the original variables [18,19]. erefore, we take the principal component with a cumulative contribution rate greater than 85% as the comprehensive index to best reflect the information of the sample data.

K-Means Clustering Algorithm.
e K-means clustering algorithm is an unsupervised learning method. In this method, the sum of the distances from the data points to the initial clustering center is taken as a clustering criterion function. e adjustment rule of the iterative operation is obtained by solving for the extreme values of the function. Similar data points are clustered into corresponding K clusters through repeated iterative operations. e optimal classification of the sample data is ultimately obtained [20].
Each sample meets the condition x (i) ∈ R n in the m-dimensional sample space X � x (1) , x (2) , . . . , x (m) . K-cluster centroid points μ 1 , μ 2 , . . . , μ K ∈ R n are randomly selected as the initial cluster centers. e squared errors of the algorithm and the criterion function are defined in where Γ ik is 1 when data point i is classified into the Kth cluster. Otherwise, Γ ik is 0.

Advances in Civil Engineering
As J c becomes smaller, the clustering process gradually reaches a steady state. e iterative operation ends when J c reaches a minimum value. In order to ensure that the true minimum value of J c is reached, an iterative operation is adopted and μ k is initially fixed to select the optimal Γ ik . As long as the data point is classified to the nearest cluster center, J c can be guaranteed to be the minimum. en, Γ ik is fixed to solve for the optimal μ k . It is a straightforward operation to obtain the minimum value of J c when μ k is the average of the data points in all clusters by differentiating J c with respect to μ k and setting the derivative equal to zero. e iterative operation formula is shown in

Probabilistic Neural Network.
e PNN is a parallel algorithm based on the Bayesian classification rules and probability density function estimation methods. PNNs can realize the functionality of a nonlinear learning algorithm by using a linear learning algorithm [21] and exhibits the high precision of nonlinear algorithms. e PNN algorithm includes an input layer, a pattern layer, a summation layer, and an output layer [22]. e structure of the PNN used in this study is shown in Figure 1. e input layer passes the input eigenvector to each node of the pattern layer through a linear function. is results in a scalar product with a weighted sum, as shown in where x is the input eigenvector and w i is a weighted coefficient corresponding to the eigenvector. e degree of matching between the input eigenvector and each pattern in training set is calculated in the pattern layer. And then the input eigenvector is sent to a Gaussian distance function to obtain the output result. e radial basis operation exp[(Z i − 1)/σ 2 ] is used as the activation function in the pattern layer. e input and output relations corresponding to the jth neuron of the ith mode are determined by where, i � 1, 2, . . . , M, M is the total number of training samples; j � 1, 2, . . . , N i , N i is the number in training samples of category i; d is the dimension of the sample space data; x ij is the jth center of the samples of in category i; and σ is the smoothing factor. e summation layer is responsible for connecting the neurons belonging to the same category in the pattern layer and finding the weighted average in accordance with where v i is the output of category i and N i is the number of neurons in category i. e output layer is responsible for producing the category with the highest weighted average in the summation layer, as given by

Support Vector
Machine. SVM is a generalized linear classifier used for the binary classification of data with good separation effects. e climate regionalization challenge in this work was to solve a multiclassification problem. erefore, it was necessary to combine multiple binary classifiers to construct a multiclassifier. e "one-on-one" method has a higher test accuracy when there are fewer categories [23]. erefore, this study utilized a "one-on-one" method to build a multivariate classification SVM. e calculation process is shown in Figure 2.
e "one-to-one" SVM makes use of a voting methodology. First, a decision boundary is established for any two of n samples, meaning that there are a total of C 2 n decision boundaries. Each decision boundary produces a classification result and a vote for a category. en, the number of votes for each type is determined. Finally, the category with the highest number of votes is defined as the classification of the sample. For the decision boundary consisting of the samples in category i (labeled "+1") and category j (labeled "−1"), the optimal decision function is found by using s.t.

Advances in Civil Engineering
where C is a penalty parameter that controls the degree of penalty applied for the wrong classification sample and x i is the input sample. w is a weight vector. ξ ij i is a slack variable that relaxes the hard boundary constraint.

Selection of Major Climate Variables.
To obtain accurate and reasonable results, we established climate regionalization by considering the effects of major climatic variables on the service performance of asphalt pavements. e major climate variables included high temperature, low temperature, precipitation level, and solar radiation level.

High Temperature.
High temperatures have a significant impact on the service life of asphalt pavements and are often the basis of determining the design of the pavement mixture and section. Asphalt becomes a plastomer that deforms under high-temperature conditions, and its shape cannot be completely recovered. Hence, asphalt pavements are prone to permanent deformation under the influence of repeated loading [24], and high temperatures are one of the main causes of the permanent deformation of asphalt pavements. High temperatures also directly affect the smoothness of the pavement and reduce the performance and service life of the road. ere is also a direct positive correlation between rut depth, asphalt penetration, and temperature [25]. In addition, high-temperature conditions have a significant effect on the thermal-oxidative aging of asphalt [26]. erefore, the influence of high temperatures should be considered when implementing climate regionalization. e maximum high temperature over a span of n years, the annual average maximum temperature, and the annual average temperature were used as the main high-temperature variables. e formulas for each of these parameters are shown in where T max,year is the maximum high temperature over a span of n years (°C), T max ,year,i is the maximum high temperature in year i (°C), T average, max ,year is the annual average maximum temperature (°C), T max ,month,i is the maximum temperature in month i (°C), T average,year is the annual average temperature (°C), and T average,month,i is the annual average temperature in month i (°C).

Low Temperature.
Asphalt pavements are prone to temperature stresses which often result in the formation of cracks in low-temperature environments [27]. Generally, asphalt mixtures have good low-temperature crack resistance [28]. But cracking tends to become more severe with repeated fatigue due to seasonal alternation and temperature cycling. A large area of the alligator cracking will often appear. ese cracks will shorten the service life of roads dramatically. is is one of the main forms of damage to asphalt pavement. Liaoning Province is a large province in China with frequent freeze-thaw cycles. e low-temperature degradation of asphalt pavements is particularly prominent. erefore, low-temperature factors must also be considered when implementing climate regionalization. e minimum low temperature over a span of n years and annual average minimum temperature were used as the main low-temperature variables for the calculation. e formulas for those parameters are shown in T average, min ,year � 1 12 where T min,year is the minimum low temperature over a span of n years (°C), T min,year,i is the minimum low temperature in year i (°C), T average, min ,year is the annual average low temperature (°C), and T max,month,i is the minimum low temperature in month i (°C).

Precipitation.
e infiltration of water is another cause of early damages to asphalt pavements. Water damage is characterized by pit slot and pavement pumping. Water first penetrates into the asphalt mixture. en, water will cause the adhesion between asphalt and aggregate to decrease. Water infiltration into the asphalt mixture reduces the adhesion of the asphalt to the aggregate. e damage is more severe especially under vehicle loads and hydrodynamic pressure.
e strength and durability of the mixture are reduced and the pavement is often destroyed [29][30][31]. In addition, snow melting cycles in cold regions will reduce the strength of the bond between the asphalt and skeleton structure. e stability of the pavement structure is reduced and it can be easily destroyed by the impact of continuous external loads [32,33]. erefore, the precipitation factor should be considered when implementing climate regionalization.
e maximum precipitation over a span of n years and the annual average precipitation were used as the main precipitation variables for the calculations. e formulas for those parameters are shown in where P max ,year is the maximum precipitation over a span of n years (mm), P average,year is the annual average precipitation (mm), P max ,year,i is the maximum precipitation in yeari (mm), and P average,month,i is the average precipitation in month i (mm).

Solar
Radiation. Asphalt pavements exposed to the atmosphere absorb a considerable amount of solar radiation. It directly leads to the increase of asphalt pavement temperature and the aging of pavement materials. e maximum temperature of asphalt pavements subjected to solar radiation can reach 60°C-70°C [34]. At this time, the antideformation ability of asphalt pavement is greatly reduced, and rut is likely to occur under vehicle loads [24]. Moreover, the intensive solar radiation will also accelerate the thermal oxidation and photooxidation of asphalt materials [35][36][37]. e long-term ultraviolet radiation causes serious degradation of asphalt performance and attenuates the durability of pavements [38]. erefore, solar radiation intensity must be considered when implementing climate regionalization. e maximum radiation over a span of n years and annual average radiation were used as the main radiation variables. e formulas are shown as where R max ,year is the maximum radiation for a span of n years (10 8 J/m 2 ), R average,year is the annual average radiation (10 8 J/m 2 ), P max ,year,i is the maximum radiation in yeari(10 8 J/ m 2 ), and P average,month,i is the average radiation in month i (10 8 J/m 2 ).

Sample Clustering.
A PCA was implemented via a method based on the eigenvalue decomposition of the covariance matrix. Many correlated climate variables were transformed into a few uncorrelated variables and nine major climate variables from 56 weather stations were used as the input data. First, each major climate variable was subtracted from its respective group average to meet the initial conditions of the PCA. e covariance matrix was calculated according to (1). en, the eigenvalues and eigenvectors of the covariance matrix were calculated using the SVD function and the eigenvalues were sorted from large to small. e eigenvalues in the Kth cluster and the maximum range were chosen. e appropriate principal components were then selected as the clustering factors using (2) and the condition that the eigenvalue was greater than one. Finally, the clustering factors determined by the PCA were taken as the input data space. e seed point of the initial cluster was determined by the sequential method. ese initial clusters seed points were analyzed iteratively according to (3). e optimal classification of the sample data was determined when the iteration caused J c to reach a minimum value.

Research Region.
e Liaoning Province, located between 118°53′ and 125°46′ east longitude and between 38°43′ and 43°26′ north latitude, belongs to the North Temperate Zone. e northern latitude of 42°experiences a high percentage of annual solar radiation and runs through the central part of Liaoning Province. Liaoning Province is the region with the strongest radiation, the most illumination, the highest temperature, and the shortest frost period in the Northeast. It is located on the east coast of the Eurasian continent. e coastal location results in a temperate continental monsoon climate. Precipitation is concentrated in the summer with relatively constant and abundant rainfall.
Advances in Civil Engineering e spatial distribution of the precipitation decreases from southeast to northwest.

Primary Types and Causes of Degradation of the Main Roads.
A representative road section in Liaoning Province was selected for an on-site investigation. Core samples were taken to determine the main causes and locations of typical degraded conditions. e primary types and proportions of the pavement degradations in Liaoning Province are shown in Figure 3.
It can be seen from Figure 3 that cracking of the pavement is common in the main roads in Liaoning Province. Alligator cracking, longitudinal cracking, and transverse cracking account for 84.24% of the total degradation area. Alligator cracking is the most serious, accounting for 32.68%, and is mainly caused by the fatigue of the asphalt layer or inorganic binder layer. Longitudinal cracking accounts for 26.62% and is mainly caused by insufficient strength of the base layer and vehicle loads shearing.
e amount of transverse cracking is roughly equivalent to the amount of longitudinal cracking and accounts for 24.94% of the degradation area. Transverse cracking is mainly caused by the shrinkage of the asphalt mixture under dry or low-temperature conditions. e cracking of the pavement provides a basis for water damages, which can easily cause deterioration, and accounts for 2.18% of the degradation area. Of the types of degradation observed at high temperatures, rutting is predominant and accounts for 6.13% of the degradation area. Rutting is mainly caused by a combination of vehicle loads and environmental factors. e remaining types of degradations are related to construction factors and account for 7.45% of the degradation area. Climatic conditions are the most important factors and leading causes of these conditions. Overall, vehicle loads and climatic conditions are the main factors that directly affect the performance of asphalt pavements.

Major Climate Variables.
According to meteorological data from 2001 to 2014 in Liaoning Province, the major climate variables were determined to be as follows.
e high-temperature variables were calculated with equations (10)- (12), and the statistical results are shown in Figure 4. e annual average temperature in Liaoning Province ranged from about 6°C to 13°C. e maximum temperature was in the range of 33.4°C to 41.8°C, with the maximum temperature of 41.8°C in Jinzhou. e low-temperature variables are calculated with equations (13) and (14), and the statistical results are shown in Figure 5. e minimum temperature in winter is in the range −41.1°C to −19°C and is −41.1°C in Xifeng. e precipitation variables were calculated with equations (15) and (16), and the statistical results are shown in Figure 6. e average annual precipitation in Liaoning Province reached a maximum value of 1023.4 mm. e upper limit of the annual maximum precipitation over a span of n years was 1765.5 mm. e solar radiation variables were calculated with equations (17) and (18), and the statistical results are shown in Figure 7. e annual average amount of radiation was within the range of 4.5 × 10 8 J/m 2 to 5.5 × 10 8 J/m 2 . e upper limit of the solar radiation over a span of n years was within the range of 48 × 10 8 J/m 2 to 57 × 10 8 J/m 2 , which is more far than the annual average.

Cluster Factor.
e common factor between the nine major climate variables was extracted by means of the main axis method, as shown in Table 1. e second column contains the eigenvalue of each factor. e third column lists the variance contribution rate for each factor. e fourth column is the cumulative variance contribution rate of all factors. e larger the three values in those columns are, the more variable information the factor contains. A cumulative  variance contribution rate of 85% meets the general requirements of this study. Factors with eigenvalues greater than one should be retained and three factors were retained according to this criterion. e cumulative variance contribution rate of these factors reached 88.415%, exceeding the general requirement of 85%. A steep slope test of the sample factors was performed in MATLAB. e eigenvalues of each factor were ranked from largest to smallest and resulted in a slope according to the theory of PCA. If a factor reduces the eigenvalue of the slope, the eigenvalue of this factor should be defined as a threshold. e factors with eigenvalues at or above the threshold should be retained. In this study, the first three factors were retained, as shown in Figure 8. e three retained factors were analyzed by the maximum variation method of orthogonal rotation, as shown in Table 2. e load value represents the relationship between each factor and the principal component. e greater the load value, the closer the correlation between the factor and The maximum precipitation for a span of n years Average annual precipitation    The minimum low temperature for a span of n years Average annual minimum temperature   Advances in Civil Engineering principal component. e extreme maximum precipitation and the annual average precipitation had higher load values for the first principal component. Each temperature variable had a higher load value for the second principal component, and the annual average radiation and extreme maximum radiation over many years had higher load values for the third principal component. e precipitation factors, temperature factors, and radiation factors were able to explain the original climate variables comprehensively. erefore, the first, the second, and the third principal components were defined as the precipitation factor, temperature factor, and radiation factor, respectively. e three principal components were then used as clustering factors for the climate regionalization.
All retained factor variables must be uncorrelated, or orthogonal, according to PCA theory. ere is no correlation between the factors retained in this study, because the PCA eliminated the correlation between the major climate variables. e covariance matrix of the principal component factors is shown in Table 3. e scores of the three principal component factors were used as the scores of each clustering factor in the climate regionalization. e clustering factor scores are shown in Figure 9.

Climate Regionalization.
e 56 meteorological stations in the research area were divided into four station groups according to the clustering factors obtained with MATLAB. e clustering results are shown in Figure 10.
Finally, ArcGIS was used to delineate the four climate zones resulting from consideration of the asphalt pavement environment, together with regional administrative boundaries. e specific regional types are shown in Figure 11. e results of the climate regionalization are as follows.
Zone I: this zone is significantly affected by high temperatures, low temperatures, and precipitation. is zone includes Fushun, Tieling, Faku, Kangping, and the central and northern parts of Liaoning Province. e extreme minimum temperature reached −43.4°C in the winter. e key in this zone is to improve the lowtemperature crack resistance of asphalt pavements and prevent the occurrence of reflective cracks. An asphalt mixture with a small void ratio, high strength, and good resistance to low-temperature cracking should be used. It is also necessary to ensure stability of the asphalt pavement regarding exposure to high temperatures and water. Zone II: this inland zone is characterized by high temperatures, little rainfall, and high radiation levels.
is zone includes Yuzhong, Fuxin, Yixian, Yangshan, and Kazou in the northwest. e extreme maximum temperature reached 41.3°C and the key in this area is to focus on ensuring the high-temperature stability of the asphalt pavement. erefore, asphalt mixture grading should be strictly controlled. A skeleton compact    structure should be adopted. Asphalt modified or high modulus additives should be used to improve the rut resistance of the road. Low-temperature crack resistance and high-temperature rut resistance must also be taken into account when designing asphalt mixtures in this zone. Zone III: this coastal zone is characterized by high temperatures and abundant rainfall. It includes Dalian, Yingkou, Panjin, and Gaizhou in the southwest, and the extreme maximum temperature reached 39.3°C. e extreme maximum rainfall was 1132.5 mm. Asphalt pavements in this zone should be guaranteed to have sufficient stability under high temperature and abundant rainfall conditions. At the same time, the pavement must have good antislip performance in view of the large amount of rainfall. erefore, it is advisable to use stone mastic asphalts (SMA) and modified asphalts as much as possible. Zone IV: this is a rainy and mountainous zone and includes Zhuanghe, Yuyan, Caohekou, Benxi, and Xinbin in the east. e area is subjected to considerable rainfall, with a maximum precipitation of 1765.5 mm.
Asphalt mixtures in this zone should be SMA, if necessary. e base material should have good erosion and corrosion resistance, and a dense framework structure should be adopted.

Probabilistic Neural Network Classification Prediction.
A PNN was established in this study through the use of the "newpnn" function in MATLAB. e input layer included the precipitation factor, temperature factor, and radiation factor. e output layer provided the four types of climate zones. After selecting data sets 1 to 36 as a training sample, the expected category was converted to a vector. Similarly, data sets 37 to 56 were used as test samples. e classification effect and error of the training data after training of the PNN are shown in Figure 12. e training samples were input into the established PNN network. e results showed that there was no error in the process of judging the sample data. erefore, the PNN can perform additional sample classification predictions. e classification of the predicted data after the test is shown in Figure 13. Only two sets of data were handled incorrectly when 20 sets of samples were predicted. erefore, the accuracy rate of the classification prediction was 90%.

Support Vector Machine Classification Prediction.
In this paper, an SVM climate regionalization model was established using the LIBSVM toolbox in the MATLAB software.
irty-six sets of data were used as training sets and twenty sets of data were used as test sets, similarly to the PNN classification. In order to find the best parameter pair, the SVM parameters were selected by the grid search method. e generalization ability of each set of parameter pairs (c, g) on the grid was evaluated. A cross-validation method called coefficient of variation was used to assess generalization capabilities. A three-dimensional view of the predicted accuracy was obtained through the LIBSVM toolbox and is shown in Figure 14.
e prediction model has the highest classification accuracy when the model parameters were c � 1.4 and g � 4.0, in which case the accuracy reached 94.4%. erefore, the optimal parameters of the model were determined by the cross-validation parameter search. en, the SVM classification model was established. e classification prediction is shown in Figure 15. e similarity between the prediction results of the SVM and K-means clustering analysis of the climate regionalization was 95% which indicates that the K-means clustering analysis results exhibited a reasonable degree of accuracy and reliability.

Verification of Physical Engineering.
To further verify the reliability of the proposed regionalization method, this study compared the degradation of the pavement in the research regions with the results of the climate regionalization. Transverse cracks and longitudinal cracks were common in all studied roads. However, the situation in Zone I was the most serious, with cracking frequencies that reached 13.2% and 20.7% for the transverse and longitudinal cracks, respectively. ose same two conditions occurred with similar frequency in the other zones and were not representative; therefore, transverse and longitudinal cracks were not considered in other zones. e percentage of the damaged area on each road is shown in Figure 16.
A regionalization plan was generated for each road according to the method proposed in this study:  Figure 11: e specific regional types.
Advances in Civil Engineering cracks. ese types of cracks occur when water 3D view of parameter selection infiltrates along the cracks during a precipitation event, where it readily causes the loss of cementitious materials. In addition, water in the aggregate can cause peeling under the influence of vehicle loads and freeze-thaw cycles. Longitudinal cracks penetrate the base layer, which is weakly connected to the surface layer. ese cracks are mainly the result of uneven settlement or frost heaving of the roadbed. Alligator cracks and rutting are the result of aging of the asphalt surface layer and fatigue of the inorganic binder layer. Settlement is mainly caused by insufficient strength of the soil base when the base layer is eroded by water. Overall, the pavement is affected by high temperatures, low temperatures, and precipitation in Zone I. Zone II includes five routes: Jinchi, Jingshen, Xiaoye, Anyang, and Binhai Road. e pavement degradation consists mainly of alligator cracks, rutting, shoving, and bleeding. e causes of cracks, rutting, and shoving in this zone are similar to those in Zone I and mainly arise as a result of insufficient overall road strength and the combined effect of vehicle loads and high temperatures. Overall, the pavement is most affected by high temperatures in Zone II. Zone III includes a total of four routes: Haipi, Kuer, Zhuanglin, and Heida (Dalian). e pavement problems are mainly rutting and bleeding in this zone. e rut deformation occurs primarily in the asphalt mixture layer. Judging from its deformation mode, the asphalt mixture undergoes lateral displacement and produces ruts under the repeated influence of high-temperature conditions and vehicle loads. Overall, the pavement is most primarily affected by high temperatures in Zone III. Zone IV includes a total of six routes: Xiaoxiao, Heida (Liaoyang), Danji, Huanghai (Dandong), Huanghai (Anshan), and Benhuan. e pavement problems are mainly the result of alligator cracks, rutting, and bleeding. e road surfaces in Zone IV do not typically exhibit water-induced damage such as deterioration or pit slot. However, the cracks provide an opportunity for water damage. e cracks develop from top to bottom until they finally penetrate the surface layer. is results in a poor connection between the top layer and base layer. e core sample taken of the surface layer was significantly fragmented and had a low cement content and low strength.
e process whereby these cracks develop is accelerated when water infiltrates the crack during a precipitation event and then remains in the base. In these circumstances, the bottom of the surface layer and the top of the base layer at the crack site are saturated for long periods of time. Vehicle loads generate hydrodynamic pressure, and temperature changes produce freeze-thaw cycles. e fine aggregate and cementitious materials peel off, and the adhesion between the surface layer and base layer is diminished. As the water has not penetrated the full depth, the base layer maintains good strength and a high cement content. Overall, precipitation is the main cause of road degradation in Zone IV.
In summary, the causes of pavement degradation in each road are consistent with the results of regionalization. erefore, the method proposed in this study exhibits considerable reliability.

Conclusions
Based on the study of the climate regionalization of asphalt pavement, the following conclusions can be obtained: (1) e four climate factors of high temperature, low temperature, precipitation, and solar radiation are analyzed in detail. ese factors are the main causes of pavement damage. Nine major climate variables related to the four climate factors were proposed and a PCA was used to extract the three principal component factors of precipitation, temperature, and radiation from the main climate variables. (2) e asphalt pavement in Liaoning Province was divided into four climate zones based on the K-means clustering analysis method. Zone I is consistently affected by high temperatures, low temperatures, and precipitation; Zone II is an inland zone with high temperatures, little rainfall, and radiation; Zone III is a coastal zone with high temperatures; Zone IV is a rainy, mountainous zone.    (3) e accuracy of the climate regionalization results was evaluated by the developed PNN and SVM prediction models. e results showed that the accuracy of the regionalization obtained from the K-means clustering algorithm was greater than 90%. us, this method has clear accuracy and reliability. (4) e method for the climate regionalization of asphalt pavements only analyzed the meteorological data of Liaoning Province in this paper. However, this method can serve as a reference and guidance for the design and construction of highway projects with similar climatic concerns around the world.
Data Availability e data in this study can be obtained from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.