Radiosity from Individual Urban Landscape Elements Measured Using a Modified Low-Cost Temperature Sensor

: Loss of green space in our suburban environment is contributing to increased urban heat. The material properties of surface treatments or landscape elements (LEs) are a determining factor in the amount, timing, and type of radiation present in the local environment. Landscape designers can use this information to better design for urban heat management, as emitted and reflected radiation (radiosity) from LEs can affect pedestrians via heat stress and glare and affect energy usage in buildings and houses if the landscape sky view factor is low. Low-cost black painted iButton temperature sensors were successfully used as radiometers to concurrently measure the daytime radiosity from 19 LEs samples located on an oval in the warm temperate climate of Perth, (Australia). Normalisation against gloss white paint on polystyrene removed the effect of varying weather conditions. Each LE had the same normalised average radiosity ( DR av ) between seasons (within ±5%), meaning the relative radiosity of new LEs can be measured on any day. White and lighter coloured LEs had the highest DR av and would have the most detrimental effect on nearby objects. Plants and moist LEs had the least DR av and would be most beneficial for managing local daytime urban heat. Measuring relative radiosity with iButtons presents a new way to examine the effect of LEs on the urban environment.


Introduction
The urban heat island (UHI) effect is the increased night−time ambient temperatures in urban areas when compared with nearby rural areas [1]. Superimposed on this are the increasing ambient temperatures due to global warming, causing particular concern in urban areas. One of the main drivers of the UHI is land use change and the replacement of vegetation with hardscape such as concrete, asphalt and pavers [2][3][4][5]. The landscapes surrounding buildings and houses in the urban and suburban environment may thus play an important role in subsequently managing this accumulated heat.
Landscapes can consist of a number of different individual surface treatments defined here as landscape elements (LEs), such as turf grass, concrete paving, crushed rock, mulch, trees, and shrubs. The amount and timing of urban heat is influenced by the material properties of these LEs, in particular albedo, moisture availability, surface roughness, and thermal inertia [1,6]. On reaching the surface of the LE, a portion of the solar radiation is reflected (albedo) and the remainder is absorbed. The surface temperature rises until there is a balance between heat coming in and heat which is convected, conducted (thermal inertia), emitted and evaporated (if moisture is present) (Figure 1). These properties affect urban heat through the parameters of surface temperature, sensible heat flux (convective heat), and total radiant energy (radiosity). Each of these properties and parameters are relevant to urban heat in different ways.
The albedo of a surface is the fraction of solar radiation which is reflected. White and/or reflective materials generally have lower surface temperatures due to greater reflection of radiation [7][8][9]. However, studies have shown that increased reflected radiation can significantly affect human thermal comfort within that landscape [4,[9][10][11][12].
The presence of water at the surface of a LE will lead to some thermal energy being lost through evaporation. Wet permeable pavement was studied by Li [13] and the increased rate of evaporation was found to decrease surface temperatures. The caveat to this was that evaporation rates were only high whilst water was available close to the surface. The surface temperature of plants has been found to be lower than other materials by many researchers [14][15][16][17], and this has been associated with higher rates of evapotranspiration [18,19]. Water availability is crucial for this process to occur properly, which may pose problems in drier climates.
The heat storage of a LE is affected by both its heat capacity and thermal conductivity, together called the thermal admittance or thermal inertia [1]. The amount of heat absorbed is a function of its coefficient of thermal conductivity and the temperature difference between the surface and subsurface layer [20]. Higher thermal inertia of pavements has been found to reduce daytime heat through lower surface temperatures [4], but found to increase heat at night−time [6]. Hence, as well as being a measure of how uncomfortable a surface may be for pedestrians, surface temperature also affects the amount of heat available for storage in the LE. Many studies of urban LE surface temperatures have been undertaken [4,13,15,[21][22][23][24][25][26][27][28]. In general, these studies found that manufactured, low albedo and dry materials were hotter than plants or than LEs that were moist or had higher albedos.
Surface temperature also directly affects the amount of heat convected into the air from that surface. From Newton's law of cooling for convection, the amount of convected heat is proportional to the difference between the surface and ambient temperatures, and is also proportional to the coefficient of convection for that surface [29]. The coefficient of convection comprises natural and forced components. When wind speed increases, the convection increases due to the forced convection component. Large areas of a LE with a high surface temperature in urban spaces thus have the potential to significantly increase urban heat.
Longwave radiation (LWR) is emitted from all bodies at temperatures above absolute zero. The higher the surface temperature of a LE, the greater the amount of emitted LWR through the application of Stefan-Boltzmann's law [20]. The amount of emitted LWR is also dependent on the emissivity (ε) of the body, a measure of how effective that body is at emitting heat (0 for a perfect reflector, 1 for a perfect emitter). The emissivity of most LEs is very high, ranging from ~0.76 to ~0.96, with the exception of shiny metals, which can be as low as 0.04 [30], and it varies depending on material surface properties [31].
The radiosity from a LE is the combination of this emitted LWR as well as the reflected solar radiation ( Figure 1). Whilst radiosity from LEs does not directly heat the air, it can increase urban heat by striking nearby objects, causing them to increase in temperature [12]. The radiosity of LEs has previously been measured using commercial net radiometers over in situ LEs [12,32], or over specifically constructed LE samples [33]. Other researchers have used pyranometers in conjunction with surface temperature sensors to effectively determine the radiosity over constructed samples [15,34]. As these instruments are relatively expensive, measurements were usually made of one or two LEs at a time, by either alternately facing the instrument towards the LE and then to the sky in the case of a pyranometer, or by leaving the net radiometer above one or two LEs for a period of time, before moving on to the next. This process drastically reduces the number of measurable LEs. Furthermore, the radiosity of smaller LEs cannot be measured using these instruments, as with a diameter of ~ 0.1 m, the shadow cast by the instrument becomes significant, especially as the height above the LE needs to be lowered to increase the amount of radiation received from the LE compared with the surroundings.
Diurnal variations in radiosity are important, as the thermal properties of LEs lead to energy being released back to the surroundings at different times. The timing of this energy release could impact how and when landscapes are used by people. Seasonal variation of radiosity could also be important but does not appear to have been well studied. Li [13], however, measured the albedo and surface temperature of concrete, asphalt and pavers across seasons. Whilst these factors were not combined into a single radiosity measurement, the surface temperature was found to be higher in summer, but the albedo did not change with season.
To address the glare of high albedo materials whilst retaining the benefits of cooler surface temperatures, new materials named cool-coloured materials have been engineered. These increase the reflectance of solar radiation in the near-infrared (NIR) band, offering non-white colours which reflect a significant amount of solar radiation. A laboratory study by Xie et al. [35] found that increased NIR reflectance reduced surface temperatures, while Kyriakodis and Santamouris [36] found the application of these materials to roofs, roads and pavements also reduced surface temperatures when compared with conventional materials (although cool asphalt reflectance halved after deposition of tyre rubber over time). However, although reduced glare for pedestrians is a major benefit of these materials, the effect of the reflected NIR on people and buildings has not been investigated. Hence, as new materials are developed, radiosity testing will be required in the laboratory, but more importantly in the field, where other objects and weathering will influence their radiosity.
Measurement of the radiosity from different LEs will help identify those which would have the most impact on heat management. To reduce the influence of fluctuations in weather conditions on the relative radiosity of different LEs in the field, they should be measured concurrently. Taking a thermal image to find surface temperature is a fairly simple process; however, measuring, and then re-locating, a single radiometer above each LE for each measurement is neither practical, nor time efficient. Furthermore, the high cost of a radiometer and its associated data logger is a barrier to the purchase of multiple units. Unlike other studies of radiosity, a reference material of white painted polystyrene was used in this research to take into account the effects of changing weather conditions (ambient temperature, wind speed and sky temperature).
Akbari et al. [37] developed a standardised procedure for ranking roofing materials according to both their emitted heat and their solar reflectivity, called the Solar Reflectance Index (SRI). It is generated by measuring a material under two controlled laboratory measurement procedures of reflectance [38] and emittance [39], and calculates SRI based on one assumed set of standard ambient conditions and a convective coefficient for the material. This method could be used for LEs; however, it does not take into account diurnally and seasonally varying ambient temperatures and solar radiation. The assumed convective coefficient would not be applicable to more massive or rough surfaces as are found in LEs as opposed to roofing materials. However, if the reflected and emitted radiant energy from LEs could be measured concurrently and diurnally using a single instrument in the field, then this may provide a simpler and more efficient way of assessing LEs under real world conditions. Aoki and Mizutani [40] used a low-cost instrument to investigate the vertical distribution of reflected SWR from three different LEs. They constructed a polystyrene post consisting of numerous thin matte black painted aluminium plates mounted at set intervals and monitored the plates' surface temperatures using thermocouples. A similar principle to this has been explored in this paper by using small modified temperature sensors or thermochrons called iButtons [41] to detect radiation across both the visible and thermal spectra. The purpose of this paper was to investigate a low-cost (~AU$50) temperature sensor for use as a radiometer, and to use multiple units concurrently to measure the radiosity from a number of different LEs typically used in Perth urban and suburban areas. Small LEs were tested across the four seasons, whilst larger, ground coupled LEs were tested in summer. LEs were ranked in order of summer radiosity and the implications of the LEs material properties on this ranking are discussed.

Materials and Methods
Obtaining sufficient measurement accuracy in outdoor urban environments is challenging due to variable weather conditions and interactions between intended measurement materials and surrounding objects such as buildings, vehicles and pedestrians. However, the benefit of measuring in the outdoor environment is the ability to capture the response of the tested LEs to the natural varying diurnal solar radiation (SR) across seasons. Using samples of LEs is beneficial, as they can be placed in an open area with a high sky view factor (SVF) [42] in order to minimise the radiative interference from surrounding objects. The LEs would be exposed to full SR throughout the day. In this study, measurements were performed on samples of LEs located outdoors under a high SVF.
This study consisted of two phases. The first study (phase 1) tested LEs which were smaller in size and placed directly on to the turf grass on the oval. These were measured across each of the four seasons. The second study (phase 2) used larger LEs which were coupled to the ground through a bed of sand placed on the turf grass. Due to their size and subsequent increase in cost, only a subset of the phase 1 LEs were used in phase 2. These were measured in summer only. Loveday et al. [28] indicated that the size and coupling of the LEs to the ground affected the quantity of thermal radiation. Hence, the purpose of phase 2 was to investigate the effect of these factors on radiosity. LEs were compared with each other to determine a relative radiosity ranking.

Materials
Aluminium foil, white painted polystyrene (WPP), black painted cardboard on polystyrene and 19 different LEs (see Table 1) were positioned in the centre of a turf grassed oval at the Booragoon Primary School in Perth, Western Australia (coordinates: −32.035918°, 115.826655°). Perth has a temperate climate with a dry summer, classified as Csa under the Köppen-Geiger climate classification system [43]. The foil was first crumpled and then straightened out (but not smoothed) and mounted shiny side up on a flat wooden board. This provided a diffuse reflector with a reflectivity as close to 1.0 as possible and was used to measure reflected sky temperature (as per the reflector method used in [44]). The black and white painted polystyrene were used as reference materials, with the polystyrene insulating these from the ground. Most LEs were placed directly onto the turf grass ( Figure 2). These were the red, sandstone, and grey pavers, both concrete slabs, artificial turf grass, all of the shade cloths, the decking, and the limestone block. Some LEs were placed onto a white polypropylene (PP) bag on the turf grass. These were the white and polished black stones, and the moist soil. The remaining LEs (soil, asphalt, mulch, and white and grey sand), were also placed onto a PP bag but were then placed within a low-sided plastic crate which allowed the samples to be transported between seasons whilst retaining the same surface area and depth of the LE. The artificial grass was set up in two ways to observe if there were any differences related to the installation method. One sample of artificial grass was placed directly onto the oval grass, while a second sample was first laid onto a bed of crushed rock and the surface sprinkled with white sand, as per a typical installation procedure, and then placed on to a PP bag within a plastic crate. To prevent them from blowing away, all light mass LEs (artificial grass, shade cloths, WPP and BPP), as well as the edges of any PP (including the PP in the plastic crates), were pinned to the ground using four small metal tent pegs. The pegs were hammered straight into the ground next to the LEs, with the short end of the peg unobtrusively holding down each of the four corners of the LE or PP. The turf grass on the oval was Pennisetum clandestinum, (the common name being kikuyu), and was managed and irrigated. LEs were spaced between 200 and 600 mm apart to reduce their effect on each other. To prevent any interference from the public, safety tape was erected around the site.
Ideally, all the LEs would have the same surface area in order to minimise confounding factors. In reality, both the standard sizes of the each LE type, and the way they were constructed for the test work, varied ( Table 1). The different shapes were accounted for by calculation of view factors, described in Section 2.3.4.
Small wireless temperature sensors called Themochron iButtons (Maxim Integrated, San Jose, CA, USA) [42], capable of logging data over long time periods, were used as low-cost radiometers by painting one surface with matt black paint [45]. These painted sensors measure 17 mm in diameter and are 6 mm high (Figure 3a). The undercoat was a spray of matt black (brand Fiddly Bits), with a topcoat of Rust Guard satin black spray (brand White Knight), applied to give a visibly darker finish. The types of iButtons used were the DS1921G low temperature (−40 to +85 °C, accuracy ±1.0 °C), and the DS1922T, high temperature (0 to +125 °C, accuracy ±0.5 °C) thermochrons. The response time of the iButtons is given as 130 s, although validation studies performed by van Marken Lichtenbelt et al. [46], albeit on the DS1291H (range 15-45 °C), indicated a better accuracy and response time in practice than given by the manufacturer (on average 0.09 °C less than a standard thermometer, and a response of only 19 s). The black paint of the iButton is expected to absorb both longwave radiation (LWR) and short-wave radiation (SWR), with absorptivities of 0.88 [30] and 0.95 [47], respectively. It will also respond to convected heat from the LE, although this depends on its proximity to the LE and the strength of the wind currents.  Each LE had a black painted iButton mounted in a polystyrene block (described in [45]), suspended 0.2 m above it, protruding 0.15 m towards the centre of the LE from the middle of one edge of the LE (Figure 3b). The top side of the polystyrene was covered in aluminium foil to reduce solar heat gain to the iButton. A thermal imaging camera (Testo 876, Testo Pty. Ltd, Croydon South, Australia) with range of 8-14 µm and accuracy ±2 °C was used to take images from directly above each LE from a height of approximately 1 m.

Field Measurement Conditions
Measurements were performed over one day during each of the four seasons: 9th April 2017 (autumn); 24th August 2017 (winter); 5th October 2017 (spring); and 9th January 2018 (summer). Measurements commenced a few hours before sunrise and continued until a few hours after sunset. The chosen days had no visible cloud cover and had average wind speed for that time of year. Weather parameters logged every 600 s were accessed via a nearby weather station at Murdoch University [48], 3.6 km away. Conditions were close to the Bureau of Meteorology averages for the measured months. LEs were placed on the oval the day previous to each test day in order to equilibrate to outdoor conditions. The sampling rate of each iButton was set to 600 s by connecting it to a computer through a USB interface (LinkUSBi, Maxim Integrated, San Jose, CA, USA) and using the manufacturer's software (OneWireViewer, Maxim Integrated, San Jose, CA, USA). Thermal images were taken every half an hour.

Materials
The LEs used for phase 2 are described in Table 2. They were all the same product type as those used in the seasonal measurements, except for decking and artificial turf grass, where the original types were unavailable. Although the same product types were purchased for the other LEs, there were some dissimilarities in colour and size between some of the phase 1 and phase 2 LEs, particularly for the more natural materials such as crushed white rocks, pine bark mulch and the soils. The new artificial grass was a slightly darker colour, whilst the new decking was a different shade of brown ( Figure 4). In addition, an annual non-native ground cover of petunias (Petunias x hybrida) and some native saltbush seedlings (Atriplex nummularia) were tested as a comparison to turf grass. The ground cover (petunias) were in 100 mm pots closely packed together to give nearly 100% coverage of the sand. The saltbush seedlings were in 50 x 50 mm containers within seedling trays and similarly closely packed together. The planting densities are given in Table 2. The LEs used in this test were all 1200 mm x 1200 mm square, and each were separated by 1 m of managed and irrigated turf grass (kikuyu) ( Figure 5). As per phase 1, reference materials of WPP and black painted polystyrene (BPP) were included in these tests. Note that the BPP did not include cardboard as per phase 1, and, instead, the black paint was applied directly onto the polystyrene. Crumpled aluminium foil was again used to measure the reflected sky temperature. Ten of the LEs were laid onto a bed of a 1200 mm x 1200 mm square of yellow sand with a thickness of 50 mm, (seen around the edges of the plants in Figure 6). This sand provided a good thermal coupling between the LE and the ground, similar to how the LE would be installed in situ. The reasons for not placing the a) b) remaining LEs on the sand were as follows: the decking simulated typical decking, which is not directly coupled to the ground; the reference LEs (BPP and WPP) were required to be insulated from the turf grass; the soils were naturally effectively coupled to the ground through the infiltration of their particles through the turf grass to the soil below; and the turf grass LE was simply the existing turf grass on the oval.

Field Measurement Conditions
Phase 2 was performed on the same oval as the phase 1 tests. The LEs were set up on the oval on the 19th January 2019 to allow them to equilibrate with the conditions. The test day had clear skies only. Thermal images were taken every 30 min and iButton data were logged every 10 min from 06:00 to 21:00 on the 21st January 2019 (summer).

Energy Balance
The energy balance between the iButtons and the LEs was modelled to determine the unknown parameters of emissivity (described below) and coefficients of convection for the iButtons and the LEs. This would allow for the calculation of the radiosity for each LE. However, the model was not sufficiently robust to provide consistent convective coefficients, possibly due to unaccounted for local weather fluctuations, which varied the amount of convected heat reaching the iButtons. A simpler method was developed which involved using the iButton temperatures (TiB) alone, and which provided a relative, rather than absolute, measure of the radiosity from the LEs.

iButton Temperatures and Emissivity
The paint on the iButtons was applied manually to the surface of each iButton, resulting in a natural variation in quality and thickness. This variation would result in a slightly different emissivity for each iButton (εiB), ultimately affecting the TiB. The purpose of this section is to explain why the variation in emissivity was not significant.
Approximate values of the convective coefficients were used from the previous modelling work to determine how much the variation in the black iButton paint would affect TiB. It was found that changing TiB by up to ±1.0 °C changed εiB by a maximum of ±0.2. As the emissivities are unlikely to vary by this large amount, the effect of emissivity variation could be ignored, as it was negligible when compared with the uncertainty in the iButton temperature loggers themselves (given by the manufacturer as ±0.5 °C and ±1.0 °C for high and low temperature iButtons, respectively). The fact that TiB was so insensitive to the εiB was perhaps another reason why the initial energy balance modelling described above was unsuccessful.

iButton Temperatures and Convection
The resultant iButton temperatures TiB are likely to be more sensitive to the convective heat transfer from their surface than to the emissivity. The iButtons were used in the downward facing position, and because of this, the heat loss due to forced convection (wind speed), was a more significant factor than the heat loss due to natural convection, which is of the order of 1 W·m −2 ·K −2 under calm conditions [49]. However, as the LEs were all measured concurrently, each iButton would have responded similarly to the changes in wind speed (and in fact, to all weather conditions). The relative iButton temperatures, and hence the relative radiosity measured by the iButtons, would thus be comparable.

View Factors
The height of the iButton above the LE was chosen so that its shadow did not impact significantly on the LE, but such that the majority of the radiation it received would be from the LE rather than from the surrounding grass. The optimisation of height was discussed in terms of the view factor (VF) by Cueto et al. [33], who measured net radiation from various LEs using a net radiometer. The VF of the LE to the measurement instrument is defined as the fraction of radiation received by the instrument from the LE as a percentage of the total radiation received. In this case, the temperature of each iButton was affected by the field of view (FOV) of its painted face which was pointing directly down at the LE and estimated to be ~150°. The assumption was that the sides and back of the iButton were well insulated within the polystyrene block.
The initial experimental design using phase 1 LEs was developed for measuring surface temperature, where interference from nearby LEs was not considered to be significant as long as they were all low-set and located a reasonable distance apart on the ground. However, when using iButtons, this set up causes the iButton to receive radiosity not only from the LE plus a smaller fraction from the surrounding turf grass, but also from the PP bag, the plastic crate, and the nearby LEs. Hence, more complex analyses are required to determine the radiosity purely from the LE alone.
To calculate the approximate VF of each iButton/LE configuration, each LE was assessed to be either a circular disc, or a rectangle. The iButtons above circular LEs were treated as a circle of radius 8 mm; however, for ease of calculation, the iButtons above the rectangular and square LEs were considered to be square with side lengths 7 mm. An online calculator by Howell [50] was used to find the VFs for each iButton/LE pair plus their surrounding grass, and, where required, the edges of the PP bag. The rectangle-to-rectangle in a parallel plane function was used to determine the VFs for the rectangular LEs, and the disk-to-parallel coaxial disk of unequal radius function was used for the circular LEs. Anything other than LE, grass, and PP, was termed Other, and the VF was denoted as VFOther. Equation 1 relates the combination of iButton temperatures and VFs.
where TiB,LE+grass+PP+Other is the iButton temperature, TiB,LE is the iButton temperature due to the LE only, VFLE is the VF of the LE only, with similar nomenclature applying to grass, PP, and Other. Equation 1 was used to recalculate TiB,LE for each LE, by removing the fraction of energy from the visible turf grass (using the TiB,grass measured at the same time), and, for some LEs, by also removing the TiB,PP, which was assumed to have a similar TiB as white shade cloth. VFOther was the remaining fraction of the total FOV once the VFs of the LE, grass and PP were summed. As there were no data for TiB,Other, in these analyses it was assumed to be equal to TiB,grass, i.e., everything other than LE and PP was assumed to be turf grass. Considering this, Equation 1 became Equation 2:

Material Properties Affecting iButton Temperatures
The expected values of each TiB,LE can be estimated based on the material properties of thermal inertia, colour (albedo), roughness, and moisture content. As discussed in the introduction, darker but higher thermal inertia LEs (such as asphalt and black polished stones) have the capacity to absorb more incident energy during the day, hence their TiB,LE would be cooler than that of dark LEs with lower thermal inertia (such as black shade cloth, BPP and artificial grass). Convected heat from the LEs can be removed by wind currents [29], hence the convected heat received by the iButtons will vary with local wind conditions. In comparison, the surface of a white LE reflects most of the incident radiation onto the iButton and is less affected by wind conditions. Hence, the TiB,LE of dark LEs would tend to be lower, whilst that of white LEs would tend to be higher. LEs with rough surfaces have a higher specific surface area, i.e., more surface area available for cooling through forced convection, and would thus have a lower TiB,LE than smooth surfaces under windier conditions due to less radiant LWR. A higher specific surface area also has a lower albedo, as light has more chances to be internally reflected before escaping from the surface. LEs with a high moisture content are able to use that moisture for cooling through evaporation (or evapotranspiration in the case of plants), and would have a lower TiB,LE.

Relative Radiosity Ranking
Each LE was compared with the reference material of WPP in order to find its relative ranking with respect to the other LEs. The temperature of the iButton over the WPP (TiB,WPP,) is always going to be the hottest during the daytime as it will have the most reflected incident radiation (due to WPP having the highest albedo), and because the WPP is insulated, no heat will be conducted away into the ground. The other LEs absorb more radiation (than the WPP), which heats their surface. Some of this is re-emitted as LWR, some is convected away, and some is conducted into the body of the LE and then into the ground. Consequently, the black iButton above these LEs (TiB,LE) will never be as hot as that above the WPP. Thus, WPP was an ideal material to be used as an upper reference point and TiB,WPP was used to normalise each TiB,LE.
A direct ratio of TiB,LE over the black iButton temperature of the WPP (TiB,WPP) would not be consistent across the day, because the iButtons are also affected by the varying ambient temperature. Removing Tamb from each iButton temperature removes this diurnal variability from the data. In the evening, the WPP temperature drops below ambient temperature because it has no thermal inertia and is radiating heat to the cold night sky. This method was therefore not used after sunset. All TiB were first corrected for VFs. Equation 3 was used to determine the daytime ratio (DR) for ranking the radiosity of the LEs relative to the WPP.
The uncertainties in TiB,WPP were ±1.0 °C, whilst the uncertainties in the other iButtons were either ±1.0 or ±0.5 °C, depending on their type. The uncertainty in Tamb was estimated to be ±0.5 °C. The standard errors (SE) were calculated assuming random errors and a normal distribution [51]). DR was averaged (DRav) over the time periods where the SR was greater than 70% of that occurring at solar noon on that day, in line with the standard for measuring reflected radiation from objects in the field [38]. For autumn this was 09:30-15:00; for winter, 09:20-15:00; for spring, 09:00-15:00; and for both summer measurements, 09:00-15:40.
Whilst the iButton and Tamb errors are relatively small, the uncertainty in the mean of the data due to the methodology used to derive the day−time ranking was expected to be higher than that calculated from the instrumentation alone. This is because the factors affecting the measurement and calculation of DRav relied on physical properties of materials, convection coefficients, and on varying weather parameters. Hence, the DRav uncertainty was found by calculating the standard error (SE) of the actual DRav data for each LE during the daytime hours, where the SE is equal to the standard deviation divided by the square root of the number of sampling points.
Microsoft Excel's Data Analysis package (Microsoft Corporation, Makati, WA, USA) was used to perform multiple two-sample t-tests (assuming equal variances) on both the seasonal DR data from phase 1, and on the summer DR data from phase 2. The summer 2018 and summer 2019 DR data was also compared using this t-test. This analysis tested the null hypothesis that the means between each LE whose rankings placed them next to each other, were the same (at a significance level α = 0.05).

Delimitations
This study was performed under clear sky conditions with a high SVF and with minimal wind. The results may not be applicable to measurements made under cloudy skies or where wind gusts are strong enough to significantly cool LEs and iButtons. Ambient conditions would be best measured using a local weather station; however, the local weather station malfunctioned during some of the test work and data was used from the nearby meteorological station. It is recognized that in reality, some of the hardscape LEs are coupled to the ground using more complex substrates than the sand used in this study. Furthermore, mechanical compaction of the asphalt LE was not performed according to standard road surface preparation. The number and type of LEs tested in phase 2 was limited by the time available for both construction and measurement (LEs were on a school oval), as well as by the cost of these LEs. Hence, replicate LEs were not tested. The moisture content of the LEs and of the turf grass surrounding them was not measured as it was outside the scope of this study; however, moisture content contributes to surface temperature through evaporative cooling.

Weather Conditions
Apart from 21st January, 2019, which had slightly cooler ambient temperatures and a higher relative humidity than typical for this season, the weather conditions for each of the days measured were close to typical for that season, according to data from the Bureau of Meteorology [52]. The ranges and conditions for both studies are given in Table 3, and shown in Figure 7. Smoke haze from a bush fire was present from around 13:30 to the early evening on the spring day (5 th October, 2017) and this can be seen in Figure 7 as a reduction in the SR at these times. This may have led to lower than expected LE temperatures for this season. In summer 2018 at ~13:30, the ambient temperature dropped whilst the relative humidity and the wind speed increased suddenly, possibly due to a weather front moving in. The effects of this weather change are seen in results presented later.

View Factors
For phase 1 LEs, VFOther ranged from ~1% up to ~8% of the total FOV of the iButton depending on the LE. Hence the assumption made that TiB,Other was equal to TiB,grass, introduced an uncertainty in the calculation of TiB,LE, which was dependent on the difference between the radiosity of turf grass and the radiosity of Other. As the VFs of Other were relatively small, this uncertainty was also expected to be small. The VFs of the phase 1 LEs, where most LEs had an approximate surface area of ~0.13m 2 , were between 0.43 and 0.75 (Table 4). VFs varied slightly with season as the LEs were re-installed each season and, whilst the order of installation was kept the same, the distances between LEs varied. The large VFgrass for these LEs (Table 4) introduces further uncertainty due to the natural variability of the turf grass across the oval and hence around each LE. Whilst the variation in VFs between phase 1 LEs was not ideal, a more significant factor was the variation in the surrounding turf grass quality and the distance to the nearest LEs. To reduce these uncertainties, the height of the iButton over the smaller LEs would have to be reduced in order to increase the VFLE. However, this would introduce new uncertainties as the size of the shadow of the iButton would increase, thus reducing the radiosity from the LE.
For the larger phase 2 LEs, VFLE were all 0.92, as these LEs all had the same surface area of 1.44 m 2 . Calculations from [53] estimated VFOther to be ~1% of the FOV of the iButton. Hence, everything other than the LE was assumed to be turf grass, with a VF of 0.08.
In summary, the uncertainty in the TiB,LE from the phase 1 LEs is much greater than for the phase 2 LEs due to the assumptions made regarding the VFs and the variability in the turf grass. Table 4. Seasonal radiosity view factors to the iButton surfaces from the small phase 1 landscape elements (LEs), and from the surrounding polypropylene (PP) and turf grass.

Seasonal Comparison-Phase 1
Equation 2 was used to calculate TiB,LE for each LE (data given in Appendix A) based on the VFs given in Table 4. Seasonal data were calculated by normalising all of the TiB,LE with TiB,WPP as per Equation 3. Both seasonal data and the summer 2019 data were analysed using this method. Examples of the temporal daytime ratio data from summer 2018 are shown below, with data separated for clarity (Figures 8-10). The standard errors (SE) given in these figures were calculated assuming random errors and a normal distribution [51]. Temporal data from other seasons were similar but are not presented here. The tails of these figures are extremely unstable due to the nature of the ratio calculation. For example, as the SR decreases to zero at around 18:30, TiB,WPP approaches Tamb, so the denominator (TiB,WPP−Tamb) approaches zero, and then becomes negative as TiB,WPP begins to reflect the cold sky temperature. This results in a very large positive ratio just on sunset, followed by a very large negative ratio, as seen in the edges of the graph. Hence this method is only valid during the day. Figures 8-10 illustrate that the normalisation process is able to take account of the weather variability and varying solar radiation as seen by the reasonably consistent ratio during the daytime from 09:00 to 15:40 (standard errors ≤0.02). This is despite the larger uncertainties in this data due to the experimental design and the corresponding assumptions made around the VFs of these phase 1 LEs.
However, some generalisations can be made. The DRs of some of the low albedo, higher thermal inertia LEs (red and grey pavers, and the polished black stones) ( Figure 9a and Figure 10a, respectively) have a slight general trend upwards across the day. However, the high albedo, higher thermal inertia LEs (concrete, sandstone pavers, and limestone) (Figure 8b, Figure 9a and Figure 10a, respectively) have a fairly consistent DR through the day. It seems the daytime radiosity may be driven more by albedo than by LWR. This is consistent with results by [12] who found that increases in albedo significantly increased the upwards SWR whilst delivering only a minimal difference in the upwards LWR. In the later afternoon, the accumulated heat from the day is released at a slower rate in the higher thermal inertia LEs when compared with WPP, and hence the DR increases slowly for these LEs. In

Relative Average Radiosity
The DR data from phase 1 were averaged for each season over the times given in Table 3, to give DRav (data are given in Appendix B). DRav gave a relative ranking for each LE for each season, based on its radiosity ( Figure 11). The cross-seasonal average DRav was calculated and is shown in Figure  11 as the thick black line, with the data sorted from lowest to highest according to this average. The standard errors of the measured data for each season and LE were small (around 1%) and are shown as error bars in Figure 11. Hence, the differences seen between the seasonal DRav of each LE are likely to be significant. However, the uncertainties due to the experimental design are likely to be more significant, as discussed previously. On inspection of the DRav data, each LE was observed to lie within a band of around ±0.05 across all seasons (an uncertainty of ±5%). As 5% is not a large uncertainty for real world data, the seasonal differences would not be very significant for all practical purposes. Hence, the relative radiosity from the LEs could be considered to be similar across seasons. When calculating the cross-seasonal average for each LE, the autumn data points for limestone block and for asphalt were not included. Both of these LEs have the same DRav in summer, winter, and spring, but their autumn values are quite different and appear to be outliers. The possible reasons for this were that the limestone block was bought new for autumn but then stored outside between seasons. Aging of the new surface is suspected to have occurred, reducing the albedo for subsequent seasons and causing DRav to reduce. Conversely, asphalt had a lower DRav in autumn, as it was newer in this season and was therefore a darker black, becoming lighter in colour after aging outside between seasons. The winter DRav of red pavers and white shade cloth are outside the 5% band, but these data points were retained for the cross-seasonal average calculation as they were not significantly outside the band and their DRav from other seasons were also not consistent. The variance in the red pavers data may have been due to different arrangements of the six unevenly speckled pavers between seasons. Within the 5% band, the dry and moist soils are also slightly variable across seasons due to their moisture content, which is assumed to have been different, although this was not measured.
As discussed in Section 2.3.5, the effects of the material properties of the LEs can be seen in Figure 11, where white and lighter coloured LEs cause the highest DRav (white shade cloth, white crushed rock, white sand, sandstone pavers, and cream shade cloth). The lowest DRav were from moist soil, pine bark mulch, polished black stones, asphalt, and laid artificial grass, all of which could be classified as having higher surface roughness and/or moisture, and/or having a darker colour.
In general, white and lighter LEs had a higher radiosity than the darker LEs, indicating that during the daytime, SR is the dominant driver of the radiosity from the LEs. Based on the results shown in Figure 11, LEs can be measured at any time of the year using WPP as a reference, and the DRav will provide the relative radiosity of that LE to within ±5%. Table 5 shows the results from the two-sample t-tests on the LEs ranked next to each other in Figure 11. All LEs, apart from six pairs (p-values given in bold), show a statistically significant difference in their means. The negative value of t-Stat indicates that, moving down the table, each mean is larger than the next. Where the p-value >0.05, further t-tests were performed, comparing the LE with the two rows beneath it in the table, e.g., asphalt was compared with soil (dry), and soil (dry) was compared with turf grass. These further t-tests all had p-values <0.01, indicating that these means were all significantly different, and that the position of each LE in the ranking table would only vary by one place at most.
Whilst the focus of the phase 1 test work was a cross-seasonal comparison, the absolute ranking position of each LE was potentially affected by the different ways these LEs were placed in contact with the ground (Section 2.1.1), and the consequential variation in their effective thermal inertia. Although all LEs were buffered from the underlying soil by the turf grass, it is assumed that the air volume underneath each LE would be dependent on its mass, e.g., heavier LEs such as concrete and pavers would have had the least air volume, whilst lighter LEs such as artificial grass and the shade cloths would have had a greater air volume. Similarly, LEs mounted on the PP alone would also compress the turf grass by an amount dependent on their mass, but with the added thin layer of PP. LEs on PP and then on plastic crates had the greatest air volume (apart from decking). The effect of this changing air volume beneath the LEs was not able to be quantified.
As each LE was in contact with the ground in the same manner across each season, this varying air volume does not affect the cross-seasonal result that the relative radiosity of a particular LE will be the same within ±5% regardless of season. However, the ranking order in Figure 11 may be slightly different than if all of the LEs were coupled in the same manner. To remove the variability associated with ground coupling and to provide a more accurate ranking, the phase 2 investigation was undertaken.  Figure 12 shows the temporal daytime ratio data for the larger, ground-coupled LEs measured in summer 2019. Similarly to the higher thermal inertia LEs for the phase 1 LEs, the DR of the pavers, asphalt, and concrete showed a slight upwards trend throughout the time period 09:00 to 15:40 ( Figure 12a). The slope of this trend appeared to be lower for the phase 2 LEs; however, the increased ground coupling of these LEs meant they were warmer than the smaller LEs at the start of the day. This was confirmed from surface temperature measurements, which showed that these phase 2 LEs were around 4 °C warmer than the phase 1 equivalent LEs just before sunrise (relative to ambient temperature). In contrast, the softscape LEs ( Figure 12b) have a fairly consistent DR across the day, similar to the equivalent phase 1 LEs. Figure 13 shows the summer 2019 DRav from phase 2 LEs ordered from lowest to highest. Although the LEs from phase 1 are not all identical to those used in phase 2, they are shown in this figure to give some idea of how they compared. The uncertainties shown are the standard errors for each LE, based on the measured DRav data. However, the uncertainties in the phase 1 LEs are likely to be higher, as discussed previously in the View Factors section.  Although the data are difficult to compare given the variability in product types, the phase 1 LEs generally seem to have a lower DRav than the corresponding similar phase 2 LEs. The p-values for this data are given in Table 6. Only the means of grey and red pavers were the same (p-values > 0.05). The positive values of t-Stat for moist soil, sandstone pavers, and white rocks indicated that for these LEs, the summer 2018 means were higher than the summer 2019 means. One reason why the 2018 DRav data for most of the LEs are generally lower than 2019 data is that the smaller LEs will have less thermal inertia than the ground-coupled LEs. This means that their response to changes in ambient conditions will be faster and they will follow the WPP changes more closely than the larger LEs will. Hence, the DRav for the smaller LEs would be lower. Another reason is that the LEs are warmer in their central area, but become cooler towards their edges [28]. The surface area of the cooler edges makes up a larger proportion of the total surface area of the small LEs than it does of the larger LEs. The iButtons used in this study measured the whole LE, including the edges. Thus, the thermal radiation received by the iButtons from the smaller LEs will be less than that received by the larger LEs, resulting in a lower DRav of smaller LEs. Note that for all LEs: Number of observations = 41, t Critical two-tail = 1.99.

Summer Comparison of LE Size and Ground coupling-Phase 2
The extremes of the phase 2 LEs presented in Figure 13 show that crushed white rock had the highest DRav, followed by sandstone pavers and decking, whilst the moist soil and the plants had the lowest. Plants use the incident radiation to evaporate water, which reduces the amount of radiation available for heating, hence the plants and similarly the moist soil are ranked lowest. Apart from moisture content, the DRav seems to be largely dominated by the albedo of the LEs rather than other material properties.
It is believed that the radiosity in the evenings, where no SR is present, would be dominated by the thermal inertia. Evening measurements were taken for all of the LEs across all seasons for this study; however, the analysis was more complex, and initial results indicate that data should be collected over at least three or more evenings in order to account for the faster changing evening ambient temperature, sky temperature, and wind speed. Hence, this data has not been presented here.
In summary, the larger, ground-coupled LEs provided a larger VF of the LE for the iButtons than did the small LEs, and thus assumptions about the background material had an insignificant effect on the accuracy of phase 2 DRav. For this reason, the larger LEs are recommended for future radiosity tests, given an iButton height of 200 mm.

Ranking LEs
The summer 2019 values of DRav are given in Table 7 and indicate the relative radiosity ranking of the LEs. The results from the two-sample t-tests on the LEs ranked next to each other in Figure 13 are also shown. Six pairs of LEs do not have statistically different means (p-values given in bold). However, similar to the phase 1 seasonal tests, the negative value of t Stat indicates that, moving down the table, each mean is larger than the next. Further t-tests were performed on the six LEs, comparing each of these with the LE two rows beneath it in the table. These t-tests all had p-values <0.01, apart from dry soil and red pavers. These LEs were further tested against LEs lower in the table. The mean of the dry soil was not significantly different to red pavers or BPP; however, it was significantly different to concrete (x40mm), the next LE in the table (p-value = 0.03). Similarly, the result for red pavers was not significantly different to BPP or concrete (x40mm), but was significantly different to decking (p-value <0.01). These results indicate that there is some overlap in the ranking order of the LEs in Table 7, particularly in the mid-range of relative radiosities. The ranking order of the LEs in Table 7 is found to be the same as similar LEs tested by other authors. Cueto et al. [33] found day-time upwelling radiation from six surfaces were ordered from least to most as: grass, asphalt, concrete, clay, and a white painted polystyrene. Similarly, Montague and Kjelgren [15] found their six surfaces to be ordered from least to most as: turf, lava rock mulch, asphalt, pine bark mulch, gravel, and concrete. Table 7 does not give absolute values for radiosity, hence it is unable to be used for calculating human thermal comfort (HTC) levels or indexes such as the physiological equivalent temperature (PET), which are also affected by other parameters such as wind speed, air temperature, age, and clothing type [54]. However, DRav does provide guidance for choosing a more comfortable LE for a particular location or situation. For example, in areas of high daytime pedestrian use, high radiosity LEs would decrease HTC levels [10,12], and LEs such as plants, mulch or asphalt would be better (depending on the type of usage the surface is required to withstand). Similarly, LEs with high radiosities would not be suitable for use close to buildings, where the higher radiant energy would cause increased building envelope temperatures during the time of occupation [40].
Decking seems to have an unusually high DRav, perhaps because the surface of the decking, although unfinished, appears to be quite reflective. As the surface ages, or with application of a surface finish, the reflectivity of decking would change, indicating the importance of being able to easily measure DRav at any time. Although moist soil would not be recommended for use as a LE as such, it does provide a good contrast with the dry soil, indicating the important role of evaporative cooling in reducing radiosity. This supports findings by Lindberg et al. [55] that mean radiant temperature in full sun was less over moist grass than over asphalt, likely caused by the presence of moisture, as the albedos were similar. Table 7 shows that plants are the most suitable LEs to use if radiosity is to be minimised in a landscape, with the assumption that they are watered as required.

Limitations
The repeatability of the results would be improved by measurements over longer time periods during each season. Similarly, analysis of evening radiosity measurements indicated that more data should be collected over a number of evenings to improve accuracy. The lower VFs of the phase 1 LEs compared with their VF of the surrounding grass, PP, plastic crate, and nearby LEs, led to greater uncertainties in their DRav, uncertainties which were decreased when using the larger LEs.

Conclusions
Black painted iButton temperature sensors, along with the normalised methodology, performed well as relative low-cost radiometers for measuring multiple LEs concurrently. Between seasons, the range in values of the normalised average radiosity varied by up to 10% absolutely. Hence, LEs can be measured at any time of the year using WPP as a reference, and the DRav will provide the relative radiosity of that LE to within ±5%. Future work could involve exploring sources of the slight variations between seasons to determine if they are caused by the method (e.g., the size of the LEs), the weather conditions, or by actual physical differences in the LEs between seasons.
When comparing similar material type phase 1 (small) with phase 2 (larger, ground-coupled) LEs, most of the small LEs had lower relative radiosities than the large LEs.
Accuracy of the measured radiosity was improved when using larger, ground-coupled LEs due to the increased VFLE (0.92). It is therefore recommended that larger, ground-coupled LEs be used in future measurements.
Daytime radiosity appears to be driven more by albedo than by LWR, consistent with results by [12]. During the daytime, white and lighter coloured LEs caused the highest iButton temperatures and would therefore have the most detrimental effect on nearby objects such as buildings and/or people in hotter periods, potentially increasing energy costs [56,57], and/or heat stress issues, respectively [10,12,55]. Conversely, white and lighter coloured LEs would reduce urban heat if they were positioned in open spaces or on rooftops, with a high SVF [58,59]. Plants and moist LEs caused the smallest rise in daytime iButton temperature and would be most beneficial for managing daytime urban heat [59].
The use of low-cost iButtons to measure relative radiosity from numerous LEs concurrently presents a new way to examine the relative effect of different LEs on the urban environment. They may be a useful tool when using landscape design for urban heat mitigation. Acknowledgments: Appreciation to Bruce Ivers for his assistance in purchasing some LEs and for the use of his vehicle and trailer. Thanks to Grant Bayne for his aerial photography, and to the Booragoon Primary school for the use of their oval.

Conflicts of Interest:
The authors declare no conflict of interest.