Next Article in Journal
2-D Selector Simulation Studies on Grain Selection for Single Crystal Superalloy of CM247LC
Next Article in Special Issue
Mechanical Characterization of Gypsum Composites Containing Inert and Insulation Materials from Construction and Demolition Waste and Further Application as A Gypsum Block
Previous Article in Journal
Heat Treatment, Impact Properties, and Fracture Behaviour of Ti-6Al-4V Alloy Produced by Powder Compact Extrusion
Previous Article in Special Issue
Performance and Durability of Porous Asphalt Mixtures Manufactured Exclusively with Electric Steel Slags
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Temperature and Moisture Content on Pavement Bearing Capacity with Improved Subgrade

by
Audrius Vaitkus
1,
Laura Žalimienė
1,*,
Jurgita Židanavičiūtė
1 and
Daiva Žilionienė
2
1
Road Research Institute, Vilnius Gediminas Technical University, Linkmenų g. 28, 08217 Vilnius, Lithuania
2
Department of Roads, Vilnius Gediminas Technical University, Saulėtekio al. 11, 10223 Vilnius, Lithuania
*
Author to whom correspondence should be addressed.
Materials 2019, 12(23), 3826; https://doi.org/10.3390/ma12233826
Submission received: 31 October 2019 / Revised: 15 November 2019 / Accepted: 18 November 2019 / Published: 21 November 2019
(This article belongs to the Special Issue Recent Research in the Design of New Sustainable Building Materials)

Abstract

:
Environmental conditions (temperature, moisture and the intensity of the sun) influence variation in asphalt pavement strength during the year. Lithuania is situated in a zone by average warm summers and average cold winters, and the most important climatic factor is the variation of the air temperature. This study presents the influence of temperature (of asphalt concrete (AC) and subgrade layers) and moisture content (of subgrade layers) to the pavement bearing capacity. The experimental research was obtained in five pavement sections of the experimental road. This experimental road was constructed in 2007 in Lithuania and is operated for more than 12 years. This paper presents a statistical analysis between the bearing capacity and the thickness of the asphalt concrete layers, the temperature and moisture content of different pavement layers, among sections, loaded and unloaded lanes (right and left wheel paths and tracks). The bearing capacity was evaluated by a falling weight deflectometer (FWD), temperature and moisture content by electronic sensors and thickness of AC layers by Georadar. Analysis of overall E0 and E0,h9 (bearing capacity at a depth of 9 cm from asphalt surface) data declares that seasonal impact on pavement structural strength due to a change of subgrade bearing capacity remains after correction of asphalt stiffness dependent on the temperature in the layer. However, it was detected that neither E0 nor E0,h9 are related to moisture content at a depths of 100 cm, 130 cm and 150 cm.

1. Introduction

Designing pavements according to the present design method assuring and fulfilling construction conditions provides a strength pavement structure (PS) for a given traffic load during a design period under a given climate condition [1].
In addition, Fonod [1] mainly characterizes the serviceability of the pavement:
  • The roughness of the pavement affects driving comfort;
  • Skid resistance affects safe driving;
  • Bearing capacity of the road pavement affects road pavement service life;
  • Structural condition (cracking and patching) affects driving comfort in addition to aesthetic effects.
The bearing capacity of pavement is defined as a characteristic of its structural state, related to the capacity to support the expected traffic. The evaluation of this capacity can be made through parameters such as residual life, required thickness for reinforcement, stiffness modulus and deflection [2].
According to Straube and Jansen [3], the bearing capacity of AC pavements is evaluated by a FWD (the equivalent elastic modulus—E0). Therefore, the elastic modulus of the AC layers is influenced by the load frequency and AC layer temperature. Whereas, the elastic modulus of the subgrade is influenced by the current water content of the unbound layers and the temperature of the AC layer in case of air temperature below 0 °C.
By Salour [4], solar radiation of the AC surfacing and geothermal heat flux contribute to pavement heat intake, and surface convection and radiation emissions contribute to heat extraction [5].
Alawi and Helal [6] results show that the maximum pavement temperature is not on the surface, but rather at a full depth from the surface. This type of temperature change explained due to the heat absorption by the asphalt concrete inside the pavement, and due to air temperature and wind action on the pavement surface. Pavement temperatures remain in the range of 35–50 °C for most of the year. They noticed that temperatures (>68 °C) are attained at a depth of 20 mm and then the temperatures decrease in each layer to its normal distribution. Therefore, the thickness of the asphaltic pavement must be taken into consideration, increasing indirect tensile stiffness modulus for the upper layers of the asphaltic pavement and using materials less affected by high temperature.
According to Motiejūnas et al. [7] based on the analysis of experimental investigation results it is recommended to measure pavement deflections by the falling weight deflectometer when the temperature of asphalt pavement layers varies within the interval +5 °C to +25 °C, except cases when special investigations are carried out during winter freeze and spring thaw periods. Under higher temperature +30 °C to +35 °C the layers of asphalt pavement lose part of their smoothness.
García and Castro [8] studied the effect of temperature on the pavement deflection based on seven deflection tests, at different temperatures, in a section of road with the flexible pavement. The statistical analysis showed a good agreement between the experimental and theoretical obtained factors. However, the influence of temperatures on pavement deflections through experimental factors (with flexible pavements composed of 250 mm of asphalt mixes in a test that representative asphalt temperature range between 15° and 25 °C) is slightly greater in comparison with theoretical factors.
Zheng et al.’s [9] study showed that there is a particular relationship between the pavement thickness and the temperature correction coefficient of the asphalt pavement deflection. With the increase in asphalt pavement thickness, the correction coefficient change rate speeds up. The correction coefficient of 70 mm asphalt pavement thickness is 1.11 at a low temperature of 0 °C and 0.71 at high temperature of 50 °C, while the values are 1.29 and 0.40, respectively for 250 mm asphalt pavement thickness, which indicates that high attention should be given to the temperature correction of asphalt pavement thickness deflection. With the increase in asphalt pavement thickness, the changing trend among temperature correction coefficients becomes stable, for example, the difference of temperature correction coefficients between 150 and 250 mm asphalt pavement thickness deflection is minimal.
As surfacing temperature depends significantly on the weather, typically changing hourly and daily, the physical process never reaches a steady-state, i.e., there would be only a long-term equilibrium [10].
It is known that high temperatures of pavement allow the development of rutting and structural damage during hot and cold periods [11,12].
In seasonal frost conditions, load restrictions are commonly imposed during the spring-thaw period to prevent severe pavement deterioration. During this period, the PSs are usually exposed to excess moisture content, which results in reduced bearing capacity of unbound layers and high resilient and permanent deformations [4].
Precipitation (rain) is the primary source of water in pavements. Rainwater enters the road PS through cracks, infiltration from shoulders and drainage channel, high groundwater. The water gets inside the structure by an energy gradient, such as gravity, capillary forces, osmotic forces and temperature or pressure differences.
Moisture content can alter the resistance of permanent deformation in unbound granular materials dramatically [13,14]. The repeated load may cause positive pore water pressure to increase in unbound granular layers, which will consequently reduce the effective stress in unbound granular layers, which in turn decreases the stiffness and resistance to permanent deformation [15].
The permeability of the material that constitutes a road increases typically from the top of the pavement downward until about 0.7 m depth. It is estimated that about 80% of the problems encountered in pavements is related to the presence of water [16]. Otherwise, water would accumulate onto the low permeability layer and keep the upper layer wet; freezing of the accumulated water might then unbind the upper layer. Water would decrease the AC pavement bearing capacity and service life [17,18].
There is some evidence to suggest that water has less impact on the ticker and well-constructed pavements than it does on thinner ones [19]. It appears that in thicker pavements, the effect of water may be more indirect than in thinner ones, reducing material stiffness leading to later distress [17].
Charlie et. al [20] modeled the response of pavement at different water contents and found that an increase of moisture content from 2.3% to 4.8% led to an increase of the strains by about 60%. Chen, Sun and Yao [21] detected that base’s resilient modulus reduced significantly (approximately 30%) with each 1% increase in moisture content.
According to Salour and Erlingsson [22], collected data showed evident correlation between the moisture content and the back-calculated unbound layer stiffness. Back-calculation of the FWD deflection data revealed 63% reduction in subgrade and 48% reduction in granular layer during the spring thaw period compared to the recovered summer values. They observed a rapid increase in basin deflection during the thawing period, with a 4.7 times increase in maximum deflection value from a fully frozen structure to a thawed base and subbase layer over a 60 cm frozen subgrade.
This study aims to evaluate how temperature (of AC and subgrade layers) and moisture content (of subgrade layers) influence pavement bearing capacity. To achieve this goal, five pavement sections were selected of the road of experimental pavement structures (furthermore the Test Road) that was constructed in 2007. In this way, it was assessed the difference of the bearing capacity of the road PS among sections, loaded and unloaded lanes (right and left wheel paths and tracks). Further, the statistical analysis was done between these factors: the thickness of the AC layers, the temperature on the top of the AC surfacing and in the AC layers, and the temperature and moisture content of the soil in the subgrade and the bearing capacity. Finally, conclusions according to the results of the experimental study are represented at the end of the paper.

2. Main Characteristics of the Test Road

The site in Pagiriai (about 20 km from Vilnius the capital of Lithuania) was selected for constructing the Test Road. This place fulfilled all the conditions required for the experiment: it has a sufficient traffic volume of heavy vehicles, located in open terrain, has no horizontal plane curves or vertical curves in longitudinal section and could be distinguished by the same irrigation conditions within the route of the road section. The Test Road was performed on the road to the gravel quarry where the right lane is used by loaded traffic (herein loaded lane) and the left lane by unloaded traffic (herein unloaded lane; Figure 1). The total length of the Test Road is 710 m and consists of 23 sections of 30 m, one section of 20 m and three sections of 15 m in length (in total 27). The width of the Test Road carriageway is 7 m, i.e., two lanes of 3.0 m, and two shoulders of 0.5 m. The road pavement structures (PSs) are designed for 100 kN axle load ESAL100 = (0.8–3.0) million by Construction Technical Regulation STR 2.06.03:2001 Automobile Roads [23,24,25,26,27]. It has all necessary electronic sensors (loop profilers, temperature and moisture sensors and stress and strain transducers). The bearing capacity of subgrade was improved >100 MPa (Ev2) during the construction period in each section (the bearing capacity of subgrade 459 MPa (PS No 4), 205 MPa (No 12), 150 MPa (No 18), 177 MPa (No 24) and 137 MPa (No 25)). The paper presented experimental researches in PS No 4, No 12, No 18 and No 24–25 (Table 1) because these sections are equipped with temperature and humidity sensors. The numbering of structures in this paper was accepted by the serial number of the Test Road sections.
The temperature sensors and moisture content sensors were installed in the soil of subgrade at a different depth (Table 1), and measurements were carried out continuously throughout the study period (every hour) in each section [25].
Every year in the sections, these measurements were made:
  • Traffic monitoring;
  • Temperature and moisture content at a different depths of the road PS;
  • Roughness;
  • Rutting;
  • Longitudinal gradient and cross fall;
  • Visual assessment of distress on the road PS;
  • Road pavement deflection (measured with a FWD, Dynatest, Vinci, Italy and Benkelman Beam, Infratest, Brackenheim, Germany);
  • Skid resistance (measured with a Pendulum device).
Five road PSs with various materials were selected for the study (Figure 1). AC layers of these sections were constructed from AC11WL, AC16BnL and AC32BsL. Herein AC11WL was an asphalt concrete wearing layer, with a mesh size of top sieve of 11 mm and affected by heavy loads; AC16BnL was an asphalt concrete binder layer, with a mesh size of top sieve of 16 mm and affected by heavy loads and AC32BsL was an asphalt concrete base layer, with a mesh size of top sieve of 16 mm and affected by heavy loads. Except one section, No 12, where surfacing performed from the SMA11 PMB WL (stone mastic asphalt wearing layer, with a mesh size of the top sieve of 11 mm, with polymer-modified bitumen and affected by heavy loads). The base course in all sections performed from crushed dolomite 0/56 and frost-blanket layer performed from sand 0/11. The different type of geosynthetic materials installed in sections No 24 and No 25.

3. Test Plan and Methods

Sixty measurements of the deflection were obtained in each section at some time. These measurements were performed after ten measurements on the left wheel path and right wheel path, including the track and repeated measurements on the loaded and unloaded lanes (Figure 2). The measurement repeated at eight times (Table 2). In total, 1920 E0 values were obtained.
This paper presented measurements data in sections and results of a statistical analysis based on linear mixed models (LMMs). These indicators were measured and used in statistical analysis:
  • The thickness of the asphalt concrete layer (measured by Georadar in 2015, Figures 4 and 5);
  • The temperature on the top of the road PS (Table 1, Figure 6);
  • The temperature at a different depth in the AC layers of the road PS (Table 1, Figure 6);
  • Temperature and moisture content in the soil of subgrade at a different depth (Table 1, Figures 6 and 7);
  • Surface modulus (E0, MPa) of the road PS (Figure 8);
  • Periods (Table 2).

4. Statistical Analysis and Evaluation of Results

The statistical analysis (Figure 3) consisted of two parts:
  • The first part was descriptive and contained some statistical characteristics of the factors under study (mean, median, minimum and maximum values, standard deviation and quartile) in road PS (furthermore RPS) No 4, No 12, No 18 and No 24–25. The data was visualized in various layers in box-plot diagrams;
  • In the second part, applying LMMs [28], it was attempted to assess whether the bearing capacity of the road PS has a statistically significant difference among sections and loaded and unloaded lanes (right and left wheel paths and tracks). The effect to the bearing capacity of the road PS of these factors was analyzed: the thickness of the AC layers, the temperature on the top of the AC surfacing and in the AC layer and the temperature and moisture content of the soil in the subgrade.

4.1. Descriptive Statistics

On all sections, the average asphalt concrete layer was 168.80 ± 7.57 mm. A slightly bigger thickness was recorded in the road PS No 12, it was 171.90 ± 8.66 mm, while the smaller thickness was distinguished by the road PS No 4 and No 18, with a thickness of the asphalt concrete layer 167.30 ± 7.35 mm and 166.60 ± 5.91 mm, respectively (Figure 4 and Figure 5).
When comparing the thickness of the AC layer on loaded and unloaded lanes, a thicker AC layer was observed on the loaded line than in the unloaded lane, respectively 171.80 ± 7.83 mm and 165.80 ± 5.95 mm (Figure 4). When comparing the thickness of the AC layer in different lanes (Figure 5), it was seen that the right wheel path was distinguished by a wider or more varied range of AC layer thickness values and these values continued to decrease towards the left wheel path. The average thickness of the AC layer in the right wheel path was 168.30 ± 8.46 mm, in the track was 170.50 ± 7.89 mm and in the left wheel path 167.70 ± 5.82 mm. The analysis of the AC layer thickness on lanes of individual sections (Figure 5) shows that PS No 4 and No 12 of the unloaded lane had the opposite tendency, i.e., the thickness of the AC layer was observed thinner in the right wheel path of PS No 4, and the thinner AC layer was observed on the left wheel path in the PS No 12.
Figure 6 presents the data of temperature on the top of the road PS and in the AC layer at a depth of 9 cm. Results of measured temperature at a depth 8 cm and 10 cm in the AC layer did not show the results because the temperature at a depth of 8 cm, 9 cm and 10 cm varied very slightly (based on the strong correlation (above 0.9)). Therefore, the temperature at a depth of 9 cm was used to recalculate E0 values. Figure 6 presents data of temperature, and Figure 7 presents data of moisture content of soil in the subgrade of road PS in sections, respectively. On all sections, the average temperature of asphalt concrete layers was 14.40 ± 7.9 °C; the average temperature of sub-grade was 14.40 ± 5.9 °C. The average moisture content of sub-grade was 8.86% ± 4.52%.
Surface modulus E0 was calculated according to the formula that was presented together with the earlier results by Motiejūnas et al. [7]. E0 was calculated by a measured deflection that was reduced to 50 kN load and to a standard temperature of +20 °C. Bearing capacity of the asphalt layers at depths of 8, 9 and 10 cm (herein E0,h8, E0,h9 and E0,h10) were calculated by Equation (1):
E 0 , h x = E 0 10 0.000221 · h a s f 1.0229 · ( T 20 )
where hx—depth of asphalt layer used for calculation of E0, cm; hasf—thickness of asphalt layers, cm and T—a temperature of the asphalt layers measured at a certain depth, °C.
When taking measurements from all sections, the mean of E0 of the road PS was 984.40 ± 100.91 MPa, and the means of the recalculated values E0,h8, E0,h9 and E0,h10 were approximately the same and equal to 933.00 ± 77.87 MPa, 926.80 ± 78.34 MPa and 935.10 ± 77.92 MPa, respectively (Figure 8).
There was no significant difference between the values of E0,h8, E0,h9 and E0,h10 because of the strong correlation (above 0.9) between the temperatures measured in AC layers at different depths (top of the road PS, 8 cm, 9 cm and 10 cm (Table 3)) and used for E0 values recalculation. Given that, only E0 and E0,h9 values were used in the analysis below.
Table 3 shows a matrix of correlations to evaluate how the factors used for analysis correlate not only with E0 but also how they relate to each other. In this correlation graph, the intensity of the color corresponds to the correlation strength, i.e., the darker blue color shows a stronger positive correlation and the darker red color shows a stronger negative correlation. Table 3 shows that E0 has a similar negative correlation (about –0.5) with all temperatures. A negative correlation means that E0 is decreasing with increasing temperature.
The comparison of E0 and E0,h9 values of the road PS in the sections (Figure 9) show that E0 and the recalculated E0,h9 values in section No 12 were on average 3−4% lower than in other sections.
The compared E0 modulus average value in the loaded or unloaded traffic lane (Figure 10 and Figure 11) changed slightly on average 5% (on average 932–976 MPa of the loaded lane and 990–1034 MPa of the unloaded lane). It is observed that the bearing capacity of the loaded lane was lower on average up to 15% (967 MPa of the loaded lane and 1108 MPa of the unloaded lane) in a loaded lane than in an unloaded lane (Figure 10).
When analyzing loaded and unloaded lanes of the road PS in the sections separately (Figure 11), it was seen that the road PS No 12 on the loaded lane was also characterized by a relatively lower E0 value than the other road PS on the loaded lane. However, this difference was not so high on the unloaded lane and not only the road PS No 12 but also the road PS No 18 on the unloaded lane were distinguished by somewhat lower E0 values than other sections, and the values of bearing capacity of these two sections No 12 and No 18 on the unloaded lane were very similar.
Additionally, E0 and E0,h9 values by the lanes (Figure 11) show that the right wheel path on the unloaded lane of all tested PS is distinguished by a lower E0, which is gradually increasing to the track and the left wheel path. The E0 value of the right wheel path in the road PS No 12 on the unloaded lane is the lowest compared to the other road PS. On the loaded lane, due to the lower loading all tracks of the road PS are distinguished by the lower E0, but there is no clear tendency between the left and right wheel paths on the loaded lane. The road PS No 12 again showed a significant tendency on the loaded lane. In this road PS, E0 was much lower on the loaded left wheel path than on the loaded right wheel path. There was a reverse tendency on this road PS on the unloaded lane, i.e., the lower E0 was observed on the unloaded right wheel path than on the unloaded left wheel path. Thus, it is to be seen that the road PS No 12 distinguished by its characteristics from another road PS.
During the exploitation of the road, negative factors of load are found in the path of the motor vehicle wheel. These factors influence the accelerated degradation of the road surface and the lower bearing capacity of the road PS compared to the unloaded lane, i.e., the track. The lower bearing capacity is particularly noticeable on the right wheel path of the loaded lane (up to 10%, the average E0 was 911 MPa at the right wheel path and the average E0 was 984 MPa at the left wheel path), because the right wheel path of the cross fall of the road PS was more loaded. By it, a recommendation is to install the crashed dolomite course up to the edge of the slope, thereby to strengthen with the unbound mineral mixtures the area of the shoulder.
The bearing capacity of the road PS was mainly influenced by the bearing capacity and sensitivity to the temperature and moisture content of the soils in the subgrade, the overall thickness of the road PS and the overall thickness of the AC layers with bituminous binders (Table 3). It can also be concluded that the use of fine-grained sand 0/4 for the frost-blanket layer (road PS No 4, the average E0 was 984 MPa) did not affect the bearing capacity and the temperature and moisture content of the soils in the subgrade compared to the road PS in which sand 0/11 was used (PS No 12, No 18, No 24, and No 25; the average E0 was 931 MPa). A difference of bearing capacity was about 6%. The bearing capacity of the insignificant impact is the type of unbound materials of the base course and the type of mineral mixtures with bituminous binders of base course and surfacing.
The variation of E0 and E0,h9 values were analyzed in time at each point of measurements separately in each section, i.e., in each section were 60 points where E0 was measured (Figure 12). Graphs show that the values of E0 and E0,h9 in all 60 points in each section had similar dynamics in time. The E0 was characterized by seasonal change in time, while the E0,h9 was characterized by a decrease in time, which was more noticeable at the last moments of measurement. However, these graphs also show a tendency that both E0 and E0,h9 values were different enough at different points in the same section. This tendency indicates that the strength of the road PS varies from one point to another.
In the second part of the statistical analysis, it was assessed what has the most significant influence on the distribution of E0 and E0,h9 values among different points of measurements, distinguishing the main factors (load, lane, right wheel path, left wheel path and track). However, the difference among points is still to be noticed, even under identical conditions. Consequently, besides the known conditions, there are specific, not observed and not included in the statistical analysis, factors, which influence the differences of E0 at different points of measurements.
Figure 12 shows that the decrease of v was observed in all the last measurements of all sections, and in section No 12, this decrease was sharper than in others. This difference of section No 12 could be influenced by a lower factor of the bearing capacity, as it was already noticed above that the average E0 of section No 12 was lower than of the other investigated sections. It can also be assumed that the significantly lower E0 values of the last measurement (April 2015) indicate that the bearing capacity of the road PS was influenced by season (hydrothermal mode). This measurement period was just a reflection of spring thaw and increased moisture content in the subgrade, but such a period did not occur during the three years of research. In this regard, it is necessary to carry out additional measurements of the bearing capacity of the hydrothermal mode in March, April and May.
Figure 13 shows the variation of E0 and E0,h9 values in each of the four sections, excluding the dynamics of each point separately. The tendency of E0 and E0,h9 were similar in all tested PS. The seasonal increase and decrease in the E0 values over time were observed. It can be assumed that this seasonal impact in time was influenced by change of the temperature and moisture content in the soil of the subgrade during the year. When comparing the values of E0 and E0,h9, it was observed that even after recalculating E0 to E0,h9 it was not always possible to avoid seasonal change in the data. The lowest average E0 value measured in July 2014 was 871 MPa (air temperature was 25 °C, Table 2) and the highest was 1112 MPa in November 2013 (air temperature was 5 °C, Table 2) and the lowest average E0,h9 value measured in April 2015 was 834 MPa (air temperature was 9 °C, Table 2) and the highest was 969 MPa in November 2013 (air temperature was 5 °C Table 2; Figure 14 and Figure 15). That also reflects in the road PS No 12, where the deviation from the mean of the E0,h9 values was almost 5% in July 2014 and 8% in April 2015. At the same time of measurement, the similar deviations from the mean were observed in the road PS No 4, the deviations from the mean were 3% and 9% (Figure 14 and Figure 15), respectively. The type of improvement of the subgrade could influence this difference because in different sections the soil of the subgrade was improved with appropriately used asphalt granules or mineral materials of different fractions.
By analyzing the data of Figure 12 and Figure 13, both the dynamics of individual values of E0 and E0,h9, as well as the mean values, it was found that by neutralizing the influence of asphalt temperature on the bearing capacity of the road PS, the effect of seasonal impact on the weakening of the subgrade remained. Due to the lower ambient temperature, the more strength asphalt in the autumn and spring seasons was not sufficient to compensate for the decrease in the bearing capacity of the subgrade due to the increased moisture content of the soil in the subgrade. Even a subgrade of no less than 100 MPa, when such a bearing capacity was only achieved by improving the soil in the subgrade, did not ensure the complete elimination of seasonal effects. Therefore, only reinforcement of the soil in the subgrade should be applied in the road PS of Class I (where ESAL100 = (10.0–32.0) million) and Class SV (where ESAL100 is more than 32.0 million), and the improvement or reinforcement of the soil in the subgrade should be applied in the road PS of Class III (where ESAL100 = (0.8–3.0) million) and Class II (where ESAL100 = (3.0–10.0) million). Analyzing data more in details was decided that this is a seasonal effect on the change in bearing capacity of the road PS.
The values of E0 and E0,h9 of all five road PS in sections measured at different measuring periods were combined to conclude the entire Test Road. Therefore, an analysis in a time of the variation of E0 and E0,h9 values together was carried out (Figure 16). According to Figure 16, it was seen that the same tendency as in the analysis of data for each section separately, E0 and E0,h9 mean values changed over time (increasing and decreasing) without seasonal adjustment.
By analyzing the dynamics of E0 and E0,h9 values in each individual point (left wheel path, right wheel path and track) on loaded and unloaded lanes of each section were found that in the loaded lane, the dispersion of the measuring values of individual points was higher than in the unloaded lane (Figure 17, Figure 18 and Figure 19). In unloaded lanes, the E0 values at each point had small dispersion, i.e., the curves representing the E0 dynamics at each point of the PS were somewhat flattened, while in loaded lanes at different points of the same section E0 was entirely different, i.e., the curves depicting the E0 dynamics at each point of the road PS seemed to be separated from each other.
Therefore, it was assumed that in the loaded lane, the values of E0 and E0,h9 differed at individual points due to the uneven loading from the wheels of the motor vehicle. The uneven load was due to the movement of the motor vehicle wheel along the driving trajectory. As well as during the measurements, the FWD might have a partial deviation in time from the point of measurement, although the latter assumption was denied by the glare of the measured values in the unloaded lane.
Observed differences in distribution patterns between measured data on loaded and unloaded lanes were further analyzed at a more aggregated level. The graphs of Figure 19 show the dynamics of the E0 and E0,h9 values at each measuring point by analyzing the data of all tested PS together but separating loaded and unloaded lanes. It can be seen that the values of the strength of the road PS in the unloaded lanes were higher and with less dispersion than in loaded lanes. From the given dynamics E0,h9 on different lanes, it could also be seen that the right wheel path of the unloaded lane had a more pronounced downward tendency of E0,h9 than the left wheel path or track. It could be suspected that in the right wheel path of the unloaded lane the strength of the road PS for some reason decreased faster. In the loaded lane, it was hard to see differences, as the point dynamics were quite chaotic.
Therefore, the effect of heavy vehicle loads had the most effect on the path trajectory of the right wheel, otherwise defined, on the right wheel path, which was closer to the shoulder. In this regard, it was recommended that during the construction or reconstruction of the road PS of Class III (where ESAL100 = (0.8–3.0) million) and a higher class, the width of the crushed dolomite course should be extended to the slope. Thus ensuring the stability and bearing capacity of the edge of the road and the layers of the road PS in the area of the shoulder.
The bearing capacity during measuring autumn and spring periods of the loaded lane (Figure 20) changed respectively up to 6.3% (in autumn) and up to 8.3% (in spring; in PS No 4), up to 4.9% and up to 8.7% (No 12), up to 0.8% and up to 2.1% (No 18) and up to 1.9% and up to 5.3% (No 24–25). The research results of the bearing capacity of the unloaded lane (Figure 19) changed respectively up to 0.5% (in autumn) and up to 12.5% (in spring; No 4), up to 5.5% and up to 9.3% (No 12), up to 2.2% and up to 1.5% (No 18) and up to 4.8% and up to 7.9% (No 24–25). The temperature of asphalt pavement surface was on average about 14–18 °C (in autumn) and 3–10 °C (in spring; Table 2). The bearing capacity results of summer measurement (when temperature of asphalt pavement surface was in average about 20–31 °C, Table 2) were on average up to 12.5% lower than in autumn and up to 6.1% lower than in spring (Figure 20 and Figure 21).

4.2. The Effect of the Thickness of the Asphalt Concrete Layers for the E0 and E0,h9

The higher total thickness of the AC layers in the loaded lane than in the unloaded lane, along with other factors, affected the bearing capacity of the tested PS, i.e., E0 numerical values. The numerical values of E0 varied from 757 to 1287 MPa (average 1012.00 ± 96.97 MPa) in the unloaded lane and from 870 to 1228 MPa (average 956.50 ± 97.10 MPa) in the loaded lane. The numerical values of E0,h9 varied from 703.60 to 1134.00 MPa (average 952.40 ± 69.51 MPa) in the unloaded lane and from 606 to 1135 MPa (average 901.20 ± 78.32 MPa) in the loaded lane.
By analyzing the recalculated values E0,h9, there was a significantly lower distribution of these values than the E0 (Figure 22 and Figure 23). The tendency that the mean value Ei of the PS in an unloaded lane in which the total average thickness of the AC layers was 165.80 ± 5.95 mm was more than the loaded lane in which the total average thickness of the AC layers was about 171.80 ± 7.83 mm that remained.
It has been observed in studies that, regardless of the overall thickness of the AC layers, the average value Ei of the tested PS in the loaded or unloaded lanes varies slightly, but when compared to the average values Ei of the loaded and unloaded PS, this tendency is not valid. In this way, it was concluded that the total heavy vehicle loads had a significant effect on the E of the tested PS.
Figure 24 shows that the temperatures and dispersion graphs of E0 and E0,h9 indicate that there was an E0 linear negative temperature dependence at a given depth (the higher the temperature, the lower the E0).
During the study, the average moisture content at a depth of 100 cm was 6.344% ± 3.68%, at a depth of 130 cm, 8.36% ± 4.21% and at a depth of 150 cm, 11.88% ± 3.67%. From moisture content at a given depth, E0 and at a given depth of moisture content and the E0,h9 dispersion graph (Figure 25), it can be seen that neither E0 nor E0,h9 were related to moisture content.

4.3. Fitting Linear Mixed Models (LMMs)

In this statistical analysis, applying LMMs to clustered data, in which the units of analysis (points in the road) were within clusters (sample of four different road sections), and repeated measures in these different points were collected over time with unequally spaced time intervals. The data were balanced—each cluster had the same number of units, which were measured at the same number of time points. Points in the section of the road PS were dependent, as this dependence was due to the different and specific characteristics of the particular road PS. Due to the data and different clusters in which data could be correlated, linear regression models (LRMs) were not appropriate here, and therefore LMMs were applied. LMMs, like LRMs, are linear in the parameters, but independent variables in LMMs involve not only fixed but also random effects [28]. Fixed effects are treated the same as regression coefficients in LMMs—they indicate the relationships of the covariates with the continuous outcome variable. Random effects are used in modeling the random variation in the dependent variable at different levels of the data. This random variation occurs because each point in the section of the road PS is very individual and examining even the same characteristics of these points; these characteristics depend not only on the features included in the research but also on the individual, unknown peculiarities of the researchers. Thus, the expression of the LMM for the ith subject (n is the number of all survey subjects), whose measurements of the test variable values are repeated k times, is [30] (Equations (2)–(3)):
Y i = X i β + Z i u i + ε i , i = 1 , , n , u i ~ N ( 0 , D ) , ε i ~ N ( 0 , R i ) ,
Y i = ( Y 1 i Y 2 i Y k i ) ,   X i = ( X 1 i ( 1 ) X 2 i ( 1 ) X k i ( 1 ) X 1 i ( p ) X 2 i ( p ) X k i ( p ) ) ,   β = ( β 1 β 2 β p ) ,   Z i = ( Z 1 i ( 1 ) Z 2 i ( 1 ) Z k i ( 1 ) Z 1 i ( q ) Z 2 i ( q ) Z k i ( q ) ) ,   u i = ( u 1 i u 2 i u q i ) ,   ε i = ( ε 1 i ε 2 i ε k i ) ,
where Yi—dependent variable. The dependent variable in this study is the surface modulus E0 or, according to the method proposed by previous researchers, recalculated E0,h9;
Xi—a matrix of explanatory variables in which the number of variables p corresponds to the number of fixed parameters and their interactions, included in the model. The explanatory variables of this study were of two types: (1) road section, lane, factor of the load on the road PS and thickness of the AC layers are indications whose values are unchanging over time; and (2) the temperature of different types and moisture content, as well as the values of the dependent variable, vary depending on weather conditions, during a particular measurement (Table 4);
Zi—matrix of random effects q. For example, if the model contains only two random parameters, one of which is for estimating the randomness of the intercept and the other one for evaluating the randomness of slope between the subjects, then the random effects matrix Zi consists of two columns, i.e., the unit column and column of the time moments, respectively: Z i = ( 1 1 1 1 2 k ) ,
with the vector of parameters u i = ( u 1 i u 2 i ) .
Assuming that the mean of the random parameters is zero, the random variables are estimated with a covariance matrix D = ( V a r ( u 1 i ) cov ( u 1 i , u 2 i ) cov ( u 2 i , u 1 i ) V a r ( u 2 i ) ) , whose diagonal elements correspond to the variance of the parameters u 1 i and u 2 i , and non-diagonal elements correspond to the covariance among these random parameters.
εi—model residuals, where the covariance matrix Ri defines residuals dependence due to repeated measurements. This matrix could be diagonal, R i = V a r ( ε i ) = ( σ 2 σ 2 ) , k × k , which shows the dispersion of residuals of the ith subject at each measurement moment, but it can also be of a different structure [30].
LMMs are implemented in many software packages. Appropriate procedures could be used in most popular statistical packages: SAS, SPSS and R [28]. This statistical analysis used the R program package and its library nlme [31].
This study used a top-down LMM model specification strategy [28]. First of all the model with all fixed parameters were estimated (time, lane, load and interactions of these variables), in the second step possible random effects were included in the model and finally appropriate covariance matrix structure for model residuals were estimated. To evaluate which of the two models was more appropriate the Likelihood ratio test was used [28].

4.3.1. M1 Model for Dependent Variable E0,h9

In the model M1, E0,h9 was the dependent variable, whose values were recalculated from E0 removing the effect of the recent temperature on the AC layer at a depth of 9 cm and taking into account the thickness of the AC layer. This correction of the bearing capacity should smooth the dynamics of the initial E0 values in time, eliminating the apparent seasonal impact. However, Table 3 shows that seasonal change had not been successfully removed in all cases. Therefore, it was assumed that E0 values depended not only on the temperature of the AC layer but also on other seasonal factors that were not included in the study. Another reason might be that at a specific temperature, the bearing capacity was not very typical for the formula used for recalculation. Therefore, a slight seasonal impact was seen in the recalculated data.
Since E0,h9 had a slightly different level among the different pavement sections (Figure 12), the random intercept parameter between the pavement sections in the M1 was included. Figure 12 shows that in all of the studied pavement sections there was a decrease in the E0,h9 values in time, but this did not happen equally quickly across all sections (the slope was different). Due to this characteristic, the random trend among the sections was included in model M1. By taking all the points in each pavement section, their E0,h9 values were also scattered widely, and this could lead to the assumption that there was a certain randomness among the points that were not described by the fixed parameters. To estimate this randomness random intercept for individual points in the sections was included in model M1. Since the dynamics of E0,h9 in single points on individual sections were quite similar in time, the random trend of single points in the section was not included in M1.
In the final M1 model fixed variables (lane, load and time) and the random effects associated to pavement sections (clusters) and pavement points in road sections were included. To estimate residuals dependence in the model a diagonal covariance matrix with different dispersions of residuals at each time point was chosen The final model for E0,h9 in the i-point (i = 1, …, 60) on the j-section (j = 1, …, 4) at a given time t (t = 1, …, 8) had this form [29] Equation (4):
E 0 , h 9 t i j = β 0 + β 1 L t + β 2 K i j + β 3 T i j + β 4 A i j + β 5 L t × A i j + β 6 L t × K i j + β 7 L t × T i j + + u 0 j + u 1 j × L + u 0 i j + ε t i j ,
where L—time (measurement at a given time), K—indicator variable indicating the left wheel path; T—indicator variable indicating the track, A—indicator variable indicating loaded lane, L × A—interaction of time and load variables, which estimates the differences between the dynamics of the bearing capacity of the road PS on loaded and unloaded lanes; L × K—interaction of time and left wheel path variables, which estimates the differences between the dynamics of the bearing capacity of the road PS in the left and right wheel paths; L × T—interaction of the time and track variables, which estimates the differences between the dynamics of the bearing capacity of the road PS in the track and right wheel paths; β0 – β7—parameters that correspond to fixed effects in the model, i.e., intercept, time, lane, load and interactions of these effects; u 0 i j is a random effect associated to the difference of the measured points in the separate section of the road PS and ε i j is the model residuals.
The random parameters u 0 j and u 1 j describe the remaining differences of E0,h9, which have not been estimated by the fixed parameters in the M1 model. These random parameters correspond to the differences between the intercept and the direction of the trend between the sections. The mean of these random parameters is zero, and their variance and covariance are given in the D matrix. In this matrix, the diagonal elements correspond to the dispersion of the random effects u 0 j and u 1 j . Non-diagonal elements are the covariance among these random parameters. This model assumes that the vector ε ij for a given subject is a random variable that follows a multivariate normal distribution with a zero mean and a diagonal dispersion matrix R i j (Equation (5)):
ε i j = ( ε 1 i j ε 8 i j ) ~ N ( 0 , R i j ) , R i j = ( r 11 0 0 0 0 r 22 0 0 0 0 0 0 0 0 r 88 ) .
All fixed parameters were statistically significant in the M1 model (Table 5). Negative parameters, which are related to time and load variables, show correspondingly the significant decrease of the bearing capacity of the road PS in time and the negative influence of load on the bearing capacity. At loaded points E0,h9 was an average of 88.5 MPa lower than unloaded points. The dynamics of the bearing capacity of loaded and unloaded points varied in time differently, i.e., at loaded points bearing capacity values decreased slower than the at unloaded points (parameter β5; Table 5). Eliminating the time effects, there were still significant differences of E0,h9 between different ruts: The bearing capacity of the road PS on the left wheel path and the track was less than on the right wheel path.

4.3.2. M2 Model for Dependent Variable E0

This model estimates the same linear mixed model as M1 (with the same fixed and random parameters), but the dependent variable in it is E0 values from which are not eliminated the influence of temperature and thickness of the AC layers. Therefore additional fixed parameters related to the temperature, thickness of the AC layers and interaction between temperature and moisture content were included in the M2 model respectively: β8, β9 and β10 parameters. The negative correlation between E0 and different type of temperatures was observed. Since the various types of temperature strongly correlated with each other (correlation >0.9), only the temperature in the AC layer at a depth of 9 cm was included in the M2 model. This temperature was used to calibrate the E0,h9 values that were used in the M1 model (Equation (2)). Thus, the expression for model M2 is (Equation (6)):
E 0 t i j = β 0 + β 1 L t + β 2 K i j + β 3 T i j + β 4 A i j + β 5 L t × A i j + β 6 L t × K i j + β 7 L t × T i j + β 8 T E M P + β 9 A S P H + β 10 T E M P × M O I S T U R E + u 0 j + u 1 j × L + u 0 i j + ε t i j .
Almost all fixed parameters in the M2 model were statistically significant with similar parameters estimates as in the M1 model (Table 5). The additional variables in the M2 model, temperature and the thickness of the AC layer had accordingly the negative, and the positive effect on the bearing capacity of the road pavement construction. When the values of all other parameters in the model were fixed, then increasing the temperature by 1 °C at a depth of 9 cm in the AC layer influenced E0 decreases by an average of 8.41 MPa, and increasing thickness of the AC layer by 10 mm, E0 increases by an average of almost 44 MPa. The interaction of the temperature and the moisture content was statistically non-significant (p = 0.391) in the M2 model.
The estimation of random variables parameters (dispersion and covariance) of the M1 and M2 models are also presented in Table 5. It can be seen that the random parameters of a section were interdependent; their correlation was −0.594 in the M1 model and −0.34 in the M2 model. Therefore, if a higher intercept is estimated in a particular section, then there is a lower slope in this section estimated. In other words, if the bearing capacity of a road PS is higher in a particular section, then values in time decrease slower (or vice versa. In this case, the correlation does not indicate causality). Statistically significant random parameters between measured points in each section could be related to different loading rate of each measured point in the section, which is not included in the survey data. It is known only that the point on the road pavement is loaded or not, but how much it is loaded is not known.
The residual values of the models presented in Figure 26 could be used to determine how adequate the fitted models are. Figure 26 shows that the residual values of the M1 model were more concentrated than the M2 model. Such dispersion might be affected by the fact that for the M1 model recalculated E0,h9 values, with a smaller dispersion than initial E0 values were used.
On the other hand, at this stage, the purpose of the LMM models was not the prediction of the strength of the road PS, but the assessment of the strength differences between the different places in the road PS and properties of the road PS. Further research needs to include additional variables that describe the road PS to predict the strength of the road PS.

5. Conclusions

  • The possibilities of the LMM model allowed evaluating the dependence of the bearing capacity of the road PS on various fixed characteristics: the lane (unloaded or loaded) and their load mode, and these interactions, by finding specific random differences between points. These differences were not explained by the fixed characteristics included in the model. One of the reasons was that, even with the evaluation of fixed parameters, the bearing capacity of the road PS was different at points, there might be certain factors not included in the survey, i.e., structural characteristics of the section (cluster) or load size for each point at section (cluster). In this analysis, the load characteristic had only two values, i.e., loaded and unloaded. It was not known how much each point was loaded.
  • Heavy vehicle loads affected the performance of road pavement and this effect was more significant (aggressive) in the trajectory of the path of the right wheels. Due to the lower bearing capacity of roadside up to 10% (the average E0 was 911 MPa at the right wheel path, and the average E0 was 984 MPa at the left wheel path) and ending of unbound base layer at the edge of asphalt and higher wheel stress asphalt wearing and binder layers rutted much more in right wheel path. Considering this, it was recommended for roads with ESAL’s of 1.0 mln and higher to prolong unbound base layer until slope of the roadside; this would increase structural strength of PS and decrease the depth of wheel stress.
  • The usage of fine-grained sand 0/4 (the average E0 was 984 MPa) for the frost-blanket layer did not influence the PS bearing capacity and hydrothermal conditions, comparing with PS with the frost-blanket layer of sand 0/11 (the average E0 was 984 MPa). A difference of the bearing capacity was about 6%.
  • Independently from asphalt layers overall thickness, E0 modulus average value in loaded or unloaded traffic lane changed slightly on average 5% (on average 932–976 MPa of the loaded lane and 990–1034 MPa of the unloaded lane). However, comparing the E0 modulus average values of loaded and unloaded traffic lanes, this tendency was not valid (the difference was up to 15%, 967 MPa of loaded lane and 1108 MPa of the unloaded lane). This way it was concluded that total traffic loads had a significant impact to the E modulus of PS.
  • The most significant influence for pavement structural strength had the bearing capacity of the subgrade, susceptibility to hydrothermal impact, the thickness of asphalt layers and thickness of the whole PS. The unbound and bound layers only slightly influenced pavement structural strength.
  • Analysis of overall E0 and E0,h9 data declared that seasonal impact on pavement structural strength due to a change of subgrade bearing capacity remained after correction of asphalt stiffness dependent on layer temperature. This assigned that stiffer asphalt layers in autumn and springtime were insufficient to compensate for the decrease of pavement structural strength due to decreased bearing capacity of subgrade. The bearing capacity during measuring autumn and spring periods of loaded lane changed respectively up to 6.3% (in autumn) and up to 8.3% (in spring; in PS No 4), up to 4.9% and up to 8.7% (No 12), up to 0.8% and up to 2.1% (No 18) and up to 1.9% and up to 5.3% (No 24–25). The temperature of asphalt pavement surface was on average about 14–18 °C (in autumn) and 3–10 °C (in spring). The bearing capacity results of the summer measurement (when temperature of asphalt pavement surface was on average about 20–31 °C) were one average up to 12.5% lower than in autumn and up to 6.1% lower than in spring.
  • A subgrade improvement with the aggregates mix reaching deformation modulus EV2 higher than 100 MPa did not ensure the complete elimination of seasonal impact. Due to this, it was recommended for roads with ESAL’s of 3.0 mln and higher to use subgrade stabilization and for roads with ESAL’s of lower than 3.0 mln to use subgrade improvement or stabilization.
  • It was detected that neither E0 nor E0,h9 were related to moisture content at a depths of 100 cm, 130 cm and 150 cm.
  • This analysis was an initial step towards more detailed testing of road PS in the Test Road, it found some of the characteristics of the road PS and allowed for the adoption of appropriate engineering solutions.

Author Contributions

Conceptualization, A.V.; methodology, A.V., L.Ž., and J.Ž.; software, J.Ž. and L.Ž.; validation, J.Ž. and A.V.; formal analysis, J.Ž. and L.Ž.; investigation, L.Ž. and A.V.; resources, A.V.; data curation, A.V.; writing—original draft preparation, L.Ž. and J.Ž.; writing—review and editing, A.V. and D.Ž.; visualization, A.V.; supervision, A.V.; project administration, A.V.; funding acquisition, A.V.

Funding

This research were done at the Test Road that was funded by Lithuanian Road Administration under the Ministry of Transport and Communications fund.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fonod, A. Degradation of the bearing capacity of asphalt pavements. Slovak J. Civ. Eng. 2005, 2, 35–39. [Google Scholar]
  2. Amorim, S.I.R.; Pais, J.C.; Vale, A.C.; Minhoto, M.J.C. A model for equivalent axle load factors. Int. J. Pavement Eng. 2015, 16, 1, 881–893. [Google Scholar] [CrossRef]
  3. Straube, E.; Jansen, D. Temperature correction of falling-weight-deflectometer measurements. In Proceedings of the 8th International Conference Bearing Capacity of Roads, Railways and Airfields, Two Volume Set (BCR2A’09), Champaign, IL, USA, 29 June–2 July 2009; Tutumluer, E., Al-Qadi, I.L., Eds.; CRC: London, UK, 2009; pp. 789–798. [Google Scholar]
  4. Salour, F. Moisture Influence on Structural Behaviour of Pavements: Field and Laboratory Investigations. Ph.D Thesis, KTH Royal Institute of Technology, Stockholm, Sweden, 2015. [Google Scholar]
  5. Hermansson, Å. Simulation model for calculating pavement temperatures including maximum temperature. Transp. Res. Rec. J. Transp. Res. Board 2000, 1699, 134–141. [Google Scholar] [CrossRef]
  6. Alawi, M.H.; Helal, M.M. A mathematical model for the distribution of heat through pavement layers in Makkah roads. J. King Saud Univ. Eng. Sci. 2014, 26, 41–48. [Google Scholar] [CrossRef]
  7. Motiejūnas, A.; Paliukaitė, M.; Vaitkus, A.; Čygas, D.; Laurinavičius, A. Research on the Dependence of Asphalt Pavement Stiffness upon the Temperature of Pavement Layers. Baltic J. Road Bridge Eng. 2010, 5, 50–54. [Google Scholar] [CrossRef]
  8. García, J.A.R.; Castro, M. Analysis of the temperature influence on flexible pavement deflection. Constr. Build. Mater. 2011, 25, 3530–3539. [Google Scholar] [CrossRef]
  9. Zheng, Y.; Zhang, P.; Liu, H. Correlation between pavement temperature and deflection basin form factors of asphalt pavement. Int. J. Pavement Eng. 2019, 20, 874–883. [Google Scholar] [CrossRef]
  10. Hermansson, Å.; Charlier, R.; Collin, F.; Erlingsson, S.; Laloui, L.; Sršen, M. Heat transfer in soils. In Water in Road Structures; Springer: Dordrecht, The Netherlands, 2009; pp. 69–79. [Google Scholar] [CrossRef]
  11. Bobes-Jesus, V.; Pascual-Muñoz, P.; Castro-Fresno, D.; Rodriguez-Hernandez, J. Asphalt solar collectors: A literature review. Appl. Energy 2013, 102, 962–970. [Google Scholar] [CrossRef]
  12. Javilla, B.; Mo, L.; Hao, F.; Shu, B.; Wu, S. Significance of initial rutting in prediction of rutting development and characterization of asphalt mixtures. Constr. Build. Mater. 2017, 153, 157–164. [Google Scholar] [CrossRef]
  13. Korkiala-Tanttu, L.; Dawson, A. Relating full-scale pavement rutting to laboratory permanent deformation testing. Int. J. Pavement Eng. 2007, 8, 19–28. [Google Scholar] [CrossRef]
  14. Wiman, L.G. Accelerated Load Testing of Pavements: HVS-NORDIC Tests at VTI Sweden 2003–2004. VTI Rapport 544A. 2006. Available online: http://www.diva-portal.org/smash/get/diva2:670625/FULLTEXT01.pdf (accessed on 8 April 2019).
  15. Qiao, Y.; Dawson, A.; Huvstig, A.; Korkiala-Tanttu, L. Calculating rutting of some thin flexible pavements from repeated load triaxial test data. Int. J. Pavement Eng. 2015, 16, 467–476. [Google Scholar] [CrossRef]
  16. Birgisson, B.; Ruth, B. Improving performance through consideration of terrain conditions: Soils, drainage, and climate. Transp. Res. Rec. J. Transp. Res. Board 2003, 1819, 369–377. [Google Scholar] [CrossRef]
  17. Charlier, R.; Hornych, P.; Sršen, M.; Hermansson, Å.; Bjarnason, G.; Erlingsson, S.; Pavšič, P. Water influence on bearing capacity and pavement performance: Field observations. In Water in Road Structures; Springer: Dordrecht, The Netherlands, 2009; pp. 175–192. [Google Scholar] [CrossRef]
  18. Dawson, A.; Kringos, N.; Scarpas, T.; Pavšič, P. Water in the pavement surfacing. In Water in Road Structures; Springer: Dordrecht, The Netherlands, 2009; pp. 81–105. [Google Scholar] [CrossRef]
  19. Hall, K.T.; Crovetti, J.A. Effects of Subsurface Drainage on Pavement Performance: Analysis of the SPS-1 and SPS-2 Field Sections. Trans. Res. Board 2007, 583, 190. [Google Scholar] [CrossRef]
  20. Charlier, R.; Laloui, L.; Brenčič, M.; Erlingsson, S.; Hansson, K.; Hornych, P. Modelling coupled mechanics, moisture and heat in pavement structures. In Water in Road Structures; Springer: Dordrecht, The Netherlands, 2009; pp. 243–282. [Google Scholar] [CrossRef]
  21. Chen, D.H.; Sun, R.; Yao, Z. Impacts of aggregate base on roadway pavement performances. Constr. Build. Mater. 2013, 48, 1017–1026. [Google Scholar] [CrossRef]
  22. Salour, F.; Erlingson, S. Investigation of a pavement structural behaviour during spring thaw using falling weight deflectometer. Road Mater. Pavement Des. 2013, 14, 141–158. [Google Scholar] [CrossRef]
  23. Čygas, D.; Laurinavičius, A.; Vaitkus, A.; Paliukaitė, M.; Motiejūnas, A. A test road section of experimental pavement structures in Lithuania (II). In Proceedings of the 8th International Conference “Environmental Engineering”, Vilnius, Lithuania, 19–20 May 2011; Technika Press: Vilnius, Lithuania; pp. 1205–1209, ISBN 9789955288299. [Google Scholar]
  24. Vaitkus, A.; Laurinavičius, A.; Oginskas, R.; Motiejūnas, A.; Paliukaitė, M.; Barvidienė, O. The road of experimental pavement structures: Experience of five years operation. Baltic J. Road Bridge Eng. 2012, 7, 220–227. [Google Scholar] [CrossRef]
  25. Čygas, D.; Laurinavičius, A.; Paliukaitė, M.; Motiejūnas, A.; Žiliūtė, L.; Vaitkus, A. Monitoring the mechanical and structural behavior of the pavement structure using electronic sensors. Comput.-Aided Civ. Infrastr. Eng. 2015, 30, 317–328. [Google Scholar] [CrossRef]
  26. Vaitkus, A.; Laurinavičius, A.; Čygas, D.; Žiliūtė, L.; Šernas, O.; Židanavičiūtė, J. Asphalt Mix Composition Influence on Layers Performance—Test Road Experience. In Proceedings of the XXVth World Road Congress Roads and Mobility Creating New Value from Transport, Coex, Seoul, Korea, 2−6 November 2015. [Google Scholar]
  27. Vaitkus, A.; Paliukaitė, M. Evaluation of time loading influence on asphalt pavement rutting. Procedia Eng. 2013, 57, 1205–1212. [Google Scholar] [CrossRef]
  28. West, B.T.; Welch, K.B.; Galecki, A.T. Linear Mixed Models: A Practical Guide Using Statistical Software, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2015; p. 42. [Google Scholar]
  29. Židanavičiūtė, J.; Vaitkus, A. Tiesinių mišriųjų modelių taikymas kelio dangos savybių pakartotiniuose matavimuose. Lith. J. Stat. 2015, 54, 101–109. [Google Scholar]
  30. Littell, R.C.; Pendergast, J.; Natarajan, R. Tutorial in biostatistics: Modelling covariance structure in the analysis of repeated measures data. Stat. Med. 2000, 19, 1793–1819. [Google Scholar] [CrossRef]
  31. Pinheiro, J.; Bates, D.; DebRoy, S.; Sarkar, D.; R Core Team. Linear and Nonlinear Mixed Effects Models. R Package Version 3.1-137. 2018. Available online: https://CRAN.R-project.org/package=nlme (accessed on 9 July 2019).
Figure 1. PSs of the Test Road: No 4, No 12, No 18 and No 24–25.
Figure 1. PSs of the Test Road: No 4, No 12, No 18 and No 24–25.
Materials 12 03826 g001
Figure 2. The layout of points of measured deflection in the tested PS No 4 (where: R/L—right wheel path of the loaded lane, T/L—track of the loaded lane, L/L—left wheel path of the loaded lane, R/U—right wheel path of the unloaded lane, T/U—track of the unloaded lane and L/U—left wheel path of the unloaded lane. Distances between measuring points are given in meters (m)
Figure 2. The layout of points of measured deflection in the tested PS No 4 (where: R/L—right wheel path of the loaded lane, T/L—track of the loaded lane, L/L—left wheel path of the loaded lane, R/U—right wheel path of the unloaded lane, T/U—track of the unloaded lane and L/U—left wheel path of the unloaded lane. Distances between measuring points are given in meters (m)
Materials 12 03826 g002
Figure 3. Statistical analysis of the experiment.
Figure 3. Statistical analysis of the experiment.
Materials 12 03826 g003
Figure 4. The box-plot of the total thickness of asphalt concrete layers in loaded and unloaded lanes.
Figure 4. The box-plot of the total thickness of asphalt concrete layers in loaded and unloaded lanes.
Materials 12 03826 g004
Figure 5. The box-plot of the total thickness of the asphalt concrete (AC) layers in different PSs on different lanes.
Figure 5. The box-plot of the total thickness of the asphalt concrete (AC) layers in different PSs on different lanes.
Materials 12 03826 g005
Figure 6. The box-plot of the measured temperature.
Figure 6. The box-plot of the measured temperature.
Materials 12 03826 g006
Figure 7. The box-plot of moisture content of the soils in the subgrade.
Figure 7. The box-plot of moisture content of the soils in the subgrade.
Materials 12 03826 g007
Figure 8. The box-plot of the surface and separate layer modulus of the road PS.
Figure 8. The box-plot of the surface and separate layer modulus of the road PS.
Materials 12 03826 g008
Figure 9. The box-plots of E0 and E0,h9 values in each tested road section.
Figure 9. The box-plots of E0 and E0,h9 values in each tested road section.
Materials 12 03826 g009
Figure 10. The box-plots of E0 and E0,h9 values in loaded and unloaded lanes.
Figure 10. The box-plots of E0 and E0,h9 values in loaded and unloaded lanes.
Materials 12 03826 g010
Figure 11. The box-plots of E0 and E0,h9 values in each loaded and unloaded lane by distinguishing the left wheel path, track and right wheel path.
Figure 11. The box-plots of E0 and E0,h9 values in each loaded and unloaded lane by distinguishing the left wheel path, track and right wheel path.
Materials 12 03826 g011
Figure 12. E0 and E0,h9 at individual measuring points [29].
Figure 12. E0 and E0,h9 at individual measuring points [29].
Materials 12 03826 g012
Figure 13. The distribution of individual E0 and E0,h9 values for each of the road PS.
Figure 13. The distribution of individual E0 and E0,h9 values for each of the road PS.
Materials 12 03826 g013
Figure 14. E0 and E0,h9 mean values and their fluctuations regarding the total average over time.
Figure 14. E0 and E0,h9 mean values and their fluctuations regarding the total average over time.
Materials 12 03826 g014
Figure 15. E0 and E0,h9 mean values and their fluctuations regarding the total average in each section.
Figure 15. E0 and E0,h9 mean values and their fluctuations regarding the total average in each section.
Materials 12 03826 g015
Figure 16. The box-plots of the distribution of individual E0 and E0,h9 values of the tested PS.
Figure 16. The box-plots of the distribution of individual E0 and E0,h9 values of the tested PS.
Materials 12 03826 g016
Figure 17. The dynamic of the individual E0 values of each section for each point of measurement on the right and left wheel paths and in the loaded and unloaded lanes.
Figure 17. The dynamic of the individual E0 values of each section for each point of measurement on the right and left wheel paths and in the loaded and unloaded lanes.
Materials 12 03826 g017
Figure 18. The dynamic of the individual E0,h9 values of each section for each point of measurement on the right and left wheel paths and in the loaded and unloaded lanes.
Figure 18. The dynamic of the individual E0,h9 values of each section for each point of measurement on the right and left wheel paths and in the loaded and unloaded lanes.
Materials 12 03826 g018
Figure 19. The dynamics of E0 and E0,h9 values of each point of measurement in unloaded and loaded lanes by analyzing the data of all road PS together.
Figure 19. The dynamics of E0 and E0,h9 values of each point of measurement in unloaded and loaded lanes by analyzing the data of all road PS together.
Materials 12 03826 g019
Figure 20. The average values of bearing capacity in different periods of the loaded lane of PS No 4, No 12, No 18 and No 24–25.
Figure 20. The average values of bearing capacity in different periods of the loaded lane of PS No 4, No 12, No 18 and No 24–25.
Materials 12 03826 g020
Figure 21. The average values of bearing capacity in different periods of the unloaded lane of PS No 4, No 12, No 18 and No 24–25.
Figure 21. The average values of bearing capacity in different periods of the unloaded lane of PS No 4, No 12, No 18 and No 24–25.
Materials 12 03826 g021
Figure 22. Influence of the total thickness of asphalt concrete layers on E0 by analyzing each section individually in loaded and unloaded lanes.
Figure 22. Influence of the total thickness of asphalt concrete layers on E0 by analyzing each section individually in loaded and unloaded lanes.
Materials 12 03826 g022
Figure 23. Influence of total thickness of asphalt concrete layers on E0 and E0,h9 values in loaded and unloaded lanes.
Figure 23. Influence of total thickness of asphalt concrete layers on E0 and E0,h9 values in loaded and unloaded lanes.
Materials 12 03826 g023
Figure 24. Scatter plots for E0 (a) and E0,h9 (b) values and moisture content at a given depth.
Figure 24. Scatter plots for E0 (a) and E0,h9 (b) values and moisture content at a given depth.
Materials 12 03826 g024aMaterials 12 03826 g024b
Figure 25. Scatter plots for E0 (a) and E0,h9 (b) values and temperature content at a given depth.
Figure 25. Scatter plots for E0 (a) and E0,h9 (b) values and temperature content at a given depth.
Materials 12 03826 g025
Figure 26. Residual values of models M1 and M2: (a) diagnostic of estimated M1 model: predicted values of E0,h9 versus standardized residuals; (b) diagnostic of estimated M2 model: predicted values of E0 versus standardized residuals; (c) diagnostic of estimated M1 model: residuals versus quantiles of standard normal distribution; (d) diagnostic of estimated M2 model: residuals versus quantiles of standard normal distribution.
Figure 26. Residual values of models M1 and M2: (a) diagnostic of estimated M1 model: predicted values of E0,h9 versus standardized residuals; (b) diagnostic of estimated M2 model: predicted values of E0 versus standardized residuals; (c) diagnostic of estimated M1 model: residuals versus quantiles of standard normal distribution; (d) diagnostic of estimated M2 model: residuals versus quantiles of standard normal distribution.
Materials 12 03826 g026
Table 1. Positions of temperature and moisture content sensors in road pavement structures (PSs).
Table 1. Positions of temperature and moisture content sensors in road pavement structures (PSs).
No of the Road PS Road PSSoil of Subgrade
(Continuous Measurements)
TemperatureTemperatureMoisture Content
No 4on the top;
at a depth of 9 cm
at a depth of 100 cm; 130 cm; 150 cm
No 12on the top
at a depth of
8 cm; 9 cm; 10 cm
at a depth of 100 cm; 130 cm; 150 cm
No 18on the top;
at a depth of
8 cm; 9 cm; 10 cm
at a depth of
100 cm; 130 cm; 150 cm
no sensor
No 24–25on the top;
at a depth of 9 cm
at a depth of 100 cm; 130 cm; 150 cm
Table 2. Date of deflection measurements in sections.
Table 2. Date of deflection measurements in sections.
Measurement NoDateESAL100Average Air Temperature during Measurement Time, °CAverage Pavement Surface Temperature during Measurement Time, °CAverage Asphalt Concrete Layer Temperature at a Depth of 9 cm during Measurement Time, °C
122.08.2013484,000222021
210.10.2013508,000131012
314–15.11.2013537,000535
407.05.2014596,000121412
514.07.2014621,000253124
616.09.2014637,000192318
713.11.2014654,000776
822–23.04.2015678,00091810
Table 3. Correlation of temperatures and moisture content at different depths of the asphalt concrete layer.
Table 3. Correlation of temperatures and moisture content at different depths of the asphalt concrete layer.
E0E0,h9TemperatureMoisture ContentThe Thickness of the Asphalt Layer
At a Depth of 100 cmAt a Depth of 130 cmAt a Depth of 150 cmOn the Top of RPS (Measured Manually)At a Depth of 9 cm of RPS (Measured Manually)At a Depth of 100 cmAt a Depth of 130 cmAt a Depth of 150 cm
E01.000.75−0.38−0.39−0.37−0.53−0.54−0.09−0.120.080.07
E0,h90.751.000.250.240.240.030.140.07−0.100.010.12
TemperatureAt a depth of 100 cm−0.380.251.000.990.990.590.890.09−0.17−0.03−0.03
At a depth of 130 cm−0.390.240.991.001.000.580.890.04−0.21−0.100
At a depth of 150 cm−0.370.240.991.001.000.550.870.03−0.23−0.110
On the top of RPS (measured manually)−0.530.030.590.580.551.000.840.380.3800.03
At a depth of 9 cm of RPS (measured manually)−0.540.140.890.890.870.841.000.230.07−0.120.03
Moisture contentAt a depth of 100 cm−0.090.070.090.040.030.380.231.000.760.26−0.02
At a depth of 130 cm −0.12−0.10−0.17−0.21−0.230.380.070.761.00−0.020.03
At a depth of 150 cm0.080.01−0.03−0.10−0.110−0.120.26−0.021.00−0.16
The thickness of the asphalt layer0.070.12−0.03000.030.03−0.020.03−0.161.00
Table 4. Types of variables in the study.
Table 4. Types of variables in the study.
Dependent VariableExplanatory Variables
Tested PSVariables Related to the Features of Measured Points in the PSTime-Varying Properties
E0 or E0,h9Identification number of the tested PS (section)Identification number of the point in the tested PS.
Lane (left wheel path, track and right wheel path).
Loading feature (Yes/No).
The thickness of the AC layer.
Time (of repeated measurements).
Temperature (°C) and moisture content (%) in the soil of the subgrade at a depth of 100 cm, 130 cm and 150 cm.
Temperature (°C) on the top of the road PS.
Temperature (°C) of the AC layer at a depth of 9 cm.
Table 5. Estimated parameters of models M1 and M2.
Table 5. Estimated parameters of models M1 and M2.
Model ParametersEstimated Parameters
M1 ModelM2 Model
EstimateStandard Errorp-ValueValueStandard Errorp-Value
β0 (intercept)1002.6511.19<0.001454.9679.89<0.001
β1 (time)−17.271.76<0.001−18.721.99<0.001
β2 (left wheel path)−21.9710.040.029−21.3110.080.0356
β3 (track)−21.6410.060.0324−32.1310.140.0017
β4 (loaded lane)−88.458.21<0.001−115.118.71<0.001
β5 (time × loaded lane)9.331.01<0.00110.031.08<0.001
β6 (time × left wheel path)14.161.24<0.00114.991.32<0.001
β7 (time × track)15.581.24<0.00117.091.32<0.001
β8 (temperature)---−8.410.17<0.001
β9 (thickness of AC layer)---4.370.48<0.001
β10 (temperature of the AC layer at a depth of 9 cm × moisture content at a depth of 100 cm)---0.0270.030.391
Standard deviation of random parameters σ 2 and correlation (elements of matrix D)σintercept among sections = 14.83
σtrend among sections = 2.57
corintercept and trend among sections = −0.516
σintercept among points in the section = 48.55
σintercept among sections = 18.89
σtrend among sections = 3.30
corintercept and trend among sections = −0.428
σintercept among points in the section = 52.19
Note. β0β10—the value of the parameters estimated without the insignificant interaction of the temperature and the moisture content.

Share and Cite

MDPI and ACS Style

Vaitkus, A.; Žalimienė, L.; Židanavičiūtė, J.; Žilionienė, D. Influence of Temperature and Moisture Content on Pavement Bearing Capacity with Improved Subgrade. Materials 2019, 12, 3826. https://doi.org/10.3390/ma12233826

AMA Style

Vaitkus A, Žalimienė L, Židanavičiūtė J, Žilionienė D. Influence of Temperature and Moisture Content on Pavement Bearing Capacity with Improved Subgrade. Materials. 2019; 12(23):3826. https://doi.org/10.3390/ma12233826

Chicago/Turabian Style

Vaitkus, Audrius, Laura Žalimienė, Jurgita Židanavičiūtė, and Daiva Žilionienė. 2019. "Influence of Temperature and Moisture Content on Pavement Bearing Capacity with Improved Subgrade" Materials 12, no. 23: 3826. https://doi.org/10.3390/ma12233826

APA Style

Vaitkus, A., Žalimienė, L., Židanavičiūtė, J., & Žilionienė, D. (2019). Influence of Temperature and Moisture Content on Pavement Bearing Capacity with Improved Subgrade. Materials, 12(23), 3826. https://doi.org/10.3390/ma12233826

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop