Next Article in Journal
An Integrated Evaluation Method of the Wind Power Ramp Event Based on Generalized Information of the Source, Grid, and Load
Previous Article in Journal
The Use of Hypergeometric Functions in Hysteresis Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Solar Resource Based on Meteorological and Geographical Data: Sonora State in Northwestern Territory of Mexico as Case Study

by
Enrique A. Enríquez-Velásquez
1,
Victor H. Benitez
2,*,
Sergey G. Obukhov
3,
Luis C. Félix-Herrán
4 and
Jorge de-J. Lozoya-Santos
5
1
School of Engineering, The University of Edinburgh, Sanderson Building, Robert Stevenson Road, The King’s Buildings, Edinburgh EH9 3FB, UK
2
Department of Industrial Engineering, Universidad de Sonora, Hermosillo 83000, Sonora, Mexico
3
Department of Electric Power and Electrical Engineering, National Research Tomsk Polytechnic University, 30, Lenin Avenue, 634050 Tomsk, Russia
4
Tecnologico de Monterrey, School of Engineering and Sciences, Blvd. Enrique Mazón López 965, Hermosillo 83000, Sonora, Mexico
5
Tecnologico de Monterrey, School of Engineering and Sciences, Ave. Eugenio Garza Sada 2501, Monterrey 64849, Nuevo Leon, Mexico
*
Author to whom correspondence should be addressed.
Energies 2020, 13(24), 6501; https://doi.org/10.3390/en13246501
Submission received: 5 November 2020 / Revised: 26 November 2020 / Accepted: 30 November 2020 / Published: 9 December 2020
(This article belongs to the Section A2: Solar Energy and Photovoltaic Systems)

Abstract

:
A model developed at the University of Tomsk, Russia, for high latitudes (over 55° N) is proposed and applied to the analysis and observation of the solar resource in the state of Sonora in the northwest of Mexico. This model utilizes satellite data and geographical coordinates as inputs. The objective of this research work is to provide a low-cost and reliable alternative to field meteorological stations and also to obtain a wide illustration of the distribution of solar power in the state to visualize opportunities for sustainable energy production and reduce its carbon footprint. The model is compared against real-time data from meteorological stations and satellite data, using statistical methods to scrutinize its accuracy at local latitudes (26–32° N), where a satisfactory performance was observed. An annual geographical view of available solar radiation against maximum and minimum temperatures for all the state municipalities is provided to identify the photovoltaic electricity generation potential. The outcomes are proof that the model is economically viable and could be employed by local governments to plan solar harvesting strategies. The results are generated from an open source model that allows calculating the available solar radiation over specific land areas, and the application potential for future planning of solar energy projects is evident.

Graphical Abstract

1. Introduction

The potential of renewable energy in Mexico is well known due to the great variety of natural resources and climates in this country. One of these resources is solar radiation, which, according to several satellite models [1] and studies from international agencies such as National Aeronautics and Space Administration (NASA) or Global Solar Atlas (SOLARGIS), is abundant in Mexico [2,3]. A high solar potential region is presented in the northwest of the country, where Sonora state is located. An advantage of Sonora for solar power development is its proximity to the US which allows technological and economic interchange. This could prompt an international renewable network collaboration. The state is mainly a desert and presents extremely high and low temperatures during the year. In addition, the promotion of solar harvesting projects for energy production would reduce greenhouse gas emissions (GHGs) in the region, where electricity production represents the majority of the state CO2 emissions, with 12.2 million metric tons forecasted for the year 2020 [4]. This sector represented 34.5% of the total GHGs in Sonora in 2005 and maintains a tendency of growing each year, according to the Commission of Ecology and Sustainable Development of the State of Sonora (CEDES) and the Center for Climate Strategies (CCS) [4]. This fact, together with the potential of the region for solar energy, promotes the expansion of solar harvesting projects for energy production as an economic and ecological way of development and municipal energy planning. A key point to the design and installation of successful solar power projects in the region is to utilize a reliable source of data to illustrate the solar resource availability in the territory.
For many years, a national network for solar measurement has been promoted by the government and several institutions such as the World Meteorological Organization (WMO), the National Autonomous University of Mexico (UNAM), the Mexican Center of Innovation in Solar Energy (CEMIE-SOL), the National Council for Science and Technology (CONACYT), the Secretary of Energy (SENER), and the Mexican National Weather Service (SMN) [5,6]. The national network consists of meteorological stations distributed throughout the national territory, which measure solar radiation, using devices like pyranometers and pyrheliometers. These weather stations are classified into two types: automatic meteorological stations (EMAs) and synoptic meteorological stations (ESIMEs). EMAs consist of a system of electrical and mechanical sensors that automatically measure, record, and transmit weather variables (wind speed and direction, atmospheric pressure, temperature, relative humidity, solar radiation, and precipitation). These measurements are taken every ten minutes, and the records are transmitted every one to three hours to the national weather system. On the other hand, ESIMEs utilize electrical sensors that automatically measure the same weather variables as EMAs, but also visibility and temperature ten centimeters above ground level. In addition, ESIMEs create and send synoptic report records to the national meteorological system every three hours. Both types of weather stations of the SMN use Coordinated Universal Time (UTC) and cover a representative area of five kilometers [7]. The main weakness of the Mexican national weather network is the limited number of meteorological stations nationwide (less than 300) [7]. This leaves several locations in the country without coverage. Furthermore, these stations require constant monitoring and maintenance in order to obtain reliable data. This requires a constant supply of qualified experts and financial resources [8,9].
Currently, according to the SMN, there are 15 meteorological stations in the state of Sonora [7], which are not enough to cover several regions of the state for solar radiation measurement. Therefore, these stations provide a general view of the solar resource in the state, but not specific data for every point in the whole territory. The Mexican network of meteorological stations of the SMN is shown in Figure 1.
Solar power plants require specific and reliable long-term data to know their viability of implementation and success of including them into the electric grid. Examples of parameters required to implement solar projects are temperature, measurements of solar radiation, geographical analysis, humidity, and meteorological data. This implies solar energy is a variable power source affected by climate, daily weather, and the day–night cycle. This illustrates the importance of specific data according to the geographical location of the intended solar power system [11]. An alternative to the use of meteorological stations is the use of mathematical models based on climatological and geographical data. These models offer a low-cost alternative compared to field stations. Moreover, mathematical models provide a good level of reliability and accuracy in results.
In recent decades, countries have researched and implemented mathematical models to predict solar radiation using different approaches. For example, a cloud observation-based model research in Gothenburg, Sweden [12], calculated global solar radiation based on the Oktas cloud scale. According to [12], this model can be used without any geographical restriction and no previous local radiation values are needed. However, this model presented limitations in resolution, which depends on long-term cloud observation records and a constant cloud observation system for higher resolution. Another case is the research conducted in Chad [13] calculating solar radiation of several cities located in different climatic zones in the country. In this study, models based on meteorological and astronomical parameters were used. These models are the Angstrom–Prescott model, the Allen model, and the Sabbagh model. These models were compared to data from the General Directorate of the National Meteorology of Chad to determine their performance and accuracy using a series of statistical metrics. The results obtained were discussed and evaluated to determine the model which predicts the best solar radiation for each city. Another interesting analysis is [14], which presented several solar radiation models developed for specific conditions in several countries and tested them in several cities of India. The results were then compared with local measurements of solar irradiance to determine which of the models has better performance in India.
The review presented in [9] analyzed several empirical models based on various approaches calculating global monthly solar radiation using meteorological parameters (sunlight, clouds, temperature, and other variables). These models were tested in the city of Yazd, Iran. The sunlight models, according to the statistical analysis in [9], were the most accurate to predict solar radiation in Iran. However, this approach requires sunlight records which are not always available. They concluded that other approaches like cloud-based or temperature can be used instead with acceptable accuracy. According to [9], the advantages of these approaches are inputs easily available everywhere. Other interesting research studies were presented in [15,16], where various mathematical models were compared to determine the optimal tilt angles to capture solar radiation on tilt surfaces for several locations worldwide. The performance of these models was analyzed for several latitudes and parameters to determine a tilt angle that provides the best performance for solar radiation for inclined surfaces. According to [16], the Koronakis and Liu–Jordan models were among the most accurate models to determine optimal surface inclination angles throughout the year. These studies are important for solar projects, which are searching for the best performance possible in their solar power installations.
A different approach to calculate daily solar radiation was presented by [17], where a model based on trigonometrical correlations was presented. This model was examined using statistical metrics against three existing models implemented on nine meteorological stations in China as reference data. Then, the proposed model was compared against 70 other meteorological stations in China to determine its errors and performance. According to [17], the proposed model resulted as simple, effective, and reliable to determine daily solar radiation in China. A similar research was presented by [18], where a series of new daily and monthly solar radiation models were presented. These models were examined in six cities in Iran against four existing models for several parameters to determine the accuracy of the new models using statistical metrics. Further, four extra cities in other countries were used to prove the accuracy of the new proposed models in latitudes worldwide. According to [18], the new models presented excellent results in general for all calculations. In addition, a highlight finding was that the new models presented a general better performance than the existing reference models.
A cloud cover model developed for Yakutia, Russia, was presented by [19], where solar radiation on tilt and horizontal surfaces was analyzed for production of hot water in permafrost regions in the Russian Arctic. The model utilized a five-year database of cloudiness in the city of Yakutsk. Hourly, daily, monthly, and annual radiations were calculated using the model. The results were compared against field testing measures in the region. In addition, the results obtained from the cloud cover model were analyzed for hot water production and storage using performance models for solar thermal systems. The objective was determining the cost efficiency of these systems throughout the year in Yakutia. The conclusion of [19] was that the information provided by the research in the permafrost region of Russia can be used to develop efficient thermal solar systems for reducing the carbon footprint of the region by using solar energy for water heating.
The review of cases provided in the previous paragraphs demonstrates the reliability and economic viability of mathematical models for the prediction of solar radiation worldwide. For the current inquiry, a model developed by the University of Tomsk [20] calculating solar global radiation on tilt surfaces was used as a base. This model was adapted to determine global solar radiation on a horizontal surface (ground-level solar radiation) in the state of Sonora in Mexico. This model was originally designed for high latitudes in Russia (over 55° N), hence a performance analysis at local latitudes in Sonora is required to ensure its accuracy. The main intention of the current research is to provide a wide distribution view of solar radiation in the state. This research work proposes this model as a useful tool to provide specific, low-cost, and reliable data for a wide and detailed view of solar radiation available in the state. Furthermore, it complements field measures from meteorological stations and satellite sources to guide the implementation of future solar power projects and promote the carbon footprint reduction due to electricity generation in the region.
The remainder of this paper is organized as follows. Section 2 presents the data sources utilized for calculations and for the determination of the accuracy in the model results. In addition, the criterion of the selection of field data is explained. Section 3 explains the methodology of the model and the strategies utilized to determine its performance. This section also describes a proposal of a geographic view of the results obtained from the model for later analysis and discussion. In Section 4, the results obtained are analyzed and discussed. Additionally, the utility, opportunities, and performance of the results in future solar projects are highlighted. Finally, Section 5 is the conclusion, where the main points of the research are presented. This section provides a global view of findings, key ideas, scopes, limitations, and suggestions about the current research work. In addition, a brief suggestion for future research is provided.

2. Data Sources

Two data sources were utilized in this research work to calculate the results obtained by the proposed model. Moreover, these sources were used as references to compare against the model and determine its performance and accuracy. These data resources are as follows: real-time measurements of the SMN meteorological stations taken from the Mexican National Water Commission (CONAGUA) [7], and satellite measurements from the NASA Surface Meteorology and Solar Energy (SSE) database [2].

2.1. SMN

The meteorological stations from the SMN selected were determined utilizing a selection criterion to avoid data gaps in the record and ensure proper data to analyze and compare against the model results. This criterion is based on the restrictions discussed in [8]. The selection conditions are as follows:
  • Stations with less than 85% of the data available are discarded.
  • Consistency in the data record; avoiding data gaps which may mean errors in measures reported by the stations. The previous fact may be due to the need for calibration of the station equipment [6] or due to shadows projected on the sensor because of obstacles that trigger a consistent daily loss of readings in data.
  • A range of G m < 1300 W/m2 for readings in stations over 1000 m above the sea level and G m < 1100 W/m2 for the rest of the stations.
Meteorological stations in the Mexican network (SMN) provide data every 10 min, and a daily average of total solar radiation in W/m2 is taken. Then, monthly and yearly averages are calculated and converted to kWh/m2. A sample of 5 years (from 2015 to 2019) was analyzed to calculate mean values for monthly and yearly data for comparison.
Five meteorological stations were considered for analysis. These stations represent 33.3% of the total of fifteen stations in the state. The selected stations are described in Table 1.
During the analysis and selection of the stations, some unexpected high values were found (measures of more than 1100 kWh/m2) in certain stations. This was especially so in El Pinacate, but it was selected for comparison due to its data availability.

2.2. NASA SSE Power Project

The NASA SSE Power project is an international satellite database available online which includes data from 1981 to present. This information includes meteorological and renewable resources such as measurements for solar, wind, temperature, moisture, and cloud index. This resource is reliable, and it is based on satellite observation and mathematical models for information worldwide [2,21]. The current research work utilized from this data source yearly averages (39 years) of the following parameters: all-sky solar radiation arriving at a horizontal surface, clearness index, and surface albedo.

3. Methodology

The methodology of the research is based on three phases. The first phase is the description of the mathematical model utilized. The second is the analysis of the model accuracy and performance based on statistical metrics for a case of study (Sonora), where the model is compared against real-time measurements of solar irradiance gathered from meteorological stations from the SME and the NASA SSE database. The last phase is the comparison and analysis of data obtained from the model for municipalities of Sonora, using heat maps to provide a geographical view of the solar resource distribution and viability of future solar power projects in the region.

3.1. Mathematical Model

The model is based on clustering the total radiation arriving at an inclined surface in three elements: direct, scattered, and reflected. The model inputs are surface albedo and clearness index obtained from the NASA SSE Power project database for the geographical location of interest [20]. All of these data are available using online tools such as the NASA SSE Power database [2] and Google Earth [10]. This allows the model to avoid investment in field equipment to collect solar irradiance measurements and provides great flexibility in data collection [20]. The model is explained in detail below.
Solar radiation arriving at a receiving surface depends upon factors such as geographical location (latitude and longitude), the day–night cycle, the date in the year, and the orientation of the receiving surface with respect to the Sun. These factors are called deterministic factors, and they vary yearly or daily and can be calculated by analytical formulas based on solar geometry [12,20].
According to [22], the day number ( N ) can be simplified from 365 days (full year) to 12 days to provide a representative monthly average, which presents an accurate estimation for solar radiation analysis [23]. This monthly average relationship is given in Table 2.
The declination angle of the Sun ( δ ) is the angular relationship between the solar position at solar noon and the equatorial plane, where: northwards (+) and southwards (−). This angle varies throughout the year ( ± 23.45 ° ) [24,25]. According to [15,26], it is given by the following equation:
δ = [ 23.45 sin ( ( 360 ) 284 + N 365 ) ] ( π 180 ) ,   Radians
The solar hour angle ( ω ) depicts the rotation of Earth around its axis, and each hour is equal to 15°, where: solar noon = 0, morning hours (−), and afternoon hours (+) [16,27]. According to [23], it can be obtained from
ω = { 15 [ ( T i m e   h r ) 12 ( D i f   G M T ) E o T ] } + Ψ ( π 180 ) ,   Radians
where: T i m e   h r is the time of the day in hours; D i f   G M T is the time zone hour difference with respect to the Greenwich meridian as the origin, for Sonora = −7; E o T is the equation of time; and Ψ is the longitude in degrees.
The equation of time ( E o T ) reflects the speed variation on Earth’s orbital movement throughout the year [11], and it is calculated from [23]
E o T = ( 9.87 sin 2 B ) ( 7.53 cos B ) ( 1.5 sin B ) 60 ,   Hours
with B determined by the equation
B = ( 360 365 ) ( N 81 ) ( π 180 )
The incidence angle ( θ ) portrays the angle between sunbeams and the normal on a receiving surface [11]. It can be obtained by
θ = cos 1 [ ( sin δ · sin φ · cos β ) + ( cos δ · sin φ · sin β · cos γ · cos ω ) ( sin δ · cos φ · sin β · cos γ ) + ( cos δ · sin β · sin γ · sin ω ) + ( cos δ · cos φ · cos β · cos ω ) ] ,   Radians
where: φ is the latitude given in radians; β is the slope angle of the receiving surface with respect to the horizontal; and γ is the surface azimuth angle, illustrating the variation in the projection on a horizontal plane of the normal to the surface from the local meridian, with south-facing considered as γ = 0, eastwards (−), and westwards (+), ( ± 180 ° ) [24]. In this paper, angles β and γ are considered 0 for this paper’s calculations, due to the studied surface being the ground itself and there being no inclination difference to consider. All the angles are considered in radians in the previous formula.
The solar zenith angle ( θ z ) is the angular relationship between the Sun’s rays and the vertical. For horizontal surfaces, θ Z = θ [11]. It is given by
θ z = cos 1 [ ( sin φ · sin δ ) + ( cos φ · cos δ · cos ω ) ] ,   Radians
The solar altitude angle ( h ) is the angle between the sunbeams and the horizontal [24]. It keeps the following relation:
h = ( π 2 ) θ Z ,   Radians
utilizing trigonometric relationships, the following equation can be obtained:
h = sin 1 [ ( sin φ · sin δ ) + ( cos φ · cos δ · cos ω ) ] ,   Radians
The solar azimuth angle ( A Z ) reflects the relationship between sunbeams and the horizontal plane, where east direction (−) and west (+) [11,28]. According to [23], this angle is determined by
A z = cos 1 [ ( sin h · sin φ ) sin δ cos h · cos φ ] ,   Radians
According to [29], the hourly diffuse coefficient ( r d ) is
r d = ( π 24 ) ( cos ω cos ω s s sin ω s s [ ω s s ( cos ω s s ) ] )
The sunset hour angle ( ω S S ) is taken from [11]
ω s s = cos 1 ( tan φ · tan δ ) ,   Radians
The hourly transparency coefficient ( r t ), according to [30], can be calculated by
r t = ( a + b cos ω ) r d
where the values of a and b can be obtained utilizing the following formulas:
a = 0.409 + 0.5016 [ sin ( ω s s π 3 ) ]
b = 0.6609 0.4767 [ sin ( ω s s π 3 ) ]
There are other factors called stochastic factors that also influence the solar radiation arriving at a receiving surface on Earth. These factors represent the state of the atmosphere and they are estimated based on an analysis of statistical models for long periods of observations and simulations. Examples of these factors are as follows: distribution of atmospheric gases, the density of particles in the air, air pressure, or cloud type at the site under study [12,20]. The clearness index ( K T ) and diffusion index ( K D ) are stochastic factors.
The clearness index ( K T ), according to [29], maintains the following relationship:
K T = H H 0
with K T taken from the NASA SSE database [2]; H 0 is the average daily extra-atmospheric insolation arriving at a horizontal surface [20,23]; and H is the average daily radiation arriving at a horizontal surface [11]. They can be calculated from the following formulas:
H 0 = G s c ( 24 π ) { 1 + [ 0.033 cos [ ( 360 N 365 ) ( π 180 ) ] ] } [ ( sin φ · sin δ · ω s s ) + ( cos φ · cos δ · sin ω s s ) ] ,   Wh / m 2
where: G S C is the solar constant = 1367 W/m2, representing the energy acting on a unit area per unit of time outside of the atmosphere at halfway from Earth to Sun [24]. All angles in the previous formula are in radians.
H = K T ( H 0 ) ,   Wh / m 2
The diffusion index ( K D ), as stated by [29], keeps the relationship below:
K D = H D H
The average daily diffuse radiation arriving at a horizontal surface H D is calculated from
H D = K D ( H ) ,   Wh / m 2
According to [21], K D can be determined by a series of polynomial equations conditioned by angles ( φ and ω s s ), where K D is determined based on K T .
K D = f ( K T )
Table 3 shows the equations and conditions to determine the diffusion index K D . All the angles shown in Table 3 are expressed in degrees.
The hourly total radiation arriving at a horizontal surface is divided into two components: direct and diffuse. The reflection part is not included for horizontal surfaces because there are not any angles that allow reflection between the surface and the horizontal. These components can be calculated as follows:
Total ( G H ):
G H = r t ( H ) ,   W / m 2
Diffuse ( G D H ):
G D H = r d ( H D ) ,   W / m 2
Direct ( G D ):
G D = G H G D H ,   W / m 2
According to [24], the hourly extra-atmospheric radiation arriving at a horizontal surface ( G 0 ) is given by
G 0 = G s c { 1 + [ 0.033 cos [ ( 360 N 365 ) ( π 180 ) ] ] } sin h ,   W / m 2
Finally, as specified by [20,31], the total radiation arriving at an inclined surface ( G ) can be calculated by the following equation:
G = [ G D ( cos θ cos θ Z ) ] + G D H { [ A i ( cos θ cos θ Z ) ] + [ ( 1 A i ) ( 1 + cos β 2 ) ] } + [ ( G H ) ( ρ ) ( 1 cos β 2 ) ] ,   W / m 2
where: ρ is the surface albedo obtained from the NASA SSE database, expressing the reflective capacity of the ground and varying according to different terrains (snow, sand, grass, etc.) and the seasons of the year [25]; and, as determined by [20,31], A i is the anisotropic index and can be obtained by the following formula:
A i = G D G 0

3.2. Determination of Accuracy and Performance

The accuracy and performance of the model are determined by comparing data from several meteorological stations in Sonora. The proposed model is compared against two data sources: real-time measurements of the SMN taken from CONAGUA [7], and the NASA SSE database [2].

3.2.1. Statistical Metrics

Eight statistical metrics are proposed for the analysis and determination of the accuracy and performance of the model. These metrics are described as follows:
Mean absolute error ( M A E ) provides the average error magnitude, but without the error direction. The smaller, the better [32]. It can be calculated by
M A E = 1 n i = 1 n | X i Y i |
where: n is the number of months in the year; i is the number of the analyzed month; Y is the reference value (from field data or from the NASA SSE database); and X is the calculated value (from the model). Both X and Y are G values in kWh/m2.
Mean bias error ( M B E ) provides the trend of the average error, which reflects the performance of the analyzed model, where underestimation (−), overestimation (+), and exactitude (0). The closer to 0, the better [32]. It is given by
M B E = 1 n i = 1 n ( X i Y i )
Root mean square error ( R M S E ) represents the standard deviation of the calculated errors. The smaller the value, the better the accuracy [8,13]. It can be obtained by
R M S E = 1 n i = 1 n ( X i Y i ) 2
Mean percentage error ( M P E ) determines the error performance, where the range ± 10 % is considered acceptable [13]. It is provided by
M P E = [ 1 n i = 1 n ( X i Y i Y i ) ] ( 100 )
Relative percentage error ( R P E ), with a range of ± 10 % considered as acceptable [9]. It is calculated from
R P E = ( X i Y i Y i ) ( 100 )
Correlation coefficient ( r ) measures the linear correlation between two variables in a range of ± 1 , where: total positive (1), total negative (−1), and no linear correlation (0) [17,33]. It is given by
r = i = 1 n ( X i X ¯ ) ( Y i Y ¯ ) i = 1 n ( X i X ¯ ) 2 i = 1 n ( Y i Y ¯ ) 2
with X ¯ and Y ¯ as yearly values of total radiation for calculated data using (25) and reference data, respectively.
Coefficient of determination ( R 2 ) determines the proximity between the regression line of the calculated values and the real values. The closer to 1, the more accuracy [9,32]. It can be calculated using the following formula:
R 2 = 1 ( i = 1 n ( Y i X i ) 2 i = 1 n ( Y i Y ¯ ) 2 )
The t-Student distribution test ( t ) is a common statistical analysis of data distribution. It can be used to determine if the values obtained from the model are statistically significant or not. The smaller the t-value, the better the performance. The statistical significance of the model is considered based on a t-distribution table of critical values, where the critical value is obtained using a confidence level ( α ) and a degree of freedom ( d f ). The t-value calculated must be less than the critical value to be statistically significant [9,34]. The t-value can be calculated by
t = ( d f ) ( M B E ) 2 ( R M S E ) 2 ( M B E ) 2
where:
d f = n 1

3.2.2. Performance and Accuracy Visualization Tools

As a visual comparison tool, a series of graphs were utilized to exhibit the behavior of monthly solar irradiance data throughout the year between the three sources of data: the SMN, the NASA SSE database, and the model for the selected meteorological stations in Section 2. This visualization provides a clear analysis to evaluate the performance and accuracy of the model.

3.3. Geographical Analysis

A series of monthly heat maps that show the relation of solar irradiance against average maximum and minimum temperatures during the year for the 72 municipalities in Sonora are proposed. Municipalities were chosen and analyzed as points of interest because they can provide an idea of the potential for solar residential projects in urban areas and a general view of solar potential across the state. In addition, the temperature was selected as a comparison parameter because it strongly affects the efficiency of most photovoltaic panels in the market, as well as the design of PV systems. Important to note is that hot temperatures decrease the efficiency in the energy production of photovoltaic modules, as the higher the temperature, the higher the loss [35,36]. This has been demonstrated by research such as the Solar Advisor Model (SAM) [37]. On the other hand, minimum temperatures determine the design and sizing of a photovoltaic array. The voltage produced by a solar panel increases with a decreasing temperature [38]. The voltage of every panel is summed depending on the configuration (series or parallel) of the string connected to an inverter (component to convert to AC). The maximum DC voltage input of the inverter limits the number of panels connected per string ( V D C   s t r i n g < V DC   max i n v e r t e r ) [11,37]. The temperatures utilized for this analysis were taken from the NASA SSE database for the coordinates of each municipality.
The proposed heat maps were used as a visual tool to show the potential areas to develop solar power projects in Sonora, and to discuss their viability in power production. Moreover, monthly maximum, minimum, and average temperatures together with the values of the model and the NASA SSE database for solar irradiance averages per month and annually by the municipality are shown in Appendix A.

4. Results and Discussion

The total solar irradiance in kWh/m2 is presented and analyzed for three data sources (the results calculated by the proposed model, the on-site measurements from the stations of the Mexican meteorological network “SMN”, and the satellite measurements from the solar irradiance records from the NASA SSE database) for each of the five stations selected. In addition, the statistical metrics previously described in Section 3 were calculated for each station to determine the performance among the data. Among these statistical metrics, namely mean absolute and bias errors (MAE and MBE), together with root mean square error (RMSE), the smaller their value, the better the performance. For the mean and relative percentage errors (MPE and RPE), the acceptable range is ± 10 % . The correlation coefficient (r) establishes the linear correlation of the data (0 = no correlation, and −1 or +1 as negative or positive linear correlations, respectively). The coefficient of determination (R2) indicates the proximity between the calculated and reference values analyzed. In this metric, the closer to 1, the better. For the last statistic metric, the t-Student distribution test (t), the smaller, the better, and it must respect the rule of being smaller than the critical value selected for the analysis to be statistically significant. In this research, a critical value of 4.025 was taken for a confidence level of 99% (0.001) and 11 degrees of freedom. Three comparisons were made between the data sources to show the statistical behavior analysis (NASA SSE-Model, SMN-Model, and SMN-NASA SSE).

4.1. NASA SSE-Model

This comparison tests the accuracy of the proposed model with respect to the solar irradiance records of NASA SSE used as reliable reference data. In this case, the results showed acceptable values for M A E , M B E , and R M S E . All of them were close to 0, with a slightly higher value in Yécora station. For R P E and M P E , all the values were acceptable (in the range of ± 10 % ). The metric r showed a positive linear correlation in all the stations. For R 2 , the values were considerably close to 1, which means good performance in the series of the analyzed values. This metric was slightly lower in the case of Yécora. According to the selected t-Student critical value (4.025), all the stations presented a statistical significance, except for Yécora station. With these results, a notably accurate approximation was shown between the NASA SSE and the model values. A slightly higher difference was shown in the case of Yécora station which resulted in a t-test parameter out of tolerance. The results of this analysis are presented in Table 4.

4.2. SMN-Model

The analysis in this comparison examines the proposed model against the on-site measurements of the stations in the Mexican meteorological network (SMN) utilized as valid data references in this case. In this comparison, the model showed good results for M A E , M B E , and R M S E in the cases of Hermosillo and Yécora, and higher errors for the other stations, with El Pinacate as the worst of the analyzed stations. For R P E and M P E , Hermosillo and Yécora were acceptable in all their values. Yécora had a better average performance than Hermosillo. On the other hand, the rest of the stations showed several months out of tolerance (less than −10%) for R P E . The M P E value was out of tolerance for the El Pinacate case, but acceptable in the cases of Caborca and Nogales. All the stations showed a positive linear correlation according to the r parameter. The metric R 2 exhibited good performance for Hemosillo and Yécora, slightly less for Caborca and Nogales, and the worst for El Pinacate. For the t-test, just Hermosillo and Yécora were statistically significant, and the rest of the stations were not. The statistical metrics exhibited a good approximation between the values of the SMN and the model in the cases of Hermosillo and Yécora. On the other hand, Caborca and Nogales had a less accurate performance. The highest error was found in El Pinacate. Table 5 provides the results of the statistical analysis.

4.3. SMN-NASA SSE

The last comparison evaluates the two data references used in the two previous comparisons (the on-site measures from the meteorological stations of the SMN and the satellite solar irradiance records in the NASA SSE database). The comparison between the two references seeks to clarify and find discrepancies in the results of the previous two cases, where they were compared against the proposed model. The results presented in this comparison were good for the M A E , M B E , and R M S E values for Hermosillo and Yécora, with slightly higher errors in Caborca and Nogales, and the highest error in El Pinacate. For R P E , Caborca, Hermosillo, and Yécora showed a good performance, except for November in Hermosillo and August in Caborca, which were slightly out of tolerance. The rest of the stations presented several months with unacceptable values, with El Pinacate as the worst on average. The M P E was acceptable for four of the stations and just El Pinacate was out of tolerance. The parameter r revealed a linear correlation in all the cases. R 2 showed a slightly acceptable dispersion in four of the stations but showed a higher dispersion in the case of El Pinacate. The t-test exhibited statistical significance for Hermosillo and Yécora, but not for the rest of the stations. The results of this comparative analysis are shown in Table 6.

4.4. General Discussion and Analysis

The results obtained reveal a close approximation in all the cases for the comparison between the model and the NASA SSE values. This means a good performance and small errors involved between both. On the other hand, the analysis also reveals slight errors in the cases of Hermosillo and Yécora for the comparisons SMN-Model and SMN-NASA SSE, with slightly wider errors in the cases of Caborca and Nogales for certain months, and the worst of the cases was El Pinacate with several months with wider errors. Further, in El Pinacate, the difference was consistent and may mean a perturbation in the values obtained from the station. However, further research of the station equipment is recommended to provide a more precise explanation. The previous trends were consistent for both comparisons, which means a significant difference between the model and the NASA SEE with respect to the SMN.

4.4.1. Visualization of Performance and Accuracy

The performance of the model is presented utilizing a visual comparison of the data obtained for solar irradiance in kWh/m2 from each data source (the SMN, the NASA SSE, and the model) for each meteorological station analyzed. Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 present this analysis.
In Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6, the statement described by the statistical analysis was visually reinforced. A close trajectory for the values of NASA SSE (red) and the model (green) is visible in all the stations. However, the SMN (blue) maintains a close trajectory with small errors for the stations of Hermosillo and Yécora. These were the two stations with better performance in the model approximation for the other two references (NASA SSE and SMN). In the graphs, it can be seen in the case of Yécora that the model is closer to the SMN than NASA SEE. This explains why, in this case, the statistical metrics presented a better performance of the model with respect to the SMN compared to NASA SSE. The errors were wider in the case of Caborca and Nogales, with significant gaps between the red and green lines with respect to the blue line. The highest error was in El Pinacate with separate trajectories between the SMN and the other two sources of data. This means a significant standard deviation between NASA SSE and the model with respect to the SMN values in this station. Additionally, it merits noticing that the errors described by R P E and M B E in the statistical analysis described the direction of the error and magnitude for each month, where the SMN was taken as the reference for both. When the NASA SSE or the model values were smaller than the SMN, R P E and M B E were negatives, otherwise, they were positives. The magnitude of these parameters represents the gaps between the trajectories of the three data sources.
As a result of the obtained values in the comparative analysis, it can be stated that the model was quite accurate as an approximation of the values obtained from the NASA SSE database for latitudes of 32–26° N. The model had originally been tested at high latitudes in Russia at over 50° N. These results may imply that the model can be applied to any geographical point and obtain a close approximation of the total solar radiation available. However, more research in other world regions for latitudes further south is needed to ensure this statement.

4.4.2. Geographical Overview of Solar Resource and Opportunities of Solar Power in Sonora

A geographical information system (GIS) comparison analysis between solar total radiation calculated from the model and temperatures (maximum/minimum) for all municipalities in Sonora based on a series of heat maps is presented as a visualization tool of solar potential in Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18. In these figures, the months of the year are represented in order from January to December, where Figure 7 represents January, Figure 8 represents February, and following this format until Figure 18 representing December. Further, each of the figures is separated into three GIS representations of the state of Sonora, where the first to the left is the average solar irradiance (kWh/m2), calculated by the model. The central heat map represents the average maximum temperatures (°C). Finally, the heat map on the right of each figure illustrates the average minimum temperatures (°C). A scale of colors was utilized to represent the values for each municipality in each GIS representation in Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18. The coordinates together with the results of solar irradiance calculated by the model and from the NASA SSE database (for comparison), temperatures, surface albedo, and transparency index throughout the year for each municipality in the state can be found in Appendix A.
The geographical analysis showed high solar irradiance in Sonora throughout the year, with values over 2.99 kWh/m2 on average. The regions with constant higher values were the north-central, central, east, and southern regions of the state. However, some of these areas were also the hottest with average maximum temperatures over 40 °C during April–September, which may mean a significant economic cost due to the need for more modules or cooling systems to compensate for the loss of efficiency due to high temperatures. On the other hand, average minimum temperatures under 0 °C were presented in the northeast and eastern regions during December–February. This represents a good potential for solar development in the eastern region of the state because of less severe maximum temperatures and good solar irradiance on average. However, extreme average minimum temperatures must be considered during the design and sizing of solar photovoltaic projects. The main problem in the eastern region is the mountainous terrain, where any solar project must be carefully analyzed due to the space required. Finally, the western and northwest regions of the state were also worth consideration because, even though there are fewer solar resources compared to the central region, the solar irradiance in these regions is still high on average and the terrain is flatter than in other regions of the state, which facilitates the planning and installation of large-scale solar harvesting projects. Temperatures are also not so extreme, so better efficiency values in photovoltaic energy production are expected.
The results obtained in this GIS analysis provide a view and perspective of a great zone of opportunity in the northwest territory of Mexico for energy production based on the abundant solar resource available. This may represent a tremendous potential of local development of municipalities in the region and their decarbonization for electricity and heating production, which may mean a big reduction in its carbon footprint. This would imply a great economic development and possible energy exportation out of the region either to the United States or other regions of Mexico.

5. Conclusions

The accuracy of the model for the latitudes of the state of Sonora was proved and analyzed under several statistical metrics comparing it against real-time measurements from the stations of the Mexican meteorological network (SMN) and the satellite solar irradiance records from the NASA SSE database. The SMN meteorological stations available for a comparative analysis in the state were only 33.3% of the total. The rest of the stations presented important data gaps for the period 2015–2019 which discarded them as reliable reference sources. The performance obtained in the statistical analysis demonstrated satisfactory accuracy of the model for latitudes of 32–26° N. However, some differences were observed between NASA SSE and the model compared to real-time data from the SMN for certain meteorological stations. This may imply measurement errors or a need for calibration, but more research of equipment in the SMN stations must be conducted to ensure this statement.
A series of heat maps was utilized for comparing solar radiation against maximum and minimum temperatures in the municipalities in Sonora state. This geographical analysis highlighted that Sonora is a region with high radiation values in a range from 2.99 kWh/m2 to 8.07 kW/m2 during the year for all the municipalities calculated. The advantages and disadvantages were discussed and identified based on the range of annual temperatures for each region in the state. This analysis revealed regions of opportunities for decision-makers and stakeholders in different organizations, both public and private, interested in reducing the carbon footprint, providing an acceptable and wide vision for planning photovoltaic power projects in the region.
The extensive set of tests and results remark the non-complexity nature and accuracy of the applied mathematical model. It requires three inputs (surface albedo, clearness index, and geographical coordinates) for the calculation of solar radiation in a convenient and economically viable way. This research work concludes that the applied model could be used as a support or reference for the calculation of solar resources in a county, a state, or large areas of countries for decision making related to sustainable development policies aimed at solar energy harvesting systems and carbon footprint reduction.
Finally, as future work, it is proposed to continue calculating the amount of solar radiation, but in other latitudes of the planet to validate the model and the herein results. The outcomes of this work reinforce the potential of this approach as a tool to be applied during solar energy planning by local governments.

Author Contributions

Conceptualization, E.A.E.-V. and V.H.B.; data curation, E.A.E.-V.; formal analysis, E.A.E.-V.; funding acquisition, V.H.B., L.C.F.-H. and J.d.-J.L.-S.; investigation, E.A.E.-V. and V.H.B.; methodology, E.A.E.-V. and V.H.B.; project administration, E.A.E.-V. and V.H.B.; resources, E.A.E.-V. and V.H.B.; software, E.A.E.-V. and S.G.O.; supervision, V.H.B., L.C.F.-H. and J.d.-J.L.-S.; validation, E.A.E.-V., V.H.B. and S.G.O.; visualization, E.A.E.-V., V.H.B. and S.G.O.; writing—original draft, E.A.E.-V.; writing—review and editing, E.A.E.-V., V.H.B., S.G.O., L.C.F.-H. and J.d.-J.L.-S. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by CampusCity Initiative from the School of Engineering and Sciences at the Tecnologico de Monterrey.

Acknowledgments

A sincere thanks goes to the SMN for facilitating the data measured at its meteorological stations to carry out the statistical analysis presented in this research paper.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Nomenclature

Symbols
aHourly transparency coefficient first variable (dimensionless)
AiAnisotropic index (dimensionless)
AzSolar azimuth (angle)
BEquation of time coefficient (dimensionless)
bHourly transparency coefficient second variable (dimensionless)
dfT-Student degrees of freedom (dimensionless)
Dif GMTTime zone hour difference with respect to Greenwich meridian (hours)
EoTEquation of time (hours)
GTotal radiation arriving at an inclined surface (W/m2)
G0Hourly extra-atmospheric radiation arriving at a horizontal surface (W/m2)
GDHourly direct radiation arriving at a horizontal surface (W/m2)
GDHHourly diffuse radiation arriving at a horizontal surface (W/m2)
GHHourly total radiation arriving at a horizontal surface (W/m2)
GmTotal solar irradiance measured on-site (W/m2)
GSCSolar constant = 1367 W/m2
hSolar altitude angle
HAverage daily radiation arriving at a horizontal surface (Wh/m2)
H0Average daily extra-atmospheric insolation arriving at a horizontal surface (Wh/m2)
HDAverage daily diffuse radiation arriving at a horizontal surface (Wh/m2)
iNumber of the analyzed month
KDDiffusion index (dimensionless)
KTClearness index (dimensionless)
Max. TMaximum average temperature (°C)
Mean TMean average temperature (°C)
Min. TMinimum average temperature (°C)
nNumber of months in the year
NDay number
rCorrelation coefficient (dimensionless)
R2Coefficient of determination (dimensionless)
rdHourly diffuse coefficient (dimensionless)
rtHourly transparency coefficient (dimensionless)
tT-Student distribution test (dimensionless)
Time hrTime of the day in hours
VDC max inverterMaximum DC voltage input of an inverter (Volts)
VDC stringDC voltage per string connected to an inverter (Volts)
XiReference value for statistical analysis
YiCalculated value for statistical analysis
Greek letters
ω S S Sunset hour angle
θ z Solar zenith angle
δDeclination angle
ΨLongitude
ωSolar hour angle
β Slope angle of the receiving surface with respect to the horizontal
γ Surface azimuth angle
θ Incidence angle
ρ Surface albedo (dimensionless)
φ Latitude
Abbreviations
ACAlternating current
CCSCenter for Climate Strategies
CEDESCommission of Ecology and Sustainable Development of the State of Sonora
CEMIE-SOLMexican Centre of Innovation in Solar Energy
CONACYTNational Council for Science and Technology
CONAGUAMexican National Water Commission
DCDirect current
EMAsAutomatic meteorological stations
ESIMEsSynoptic meteorological stations
GHGsGreenhouse gas emissions
GISGeographical information system
MAEMean absolute error (dimensionless)
MBEMean bias error (dimensionless)
MPEMean percentage error (%)
NASANational Aeronautics and Space Administration
NASA SSENASA Surface Meteorology and Solar Energy
PVPhotovoltaic
RMSERoot mean square error (dimensionless)
RPERelative percentage error (%)
SAMSolar Advisor Model
SENERSecretary of Energy
SMNMexican National Weather Service
SOLARGISGlobal Solar Atlas
UNAMNational Autonomous University of Mexico
UTCCoordinated Universal Time
WMOWorld Meteorological Organization

Appendix A. Calculations of Solar Irradiance Parameters for All Municipalities in Sonora

Table A1 presents the solar irradiance in kWh/m2 calculated using the proposed model and from the satellite NASA SSE database for all the 72 municipalities in the state of Sonora. Further, surface albedo ( ρ ) and transparency index (Kt) are shown. These two parameters are dimensionless. Other important parameters illustrated are monthly maximum, mean, and minimum average temperatures in °C for each municipality.
Table A1. Solar irradiance in kW/m2 from the model, the NASA SSE database, temperatures, surface albedo, and transparency index for each municipality in Sonora.
Table A1. Solar irradiance in kW/m2 from the model, the NASA SSE database, temperatures, surface albedo, and transparency index for each municipality in Sonora.
MunicipalityLatitudeLongitudeParameterMonthAnnual Average
JanFebMarAprMayJunJulAugSeptOctNovDec
Aconchi29.82° N110.26° WS. Albedo0.170.170.170.170.180.190.180.170.170.160.160.16
T. index (KT)0.630.640.690.710.70.680.60.590.620.650.650.63
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
Model3.74.586.017.217.747.746.726.185.714.944.013.465.67
Max. T (°C)25.0527.1232.2837.5742.4746.5540.6437.8538.3835.7729.7824.57
Mean T (°C)11.613.2417.0121.6526.5731.3929.7528.0926.7122.116.0711.44
Min. T (°C)3.444.065.999.0513.2718.6521.5120.8518.4813.017.583.69
Agua Prieta31.31° N109.66° WS. Albedo0.150.160.160.160.170.180.170.160.160.150.150.14
T. index (KT)0.620.630.680.710.70.670.550.560.620.650.650.62
NASA SSE3.534.365.927.27.737.586.185.775.614.823.893.275.49
Model3.494.385.817.167.757.666.175.855.634.813.853.255.48
Max. T (°C)18.4521.3126.4131.9237.8943.6639.9837.6536.1630.6623.418.01
Mean T (°C)5.677.8611.7416.6222.4628.2627.826.0923.4217.1110.155.5
Min. T (°C)−2.28−1.330.84.159.315.0518.4817.5714.27.811.58−2.11
Álamos27.03° N109.03° WS. Albedo0.120.120.120.120.130.130.140.130.130.110.110.12
T. index (KT)0.640.650.70.70.70.670.590.580.620.640.660.63
NASA SSE4.074.946.397.267.87.566.596.095.885.184.383.775.83
Model4.044.96.287.197.727.556.566.15.845.084.353.755.78
Max. T (°C)28.8130.7134.9939.3643.2345.1838.6236.2436.6136.7632.8228.6
Mean T (°C)16.2617.6720.6824.5828.5732.2230.0828.928.0625.2720.4616.5
Min. T (°C)8.358.9410.4513.2316.8621.9423.8823.4621.9617.6212.649.01
Altar30.74° N111.85° WS. Albedo0.170.180.170.180.190.20.190.180.180.170.170.17
T. index (KT)0.630.640.680.70.690.660.580.580.620.640.640.63
NASA SSE3.654.546.017.27.557.486.495.935.664.883.993.45.57
Model3.64.55.857.087.637.536.56.065.664.783.853.375.53
Max. T (°C)24.9427.0131.9736.9742.6447.4944.4342.4541.5536.5529.8524.4
Mean T (°C)11.8213.4117.0721.3626.6531.731.9230.7228.5322.6316.3411.55
Min. T (°C)3.984.586.589.4113.7618.7322.722.3919.613.558.174.04
Arivechi29° N109.47° WS. Albedo0.140.140.140.150.150.160.150.140.140.140.130.13
T. index (KT)0.640.650.70.720.730.690.60.610.640.660.660.64
NASA SSE3.884.726.327.528.127.966.786.365.965.154.243.65.88
Model3.844.736.157.348.077.836.76.45.935.084.153.65.82
Max. T (°C)24.5226.3631.3636.4741.3244.7137.9835.2536.1634.4429.224.2
Mean T (°C)11.3512.9816.6221.1926.0830.5928.2526.6225.6521.5915.8911.39
Min. T (°C)3.033.815.648.7513.1318.4820.7619.9517.9512.847.393.49
Arizpe30.35° N110.2° WS. Albedo0.160.160.160.160.170.180.170.160.160.150.150.15
T. index (KT)0.630.640.690.720.720.680.580.580.630.650.660.63
NASA SSE3.714.586.157.397.977.856.546.065.84.994.093.465.72
Model3.644.535.977.297.967.756.56.075.774.894.013.415.65
Max. T (°C)22.7224.9730.0135.3740.5945.2940.2537.7137.5533.7727.4522.28
Mean T (°C)9.4811.2114.9419.6224.8630.0328.8227.2225.4320.1313.889.32
Min. T (°C)1.3523.947.0611.5617.0720.1919.5116.8110.975.361.58
Átil30.84° N111.6° WS. Albedo0.170.180.170.180.190.20.190.180.180.170.170.17
T. index (KT)0.630.640.680.70.690.660.580.580.620.640.640.63
NASA SSE3.654.546.017.27.557.486.495.935.664.883.993.45.57
Model3.594.495.857.087.637.536.56.065.664.783.843.365.53
Max. T (°C)24.9427.0131.9736.9742.6447.4944.4342.4541.5536.5529.8524.4
Mean T (°C)11.8213.4117.0721.3626.6531.731.9230.7228.5322.6316.3411.55
Min. T (°C)3.984.586.589.4113.7618.7322.722.3919.613.558.174.04
Bacadéhuachi29.8° N109.16° WS. Albedo0.140.140.140.150.150.160.150.140.140.140.130.13
T. index (KT)0.640.650.70.720.730.690.60.610.640.660.660.64
NASA SSE3.884.726.327.528.127.966.786.365.965.154.243.65.88
Model3.764.666.097.318.077.856.726.395.895.014.073.525.78
Max. T (°C)2324.8929.8234.9640.0744.0738.3635.936.1533.327.6622.6
Mean T (°C)10.0111.6715.2919.912529.7827.9726.425.0420.4114.59.98
Min. T (°C)1.732.544.367.5212.0417.4319.9919.2116.9211.545.972.13
Bacanora28.98° N109.42° WS. Albedo0.140.140.140.150.150.160.160.150.140.130.130.13
T. index (KT)0.640.650.70.720.720.680.590.590.620.650.660.63
NASA SSE44.836.387.458.017.736.656.25.95.184.323.75.86
Model3.844.736.157.347.967.726.596.195.755.014.153.555.75
Max. T (°C)24.5726.3231.0235.8940.3943.4136.5933.8534.6933.6929.0224.34
Mean T (°C)11.813.3216.7121.0625.6629.9627.5225.9325.0521.4816.1711.92
Min. T (°C)3.774.446.079.0713.1818.4620.5519.7817.9513.188.014.28
Bacerac30.35° N108.94° WS. Albedo0.140.140.140.140.150.160.150.140.140.140.130.14
T. index (KT)0.630.640.690.710.70.660.550.560.620.650.660.63
NASA SSE3.74.516.127.37.827.526.215.835.74.954.073.425.6
Model3.644.545.977.197.747.526.165.865.684.894.013.415.55
Max. T (°C)17.5419.5123.7328.3233.4138.0734.6633.232.1127.9522.0617.3
Mean T (°C)6.067.7210.9315.2220.2225.2124.5223.3621.3316.2110.316.08
Min. T (°C)−1.42−0.650.863.938.4213.6416.4715.8813.297.832.32−1.1
Bacoachi30.63° N109.98° WS. Albedo0.150.150.150.150.160.160.150.150.140.140.140.14
T. index (KT)0.640.640.70.720.720.680.570.580.620.650.660.63
NASA SSE3.744.556.177.448.017.816.456.055.7854.123.455.71
Model3.674.516.037.297.977.766.396.065.674.873.983.385.63
Max. T (°C)20.0322.6627.6833.1638.8644.340.0737.7736.6731.7224.8719.63
Mean T (°C)6.999.0112.817.6423.2928.8628.0326.3924.0418.111.446.86
Min. T (°C)−1.19−0.321.6959.9715.6818.9318.114.988.782.77−0.95
Bácum27.54° N110.12° WS. Albedo0.070.070.070.080.090.090.090.080.080.070.070.07
T. index (KT)0.650.660.710.710.730.710.650.630.650.670.660.63
NASA SSE4.144.996.57.448.128.167.366.696.255.444.443.796.11
Model4.054.936.347.288.058.027.246.236.15.284.33.75.96
Max. T (°C)29.0430.8235.0839.242.9546.0141.0738.939.5438.5933.5628.88
Mean T (°C)17.3418.6621.6325.2829.1933.2332.2331.1630.4527.2421.9217.65
Min. T (°C)10.1710.7612.3915.0218.5823.4926.0525.8624.4919.9914.7410.88
Banámichi30.01° N110.24° WS. Albedo0.160.160.160.160.170.180.170.160.160.150.150.15
T. index (KT)0.630.640.690.720.720.680.580.580.630.650.660.63
NASA SSE3.714.586.157.397.977.856.546.065.84.994.093.465.72
Model3.684.575.997.317.967.746.496.075.794.924.053.445.67
Max. T (°C)22.7224.9730.0135.3740.5945.2940.2537.7137.5533.7727.4522.28
Mean T (°C)9.4811.2114.9419.6224.8630.0328.8227.2225.4320.1313.889.32
Min. T (°C)1.3523.947.0611.5617.0720.1919.5116.8110.975.361.58
Baviácora29.7° N110.19° WS. Albedo0.170.170.170.170.180.190.180.170.170.160.160.16
T. index (KT)0.630.640.690.710.70.680.60.590.620.650.650.63
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
Model3.714.596.017.217.747.746.716.185.714.954.023.485.67
Max. T (°C)25.0527.1232.2837.5742.4746.5540.6437.8538.3835.7729.7824.57
Mean T (°C)11.613.2417.0121.6526.5731.3929.7528.0926.7122.116.0711.44
Min. T (°C)3.444.065.999.0513.2718.6521.5120.8518.4813.017.583.69
Bavispe30.48° N108.92° WS. Albedo0.140.140.140.140.150.160.150.140.140.140.130.14
T. index (KT)0.630.640.690.710.70.660.550.560.620.650.660.63
NASA SSE3.74.516.127.37.827.526.215.835.74.954.073.425.6
Model3.634.525.967.197.747.536.165.865.674.8843.395.54
Max. T (°C)17.5419.5123.7328.3233.4138.0734.6633.232.1127.9522.0617.3
Mean T (°C)6.067.7210.9315.2220.2225.2124.5223.3621.3316.2110.316.08
Min. T (°C)−1.42−0.650.863.938.4213.6416.4715.8813.297.832.32−1.1
Benito Juárez27.28° N110.01° WS. Albedo0.070.070.070.080.090.090.090.080.080.070.070.07
T. index (KT)0.650.660.710.710.730.710.650.630.650.670.660.63
NASA SSE4.144.996.57.448.128.167.366.696.255.444.443.796.11
Model4.084.966.367.288.058.017.236.636.115.34.333.736.01
Max. T (°C)25.8526.9329.9833.2236.6439.5637.2236.2536.5435.4730.8426.46
Mean T (°C)18.1518.8721.0223.9927.5431.0931.2931.0530.5627.9923.1819.12
Min. T (°C)13.4813.7415.0417.4120.7424.8727.3127.6626.7923.3518.5114.72
Benjamín Hill29.88° N111.35° WS. Albedo0.150.160.160.160.170.180.170.160.160.150.150.15
T. index (KT)0.620.640.680.690.690.660.590.590.620.640.640.62
NASA SSE3.724.66.097.097.547.476.66.045.784.944.073.475.62
Model3.634.585.917.017.637.516.66.185.74.863.943.45.58
Max. T (°C)26.4228.5233.4338.4243.3447.6842.8440.4440.3136.7630.7725.83
Mean T (°C)13.1314.6818.2822.5827.3432.1731.1229.6828.0623.0817.2812.87
Min. T (°C)5.35.817.710.4814.3919.522.5822.1319.7914.229.185.42
Caborca30.34° N112.55° WS. Albedo0.10.110.110.120.120.120.120.110.110.10.10.1
T. index (KT)0.620.630.670.690.690.670.610.610.630.630.620.6
NASA SSE3.594.475.937.087.617.686.886.315.784.793.873.295.61
Model3.594.475.86.997.637.646.836.385.774.743.873.255.58
Max. T (°C)22.7823.927.1930.634.5238.4139.0938.8838.334.1828.2923.23
Mean T (°C)14.8415.7118.221.3125.1329.0731.4831.6730.325.6420.1515.59
Min. T (°C)10.0710.3911.9814.3517.6921.5325.926.632520.1915.1911.04
Cajeme27.58° N110.38° WS. Albedo0.070.070.070.080.090.090.090.080.080.070.070.07
T. index (KT)0.650.660.710.710.730.710.650.630.650.670.660.63
NASA SSE4.144.996.57.448.128.167.366.696.255.444.443.796.11
Model4.054.936.347.288.058.027.246.626.095.284.33.75.99
Max. T (°C)29.0430.8235.0839.242.9546.0141.0738.939.5438.5933.5628.88
Mean T (°C)17.3418.6621.6325.2829.1933.2332.2331.1630.4527.2421.9217.65
Min. T (°C)10.1710.7612.3915.0218.5823.4926.0525.8624.4919.9914.7410.88
Cananea31° N110.35° WS. Albedo0.160.160.160.170.180.180.170.160.160.150.160.15
T. index (KT)0.620.630.680.710.710.680.550.560.620.640.650.62
NASA SSE3.554.45.977.257.97.776.285.785.644.793.933.295.55
Model3.524.415.847.177.867.776.175.855.654.763.893.295.52
Max. T (°C)18.7121.3726.3231.6837.6243.2739.6537.4436.0830.6423.5718.26
Mean T (°C)5.937.9211.6616.3822.1427.8827.5325.9223.3317.0510.295.73
Min. T (°C)−1.94−1.170.834914.6718.317.4914.157.81.82−1.81
Carbó29.68° N110.98° WS. Albedo0.170.170.170.170.180.190.180.170.170.160.160.16
T. index (KT)0.630.640.690.710.70.680.60.590.620.650.650.63
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
Model3.714.66.027.217.747.746.716.185.714.954.023.485.67
Max. T (°C)25.9328.0233.0338.0842.846.9641.8139.4139.5936.5330.5425.43
Mean T (°C)12.6514.2317.9122.3427.0631.8830.6329.1727.6722.9317.0112.46
Min. T (°C)4.765.297.1910.0714.0319.2822.2721.7519.43148.84.97
Cucurpe30.33° N110.72° WS. Albedo0.160.160.160.160.170.180.170.160.160.150.150.15
T. index (KT)0.630.640.690.720.720.680.580.580.630.650.660.63
NASA SSE3.714.586.157.397.977.856.546.065.84.994.093.465.72
Model3.654.545.977.37.967.756.56.075.774.894.013.415.65
Max. T (°C)23.6725.8330.823641.0745.6640.9438.6638.4734.6328.3623.22
Mean T (°C)10.5312.1515.8420.3425.4330.4929.4628.0226.2621.0114.8810.35
Min. T (°C)2.613.155.088.0412.3417.6620.8320.3117.7111.976.612.8
Cumpas30.34° N110.5° WS. Albedo0.160.160.160.160.170.180.170.160.160.150.150.15
T. index (KT)0.630.640.690.720.720.680.580.580.630.650.660.63
NASA SSE3.714.586.157.397.977.856.546.065.84.994.093.465.72
Model3.654.545.977.297.967.756.56.075.774.894.013.415.65
Max. T (°C)22.7224.9730.0135.3740.5945.2940.2537.7137.5533.7727.4522.28
Mean T (°C)9.4811.2114.9419.6224.8630.0328.8227.2225.4320.1313.889.32
Min. T (°C)1.3523.947.0611.5617.0720.1919.5116.8110.975.361.58
Divisaderos29.62° N109.48° WS. Albedo0.140.140.140.150.150.160.150.140.140.140.130.13
T. index (KT)0.640.650.70.720.730.690.60.610.640.660.660.64
NASA SSE3.884.726.327.528.127.966.786.365.965.154.243.65.88
Model3.784.676.117.328.077.856.716.395.95.034.093.545.79
Max. T (°C)2324.8929.8234.9640.0744.0738.3635.936.1533.327.6622.6
Mean T (°C)10.0111.6715.2919.912529.7827.9726.425.0420.4114.59.98
Min. T (°C)1.732.544.367.5212.0417.4319.9919.2116.9211.545.972.13
Empalme28° N110.85° WS. Albedo0.160.160.160.160.170.180.170.170.170.160.150.15
T. index (KT)0.620.630.680.680.680.650.580.570.610.630.640.62
NASA SSE3.844.686.146.997.417.226.445.895.684.964.163.575.58
Model3.824.676.046.967.517.356.475.995.74.934.123.595.60
Max. T (°C)27.8529.633.7437.9642.0145.4541.3339.3139.6237.2831.9827.42
Mean T (°C)16.117.4820.4424.1728.2332.4832.0130.9629.9625.8820.3816.15
Min. T (°C)8.899.4911.0913.6717.2922.2725.3825.2323.6118.4613.189.36
Etchojoa26.9° N109.72° WS. Albedo0.070.070.070.080.090.090.10.090.080.070.070.07
T. index (KT)0.650.660.70.710.730.70.640.620.640.660.660.63
NASA SSE4.245.086.517.398.087.947.186.596.155.464.523.896.09
Model4.124.996.297.298.057.897.116.536.035.254.363.775.97
Max. T (°C)24.8325.7728.4131.534.8437.2835.1834.5234.5933.8529.7825.61
Mean T (°C)18.218.7820.6823.5226.9330.1830.4230.3929.9127.6923.2719.32
Min. T (°C)14.0814.2515.3817.7320.9724.9127.1827.6226.8523.7919.2115.45
Fronteras30.91° N109.6° WS. Albedo0.150.150.150.150.160.160.150.150.140.140.140.14
T. index (KT)0.640.640.70.720.720.680.570.580.620.650.660.63
NASA SSE3.744.556.177.448.017.816.456.055.7854.123.455.71
Model3.644.486.017.287.977.766.396.065.654.853.953.355.62
Max. T (°C)20.0322.6627.6833.1638.8644.340.0737.7736.6731.7224.8719.63
Mean T (°C)6.999.0112.817.6423.2928.8628.0326.3924.0418.111.446.86
Min. T (°C)−1.19−0.321.6959.9715.6818.9318.114.988.782.77−0.95
General Plutarco Elías Calles32.01° N113.16° WS. Albedo0.210.210.210.220.230.230.230.230.230.210.210.21
T. index (KT)0.60.610.640.660.680.680.620.60.60.610.620.6
NASA SSE3.294.125.416.597.537.736.976.165.334.413.593.035.35
Model3.314.185.436.637.537.796.976.255.424.463.613.085.39
Max. T (°C)23.8326.4731.7836.943.0148.2548.6147.345.2838.2329.6223
Mean T (°C)11.3613.4217.4421.9827.6132.6835.1634.3531.2223.9816.3710.83
Min. T (°C)3.924.867.2710.3915.0319.7824.8524.9221.5414.658.333.67
Granados29.86° N109.33° WS. Albedo0.140.140.140.150.150.160.150.140.140.140.130.13
T. index (KT)0.640.650.70.720.730.690.60.610.640.660.660.64
NASA SSE3.884.726.327.528.127.966.786.365.965.154.243.65.88
Model3.754.656.097.318.077.856.726.395.895.014.063.515.78
Max. T (°C)2324.8929.8234.9640.0744.0738.3635.936.1533.327.6622.6
Mean T (°C)10.0111.6715.2919.912529.7827.9726.425.0420.4114.59.98
Min. T (°C)1.732.544.367.5212.0417.4319.9919.2116.9211.545.972.13
Guaymas28° N111.07° WS. Albedo0.120.120.110.130.140.140.140.140.140.120.120.12
T. index (KT)0.640.650.690.710.710.690.630.610.640.650.650.63
NASA SSE3.914.796.287.337.887.886.996.3765.124.223.635.87
Model3.944.826.137.267.847.817.026.415.985.094.193.655.85
Max. T (°C)24.4925.728.7732.136.0940.1438.6437.4137.3934.629.1724.67
Mean T (°C)16.2517.2419.4922.5826.530.7431.6131.1830.2926.3920.9616.79
Min. T (°C)11.0511.5112.8715.2418.823.3426.6926.9925.6821.0315.7711.89
Hermosillo29.17° N111.03° WS. Albedo0.150.160.160.160.170.180.170.160.160.150.150.15
T. index (KT)0.620.640.680.690.690.660.590.590.620.640.640.62
NASA SSE3.724.66.097.097.547.476.66.045.784.944.073.475.62
Model3.714.645.967.037.637.56.596.195.744.914.013.485.62
Max. T (°C)28.0430.1134.8439.5544.1848.4843.9341.4441.4938.1432.2327.38
Mean T (°C)14.8516.4319.924.0328.5233.3132.6131.1929.7424.8418.9814.61
Min. T (°C)7.047.629.4512.1315.8521.0124.392421.8416.21117.23
Huachinera30.21° N108.97° WS. Albedo0.140.140.140.140.150.160.150.140.140.140.130.14
T. index (KT)0.630.640.690.710.70.660.550.560.620.650.660.63
NASA SSE3.74.516.127.37.827.526.215.835.74.954.073.425.6
Model3.664.555.987.27.747.526.165.865.694.94.033.425.56
Max. T (°C)17.5419.5123.7328.3233.4138.0734.6633.232.1127.9522.0617.3
Mean T (°C)6.067.7210.9315.2220.2225.2124.5223.3621.3316.2110.316.08
Min. T (°C)−1.42−0.650.863.938.4213.6416.4715.8813.297.832.32−1.1
Huásabas29.91° N109.32° WS. Albedo0.140.140.140.150.150.160.150.140.140.140.130.13
T. index (KT)0.640.650.70.720.730.690.60.610.640.660.660.64
NASA SSE3.884.726.327.528.127.966.786.365.965.154.243.65.88
Model3.754.656.097.318.077.856.726.395.8954.063.515.77
Max. T (°C)2324.8929.8234.9640.0744.0738.3635.936.1533.327.6622.6
Mean T (°C)10.0111.6715.2919.912529.7827.9726.425.0420.4114.59.98
Min. T (°C)1.732.544.367.5212.0417.4319.9919.2116.9211.545.972.13
Huatabampo26.86° N109.71° WS. Albedo0.070.070.070.080.090.090.10.090.080.070.070.07
T. index (KT)0.650.660.70.710.730.70.640.620.640.660.660.63
NASA SSE4.245.086.517.398.087.947.186.596.155.464.523.896.09
Model4.124.996.297.298.047.897.116.536.035.264.373.775.97
Max. T (°C)24.8325.7728.4131.534.8437.2835.1834.5234.5933.8529.7825.61
Mean T (°C)18.218.7820.6823.5226.9330.1830.4230.3929.9127.6923.2719.32
Min. T (°C)14.0814.2515.3817.7320.9724.9127.1827.6226.8523.7919.2115.45
Huépac29.92° N110.22° WS. Albedo0.170.170.170.170.180.190.180.170.170.160.160.16
T. index (KT)0.630.640.690.710.70.680.60.590.620.650.650.63
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
Model3.694.5767.217.747.746.726.185.74.9343.455.66
Max. T (°C)25.0527.1232.2837.5742.4746.5540.6437.8538.3835.7729.7824.57
Mean T (°C)11.613.2417.0121.6526.5731.3929.7528.0926.7122.116.0711.44
Min. T (°C)3.444.065.999.0513.2718.6521.5120.8518.4813.017.583.69
Ímuris30.82° N110.9° WS. Albedo0.160.160.160.160.170.180.170.160.160.150.150.15
T. index (KT)0.630.640.690.720.720.680.580.580.630.650.660.63
NASA SSE3.714.586.157.397.977.856.546.065.84.994.093.465.72
Model3.64.495.937.287.977.766.56.065.754.853.963.365.63
Max. T (°C)21.6423.9228.7933.9539.4344.4940.6238.5937.7232.8326.3121.19
Mean T (°C)8.6510.3313.9618.4523.8629.2128.6427.2525.0619.2112.928.45
Min. T (°C)0.761.323.216.1610.7316.1519.6619.1216.1610.14.620.91
La Colorada28.81° N110.58° WS. Albedo0.160.160.160.160.170.180.170.170.170.160.150.15
T. index (KT)0.620.630.680.680.680.650.580.570.610.630.640.62
NASA SSE3.844.686.146.997.417.226.445.895.684.964.163.575.58
Model3.744.65.996.947.517.376.485.985.664.874.043.515.56
Max. T (°C)28.5930.6435.5140.2844.7548.4843.4140.9241.1438.3232.6927.89
Mean T (°C)15.1116.7620.2624.4328.8333.432.4631.0229.725.0619.2814.88
Min. T (°C)7.027.699.5212.2416.0221.2824.5124.0722.0416.4311.17.25
Magdalena30.62° N111.06° WS. Albedo0.170.180.170.180.190.20.190.180.180.170.170.17
T. index (KT)0.630.640.680.70.690.660.580.580.620.640.640.63
NASA SSE3.654.546.017.27.557.486.495.935.664.883.993.45.57
Model3.624.515.867.087.637.536.56.065.674.793.863.385.54
Max. T (°C)24.0126.0431.0136.1441.7446.6142.6340.4739.835.0128.6523.49
Mean T (°C)10.7712.3216.0120.4325.8130.9730.3628.9726.9121.1815.0710.52
Min. T (°C)2.893.395.368.2612.7217.8721.2920.8118.0312.126.872.98
Mazatán29.01° N110.15° WS. Albedo0.170.170.170.170.180.190.180.170.170.160.160.16
T. index (KT)0.630.640.690.710.70.680.60.590.620.650.650.63
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
Model3.784.666.067.247.747.726.76.195.7554.093.555.71
Max. T (°C)26.928.9634.1339.2644.0547.7341.6238.9139.4637.1431.4226.36
Mean T (°C)13.3214.9918.7523.2628.0232.6230.9129.328.0223.617.7213.16
Min. T (°C)5.135.87.7410.6714.7720.122.8922.2520.0714.679.315.4
Moctezuma30.09° N110.16° WS. Albedo0.160.160.160.160.170.180.170.160.160.150.150.15
T. index (KT)0.630.640.690.720.720.680.580.580.630.650.660.63
NASA SSE3.714.586.157.397.977.856.546.065.84.994.093.465.72
Model3.674.565.997.37.967.756.56.075.794.914.043.445.67
Max. T (°C)22.7224.9730.0135.3740.5945.2940.2537.7137.5533.7727.4522.28
Mean T (°C)9.4811.2114.9419.6224.8630.0328.8227.2225.4320.1313.889.32
Min. T (°C)1.3523.947.0611.5617.0720.1919.5116.8110.975.361.58
Naco31.34° N109.98° WS. Albedo0.150.160.160.160.170.180.170.160.160.150.150.14
T. index (KT)0.620.630.680.710.70.670.550.560.620.650.650.62
NASA SSE3.534.365.927.27.737.586.185.775.614.823.893.275.49
Model3.494.385.817.167.757.666.175.845.634.813.853.255.48
Max. T (°C)18.4521.3126.4131.9237.8943.6639.9837.6536.1630.6623.418.01
Mean T (°C)5.677.8611.7416.6222.4628.2627.826.0923.4217.1110.155.5
Min. T (°C)−2.28−1.330.84.159.315.0518.4817.5714.27.811.58−2.11
Nácori Chico29.69° N109.02° WS. Albedo0.140.140.140.150.150.160.150.140.140.140.130.13
T. index (KT)0.640.650.70.720.730.690.60.610.640.660.660.64
NASA SSE3.884.726.327.528.127.966.786.365.965.154.243.65.88
Model3.774.676.17.328.077.856.716.395.95.024.083.535.78
Max. T (°C)2324.8929.8234.9640.0744.0738.3635.936.1533.327.6622.6
Mean T (°C)10.0111.6715.2919.912529.7827.9726.425.0420.4114.59.98
Min. T (°C)1.732.544.367.5212.0417.4319.9919.2116.9211.545.972.13
Nacozari de García30.38° N109.71° WS. Albedo0.150.150.150.150.160.160.150.150.140.140.140.14
T. index (KT)0.640.640.70.720.720.680.570.580.620.650.660.63
NASA SSE3.744.556.177.448.017.816.456.055.7854.123.455.71
Model3.74.536.057.297.967.756.396.075.684.894.013.45.64
Max. T (°C)22.1524.4329.4834.940.3245.1139.9837.4237.1133.1926.8821.7
Mean T (°C)8.8510.6714.4119.224.5829.8128.4726.8124.9619.6113.38.72
Min. T (°C)0.51.293.246.511.2216.819.7919.0116.2610.384.60.79
Navojoa27.11° N109.51° WS. Albedo0.120.120.120.120.130.130.140.130.130.110.110.12
T. index (KT)0.640.650.70.70.70.670.590.580.620.640.660.63
NASA SSE4.074.946.397.267.87.566.596.095.885.184.383.775.83
Model4.044.96.287.197.727.566.566.15.845.084.343.745.78
Max. T (°C)28.8830.634.6738.6942.4144.9539.7637.7138.1337.8333.4128.91
Mean T (°C)17.5218.7421.5625.1529.0132.7831.5630.6129.8927.1222.1618.02
Min. T (°C)10.491112.5315.1818.7223.5625.8525.6824.3820.2415.1711.37
Nogales31.33° N111° WS. Albedo0.160.160.160.170.180.180.170.160.160.150.160.15
T. index (KT)0.620.630.680.710.710.680.550.560.620.640.650.62
NASA SSE3.554.45.977.257.97.776.285.785.644.793.933.295.55
Model3.494.385.817.167.867.776.175.855.634.743.853.255.50
Max. T (°C)19.8422.2327.132.2438.0743.5640.1838.1336.8631.4124.5619.36
Mean T (°C)7.028.7912.4216.9422.5728.1727.9726.5224.0217.811.296.79
Min. T (°C)−0.79−0.181.734.729.5114.9918.7518.1314.928.642.98−0.68
Ónavas28.46° N109.53° WS. Albedo0.140.140.140.150.150.160.160.150.140.130.130.13
T. index (KT)0.640.650.70.720.720.680.590.590.620.650.660.63
NASA SSE44.836.387.458.017.736.656.25.95.184.323.75.86
Model3.94.786.197.367.957.716.586.25.775.054.213.65.78
Max. T (°C)28.0830.0234.8939.7843.9946.8339.6136.5337.4137.0132.3527.69
Mean T (°C)14.9216.4819.9224.1428.4532.6630.2628.627.7724.4819.2114.94
Min. T (°C)6.727.399.1311.9815.7921.0723.2322.5520.8416.0810.947.16
Opodepe30.4° N111.56° WS. Albedo0.170.180.170.180.190.20.190.180.180.170.170.17
T. index (KT)0.630.640.680.70.690.660.580.580.620.640.640.63
NASA SSE3.654.546.017.27.557.486.495.935.664.883.993.45.57
Model3.74.65.887.097.637.526.56.075.684.813.893.45.56
Max. T (°C)25.5327.6132.637.4842.8547.5143.9341.9341.3536.8530.3224.95
Mean T (°C)12.313.8917.5421.7426.7831.7331.6730.5128.5322.9116.7612.04
Min. T (°C)4.395.0279.7513.8818.8522.6622.419.813.848.544.49
Oquitoa30.74° N111.74° WS. Albedo0.170.180.170.180.190.20.190.180.180.170.170.17
T. index (KT)0.630.640.680.70.690.660.580.580.620.640.640.63
NASA SSE3.654.546.017.27.557.486.495.935.664.883.993.45.57
Model3.64.55.857.087.637.536.56.065.664.783.853.375.53
Max. T (°C)24.9427.0131.9736.9742.6447.4944.4342.4541.5536.5529.8524.4
Mean T (°C)11.8213.4117.0721.3626.6531.731.9230.7228.5322.6316.3411.55
Min. T (°C)3.984.586.589.4113.7618.7322.722.3919.613.558.174.04
Pitiquito30.72° N112.11° WS. Albedo0.10.110.110.120.120.120.120.110.110.10.10.1
T. index (KT)0.620.630.670.690.690.670.610.610.630.630.620.6
NASA SSE3.594.475.937.087.617.686.886.315.784.793.873.295.61
Model3.554.435.776.987.637.656.846.385.754.713.733.215.55
Max. T (°C)25.2827.1931.9136.5441.8846.6845.0743.4942.4537.2930.424.79
Mean T (°C)12.7814.3117.7821.8126.7931.6732.9332.0429.8223.917.5412.62
Min. T (°C)5.385.997.9610.714.8119.5223.9923.9621.3115.279.785.51
Puerto Peñasco31.38° N113.54° WS. Albedo0.190.190.190.20.220.220.210.20.20.190.190.19
T. index (KT)0.570.580.620.640.670.660.60.590.580.580.570.57
NASA SSE3.153.935.326.47.257.46.625.925.134.183.342.925.13
Model3.24.035.36.457.427.556.736.165.274.293.372.995.23
Max. T (°C)19.2320.0522.6125.528.932.0834.6235.434.7830.6325.1820.38
Mean T (°C)15.2815.8717.9620.6823.9927.1430.4231.4230.3926.1721.0416.6
Min. T (°C)12.8713.1414.7317.0820.1823.2427.3628.6827.5123.3218.5114.31
Quiriego27.53° N109.26° WS. Albedo0.120.120.120.120.130.130.140.130.130.110.110.12
T. index (KT)0.640.650.70.70.70.670.590.580.620.640.660.63
NASA SSE4.074.946.397.267.87.566.596.095.885.184.383.775.83
Model3.994.866.257.177.727.576.576.15.825.044.33.75.76
Max. T (°C)27.1628.9733.3137.8641.9344.1937.4534.8735.3135.1631.126.92
Mean T (°C)14.5215.9619.0523.0927.2631.1528.8127.4326.5423.5318.6114.67
Min. T (°C)6.527.138.6411.4515.1820.3722.321.6920.0815.6510.667.06
Rayón29.72° N110.58° WS. Albedo0.170.170.170.170.180.190.180.170.170.160.160.16
T. index (KT)0.630.640.690.710.70.680.60.590.620.650.650.63
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
Model3.714.596.017.217.747.746.716.185.714.954.023.475.67
Max. T (°C)25.9328.0233.0338.0842.846.9641.8139.4139.5936.5330.5425.43
Mean T (°C)12.6514.2317.9122.3427.0631.8830.6329.1727.6722.9317.0112.46
Min. T (°C)4.765.297.1910.0714.0319.2822.2721.7519.43148.84.97
Rosario28.42° N109.72° WS. Albedo0.140.140.140.150.150.160.160.150.140.130.130.13
T. index (KT)0.640.650.70.720.720.680.590.590.620.650.660.63
NASA SSE44.836.387.458.017.736.656.25.95.184.323.75.86
Model3.94.786.197.357.957.76.586.25.775.054.213.615.77
Max. T (°C)28.0830.0234.8939.7843.9946.8339.6136.5337.4137.0132.3527.69
Mean T (°C)14.9216.4819.9224.1428.4532.6630.2628.627.7724.4819.2114.94
Min. T (°C)6.727.399.1311.9815.7921.0723.2322.5520.8416.0810.947.16
Sahuaripa29.07° N109.26° WS. Albedo0.140.140.140.150.150.160.150.140.140.140.130.13
T. index (KT)0.640.650.70.720.730.690.60.610.640.660.660.64
NASA SSE3.884.726.327.528.127.966.786.365.965.154.243.65.88
Model3.834.726.147.348.077.836.76.45.935.084.143.65.82
Max. T (°C)24.5226.3631.3636.4741.3244.7137.9835.2536.1634.4429.224.2
Mean T (°C)11.3512.9816.6221.1926.0830.5928.2526.6225.6521.5915.8911.39
Min. T (°C)3.033.815.648.7513.1318.4820.7619.9517.9512.847.393.49
San Felipe de Jesús29.87° N110.24° WS. Albedo0.170.170.170.170.180.190.180.170.170.160.160.16
T. index (KT)0.630.640.690.710.70.680.60.590.620.650.650.63
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
Model3.694.5867.217.747.746.726.185.74.9343.465.66
Max. T (°C)25.0527.1232.2837.5742.4746.5540.6437.8538.3835.7729.7824.57
Mean T (°C)11.613.2417.0121.6526.5731.3929.7528.0926.7122.116.0711.44
Min. T (°C)3.444.065.999.0513.2718.6521.5120.8518.4813.017.583.69
San Ignacio Río Muerto27.41° N110.25° WS. Albedo0.070.070.070.080.090.090.090.080.080.070.070.07
T. index (KT)0.650.660.710.710.730.710.650.630.650.670.660.63
NASA SSE4.144.996.57.448.128.167.366.696.255.444.443.796.11
Model4.074.956.357.288.058.027.236.636.15.294.313.716.00
Max. T (°C)25.8526.9329.9833.2236.6439.5637.2236.2536.5435.4730.8426.46
Mean T (°C)18.1518.8721.0223.9927.5431.0931.2931.0530.5627.9923.1819.12
Min. T (°C)13.4813.7415.0417.4120.7424.8727.3127.6626.7923.3518.5114.72
San Javier28.59° N109.74° WS. Albedo0.140.140.140.150.150.160.160.150.140.130.130.13
T. index (KT)0.640.650.70.720.720.680.590.590.620.650.660.63
NASA SSE44.836.387.458.017.736.656.25.95.184.323.75.86
Model3.884.776.187.47.957.716.596.195.775.044.193.595.77
Max. T (°C)26.8628.8133.8839.0143.6446.7839.5836.4337.3736.3131.326.43
Mean T (°C)13.5215.1518.7823.2327.932.2829.8828.1227.1923.4417.913.49
Min. T (°C)5.265.977.8110.7814.8920.2322.5221.7519.8514.799.515.67
San Luis Río Colorado32.49° N114.86° WS. Albedo0.20.20.20.210.230.230.220.220.220.20.20.2
T. index (KT)0.590.610.630.650.680.680.630.60.60.610.610.59
NASA SSE3.254.095.356.457.447.747.046.185.334.413.522.985.32
Model3.214.145.316.527.537.87.086.245.394.423.52.985.34
Max. T (°C)23.2526.6732.337.6243.2648.449.9549.3646.4338.3729.2122.36
Mean T (°C)11.3613.8818.0922.6727.9732.8236.135.7732.0924.316.4210.79
Min. T (°C)4.155.568.0311.2515.4719.7425.0925.6822.0515.078.683.88
San Miguel de Horcasitas29.49° N110.74° WS. Albedo0.170.170.170.170.180.190.180.170.170.160.160.16
T. index (KT)0.630.640.690.710.70.680.60.590.620.650.650.63
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
Model3.734.616.037.227.747.736.716.185.724.964.043.55.68
Max. T (°C)27.8629.9534.8939.7544.3348.2443.2740.9741.138.0332.2127.26
Mean T (°C)14.4216.0619.6923.9828.533.1732.1230.7129.2524.5418.714.2
Min. T (°C)6.487.088.9511.7215.5520.7723.8923.4121.2315.7710.546.69
San Pedro de la Cueva29.29° N109.76° WS. Albedo0.140.140.140.150.150.160.150.140.140.140.130.13
T. index (KT)0.640.650.70.720.730.690.60.610.640.660.660.64
NASA SSE3.884.726.327.528.127.966.786.365.965.154.243.65.88
Model3.814.76.137.338.077.846.716.45.925.064.123.575.81
Max. T (°C)26.0228.0133.2738.5643.4746.9840.0237.0337.9736.2230.6925.56
Mean T (°C)12.4414.1117.9322.5827.5132.0629.7828.0426.9922.8116.9512.36
Min. T (°C)4.034.786.769.8514.1819.5722.0321.2419.1313.838.354.39
Santa Ana30.55° N111.21° WS. Albedo0.170.180.170.180.190.20.190.180.180.170.170.17
T. index (KT)0.630.640.680.70.690.660.580.580.620.640.640.63
NASA SSE3.654.546.017.27.557.486.495.935.664.883.993.45.57
Model3.624.525.877.097.637.536.56.075.674.83.873.395.55
Max. T (°C)24.0126.0431.0136.1441.7446.6142.6340.4739.835.0128.6523.49
Mean T (°C)10.7712.3216.0120.4325.8130.9730.3628.9726.9121.1815.0710.52
Min. T (°C)2.893.395.368.2612.7217.8721.2920.8118.0312.126.872.98
Santa Cruz31.232° N110.597° WS. Albedo0.160.160.160.170.180.180.170.160.160.150.160.15
T. index (KT)0.620.630.680.710.710.680.550.560.620.640.650.62
NASA SSE3.554.45.977.257.97.776.285.785.644.793.933.295.55
Model3.54.395.827.167.867.776.175.855.644.743.863.265.50
Max. T (°C)19.8422.2327.132.2438.0743.5640.1838.1336.8631.4124.5619.36
Mean T (°C)7.028.7912.4216.9422.5728.1727.9726.5224.0217.811.296.79
Min. T (°C)−0.79−0.181.734.729.5114.9918.7518.1314.928.642.98−0.68
Sáric31.11° N111.39° WS. Albedo0.160.180.170.180.190.20.190.180.170.170.170.16
T. index (KT)0.620.640.680.70.70.680.580.570.620.640.650.62
NASA SSE3.544.435.927.227.747.766.565.95.664.793.93.275.56
Model3.514.475.837.077.757.776.515.955.644.753.873.285.53
Max. T (°C)2224.1529.0334.1439.9545.1641.8839.7638.6833.2826.6721.49
Mean T (°C)9.110.714.3218.7524.329.6729.4928.0725.7319.6213.378.83
Min. T (°C)1.341.873.816.7311.3516.5920.2819.7316.7310.65.261.42
Soyopa28.77° N109.64° WS. Albedo0.140.140.140.150.150.160.160.150.140.130.130.13
T. index (KT)0.640.650.70.720.720.680.590.590.620.650.660.63
NASA SSE44.836.387.458.017.736.656.25.95.184.323.75.86
Model3.874.756.177.347.957.716.596.195.765.024.183.575.76
Max. T (°C)26.8628.8133.8839.0143.6446.7839.5836.4337.3736.3131.326.43
Mean T (°C)13.5215.1518.7823.2327.932.2829.8828.1227.1923.4417.913.49
Min. T (°C)5.265.977.8110.7814.8920.2322.5221.7519.8514.799.515.67
Suaqui Grande28.4° N109.93° WS. Albedo0.140.140.140.150.150.160.160.150.140.130.130.13
T. index (KT)0.640.650.70.720.720.680.590.590.620.650.660.63
NASA SSE44.836.387.458.017.736.656.25.95.184.323.75.86
Model3.94.786.197.367.957.76.586.25.785.054.213.615.78
Max. T (°C)28.0830.0234.8939.7843.9946.8339.6136.5337.4137.0132.3527.69
Mean T (°C)14.9216.4819.9224.1428.4532.6630.2628.627.7724.4819.2114.94
Min. T (°C)6.727.399.1311.9815.7921.0723.2322.5520.8416.0810.947.16
Tepache29.63° N109.93° WS. Albedo0.140.140.140.150.150.160.150.140.140.140.130.13
T. index (KT)0.640.650.70.720.730.690.60.610.640.660.660.64
NASA SSE3.884.726.327.528.127.966.786.365.965.154.243.65.88
Model3.784.676.117.328.077.856.716.395.95.034.093.545.79
Max. T (°C)24.4126.4431.6637.0542.1846.2139.9137.0137.635.0729.1623.94
Mean T (°C)10.912.5916.3921.1426.2631.0829.1427.426.0921.4815.4310.78
Min. T (°C)2.473.235.198.3812.8618.320.9820.1917.8212.336.722.79
Trincheras30.41° N111.57° WS. Albedo0.170.180.170.180.190.20.190.180.180.170.170.17
T. index (KT)0.630.640.680.70.690.660.580.580.620.640.640.63
NASA SSE3.654.546.017.27.557.486.495.935.664.883.993.45.57
Model3.644.535.887.097.637.526.56.075.684.813.893.45.55
Max. T (°C)25.5327.6132.637.4842.8547.5143.9341.9341.3536.8530.3224.95
Mean T (°C)12.313.8917.5421.7426.7831.7331.6730.5128.5322.9116.7612.04
Min. T (°C)4.395.0279.7513.8818.8522.6622.419.813.848.544.49
Tubutama30.89° N111.48° WS. Albedo0.170.180.170.180.190.20.190.180.180.170.170.17
T. index (KT)0.630.640.680.70.690.660.580.580.620.640.640.63
NASA SSE3.654.546.017.27.557.486.495.935.664.883.993.45.57
Model3.594.495.847.077.647.546.516.065.654.773.843.355.53
Max. T (°C)24.0126.0431.0136.1441.7446.6142.6340.4739.835.0128.6523.49
Mean T (°C)10.7712.3216.0120.4325.8130.9730.3628.9726.9121.1815.0710.52
Min. T (°C)2.893.395.368.2612.7217.8721.2920.8118.0312.126.872.98
Ures29.43° N110.41° WS. Albedo0.170.170.170.170.180.190.180.170.170.160.160.16
T. index (KT)0.630.640.690.710.70.680.60.590.620.650.650.63
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
Model3.744.626.037.227.747.736.716.185.734.974.053.55.69
Max. T (°C)26.928.9634.1339.2644.0547.7341.6238.9139.4637.1431.4226.36
Mean T (°C)13.3214.9918.7523.2628.0232.6230.9129.328.0223.617.7213.16
Min. T (°C)5.135.87.7410.6714.7720.122.8922.2520.0714.679.315.4
Villa Hidalgo30.17° N109.35° WS. Albedo0.150.150.150.150.160.160.150.150.140.140.140.14
T. index (KT)0.640.640.70.720.720.680.570.580.620.650.660.63
NASA SSE3.744.556.177.448.017.816.456.055.7854.123.455.71
Model3.724.556.077.37.967.756.386.075.694.914.033.435.66
Max. T (°C)2123.1928.0533.2438.6643.4338.8236.6636.0231.8925.720.6
Mean T (°C)8.139.9413.5718.2523.628.7727.5626.0524.1318.7912.588.06
Min. T (°C)−0.110.772.615.8810.6416.131918.2415.539.7740.24
Villa Pesqueira29.5° N110.08° WS. Albedo0.170.170.170.170.180.190.180.170.170.160.160.16
T. index (KT)0.630.640.690.710.70.680.60.590.620.650.650.63
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
Model3.734.616.037.227.747.736.716.185.724.964.043.55.68
Max. T (°C)25.0527.1232.2837.5742.4746.5540.6437.8538.3835.7729.7824.57
Mean T (°C)11.613.2417.0121.6526.5731.3929.7528.0926.7122.116.0711.44
Min. T (°C)3.444.065.999.0513.2718.6521.5120.8518.4813.017.583.69
Yécora28.38° N108.96° WS. Albedo0.130.130.120.130.140.140.140.130.130.120.120.12
T. index (KT)0.650.660.710.730.730.660.560.580.630.660.680.65
NASA SSE4.14.926.497.548.147.646.376.15.955.284.433.785.89
Model3.974.866.287.468.067.486.256.095.875.134.343.735.79
Max. T (°C)20.5421.9225.4729.5133.636.5532.0530.363028.4224.4720.56
Mean T (°C)9.2810.5313.1516.9721.1725.2923.7922.7421.5118.0213.279.57
Min. T (°C)1.792.33.396.159.9415.1817.3716.9115.1510.585.792.4

References

  1. Callejas-Cornejo, J.M.; Peña-Cruz, M.I.; Valentín-Coronado, L.M. Satellite Image Processing for Retrieving Historical Solar Irradiance Data Within the Mexican Territory. In Progress in Optomechatronic Technologies; Springer: Singapore, 2019; pp. 117–126. [Google Scholar] [CrossRef]
  2. National Aeronautics and Space Administration. NASA POWER—Prediction of Worldwide Energy Resources. 2019. Available online: https://power.larc.nasa.gov/ (accessed on 25 April 2020).
  3. Energy Sector Management Assistance Program. SOLARGIS, WB, IFC, Global Solar Atlas. 2019. Available online: https://globalsolaratlas.info/?m=sg:ghi (accessed on 25 April 2020).
  4. Guzmán, Á.L.; Lemus, R.A.C.; Campas, M.E.C.; Quintero, J.M.; Maldonado, J.A.; Roe, S.M.; Anderson, R. Greenhouse Gas Emissions in Sonora and Reference Case Projections 1990–2020; CEDES: Hermosillo, Sonora, Mexico; Center for Climate Strategies: Harriburg, PA, USA, 2010; Available online: www.climatestrategies.us/library/library/download/927 (accessed on 12 November 2020).
  5. Gonzalez-Cabrera, A.E.; Riveros-Rosas, D.; Valdes-Barrón, M.G.; Bonifaz-Alfonzo, R.; Velasco-Herrera, V.M.; Estevez-Perez, H.R.; Carabali, G. New Reference Solarimetric Network for Mexico. In AIP Conference Proceedings; AIP Publishing LLC: Melville, NY, USA, 2018. [Google Scholar] [CrossRef]
  6. Valdes-Barrón, M.; Riveros-Rosas, D.; Arancibia-Bulnes, C.A.; Bonifaz, R. The Solar Resource Assessment in Mexico: State of the Art. Energy Procedia 2014, 57, 1299–1308. [Google Scholar] [CrossRef] [Green Version]
  7. Mexican National Water Commission. Estaciones Meteorológicas Automáticas (EMA’s). 2020. Available online: https://smn.conagua.gob.mx/es/observando-el-tiempo/estaciones-meteorologicas-automaticas-ema-s (accessed on 25 April 2020).
  8. Riveros-Rosas, D.; Arancibia-Bulnes, C.A.; Bonifaz, R.; Medina, M.A.; Peón, R.; Valdes, M. Analysis of a Solarimetric Database for Mexico and Comparison with the CSR Model. Renew. Energy 2015, 75, 21–29. [Google Scholar] [CrossRef]
  9. Besharat, F.; Dehghan, A.A.; Faghih, A.R. Empirical Models for Estimating Global Solar Radiation: A Review and Case Study. Renew. Sustain. Energy Rev. 2013, 21, 798–821. [Google Scholar] [CrossRef]
  10. Google. Google Earth, (n.d.). Available online: https://earth.google.com/web/ (accessed on 14 March 2020).
  11. Kalogirou, S.A. Solar Energy Engineering: Processes and Systems, 2nd ed.; Academic Press: Waltham, MA, USA, 2014; ISBN 9780123972705. [Google Scholar]
  12. Ehnberg, J.S.G.; Bollen, M.H.J. Simulation of Global Solar Radiation Based on Cloud Observations. Sol. Energy 2005, 78, 157–162. [Google Scholar] [CrossRef]
  13. Soulouknga, M.H.; Falama, R.Z.; Ajayi, O.O.; Doka, S.Y.; Kofane, T.C. Determination of a Suitable Solar Radiation Model for the Sites of Chad. Energy Power Eng. 2017, 9, 703–722. [Google Scholar] [CrossRef] [Green Version]
  14. Katiyar, A.K.; Pandey, C.K. Solar radiation models. Part 1: A review. Hindawi J. Renew. Energy 2013, 2013, 1–11. [Google Scholar] [CrossRef] [Green Version]
  15. Stanciu, C.; Stanciu, D. Optimum Tilt Angle for Flat Plate Collectors all Over the World—A Declination Dependence Formula and Comparisons of Three Solar Radiation Models. Energy Convers. Manag. 2014, 81, 133–143. [Google Scholar] [CrossRef]
  16. Maleki, S.A.M.; Hizam, H.; Gomes, C. Estimation of Hourly, Daily and Monthly Global Solar Radiation on Inclined Surfaces: Models re-visited. Energies 2017, 10, 134. [Google Scholar] [CrossRef] [Green Version]
  17. Li, H.; Ma, W.; Lian, Y.; Wang, X. Estimating Daily Global Solar Radiation by Day of Year in China. Appl. Energy 2010, 87, 3011–3017. [Google Scholar] [CrossRef]
  18. Mirzabe, A.H.; Hajiahmad, A.; Keyhani, A.; Nejad, B.F. Developing New Mathematical Model Formats for Estimating Daily and Monthly Solar Radiation: Case Study: Iran. Int. J. Sustain. Energy 2020, 39, 896–926. [Google Scholar] [CrossRef]
  19. Tsvetkov, N.A.; Krivoshein, U.O.; Tolstykh, A.V.Y.; Khutornoi, A.N.; Boldyryev, S. The Calculation of Solar Energy Used by Hot Water Systems in Permafrost Region: An Experimental Case Study for Yakutia. Energy 2020, 210, 118577. [Google Scholar] [CrossRef]
  20. Obukhov, S.; Plontnikov, I.A.; Masolov, V.G. Mathematical Model of Solar Radiation Based on Climatological Data from NASA SSE. In IOP Conference Series: Materials Science and Engineering; IOP Publishing: Moscow, Russia, 2018; Volume 363. [Google Scholar] [CrossRef] [Green Version]
  21. Stackhouse, P.W.; Chandler, W.S.; Zhang, T.; Westberg, D.; Barnett, A.J.; Hoell, J.M. Surface Meteorology and Solar Energy (SSE) Release 6.0 Methodology Version 3.2.0. 2016. Available online: https://power.larc.nasa.gov/documents/SSE_Methodology.pdf (accessed on 8 May 2020).
  22. Klein, S.A. Calculation of Monthly Average Insolation on Tilted Surfaces. Sol. Energy 1976, 19, 325–329. [Google Scholar] [CrossRef] [Green Version]
  23. Obukhov, S.G.; Plotnikov, I.A. Mathematical Model of Solar Radiation Incident on an Arbitrarily Oriented Surface for any Region in Russia. Int. Sci. J. Altern. Energy Ecol. 2017, 16–18, 1–14. [Google Scholar] [CrossRef]
  24. Duffie, J.A.; Beckman, W.A. Solar Engineering of Thermal Processes, 3rd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2006; ISBN 978-0471698678. [Google Scholar]
  25. Munner, T. Solar Radiation and Daylight Models, 2nd ed.; Elsevier: Edinburgh, UK, 2004; ISBN 0750659742. [Google Scholar]
  26. Cooper, P.I. The Absorption of Radiation in Solar Stills. Sol. Energy 1969, 12, 333–346. [Google Scholar] [CrossRef]
  27. Sayigh, A.A.M. Solar Energy Engineering, 2nd ed.; Elsevier: Riyadh, Saudi Arabia, 2013; Available online: http://proquest.safaribooksonline.com/9780123972705 (accessed on 21 July 2020).
  28. Messenger, R.A.; Ventre, J. Photovoltaic Systems Engineering, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]
  29. Liu, B.Y.H.; Jordan, R.C. The Interrelationship and Characteristic Distribution of Direct, Diffuse and Total Solar Radiation. Sol. Energy 1960, 4, 1–19. [Google Scholar] [CrossRef]
  30. Collares-Pereira, M.; Rabl, A. The Average Distribution of Solar Radiation-Correlations Between Diffuse and Hemispherical and Between Daily and Hourly Insolation Values. Sol. Energy 1979, 22, 155–164. [Google Scholar] [CrossRef]
  31. Obukhov, S.G.; Plotnikov, I.A. Selection of Parameters and Analysis of Efficiency of Application of Sun Tracking Systems. Bull. Tomsk Polytech. Univ. Geo Assets Eng. 2018, 329, 95–106. [Google Scholar]
  32. Despotovic, M.; Nedic, V.; Despotovic, D.; Cvetanovic, S. Review and Statistical Analysis of Different Global Solar Radiation Sunshine Models. Renew. Sustain. Energy Rev. 2015, 52, 1869–1880. [Google Scholar] [CrossRef]
  33. Nottrott, A.; Kleissl, J. Validation of the NSRDB-SUNY Global Horizontal Irradiance in California. Sol. Energy 2010, 84, 1816–1827. [Google Scholar] [CrossRef]
  34. Tamhane, A.C. Statistical Analysis of Designed Experiments: Theory and Applications; John Wiley & Sons, Inc.: Hoboken, NY, USA, 2009. [Google Scholar] [CrossRef]
  35. Rahman, M.M.; Hasanuzzaman, M.; Rahim, N.A. Effects of Various Parameters on PV-Module Power and Efficiency. Energy Convers. Manag. 2015, 103, 348–358. [Google Scholar] [CrossRef]
  36. Nishioka, K.; Hatayama, T.; Uraoka, Y.; Fuyuki, T.; Hagihara, R.; Watanabe, M. Field-Test Analysis of PV System Output Characteristics Focusing on Module Temperature. Sol. Energy Mater. Sol. Cells 2003, 75, 665–671. [Google Scholar] [CrossRef]
  37. Luque, A.; Hegedus, S. Handbook of Photovoltaic Science and Engineering, 2nd ed.; John Wiley & Sons: Chennai, India, 2011; ISBN 9780470721698. [Google Scholar]
  38. Cuce, E.; Cuce, P.M.; Bali, T. An Experimental Analysis of Illumination Intensity and Temperature Dependency of Photovoltaic Cell Parameters. Appl. Energy 2013, 111, 374–382. [Google Scholar] [CrossRef]
Figure 1. Location of meteorological stations from the SMN in Mexico, highlighting the state of Sonora with a red outline [7,10].
Figure 1. Location of meteorological stations from the SMN in Mexico, highlighting the state of Sonora with a red outline [7,10].
Energies 13 06501 g001
Figure 2. Comparison between the SMN, NASA SSE, and the model for the station of Caborca.
Figure 2. Comparison between the SMN, NASA SSE, and the model for the station of Caborca.
Energies 13 06501 g002
Figure 3. Comparison between the SMN, NASA SSE, and the model for the station of El Pinacate.
Figure 3. Comparison between the SMN, NASA SSE, and the model for the station of El Pinacate.
Energies 13 06501 g003
Figure 4. Comparison between the SMN, NASA SSE, and the model for the station of Hermosillo.
Figure 4. Comparison between the SMN, NASA SSE, and the model for the station of Hermosillo.
Energies 13 06501 g004
Figure 5. Comparison between the SMN, NASA SSE, and the model for the station of Nogales.
Figure 5. Comparison between the SMN, NASA SSE, and the model for the station of Nogales.
Energies 13 06501 g005
Figure 6. Comparison between the SMN, NASA SSE, and the model for the station of Yécora.
Figure 6. Comparison between the SMN, NASA SSE, and the model for the station of Yécora.
Energies 13 06501 g006
Figure 7. GIS analysis of Sonora for January: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 7. GIS analysis of Sonora for January: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g007
Figure 8. GIS analysis of Sonora for February: (left) average solar total radiation (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 8. GIS analysis of Sonora for February: (left) average solar total radiation (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g008
Figure 9. GIS analysis of Sonora for March: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 9. GIS analysis of Sonora for March: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g009
Figure 10. GIS analysis of Sonora for April: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 10. GIS analysis of Sonora for April: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g010
Figure 11. GIS analysis of Sonora for May: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 11. GIS analysis of Sonora for May: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g011
Figure 12. GIS analysis of Sonora for June: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 12. GIS analysis of Sonora for June: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g012
Figure 13. GIS analysis of Sonora for July: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 13. GIS analysis of Sonora for July: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g013
Figure 14. GIS analysis of Sonora for August: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 14. GIS analysis of Sonora for August: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g014
Figure 15. GIS analysis of Sonora for September: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 15. GIS analysis of Sonora for September: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g015
Figure 16. GIS analysis of Sonora for October: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 16. GIS analysis of Sonora for October: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g016
Figure 17. GIS analysis of Sonora for November: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 17. GIS analysis of Sonora for November: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g017
Figure 18. GIS analysis of Sonora for December: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Figure 18. GIS analysis of Sonora for December: (left) average total solar irradiance (kWh/m2), (center) maximum average temperatures (°C), and (right) minimum average temperatures (°C).
Energies 13 06501 g018
Table 1. Selected meteorological stations in Sonora (location and altitude) [7].
Table 1. Selected meteorological stations in Sonora (location and altitude) [7].
StationLatitudeLongitudeAltitudeData%
Caborca30.775° N112.435° W188 m100
El Pinacate31.68° N113.305° W99 m93.3
Hermosillo29.078° N110.93° W211 m93.3
Nogales31.297° N110.914° W1269 m98.3
Yécora28.36° N108.93° W1531 m100
Table 2. Day numbers per month according to [22,23].
Table 2. Day numbers per month according to [22,23].
MonthAverage Day of the Month
Day of the MonthDay Number of the Year (N)
January1717
February1647
March1675
April15105
May15135
June11162
July17198
August16228
September15258
October15288
November14318
December10344
Table 3. Conditional table to determine the diffusion index ( K D ) [21].
Table 3. Conditional table to determine the diffusion index ( K D ) [21].
ConditionEquation
i f   φ < 45 ° K D = 0.96268 1.452 K T + 0.27365 K T 2 + 0.04279 K T 3 + 0.000246 ω s s + 0.001189 [ ( 90 φ ) + δ ]
e l s e   i f   ω s s > 150 ° K D = 0.6563 2.893 K T + 4.594 K T 2 3.23 K T 3 + 0.004 ω s s 0.0023 [ ( 90 φ ) + δ ]
e l s e   i f   ω s s > 125 ° K D = 1.6586 4.412 K T + 5.8 K T 2 3.1223 K T 3 + 0.000144 ω s s 0.000829 [ ( 90 φ ) + δ ]
e l s e   i f   ω s s > 100 ° K D = 0.3498 + 3.8035 K T 11.765 K T 2 + 9.1748 K T 3 + 0.001575 ω s s 0.002837 [ ( 90 φ ) + δ ]
e l s e   i f   ω s s > 81.4 ° K D = 1.6821 2.5866 K T + 2.373 K T 2 0.5294 K T 3 0.00277 ω s s 0.004233 [ ( 90 φ ) + δ ]
e l s e K D = 1.441 3.6839 K T + 6.4927 K T 2 4.147 K T 3 + 0.0008 ω s s 0.008175 [ ( 90 φ ) + δ ]
Table 4. Comparison between the NASA SSE database and the model.
Table 4. Comparison between the NASA SSE database and the model.
StationParameterMonthAnnual AverageStatistical Parameters
JanFebMarAprMayJunJulAugSeptOctNovDecMAEMBERMSEMPErR2t
CaborcaNASA SSE3.594.475.937.087.617.686.886.315.784.793.873.295.610.057−0.0570.082−1.2250.9990.9973.152
Model3.544.435.776.987.637.656.846.375.754.713.733.25.55
RPE−1.393−0.895−2.698−1.4120.263−0.391−0.5810.951−0.519−1.670−3.618−2.736
El PinacateNASA SSE3.153.935.326.47.257.46.625.925.134.183.342.925.130.0850.0850.1121.5430.9990.9953.839
Model3.1745.286.447.427.556.746.155.254.273.352.965.22
RPE0.6351.781−0.7520.6252.3452.0271.8133.8852.3392.1530.2991.370
HermosilloNASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.730.032−0.0320.061−0.6160.9990.9982.035
Model3.774.656.067.237.747.726.76.195.7454.083.545.70
RPE−0.789−0.215−2.100−1.0940.2590.1300.1490.814−1.205−1.186−2.1580.000
NogalesNASA SSE3.554.45.977.257.97.776.285.785.644.793.933.295.550.048−0.0480.073−0.9260.9990.9982.875
Model3.494.385.827.167.867.776.175.855.634.743.853.265.50
RPE−1.690−0.455−2.513−1.241−0.5060.000−1.7521.211−0.177−1.044−2.036−0.912
YécoraNASA SSE4.14.926.497.548.147.646.376.15.955.284.433.785.890.102−0.1020.114−1.7790.9990.9936.433
Model3.974.866.287.468.067.486.256.095.875.134.343.735.79
RPE−3.171−1.220−3.236−1.061−0.983−2.094−1.884−0.164−1.345−2.841−2.032−1.323
Table 5. Comparison between the SMN and the model.
Table 5. Comparison between the SMN and the model.
StationParameterMonthAnnual AverageStatistical Parameters
JanFebMarAprMayJunJulAugSeptOctNovDecMAEMBERMSEMPErR2t
CaborcaSMN3.74.656.097.588.047.887.087.136.035.194.013.345.890.343−0.3430.387−5.7800.9950.9436.358
Model3.544.435.776.987.637.656.846.375.754.713.733.25.55
RPE−4.324−4.731−5.255−7.916−5.100−2.919−3.390−10.659−4.643−9.249−6.983−4.192
El PinacateSMN3.674.686.177.578.218.27.637.386.195.1843.356.020.804−0.8040.838−13.6690.9940.76411.380
Model3.1745.286.447.427.556.746.155.254.273.352.965.22
RPE−13.624−14.530−14.425−14.927−9.622−7.927−11.664−16.667−15.186−17.568−16.250−11.642
HermosilloSMN3.694.55.927.187.87.546.916.665.614.973.783.245.650.0520.0520.2141.6830.9920.9800.825
Model3.774.656.067.237.747.726.76.195.7454.083.545.70
RPE2.1683.3332.3650.696−0.7692.387−3.039−7.0572.3170.6047.9379.259
NogalesSMN4.014.826.247.638.317.846.486.456.065.364.293.495.920.417−0.4170.442−7.5740.9950.9139.438
Model3.494.385.827.167.867.776.175.855.634.743.853.265.50
RPE−12.968−9.129−6.731−6.160−5.415−0.893−4.784−9.302−7.096−11.567−10.256−6.590
YécoraSMN4.184.916.17.558.317.455.995.925.675.144.173.685.760.0380.0380.1630.7580.9940.9860.784
Model3.974.866.287.468.067.486.256.095.875.134.343.735.79
RPE−5.024−1.0182.951−1.192−3.0080.4034.3412.8723.527−0.1954.0771.359
Table 6. Comparison between the SMN and the NASA SSE database.
Table 6. Comparison between the SMN and the NASA SSE database.
StationParameterMonthAnnual AverageStatistical Parameters
JanFebMarAprMayJunJulAugSeptOctNovDecMAEMBERMSEMPErR2t
CaborcaSMN3.74.656.097.588.047.887.087.136.035.194.013.345.890.287−0.2870.354−4.5930.9940.9524.586
NASA SSE3.594.475.937.087.617.686.886.315.784.793.873.295.61
RPE−2.973−3.871−2.627−6.596−5.348−2.538−2.825−11.501−4.146−7.707−3.491−1.497
El PinacateSMN3.674.686.177.578.218.27.637.386.195.1843.356.020.889−0.8890.930−14.9720.9920.70910.803
NASA SSE3.153.935.326.47.257.46.625.925.134.183.342.925.13
RPE−14.169−16.026−13.776−15.456−11.693−9.756−13.237−19.783−17.124−19.305−16.500−12.836
HermosilloSMN3.694.55.927.187.87.546.916.665.614.973.783.245.650.0830.0830.2542.3420.9900.9731.154
NASA SSE3.84.666.197.317.727.716.696.145.815.064.173.545.73
RPE2.9813.5564.5611.811−1.0262.255−3.184−7.8083.5651.81110.3179.259
NogalesSMN4.014.826.247.638.317.846.486.456.065.364.293.495.920.369−0.3690.402−6.7070.9950.9287.709
NASA SSE3.554.45.977.257.97.776.285.785.644.793.933.295.55
RPE−11.471−8.714−4.327−4.980−4.934−0.893−3.086−10.388−6.931−10.634−8.392−5.731
YécoraSMN4.184.916.17.558.317.455.995.925.675.144.173.685.760.1390.1390.2192.5880.9930.9752.727
NASA SSE4.14.926.497.548.147.646.376.15.955.284.433.785.89
RPE−1.9140.2046.393−0.132−2.0462.5506.3443.0414.9382.7246.2352.717
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Enríquez-Velásquez, E.A.; Benitez, V.H.; Obukhov, S.G.; Félix-Herrán, L.C.; Lozoya-Santos, J.d.-J. Estimation of Solar Resource Based on Meteorological and Geographical Data: Sonora State in Northwestern Territory of Mexico as Case Study. Energies 2020, 13, 6501. https://doi.org/10.3390/en13246501

AMA Style

Enríquez-Velásquez EA, Benitez VH, Obukhov SG, Félix-Herrán LC, Lozoya-Santos Jd-J. Estimation of Solar Resource Based on Meteorological and Geographical Data: Sonora State in Northwestern Territory of Mexico as Case Study. Energies. 2020; 13(24):6501. https://doi.org/10.3390/en13246501

Chicago/Turabian Style

Enríquez-Velásquez, Enrique A., Victor H. Benitez, Sergey G. Obukhov, Luis C. Félix-Herrán, and Jorge de-J. Lozoya-Santos. 2020. "Estimation of Solar Resource Based on Meteorological and Geographical Data: Sonora State in Northwestern Territory of Mexico as Case Study" Energies 13, no. 24: 6501. https://doi.org/10.3390/en13246501

APA Style

Enríquez-Velásquez, E. A., Benitez, V. H., Obukhov, S. G., Félix-Herrán, L. C., & Lozoya-Santos, J. d. -J. (2020). Estimation of Solar Resource Based on Meteorological and Geographical Data: Sonora State in Northwestern Territory of Mexico as Case Study. Energies, 13(24), 6501. https://doi.org/10.3390/en13246501

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