<figure><img src="image_1.png" /></figure>

HAL
open science

# Modeling Mean Radiant Temperature Distribution in Urban Landscapes Using DART

Maria Angela Dissegna, Tiangang Yin, Hao Wu, Nicolas Lauret, Shanshan Wei, Jean-Philippe Gastellu-Etchegorry, Adrienne Grêt-Regamey

► To cite this version:

Maria Angela Dissegna, Tiangang Yin, Hao Wu, Nicolas Lauret, Shanshan Wei, et al. Modeling Mean Radiant Temperature Distribution in Urban Landscapes Using DART. Remote Sensing, 2021, 13 (8), pp.1443. 10.3390/rs13081443. hal-04644120

HAL Id: hal-04644120

https://hal.science/hal-04644120v1

Submitted on 12 Jul 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

# Modeling Mean Radiant Temperature Distribution in Urban Landscapes Using DART

Maria Angela Dissegna<sup>1,2,*</sup>,<sup>1</sup> Tiangang Yin<sup>3</sup>,<sup>1</sup> Hao Wu<sup>1</sup>, Nicolas Lauret<sup>4</sup>, Shanshan Wei<sup>5</sup>,
Jean-Philippe Gastellu-Etchegorry<sup>4</sup> and Adrienne Grêt-Regamey<sup>1,2</sup>

<figure><img src="image_4.png" /></figure>

check for
updates

**Citation:** Dissegna, M.A.; Yin, T.; Wu, H.; Lauret, N.; Wei, S.; Gastellu-Etchegorry, J.-P.; Grêt-Regamey, A. Modeling Mean Radiant Temperature Distribution in Urban Landscapes Using DART.

*Remote Sens.* **2021**, *13*, 1443. https://doi.org/10.3390/rs13081443

**Academic Editors:**
Panagiotis Bismanidis and
Zina Mitraka

Received: 31 January 2021

Accepted: 29 March 2021

Published: 8 April 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** The microclimatic conditions of the urban environment influence significantly the thermal comfort of human beings. One of the main human biometeorology parameters of thermal comfort is the Mean Radiant Temperature (Tmrt), which quantifies effective radiative flux reaching a human body. Simulation tools have proven useful to analyze the radiative behavior of an urban space and its impact on the inhabitants. We present a new method to produce detailed modeling of Tmrt spatial distribution using the 3-D Discrete Anisotropic Radiation Transfer model (DART). Our approach is capable to simulate Tmrt at different scales and under a range of parameters including the urban pattern, surface material of ground, walls, roofs, and properties of the vegetation (coverage, shape, spectral signature, Leaf Area Index and Leaf Area Density). The main advantages of our method are found in (1) the fine treatment of radiation in both short-wave and long-wave domains, (2) detailed specification of optical properties of urban surface materials and of vegetation, (3) precise representation of the vegetation component, and (4) capability to assimilate 3-D inputs derived from multisource remote sensing data. We illustrate and provide a first evaluation of the method in Singapore, a tropical city experiencing strong Urban Heat Island effect (UHI) and seeking to enhance the outdoor thermal comfort. The comparison between DART modelled and field estimated Tmrt shows good agreement in our study site under clear-sky condition over a time period from 10:00 to 19:00 ($R^2$ = 0.9697, RMSE = 3.3249). The use of a 3-D radiative transfer model shows promising capability to study urban microclimate and outdoor thermal comfort with increasing landscape details, and to build linkage to remote sensing data. Our methodology has the potential to contribute towards optimizing climate-sensitive urban design when combined with the appropriate tools.

**Keywords:** mean radiant temperature; DART; 3-D urban landscapes; urban vegetation; urban microclimate; outdoor thermal comfort

## 1. Introduction

The increased number of heat waves due to global climate change experienced in growing urban areas across the world directly and indirectly affects thermal comfort and health of urban residents, particularly in the tropics. The environmental factors controlling outdoor thermal comfort are the Mean Radiant Temperature (Tmrt), wind speed, air temperature, and humidity [1]. Outdoor thermal comfort and heat related mortality are more influenced by the Tmrt and less by the other three factors [2]. Tmrt is considered an essential bio-meteorological variable having a strong influence on thermal comfort indices

<footnote n="1">¹ Future Cities Laboratory, Singapore ETH Center, 1 CREATE Way, Singapore 138602, Singapore; hao.wu@sec.ethz.ch (H.W.); gret@ethz.ch (A.G.-R.)</footnote>

<footnote n="2">² Chair of Planning of Landscape and Urban Systems, ETH Zurich, 8093 Zurich, Switzerland</footnote>

<footnote n="3">³ Earth System Science Interdisciplinary Center, University of Maryland, College Park, MD 20740, USA; tiangang.yin@nasa.gov</footnote>

<footnote n="4">⁴ Centre d'Etudes Spatiales de la BIOsphère (CESBIO)—UPS, CNES, CNRS, IRD, Université de Toulouse, CEDEX 9, 31401 Toulouse, France; nicolas.lauret@cesbio.cnes.fr (N.L.); jean-philippe.gastellu@iut-tlse3.fr (J.-P.G.-E.)</footnote>

<footnote n="5">⁵ Singapore–MIT Alliance for Research and Technology, Singapore 138602, Singapore; shanshan@smart.mit.edu</footnote>

<footnote>* Correspondence: dissegna@arch.ethz.ch</footnote>
such as the Physiological Equivalent Temperature (PET) [3], Predicted Mean Vote (PMV) [1], Universal Thermal Climate Index (UTCI) [4], the Perceived Temperature (PT) [5], and Global Outdoor Comfort Index (GOCI) [6]. A comprehensive review of thermal comfort studies in urban spaces can be found in Lai et al. 2020 [7].

In solar exposure conditions, the energy loss or gain by radiation is the most important
heat flux in the energy balance; followed by convective fluxes of sensible and latent heat
and at last, the conductive heat flux [8]. The concept of Tmrt is based on the fact that
the net exchange of radiant energy between two objects is approximately proportional to
their temperature difference multiplied by their ability to emit and absorb heat [1]. This is
valid as long as the absolute temperatures of objects in question are large compared to the
temperature differences, allowing linearization of the Stefan–Boltzmann Law in the relevant
temperature range [9,10]. Different city structures in the same thermal region differ mainly
in their radiation temperatures [11]. This is due to the geometrical complexity of buildings,
their surface radiative properties and their wavelength dependency. These properties play
an important role in the radiation exchange and therefore, on its modelling. Numerical
tools such as Rayman [9], ENVI-met [12], SOLWEIG [13], TUF-3D [14], CityComfort+ [15],
and VTUF 3D [16] are increasingly being used by planners and architects to assess the bio-
meteorological performance of urban scenarios. Such tools vary in terms of dimensionality,
i.e., 1-D/3-D, accuracy in predicting the radiation fluxes, consideration of surface material
properties, distribution of surface temperatures, physical representation of the vegetation,
scale of analysis, speed of computation, and capability to assimilate remote sensing data as
input for simulations.

Remote sensing technologies have great potential to study the thermal behavior
of cities in space and time [17], particularly when combining information of urban sur-
faces and their spectral characteristics using physically-based radiative transfer models
(RTMs). Recent studies have successfully combined a variety of remotely acquired data
with RTMs to derive urban surface properties (albedo and temperature), energy fluxes, and
Tmrt [18–23]. These studies highlight an increasing demand for advanced tools capable
of 3-D analysis of detailed urban data including the vegetation. RTMs used in the field of
remote sensing were developed to compute accurate canopy spectral radiance, to invert
and evaluate existing data, and for the configuration of future satellite missions [24]. Only
a few RTMs simultaneously simulate the canopy spectral radiance and the canopy 3-D
radiative budget (RB). To our knowledge, advanced physically-based RTM models have
not been fully employed for modelling Tmrt. Therefore, this study seeks to (1) exploit the
advanced capabilities of a physically-based 3-D RTM for the estimation of Tmrt of complex
urban scenes, and (2) introduce its use for a broader application by planners and urban
climate scientists.

In this study we propose a new method for detailed 3-D modelling of Tmrt in com-
plex urban scenes by improving the radiative budget module of the Discrete Anisotropic
Radiation Transfer model DART [25]. The characteristics that make DART a powerful tool
to model radiation are summarized hereafter. (1) DART considers the full spectral domain
from ultraviolet to thermal infrared in solar/atmosphere radiation and optical properties.
Therefore, it provides more accurate results than models that consider that radiation only
belongs to two spectral domains (i.e., short and long wave domains) and/or that the canopy
optical property is only defined by its albedo and emissivity. (2) DART considers all types
of surface optical properties, including spectral reflectance, transmittance, and absorbance;
spectral emissivity of opaque material is equal to one minus spectral reflectance of this
material. Surface reflectance can be isotropic (Lambertian) or anisotropic using Fresnel
equations or predefined bidirectional reflectance distribution function, which allows DART
to consider any material. (3) DART can simulate vegetation either as a set of facet-like
leaves with specific locations, orientations, and optical properties, or as a 3-D distribu-
tion of turbid medium (i.e., infinite number of infinitely small plane elements) defined
by specific clumping, leaf angle distribution, Leaf Area Index (LAI)/Leaf Area Density
(LAD), and optical properties. Therefore, DART simulates light propagation below trees
conversely to models that represent trees as opaque volumes. The turbid representation of vegetation is adopted in this work. (4) DART simulates 3-D multiple scattering within vegetation and between vegetation and other materials (e.g., walls). This is very important in presence of high reflectance scene elements (e.g., green vegetation in the near infrared domain), and also to simulate the grey body (emissivity less than one) in long wave domain. (5) In addition to the radiative budget of canopies, DART simulates their remote sensing signal [21,24,26]. DART is considered one of the most advanced 3-D RTMS and that it has been positively evaluated in the framework of the Radiative transfer inter-comparison project (RAMI) [27,28].

Our methodology for modelling Tmrt can be applied at different scales and under a range of parameters including the urban pattern, surface material of ground, walls, roofs, and vegetation properties (coverage, shape, spectral properties, Leaf Area Index, and Leaf Area Density). Three dimensional scenes can be generated from multi-source remote sensing data as shown on our previous study on 3-D Reconstruction of urban landscapes from satellite data, where a detailed workflow for the derivation of 3-D urban scenes was presented, together with an evaluation of 3-D radiative budget over different urban typologies found in Singapore [21]; additionally, scenes can be generated from Aerial LiDAR Scanning (ALS) point cloud data [29] or using 3-D modelling software. Depending on the desired scale of analysis, vegetation turbid plots (i.e., volume of turbid material) can be generated from satellite data [21], from ALS point clouds [29] or using realistic 3-D models. Moreover, vegetation properties (e.g., LAI/LAD) can be retrieved from terrestrial laser scanning data (TLS) at plant scale [30,31], from ALS [32,33] at neighborhood scale or from other remote sensing systems. Finally, the spatial distribution of vegetation can be modelled from hybrid point cloud data [29]. Information on suggested data sources for scene creation can be found in Tables A1 and A3.

We illustrate our approach in Singapore, a city experiencing strong UHI and seeking to enhance outdoor thermal comfort. We present a field comparison and perform a sensitivity analysis to examine how building materials and vegetation properties influence Tmrt. Then, we discuss limitations of our approach, possible areas of improvement and future research directions.

## 2. Materials and Methods

### 2.1. *Scene Generation and Parameter Setting*

A 3-D model of the Tanjong Pagar District, Singapore, was reconstructed. The building height, vegetation shape, and distribution was obtained from the Aerial LiDAR Scanning dataset of Singapore [34]. This dataset was acquired by the Optech Pegasus HA500 sensor with a planar density of about 30 points/m² which proved to be well adapted to map vegetation. The buildings were reconstructed using 3DsMax software, obtaining a Level of Detail (LOD) of 1.2. The north part of the scene corresponds to the existing urban setting, whereas its southern part (former Tanjong Pagar Port) corresponds to a design proposal with a high-density urban form [29]. The scene elements were grouped by surface material, allowing them to be linked to their respective surface temperatures and optical spectral properties in the DART database. All trees were assumed to have the same LAI value of 1.5 (single-side leaf area divided by the tree projection area). This value falls within the range of measured LAI of urban trees found in Singapore [35]. Surface temperature is an essential term for simulations in the longwave domain. It can be derived from field observations of a calibrated thermal camera [36] or from an energy balance model, i.e., DART EB. It is provided as an input to DART either as thermal functions (i.e., mean and standard deviation of hourly modelled temperatures) or as a 3-D temperature distribution. DART uses these thermal functions to distribute the temperature values over the scene using a virtual sun illumination of the scene: the most irradiated surfaces are given the highest temperature values. In our study, the mean and standard deviation of hourly surface temperatures were estimated using the Urban Eco-hydrological model UT&C [37]. UT&C is a fully coupled energy and water balance model which accounts for the biophysical and eco-physiological characteristics of urban trees. UT&C
is a 1D model, providing as output the mean and standard deviations of surface temperatures of impervious roof, impervious ground, ground covered with grass, tree, sunlit wall, and shaded wall for each of the parametrized areas of interest. Our study area was classified in urban "typologies" to extract the neighborhood parameters such as building height, height to width ratio and fraction of vegetation. The model requires meteorological data of incoming shortwave and longwave radiation, rainfall, pressure, wind speed, air temperature, and humidity at forcing height above the urban canyon. The model was run over the time period from 1 May 2013 to 30 April 2014. The mean and standard deviation of the temperatures were obtained considering all the hours with a cloudiness of less than 50% with the aim to approximate the average conditions on rather sunny days in Singapore. Information on the values such as volumetric heat capacity and thermal conductivity are reported in the supplementary information on the UT&C model development document [37].

## 2.2. Mean Radiant Temperature Computation

To compute the Tmrt of Tanjong Pagar District, we carried out two sets of DART radiative budget simulations, one in the shortwave domain (0.3–2.5 µm) and one in the longwave domain (3–50 µm). Simulations were run at a resolution of 2 m horizontally and 1 m vertically. The use of "repetitive scene mode" in the longwave domain reduces the computation time without producing substantial differences on the resulting radiative budget. However, in the shortwave domain, the use of repetitive scene mode would cause the emergence of shadows which do not correspond to the reality. Therefore, we used "isolated scene mode" when running simulations in the shortwave domain. During the simulation process, DART converts the scene into voxels at a defined dimension. These voxels contain information on the material reflectance, surface temperature, and additional properties such as LAI/LAD for the case of turbid vegetation plots. DART stores the resulting radiation that has entered each voxel. The previous versions of DART did not compute the radiative budget on empty cells; therefore, it was not suitable for estimating and mapping Tmrt. The recently improved releases of DART's radiative budget 3-D allow one to store separately the direct, diffuse, and emitted irradiance reaching each and every voxel of the scene from six directions (N-S, E-W, and Upwelling-Downwelling). This allows to virtually transform the voxels into equivalent cylinders [38] on which the absorption coefficient of the human surface to solar radiation (1 − α) and emissivity of clothed body (ε) can be considered and be varied at the Tmrt computation stage. In the short waves (i.e., solar radiation), the absorption coefficient of the skin (1 − α) ranges between 0.55 and 0.85 depending on the skin color [39]. In the long waves (i.e., terrestrial radiation), the emissivity of the human skin is assumed to be 0.99 and 0.95 for normal clothing surfaces [40]. The complete list of the simulation parameters can be found in Tables A4 and A5. The surface temperatures used for each time step are reported in Tables A6 and A7.

The computation of Tmrt is done according to the method proposed by Peter Höppe [41], with further considerations [38]. These include (1) the representation of the human body as a cylinder as opposed to a box, (2) separated treatment of direct and diffuse short-wave radiation as well as for long-wave radiation, (3) consideration of solar elevation and azimuth for every time step. We used a cylinder with the measurements reported in Table 1 to derive the surface fractions $w_i$ of the human body. The cylindrical representation of the human body is centered at a height of 1.5 m from the ground.

<table><caption>**Table 1.** Characteristics of the cylindrical representation of the human body.</caption><thead><tr><th colspan="2">Cylinder Dimensions</th></tr></thead><tbody><tr><td>Diameter:</td><td>0.28 m</td></tr><tr><td>Height:</td><td>1 m</td></tr><tr><td>Total area:</td><td>1 m<sup>2</sup></td></tr><tr><td>roof area:</td><td>0.06 m<sup>2</sup></td></tr><tr><td>wall area:</td><td>0.88 m<sup>2</sup></td></tr><tr><td>Area of frontal view:</td><td>0.28 m × 1 m = 0.28 m<sup>2</sup></td></tr><tr><td>Surface fraction w<sub>i</sub> of Cylroof:</td><td>0.0616 m<sup>2</sup>/1 m<sup>2</sup> = 0.0616</td></tr><tr><td>Surface fraction w<sub>i</sub> of Cylwall:</td><td>0.28 m<sup>2</sup>/1 m<sup>2</sup> = 0.28</td></tr></tbody></table>
The mean radiant flux (Sstr) reaching a rotationally symmetrical (cylindrical) representation of the human body is computed at a specified layer height using Equation (1).

$$S_{\text{str, cyl}} = (1 - \alpha) \cdot \left[ \mathbf{w}_{\text{Cylwall}} \cdot \sum_{1-4} K_{\text{dir,hor}} + \mathbf{w}_{\text{Cylroof}} \cdot (K_{\uparrow\text{dir}} + K_{\downarrow\text{dir}} + K_{\uparrow\text{diff}} + K_{\downarrow\text{diff}}) + \sum_{1-4} K_{\text{diff,hor}} \cdot \mathbf{w}_{\text{Cylwall}} \right] + \varepsilon \cdot \sum_{1-6} \mathbf{w}_i \cdot \mathbf{L}_i \quad (1)$$

where

* α is Albedo of clothed human body (0.37)

* K<sub>dir, hor</sub> is Horizontal direct shortwave radiation (KE, KS, KW, KN)

* K<sub>↑dir</sub>, K<sub>↓dir</sub> is Vertical direct shortwave radiation

* K<sub>↑diff</sub>, K<sub>↓diff</sub> is Vertical diffuse shortwave radiation

* K<sub>diff,hor</sub> is Horizontal diffuse shortwave radiation (KE, KS, KW, KN)

* ε is Emissivity of clothed human body (0.97)

* w<sub>i</sub> is Surface fractions of the standing "cylinder man" (w<sub>Cylwall</sub> = 0.28 and w<sub>Cylroof</sub> = 0.0616)

* L<sub>i</sub> is Longwave radiation

The mean radiation temperature Tmrt [°C] results from the Stefan-Boltzmann law using Equation (2).

$$T_{\text{mrt}} = \sqrt[4]{\frac{S_{\text{str,cyl}}}{\varepsilon \cdot \sigma}} - 273.15 \quad (2)$$

where

* S<sub>str, cyl</sub> is the mean radiant flux

* ε is Emissivity of clothed human body

* σ is Stefan-Boltzmann constant 5.670374419 × 10<sup>-8</sup> W·m<sup>-2</sup>·K<sup>-4</sup>

## 2.3. Study Area and Data Collection

Our study area is located in Tanjong Pagar district, Singapore (Figure 1a). Singapore has a tropical humid climate with no distinctive seasons. Near-surface air temperature usually ranges from 23 °C to 32 °C. The "Cantonment Towers" site, a modern social housing compound with a playground surrounded by high-density housing blocks was selected to carry out the field evaluation (Figure 1b). Three net radiometers Kipp and Zonen CNR4 [42] were mounted on a mobile platform at 1.5 m from the ground. The minimum height suggested by manufacturer to avoid interference from the mounting structure to the readings is 1.5 m. Additionally, a Vaisala WXT536 weather station [43] and a Campbell Scientific 152 mm black globe thermometer [44] were mounted on the platform to record wind speed, precipitation, air temperature, and globe temperature, respectively (Figure 1c). Special care was made to keep the tower leveled and oriented with the geodetic north [38]. Figure 1b shows the location of the measurement station in the study site.

A time series composite was generated by collecting observations of the maximum incoming radiation from 20 February 2020 to 2 March 2020 between 9:30 and 19:00 with a time step of one minute. This was done with the aim to approximate the maximum attainable Tmrt values for the site over a day with clear sky Tmrt was estimated using the integral radiation measurement technique [10,41]. Calculations of Sstr,cyl were based on angular factors for a rotationally symmetric standing person using Equation (3).

$$S_{\text{str, cyl}} = (1 - \alpha) \cdot \left[ \mathbf{w}_{\text{Cylwall}} \cdot \sum K_{\text{dir,tot}} + \mathbf{w}_{\text{Cylroof}} \cdot (K_{\uparrow} + K_{\downarrow}) + 0.88 \cdot K_{\text{diff}} \right] + \varepsilon \cdot \sum \mathbf{w}_i \cdot \mathbf{L}_i \quad (3)$$

where

* K<sub>dir, tot</sub> is Direct and reflected horizontal shortwave radiation (KE, KS, KW, KN)

* K<sub>↑</sub>, K<sub>↓</sub> is Vertical shortwave radiation

* K<sub>diff</sub> is Diffuse radiation = min (KE, KS, KW, KN)

* ε is Emissivity of clothed human body

* w<sub>i</sub> is Surface fractions of the standing cylinder man (w<sub>Cylwall</sub> = 0.28 and w<sub>Cylroof</sub> = 0.0616)

* L<sub>i</sub> is Longwave radiation

Then, the Tmrt is determined using Equation (2).
<figure><img src="image_6.png" /><figcaption>**Figure 1.** (a) Location of the study area. (b) Location of the measurement station within the study site marked in red. (c) Bio-climate measuring station for determining the Tmrt.</figcaption></figure>

# 3. Results

## 3.1. Modelled Tmrt over Tanjong Pagar, Singapore

The Figure 2 shows the results of Tmrt simulation over the Tanjong Pagar district, Singapore, at 15:00 on 29 February 2020. In sunlit areas, the Tmrt at pedestrian level reaches 68 °C. The highest Tmrt values are observed near the sunlit walls of buildings. Tmrt decreases as the distance to the sun-exposed surfaces increases. Under trees, Tmrt varies between 40 °C and 28 °C. The combined effect of building shade and tree shade results in Tmrt values that are below 28 °C. This suggests a slight underestimation of Tmrt in areas without direct solar radiation.

<figure><img src="image_7.png" /><figcaption>**Figure 2.** Visualization of pedestrian level Tmrt of the Tanjong Pagar South district, Singapore, at 15:00 on 29 February 2020. The results are overlaid into a Cloud Compare display of the original 3-D scene used for simulations (http://cloudcompare.org/ (accessed on 13 July 2020)). The north part of the scene corresponds to the existing urban setting while the south is a proposed urban design for the former Tanjong Pagar Port [29]. The area marked in red corresponds to the Cantonment Towers site used for field evaluation.</figcaption></figure>
## 3.2. Field Evaluation

A smaller subset corresponding to the Cantonment towers was generated for the purpose of field evaluation. Simulations were carried out assuming absence of clouds over the scene and using hourly values of Aerosol Optical Depth (AOD) reported by NASA-AERONET Singapore. They were run at 1m resolution. Then, the values of DART cells with the center 1.5 m above the ground were extracted to compute Tmrt maps at different times of the day (Figure 3). The highest Tmrt values are mostly observed close to the walls that receive direct solar radiation, prominently at 10:00 on the east facing walls and at 16:00 on the west facing walls. Generally, higher Tmrt occurs at 16:00.

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Selected Tmrt maps for Cantonment towers site at different simulated times of the day.</figcaption></figure>

The Figure 4 shows the comparison of modelled and field estimated Tmrt. The green line shows the field estimated Tmrt during the sunniest day of the study period, 29 February 2020. The downward peaks indicate reduction of Tmrt due to the passage of clouds over the site. The blue line presents the maximum Tmrt value of each time step obtained from the net radiometers during the entire study period (20 February 2020 to 2 March 2020). The black asterisks indicate the DART-modelled Tmrt for 29 February 2020 assuming clear sky conditions. The results show good agreement between 10:00 and 19:00 h ($R^2$ = 0.9697, RMSE = 3.3249) when comparing DART-modelled Tmrt with the maximum Tmrt values reported over the entire study period (blue line). In the morning at 9:00, the modelled Tmrt was overestimated in 19.54 °C. This explained by the fact that the DART simulation for Cantonment towers was run on a small subset of the scene, therefore neglecting the shadows that are casted by buildings outside the scene, particularly at low solar angles. Around 13:00, we observe a slight overestimation of 2.16 °C in the modelled Tmrt. From 17:15 onwards, our measuring site is affected by the shadow cast by buildings within the scene. Tmrt decreases considerably until the sunset just after 19:00. Then, Tmrt is only determined by the longwave radiation fluxes and remains relatively constant throughout the night. From 18:00 to 19:00, in absence of direct solar radiation, we observe an underestimation of 5.55 °C and 6.395 °C respectively on our modelled Tmrt.

## 3.3. Radiation Fluxes Recorded at the Study Site

The diurnal longwave and shortwave radiation fluxes received from each of the six directions (i.e., North ↔ South, East ↔ West, Up ↔ Down) on 29 February 2020 are illustrated in Figures 5 and 6. In terms of longwave radiation, we observe a steady behavior from the six directions throughout the day ranging between 450 and 600 W/m². The
maximum values are observed around 14:00. Before midday, the upwelling and West directions have the strongest signal. After midday, the downwelling direction becomes more prominent.

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Comparison between DART Tmrt and field-estimated Tmrt. The green line indicates records from 29 February 2020. The blue line corresponds to a composite of maximum values attained at each time step over the study period. The black asterisks correspond to each time step of DART-modelled Tmrt.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 5.** Diurnal longwave (L) radiation fluxes measured from six directions (9:00 to 19:00) on 29 February 2020.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 6.** Diurnal shortwave (K) radiation fluxes measured from six directions (9:00 to 19:00) on 29 February 2020.</figcaption></figure>

The downwelling shortwave radiation is the predominant flux with a maximum around 13:00 with 980W/m² and changes considerably during the day. This flux is primarily determined by the sun's position and sky view factor of the site. The upwelling shortwave is due to the reflection of solar radiation by the ground. It peaks around 13:00 with 120 W/m² with a slight variation during the day. The radiation from the South direction is greater than that from the North, hovering around 250 W/m² at 1:00 p.m. and varying gently throughout the day. This is because Singapore geographically located at 1° North. The downward peaks correspond to the effect of passing clouds over the study site.

## 3.4. Sensitivity Analysis

A sensitivity analysis was performed to demonstrate (1) how much the variation of a biophysical vegetation property (LAI) impacts on Tmrt, (2) how much does ground surface temperature impacts Tmrt, (3) how do spectral properties of wall material impacts spatially the Tmrt at pedestrian level.

### 3.4.1. Effect of Variying LAI on Tmrt under the Tree Canopy

Urban trees are an important component that contributes to the urban microclimate because of their potential to attenuate solar radiation, evapotranspiration, and control the wind speed. In tropical humid regions, the cooling effect by trees is caused mainly by the reduction of Tmrt due to shading. The benefits provided by trees are constrained by several factors such as spatial arrangement, type, age, height, phenology, crown shape, characteristics of trunk and twigs, leaf size, LAI, and leaf reflectance. Strategic placement and optimal selection of vegetation are essential to obtain the desired Tmrt attenuation. Biophysical properties of vegetation such as LAI/LAD are required to obtain realistic Tmrt simulations and to assess the differences between species and planting regimes. The effect of varying LAI on the Tmrt under the tree canopy of three common tree species of Singapore namely Albizia Saman, Khaya Senegalensis and Tabebuia Rosea was explored and reported in Appendix E. A scene was created using 3-D models of the tree species, assuming the ground material to be grass. Simulations were run at 13:00 and 16:00 using the respective surface temperatures reported in Table A6. The results show that for the three species, the impact on Tmrt due to increasing LAI is large if LAI is low. The impact
decreases when LAI increases, and becomes minimal when LAI is larger than 5.5. Khaya Senegalensis has the highest potential for reducing Tmrt. For instance, with small LAI values (LAI = 0.1), we observed a difference in Tmrt up to 16.19 °C when comparing Khaya Senegalensis against Tabebuia Rosea. This is primarily due to the amount of twigs and branches and their homogeneous distribution throughout the canopy; this implies a reduction of Tmrt even when LAI = 0.

### 3.4.2. Effect of Surface Temperature on Tmrt for Different Ground Materials

To provide an idea of the effect of surface temperature of different ground materials on Tmrt, we simulated a small scene where the surface temperatures of grass, concrete and three types of asphalt varied between 20 °C and 60 °C. In the case of grass, the range was assumed to be 20 °C to 40 °C This is because a surface temperature of grass higher than 40 °C is unrealistic. The overall response in terms of Tmrt of the studied materials is almost linear to the change in surface temperature. It was observed that light materials tend to give higher Tmrt, due to their higher reflectance. The results of the sensitivity analysis and additional details are reported in Appendix F.

### 3.4.3. Effect of Spectral Properties of Wall Material on Tmrt

The radiative response of the urban surface materials plays an important role in the radiative budget and thus, on Tmrt and in the overall urban microclimatic conditions. To assess the influence of different wall materials on Tmrt, we recreated a small scene with three buildings of dimensions (L:10 m, W:24 m, H:18 m). The first building has 100% white walls, the second has 66.6% white walls, and 33.3% glass and the third is 100% glass walls. The simulation was run at 16:00 with a cell size (x, y, z) of 1 m. The results indicate values of 79.39 °C, 68.59 °C, and 64.47 °C Tmrt in the first cell adjacent to the respective building. A detailed representation of results together with the surface temperatures assumed for this exercise is reported in Appendix G.

## 4. Discussion

This work introduced a new method, which uses for the first time the DART model to estimate and map Tmrt at different scales. Our motivation to use the DART model to estimate Tmrt stems from its capability to assimilate a large variety of 3-D data derived from state-of-the-art remote sensing techniques; providing the possibility to better analyze how Tmrt is influenced by vegetation and its biophysical properties, by the optical properties of surface materials, by surface temperature, and also by local atmospheric conditions. The field evaluation indicated good agreement between DART-simulated Tmrt and field estimated Tmrt at the Cantonment towers measurement site. However, a substantial underestimation is observed in areas with absence of solar radiation and in the late afternoon. This might be explained by a number of considerations in which our method is based. Firstly, we assumed clear sky conditions throughout the day. This is valid as long as one seeks to map temperature extremes. However, in a tropical city such as Singapore, the presence of clouds and episodes of rain have a major impact in the actual mean radiant temperature. To give an idea of this, we plotted the field observed Tmrt under different meteorological conditions over the study period in Figure 7. The inclusion of clouds on DART simulations has been recently explored [45]. However, this has not been implemented in our study. Atmospheric conditions such as water content and aerosols are highly variable over the course of a day. They influence the amount of radiation reaching the ground and consequently, Tmrt. However, in the case of Singapore there is not much variation in terms of water content. The annual average of relative humidity is 84.0%. Therefore, we only included real time (local) AOD observations obtained from the NASA AERONET website for each time step simulated.
<figure><img src="image_12.png" /><figcaption>**Figure 7.** Observed Tmrt under different meteorological conditions in Singapore over a two weeks period.</figcaption></figure>

Secondly, the designation of optical properties was merely based on site inspection and matching with the closest material existing within the DART database, this, however, is a source of uncertainty since the spectral library of DART might not well represent the actual characteristics of the surfaces of our study area. An overview of the spectral signatures of construction materials and vegetation used on our simulations can be found in Appendix D. The correct consideration of optical properties of surface materials play an important role in the radiative budget and thus, on Tmrt. To obtain a more accurate representation of optical properties, spectral signatures could be collected using a spectroradiometer and measuring the emissivity of the surfaces in the actual site. Unfortunately, we did not have access to such equipment. A detailed methodology on the derivation spectral signatures of urban materials through emittance and reflectance spectroscopy can be found in Kotthaus et al. 2014 [46].

Thirdly, the variation of surface temperatures impacts pedestrian Tmrt, particularly in sunny conditions. This is expected since in sun hours the intensity of Tmrt is primarily driven by shortwave radiation. The influence of longwave emission tends to decrease with the distance between the emitting surface and the absorbing human body. The sensitivity analysis to changing the surface material of ground showed that the ground material greatly influences Tmrt. For instance, the maximum difference of Tmrt observed between the five materials was 7.5 °C when surface temperature = 20 °C. This difference decreases to 5 °C if surface temperature equals 60 °C. A possible reason of underestimation in the areas with no direct shortwave radiation is that the surface temperatures obtained from UT&C energy balance model were obtained considering all the hours with a cloudiness of less than 50% over the time period from 1 May 2013 to 30 April 2014. This resulted in relatively lower surface temperatures compared to the actual surface temperatures occurring on a day with no cloud cover.

Fourthly, the accuracy of LAI retrievals from indirect methods relies on the integrity of the data used and on the retrieval technique. The retrieval of biophysical properties of vegetation has received particular attention in recent years, nevertheless, important gaps still exist for example in the estimation of LAI of heterogeneous tropical vegetation and for individual urban trees [35]. Further field investigations on the relationship between LAI and under canopy Tmrt would help to validate and improve the predictions of our proposed method. In our study over Tanjong Pagar, Singapore, a uniform LAI value was assumed for all the trees in the site. If local LAI/LAD data becomes available, for instance retrieved from ALS [32], this information could be assimilated. Additionally, more
evaluation sites are required to further assess the robustness of our method under different vegetation and urban settings, and in different climatic zones.

Finally, our method currently provides Tmrt at pedestrian level (1.5 m above the ground). Mapping 3-D Tmrt for instance on facades, balconies and elevated terraces could be possible since Tmrt is computed per vertical column of voxels of the scene. This would require to select the voxels in proximity to walls and roofs and to ignore the empty spaces where a person could not stand. Despite the advantages of using DART to perform analysis at different scales, a high level of detail over large scenes will result in long computation times. Therefore, it is crucial for the user to balance the trade-offs between scene size, level of detail and computation time. A table with examples of computation times is provided in Table A2. With the continuous increase in computing power, this technical limitation will continue to decrease in the coming years.

## 5. Conclusions

The increasing availability of data of urban areas opens up new possibilities for detailed analysis in terms of urban microclimate and thermal comfort. Affordable laser scanners are being developed, opening opportunities for wide range of research on the processing and data retrieval for urban 3-D mapping. This unprecedented increase of data and computation capacity highlights the need for more detailed simulation tools and methods for planners and designers to evaluate the performance of existing areas and to assess future designs under different climatic scenarios. This study explored the potential of using the 3-D physically based RTM DART to model Tmrt. We presented a new method for the estimation and mapping Tmrt at different scales with a detailed consideration of surface materials and vegetation. The field evaluation showed good agreement between modelled and field estimated Tmrt. The impact of LAI on pedestrian Tmrt was explored as part of the sensitivity analysis. The results indicate a reduction of up to 38.2 °C Tmrt under canopy when LAI = 5.5 compared to an "exposed" setting. This highlights the importance of using accurate vegetation properties for simulations. The sensitivity analysis showed consistency across different resolutions and changing simulation parameters. We demonstrated a potential first application of our proposed method. Our approach can be used to visualize locations in need of interventions, to help to optimize climate sensitive urban design, and to support urban microclimate and outdoor thermal comfort studies when combined with adequate simulation tools for wind analysis and energy balance models.

A correct designation of surface temperatures is crucial to avoid erroneous simulations and wrong conclusions of the study. In future work, we will explore coupling our method with 3-D energy balance model such as DART EB [47] for a more precise designation of surface temperatures and a detailed calibration of input parameters. Additional validation exercises in different urban settings and other climatic regions would contribute to improve the robustness and applicability of this newly proposed method.

Quantifying the influence of urban design, construction materials, and vegetation on Tmrt can be very helpful in order to evaluate urban planning scenarios to reduce heat stress in existing urban areas, to promote enhanced thermal comfort in future developments and to contribute in reducing the energy consumption for cooling systems. Future research could, for instance, explore the effect of building morphology and urban patterns; which species to plant and their optimal spatial arrangement such as in street-tree scenarios; as well as the effect of construction materials on Tmrt by using the existing DART spectral database or by adding the actual spectral information of the materials found on the site of interest.

**Author Contributions:** Conceptualization, M.A.D. and T.Y.; methodology, M.A.D., T.Y., J.-P.G.-E. and N.L.; software, H.W., J.-P.G.-E., N.L. and T.Y.; validation, M.A.D.; formal analysis, M.A.D.; investigation, M.A.D. and S.W.; resources, T.Y. and A.G.-R.; data curation, M.A.D. and H.W.; writing—original draft preparation, M.A.D.; writing—review and editing, T.Y., S.W., J.-P.G.-E. and A.G.-R.; visualization, M.A.D. and H.W.; supervision, T.Y. and A.G.-R.; project administration, A.G.-
R.; funding acquisition, A.G.-R. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was conducted at the Future Cities Laboratory, Singapore-ETH Centre, which was established collaboratively between ETH Zurich and Singapore's National Research Foundation (FI 370074016) under its Campus for Research Excellence and Technological Enterprise Programme. Shanshan Wei is funded by the Singapore National Parks Board (Nparks) project "Remote Sensing for Urban Tree Management: Species Classification and Health Monitoring".

**Institutional Review Board Statement:** Not applicable.

**Acknowledgments:** The authors would like to thank the DART team of CESBIO (CNES, CNRS, IRD, and University of Toulouse) for the consistent support with DART model. Special thanks to Professor Wong Nyuk Hien (Building Department, National University of Singapore) for providing the Net radiometer devices used for field evaluation. We acknowledge the Singapore Land Authority for providing the LiDAR data and Naika Meili for providing the surface temperatures derived from the UT&C model the authors appreciate the valuable remarks and suggestions from anonymous reviewers.

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

## Appendix A

<table><caption>**Table A1.** Remote Sensing Inputs, Level of Detail and Recommended Scene Settings for Simulation.</caption><thead><tr><th>Data Source</th><th>Scale of Study</th><th>LOD Buildings</th><th>LOD Vegetation</th><th>Scene Dimensions</th><th>Cell Size (x,y)</th><th>Cell Size (z)</th><th>Reference</th></tr></thead><tbody><tr><td>VHR satellite imagery (WorldView2)</td><td>Town</td><td>1.1-1.2</td><td>location, approx. height and volume (cube)</td><td>1000 m * 1000 m</td><td>2 m-5 m</td><td>1 m or 2 m</td><td>i.e. Dissegna et al. 2019</td></tr><tr><td>Aerial</td><td>Neighbourhood</td><td>1.1-1.2</td><td>location, accurate height and volume (rounded)</td><td>500 m * 500 m</td><td>1 m-2 m</td><td>1 m</td><td></td></tr><tr><td>LiDAR</td><td>Town</td><td>2</td><td>location, accurate height and volume (rounded)</td><td>1000 m * 1000 m</td><td>2 m-5 m</td><td>1 m or 2 m</td><td></td></tr><tr><td>Scanning</td><td>Neighbourhood</td><td>2.1</td><td>location, accurate height and volume (rounded)</td><td>500 m * 500 m</td><td>1 m-2 m</td><td>1 m</td><td>i.e. Urech et al. 2020</td></tr><tr><td>3-D Objects of architectural designs</td><td>Neighbourhood</td><td>2</td><td>location, accurate height and volume (rounded)</td><td>500 m * 500 m</td><td>1 m-2 m</td><td>1 m</td><td>Created with any 3-D modelling software</td></tr><tr><td>Hybrid</td><td>Building</td><td>3.1</td><td>n/a</td><td>&lt;100 m * 100 m</td><td>0.5 m-1 m</td><td>0.5 m-1 m</td><td>Laubwerk</td></tr><tr><td>ALS/3-D designs</td><td>Tree</td><td>n/a</td><td>Explicit shape</td><td>50 m * 50 m</td><td>0.25 m-0.5 m</td><td>0.25 m-0.5 m</td><td></td></tr><tr><td>Neighbourhood</td><td>Town</td><td>1.1-1.2</td><td>location, accurate height and volume (rounded)</td><td>1000 m * 1000 m</td><td>2 m-5 m</td><td>1 m or 2 m</td><td></td></tr><tr><td>Building</td><td>Neighbourhood</td><td>up to 1.3</td><td>location, accurate height and volume (rounded)</td><td>500 m * 500 m</td><td>1 m-2 m</td><td>1 m</td><td>i.e. Urech et al. 2020</td></tr><tr><td>Building</td><td>Building</td><td>up to 3.2</td><td>location, accurate height and volume (rounded)</td><td>&lt;100 m * 100 m</td><td>0.5 m-1 m</td><td>0.5 m-1 m</td><td></td></tr></tbody></table>

Remarks: regardless on the source, the different surface materials have to be grouped in order to be linked to their respective spectral properties found within the DART database. For more information on LOD visit: https://osmbuildings.org/blog/2018-02-28_level_of_detail/ accessed on 31 January 2021.

<table><caption>**Table A2.** Estimated Computation Times for One Time-Step Simulations.</caption><thead><tr><th>Scale</th><th>Scene Dimensions (m)</th><th>Cell Size (m)</th><th>Shortwave (DART)</th><th>Longwave (DART)</th><th>Tmrt Calculation (Python Code)</th><th>Total</th></tr></thead><tbody><tr><td>Tree</td><td>30 × 30 × 20</td><td>0.5 (x,y,z)</td><td>9'</td><td>100'</td><td>2'</td><td>111'</td></tr><tr><td>Building</td><td>100 × 100 × 20</td><td>1 (x,y,z)</td><td>12'</td><td>138'</td><td>3'</td><td>153'</td></tr><tr><td>Neighbourhood</td><td>1000 × 1000 × 30</td><td>2 (x,y) 1 (z)</td><td>450'</td><td>5175'</td><td>112'</td><td>5737'</td></tr></tbody></table>

* PC characteristics: 128 GB RAM, CPU Intel®Xeon®E5-1650 v3 @ 3.50GHz. 6 Cores, 12 Logical processors.

<table><caption>**Table A3.** Characteristics of Remote Sensing Inputs.</caption><thead><tr><th>Remote Sensing Products</th><th>Provider</th><th>Characteristics</th><th>Reference</th></tr></thead><tbody><tr><td>WorldView2 imagery</td><td>Digital Globe</td><td>Ten high-resolution WorldView2 images (8-band, 2 m pixel size). Used to derive land cover map.</td><td>Dissegna et al. 2019</td></tr><tr><td>Building footprint data</td><td>Openstreetmap</td><td></td><td></td></tr></tbody></table>
<table><caption>Table A3. *Cont.*</caption><thead><tr><th>Remote Sensing Products</th><th>Provider</th><th>Characteristics</th><th>Reference</th></tr></thead><tbody><tr><td>Digital surface model (DSM)</td><td>AW3D</td><td>Derived from photogrammetric reconstruction of Digital Globe stereo satellite images. Nominal resolution of the DSM dataset is 1m for the urban areas and 5m for the densely vegetated area. Used to derive a digital terrain model (DTM), buildings and vegetation heights.</td><td></td></tr><tr><td>2014 Aerial LiDAR of Singapore</td><td>Singapore Land Authority</td><td>Obtained using the Optech Pegasus HA500 sensor with a planar density of 30 points/m².</td><td>Urech et al. 2020</td></tr><tr><td>Terrestrial LiDAR Scanning (TLS)</td><td>SEC/SMART field campaign</td><td>Field measurements of leaf area density derived from TLS pointcloud collected using Leica P40 sensor.</td><td>Wei et al. 2020</td></tr></tbody></table>

## Appendix B

<table><caption>Table A4. DART Parameters for Longwave Simulations.</caption><tbody><tr><td colspan="4"><strong>Flux-tracking</strong></td></tr><tr><td></td><td>Radiation-Radiative method</td><td>Flux Tracking</td><td></td></tr><tr><td></td><td>Atmosphere radiative transfer TOA<->BOA</td><td>Radiative transfer simulation</td><td></td></tr><tr><td>Flux tracking parameters</td><td></td><td></td><td></td></tr><tr><td>Products</td><td></td><td></td><td>Sparse voxel acceleration: NO</td></tr><tr><td>Spectral interval</td><td>Spectral band</td><td>Longwave</td><td>Radiation budget products only<br>Central Wavelength 25, Bandwidth 50</td></tr><tr><td colspan="4"><strong>Products</strong></td></tr><tr><td></td><td>Radiative Budget</td><td>...</td><td>3D INTR, ABS, SCAR<br>W/m²</td></tr><tr><td></td><td></td><td>Radiative budget unit</td><td>Irradiance—Six directions</td></tr><tr><td></td><td></td><td>Radiative budget components</td><td></td></tr><tr><td></td><td></td><td>Scene element components</td><td></td></tr><tr><td colspan="4"><strong>Direction input parameters</strong></td></tr><tr><td></td><td>Sun angles or date</td><td></td><td>Exact date</td></tr><tr><td></td><td>Year</td><td></td><td>20xx</td></tr><tr><td></td><td>Month</td><td></td><td>xx</td></tr><tr><td></td><td>Day</td><td></td><td>xx</td></tr><tr><td></td><td>Hours</td><td></td><td>* To be set in sequence launcher tool</td></tr><tr><td></td><td>Local time</td><td></td><td>yes</td></tr><tr><td></td><td>Time Zone (UTC)</td><td></td><td>8 (For Singapore)</td></tr><tr><td colspan="4"><strong>Optical and temperatures properties</strong></td></tr><tr><td>Lambertian</td><td>Building walls</td><td>Property name<br>2D lambertian database<br>2D lambertian model</td><td>Building walls<br>Lambertian_mineral.db<br>Sand white</td></tr><tr><td>Lambertian</td><td>Glass</td><td>Property name<br>2D lambertian database<br>2D lambertian model</td><td>Glass<br>Lambertian_mineral.db<br>glass_windows</td></tr><tr><td>Lambertian</td><td>Aluminium Roof</td><td>Property name<br>2D lambertian database<br>2D lambertian model</td><td>Aluminium roof<br>Lambertian_mineral.db<br>aluminium_windows</td></tr><tr><td>Lambertian</td><td>Badminton Court</td><td>Property name<br>2D lambertian database<br>2D lambertian model</td><td>Badminton Court<br>Lambertian_mineral.db<br>concrete_brown</td></tr></tbody></table>
<table><caption>**Table A4. Cont.**</caption><tbody><tr><td rowspan="3">Lambertian</td><td rowspan="3">Stem</td><td>Property name</td><td>Stem</td></tr><tr><td>2D lambertian database</td><td>Lambertian_mineral.db</td></tr><tr><td>2D lambertian model</td><td>Bark_deciduous</td></tr><tr><td rowspan="3">Lambertian</td><td rowspan="3">Grass</td><td>Property name</td><td>Grass</td></tr><tr><td>2D lambertian database</td><td>Lambertian_mineral.db</td></tr><tr><td>2D lambertian model</td><td>Grass_rye</td></tr><tr><td rowspan="3">Vegetation</td><td rowspan="3">Tree leaf</td><td>Vegetation property name</td><td>Tree leaf</td></tr><tr><td>3D vegetation database</td><td>Lambertian_Vegetation.bd</td></tr><tr><td>3D vegetation model</td><td>leaf_deciduous</td></tr><tr><td>Temperature</td><td>Thermal function</td><td>Temperature (Mean and Delta)</td><td>Hourly mean and deltas obtained from energy balance model for each scene component (Table A6).</td></tr><tr><td rowspan="5">Atmosphere</td><td>Gas optical properties</td><td>Tropical</td></tr><tr><td>Gas temperature profile</td><td>Tropical</td></tr><tr><td>Gas O<sub>3</sub> and other gases</td><td>Tropical</td></tr><tr><td>vertical profile</td><td></td></tr><tr><td>Redefine temperature profile</td><td>✓</td></tr><tr><td>Altitude</td><td>1.5 m</td></tr><tr><td rowspan="3">Aerosol</td><td rowspan="3">Aerosol properties</td><td>Temperature [K]</td><td>Specific temperature at each time step (Table A6)</td></tr><tr><td>Aerosol optical properties</td><td>Tropical_urbav5</td></tr><tr><td>Aerosol vertical profile</td><td>Urbanav5</td></tr><tr><td></td><td></td><td>Aerosol Henyey Greenstein parameters</td><td>Urbanav5</td></tr></tbody></table>

<table><caption>**Table A5. DART Parameters for Shortwave Simulations.**</caption><tbody><tr><td colspan="3"><strong>Flux-Tracking</strong></td><td>Radiation-Radiative method</td><td>Flux Tracking</td></tr><tr><td colspan="3">Flux-tracking parameters</td><td>Atmosphere radiative transfer TOA&larr;BOA</td><td>Analytic model</td></tr><tr><td colspan="3">Products</td><td></td><td>Sparse voxel acceleration: NO</td></tr><tr><td colspan="3">Spectral interval</td><td></td><td>Radiation budget products only</td></tr><tr><td colspan="3"></td><td>Spectral band</td><td>Central Wavelength 1.4, Bandwidth 2.2</td></tr><tr><td colspan="3"><strong>Products</strong></td><td>Shortwave</td><td></td></tr><tr><td colspan="3"></td><td>Radiative Budget</td><td>3D INTR, ABS, SCAR<br>W/m<sup>2</sup></td></tr><tr><td colspan="3"></td><td>...</td><td>Irradiance—Six directions</td></tr><tr><td colspan="3"></td><td>Radiative budget unit</td><td></td></tr><tr><td colspan="3"></td><td>Radiative budget components</td><td></td></tr><tr><td colspan="3"></td><td>Scene element components</td><td></td></tr><tr><td colspan="3"><strong>Direction input parameters</strong></td><td></td><td></td></tr><tr><td colspan="3"></td><td>Sun angles or date</td><td>Exact date</td></tr><tr><td colspan="3"></td><td>Year</td><td>20xx</td></tr><tr><td colspan="3"></td><td>Month</td><td>xx</td></tr><tr><td colspan="3"></td><td>Day</td><td>xx</td></tr><tr><td colspan="3"></td><td>Hours</td><td>* To be set in sequence launcher tool</td></tr><tr><td colspan="3"></td><td>Local time</td><td>yes</td></tr><tr><td colspan="3"></td><td>Time Zone (UTC)</td><td>8 (For Singapore)</td></tr><tr><td colspan="3"><strong>Optical and temperatures properties</strong></td><td></td><td></td></tr><tr><td rowspan="3">Lambertian</td><td rowspan="3">Building walls</td><td>Property name</td><td>Building walls</td></tr><tr><td>2D lambertian database</td><td>Lambertian_mineral.db</td></tr><tr><td>2D lambertian model</td><td>Sand white</td></tr></tbody></table>
<table><caption>**Table A5. Cont.**</caption><tbody><tr><td>Lambertian</td><td>Glass</td><td>Property name<br>2D lambertian database<br>2D lambertian model</td><td>Glass<br>Lambertian_mineral.db<br>glass_windows</td></tr><tr><td>Lambertian</td><td>Aluminium Roof</td><td>Property name<br>2D lambertian database<br>2D lambertian model</td><td>Aluminium roof<br>Lambertian_mineral.db<br>aluminium_windows</td></tr><tr><td>Lambertian</td><td>Badminton Court</td><td>Property name<br>2D lambertian database<br>2D lambertian model</td><td>Badminton Court<br>Lambertian_mineral.db<br>concrete_brown</td></tr><tr><td>Lambertian</td><td>Stem</td><td>Property name<br>2D lambertian database<br>2D lambertian model</td><td>Stem<br>Lambertian_mineral.db<br>Bark_deciduous</td></tr><tr><td>Lambertian</td><td>Grass</td><td>Property name<br>2D lambertian database<br>2D lambertian model</td><td>Grass<br>Lambertian_mineral.db<br>Grass_rye</td></tr><tr><td>Vegetation</td><td>Tree leaf</td><td>Vegetation property name<br>3D vegetation database<br>3D vegetation model</td><td>Tree leaf<br>Lambertian_Vegetation.bd<br>leaf_deciduous</td></tr><tr><td>Temperature</td><td>Thermal function</td><td>Mean temperature<br>Delta temperature</td><td>Hourly mean and deltas obtained from energy<br>balance model for each scene component (Table A6).</td></tr><tr><td colspan="4"><strong>Atmosphere</strong></td></tr><tr><td></td><td></td><td>Gas optical properties<br>Gas temperature profile<br>Gas O<sub>3</sub> and other gases<br>vertical profile</td><td>Tropical<br>Tropical<br>Tropical</td></tr><tr><td>Aerosol</td><td>Aerosol properties</td><td>Aerosol optical properties<br>Aerosol vertical profile<br>Aerosol Henyey Greenstein<br>parameters<br>AOD multiplicative factor</td><td>Tropical_urbav5<br>Urbanv5<br>Urbanv5<br>Retrieved from Aeronet for each time step of the<br>simulation date (Table A7).</td></tr></tbody></table>

## Appendix C

<table><caption>**Table A6. DART Parameters for Longwave Simulations.**</caption><thead><tr><th colspan="15">Longwave Sequence Parameters</th></tr><tr><th></th><th>seq_0</th><th>seq_1</th><th>seq_2</th><th>seq_3</th><th>seq_4</th><th>seq_5</th><th>seq_6</th><th>seq_7</th><th>seq_8</th><th>seq_9</th><th>seq_10</th><th>seq_11</th><th>seq_12</th></tr></thead><tbody><tr><td>Sequence name</td><td>7:00</td><td>8:00</td><td>9:00</td><td>10:00</td><td>11:00</td><td>12:00</td><td>13:00</td><td>14:00</td><td>15:00</td><td>16:00</td><td>17:00</td><td>18:00</td><td>19:00</td></tr><tr><td>Time</td><td>7:00</td><td>8:00</td><td>9:00</td><td>10:00</td><td>11:00</td><td>12:00</td><td>13:00</td><td>14:00</td><td>15:00</td><td>16:00</td><td>17:00</td><td>18:00</td><td>19:00</td></tr><tr><td>Badminton court</td><td>299.15</td><td>300.65</td><td>304.65</td><td>309.15</td><td>314.15</td><td>314.15</td><td>315.15</td><td>318.15</td><td>316.65</td><td>313.15</td><td>312.15</td><td>310.65</td><td>307.65</td></tr><tr><td>Delta</td><td>0</td><td>1</td><td>7</td><td>14</td><td>20</td><td>20</td><td>22</td><td>26</td><td>25</td><td>18</td><td>16</td><td>15</td><td>11</td></tr><tr><td>Glass windows</td><td>299.15</td><td>301.15</td><td>303.65</td><td>306.15</td><td>309.15</td><td>307.65</td><td>307.15</td><td>308.15</td><td>308.65</td><td>309.65</td><td>308.65</td><td>306.15</td><td>303.15</td></tr><tr><td>Delta</td><td>0</td><td>2</td><td>5</td><td>8</td><td>12</td><td>7</td><td>4</td><td>6</td><td>9</td><td>11</td><td>9</td><td>6</td><td>0</td></tr><tr><td>Walls</td><td>299.15</td><td>300.65</td><td>302.15</td><td>308.15</td><td>313.15</td><td>312.65</td><td>311.15</td><td>316.15</td><td>316.65</td><td>312.65</td><td>311.65</td><td>310.15</td><td>307.65</td></tr><tr><td>Delta</td><td>0</td><td>1</td><td>2</td><td>12</td><td>18</td><td>17</td><td>14</td><td>22</td><td>25</td><td>17</td><td>15</td><td>14</td><td>11</td></tr><tr><td>Aluminum roof</td><td>299.15</td><td>300.65</td><td>302.15</td><td>303.65</td><td>306.15</td><td>308.65</td><td>311.65</td><td>312.15</td><td>311.15</td><td>310.15</td><td>308.15</td><td>305.65</td><td>302.65</td></tr><tr><td>Delta</td><td>0</td><td>1</td><td>2</td><td>3</td><td>6</td><td>9</td><td>13</td><td>14</td><td>14</td><td>12</td><td>8</td><td>5</td><td>1</td></tr><tr><td>Grass</td><td>299.15</td><td>301.15</td><td>304.15</td><td>305.15</td><td>306.15</td><td>307.15</td><td>308.15</td><td>307.15</td><td>307.15</td><td>306.15</td><td>305.15</td><td>304.15</td><td>302.15</td></tr><tr><td>Delta</td><td>0</td><td>2</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>4</td><td>2</td><td>0</td></tr><tr><td>Leaves</td><td>299.15</td><td>301.15</td><td>304.15</td><td>305.15</td><td>306.15</td><td>307.15</td><td>308.15</td><td>307.15</td><td>307.15</td><td>306.15</td><td>305.15</td><td>304.15</td><td>302.15</td></tr><tr><td>Delta</td><td>0</td><td>2</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>4</td><td>2</td><td>0</td></tr><tr><td>Trunk</td><td>299.15</td><td>301.15</td><td>303.15</td><td>304.15</td><td>306.15</td><td>305.15</td><td>304.15</td><td>304.15</td><td>305.15</td><td>306.15</td><td>305.15</td><td>304.15</td><td>302.15</td></tr><tr><td>Delta</td><td>0</td><td>2</td><td>4</td><td>4</td><td>6</td><td>4</td><td>2</td><td>2</td><td>4</td><td>6</td><td>4</td><td>2</td><td>0</td></tr><tr><td>Air temperature *</td><td>299.51</td><td>300.06</td><td>301.01</td><td>301.48</td><td>303.5</td><td>303.75</td><td>304.55</td><td>306.28</td><td>305.31</td><td>305.08</td><td>304.58</td><td>304.25</td><td>303.15</td></tr></tbody></table>

<footnote>* Measured from a Vaisala weather station mounted next to radiometers.</footnote>
<table><caption>**Table A7.** DART Parameters for Shortwave Simulations.</caption><thead><tr><th colspan="14">Shortwave Sequence Parameters</th></tr><tr><th>Sequence name</th><th>seq_0</th><th>seq_1</th><th>seq_2</th><th>seq_3</th><th>seq_4</th><th>seq_5</th><th>seq_6</th><th>seq_7</th><th>seq_8</th><th>seq_9</th><th>seq_10</th><th>seq_11</th><th>seq_12</th></tr></thead><tbody><tr><td>Time</td><td>7:00</td><td>8:00</td><td>9:00</td><td>10:00</td><td>11:00</td><td>12:00</td><td>13:00</td><td>14:00</td><td>15:00</td><td>16:00</td><td>17:00</td><td>18:00</td><td>19:00</td></tr><tr><td>AOD multiplicative factor *</td><td>0.171</td><td>0.172</td><td>0.185</td><td>0.192</td><td>0.361</td><td>0.422</td><td>0.337</td><td>0.340</td><td>0.297</td><td>0.211</td><td>0.214</td><td>0.158</td><td>0.185</td></tr></tbody></table>

## Appendix D

<figure><img src="image_13.png" /><figcaption>**Figure A1.** Spectral Signatures of Construction Materials and Vegetation (from DART Database).</figcaption></figure>

<footnote>* Obtained from NASA AERONET for Singapore for simulation date 29 February 2020.</footnote>
## Appendix E. Sensitivity to Changing LAI

<table><thead><tr><th colspan="4">Sensitivity to Changing LAI at 16:00</th></tr><tr><th rowspan="2">LAI</th><th colspan="3">Average Tmrt under canopy [°C]</th></tr><tr><th>Khaya S.</th><th>Albizia S.</th><th>Tabebuia R.</th></tr></thead><tbody><tr><td>0.1</td><td>33.15</td><td>48.81</td><td>49.34</td></tr><tr><td>0.5</td><td>30.27</td><td>44.33</td><td>45.89</td></tr><tr><td>1.5</td><td>28.25</td><td>38.05</td><td>40.46</td></tr><tr><td>3.5</td><td>26.84</td><td>33.78</td><td>34.4</td></tr><tr><td>5.5</td><td>25.89</td><td>30.32</td><td>31.7</td></tr></tbody><tfoot><tr><td colspan="4">Tmrt outside canopy shade is 64.1 °C</td></tr></tfoot></table>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /><figcaption>**Figure A2.** (a) Sensitivity to Changing LAI at 16:00 (b) Visualization of Tmrt at pedestrian level under different tree species with changing LAI.</figcaption></figure>
# Appendix F. Sensitivity to Changing Surface Temperatures of Different Ground Material

<table><caption>**Table A8.** Sensitivity to Changing Surface Temperatures of Different Ground Materials at 16:00.</caption><caption>Sensitivity to Changing Surface Temperatures of Different Ground Materials at 16:00</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="7">312.65 K</th></tr><tr><th>20 °C</th><th>25 °C</th><th>30 °C</th><th>35 °C</th><th>40 °C</th><th>50 °C</th><th>60 °C</th></tr></thead><tbody><tr><td rowspan="7">80 m × 80 m scene with<br>one building<br>(W:10 m × L:20 m × H:15 m)</td><td>Surface temperature concrete building</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Surface temperature ground</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Surface temperature ground (K)</td><td>293.15</td><td>298.15</td><td>303.15</td><td>308.15</td><td>313.15</td><td>323.15</td><td>333.15</td></tr><tr><td>Tmrt (°C) Grass</td><td>57.53</td><td>59.37</td><td>61.20</td><td>63.08</td><td>65.03</td><td>*</td><td>*</td></tr><tr><td>Tmrt (°C) Concrete</td><td>60.01</td><td>61.67</td><td>63.38</td><td>65.16</td><td>66.99</td><td>70.84</td><td>74.91</td></tr><tr><td>Tmrt (°C) Asphalt Light</td><td>59.75</td><td>61.47</td><td>63.24</td><td>65.08</td><td>66.98</td><td>70.95</td><td>75.16</td></tr><tr><td>Tmrt (°C) Asphalt Grey</td><td>57.29</td><td>59.00</td><td>60.78</td><td>62.62</td><td>64.51</td><td>68.49</td><td>72.69</td></tr><tr><td></td><td>Tmrt (°C) Asphalt Dark</td><td>52.60</td><td>54.51</td><td>56.49</td><td>58.62</td><td>60.73</td><td>65.12</td><td>69.75</td></tr></tbody></table>

<figure><img src="image_16.png" /><figcaption>**Figure A3.** Sensitivity to Changing Surface Temperatures of Different Ground Material.</figcaption></figure>

# Appendix G. Sensitivity to Changing Wall Material

<table><caption>**Table A9.** Sensitivity to Changing Wall Material (16:00).</caption><caption>Sensitivity to Changing Wall Material (16:00)</caption><tbody><tr><td>Surface temperature ground * (Asphalt road light)</td><td>303.15 K Delta: 8</td></tr><tr><td>Surface temperature white concrete wall and roof *</td><td>303.15 K Delta: 8</td></tr><tr><td>Surface temperature glass *</td><td>303.15 K Delta: 8</td></tr><tr><td>Temperature range</td><td>26–34 °C</td></tr><tr><td>Tmrt at 0.5 m from sunlit white concrete wall</td><td>79.39 °C</td></tr><tr><td>Tmrt at 0.5 m from sunlit white concrete—glass wall</td><td>68.59 °C</td></tr><tr><td>Tmrt at 0.5 m from sunlit glass wall</td><td>64.47 °C</td></tr></tbody></table>

<footnote>* Omitted Tmrt values since a surface temperature above 40 °C for grass is unrealistic.</footnote>

<footnote>* Assumed surface temperature.</footnote>
<figure><img src="image_17.png" /><figcaption>**Figure A4.** (a,b) Sensitivity to Changing Wall Material. (c) Plot showing the values corresponding to the cells comprised in the black rectangles 1, 2, 3 on each of the three buildings.</figcaption></figure>

## References

1. Fanger, P.O. *Thermal Comfort*; Danish Technical Press: Copenhagen, Denmark, 1970.

2. Thorsson, S.; Rocklöv, J.; Konarska, J.; Lindberg, F.; Holmer, B.; Dousset, B.; Rayner, D. Mean radiant temperature—A predictor of heat related mortality. *Urban Clim.* **2014**, *10*, 332–345. [CrossRef]

3. Höppe, P. The physiological equivalent temperature—A universal index for the biometeorological assessment of the thermal environment. *Int. J. Biometeorol.* **1999**, *2466*, 71–75. [CrossRef] [PubMed]

4. Jendritzky, G.; De Dear, R.; Havenith, G. Institutional Repository UTCI—Why another thermal index? *Int. J. Biometeorol.* **2012**, *56*, 421–428. [CrossRef] [PubMed]

5. Staiger, H.; Laschewski, G.; Grätz, A. The perceived temperature—A versatile index for the assessment of the human thermal environment. Part A: Scientific basics. *Int. J. Biometeorol.* **2012**, *56*, 165–176. [CrossRef]

6. Golasi, I.; Salata, F.; de Lieto Vollaro, E.; Coppi, M. Complying with the demand of standardization in outdoor thermal comfort: A first approach to the Global Outdoor Comfort Index (GOCI). *Build. Environ.* **2018**, *130*, 104–119. [CrossRef]
7. Lai, D.; Lian, Z.; Liu, W.; Guo, C.; Liu, W.; Liu, K.; Chen, Q. A comprehensive review of thermal comfort studies in urban open spaces. *Sci. Total Environ.* **2020**, *742*, 140092. [CrossRef]

8. Modest, M.F. *Radiative Heat Transfer*; McGraw-Hill: New York, NY, USA, 1993.

9. Matzarakis, A.; Rutz, F.; Mayer, H. Modelling radiation fluxes in simple and complex environments: Basics of the RayMan model. *Int. J. Biometeorol.* **2010**, *54*, 131–139. [CrossRef]

10. Verein Deutscher Ingenieure. *Environmental Meteorology Methods for the Human Biometeorological Evaluation of Climate and Air Quality for Urban and Regional Planning at Regional Level. Part 1 Climate*; VDI 3787; Verlag des Vereins Deutscher Ingenieure: Düsseldorf, Germany, 2008; ISBN 5662422001.

11. Hoppe, P.; Mayer, H. Planungsrelevante Bewertung der thermischen Komponente des Stadtklimas. *Landsch. Stadt* **1987**, *19*, 22–30.

12. Bruse, M.; Huttner, S. Numerical modeling of the urban climate—A preview on ENVI-MET 4.0. *Int. Conf. Urban Clim.* **2009**, *29*, 1–4.

13. Lindberg, F.; Holmer, B.; Thorsson, S. SOLWEIG 1.0—Modelling spatial variations of 3D radiant fluxes and mean radiant temperature in complex urban settings. *Int. J. Biometeorol.* **2008**, *52*, 697–713. [CrossRef]

14. Krayenhoff, E.S.; Voogt, J.A. A microscale three-dimensional urban energy balance model for studying surface temperatures. *Bound.-Layer Meteorol.* **2007**, *123*, 433–461. [CrossRef]

15. Huang, J.; Cedeño-Laurent, J.G.; Spengler, J.D. CityComfort+: A simulation-based method for predicting mean radiant temperature in dense urban areas. *Build. Environ.* **2014**, *80*, 84–95. [CrossRef]

16. Nice, K.A.; Coutts, A.M.; Tapper, N.J. Urban climate development of the VTUF-3D v1: 0 urban micro-climate model to support assessment of urban vegetation in fluences on human thermal comfort. *Urban Clim.* **2018**, *24*, 1052–1076. [CrossRef]

17. Mitraka, Z.; Chrysoulakis, N.; Doxani, G.; Del Frate, F.; Berger, M. Urban surface temperature time series estimation at the local scale by spatial-spectral unmixing of satellite observations. *Remote Sens.* **2015**, *4139–4156*. [CrossRef]

18. Landier, L.; Gastellu-Etchegorry, J.P.; Al Bitar, A.; Chavanon, E.; Lauret, N.; Feigenwinter, C.; Mitraka, Z.; Chrysoulakis, N. Calibration of urban canopies albedo and 3D shortwave radiative budget using remote-sensing data and the DART model. *Eur. J. Remote Sens.* **2018**, *51*, 739–753. [CrossRef]

19. Chrysoulakis, N.; Grimmond, S.; Feigenwinter, C.; Lindberg, F.; Gastellu-Etchegorry, J.P.; Marconcini, M.; Mitraka, Z.; Stagakis, S.; Crawford, B.; Olofson, F.; et al. Urban energy exchanges monitoring from space. *Sci. Rep.* **2018**, *8*, 1–8. [CrossRef]

20. Morrison, W.; Yin, T.; Lauret, N.; Guilleux, J.; Kotthaus, S.; Gastellu-etchegorry, J.; Norford, L.; Grimmond, S. Remote sensing of environment atmospheric and emissivity corrections for ground-based thermography using 3D radiative transfer modelling. *Remote Sens. Environ.* **2020**, *237*, 111524. [CrossRef]

21. Dissegna, M.A.; Yin, T.; Wei, S.; Richards, D.; Grêt-Regamey, A. 3-D reconstruction of an urban landscape to assess the influence of vegetation in the radiative budget. *Forests* **2019**, *10*, 700. [CrossRef]

22. Landier, L.; Al Bitar, A.; Lauret, N.; Gastellu-Etchegorry, J.P.; Aubert, S.; Mitraka, Z.; Feigenwinter, C.; Parlow, E.; Heldens, W.; Kotthaus, S.; et al. 3D modeling of radiative transfer and energy balance in urban canopies combined to remote sensing acquisitions. *Int. Geosci. Remote Sens. Symp.* **2016**, *2016*, 6738–6741.

23. Chen, Y.; Chen, C.; Matzarakis, A.; Liu, J.; Lin, T. Modeling of mean radiant temperature based on comparison of airborne remote sensing data with surface measured data. *Atmos. Res.* **2016**, *174–175*, 151–159. [CrossRef]

24. Gastellu-Etchegorry, J. Discrete Anisotropic Radiative Transfer (DART 5) for modeling airborne and satellite spectroradiometer and LIDAR acquisitions of natural and urban landscapes. *Remote Sens.* **2015**, *7*, 1667–1701. [CrossRef]

25. Gastellu-Etchegorry, J.P.; Martin, E.; Gascon, F. DART: A 3D model for simulating satellite images and studying surface radiation budget. *Int. J. Remote Sens.* **2004**, *25*, 73–96. [CrossRef]

26. Gastellu-Etchegorry, J. Modeling radiative transfer in heterogeneous 3D vegetation canopies. *Remote Sens. Environ.* **1996**, *58*, 131–156. [CrossRef]

27. Widlowski, J.L.; Taberner, M.; Pinty, B.; Bruniquel-Pinel, V.; Disney, M.; Fernandes, R.; Gastellu-Etchegorry, J.P.; Gobron, N.; Kuusk, A.; Lavergne, T.; et al. Third Radiation Transfer Model Intercomparison (RAMI) exercise: Documenting progress in canopy reflectance models. *J. Geophys. Res. Atmos.* **2007**, *112*, 1–28. [CrossRef]

28. Widlowski, J.; Pinty, B.; Lopatka, M.; Atzberger, C.; Buzica, D.; Chelle, M.; Disney, M.; Gerboles, M.; Gobron, N.; Grau, E.; et al. The fourth radiation transfer model intercomparison (RAMI-IV): Proficiency testing of canopy reflectance models with ISO-13528. *J. Geophys. Res. Atmos* **2013**, *118*, 6869–6890. [CrossRef]

29. Urech, P.R.W.; Dissegna, M.A.; Girot, C.; Grêt-Regamey, A. Point cloud modeling as a bridge between landscape design and planning. *Landsc. Urban Plan.* **2020**, *203*, 103903. [CrossRef]

30. Béland, M.; Widlowski, J.L.; Fournier, R.A.; Côté, J.F.; Verstraete, M.M. Estimating leaf area distribution in savanna trees from terrestrial LiDAR measurements. *Agric. For. Meteorol.* **2011**, *151*, 1252–1266. [CrossRef]

31. Béland, M.; Widlowski, J.L.; Fournier, R.A. A model for deriving voxel-level tree leaf area density estimates from ground-based LiDAR. *Environ. Model. Softw.* **2014**, *51*, 184–189. [CrossRef]

32. Richardson, J.J.; Moskal, L.M.; Kim, S. Modeling approaches to estimate effective leaf area index from aerial discrete-return LIDAR. *Agric. For. Meteorol.* **2009**, *149*, 1152–1160. [CrossRef]

33. Zheng, G.; Ma, L.; Eitel, J.U.H.; He, W.; Magney, T.S.; Moskal, L.M.; Li, M. Retrieving directional gap fraction, extinction coefficient, and effective leaf area index by incorporating scan angle information from discrete aerial lidar Data. *IEEE Trans. Geosci. Remote Sens.* **2017**, *55*, 577–590. [CrossRef]
34. Singapore Land Authority. 2014 *Aerial LiDAR of Singapore*; SLA: Singapore, 2014.

35. Wei, S.; Yin, T.; Dissegna, M.A.; Whittle, A.J.; Ow, G.L.F.; Yusof, M.L.M.; Lauret, N.; Gastellu-Etchegorry, J.P. An assessment study of three indirect methods for estimating leaf area density and leaf area index of individual trees. *Agric. For. Meteorol.* **2020**, *292–293*, 108101. [CrossRef]

36. Morrison, W.; Kotthaus, S.; Grimmond, C.S.B.; Inagaki, A.; Yin, T.; Gastellu-Etchegorry, J.P.; Kanda, M.; Merchant, C.J. A novel method to obtain three-dimensional urban surface temperature from ground-based thermography. *Remote Sens. Environ.* **2018**, *215*, 268–283. [CrossRef]

37. Meili, N.; Manoli, G.; Burlando, P.; Bou-Zeid, E.; Chow, W.T.L.; Coutts, A.M.; Daly, E.; Nice, K.A.; Roth, M.; Tapper, N.J.; et al. An urban ecohydrological model to quantify the effect of vegetation on urban climate and hydrology (UT&C v1.0). *Geosci. Model Dev.* **2020**, *13*, 335–362.

38. Holmer, B.; Lindberg, F.; Rayner, D.; Thorsson, S. How to transform the standing man from a box to a cylinder—A modified methodology to calculate mean radiant temperature in field studies and models. In Proceedings of the 9th International Conference on Urban Climate (ICUC9), Toulouse, France, 28 May 2015; pp. 20–24.

39. Dirmhirn, I. *Das Strahlungsfeld im Lebensraum*; Akademische Verlagsgesellschaft: Frankfurt, Germany, 1964.

40. Höppe, P. *Die Energiebilanz des Menschen*; Universität Munchen, Meteorologisches Institut: Munich, Germany, 1984; Volume 49.

41. Höppe, P. Ein neues Verfahren zur Bestimmung der mittleren Strahlungstemperatur in Freien. *Wetter Leb.* **1992**, *44*, 147–151.

42. Kipp and Zonen. *Instruction Manual for Kipp & Zonen CNR 4 Net Radiometer*; Kipp and Zonen: Delft, The Netherlands, 2008.

43. Oyj, V. *Vaisala User Guide Vaisala Weather Transmitter WXT530 Series*; Vaisala: Helsinki, Finland, 2017.

44. Campbell Scientific. *BlackGlobe Temperature Sensor for Heat Stress User Manual*; Campbell Scientific: Logan, UT, USA, 2015.

45. Wang, Y.; Grimaldi, J.; Landier, L.; Chavanon, E.; Gastellu-Etchegorry, J.P. Introduction of clouds in DART model. *Int. Arch. Photogramm.* **2020**, *XLIII*, 843–848.

46. Kotthaus, S.; Smith, T.E.L.; Wooster, M.J.; Grimmond, C.S.B. ISPRS Journal of Photogrammetry and Remote Sensing Derivation of an urban materials spectral library through emittance and reflectance spectroscopy. *ISPRS J. Photogramm. Remote Sens.* **2014**, *94*, 194–212. [CrossRef]

47. Gastellu-Etchegorry, J.P. 3D modeling of satellite spectral images, radiation budget and energy budget of urban landscapes. *Meteorol. Atmos. Phys.* **2008**, *102*, 187–207. [CrossRef]