In this work an approach to the theoretical solar radiation in the South Region of Ecuador is presented, using the Hottel Model ([1] [2]), considering the cloudiness and the topography Shadows. The maps o f Global Solar Radiation are shown for the months of January, August and Annual one, considering the heights over the middle sea level extracted of the Digital Elevation Model. This work will allow specialists and engineers to know the most favorable areas for the location of photovoltaic or photothermic plants in the southern region of Ecuador. The effects have been analyzed by the orography and the methods used were discussed. The role of the orography is analyzed to improve the values’ product of the modeling of the solar radiation. Conclusions and recommendations are given.
The present work is part of the hypothesis of that it is possible to obtain maps of radiation from theoretical models that are consistent with the reality, and the objective is the preparation of maps of global radiation that allow their use to know with greater precision solar energy potentials in the Southern Region of Ecuador (RSE), as for the determination of different environmental components, such as potential evapotranspiration [
The work focused on obtaining the values of direct, diffuse and global radiation on a horizontal surface, considering a clear sky, using the “Hottel model” [
It is known that the atmosphere exerts a redistribution effect of the radiation that receives from the sun, for example, on a very clear day; a relatively small part becomes diffuse radiation, whereas most of it remains as direct. On the other hand, on a cloudy day, the redistribution of radiation is much more noticeable. Dense clouds have a very high albedo (reflected energy fraction), which makes a large part of the solar radiation reflected in the outer space on a densely cloudy day. In addition, the energy that manages to pass through the clouds is only diffuse radiation.
According to recent developments in the calibration of radiometric broadband, solar instrumentation and improvements in the accuracy of the measurements are profiled in relation to the diffuse radiation and the calibration of radiometers [
Errors between 25 and 100 w/m2 are possible in the pirheliometric data and errors of 25 w/m2 in the Pirihelio Métricos data can occur under clear sky conditions. These instrumental errors are important if it is modeling instant profiles, schedules, subtimes, or diaries. Analyses of annual totals of corrected and uncorrected data have shown differences of less than 0.1% in the outcome [
On the other hand, a reflectance map can be based on experimental data. It can mount the sample of a surface and measure its apparent brightness from a fixed display point under fixed lighting conditions, while its orientation varies. Similarly, it can take a box from a test object in a known way and calculate the orientation of the corresponding surface element for each point in the image. The reflectance map is then obtained by reading the brightness measurements at each point [
Alternatively, it can even use more detailed information about light reflection from the surface. The bidirectional Reflectance distribution function (BRDF) describes how bright a view surface can appear from a specific address when it is illuminated from another direction fixed with Advance. Through integration on a given distribution of light sources, the reflectance map can be calculated from this information.
Then, a reflectance map is similar to a “convolution” of the BRDF and the distribution of the radiation source. More commonly, reflectance maps are based on phenomenological models rather than physical reality. The so-called “Lambertianas surfaces”, or the perfect diffuser, for example, have the property of appearing equally bright from all directions from where it is observed.
In the same way they reflect all the light without absorbing anything. In fact, these two conditions are sufficient to determine in a unique way the BRDF of such a surface, and from this the reflectance map, providing the positions of the given light sources. Some reflectance maps are based on mathematical models of the interaction of light with the surface. These models tend to be very complex to obtain an analytical solution or very simple to represent effectively the real surfaces. However, some have come very close to the prediction of observed behavior of particular surfaces.
New reflectance maps were determined based on methods to produce maps of shaded reliefs that do not depend on the comprehension of the process of formation of the images or the physics of the reflection of light, but of the way in which the light at one point depends on the gradient of the underlying geographical surface.
Because reflectance maps provide apparent brilliance as a function only of a local surface gradient, it does not take into account the effects dependent on the position of the surface element. One of these effects is the illumination of one surface element on the other. Fortunately, this mutual lighting effect is small unless the surface reflectance is very high.
Another position-dependent effect on the apparent brightness is blocking the light by a portion of the surface before the light reaches another.
Likewise the scattering of light by the molecules of air and aerosols changes the apparent brilliance of an element of surface seen through the atmosphere. However, the difference between brilliance and the value of the environment decreases with the thickness of the gaseous layer through which the surface is visualized [
For the realization of this work, initially, the preparation of an Excel book was started using the methodology proposed by Passamai (2000) [
Later and due to the need to know the behavior of solar radiation, and to determine also, which sectors of the southern region of Ecuador and at what moment of the year, have the greatest solar potential, and considering the absence of sufficient series of observations of this meteorological variable, were used the theoretical methods of Hottel and PVGIS for the calculation of the variable, which are presented in the work of Molina and others [
The model of Hottel [
τ b = a 0 + a 1 ⋅ e ( − k / cos ( θ z ) ) (1)
Here a0, a1 and k are parameters adjusted empirically, used by [
θ z = a cos ( sin ( r a d i a n e s ( λ ) ) ) × sin ( r a d i a n e s ( δ ) ) + cos ( r a d i a n e s ( λ ) ) × cos ( r a d i a n e s ( δ ) ) × cos (ω)
where:
λ = Latitude;
δ = Declination of the Sun;
ω = (−0.261799 × (HS-12)) × 360/2π;
HS = Solar Time.
For estimating irradiation or diffuse irradiance τd for a horizontal surface is used the expression of Liu and Jordan [
τ d = 0.2710 − 0.2939 τ b (2)
where τd and τb are adimensional.
By integrating these equations, from the sunrise to the sunset, it can get irradiation over a day. In the Hottel model, the standard atmosphere represents a typical average latitude and pollution-free condition, independent of the thickness of the ozone.
Elevation data obtained from the digitized topographic map of the entire southern region of Ecuador were used, with corresponding level curves from the information pertaining to the Radar Topography Mision (SRTM) [
The geographical points were selected from the values of longitude and geographic latitude and the height on the middle sea level obtained from the numerical model of altitude for the project zone, calculated the time values of direct, diffuse and global radiation for each point and subsequently, the corresponding global radiation maps were made, for each month and for the yearly.
Once the corresponding maps were obtained, we proceeded to rectify those using data for reanalysis of the daytime clouding values, recalculating the maps. From the monthly and annual values of global radiation, on a horizontal surface, considering the condition of clear day, for the different months, as well as the annual average, obtained by applying corrections due to the cloudiness, obtained
from the reanalysis carried out by NASA, the values of the global radiation are obtained for the months and for the annual average [
Reanalysis data were obtained from the NASA website: Surface Meteorology and Solar Energy [
From the data of the tables obtained for each one of the points of a grid with the ends: −77.0˚W, −81.0˚W; and −3.0˚S, −6.0˚S, longitude and latitude respectively, with a cell size of 0.001˚, the average monthly and annual values were calculated for the period of 22 years included in the reanalysis, by means of an interpolation using a linear Krigging model.
We worked on the Software Envi 4.7, where we used the topographical modeling on data of the numerical model of altitude (MNA) obtained from the data SRTM and the Geological Service of the United States (USGS) [
When the color of a pixel (digital unit of an image) based on their neighbors, it is talking about a window or Kernel (usually a 3 × 3 or 5 × 5 grid); Which is nothing more than a convolution matrix that indicates which pixels should influence the original pixel and how much. A weight is made for each pixel covered by the specified Kernel. This is the basis for many filtering algorithms [
On the other hand, the possibilities offered by the ENVI geographic Information system for topographical modeling are used to obtain brilliant or shaded embossing maps. The following
In order to compare the results, data from NASA’s reanalysis for insolation were obtained, which are shown in
units, from which the corresponding maps were made. An overlay was made for the data of the monthly sum of the theoretical radiation obtained by the Hottel method for each particular month and for the annual average, by means of multiplication with the percent of sun hours obtained for each month and annual average to determine the quantity of radiation corresponding to each of the periods.
The shading maps obtained from the Digital elevation model were used, considering the following combinations of elevation angles and azimuth: (15.45); (15.135); (15.225); (15.315); (30.45); (30.135); (30.225); (30.315); (45.45); (45.135); (45.225); (45.315); (60.45); (60.135); (60.225); (60.315). From these maps were made the average maps according to the azimuth North and South and finally were averaged to obtain the map of shadows for the entire project area.
Once the shadow map was obtained, the maps of global solar radiation were multiplied for months and yearly considering the hours of sunshine for each period by the shadow map for the area. A comparison was made of the maps with respect to the maps obtained from the reanalysis of NASA’s insolation [
Figures 4-6 show the maps of the sum of the solar radiation in clear days corresponding to the months of January, August and the annual average in the study area.
Figures 7-9 show the maps of the months of January, August and annual average percent of the hours without clouding obtained from NASA reanalysis [
Daily average value maps for months considering percent of clear days per month for January, August and the annual period are shown in Figures 10-12.
It is observed, that due to the presence of the cloudiness, considerably decrease the values of solar radiation possible to obtain according to the model of clear day of Hottel, in addition the zones with greater and lower values of this meteorological factor are inverted, due to the cloudiness generated in the highest areas of the Andean cordillera.
When comparing the result of annual irradiation with NASA's irradiation values considering the annual frequency of sunny days, it can be observed that there are no major differences (
Multiplications of the maps of each month and yearly were carried out by the shaded relief map (
by CONELEC [
The maximum values of global radiation, considering the days with clear skies, are to the west of the project area while the minimums appear, mainly to the east of it, coinciding with the zones of higher percent monthly of clear days and cloudy days respectively.
There are no significant differences between the annual values obtained and those reported in the NASA reanalysis [
By applying the slope and orientation of the slopes to the present model, the shaded relief map has obtained a more representative model of the average global radiation to be obtained in the southern region of Ecuador.
This scientific work has been financed by the Prometheus Project of the National Secretariat of Science, Technology and Innovation (Ecuador).
The authors declare no conflicts of interest regarding the publication of this paper.
Álvarez, O.H. and Montaño, T.E. (2019) The Global Solar Radiation in the Southern Region of Ecuador. Open Access Library Journal, 6: e5183. https://doi.org/10.4236/oalib.1105183