Numerical Validation of the Radiative Model for the Solar Cadaster Developed for Greater Geneva

: The achievement of the targets for reducing greenhouse gas emissions set by the Paris Agreements and the Swiss federal law on the reduction of greenhouse gas emissions (CO 2 law) requires massive use of renewable energies, which cannot be achieved without their adoption by the general public. The solar cadaster developed as part of the INTERREG G2 Solar project is intended to assess the solar potential of buildings at the scale of Greater Geneva—for both industrial buildings and for individual residential buildings—at a resolution of 1m. The new version of the solar cadaster is intended to assess the solar potential of roofs, as well as that of vertical facades. The study presented here aims to validate this new version through a comparison with results obtained with two other simulation tools that are widely used and validated by the scientiﬁc community. The good accordance with the results obtained with ENVI-met and DIVA-for-Rhino demonstrates the capability of the radiative model developed for the solar cadaster of Greater Geneva to accurately predict the radiation levels of building facades in conﬁgurations with randomly distributed buildings (horizontally or vertically).


Introduction
The massive use of renewable energies and, more particularly, solar energy is necessary in order to meet the objectives of the Paris Agreements and the Swiss federal law on the reduction of greenhouse gas emissions (CO 2 law) [1]. As buildings are some of the main contributors to climate change, one of the main challenges is, therefore, turning buildings from energy consumers into energy producers. This cannot be done without the development of tools dedicated to the evaluation of solar potential.
This massive deployment involves the evaluation of solar potential on a large scale. In addition, the urban environment has many specific features that are to be taken into consideration when evaluating this potential. This can be achieved with solar cadasters.
However, these tools are usually limited to the evaluation of the solar potential on rooftops [2][3][4][5]. This limitation can be an obstacle to the development of the use of solar energy because the potential of vertical surfaces in an urban environment can prove to be predominant [6][7][8][9].
Thus, there is a need for large-scale tools that are able to take vertical facades into account, in addition to roofs. From a physical point of view, the effects of surrounding surfaces (the ground or surrounding buildings) on the vertical radiative balance are not well known and may appear to be non-negligible. From the point of view of solar production potential, the issues of the variability linked to shading effects are added. Tools that are able to consider these issues are being developed and need validation.
A solar cadaster tool is being developed at the scale of Greater Geneva, which is a trans-border agglomeration around the city of Geneva between France and Switzerland, and it totals an area of 2000 km 2 . A first solar cadaster was developed for the Canton of Geneva only (280 km 2 ), as presented by Desthieux et al. [10]. The goal now is to update and extend this solar cadaster for the whole of the Greater Geneva area in the framework of INTERREG G2 Solar. Its aim is to intensify solar energy production in the agglomeration by mobilizing a large panel of stakeholders from public institutions, municipalities, energy providers, the private sector, and universities. The difficulty lies in the fact that the solar cadaster must cover a large area (2000 km 2 for Greater Geneva) while providing accurate information at a small spatial resolution and with a reasonable computation time.
The main objective of this study is to validate a solar cadaster tool by using a comparison with the results from other numerical tools. To this purpose, the results from the solar cadaster are compared with those of two widely used tools: DIVA-for-Rhino and ENVI-met. The latter, thanks to its holistic model, makes it possible to deepen the assessment of the buildings' photovoltaic potential by evaluating their surface temperatures.
The validation of the solar cadaster involves three different tools based on different models and methods and with a varied field of applications. Thus, these tools are presented first, prior to introducing the methodology used. The results in terms of predicted shortwave irradiance on both the roofs and the facades are then detailed and discussed.

Modeling Approaches for Radiative Transfers
Two main methods are commonly used for the evaluation of radiative transfers: raytracing-based models or radiosity methods. Their principles, as well as their advantages and drawbacks, are described in the following.

Ray-Tracing Methods
Ray-tracing methods consist of following the paths taken by electromagnetic rays. There are two types of ray-tracing approaches. Forward ray-tracing methods, also called light tracing, consist of launching rays from light sources in a set of directions. Backward ray-tracing methods consist of following the light path backward. In this latter case, rays are launched in a set of directions from the element of interest (the eye of an observer or an irradiated surface). A comparison between forward and backward ray tracing is illustrated in Figure 1.

Eye
Light source Incident Ray Reflected Ray  For both of these methods, the directions can be random (probabilistic approach) or set (deterministic approach). After rays have been launched, an intersection test is carried out between each of these rays and the objects in the scene; then, new rays are cast from the points of impact while considering the properties of the materials encountered (reflection, transmission, absorption). This method makes it possible to take the effects of light into account, such as refraction of light (for example, in the case of a scene comprising glasses). The behavior of a ray at the point of intersection then depends on the refractive properties of the material encountered. This feature is particularly useful in the context of image synthesis or the study of transparent media with different refractive indexes.
This method also has the advantage of making the calculation of specular reflections possible (see Figure 2). On the other hand, it is not suitable for a calculation of diffuse reflections (which are usually the case for urban surfaces), as these require a significant calculation time because of the number of rays to be launched (in all directions) at each new intersection. This method has to be coupled with a diffuse radiation model.

Radiosity Method
The radiosity method is based on the calculation of the balance of radiosity, which is the total radiative flux exiting a surface (W m −2 ). This method allows the calculation of diffuse and reflected components of the irradiation only by solving a system of equations involving all the objects of a scene. Taking specular reflections (see Figure 2) into account is not directly possible with this method. However, unlike ray tracing, the calculations are independent of the position of the possible observer (surface of interest)-the calculation is generalized over the whole scene.
The radiosity method is based on solving the equation introduced by Kajiya [11] and detailed in the work of Sillion and Puech [12], where the exiting power is expressed as: where x is a point in space, B(x) is the exiting power (emitted, reflected, and transmitted) at point x in space (in W m −2 ), E(x) is the power emitted by point x (source term, in W m −2 ), ρ is the coefficient of reflection of the surface, L i (x, θ, φ) is the incident power at point x from the solid angle determined by the angles θ and φ, Ω is the set of directions (θ, φ) of the hemisphere, and dω is the elementary solid angle (sr). This equation can be discretized: where B i and B j are the power exiting the i-th cell and j-th cell, respectively, or the radiosity of the i-th or j-th cell (in W m −2 ); this term takes into account both the term emitted and the effect of the reflections. E i is the power initially emitted by the i-th cell, ρ i is the reflection coefficient of the i-th cell, i.e., the fraction of the energy received that the cell returns (unitless value between 0 and 1), N is the number of meshes of which the scene is composed, and F ij is the form factor between the i-th and the j-th cells.
Solving the radiosity equation requires one to first know the source term E i (known if the source emits directly), the reflexivity ρ i , which is generally known as one of the characteristics of the materials, and the matrix of the form factors F ij , which must be calculated because it is dependent on the studied configuration.

Tools Considered
Three different tools were considered in this study: the solar cadaster (CadSol) for Greater Geneva (which is the one to be validated) and two widely validated tools-DIVAfor-Rhino (DIVA) and ENVI-met (EM). Their models and applications are described in the following.

Solar Cadaster
As summarized by Desthieux et al. [10] and Freitas et al. [13], the solar cadaster tool is a geographic information system (GIS) tool. Such tools enable one to process large amounts of data and spatial analyses. These tools also provide automatic or systematic environmental analyses of urban areas, such as solar radiation calculations. These are different from tools that are classified in the category of computer-aided design (CAD). The latter process more accurate spatial and weather data in terms of spatial and time scales, but they require much more computing time and, thus, address the local scale (limited sets of buildings).
The irradiance received from the direct component I b on an element of a surface at location x at a given time t is given by: in which BN I corresponds to the direct normal component of the irradiance (also called the beam), r b is the transposition factor, and S b corresponds to the shadow cast from a neighborhood building (S b (x, t) = 0 if, at time t, the surface located at x is shaded by another building, or S b (x, t) = 1 otherwise). The transposition factor, r b , depends on the solar elevation h and the slope of the considered surface β (β = 0 corresponds to a horizontal surface and β = 90 • corresponds to a vertical surface; see Figure 3). It is calculated as r b (x, t) = sin(h(t))/cos(β(x)). The irradiance over the element of the surface is then considered as the average of the calculated irradiances at its edges. The calculation differences between the three tools, as described below, are illustrated in Figure 4. The diffuse component, I d , of the solar radiation is then calculated as follows: where GH I is the global irradiation on a horizontal surface, DH I is the diffuse irradiation on a horizontal surface, I 0 is the hourly extraterrestrial irradiation, r b is the transposition factor as defined in (3), and SVF is the sky-view factor. Finally, the reflected component is simply considered at the current stage as isotropic and is estimated based on Iqbal (1983) [15] as follows: where GH I is as defined in (4), ρ is the coefficient of reflection of the surface, as defined in (1), and β is the slope of the surface. More details about the calculation method and the computing tools used are available in [10,16].

DIVA-for-Rhino
DIVA-for-Rhino is a highly optimized daylighting and energy modeling plug-in for Rhinoceros. This software uses ray-tracing and light-backwards algorithms based on the physical behavior of light in a 3D volumetric model. For hourly solar radiation, the Daysim interface is used. Daysim has been validated by several studies to be accurate in modeling visible-wavelength natural light for multiple sky conditions [17].
The daylight coefficient approach and the all-weather sky luminance model according to [18] are used here. In this approach, the irradiance received on an element of surface x is calculated as the sum of all sky segments visible from this element of the surface. For the diffuse component, 145 sky segments are used [19] concomitantly with three ground segments [20]. For the direct component, Daysim uses 65 sun positions; therefore, at a specific time t, Daysim will use one of the 65 positions that is closest to the real sun position at time t.
To calculate the complete set of daylight coefficients, two ray-tracing runs are performed: The default Daysim simulation parameters were chosen. Up to two reflections from direct solar irradiation and one reflection from diffuse sky irradiation from the environment were considered.
Similarly to the solar cadaster, the irradiance calculated by Diva-for-Rhino for each surface element corresponds to the average of those calculated at its vertices.

ENVI-Met
ENVI-met is a software aiming at simulating the urban microclimate by taking into consideration all of the phenomena that occur in an urban environment. It is based on coupled balance equations (including those of mass, momentum, and energy). This involves taking the built and natural environment into account.
The ENVI-met model has been widely used in numerous studies dealing with different issues, including the evaluation of outdoor thermal comfort [21,22], mitigation of the urban heat island effect [22,23], or assessment of buildings' solar and photovoltaic potential [24].
Regarding radiative transfers, ENVI-met uses a hybrid approach based on a radiosity method (see Section 2.1.2) for the evaluation of the irradiance with a deterministic raytracing method for the calculation of the view factors. As ENVI-met is under a proprietary license, access to its code is limited. However, some information is available. Regarding the calculation of the diffuse solar radiation, the model is isotropic. This means that there is no difference between the different sky segments and no dependency on the actual position of the sun. Concerning the sky, it is divided into 414 segments (207 upward and 207 downward) [25].
Unlike the two other tools, the irradiances predicted by ENVI-met are given as that calculated at the center of the mesh.

Case Study
The present study focuses on two fictitious districts composed of three rows and three columns of buildings, i.e., a total of nine buildings. This ability of this kind of fictitious district has already been proven with respect to the study of the solar potential of neighborhoods [26][27][28]. In addition, it makes it possible to study the same district at different location, depending on the weather conditions used as input for the simulations.
These two districts, which are called the homogeneous and the heterogeneous district, have different arrangements, but they share some common urban morphology indicators, which are listed in Table 1. Regarding the thermo-radiative properties of the different surfaces, the coefficient of reflection is set to 0.2 for both the ground's and the buildings' envelopes (vertical facades and roofs). These values are taken as a compromise between the typical values for concrete and soil.

Study of the Influence of the Mesh Resolution on the Results
The mesh sensitivity was tested on a simple configuration-presented in Figure 5which consisted of two buildings that were 25 and 35 m high and separated by 10 m. Three grid sizes of 0.5, 1, and 2 m square meshes were chosen. The daily cumulative irradiation on the south facade of the highest building is given in Figure 6. It appears that the cumulative irradiation predicted by the solar cadaster increased along with the coarseness of the resolution.  The influence of the mesh resolution on the results is detailed in Table 2 for the months of February and August. The mesh resolution had an influence on the irradiation predicted for every tool considered. Indeed, the coarser the resolution, the higher the discrepancy in terms of cumulative irradiation. Nonetheless, DIVA, which is based on a ray-tracing method, was less influenced by the mesh resolution, unlike ENVI-met and the solar cadaster, which are based on radiosity methods. On the other hand, the influence of the mesh resolution was more important in February than in August. This was due to the sun's lower course, which led to more shadowing.
In the particular case of the solar cadaster tool, a resolution of 2 m is, therefore, coarse and is not recommended for an analysis of irradiation in an urban environment. A resolution of 1 m remains acceptable for large areas, and the deviation from 0.5 m is small. This difference is explained by the fact that the obstacles' mutual shadows between the two plots will be better detected and considered in every facade point with a finer resolution.
A mesh composed of square cells with a resolution of 1 m was retained for the following study. This constituted a good compromise between the accuracy of the results and the computation time.

Homogeneous Neighborhood
The homogeneous neighborhood was composed of nine identical buildings, which were 20 m wide and 30 m high, and each was separated by 20 m, as shown in Figure 7. In this case, the widths of the streets were the same, and so was the height-to-width ratio ( Table 3) for all of the buildings. This configuration also made it possible to have the same sky-view factor regardless of the orientation of the vertical facade, as well as a sky-view factor that was equal to 0.5 for the roofs. The homogeneous neighborhood provided a simple configuration and made it easier to study the different radiative phenomena that occurred on the buildings' facades.

Heterogeneous Neighborhood
As a counterpart used to make the study of radiative phenomena easier, the homogeneous neighborhood (see Figure 7) could not be fully representative of an actual district, as actual districts are generally composed of buildings of different heights or random positions. To bridge this gap, a heterogeneous neighborhood was considered as well. Its layout is given in Figure 8.
This district is representative of a common type of actual district because of the non-homogeneous spatial distribution of the buildings, as well as their different heights. The buildings' footprints and their total volume were kept constant, as were the urban morphology indicators given in Table 1. The urban morphology indicators specific to the heterogeneous district are given in Table 4.

Weather Conditions for the Study
The present study focused on two different days: 15 February, which had a short daylight period, and 16 August, which had a longer one. These two days are representative of a winter and a summer day for a location such as Geneva (cold and not very sunny for the first, hot and more sunny for the second). Furthermore, the sun's path corresponding to these two days was halfway between the solstices and the equinoxes.
The two days were actually those for which the sun's path was closest to the mean paths for the months. Indeed, in order to reduce the number of simulated days (two days to be simulated for the months of February and August instead of 59) while avoiding the specific conditions of a given day, this study considered two representative average days for the two considered months [10]. A representative average day (RAD) is defined, according to Equation (6), as a day for which the meteorological conditions (including irradiation level, temperature, and wind velocity and direction) for each hour are equal to the average of these conditions over all days of the month (N day ). Carrying out the study for two different months allowed us to evaluate the influence of the irradiation level on the accuracy of the predicted results.
where < · · · > stands for the time-average operator. The two considered months had a double advantage. First, they allowed the validation of the values predicted by the solar cadaster under low and high irradiation levels (February and August, respectively). Second, they made it possible to evaluate the influence of the evolution of the sun's path over the year on the irradiance profiles.
Since the solar cadaster was developed for Greater Geneva, the fictitious districts studied were located at this place. The data used as input for the simulations came from the METEONORM ® version 7.3 database for this city. The hourly averaged values issued from the METEONORM database were then averaged by month according to Equation (6) in order to reduce the computation time. Since the values used as input for the simulations were averaged by hour, the level of irradiation was considered constant between two consecutive hours for both the input and output data.
In this study, only the shortwave radiation (SW) was considered, as it is the most important energy content in the solar spectrum. Thus, the waveband considered was 2500 nm. The daily profiles of the different components of the solar radiation over Geneva for the representative days of the months of February and August are given in Figure 9. The direct shortwave radiation represents the amount of energy that comes straight from the sun, while the diffuse part is the amount of solar radiation reflected by the atmosphere prior to hitting the ground. The total shortwave radiation is the sum of the direct and the diffuse parts.
Although the overall evolution of the shortwave radiation is similar between February and August, differences are noticeable. Indeed, the maximum in terms of GHI is two times higher in August. This is mainly due to the direct part of the solar radiation, which exceeds 350 W m −2 in August, in comparison with 150 W m −2 in February. 16

Results
The validation of the solar cadaster radiation model is a multi-step process. In the first step, the focus is on the horizontal surfaces of the homogeneous district because they are subject to a smaller number of phenomena. In the second step, the mean irradiance values over the vertical facades of the central building (see building 5 in Figures 7 and 8) are studied. Finally, the spatial distribution of the irradiance is analyzed.

Mean Irradiance over the Unshaded Horizontal Surface
The mean irradiance level over the roofs is given in Figure 10. For the homogeneous district (see Figure 7), the irradiance level is the same for the roofs of the nine buildings, since they are not subject to shading effects. Regarding the heterogeneous district, the irradiance level is the that of buildings 6 and 7 (see Figure 8), which are the highest buildings.
This comparison demonstrates the ability of the solar cadaster to accurately reproduce the basics of the solar conditions (including the beam horizontal irradiance (BHI) and the diffuse horizontal irradiance (DHI)). This first preliminary analysis allows us to ensure that the results of the modeling of the global horizontal irradiance (GHI) are the same for all tools. In other words, the difference that will be observed in what follows will be the consequence of the presence of shading and the inclination of the facades.

Mean Irradiance over the Vertical Facades
The goal of this study is to demonstrate the ability of the solar cadaster to provide accurate results for vertical facades. The results for the homogeneous district are given in Figures 11 and 12, while the results for the heterogeneous district are given in Figures 13 and 14. It appears that the values predicted by the solar cadaster are in good accordance with the values predicted by DIVA-for-Rhino and ENVI-met. The different orientations of the facades make it possible to evaluate the prediction of the irradiance level for different solar conditions. Indeed, the north facade is mainly irradiated by the diffuse and reflected parts of the solar radiation, while the proportion of direct radiation is higher for the south facade. The east facade is directly irradiated in the morning, while the west facade faces the sun later in the day.
Regarding the level of irradiation, it does not have an influence on the concordance of the predicted values. Indeed, the results for the summer day are, in general, as accurate as those for the winter day. The differences between the results of the tools-quantified here with the normalized root mean squared error (NRMSE)-are not significantly impacted by the month of the year considered (see Table 5). The NRMSE values are calculated from differences between the daily profiles of the spatial average irradiance received on the facades and those predicted by the solar cadaster and ENVI-met on the one hand or by the solar cadaster and DIVA-for-Rhino on the other hand.
In both cases, the NRMSE remains lower for the roofs (less than 1 % for the roofs versus 4 % to 25 % for the vertical facades; see Table 5). This demonstrates the good accordance between the three tools with respect to the evaluation of the incident radiation on the roof, which are without a mask or reflection here. The observed differences for the vertical facades, which result in an increase in the NRMSE in the last four lines in Table 5, are then due to a difference in the consideration of masks and inter-building reflections. Nevertheless, the time of the year does not influence only the level of irradiance, but also the sun's path. Indeed, regarding the north facade (see , the results from the solar cadaster are close to the those from DIVA-for-Rhino and ENVI-met (although a bit overestimated). Nonetheless, one peak can be observed in the results at 6:00 a.m. in the month of August (see Figure 12), which is not present for the same facade in February (see Figure 11). This difference in terms of the profile over time for the irradiance level is due to the evolution of the sun's path over the year. Indeed, during the month of August, the sun rises in the north-east and sets in the north-west, while it rises in the east and sets in the west in February.
The sun's path has an influence on the level of irradiation as well. Indeed, in the month of February (Figure 11), the south facade shows a slight decrease around noon, which is not present in the GHI (see Figure 9). This decrease is actually due to the shadow cast by building 2 on the central building (see Figures 7 and 8). This does not occur in August because the sun's path is high enough for the shading effect not to occur.
The results for the heterogeneous district make it possible to evaluate the ability of the solar cadaster to accurately predict the irradiation level for a more complex city. Indeed, the buildings are randomly located in this case (see Figure 8). The results for the two months considered are given in Figures 13 and 14.
It appears that the morphology of the neighborhood, although it is more random and complex than that of the homogeneous neighborhood, does not have a significant negative impact on the accuracy of the results predicted by the solar cadaster, which are in good accordance with those of DIVA-for-Rhino and ENVI-met (see Tables 5 and 6). Nonetheless, the level of solar irradiance over the different facades is lower in the case of the heterogeneous district than in the case of the homogeneous one. This can be explained by the fact that the central building is smaller than those around it and is thus subject to more of a shading effect.

Irradiance Maps
Although the results presented so far show a good accordance between the three tools, a consideration of the spatial mean value alone is not sufficient. Indeed, as shown in Figures 15 and 16, the range of the level of irradiance received over the facade can be very important.
The range of the predicted values may be due to the method of modeling radiative phenomena (ray-tracing or radiosity) or the shading effect. The latter may have an important impact on the level of irradiance. Thus, the study of the mean irradiance over the facades needs to be complemented with the study of irradiance maps. The level of irradiance on unshaded roofs appears to be perfectly homogeneous and equal for the three tools considered, and it corresponds to the sum of the BHI and the DHI when given as the input of the simulation.
Nevertheless, the level of irradiance on the vertical facades depends on the orientation of the facade, the time of the day, and the time of the year. The levels of irradiance for the south facade of the central building of the homogeneous district are given in Figure 17. Regarding the data, not all are used for comparisons. All of the facades of the buildings are 20 m wide. However, the data at the edges of facades were not analyzed. Indeed, these data are difficult to compare because of the difference in terms of the calculation method for the three tools considered (depending on the tool, irradiance is considered as that at the center of the mesh or as the average over it); see Figure 4. This source of difference, due to the schemes specific to the tools, does not lead to a significant difference in the results, except at the edges of the facades.
Nevertheless, the level of irradiance on the vertical facades depends on the orientation of the facade, the time of the day and the time of the year. The level of irradiance for the south facade of the central building of the homogeneous district are given in Figure 17. Regarding the data, not all are used for comparisons. All the facades of the buildings are 20 m wide. However, the data at the edges of facades were not analyzed. Indeed, these data are difficult to compare because of the difference in terms of calculation method for the three tools considered (depending on the tool, irradiance is considered as that at the center of the mesh or as the average on it), see Figure 4. This source of difference, due to the schemes specific to the tools, does not lead to a significant difference in the results, except at the edges of the facades. The irradiance levels appears more homogeneous on the south facade for the month of August than for the month of February. This can be explained by a higher sun's path at this time of the year. The higher range in terms of value (see Figures 15 and 16) is then due the shading effect, more intense in February, which leads to a range of irradiance over the south facade of 50 W m −2 in August against 350 W m −2 in February.
Regarding the level of irradiance predicted by the solar cadaster, shown in Figure 17a,d, its values are similar to those predicted by ENVI-met and DIVA-for-Rhino. The irradiance levels appear to be more homogeneous on the south facade for the month of August than for the month of February. This can be explained by the sun's higher path at this time of year. The greater range in terms of value (see Figures 15 and 16) is then due to the shading effect, which is more intense in February and leads to a range of irradiance of 50 W m −2 over the south facade in August in comparison with 350 W m −2 in February.
Regarding the level of irradiance predicted by the solar cadaster, which is shown in Figure 17a,d, the values are similar to those predicted by ENVI-met and DIVA-for-Rhino.
Nevertheless, although the predicted irradiance levels are close to each other, there are two particular points. First, the shape induced by the shading is not exactly the same for the solar cadaster and ENVI-met on the one hand and for DIVA-for-Rhino on the other hand. The latter predicts a smoother evolution of the predicted irradiance level. This smoothness is due to the method used for the evaluation of the radiative phenomena; DIVA-for-Rhino uses ray-tracing while ENVI-met and the solar cadaster use the radiosity method. The second point is how the irradiance is evaluated. ENVI-met considers the irradiance to be that at center of the mesh, while DIVA-for-Rhino and the solar cadaster consider it to be the average of the irradiance over the mesh. In the case of DIVA-for-Rhino, it reinforces the smoothness of the evolution of the irradiance over the facade. Concerning the solar cadaster, this has an impact at the edge of the shade (the horizontal purple line at the height of 17 m in Figure 17a).

Discussion
The results presented so far show a general good accordance between the three tools considered in terms of both the mean values the irradiance maps. Nevertheless, two points must be noted: 1. The prediction of the irradiance is highly sensitive around sunrise and sunset; 2. A slight time offset between the results of the different tools seems to be present.
Regarding the sensitivity of the predicted irradiance, the highest discrepancies between the results of the tools appear at the first or the last hour of sunshine (see the north and the west facades in Figure 12). These discrepancies are due to a lack of precision between the calculation of the sun's path and the level of irradiation. This shift has only a slight incidence on the prediction of the irradiance during the day, but sees its influence strongly increase when the sun is close to the horizon. A sunrise that is predicted too early in the morning results in an overestimation of the predicted irradiance at this time of day, while a sunset predicted too late in the evening results in a peak of the predicted irradiance before the dusk. Figure 17 shows the existence of a time lag between the three tools studied. Indeed, the shapes of the shadow on the south facade are different among the results of the solar cadaster (Figure 17a), ENVI-met (Figure 17b), and DIVA-for-Rhino (Figure 17c). This offset has only a slight influence on the results. Nonetheless, this impact may increase incrementally, especially early in the morning or late in the evening, as mentioned before (see the north facade in Figure 12).

Conclusions
The solar cadaster has proven its ability to provide accurate results in terms of the level of irradiance on both roofs and vertical facades. Indeed, the results predicted by the new version of the solar cadaster developed for Greater Geneva are in good accordance with those predicted by ENVI-met and DIVA-for-Rhino.
Regarding the concordance of the results, the non-shaded horizontal surfaces match almost perfectly. This is not the case for the shaded facades and vertical facades, but the results are very satisfactory. The differences observed between the tools come from the differences between the models in terms of diffuse radiation, the evaluation of reflections, or shading.
Finally, the solar cadaster provides precise results at a low spatial resolution while keeping the computation time low. It ranges from a few hours for the solar cadaster and DIVA-for-Rhino to a few days for ENVI-met. It can therefore be used on a large scale and can prove to be a reliable tool for promoting and intensifying the use of solar energy on the scale of Greater Geneva. However, work is in progress to improve the modeling of reflected components and, thus, the reliability of the tool for more complex geometries.

Data Availability Statement:
The data that supports the findings of this study are available from the corresponding autho upon reasonable request.