A Study on Radiation Cooling Effect on Asphalt Concrete Pavement Using Basic Oxygen Furnace Slag to Replace Partial Aggregates

: In this study, aggregates in asphalt concrete were partially replaced by basic oxygen furnace slag (BOFS) in proportions of 45 wt.%, 55 wt.%, and 75 wt.%. The thermal performances of the specimens are discussed based on the thermal conductivity, emissivity, and the indoor and outdoor temperature measurements. Consequently, 75 wt.% of the specimen’s aggregates were replaced by BOFS, which had a high emissivity of 0.86 across the sky window. In the indoor and outdoor tests, the temperature change was recorded to estimate the thermal performance of specimens. According to the quantitative calculation, when the substitution of BOFS was higher than 55 wt.%, the specimens had a better radiation cooling ability. Among these specimens, the specimen with the BOFS substitution of 75 wt.% absorbed the most heat inside the body, contributing to less heat remaining in the environment. Furthermore, because Newton’s cooling energy accounted for about 90% of the stored energy within 7 h, the heat dissipation after the seventh hour was primarily radiation cooling, corresponding to the emission across the urban boundary layer. As for the mechanical properties, the stability value, indirect tensile strength, and British pendulum number (BPN) were in line with the speciﬁcations under the proper BOFS substitution. In conclusion, BOFS has great applicability in pavements due to its thermal performance and mechanical properties. It not only achieves the goal of urban heat island mitigation by radiation cooling, but also reﬂects the concept of resource sustainability.


Introduction
In this century with booming industries, the enormous exploitation and consumption of energy results in the ever-increasing concentration of carbon dioxide worldwide, causing the greenhouse effect and extreme climate changes. According to the statistics from the National Oceanic and Atmospheric Administration (NOAA), the average temperature has increased by up to 0.18 • C per decade since 1981. This has caused the snow in the polar regions to melt and subsequently raised the sea level [1]. Furthermore, the concentration of carbon dioxide started at about 265 ppm in 1850, reached 385 ppm in 2009, and even ascended to 407.4 ppm in 2018. The vicious circle of greenhouse gas emission has a tremendous impact not only on metropolises but also environments. Under the circumstances, building envelopes and pavements with high heat capacities bring on and even worsen the urban heat island effect.
Nowadays, to mitigate the urban heat island effect, many studies have proposed to utilize reflective coating of asphalt concrete and have recorded its ambient temperatures and internal temperatures at different depths, exploring the principle of cooling with its spectral characteristics [2,3]. One study changed the color of asphalt concrete to discuss its reflection characteristics [4]; some even improved the reflectivity across specific wavelengths, mainly aimed at visible and near infrared light [5][6][7]; others made multi-layer structure coatings and compared their cooling benefits [8]. Furthermore, adding titanium dioxide with different weight ratios and particle sizes could improve the reflectivity of coating [9,10]. Hollow glass microspheres were substituted for fillers inside the asphalt concrete, optimizing the proportion for the best cooling effect with a simulated solar heating test [11]. This replacement could improve concrete's emissivity and also reduce thermal conductivity [12]. The lower thermal conductivity materials painted on the asphalt concrete were capable of preventing heat from passing through the coating [13]. The conductive multi-layer structures contributed to diverse applications, such as a radiation-cooling pavement in summer and heat preservation on permafrost [14,15].
In terms of building envelopes, inorganic coating made from mineral powders has been painted on steel plates and concrete specimens. The temperatures and heat fluxes of the coated and uncoated surfaces and interiors were measured in the laboratory and in field tests and compared. The results showed that this coating had an obvious radiation cooling effect and great insulation ability owing to the emissivity and reflectivity. Furthermore, the adhesion and weather resistance were maintained at a certain level [16].
Apart from thermal properties, the mechanical properties of pavements need to be examined as well for practical use. Basic oxygen furnace slag (BOFS) is a byproduct of the steelmaking process. The converter is a by-product of the steelmaking process. A proper ratio for partially replacing aggregate with BOFS should be favorable for the development of bearing capacity and tensile strength when preparing asphalt or recycled concrete [17][18][19][20][21]. Furthermore, based on a two-year on-site test, it was found that the substitution of asphalt concrete entailed a longer life span with less damage in comparison with the original one [22]. Due to the roughness, multi-angularity, and rigidity of BOFS, the substitution concrete had better skid resistance, binding with bitumen, and also abrasion resistance [23]. Even when adopting the warm-mix or hot-mix asphalt approaches, this concrete still showed a great capability for resistance to deflection [24].
In this study, thermal and mechanical properties tests were done on the asphalt concrete specimens, in which aggregates were partially replaced by BOFS. Our previous research found that BOFS has a high far-infrared emissivity. After absorbing solar radiation, these specimens were capable of emitting the stored energy across the sky window (8-13 µm) to the upper sky, reducing energy accumulation and retention in the atmosphere, thereby achieving the goal of urban heat island (UHI) mitigation. The surface skid resistance is needed for the development of a radiation-cooling pavement. Ceramic particles or fine sand particles were added to increase the pavement roughness and the British pendulum number (BPN). BOFS is a byproduct of steelmaking that can not only achieve the goal of radiation cooling but also meet the needs of resource sustainability.

Materials and Their Properties
BOFS, a byproduct of steel industries, has lower thermal conductivity compared to natural aggregates. Under the same heat source and with an identical time duration, BOFS has a slower temperature-raising rate. Due to the greater hardness, BOFS applied to asphalt concrete has a bearing capacity comparable to natural aggregates.
In this study, five specimens were designed and prepared, including one general asphalt concrete specimen, and the stone aggregates of the asphalt concrete specimens were replaced by BOFS at proportions of 45 wt.%, 55 wt.%, 65 wt.%, and 75 wt.% respectively. The compositions of the aggregates used in the asphalt concrete specimens are shown in Table 1, and they are dense-graded and conform to ASTM D3515 [25]. The BOFS was also sieved with different sieve sizes. Then, we used BOFS aggregates to replace stone aggregates; the replacement percentage of every sieve size was the same. The content of asphalt to the aggregate mixture was 5.6 wt.%, meeting the specification of between 2 wt.% to 10 wt.%. The naming rules and descriptions of the specimens are shown in Table 2.

Emissivity
To identify the radiation cooling effect, the specimens (10 cm in diameter and less than 2 cm in thickness) were placed in the customized apparatus and then heated to 35 • C and 80 • C, respectively. To reflect the reality of usage, it was necessary to confirm that no aggregate was exposed on the measured surface. The environmental causes were reduced as much as possible in the measurements. When the specimens were heated, the infrared light emitted from the specimens was measured by Fourier transform infrared spectrometer (FT-IR) (Invenior, Bruker, Karlsruhe, Germany) in the range of 4000-600 cm −1 (3-16.67 µm). The raw data were the result of optical interference.
Through the Fourier transform, the emissivity of heated specimens was calculated compared to the standard blackbody. The emissivity of five specimens was classified into 3-16.7 µm and 8-13 µm, as shown in Figure 1a,b, respectively. The corresponding average emissivity for 3-16.67 µm and 8-13 µm is presented in Table 3. When the substitution of BOFS exceeded 65 wt.%, the emissivity was higher than the benchmark specimen. Among these, the BOF-75 specimen possessed the highest emissivity of 0.88 for 3-16.67 µm and 0.86 for 8-13 µm.

Thermal Conductivity
Thermal conductivity is one of the important parameters of thermal performance for materials. The conductivity was tested with a quick thermal conductivity meter (QTM-500, KEM, Tokyo, Japan) used on the specimens. After calibration, the slope equaling the temperature divided by the logarithmic time represents the thermal conductivity coefficient. The larger the slope, the faster the materials conduct heat, and vice versa. In this study, the thermal conductivity of each specimen was averaged.
As illustrated in Figure 2, the thermal conductivity showed a downward trend with the increasing replacement amount of BOFS. Among the specimens, the BOF-75 specimen had the lowest conductivity of 1.17 W/m-K, which means it had better heat-insulating properties than the others. As a consequence, BOFS was taken as a good heat-insulating material.

Thermal Conductivity
Thermal conductivity is one of the important parameters of thermal performance for materials. The conductivity was tested with a quick thermal conductivity meter (QTM-500, KEM, Tokyo, Japan) used on the specimens. After calibration, the slope equaling the temperature divided by the logarithmic time represents the thermal conductivity coefficient. The larger the slope, the faster the materials conduct heat, and vice versa. In this study, the thermal conductivity of each specimen was averaged.
As illustrated in Figure 2, the thermal conductivity showed a downward trend with the increasing replacement amount of BOFS. Among the specimens, the BOF-75 specimen had the lowest conductivity of 1.17 W/m-K, which means it had better heat-insulating properties than the others. As a consequence, BOFS was taken as a good heat-insulating material.

Thermal Conductivity
Thermal conductivity is one of the important parameters of thermal performance for materials. The conductivity was tested with a quick thermal conductivity meter (QTM-500, KEM, Tokyo, Japan) used on the specimens. After calibration, the slope equaling the temperature divided by the logarithmic time represents the thermal conductivity coefficient. The larger the slope, the faster the materials conduct heat, and vice versa. In this study, the thermal conductivity of each specimen was averaged.
As illustrated in Figure 2, the thermal conductivity showed a downward trend with the increasing replacement amount of BOFS. Among the specimens, the BOF-75 specimen had the lowest conductivity of 1.17 W/m-K, which means it had better heat-insulating properties than the others. As a consequence, BOFS was taken as a good heat-insulating material.

Experimental Measurement Methods
In this study, the experimental measurement was divided into two aspects: thermal performance and mechanical performance. The thermal performance was measured in the laboratory and also outdoors; the mechanical performance, including the stability value, indirect tensile strength, and British pendulum number, reflects the feasibility of practical application in engineering.

Indoor Temperature Measurement
To find the best radiation cooling among the benchmark and BOFS specimens, this experimental measurement was conducted in the laboratory and outdoors. The fabrication of the specimen started with a stir-fry pan and was evenly mixed, and the asphalt concrete was poured into the customized steel mold and then compacted to form a specimen. To simulate the actual pavement, the size of the specimen was molded as 50 × 50 cm and with a thickness of 5 cm. Meanwhile, the aggregate layer was 13 cm in thickness underneath the asphalt concrete.
To simulate a sunlight environment in the laboratory, halogen lamps and infrared lamps were used as the heat source in this apparatus to irradiate visible light and infrared light, which provide the majority of solar radiation. The ratios of visible light and infrared light remained at 44% and 53%. Two different radiation intensities, 623 W/m 2 and 436 W/m 2 , were tested. According to the previous results of outdoor measurements in summer at Taipei City, the radiation intensity of the heat source should be adjusted to 623 W/m 2 , composed of 279 W/m 2 from the halogen lamps erected at 30 cm above the specimen and 344 W/m 2 from the infrared lamp erected at 33 cm above the specimen. To simulate the surface temperature of pavement in winter at Taipei City, the radiation intensity was adjusted to 436 W/m 2 with the maximum surface temperature of 65 • C. Schematic diagrams of the apparatus used in the indoor temperature measurement under the radiation intensities of (a) 623 W/m 2 and (b) 436 W/m 2 are shown in Figure 3. application in engineering.

Indoor Temperature Measurement
To find the best radiation cooling among the benchmark and BOFS specimens, th experimental measurement was conducted in the laboratory and outdoors. The fabrica tion of the specimen started with a stir-fry pan and was evenly mixed, and the aspha concrete was poured into the customized steel mold and then compacted to form a spec men. To simulate the actual pavement, the size of the specimen was molded as 50 × 50 cm and with a thickness of 5 cm. Meanwhile, the aggregate layer was 13 cm in thickness un derneath the asphalt concrete.
To simulate a sunlight environment in the laboratory, halogen lamps and infrare lamps were used as the heat source in this apparatus to irradiate visible light and infrare light, which provide the majority of solar radiation. The ratios of visible light and infrare light remained at 44% and 53%. Two different radiation intensities, 623 W/m 2 and 43 W/m 2 , were tested. According to the previous results of outdoor measurements in sum mer at Taipei City, the radiation intensity of the heat source should be adjusted to 62 W/m 2 , composed of 279 W/m 2 from the halogen lamps erected at 30 cm above the spec men and 344 W/m 2 from the infrared lamp erected at 33 cm above the specimen. To sim ulate the surface temperature of pavement in winter at Taipei City, the radiation intensit was adjusted to 436 W/m 2 with the maximum surface temperature of 65 °C. Schemat diagrams of the apparatus used in the indoor temperature measurement under the radia tion intensities of (a) 623 W/m 2 and (b) 436 W/m 2 are shown in Figure 3.  The test was conducted in a shaded and airtight room; the specimen was surrounde by Styrofoam. To avoid errors from environmental factors, such as light and ambient tem perature, strict conditions were necessary. The thermal cables, which were attached o the surface or buried inside the asphalt concrete, recorded the temperature change at di ferent depths. The experimental apparatus and a schematic diagram of the sensor insta lation are shown in Figure 4a,b, respectively. The test was conducted in a shaded and airtight room; the specimen was surrounded by Styrofoam. To avoid errors from environmental factors, such as light and ambient temperature, strict conditions were necessary. The thermal cables, which were attached on the surface or buried inside the asphalt concrete, recorded the temperature change at different depths. The experimental apparatus and a schematic diagram of the sensor installation are shown in Figure 4a,b, respectively.

Outdoor Temperature Measurement
In order to measure their thermal performance, asphalt concrete specimens were exposed to solar radiation, and the outdoor test was conducted on the empty top floor, without shelter, in the summer in 2020 on the roof of the Department of Civil Engineering building, National Taipei University of Technology, Taipei City. The size of the specimens, 50 × 50 × 5 cm, was the same as for those used in the indoor test. Likewise, the aggregate layer was laid 13 cm in thickness beneath the specimens to simulate a real pavement. Also, Sustainability 2021, 13, 3708 6 of 25 thermal cables were installed to measure the temperature change, as shown in Figure 4b. The outdoor temperature measurement setup is shown in Figure 5.

Outdoor Temperature Measurement
In order to measure their thermal performance, asphalt concrete specimens were exposed to solar radiation, and the outdoor test was conducted on the empty top floor, without shelter, in the summer in 2020 on the roof of the Department of Civil Engineering building, National Taipei University of Technology, Taipei City. The size of the specimens, 50 × 50 × 5 cm, was the same as for those used in the indoor test. Likewise, the aggregate layer was laid 13 cm in thickness beneath the specimens to simulate a real pavement. Also, thermal cables were installed to measure the temperature change, as shown in Figure 4b. The outdoor temperature measurement setup is shown in Figure 5.

Mechanical Properties
The asphalt concrete specimen without BOFS, named the benchmark specimen, and the asphalt concrete with BOFS in the three ratios were prepared with a Marshall apparatus according to the ASTM D6926 standard [26,27]. First, following the mixing formula and substitution ratio, the asphalt and aggregates were mixed by stirring until uniform. The mixture was quickly poured into the steel mold. Then, through 75 times compaction on both sides and then cooling down for about one day, the asphalt concrete specimen was demolded. Its average diameter and height were measured with an electronic Vernier caliper, with an accuracy of 0.01 mm.
The stability value, indirect tensile strength, and British pendulum number are discussed below. Following AI SS-1 set by the Asphalt Institute, the stability values of specimens were measured and specimens were sunk into a constant 60 °C water tank for 30 min and tested with the computer-controlled automatic Marshall apparatus (Ye-Chance Enterprise Co., Taipei, Taiwan) at a uniform rate of 50.8 mm/min. As the pavement is

Outdoor Temperature Measurement
In order to measure their thermal performance, asphalt concrete speci posed to solar radiation, and the outdoor test was conducted on the empty t out shelter, in the summer in 2020 on the roof of the Department of Civ building, National Taipei University of Technology, Taipei City. The siz mens, 50 × 50 × 5 cm, was the same as for those used in the indoor test aggregate layer was laid 13 cm in thickness beneath the specimens to simul ment. Also, thermal cables were installed to measure the temperature cha in Figure 4b. The outdoor temperature measurement setup is shown in Fig   Figure 5. The outdoor temperature measurement equipment and setup.

Mechanical Properties
The asphalt concrete specimen without BOFS, named the benchmark the asphalt concrete with BOFS in the three ratios were prepared with a M ratus according to the ASTM D6926 standard [26,27]. First, following the m and substitution ratio, the asphalt and aggregates were mixed by stirring The mixture was quickly poured into the steel mold. Then, through 75 tim on both sides and then cooling down for about one day, the asphalt conc was demolded. Its average diameter and height were measured with an ele caliper, with an accuracy of 0.01 mm.
The stability value, indirect tensile strength, and British pendulum n cussed below. Following AI SS-1 set by the Asphalt Institute, the stability v mens were measured and specimens were sunk into a constant 60 °C wa min and tested with the computer-controlled automatic Marshall apparat Enterprise Co., Taipei, Taiwan) at a uniform rate of 50.8 mm/min. As th subjected to severe traffic loadings, the stability value must be higher tha comply with the standard.

Mechanical Properties
The asphalt concrete specimen without BOFS, named the benchmark specimen, and the asphalt concrete with BOFS in the three ratios were prepared with a Marshall apparatus according to the ASTM D6926 standard [26,27]. First, following the mixing formula and substitution ratio, the asphalt and aggregates were mixed by stirring until uniform. The mixture was quickly poured into the steel mold. Then, through 75 times compaction on both sides and then cooling down for about one day, the asphalt concrete specimen was demolded. Its average diameter and height were measured with an electronic Vernier caliper, with an accuracy of 0.01 mm.
The stability value, indirect tensile strength, and British pendulum number are discussed below. Following AI SS-1 set by the Asphalt Institute, the stability values of specimens were measured and specimens were sunk into a constant 60 • C water tank for 30 min and tested with the computer-controlled automatic Marshall apparatus (Ye-Chance Enterprise Co., Taipei, Taiwan) at a uniform rate of 50.8 mm/min. As the pavement is subjected to severe traffic loadings, the stability value must be higher than 8.006 kN to comply with the standard.
The indirect tensile strength was measured according to ASTM D6931 [28]. The specimen was attached between two load stripes and was loaded radially at a speed of 50 ± 5 mm/min. The width of stripes was 12.7 ± 0.3 mm, complying with the standard. The specimen was exposed to 25 • C water for 30-120 min, and the maximum load at fracture was measured. The indirect tensile strength can be calculated from Equation (1), as follows.
In Equation (1), S t is the indirect tensile strength (kPa), P is the maximum load (N), t is the average thickness of the specimen before test (mm), and D is the average diameter of the specimen before testing (mm). Conforming to ASTM E303, the British pendulum number was measured on the surface of the 50 × 50 × 5 cm specimen with a calibrated British pendulum anti-sliding tester (EL42-6000, ELE International, Leighton Buzzard, UK). In the end, the averaged value was regarded as the BPN value, which represents the anti-skid ability.

Indoor Temperature Measurement under Radiation Intensity of 623 W/m 2
The indoor temperature measurement of the specimens was divided into a heating period under a radiation intensity of 623 W/m 2 and cooling period without radiation intensity. The temperature-time relationships and temperature profiles of specimens were drawn and discussed.

Heating Period
The specimens were heated constantly for 24 h under the radiation intensity of 623 W/m 2 . Figure 6 shows the temperature-time relationships in the heating period under the radiation intensity of 623 W/m 2 . As seen from Figure 6, the specimens reached thermal equilibrium at 24 h. In Figure 6a, the highest temperature of 85 • C occurred in BOF-75 at the end of 24 heating hours and the temperature decreased in the order of BOF-55 and then BOF-45, while the benchmark was 70 • C. No obvious order was observed at the depths of 1 cm and 2 cm, as shown in Figure 6b,c. As seen in Figure 6d, BOF-75 possessed the lowest temperature at the depth of 3 cm. The indirect tensile strength was measured according to ASTM D6931 [28]. The specimen was attached between two load stripes and was loaded radially at a speed of 50 ± 5 mm/min. The width of stripes was 12.7 ± 0.3 mm, complying with the standard. The specimen was exposed to 25 °C water for 30-120 min, and the maximum load at fracture was measured. The indirect tensile strength can be calculated from Equation (1), as follows.
In Equation (1), St is the indirect tensile strength (kPa), P is the maximum load (N), t is the average thickness of the specimen before test (mm), and D is the average diameter of the specimen before testing (mm). Conforming to ASTM E303, the British pendulum number was measured on the surface of the 50 × 50 × 5 cm specimen with a calibrated British pendulum anti-sliding tester (EL42-6000, ELE International, Leighton Buzzard, UK). In the end, the averaged value was regarded as the BPN value, which represents the anti-skid ability.

Indoor Temperature Measurement under Radiation Intensity of 623 W/m 2
The indoor temperature measurement of the specimens was divided into a heating period under a radiation intensity of 623 W/m 2 and cooling period without radiation intensity. The temperature-time relationships and temperature profiles of specimens were drawn and discussed.

Heating Period
The specimens were heated constantly for 24 h under the radiation intensity of 623 W/m 2 . Figure 6 shows the temperature-time relationships in the heating period under the radiation intensity of 623 W/m 2 . As seen from Figure 6, the specimens reached thermal equilibrium at 24 h. In Figure 6a, the highest temperature of 85 °C occurred in BOF-75 at the end of 24 heating hours and the temperature decreased in the order of BOF-55 and then BOF-45, while the benchmark was 70 °C. No obvious order was observed at the depths of 1 cm and 2 cm, as shown in Figure 6b There were two heat sources: one was the halogen and infrared lamps, and the other was the accumulated heat underneath the asphalt concrete; both substantially affected the temperature at the depths of 1 cm and 2 cm.
The order of the temperature on the surface was totally in reverse compared to the temperature at the depths of 3 cm and deeper. The temperatures in profile taken from the thermal equilibrium condition at the end of 24 heating hours are shown in Table 4 and Figure 7. The depth-temperature curves of BOF-55 and BOF-75 were similar. In the profile, their temperatures on the surface were higher than inside and broadly decreased when going deeper. On the other hand, the depth-temperature curve of the benchmark and BOF-45 were similar in another way: the highest temperature was recorded at the depth of 3 cm. Compared to the benchmark specimen, although all the specimens with BOFS replacements had hotter surfaces, their temperatures were lower starting at 3 cm depths, especially BOF-75. There were two heat sources: one was the halogen and infrared lamps, and the other was the accumulated heat underneath the asphalt concrete; both substantially affected the temperature at the depths of 1 cm and 2 cm.
The order of the temperature on the surface was totally in reverse compared to the temperature at the depths of 3 cm and deeper. The temperatures in profile taken from the thermal equilibrium condition at the end of 24 heating hours are shown in Table 4 and Figure 7. The depth-temperature curves of BOF-55 and BOF-75 were similar. In the profile, their temperatures on the surface were higher than inside and broadly decreased when going deeper. On the other hand, the depth-temperature curve of the benchmark and BOF-45 were similar in another way: the highest temperature was recorded at the depth of 3 cm. Compared to the benchmark specimen, although all the specimens with BOFS replacements had hotter surfaces, their temperatures were lower starting at 3 cm depths, especially BOF-75.

Cooling Period
All the lamps were turned off after 24 h of heating; the temperature change subsequent 1-h cooling was recorded and shown in Figure 8. When the heat sour removed, the heat obviously dissipated from the specimen. The cooling rate at each gradually became similar, and the temperature-time relationships of all specimens depth of 3 cm were the most concentrated.

Cooling Period
All the lamps were turned off after 24 h of heating; the temperature change in the subsequent 1-h cooling was recorded and shown in Figure 8. When the heat source was removed, the heat obviously dissipated from the specimen. The cooling rate at each depth gradually became similar, and the temperature-time relationships of all specimens at the depth of 3 cm were the most concentrated.
The depth-temperature profiles of specimens within five cooling hours in the cooling period are shown in Figure 9. In the beginning of the cooling period, the sequences of the surface temperatures are the same. Notably, all the specimens containing BOFS possessed a higher temperature at a depth of 5 cm than 3 cm. As seen from Figure 9, the temperature of the benchmark specimen at the depth of 3 cm was higher than the temperature at 5 cm in the first hour, but the temperature at a depth of 3 cm was less than the temperature at a depth of 5 cm at other times.

Indoor Temperature Measurement under Radiation Intensity of 436 W/m 2
The indoor temperature measurement of the specimens was divided into a heating period under a radiation intensity of 436 W/m 2 and a cooling period without radiation intensity. The temperature-time relationships and temperature profiles of specimens were determined and discussed.
All the lamps were turned off after 24 h of heating; the temperature change in the subsequent 1-h cooling was recorded and shown in Figure 8. When the heat source was removed, the heat obviously dissipated from the specimen. The cooling rate at each depth gradually became similar, and the temperature-time relationships of all specimens at the depth of 3 cm were the most concentrated. The depth-temperature profiles of specimens within five cooling hours in the cooling period are shown in Figure 9. In the beginning of the cooling period, the sequences of the surface temperatures are the same. Notably, all the specimens containing BOFS possessed a higher temperature at a depth of 5 cm than 3 cm. As seen from Figure 9, the temperature of the benchmark specimen at the depth of 3 cm was higher than the temperature at 5 cm in the first hour, but the temperature at a depth of 3 cm was less than the temperature at a depth of 5 cm at other times.

Heating Period
The specimens were heated constantly for 24 h under the radiation intensity of 436 W/m 2 until the temperatures of the specimens reached a stable state. As shown in Figure 10a, BOF-75 had the highest temperature of 66.82 • C and BOF-55 had the next highest temperature of 64.74 • C, corresponding with the results under the radiation intensity of 623 W/m 2 . Although the order between BOF-45 and the benchmark was disrupted and did not correspond with the previous results, the temperature difference between 56.95 • C and 57.69 • C was not significant. As shown in Figure 10b,c, the temperature changes at the depths of 1 cm and 2 cm were similar to the results under the radiation intensity of 623 W/m 2 . Influenced by the upper and lower heat sources, it was difficult to find the order at these two depths. As shown in Figure 10d, the temperatures at the depth of 3 cm were in descending order of BOF-75 (63.46 • C), BOF-55 (63.53 • C), BOF-45 (64.17 • C), and then the benchmark (65.65 • C). This order happened for the results under the radiation intensity of 623 W/m 2 , which was in reverse on the surface. As shown in Figure 10e,f, the temperature-time relationship at the depths of 4 cm and 5 cm were similar to those in Figure 10d. cm depth; (e) 4 cm depth; and (f) 5 cm depth. The depth-temperature profiles of specimens within five cooling hours in th period are shown in Figure 9. In the beginning of the cooling period, the sequen surface temperatures are the same. Notably, all the specimens containing BOFS p a higher temperature at a depth of 5 cm than 3 cm. As seen from Figure 9, the tem of the benchmark specimen at the depth of 3 cm was higher than the temperatu in the first hour, but the temperature at a depth of 3 cm was less than the temp a depth of 5 cm at other times.

Indoor Temperature Measurement under Radiation Intensity of 436 W/m 2
The indoor temperature measurement of the specimens was divided into a heating period under a radiation intensity of 436 W/m 2 and a cooling period without radiation intensity. The temperature-time relationships and temperature profiles of specimens were determined and discussed.

Heating Period
The specimens were heated constantly for 24 h under the radiation intensity of 436 W/m 2 until the temperatures of the specimens reached a stable state. As shown in Figure  10a, BOF-75 had the highest temperature of 66.82 °C and BOF-55 had the next highest temperature of 64.74 °C , corresponding with the results under the radiation intensity of 623 W/m 2 . Although the order between BOF-45 and the benchmark was disrupted and did not correspond with the previous results, the temperature difference between 56.95 °C and 57.69 °C was not significant. As shown in Figure 10b,c, the temperature changes at the depths of 1 cm and 2 cm were similar to the results under the radiation intensity of 623 W/m 2 . Influenced by the upper and lower heat sources, it was difficult to find the order at these two depths. As shown in Figure 10d, the temperatures at the depth of 3 cm were in descending order of BOF-75 (63.46 °C ), BOF-55 (63.53 °C ), BOF-45 (64.17 °C ), and then the benchmark (65.65 °C ). This order happened for the results under the radiation intensity of 623 W/m 2 , which was in reverse on the surface. As shown in Figure 10e,f, the temperature-time relationship at the depths of 4 cm and 5 cm were similar to those in Figure 10d. The depth-temperature profile taken from the balanced condition at the end of 24 heating hours is shown in Table 5 and Figure 11. In the profile, the temperatures of all specimens declined constantly from 3 cm to 5 cm depths. Furthermore, it clearly shows that BOF-75 possessed the lowest temperature of 63.07 °C at the depth of 4 cm and 60.04 °C at the depth of 5 cm. Furthermore, the depth-temperature curves of BOF-55 and BOF-75 were similar. In the profile, their temperatures at the depth of 3 cm were higher than at the depth of 5 cm, and the temperatures on the surface were higher than that at the depth of 3 cm. On the other hand, the depth-temperature curves of the benchmark and BOF-45 were similar in a different way. These two specimens had the highest temperatures at the depth of 3 cm and the lowest temperatures on the surface. All the trends above were similar to the results under the radiation intensity of 623 W/m 2 .  The depth-temperature profile taken from the balanced condition at the end of 24 heating hours is shown in Table 5 and Figure 11. In the profile, the temperatures of all specimens declined constantly from 3 cm to 5 cm depths. Furthermore, it clearly shows that BOF-75 possessed the lowest temperature of 63.07 • C at the depth of 4 cm and 60.04 • C at the depth of 5 cm. Furthermore, the depth-temperature curves of BOF-55 and BOF-75 were similar. In the profile, their temperatures at the depth of 3 cm were higher than at the depth of 5 cm, and the temperatures on the surface were higher than that at the depth of 3 cm. On the other hand, the depth-temperature curves of the benchmark and BOF-45 were similar in a different way. These two specimens had the highest temperatures at the depth of 3 cm and the lowest temperatures on the surface. All the trends above were similar to the results under the radiation intensity of 623 W/m 2 .

Cooling Period
After 24-h heating, the temperature change in the subsequent 1-h cooling was recorded, as shown in Figure 12. When the heat source was removed, the heat obviously dissipated from the specimens, especially those which possessed higher conductivity. BOF-75 had the lowest thermal conductivity, leading to the slowest cooling rate. In contrast, the benchmark specimen had the highest thermal conductivity, leading to the fastest cooling rate. Also, the cooling rates for each depth gradually became similar, and the rates at the depth of 3 cm were the most concentrated.

Cooling Period
After 24-h heating, the temperature change in the subsequent 1-h cooling was recorded, as shown in Figure 12. When the heat source was removed, the heat obviously dissipated from the specimens, especially those which possessed higher conductivity. BOF-75 had the lowest thermal conductivity, leading to the slowest cooling rate. In contrast, the benchmark specimen had the highest thermal conductivity, leading to the fastest cooling rate. Also, the cooling rates for each depth gradually became similar, and the rates at the depth of 3 cm were the most concentrated.
The depth-temperature profiles of specimens within five cooling hours can be seen in Figure 13. The temperature difference between 3 cm and 5 cm depths deserves discussion. At the end of five hours, BOF-75 had the lowest temperature on the surface, and its temperature at the depth of 3 cm was lower than that at 5 cm. With regard to BOF-55, its temperature at the depth of 3 cm was lower than that at 5 cm only in the first hour, with a slight difference, and its temperature on the surface was lower than that at the depth of 3 cm all the time. Also, the temperature of the benchmark at the depth of 3 cm was higher than that at 5 cm until the fourth hour.
To sum up, compared to the benchmark asphalt concrete, no matter which radiation intensity was taken, the specimens with partial BOFS replacement showed faster cooling rates when the light was turned off, which might refer to more heat accumulating inside the body. orded, as shown in Figure 12. When the heat source was removed, the heat obviously dissipated from the specimens, especially those which possessed higher conductivity. BOF-75 had the lowest thermal conductivity, leading to the slowest cooling rate. In contrast, the benchmark specimen had the highest thermal conductivity, leading to the fastest cooling rate. Also, the cooling rates for each depth gradually became similar, and the rates at the depth of 3 cm were the most concentrated. The depth-temperature profiles of specimens within five cooling hours can be seen in Figure 13. The temperature difference between 3 cm and 5 cm depths deserves discussion. At the end of five hours, BOF-75 had the lowest temperature on the surface, and its temperature at the depth of 3 cm was lower than that at 5 cm. With regard to BOF-55, its temperature at the depth of 3 cm was lower than that at 5 cm only in the first hour, with a slight difference, and its temperature on the surface was lower than that at the depth of 3 cm all the time. Also, the temperature of the benchmark at the depth of 3 cm was higher than that at 5 cm until the fourth hour. The depth-temperature profiles of specimens within five cooling hours can be seen in Figure 13. The temperature difference between 3 cm and 5 cm depths deserves discussion. At the end of five hours, BOF-75 had the lowest temperature on the surface, and its temperature at the depth of 3 cm was lower than that at 5 cm. With regard to BOF-55, its temperature at the depth of 3 cm was lower than that at 5 cm only in the first hour, with a slight difference, and its temperature on the surface was lower than that at the depth of 3 cm all the time. Also, the temperature of the benchmark at the depth of 3 cm was higher than that at 5 cm until the fourth hour. To sum up, compared to the benchmark asphalt concrete, no matter which radiation intensity was taken, the specimens with partial BOFS replacement showed faster cooling

Outdoor Temperature Measurement
The temperature-time relationships of specimens in the outdoor test, across both daytime and nighttime for a full 24 h, were measured and are shown in Figure 14. The specimens' temperature and environmental information were recorded and are shown in Table 6. The field test results were easily affected by the environmental factors, especially after sunset. Therefore, it is more reliable to use the average ambient temperature for comparison. As shown in Table 6, during the daytime, the benchmark had a higher temperature than the ambient temperature of about 35 • C. Also, unlike the indoor test, specimens could not obtain long-term and sufficient heat from the sun. As a consequence, for BOF-75, its low thermal conductivity and large heat capacity resulted in relatively low temperatures. The calculation of the heat capacity is explained in Section 4.4.
As shown in Table 6, using the ambient temperature of about 30 • C during nighttime as the comparison basis, the temperatures of the benchmark specimen were generally higher than those of BOF-55 and BOF-75. In terms of results, the use of BOFS did have a positive effect on pavement cooling. rates when the light was turned off, which might refer to more heat accumulating inside the body.

Outdoor Temperature Measurement
The temperature-time relationships of specimens in the outdoor test, across both daytime and nighttime for a full 24 h, were measured and are shown in Figure 14. The specimens' temperature and environmental information were recorded and are shown in Table 6. The field test results were easily affected by the environmental factors, especially after sunset. Therefore, it is more reliable to use the average ambient temperature for comparison. As shown in Table 6, during the daytime, the benchmark had a higher temperature than the ambient temperature of about 35 °C. Also, unlike the indoor test, specimens could not obtain long-term and sufficient heat from the sun. As a consequence, for BOF-75, its low thermal conductivity and large heat capacity resulted in relatively low temperatures. The calculation of the heat capacity is explained in Section 4.4.

Heat Capacity, Newton's Cooling and Radiation Cooling
To confirm the radiation cooling effectiveness of specimens, the loss by radiation cooling was quantified based on several assumptions as follows. The heat which enters the asphalt concrete and then transfers to the ground has only a small impact on pedestrian and human activities; the heat below the specimens was not discussed. On the other hand, the heat on the upper surface of the asphalt concrete specimens is worth further study.
Ideally, when objects receive heat, they transfer it in three ways: reflection, transmission, and absorption. The sum of the energy from these ways must be equal to the energy incident, an illustration of which is shown in Figure 15a. However, it is hard to determine the reflection, transmission, and absorption of specimens realistically, even using a pyranometer. It is impossible to know how much heat has dissipated through radiation. We can only assume that, following the law of conversation of energy, the relationship between the energy entering and leaving the specimen is constant. In this study, the focus was on the cooling period, as shown in Figure 15b. As the temperature of the test body cooled over time under the radiation intensity of 623 W/m 2 , the time interval between data was set to 5 min.

Specimen
Benchmark BOF-45 BOF-55 BOF As shown in Table 6, using the ambient temperature of about 30 °C during nigh as the comparison basis, the temperatures of the benchmark specimen were gen higher than those of BOF-55 and BOF-75. In terms of results, the use of BOFS did h positive effect on pavement cooling.

Heat Capacity, Newton's Cooling and Radiation Cooling
To confirm the radiation cooling effectiveness of specimens, the loss by rad cooling was quantified based on several assumptions as follows. The heat which the asphalt concrete and then transfers to the ground has only a small impact on p trian and human activities; the heat below the specimens was not discussed. On the hand, the heat on the upper surface of the asphalt concrete specimens is worth f study.
Ideally, when objects receive heat, they transfer it in three ways: reflection, tran sion, and absorption. The sum of the energy from these ways must be equal to the e incident, an illustration of which is shown in Figure 15a. However, it is hard to dete the reflection, transmission, and absorption of specimens realistically, even using a nometer. It is impossible to know how much heat has dissipated through radiatio can only assume that, following the law of conversation of energy, the relationsh tween the energy entering and leaving the specimen is constant. In this study, the was on the cooling period, as shown in Figure 15b. As the temperature of the test cooled over time under the radiation intensity of 623 W/m 2 , the time interval betwee was set to 5 min. The heat stored in a specimen at various time intervals can be calculated as in tion (2). In this function, Q is the heat energy in storage (J), m is the mass of the spe The heat stored in a specimen at various time intervals can be calculated as in Equation (2). In this function, Q is the heat energy in storage (J), m is the mass of the specimen (kg), c is the heat capacity of the specimen (J/kg-K), and ∆T is the change in temperature (K).
The specimen was divided into five equal parts each with a thickness of 5 cm; the sectional drawing of the specimen is shown in Figure 16. In each part, the temperature at the center was the average of the adjoining temperature at the upper and lower part. Then, the sum of the heat stored in each part at the t minute in the cooling period represented the heat storage in the entire specimen at that time. Finally, the heat storage of (t + 5) minutes was subtracted from that at the t minute, and added up over 24 h to obtain the final heat storage of the specimen. This formula is shown in Equation (3), and the calculation results for each hour are presented in Table 7. (kg), c is the heat capacity of the specimen (J/kg-K), and ΔT is the change in temperature (K).
The specimen was divided into five equal parts each with a thickness of 5 cm; the sectional drawing of the specimen is shown in Figure 16. In each part, the temperature at the center was the average of the adjoining temperature at the upper and lower part. Then, the sum of the heat stored in each part at the t minute in the cooling period represented the heat storage in the entire specimen at that time. Finally, the heat storage of (t + 5) minutes was subtracted from that at the t minute, and added up over 24 h to obtain the final heat storage of the specimen. This formula is shown in Equation (3), and the calculation results for each hour are presented in Table 7.
Based on Newton's law of cooling, the heat stored in the pavement is also transferred into the atmosphere, which is called "Newton's cooling" in this study, as shown in Equation (4): In this function, q is the power of Newton's cooling transferred out of the specimen (watts), h is the heat transfer coefficient (W/m 2 -K), A is the heat transfer surface area (m 2 ), and ΔT is the temperature difference between the environment and specimen (K). The final value for the difference, resulting from Newton's cooling, between specimen and air can be calculated as in Equation (5). Notably, it was necessary to multiply the difference  Based on Newton's law of cooling, the heat stored in the pavement is also transferred into the atmosphere, which is called "Newton's cooling" in this study, as shown in Equation (4): In this function, q is the power of Newton's cooling transferred out of the specimen (watts), h is the heat transfer coefficient (W/m 2 -K), A is the heat transfer surface area (m 2 ), and ∆T is the temperature difference between the environment and specimen (K). The final value for the difference, resulting from Newton's cooling, between specimen and air can be calculated as in Equation (5). Notably, it was necessary to multiply the difference between the total heat of t and that at (t + 5) minutes by a time interval of 5 min. The calculation results for each hour are presented in Table 8.
The final Newton's cooling in 24 h = 288 ∑ t=0 (q 5t ) − q 5(t+1) × 5 Finally, following the above data in Tables 7 and 8, Newton's cooling energy subtracted from the stored heat energy equals the radiation cooling energy, as shown in Table 9. The overall amounts of the accumulated heat storage energy, Newton's cooling energy, and the radiation cooling energy are presented in Tables 10-12, respectively. As shown in Table 10, heat storage energy within 7 h accounted for up to 80% of heat storage energy within 24 h. We can deduce that the main cooling was completed within 7 h in the cooling period. Also, BOF-75 had the highest heat storage energy among all the specimens, as shown in Table 10. In the seventh cumulated hour, Newton's cooling energy accounted for about 90% of the stored energy, as shown in Table 11. Therefore, the main way of dissipating heat after 7 h was radiation cooling. In the above calculations, it seems that, under the same radiation intensity, the specimens with higher heat capacities could absorb more heat inside the specimen, and less heat escaped to the environment. And according to Newton's cooling law, the thermal equilibrium occurs under the condition of releasing sufficient energy for the specimen. BOF-75 had a higher surface temperature because it stored as well as absorbed more heat, which is a reasonable result. As far as radiation cooling is concerned, a larger value is considered better, which means that more heat leaves the specimen in this way. BOF-75 had the largest radiation cooling among the specimens.
Based on the results, it can be concluded that BOF-75 had better heat capacity and could dissipate less heat into the environment and absorb more heat than other specimens. Although a higher surface temperature occurred under the same radiation intensity, BOF-75 possessed the best radiation cooling ability in the cooling period among the specimens, and the heat could be transferred into long-wave radiation, which is not easily absorbed by air. In addition, the material emissivity results measured by the Fourier Infrared Spectrometer (FTIR) help to show that BOF-75 demonstrated a better performance in the radiation cooling project.

Mechanical Test Results
The stability values of the specimens were measured according to ASTM D6927 [27]. As shown in Table 13 and Figure 17, the stability values of all specimens were far larger than the standard requirement of 8.006 kN, which was set by the Asphalt Institute. As the replacement ratio of BOFS increased, the stability value increased. Among all specimens, BOF-75 had the highest average stability value of 34.54 kN. The stability values of the specimens were measured according to ASTM D6927 [27]. As shown in Table 13 and Figure 17, the stability values of all specimens were far larger than the standard requirement of 8.006 kN, which was set by the Asphalt Institute. As the replacement ratio of BOFS increased, the stability value increased. Among all specimens, BOF-75 had the highest average stability value of 34.54 kN.  Indirect tensile strength, one of the indexes of asphalt concrete, can be inferred from the resistance of rut and crack. Following ASTM D6931 [28], the indirect tensile strength was measured and calculated as shown in Table 14 and Figure 18. It implied that BOFS had a better binding ability with bitumen than natural aggregate did.  Indirect tensile strength, one of the indexes of asphalt concrete, can be inferred from the resistance of rut and crack. Following ASTM D6931 [28], the indirect tensile strength was measured and calculated as shown in Table 14 and Figure 18. It implied that BOFS had a better binding ability with bitumen than natural aggregate did.  The BPN was measured according to ASTM E303 [29]. As shown in Table 15, all specimens conformed to the BPN of 45 [30]. Specimens replaced partially by BOFS demonstrated an anti-skid ability, while the flatness was roughly the same. There is no doubt that using asphalt concrete pavement in which natural aggregates have been partially replaced by BOFS affects driving safety.

5.
The anti-skid index, the BPN value, was far larger than the standard of 45; the stability value complied with the standard of 8.006 kN, and BOF-75 could reach 34.54 kN. Consequently, BOFS can be successfully applied in pavements. 6.
All asphalt concrete with partial BOFS replacement had a greater indirect tensile strength than the traditional concrete. It is very likely that BOFS has a better binding ability with bitumen than natural aggregate does. 7.
From the test results, it can be seen that BOFS has great applicability in pavements due to its thermal performance and mechanical properties. Institutional Review Board Statement: Not applicable.