Next Article in Journal
Room Temperature NH3 Selective Gas Sensors Based on Double-Shell Hierarchical SnO2@polyaniline Composites
Next Article in Special Issue
Stingray Sensor System for Persistent Survey of the GEO Belt
Previous Article in Journal
Effective Signal Extraction Algorithm for Cerebral Blood Oxygen Based on Dual Detectors
Previous Article in Special Issue
Classification of Low Earth Orbit (LEO) Resident Space Objects’ (RSO) Light Curves Using a Support Vector Machine (SVM) and Long Short-Term Memory (LSTM)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Validation of a 3D Local-Scale Adaptive Solar Radiation Model by Using Pyranometer Measurements and a High-Resolution Digital Elevation Model

by
Eduardo Rodríguez
1,†,
Judit García-Ferrero
2,3,†,
María Sánchez-Aparicio
4,†,
José M. Iglesias
5,†,
Albert Oliver-Serra
1,†,
M. Jesús Santos
2,3,†,
Paula Andrés-Anaya
4,†,
J. Manuel Cascón
2,6,†,
Gustavo Montero García
1,†,
Alejandro Medina
2,3,†,
Susana Lagüela
4,†,
M. Isabel Asensio
2,5,*,† and
Rafael Montenegro Armas
1,†
1
Instituto Universitario de Sistemas Inteligentes y Aplicaciones Numéricas en Ingeniería (SIANI), Universidad de Las Palmas de Gran Canaria, 35017 Las Palmas de Gran Canaria, Spain
2
Instituto Universitario de Física Fundamental y Matematicas (IUFFyM), Universidad de Salamanca, 37008 Salamanca, Spain
3
Departamento de Física Aplicada, Universidad de Salamanca, 37008 Salamanca, Spain
4
Departamento de Ingeniería Cartográfica y del Terreno, Universidad de Salamanca, 37008 Salamanca, Spain
5
Departamento de Matemática Aplicada, Universidad de Salamanca, 37008 Salamanca, Spain
6
Departamento de Economía e Historia Económica, Universidad de Salamanca, 37007 Salamanca, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2024, 24(6), 1823; https://doi.org/10.3390/s24061823
Submission received: 21 February 2024 / Revised: 6 March 2024 / Accepted: 8 March 2024 / Published: 12 March 2024
(This article belongs to the Special Issue Optical Sensors for Space Situational Awareness)

Abstract

:
The result of the multidisciplinary collaboration of researchers from different areas of knowledge to validate a solar radiation model is presented. The MAPsol is a 3D local-scale adaptive solar radiation model that allows us to estimate direct, diffuse, and reflected irradiance for clear sky conditions. The model includes the adaptation of the mesh to complex orography and albedo, and considers the shadows cast by the terrain and buildings. The surface mesh generation is based on surface refinement, smoothing and parameterization techniques and allows the generation of high-quality adapted meshes with a reasonable number of elements. Another key aspect of the paper is the generation of a high-resolution digital elevation model (DEM). This high-resolution DEM is constructed from LiDAR data, and its resolution is two times more accurate than the publicly available DEMs. The validation process uses direct and global solar irradiance data obtained from pyranometers at the University of Salamanca located in an urban area affected by systematic shading from nearby buildings. This work provides an efficient protocol for studying solar resources, with particular emphasis on areas of complex orography and dense buildings where shadows can potentially make solar energy production facilities less efficient.

1. Introduction

Knowledge of local solar radiation is essential for designing optimal solar energy systems. Ideally, this knowledge should come from long-term measured data. However, the limited coverage of radiation measurement networks imposes the need to develop solar radiation models. Therefore, validation of these models against measurements is crucial [1].
Dealing with solar radiation and shadows is a common issue in Geographical Information Systems (GIS). This stems from a wide range of interests that go from visualization and rendering to many potential applications such as predicting the radiation at crucial environments like glaciers [2], the design of agroforestry exploitations [3] or, in the urban environment, project planning [4], the generation of solar radiation maps [5] essential for the deployment of concentrating solar power (CSP) plants [6], and for planning the optimal location and orientation of PV or solar thermal panels for distributed generation of electricity or heat for domestic applications on building roofs or facades [7,8,9,10,11,12,13,14,15,16]. The latter application demands reliable models for solar radiation and shadow casting at the local level that can deal with the complexities of urban architecture, which has a tremendous impact on blocking the solar radiation of nearby areas. Furthermore, this involves a high accuracy of the input so that it can guarantee to capture the mentioned architectural features. On top of all this, a contained computation cost is highly desirable.
The Sun’s radiation in the direct, diffuse and reflected forms, as well as the treatment of shadow casting, is an active matter of research for which several methodologies have been proposed. Li and Liu [10] employed the Yallop algorithm to calculate the solar position and evaluate several methods for calculating diffuse radiation, but they did not address the issue of shadows cast by nearby obstacles. In [15], Toledo et al. report on the determination of the optimal PV panel tilting angle accounting for direct and diffuse radiation calculated by using the Erbs’s correlation and the effect of shadows cast by surrounding buildings parameterized with their height and distance. In [7,8], Dorman et al. proposed an open-source tool devised for modeling and shadow casting of a 2.5D vector-based environment. The shadow calculation is based on trigonometric relations between the Sun’s rays, the ground and any intervening obstacles. In [16], Tripathy et al. report on a system for computing the shadows as the intersection of parallelograms that arise from the transformation of a polygon-defined 3D city into a cloud of points with the terrain, facades and horizontal surfaces. Stendardo and coworkers [12] employed separate digital surface/elevation models with different accuracy levels for the terrain and for a 3D city in a system powered by ray-tracing algorithms implemented through parallelized GPU computations in order to estimate the solar potential in Geneva. Liang et al. [13] proposes a GIS-integrated open source tool that includes the r.sun package for the estimation of clear-sky, diffuse and reflected solar radiation that is in turn coupled to a well-established OpenGL-based 3D-graphics engine to compute shadow casting. In [14], Toledo and coworkers employed the SOL algorithm coupled to the direct and diffuse solar radiation model from Perez [17]. The algorithm, implemented in MATLAB, employs the concepts of hyperpoints to calculate the height of shadows in facades in a rectangular grid that represents the 3D city. In [4], Kaynak et al. obtain the direct solar radiation with the elevation angle constant (EAC) method and use the local meteorological with the Anstrom–Prescott model to estimate the sunlight attenuation due to atmospheric conditions, while the finite element method, back-face detection and ray-tracing algorithms are used to make a precise calculation of cast shadows.
Regarding the accuracy level of the Digital Elevation Model (DEM), several works suggest that a finely detailed 3D representation of the geometric features is vital to accurately evaluate the shadows cast by nearby geometries. In [18], the authors used various GIS light projection models to analyze how the level of detail in a building’s 3D approximation affects the quality of its projected shadows. They found that finer details are relevant for local applications such as predicting solar panels’ potential. In [19], the solar yield was assessed based on the grid spacing of the point cloud, finding out that the error increases linearly with the grid spacing up to 4 m, which is generally larger for structured grids than for unstructured grids, and that considerably fine grids are necessary when considering roofs with complex shadowing artifacts.
In the context of GIS and the representation of architectural features in the urban environment, a common approach is the use of the so-called 2.5D polygons, which consist mainly of a simplified representation of buildings as a composition of extruded shapes. The previous studies suggest that the polygon-based 2.5D urban mapping strategies might not be able to hold the needed detail when buildings that feature more complex shapes come into play. However, some works have extended the 2.5D geometry concept through the adaptation of more complex architectural shapes by using LiDAR data to reproduce the building geometry by means of a composition of arbitrarily and vertically (walls) oriented polygons within an adaptive mesh [20]. The present work uses a similar approach; we generate a high-definition DEM that accurately captures the terrain and the architectural geometry. This DEM is created from LiDAR data with a scan spacing of 0.30 m. Computing the shadows with this high-definition DEM is too computationally expensive, so instead, we use a triangular mesh adapted to it.
In this paper, we introduce a toolbox that accurately calculates solar radiation while considering the shadow casting of complex orography and architectural features. Section 2 describes this toolbox: the MAPSol model, the construction of the high-definition DEM and the generation of the adapted triangular mesh. In Section 3, we explain the process for acquiring experimental data and how we validated our model. We use a case study where shadows originate from the intricate geometry of the tower of the Cathedral of Salamanca onto a nearby placed pyranometer. The study’s results confirm our model’s accuracy in accounting for the shadows that impact solar radiation.

2. Materials and Methods

This multidisciplinary work has required a broad methodology covering different areas of knowledge. It begins by recalling some basic concepts about solar irradiance that will allow a better understanding of the work. Then, the MAPSol 3D local-scale adaptive solar irradiance model for which the validation protocol described in this paper has been designed is presented. This model allows the detection of the shadows cast by complex orography, but this requires a high-resolution map of the study surface. Therefore, the methodology used to obtain high-resolution DEM is described. To detect the shadows, the model relies on a mesh of the surface defined by the high-resolution DEM, so the next step is to generate an adapted mesh in the most efficient way. Finally, the model is validated by means of actual solar irradiance measurements. For this reason, the corresponding procedure of experimental data acquisition is described in detail.

2.1. Solar Irradiance

Irradiance is defined as the incident energy per unit of time and surface area, i.e., it is a power density [21]. The solar irradiance received on a point of interest depends on three key factors. The first factor is the Earth’s position concerning the Sun, determined by latitude, longitude, day, and time. This factor determines the angle at which the sunbeam strikes a given spot. The second factor is the orography and buildings, which impact shadows and the orientation between the sunbeam and the point of interest. The third factor is the atmosphere, which attenuates solar irradiance due to gases, solid and liquid particles, and clouds.
Depending on the factors considered, we can obtain two types of solar irradiance: real-sky irradiance (when all factors are considered) and clear-sky irradiance (when clouds are excluded) [22]. On the other hand, due to the scattering effects generated by the atmosphere and the effect of the angle of incidence, the global irradiance can be expressed as the sum of two components: the direct irradiance coming from the solar disk, and the diffuse irradiance, due to the scattering effects in the atmosphere.
The definition of the main irradiances is presented below [23].
  • Diffuse Horizontal Irradiance, D H I : This is the solar irradiance collected on a horizontal surface from the atmospheric scattering of light, excluding circumsolar radiation.
  • Direct Normal Irradiance, D N I : It is the component of solar irradiance collected on a surface perpendicular to the Sun’s rays. The horizontal diffuse component, D H I , is neglected here. On clear days, this component is much larger than the diffuse component, while on days with high cloud cover, it is practically zero. As it is measured over the Earth’s surface, its values depend highly on atmospheric conditions and the time of the year.
  • Global Horizontal Irradiance, G H I : This is the sum of all irradiance components collected over a horizontal surface. This includes the direct and diffuse components, as well as the reflected components, which are generally neglected because of their low value. The G H I can be calculated from the following expression:
    G H I = D H I + D N I · sin α
    where α is the solar altitude angle, i.e., the complementary of the zenith angle of the Sun, D H I is the horizontal diffuse component, and D N I is the normal component.
  • Beam Horizontal Irradiance, B H I : It is the direct horizontal component of the irradiance, i.e., the direct irradiance on a plane perpendicular to the vertical of the site. It can be obtained as follows:
    B H I = G H I D H I

2.2. The MAPSol Model

This section describes the solar irradiance model MAPSol developed at the University of Las Palmas de Gran Canaria [24,25,26]. MAPSol computes the clear-sky beam irradiance in every triangle of an adaptative mesh representing the orography surface. What sets MAPSol apart from other models is that it considers the orography shadows in the irradiance computation.

2.2.1. Clear-Sky Beam Irradiance Model

In this section, the main ingredients of the clear-sky beam irradiance model implemented in MAPSol are presented; a more detailed description can be found in [22,24].
The BHI equation for a horizontal surface is
B H I = I 0 ϵ 0.8662 T L K m δ R ( m ) L f sin α .
In this equation, each of the terms can be related to the three factors impacting solar irradiance.
The first term I 0 ϵ refers to the Earth’s geometry. I 0 is the solar constant (the Sunbeam irradiance at the mean solar distance) [27], and ϵ is a correction factor due to Earth’s orbit eccentricity [22,24].
The second term 0.8662 T L K m δ R ( m ) concerns the atmospheric attenuation. Specifically, m and δ R ( m ) are related to the attenuation by gas constituents, and T L K is related to the attenuation by solid and liquid particles [22]. The parameter m is the relative optical air mass [28]; T L K is the air mass 2 Linke atmospheric turbidity factor, and δ R ( m ) is the Rayleigh optical thickness at air mass m, and both are calculated according to the improved formula in [29].
Finally, the third term L f sin α is related to the shadows of both the orography and the buildings. Parameter α represents the solar altitude angle, while L f is the lighting factor that considers cast shadows. A Solar Position Algorithm computes the solar altitude angle; the present model uses the one developed at the Plataforma Solar de Almería: the PSA algorithm [30]. On the other hand, the light factor takes values between 0 and 1; L f = 0 means complete shadowing and L f = 1 total sun exposure. The following section describes the algorithm to compute this factor.

2.2.2. Shadow Detection

The model computes the light factor using a triangular mesh of the terrain generated with the method shown in Section 2.4. Therefore, L f relates directly to finding which triangles in the mesh are shadowed. This algorithm is described in detail in [24,26,31].
In short, this geometrical problem can be solved by finding those triangles that, looking at the mesh from the Sun, have at least another triangle that covers the former. This is done using a new reference system x , y and z , where z is in the direction of the beam irradiance (see Figure 1). Then, the mesh is projected on the plane x y . This is the view of the mesh from the Sun that allows us to find out which triangles are overlapping, if any.
To transform the coordinates and align the direction of the z axis with the beam irradiance, we need the solar azimuth θ and solar altitude α (shown in Figure 1). These angles are obtained from the PSA algorithm [30].
Examining intersected triangles at each time step and applying temporal refinement to each analyzed triangle, a shading level is allocated to shadowed triangles based on the count of warning points located within other triangles (refer to Figure 2). This allocation is the light factor L f .
The drawback of this technique emerges when attempting to apply it to a sizable mesh, as the computational cost becomes high, making the process too slow for some purposes.
In response to this challenge, a technique was devised in [31]. This approach involves pre-filtering candidate triangles before examining the ultimate mesh overlap. The determination of whether a triangle is shadowed is contingent upon either self-shadows or shadows cast by other triangles.
Self-shadows occur when a surface casts shadows onto itself, for instance, when it is oriented away from the Sun. Determining this condition is straightforward by examining the incidence angle ( δ e x p ), defined as the angle between the solar vector and the normal to a triangle surface. The criterion for self-shadows is met when the incidence angle exceeds π / 2 . So, a Self Shadows Light Factor ( L f s s ) is defined with values of
L f s s = 0 | δ e x p | > ( π / 2 ) 1 ( π / 2 ) | δ e x p | 0
As the analysis of cast shadows demands substantial computational effort, an initial step involves circumventing the examination of all triangles oriented away from the Sun (self-shadowed). An effective strategy involves filtering triangles based on the following:
  • In the absence of self-shadowed triangles (those facing away from the Sun), the entire mesh is illuminated, and no shadows are present.
  • Only triangles oriented away from the Sun are capable of casting shadows. These are referred to as potential 1 triangles [31].
The shading analysis is finalized by examining the shadow projection on a specific set of points known as warning points (WPs). These points are uniformly distributed within the triangle using the Rivara 4-T algorithm [32] and correspond to the geometric centers resulting from each subdivision of the triangle (see Figure 2).
A warning point is shaded when its projection onto the x y plane falls within a triangle that can potentially cast a shadow on it, and its z coordinate is smaller than that of its projection onto the said triangle (it is further from the sun than its projection on the triangle). So, the cast shadows light factor, L f c s is computed as follows:
L f c s = n w p i n w p
where i = 0 , 1 , , n w p is the number of shadowed warning points, and n w p is the total number of warning points. The ultimate light factor, L f is
L f = L f c s · L f s s

2.3. High-Resolution DEM

There is a wide variety of technologies that allow a massive capture of terrestrial information [33]. Among the different geospatial data coming from different sensors that allow the massive capture of information, LiDAR data stand out [34]. LiDAR data, obtained from laser scanners on board different devices, are highly accurate three-dimensional data. Aerial LiDAR data are the most widely used data for massive data acquisition in urban environments [35]. Thanks to the INSPIRE (Infrastructure for Spatial Information in Europe) initiative [36] based on the Directive 2007/2 [37], many countries within the European Union provide LiDAR data to users. The National Geographic Institute (IGN) in Spain offers highly accurate and updated LiDAR data of the entire Spanish territory, captured periodically under the National Plan for Territory Observation (Plan Nacional de Observación del Territorio—PNOT) [38]. Previous studies have demonstrated the applicability of these LiDAR data in urban environments in Spain [39,40,41,42].
Although the IGN allows the user to download the LiDAR point cloud through the Download Centre [43], it also offers several three-dimensional products that are derived from LiDAR data: (i) Digital Terrain Model (DTM), (ii) Digital Elevation Model (DEM), and (iii) Digital Slope Model (DSM). All of them have one feature in common: a spatial resolution between 2 and 25 m. Although a spatial resolution of meters may be sufficient for some studies, in certain cases, a resolution of less than one meter is necessary to obtain accurate results. For the purpose of this research study, a high-resolution DEM was generated using the available LiDAR data from IGN over the study area, described in Section 3.2. The LiDAR data used belong to the first coverage of the PNOT observation program and have the characteristics described in Table 1.
From this data, a high-resolution DEM is generated by applying the grid method of triangulation, called the TIN (Triangular Irregular Network), which converts the network of individual points in a continuous surface of triangular facets. The vertices of the triangles are the points measured by the LiDAR sensor, forming the point cloud. The procedure behind the generation of the DEM is the Delaunay triangulation, where the length of the edges and the angles for the triangles can be limited in order to generate smoother surfaces with fewer peaks, avoiding sliver triangles. Thus, according to the Delaunay triangulation, triangles are generated for groups of three points, ensuring that there is no point inside the circumsphere of any triangle. Once the triangulation is generated, the triangular faces are interpolated, using the values of elevation and slope in order to generate a smoother elevation layer.

2.4. Mesh Generation

The MAPSol model described in Section 2.2 requires a discretization of the orography provided by the high-resolution DEM described in Section 2.3. The uniform triangulation associated with the DEM-raster should be discarded, as it is not computationally feasible. As an alternative, the construction of an adapted mesh is proposed, which allows capturing the singularities of the study domain and computing the beam irradiance in a more efficient way.
The mesh generation is based on a local refinement driven by the error (vertical distance) between the approximation of the orography (adapted mesh) and the actual orography (DEM). The algorithm starts with a coarse mesh and locally refines [32] those triangles that do not satisfy the error criterion. It produces graded meshes that capture the features of the domain.
The process is fully automatic and only requires an initial discretization, in this case the DEM, to capture the orography with a fixed accuracy. It has been implemented with the Neptuno library developed by L. Ferragut [44]. The entire mesh generation process, including the DEM construction described in Section 2.3, is summarized in Figure 3.
Although MAPSol does not require a high-quality mesh, the previous procedure could be combined with surface parameterization and/or smoothing techniques [45,46] to improve the mesh quality and reduce the number of triangles of the final mesh.

2.5. Experimental Measurements of Solar Irradiance with Pyranometers

The purpose of this section is to obtain a realistic picture of solar radiation in the area of the Trilingüe building of the University of Salamanca (Spain) through experimental data. The physical principles of the pyranometers used to collect solar irradiance data are presented below.
The main function of pyranometers is to convert a flux of incident photons into an electrical signal through a process of collection, selection and detection [47]. Depending on the magnitude to be measured, they require different adaptations. Here, the focus is on devices for measuring G H I and D H I . In the case of GHI, the devices for measuring global radiation register very unequal fluxes, and they also have a strong dependence on the Sun position. In order for the electrical signal obtained to be representative of the incident flux, they have to mix the radiation inside, so that the resulting signal is isotropic. This can be achieved by using a diffuser sheet or an integrating sphere as receiving elements.
On the other hand, pyranometers dedicated to measuring DHI require a complement that eliminates the direct component of the radiation. The standard procedure for determining diffuse radiation is by using a pyranometer (see Figure 4 and paragraphs below) shadowed by a shadowing ring, made of an opaque material of a calibrated width that avoids direct radiation on the receiver. The ring is positioned so that its axis is parallel to the Earth’s axis in the equatorial coordinate system and passes through the receiver. Figure 4 schematizes the positioning of a shadow ring at a latitude of 40 ° N, where the arrows on the right side represent the Sun’s rays. Once the ring has been correctly positioned, the receiver will be permanently shaded provided that the ratio of the receiver’s distance from the center of the ring to its radius is close to the tangent of the Sun’s declination. Due to the changing declination of the Sun over the course of the year, the shadow ring requires small adjustments in its position every two or three days, aiming to keep the receiver permanently shaded.
To compensate for the defect in the diffuse irradiance collection generated by the shadow ring, a correction factor is introduced. It depends on the latitude of the observation point as well as the time of the year. In reference [48], the ring manufacturer provides the full derivation of the correction factor, C, as follows:
C = 1 1 S
where S is obtained as
S = 2 · 0.185 r a d Π cos δ ( U 0 · sin λ · sin δ + sin U 0 · cos λ · cos δ )
where λ is the angle of latitude, δ the angle of declination of the Sun, and U 0 is the angle formed by the Sun’s position at sunrise and noon in the plane of the ring.
U 0 = arccos ( tan λ tan δ )
There is a wide variety of pyranometers that implement different technologies depending on the type of transducer they use. They are subdivided into two classes: indirect conversion detectors and direct conversion detectors [47].
  • Indirect conversion detectors: They work by converting the incident photon flux into another type of flux (usually heat), but they can also be a secondary photon flux. Heat flux detectors are widely used and their operation is relatively simple. To convert the photon flux into heat flux, a highly absorbing paint or varnish is applied to the detector, which causes its temperature to rise when the light beam is impinging on it. Knowing the temperature at two points and assuming that the steady state is reached, the intensity of the flux is calculated, which will be proportional to the temperature difference. Figure 5a shows a general scheme of the parts of an indirect heat flux conversion pyranometer. In the upper part there are two domes, the outer dome has the function of avoiding energy exchanges due to convective phenomena; as a whole, the domes act as an integrating sphere. As can be seen, the detector is surrounded by an anti-radiation shield to prevent radiation penetrating from anywhere other than the dome. Figure 5b shows the Pyranometer Kipp and Zonen SMP10, belonging to the Energy Optimization, Thermodynamics and Statistical Physics Group (GTFE), with which the Global Horizontal Irradiance measurements were performed.
  • Direct conversion detectors: Again, there are two types. Photoemitter cells are based on the junction of an anode and a cathode, between which there is a large potential difference (in the range of kV), and an avalanche effect is produced. On the other hand, there are detectors based on PN junctions, the photodiodes, where the current generated is proportional to the incident flux. These types of detectors have better sensitivity than avalanche detectors and work with low voltage [49].
Figure 5. Heat flux sensing pyranometer: (a) basic scheme (taken from Kipp and Zonen. Instruction Manual SMP Series [50]) and (b) pyranometer belonging to the Group of Energy Optimization, Thermodynamics and Statistical Physics (GTFE) of the University of Salamanca.
Figure 5. Heat flux sensing pyranometer: (a) basic scheme (taken from Kipp and Zonen. Instruction Manual SMP Series [50]) and (b) pyranometer belonging to the Group of Energy Optimization, Thermodynamics and Statistical Physics (GTFE) of the University of Salamanca.
Sensors 24 01823 g005
For this work, radiation data obtained in situ are used. These data are recorded with two Kipp and Zonen pyranometers, from model SMP10. These are class A pyranometers (according to ISO 9060) that implement thermopiles as heat flux detectors. Thus, they work as indirect conversion detectors. With them, different magnitudes are collected: pyranometer 1 collects G H I measurements, while pyranometer 2, equipped with a Kipp and Zonen CM 121 shadow ring, takes D H I values. Table 2 shows the main characteristics of the model in question (detailed information about the pyranometers or the shadow ring can be found in the complete manufacturer’s manual [50]).
Figure 6 shows a comparison between the wavelength response of an SMP series pyranometer and the standard solar spectrum collected at sea level. As can be seen, the pyranometer response is practically linear over the entire range of interest, from the infrared (around 2000 nm) to the ultraviolet limit of the visible (around 350 nm). In this region of the spectrum, most of the energy from the Sun is deposited.

3. Results

In this section, the case study with which the MAPSol model has been validated is precisely described, including the description and reliability of the chosen irradiance data and the generation of the high-resolution DEM of the study area, as well as the adapted grid, and the simulation with the MAPSol model. The chosen case study provides the possibility to assess the effect of shadows caused by the intricate orography generated by buildings close to the location of pyranometers.

3.1. Experimental Data Acquisition

Experimental data were obtained from two pyranometers and a weather station, shown in Figure 7, installed on the rooftop of the Trilingüe building at the Faculty of Sciences (40.96062° N, 5.67075° W) of the University of Salamanca (Spain): a pyranometer Kipp and Zonen SMP10 (Figure 7a) and a pyranometer SMP10 with shadow ring model Kipp and Zonen CM 121 (Figure 7b).
During the day, the devices perform the corresponding measurements, storing G H I and D H I values every five minutes. At the end of the day, the collected data are sent to the computer, generating a csv file containing the irradiance information. The next step is to classify and filter the files. To do this, a program developed in “Mathematica® Software v. 13.1, Licensed to Universidad de Salamanca”, identifies the files that belong to the same day and extracts the desired information from each of them. Finally, the program generates a new output file where all the information for the same day is merged.
The results presented in this work have been obtained with data between 30 July 2020 and 22 February 2023. In Figure 8, illustrative monthly averages have been performed on the irradiance, so that for each month, G H I and B H I curves are obtained as a function of time (which is represented in the UTC format). The monthly average is obtained by considering the data of all the days belonging to that month in the whole study period.
As previously mentioned, pyranometers take measurements of G H I and D H I . To obtain the D H I component, a correction factor which depends on the day of the year is required. For each of the measurements, the correction factor is calculated as a function of the solar declination angle, and hereafter, D H I is considered to include the correction factor. The B H I and D H I measurements are obtained, every 5 min, from Equations (1) and (2), respectively. During the night, G H I and D H I values should be identically zero, but it is found that negative values appear. This is usual when performing this type of measurements. However, it does not make physical sense to work with negative irradiances. To solve this problem, the following filters are introduced into the program:
  • If D H I i < 0 D H I i = 0
  • If G H I i < 0 G H I i = 0
  • If ( G H I i D H I i ) < 0 ( G H I i D H I i ) = 0
  • If sin α i < 0 D N I i = 0
The irradiance curves introduced below are intended to be as close as possible to the actual measurements, so no statistical treatment beyond averaging is carried out. Figure 8 shows the irradiance curves for four representative months of the year: March, June, September and December. In all these plots, a certain level of noise is observed within the irradiances, which evinces that the pyranometers are very sensitive to small changes. For example, on a summer day with isolated clouds, if a cloud comes between the Sun and the pyranometer for a few minutes, the direct components will become very small. When the cloud has passed, the direct component will again reach high values. When observing the set of figures, it can be noticed that the irradiance curves increase in width and height in the summer months. On the contrary, during winter, when the Sun is up for less time and at a lower height above the horizon, the curves become narrower and present lower peak values.
In order to assess the reliability of the measurements obtained, these values are compared with two independent external sources: AEMET [51] and Solargis [52] for solar irradiance. The Atlas of Solar Radiation in Spain, using data from the EUMETSAT climate FAS developed by AEMET [51], uses satellite data with a spatial resolution of 3 × 3 km. Data from the period (1983–2005) were used to elaborate the atlas. In order to compare the results in a more convenient way, Table 3 shows the results of this work compared with those of AEMET. The relative differences obtained for 9 of the 12 months of the year do not overcome 10% in the daily accumulated energy for the global component ( G H I ). This number is reduced to eight for the direct horizontal component ( B H I ).
Taking the experimental data of G H I and B H I from Table 3, it is possible to obtain the annual cumulative values and also the average value per day of the energy received. The annual cumulative values are easily obtained by performing the following operation:
E x = i = 1 12 x i N i .
where x represents the irradiance component, the subscript i is the number of the month, and N i is the number of days in month i.
The Atlas of Solar Radiation in Spain [51] using data from the EUMETSAT climate FAS provides daily average G H I and B H I values in the form of an irradiance distribution map. The values extracted from these maps can be seen in Table 4, together with the values recorded in situ.
The value for G H I obtained in this work is within the AEMET interval, touching the upper limit. The experimental B H I component is 0.01 kWh m 2 · day 1 above the maximum value of the interval. AEMET does not offer the annual average accumulated values, but they can be calculated by multiplying the daily average by the number of days in the year (the result is also shown in the same table). Again, the experimental value of G H I is within the range, while the value of B H I is slightly above it.
As it is known, a Sun chart is employed to present, at a specific location, the apparent position of the Sun, i.e., the height of the Sun at any hour of the day. From a Sun chart elaborated by the University of Oregon Program [53], together with a panoramic photo from the Trilingüe building at the University of Salamanca, Figure 9 is elaborated. This figure shows the Sun path chart from the Trilingüe building at the University of Salamanca (40.96062° N, 5.670759° W) between 21 December and 21 June, overlapped with a panoramic photo, taken from a pyranometer that registers G H I to identify shadowing sources.
One of these shadowing sources on the location of the pyranometers is the 93-m-high cathedral tower, located 296 m to the east. As can be seen in Figure 9, the cathedral tower casts its shadow in this area during the early hours of some days of the year, specifically between 22 August and 11 September, and in its symmetrical months, February and March. In these months close to spring, cloudiness is higher, so it is more difficult to obtain irradiance measurements that allow the effect of the cathedral’s shadow to be clearly seen. However, in summer, the sky is clearer, and this phenomenon can be seen more clearly in the measurements obtained.

3.2. Area Study, High-Resolution DEM and Adapted Mesh

An area of 430 m × 176 m has been selected which includes the Trilingüe building at the Faculty of Sciences of the University of Salamanca, where the pyranometers are located, and the cathedral tower located at 296 m in a straight line to the east, which casts its shadow on the pyranometers affecting the reading of the irradiance data.
Provided with the point cloud available in the area under study, the DEM generated with the process described in Section 2.3 presents a resolution of 33 cm (Figure 10). This methodology performs the generation of the Delaunay triangulation using negligible computing time (3 s for 35,228 points and 1,385,800 triangles), with a computer with Intel Core i7-6700 processor at 3.41 GHz, 64 bits, 32 Gb. RAM and 931 Gb (Dell Precision Tower 3620, Round Rock, TX, USA).
The size of the raster file is too large for the calculations involved in the MAPSol model, so the mesh is adapted using the procedure described in Section 2.4, achieving a reduction of 62 % in the number of cells. The adapted mesh for a tolerance of 1 m has 529,503 triangles and includes all the singularities of the study area, in particular the cathedral tower. This element holds significant importance in the DEM. It is not present in some of the raster files, such as those provided by the IGN (Digital Surface Model–Building and Digital Surface Model), due to their low resolution and dimensions.
The number of triangles can be further reduced without affecting the accuracy of the shadow and irradiance calculations, using a tolerance of 5 m. This means that in the mesh adaptation process, a triangle is refined if the distance from the triangle to the high-resolution DEM is greater than 5 m, or 1 m in the previous case. The number of triangles of the adapted mesh with a tolerance of 5 m is 44,313 , which represents a reduction of more than 95 % with respect to the original mesh, without, as we will show further, a significant loss in precision in the calculation of shadows and irradiance, but with a considerably lower computational cost.
Figure 11c displays the fine adapted mesh (1 m) for the entire study area and several zooms of both meshes, the fine adapted mesh (Figure 11b) and the uniform original mesh (Figure 11a), for the Cathedral building area. A 3D reconstruction of the study area is depicted in Figure 12, including the location of the pyranometers with a red point, which allows the complexity of the study area to be observed. Furthermore, this procedure demonstrates that with few free resources, it is possible to reconstruct a 3D image of a very complex area such as this one, in the old part of the city of Salamanca. For Figure 12, the most appropriate orientation has been chosen so that the orthophoto of the area projected on the fine adapted mesh provides the best visual result.

3.3. Simulation with MAPSol

Using the coarse adapted mesh (5 m), the shadows and global irradiance in each triangle have been calculated for each day for a full year, with a time step of 5 min, with the MAPSol model. Computations have been performed in a computer equipped with two AMD EPYC 7313 CPUs, 128 GB of RAM memory and Debian Linux version 11 operating system. MAPSol is written in Python and C++, using the latter for the core of the computation and the former mostly to deal with input/output. In this particular work, Python 3.9.16 version and GNU g++ 10.2.1 C++ compiler have been used. Simulations have been performed in parallel, using one core per month. The mean wall clock execution time was 3 h and 56 min, including writing results to files for later analysis.
For each time step, a VTK file representing irradiance in the whole domain is written, as well as a csv file with values of GHI in the location of the pyranometers. Monthly average values of global irradiance in the study area have been calculated, and in Figure 13, four of the most significant months of the different seasons of the year have been represented. In addition, the annual mean GHI is depicted in a 3D image in Figure 14; notice that a view from the south has been chosen so that the facades most exposed to solar radiation can be appreciated.
Optionally, additional VTK files representing shadows cast in the domain for each time step can be obtained. Figure 15 represents the computed shadows cast in the study area on 4 September 2022 at 7.00 a.m., when the shadow of the cathedral tower, located in the blue dotted square, affects the reading of the pyranometers, located in the red dot.

3.4. Comparison of Simulation Results with Experimental Data

Taking into account the availability of data, due to cloud cover, possible technical difficulties and the reduced number of days on which the effect of the cathedral tower on the location of the pyranometers can be seen, the following dates have been chosen to compare the measured and calculated G H I : 4 and 11 September 2022, when the measurements will show the effect of the shadow of the cathedral, and 15 March 2021 and 4 August 2022, when the shadow phenomenon will not occur.
Both measured and calculated G H I values have been taken every 5 min throughout the 24 h of the day, with a total of 288 data points. In the two lower graphics of Figure 16, it can be clearly appreciated that the model captures the effect of the shadow of the cathedral in the early morning. The cathedral tower shadows the pyranometers just at sunrise, causing an interruption in the increase in irradiance, either with a sharp and short decrease (Figure 16, bottom left) or with a sharp and delayed sunrise (Figure 16, bottom right). Both phenomenologies are well captured in the simulation, with the pyranometer measurement (purple line) matching the calculated G H I (green line) very accurately. The upper graphics correspond to two dates where the shadow of the cathedral does not affect the irradiance data reading. The upper left graph shows the effect of cloud cover around 9 a.m. on pyranometer measurements (purple line).
To compare the accuracy of the results and the computational cost, the selected days were also simulated using the high resolution (1 m) mesh. The results at the pyranometer’s locations are the same (with 6-digit accuracy) as those obtained using the 5 m mesh, but the mean execution time is 68 times higher, so the 5 m mesh was chosen for all irradiance simulations. However, the use of the high resolution mesh might be appropriate in case a very precise study of shadows over the whole domain is intended.
In order to evaluate the precision of the model for a clear sky, four representative statistical error indicators have been calculated for the selected dates.
  • M A E : Mean Absolute Error
    M A E = i = 1 n | G H I ^ i G H I i | n
  • N M A E : Normalized Mean Absolute Error
    N M A E = i = 1 n | G H I ^ i G H I i | n × G H I m a x · 100
  • R M S E : Root-Mean-Square Error
    R M S E = i = 1 n ( G H I ^ i G H I i ) 2 n
  • N R M S E : Normalized Root-Mean-Square Error
    N R M S E = i = 1 n ( G H I ^ i G H I i ) 2 n × G H I m a x 2
  • R 2 : Coefficient of determination
    R 2 = 1 i = 1 n ( G H I ^ i G H I i ) 2 i = 1 n ( G H I ^ i G H I ¯ ) 2
Here, G H I i and G H I ^ i represent, respectively, the measured and calculated value of the global irradiance at time t i , corresponding to the n = 288 time instants for which data are available. G H I m a x refers to the maximum measured G H I , and G H I ¯ is the mean of the measured GHI. Table 5 summarizes the statistical indicators described above for the four selected days, showing very small values, especially the normalized errors, resulting in a higher coefficient of determination ( R 2 0.99 ).

4. Discussion and Conclusions

The combination of different technologies from diverse areas of knowledge has allowed the validation of a solar radiation model that is particularly useful in complex orographies such as those with buildings in an urban area. This is especially useful for the optimized design of solar installations that are increasingly being installed in cities.
The solar radiation model MAPSol allows us to estimate the clear-sky beam irradiance and generate a very detailed shadow map, which is computed at a reasonable computational cost by pre-filtering the triangles of the surface mesh and using an adapted mesh.
For areas with complicated orography, a high-resolution DEM is required, which is calculated from LiDAR data with a resolution of 33 cm. Adapting the mesh associated with this DEM allows the number of elements to be reduced by more than 90 % without reducing the accuracy of the calculations, but at a much lower computational cost. The comparison of results and computational cost using two adapted meshes of different resolutions, one finer (1 m tolerance) and the other coarser (5 m tolerance), allows us to conclude that using the coarse mesh considerably reduces the computational cost without losing accuracy.
The accuracy of the MAPSol model has been validated by comparing the calculated G H I data with data measured by two pyranometers. The calculated G H I accuracy is close to 99 % on days with low cloud cover. It is worth noting that on some of these days, the sensors detect the shadow of the dome of the cathedral tower in the early hours of the day.
Therefore, the proposed procedure and model allow the calculation of irradiance in areas of very complex orography with a reduced computational cost and freely available data, achieving a great balance between efficiency and effectiveness.
In addition, the adapted mesh generated from the high-resolution DEM has proven to be useful for recreating 3D images of complex areas with limited resources and freely available data.
In the future, the described procedure and model will be used to elaborate a solar map of a larger area of the city of Salamanca. In addition, the comparison of the measured irradiance values, subject to variable cloud cover, and those calculated for clear skies, can be used to produce cloud cover indices, and progress in the use of the MAPSol model under real sky conditions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/s24061823/s1.

Author Contributions

E.R., A.O.-S. and G.M.G.: formal analysis, investigation, software, validation, writing—original draft and visualization concerning the solar radiation model; J.G.-F., M.J.S. and A.M.: formal analysis, investigation, software, validation, writing—original draft and visualization concerning the data acquisition; M.S.-A., P.A.-A. and S.L.: formal analysis, investigation, software, validation, writing—original draft and visualization concerning the high-resolution model; J.M.C.: formal analysis, investigation, software, validation, writing—original draft and visualization concerning the adapted mesh; J.M.I.: investigation and writing—original draft concerning the state of the art; M.I.A.: conceptualization, funding acquisition, project administration, supervision, writing—review and editing; R.M.A.: supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by the Department of Education of the regional government, the Junta de Castilla y León, grant contract SA089P20; by the European Union’s Horizon 2020-Research and Innovation Framework Program under grant agreement ID 101036926; by the Ministry of Education, Culture and Sports (Government of Spain) through human resources under grant number FPU19/06034; by the Ministry of Science and Innovation (Government of Spain) through project funding PID2021-127089OB-I00; by the European Regional Development Fund (ERDF) grant contract EIS2021-04 (IUSIANI-ODS); J.G.-F. was funded by NextGEN EU-Junta de Castilla y León, with a research contract within “Research Program”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available datasets were analyzed in this study. These data can be found here: https://centrodedescargas.cnig.es/CentroDescargas (accessed on 31 January 2024).

Acknowledgments

A. González-Sánchez and P. López Rico are acknowledged for pyranometers data visualization and the composition of some figures.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Abbreviations

The following abbreviations are used in this manuscript:
WPWarning points
DTMDigital Terrain Model
DEMDigital Elevation Model
DSMDigital Slope Model
GHIGlobal Horizontal Irradiance
DHIDiffuse Horizontal Irradiance
DNIDirect Normal Irradiance
BHIBeam Horizontal Irradiance
CSPConcentrating Solar Power
GISGeographical Information System
IGNNational Geographic Institute

References

  1. Wong, L.; Chow, W. Solar radiation model. Appl. Energy 2001, 69, 191–224. [Google Scholar] [CrossRef]
  2. Olson, M.; Rupper, S. Impacts of topographic shading on direct solar radiation for valley glaciers in complex topography. Cryosphere 2019, 13, 29–40. [Google Scholar] [CrossRef]
  3. Rosskopf, E.; Morhart, C.; Nahm, M. Modelling Shadow Using 3D Tree Models in High Spatial and Temporal Resolution. Remote Sens. 2017, 9, 719. [Google Scholar] [CrossRef]
  4. Kaynak, S.; Kaynak, B.; Özmen, A. A software tool development study for solar energy potential analysis. Energy Build. 2018, 162, 134–143. [Google Scholar] [CrossRef]
  5. Choi, Y.; Suh, J.; Kim, S.M. GIS-Based Solar Radiation Mapping, Site Evaluation, and Potential Assessment: A Review. Appl. Sci. 2019, 9, 1960. [Google Scholar] [CrossRef]
  6. Merchán, R.P.; Santos, M.J.; Medina, A.; Calvo Hernández, A. High temperature central tower plants for concentrated solar power: 2021 overview. Renew. Sust. Ener. Rev. 2022, 155, 111828. [Google Scholar] [CrossRef]
  7. Vulkan, A.; Kloog, I.; Dorman, M.; Erell, E. Modeling the potential for PV installation in residential buildings in dense urban areas. Energy Build. 2018, 169, 97–109. [Google Scholar] [CrossRef]
  8. Dorman, M.; Erell, E.; Vulkan, A.; Kloog, I. shadow: R Package for Geometric Shadow Calculations in an Urban Environment. R J. 2019, 11, 287. [Google Scholar] [CrossRef]
  9. Li, Y.; Liu, C. Estimating solar energy potentials on pitched roofs. Energy Build. 2017, 139, 101–107. [Google Scholar] [CrossRef]
  10. Toledo, C.; Gracia Amillo, A.M.; Bardizza, G.; Abad, J.; Urbina, A. Evaluation of Solar Radiation Transposition Models for Passive Energy Management and Building Integrated Photovoltaics. Energies 2020, 13, 702. [Google Scholar] [CrossRef]
  11. Li, S.Y.; Han, J.Y. The impact of shadow covering on the rooftop solar photovoltaic system for evaluating self-sufficiency rate in the concept of nearly zero energy building. Sustain. Cities Soc. 2022, 80, 103821. [Google Scholar] [CrossRef]
  12. Stendardo, N.; Desthieux, G.; Abdennadher, N.; Gallinelli, P. GPU-Enabled Shadow Casting for Solar Potential Estimation in Large Urban Areas. Application to the Solar Cadaster of Greater Geneva. Appl. Sci. 2020, 10, 5361. [Google Scholar] [CrossRef]
  13. Liang, J.; Gong, J.; Xie, X.; Sun, J. Solar3D: An Open-Source Tool for Estimating Solar Radiation in Urban Environments. ISPRS Int. J. Geo-Inf. 2020, 9, 524. [Google Scholar] [CrossRef]
  14. Brito, M.C.; Redweik, P.; Catita, C.; Freitas, S.; Santos, M. 3D Solar Potential in the Urban Environment: A Case Study in Lisbon. Energies 2019, 12, 3457. [Google Scholar] [CrossRef]
  15. Tripathy, M.; Yadav, S.; Sadhu, P.; Panda, S. Determination of optimum tilt angle and accurate insolation of BIPV panel influenced by adverse effect of shadow. Renew. Energy 2017, 104, 211–223. [Google Scholar] [CrossRef]
  16. Zhu, R.; Wong, M.S.; You, L.; Santi, P.; Nichol, J.; Ho, H.C.; Lu, L.; Ratti, C. The effect of urban morphology on the solar capacity of three-dimensional cities. Renew. Energy 2020, 153, 1111–1126. [Google Scholar] [CrossRef]
  17. Perez, R.; Seals, R.; Ineichen, P.; Stewart, R.; Menicucci, D. A new simplified version of the perez diffuse irradiance model for tilted surfaces. Sol. Energy 1987, 39, 221–231. [Google Scholar] [CrossRef]
  18. Biljecki, F.; Ledoux, H.; Stoter, J. Does a Finer Level of Detail of a 3D City Model Bring an Improvement for Estimating Shadows? In Advances in 3D Geoinformation; Lecture Notes in Geoinformation and Cartography; Springer International Publishing: Cham, Switzerland, 2016; pp. 31–47. [Google Scholar] [CrossRef]
  19. Peronato, G.; Rey, E.; Andersen, M. 3D model discretization in assessing urban solar potential: The effect of grid spacing on predicted solar irradiation. Sol. Energy 2018, 176, 334–349. [Google Scholar] [CrossRef]
  20. Zhou, Q.Y.; Neumann, U. 2.5D Dual Contouring: A Robust Approach to Creating Building Models from Aerial LiDAR Point Clouds. In Proceedings of the 11th European Conference on Computer Vision, Heraklion, Greece, 5–11 September 2010; Lecture Notes in Computer Science. Springer: Berlin/Heidelberg, Germany, 2010; pp. 115–128. [Google Scholar] [CrossRef]
  21. North, G.; Zhang, F.; Pyle, J. Encyclopedia of Atmospheric Sciences, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2014. [Google Scholar]
  22. Šúri, M.; Hofierka, J. A new GIS-based solar radiation model and its application to photovoltaic assessments. Trans. GIS 2004, 8, 175–190. [Google Scholar] [CrossRef]
  23. National Renewable Energy Laboratory. Available online: https://www.nrel.gov/grid/solar-resource/solar-glossary.html (accessed on 31 January 2024).
  24. Díaz, F.; Montero, G.; Escobar, J.; Rodríguez, E.; Montenegro, R. An adaptive solar radiation numerical model. J. Comput. Appl. Math. 2012, 236, 4611–4622. [Google Scholar] [CrossRef]
  25. Díaz, F.; Montero, G.; Escobar, J.; Rodríguez, E.; Montenegro, R. A new predictive solar radiation numerical model. Appl. Math. Comput. 2015, 267, 596–603. [Google Scholar] [CrossRef]
  26. Montero, G.; Escobar, J.; Rodríguez, E.; Montenegro, R. Solar radiation and shadow modelling with adaptive triangular meshes. Sol. Energy 2009, 83, 998–1012. [Google Scholar] [CrossRef]
  27. Page, J.K. (Ed.) Prediction of Solar Radiation on Inclined Surfaces; D. Reidel Publishing Co.: Dordrecht, The Netherlands, 1986. [Google Scholar]
  28. Kasten, F.; Young, A.T. Revised optical air mass tables and approximation formula. Appl. Opt. 1989, 28, 4735. [Google Scholar] [CrossRef]
  29. Kasten, F. The Linke turbidity factor based on improved values of the integral Rayleigh optical thickness. Sol. Energy 1996, 56, 239–244. [Google Scholar] [CrossRef]
  30. Blanco-Muriel, M.; Alarcón-Padilla, D.C.; López-Moratalla, T.; Lara-Coira, M. Computing the solar vector. Sol. Energy 2001, 70, 431–441. [Google Scholar] [CrossRef]
  31. Díaz, F.; Montero, H.; Santana, D.; Montero, G.; Rodríguez, E.; Mazorra Aguiar, L.; Oliver, A. Improving shadows detection for solar radiation numerical models. Appl. Math. Comput. 2018, 319, 71–85. [Google Scholar] [CrossRef]
  32. Rivara, M. A grid generator based on 4-triangles conforming mesh-refinement algorithms. Int. J. Numer. Methods Eng. 1987, 24, 1343–1354. [Google Scholar] [CrossRef]
  33. Toth, C.; Jóźków, G. Remote sensing platforms and sensors: A survey. ISPRS J. Photogramm. Remote Sens. 2016, 115, 22–36. [Google Scholar] [CrossRef]
  34. Deibe, D.; Amor, M.; Doallo, R. Big Data Geospatial Processing for Massive Aerial LiDAR Datasets. Remote Sens. 2020, 12, 719. [Google Scholar] [CrossRef]
  35. Gassar, A.A.A.; Cha, S.H. Review of geographic information systems-based rooftop solar photovoltaic potential estimation approaches at urban scales. Appl. Energy 2021, 291, 116817. [Google Scholar] [CrossRef]
  36. Bartha, G.; Kocsis, S. Standardization of Geographic Data: The European Inspire Directive. Eur. J. Geogr. 2022, 2. Available online: https://eurogeojournal.eu/index.php/egj/article/view/36/10 (accessed on 31 January 2024).
  37. Directive 2007/2/EC of the European Parliament and of the Council of 14 March 2007 Establishing an Infrastructure for Spatial Information in the European Community (INSPIRE). 2023. Available online: https://eur-lex.europa.eu/legal-content/EN/ALL/?uri=CELEX%3A32007L0002 (accessed on 31 January 2024).
  38. Plan Nacional de Observacion del Territorio—PNOT. 2023. Available online: https://www.ign.es/web/plan-nacional-de-observacion-del-territorio (accessed on 31 January 2024).
  39. Sánchez-Aparicio, M.; González-González, E.; Martín-Jiménez, J.A.; Lagüela, S. Solar Potential Analysis of Bus Shelters in Urban Environments: A Study Case in Ávila (Spain). Remote Sens. 2023, 15, 5189. [Google Scholar] [CrossRef]
  40. Sánchez-Aparicio, M.; Martín-Jiménez, J.A.; González-González, E.; Lagüela, S. Laser Scanning for Terrain Analysis and Route Design for Electrified Public Transport in Urban Areas. Remote Sens. 2023, 15, 3325. [Google Scholar] [CrossRef]
  41. Sánchez-Aparicio, M.; Del Pozo, S.; Martín-Jiménez, J.A.; González-González, E.; Andrés-Anaya, P.; Lagüela, S. Influence of LiDAR Point Cloud Density in the Geometric Characterization of Rooftops for Solar Photovoltaic Studies in Cities. Remote Sens. 2020, 12, 3726. [Google Scholar] [CrossRef]
  42. Martín-Jiménez, J.; Del Pozo, S.; Sánchez-Aparicio, M.; Lagüela, S. Multi-scale roof characterization from LiDAR data and aerial orthoimagery: Automatic computation of building photovoltaic capacity. Autom. Constr. 2020, 109, 102965. [Google Scholar] [CrossRef]
  43. Organismo Autónomo Centro Nacional de Información Geográfica, Centro de Descargas. 2023. Available online: https://centrodedescargas.cnig.es/CentroDescargas/index.jsp (accessed on 31 January 2024).
  44. Cascón, J.; Ferragut, L.; Asensio, M.; Prieto, D.; Álvarez, D. Neptuno ++: An adaptive finite element toolbox for numerical simulation of environmental problems. In Proceedings of the XVIII Spanish-French School Jacques- Louis Lions about Numerical Simulation in Physics and Engineering, Las Palmas de Gran Canaria, Spain, 25–29 June 2018. [Google Scholar]
  45. Escobar, J.; Rodríguez, E.; Montenegro, R.; Montero, G.; González-Yuste, J. Simultaneous untangling and smoothing of tetrahedral meshes. Comput. Methods Appl. Mech. Eng. 2003, 192, 2775–2787. [Google Scholar] [CrossRef]
  46. Montenegro, R.; Cascon, J.M.; Escobar, J.M.; Rodriguez, E.; Montero, G. An automatic strategy for adaptive tetrahedral mesh generation. Appl. Numer. Math. 2009, 59, 2203–2217. [Google Scholar] [CrossRef]
  47. de la Casinière, A.C.; Cachorro-Revilla, V.E. La Radiación Solar en el Sistema Tierra-Atmósfera; Secretariado de Publicaciones e Intercambio Universidad de Valladolid, Secretariado de Publicaciones: Valladolid, Spain, 2008. [Google Scholar]
  48. Kipp & Zonen. Instruction Manual. CM 121 Shadow Ring; OTT HydroMet: Delft, The Netherlands, 2004. [Google Scholar]
  49. Bardia, R. Dispositivos Electrónicos. Fundamentos de Electrónica, Volumen 5; Universidad Politécnica de Catalunya: Barcelona, Spain, 1999. [Google Scholar]
  50. Kipp & Zonen. Instruction Manual. SMP Series. Smart Pyranometer; OTT HydroMet: Delft, The Netherlands, 2017. [Google Scholar]
  51. Ávila, J.M.S.; Martín, J.R.; Alonso, C.J.; de Cos Escuin, M.C.S.; Cadalso, J.M.; Bartolomé, M.L. Atlas de Radiación Solar en España Utilizando Datos de SAF de Clima de EUMETSAT; Agencia Estatal de Meteorología: Madrid, Spain, 2012. [Google Scholar]
  52. Solargis. Solar Data behind the Maps. 2023. Available online: https://solargis.com/es/maps-and-gis-data/techspecs (accessed on 31 January 2024).
  53. University of Oregon. Sun Path Chart Program. 2023. Available online: http://solardata.uoregon.edu/SunChartProgram.php (accessed on 31 January 2024).
Figure 1. Reference system and Euler angles.
Figure 1. Reference system and Euler angles.
Sensors 24 01823 g001
Figure 2. Warning points for shading are equidistributed over each triangle (left graph) simply using the centers of the triangles obtained by refining with the 4-T Rivara algorithm (right graph) [32].
Figure 2. Warning points for shading are equidistributed over each triangle (left graph) simply using the centers of the triangles obtained by refining with the 4-T Rivara algorithm (right graph) [32].
Sensors 24 01823 g002
Figure 3. Scheme of the mesh generation procedure.
Figure 3. Scheme of the mesh generation procedure.
Sensors 24 01823 g003
Figure 4. Schematic of shadow-ring (positioning, taken from Kipp and Zonen. Instructions Manual CM121 Shadow Ring [48]).
Figure 4. Schematic of shadow-ring (positioning, taken from Kipp and Zonen. Instructions Manual CM121 Shadow Ring [48]).
Sensors 24 01823 g004
Figure 6. Spectral range response of a pyranometer of the SMP series manufactured by Kipp and Zonen versus the radiation spectrum at sea level (taken from Kipp and Zonen. Instruction Manual SMP Series [50]).
Figure 6. Spectral range response of a pyranometer of the SMP series manufactured by Kipp and Zonen versus the radiation spectrum at sea level (taken from Kipp and Zonen. Instruction Manual SMP Series [50]).
Sensors 24 01823 g006
Figure 7. Measurement devices installed on the rooftop of the Trilingüe building of the Faculty of Sciences (40.96062 N, 5.67075 W) of the University of Salamanca (Spain): (a) pyranometer Kipp and Zonen SMP10 and (b) pyranometer SMP10 with shadow ring model Kipp and Zonen CM 121.
Figure 7. Measurement devices installed on the rooftop of the Trilingüe building of the Faculty of Sciences (40.96062 N, 5.67075 W) of the University of Salamanca (Spain): (a) pyranometer Kipp and Zonen SMP10 and (b) pyranometer SMP10 with shadow ring model Kipp and Zonen CM 121.
Sensors 24 01823 g007
Figure 8. Averaged curves of G H I and B H I for different representative months of a year.
Figure 8. Averaged curves of G H I and B H I for different representative months of a year.
Sensors 24 01823 g008
Figure 9. Sun path chart: Apparent position of the Sun from the Trilingüe building at the University of Salamanca (40.96062° N, 5.670759° W) between 21 December and 21 June [53]. A panoramic photo, taken from pyranometer that register G H I , is overlapped aiming to identify shadowing sources.
Figure 9. Sun path chart: Apparent position of the Sun from the Trilingüe building at the University of Salamanca (40.96062° N, 5.670759° W) between 21 December and 21 June [53]. A panoramic photo, taken from pyranometer that register G H I , is overlapped aiming to identify shadowing sources.
Sensors 24 01823 g009
Figure 10. Three-dimensional view of the study area: (a) original point cloud and (b) DEM derived from the point cloud.
Figure 10. Three-dimensional view of the study area: (a) original point cloud and (b) DEM derived from the point cloud.
Sensors 24 01823 g010
Figure 11. Adapted thin mesh (1 m) of the complete area (c), zoom of the fine adapted mesh over the Cathedral area (b), and detail of the uniform original mesh over the Cathedral area (a).
Figure 11. Adapted thin mesh (1 m) of the complete area (c), zoom of the fine adapted mesh over the Cathedral area (b), and detail of the uniform original mesh over the Cathedral area (a).
Sensors 24 01823 g011
Figure 12. Three-dimensional reconstruction of the study area by simply projecting the orthophoto onto the fine adapted mesh. The pyranometers described in Section 2.5 are located at the red dot.
Figure 12. Three-dimensional reconstruction of the study area by simply projecting the orthophoto onto the fine adapted mesh. The pyranometers described in Section 2.5 are located at the red dot.
Sensors 24 01823 g012
Figure 13. Mean GHI map for January, April, July and October, computed with MAPSol and coarse adapted mesh. Notice that the calculations assume a clear sky.
Figure 13. Mean GHI map for January, April, July and October, computed with MAPSol and coarse adapted mesh. Notice that the calculations assume a clear sky.
Sensors 24 01823 g013
Figure 14. Annual mean GHI 3D map, computed with MAPSol and coarse adapted mesh. View from the south.
Figure 14. Annual mean GHI 3D map, computed with MAPSol and coarse adapted mesh. View from the south.
Sensors 24 01823 g014
Figure 15. Shadow calculated with MAPSol on 4 September 2022 at 7.00 a.m. and coarse adapted mesh (Supplementary Materials). The long shadow of the cathedral tower (blue dotted square) can be seen over the area where the pyranometers are located (red dot).
Figure 15. Shadow calculated with MAPSol on 4 September 2022 at 7.00 a.m. and coarse adapted mesh (Supplementary Materials). The long shadow of the cathedral tower (blue dotted square) can be seen over the area where the pyranometers are located (red dot).
Sensors 24 01823 g015
Figure 16. Curves of G H I measured (purple lines) and simulated (green lines) for the selected dates. The effect of the shadow of the cathedral tower at sunrise can be appreciated in the simulated and the measurement data (circled in red in bottom graphics). The accuracy of the irradiance fit calculated with the MAPSol model is very good, as can be seen in all of the graphs. The peak at 9 a.m. in the upper left graphic corresponds to a cloudy interval that affected the pyranometer readings, which cannot be simulated as the model assumes a clear sky.
Figure 16. Curves of G H I measured (purple lines) and simulated (green lines) for the selected dates. The effect of the shadow of the cathedral tower at sunrise can be appreciated in the simulated and the measurement data (circled in red in bottom graphics). The accuracy of the irradiance fit calculated with the MAPSol model is very good, as can be seen in all of the graphs. The peak at 9 a.m. in the upper left graphic corresponds to a cloudy interval that affected the pyranometer readings, which cannot be simulated as the model assumes a clear sky.
Sensors 24 01823 g016
Table 1. Main specifications of the LiDAR data of the first coverage offered by the IGN.
Table 1. Main specifications of the LiDAR data of the first coverage offered by the IGN.
FeatureFirst Coverage
Minimum point density 0.5 pt / m 2
Years of flight2009–2015
Geodetic reference systemETRS89 zones 28, 29, 30 and 31 as appropriate
Altimetric reference systemOrthometric altitudes, reference geoid EGM08
RMSE Z≤40 cm
Estimated planimetric accuracy≤30 cm
File size 2 × 2 km
File formatLAS 1.2 format 3
Table 2. Main features of the Kipp and Zonen SMP 10 pyranometers [50] used to perform the experimental measurements.
Table 2. Main features of the Kipp and Zonen SMP 10 pyranometers [50] used to perform the experimental measurements.
FeatureValue
Spectral range285–2800 nm
Response time(63%) < 0.7 s
Response time(95%) < 2 s
Non-linearity<0.2
Spectral selectivity(350–1500 nm) < 3%
Field of view 180 °
Table 3. Experimental and bibliographic values of accumulated energy by irradiance type in kWh m 2 day 1 for each month of the year.
Table 3. Experimental and bibliographic values of accumulated energy by irradiance type in kWh m 2 day 1 for each month of the year.
SourceAEMET [51]Experimental DataRelative Differences (%)
Month GHI BHI GHI BHI Δ GHI Δ BHI
January2.081.182.311.4711.0624.58
February3.091.893.091.970.00 4.23
March4.492.824.743.085.57 9.22
April5.563.505.192.896.6517.43
May6.444.086.904.657.1413.97
June7.605.457.335.133.55 5.87
July7.825.967.826.170.00 3.52
August6.845.056.955.481.61 8.51
September5.273.715.213.751.14 1.08
October3.432.143.532.322.92 8.41
November3.381.282.261.2733.14 0.78
December1.780.961.530.6714.0430.21
Table 4. Annual cumulative values (in kWh · m−2·yr−1) and average value per day (in kWh · m−2·day−1) of the energy received.
Table 4. Annual cumulative values (in kWh · m−2·yr−1) and average value per day (in kWh · m−2·day−1) of the energy received.
SourceAEMET [51]Solargis [52]Measured Records
AnnualMax.Min.Max.Min.
G H I 1708.2 1733.8 168017531733.65
B H I 1146.0 1182.6 1185.93
DailyMax.Min.Max.Min.
G H I 4.68 4.75 4.6 4.8 4.75
B H I 3.14 3.24 3.25
Table 5. Summary of errors, in terms of M A E and M R S E and the corresponding normalized indicators N M A E and N M R S E , as well as the coefficient of determination R 2 .
Table 5. Summary of errors, in terms of M A E and M R S E and the corresponding normalized indicators N M A E and N M R S E , as well as the coefficient of determination R 2 .
Date MAE NMAE MRSE NMRSE R 2
15 March 2021 9.8632 1.2207 19.5744 0.0242 0.9959
4 August 2022 7.1667 0.7465 9.1455 0.0095 0.9994
4 September 2022 14.1427 1.5802 18.7152 0.0209 0.9971
11 September 2022 7.0888 0.8224 12.0603 0.0140 0.9986
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rodríguez, E.; García-Ferrero, J.; Sánchez-Aparicio, M.; Iglesias, J.M.; Oliver-Serra, A.; Santos, M.J.; Andrés-Anaya, P.; Cascón, J.M.; Montero García, G.; Medina, A.; et al. Validation of a 3D Local-Scale Adaptive Solar Radiation Model by Using Pyranometer Measurements and a High-Resolution Digital Elevation Model. Sensors 2024, 24, 1823. https://doi.org/10.3390/s24061823

AMA Style

Rodríguez E, García-Ferrero J, Sánchez-Aparicio M, Iglesias JM, Oliver-Serra A, Santos MJ, Andrés-Anaya P, Cascón JM, Montero García G, Medina A, et al. Validation of a 3D Local-Scale Adaptive Solar Radiation Model by Using Pyranometer Measurements and a High-Resolution Digital Elevation Model. Sensors. 2024; 24(6):1823. https://doi.org/10.3390/s24061823

Chicago/Turabian Style

Rodríguez, Eduardo, Judit García-Ferrero, María Sánchez-Aparicio, José M. Iglesias, Albert Oliver-Serra, M. Jesús Santos, Paula Andrés-Anaya, J. Manuel Cascón, Gustavo Montero García, Alejandro Medina, and et al. 2024. "Validation of a 3D Local-Scale Adaptive Solar Radiation Model by Using Pyranometer Measurements and a High-Resolution Digital Elevation Model" Sensors 24, no. 6: 1823. https://doi.org/10.3390/s24061823

APA Style

Rodríguez, E., García-Ferrero, J., Sánchez-Aparicio, M., Iglesias, J. M., Oliver-Serra, A., Santos, M. J., Andrés-Anaya, P., Cascón, J. M., Montero García, G., Medina, A., Lagüela, S., Asensio, M. I., & Montenegro Armas, R. (2024). Validation of a 3D Local-Scale Adaptive Solar Radiation Model by Using Pyranometer Measurements and a High-Resolution Digital Elevation Model. Sensors, 24(6), 1823. https://doi.org/10.3390/s24061823

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

Article Metrics

Back to TopTop