Next Article in Journal
Smart Model to Distinguish Crohn’s Disease from Ulcerative Colitis
Next Article in Special Issue
Risk-Constrained Optimal Chiller Loading Strategy Using Information Gap Decision Theory
Previous Article in Journal
Effects of Velocity Profiles on Measuring Accuracy of Transit-Time Ultrasonic Flowmeter
Previous Article in Special Issue
Robust Planning of Energy and Environment Systems through Introducing Traffic Sector with Cost Minimization and Emissions Abatement under Multiple Uncertainties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Meteorological Variables’ Influence on Electric Power Generation for Photovoltaic Systems Located at Different Geographical Zones in Mexico

by
Jose A. Ruz-Hernandez
1,
Yasuhiro Matsumoto
2,
Fernando Arellano-Valmaña
1,
Nun Pitalúa-Díaz
3,*,
Rafael Enrique Cabanillas-López
3,
José Humberto Abril-García
3,
Enrique J. Herrera-López
4 and
Enrique Fernando Velázquez-Contreras
3
1
Facultad de Ingeniería, Universidad Autónoma del Carmen, Calle 56 No. 4 Esq. Avenida Concordia Col. Benito Juárez, Cd. Del Carmen C.P. 24180, Campeche, Mexico
2
Departamento de Ingeniería Eléctrica, Centro de Investigación y de Estudios Avanzados del IPN, Av. Instituto Politécnico Nacional 2508, San Pedro Zacatenco, Ciudad de México C.P. 07360, CDMX, Mexico
3
Departamento de Ingeniería Industrial, Departamento de Ingeniería Química y Metalurgia, Departamento de Investigación en Polímeros y Materiales, Universidad de Sonora, Blvd. Luis Encinas y Rosales S/N, Col. Centro, Hermosillo C.P. 83000, Sonora, Mexico
4
Biotecnología Industrial, Sublínea Bioelectrónica, Centro de Investigación y Asistencia en Tecnología y Diseño del Estado de Jalisco A.C., Camino Arenero 1227, Col. El Bajío del Arenal, Zapopan C.P. 45019, Jalisco, Mexico
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(8), 1649; https://doi.org/10.3390/app9081649
Submission received: 1 March 2019 / Revised: 4 April 2019 / Accepted: 10 April 2019 / Published: 20 April 2019
(This article belongs to the Special Issue Sustainable Energy Systems Planning, Integration and Management)

Abstract

:
In this study, the relation among different meteorological variables and the electrical power from photovoltaic systems located at different selected places in Mexico were presented. The data was collected from on-site real-time measurements from Mexico City and the State of Sonora. The statistical estimation by the gradient descent method demonstrated that solar radiation, outdoor temperature, wind speed, and daylight hour influenced the electric power generation when it was compared with the real power of each photovoltaic system. According to our results, 97.63% of the estimation results matched the real data for Sonora and 99.66% the results matched for Mexico City, achieving overall errors less than 7% and 2%, respectively. The results showed an acceptable performance since a satisfactory estimation error was achieved for the estimation of photovoltaic power with a high determination coefficient R2.

Graphical Abstract

1. Introduction

Renewables energies represent a potential alternative in the transition towards a low-carbon society, where photovoltaic sources play a key role; however, consumers are also investors and a project is implemented only if economic conditions are verified [1].
Solar technologies are characterized depending on the way they capture, convert, and distribute sunlight such as photovoltaic (PV) systems and their corresponding requirements for an energy storage arrangement [2,3]. Consequently, these technologies feed power to the electric grid by using solar panels as generators [4,5,6]. In addition, concentrating solar power plants (CSP) use mirrors to focus the energy from the sun to drive traditional steam turbines or engines to create electricity; solar heating and cooling (SHC) systems which collect the thermal energy from the sun use this heat to provide hot water, space heating, and cooling for residential, commercial, and industrial applications. These technologies displace the need to use electricity or natural gas [7,8,9].
A photovoltaic system is composed of several components: the solar panel and the inverter for grid-connected systems and additionally energy storage for stand-alone systems. The fabrication of the solar panel involves diverse stages. The first step is to define the type of panel, where the most known is monocrystalline or polycrystalline and where the solar cells are manufactured using wafers made of silicon [10]. Nowadays, the PV systems are widely used to generate electricity due its accessible cost [5,6,11,12]. Some studies aimed at reducing this cost even more, e.g., in Reference [13], a design optimization model for the residential PV systems in South Korea was proposed, where the objective function to be minimized consisted of three costs, such as the monthly electric bill, the PV-related construction costs, and the PV-related maintenance cost.
The solar radiation (photons) is responsible for the photovoltaic effect; nonetheless, some weather factors have an effect on the amount of energy generated even with the optimum radiation. A cloudy day generates a shadow on the solar panel by the time it has to capture the photons so that the incidence of these particles will be less, achieving a minor electric power. It is well-known that clouds are water steam concentrated in the air; in other words, humidity and temperature work together, and such a relation is an example of how the meteorological variables impact the generation of electric power [14,15].
Solar energy, besides wind energy, is currently the most resourceful renewable source worldwide [16,17]. Its obtainment, unlike many others currently used, does not mean any harm to the environment, and its resources overcome, by far, everyone else [12,18,19,20]. Some countries such as Germany, Italy, Spain, United States of America, and China are ahead on solar energy research; meanwhile in Mexico, being a country rich in solar radiation, with a great territorial extension, and having some solar energy studies, solar energy still does not have the necessary research compared to countries leading solar technology [4,5,6,11,12,21].
Figure 1 depicts the behavior of the horizontal solar radiation on the world. Similarly, if Mexico is compared with Europe, it can be seen that the only country with a notorious radiation incidence is Spain, achieving a maximum value between 4.8 and 5.4 kWh/m2 [22]. Mexico exceeds Spain both in incidence territory and in radiation intensity, with an average between 5.6 and 6.2 kWh/m2, as can be seen in Figure 2, exceeding even China, which mostly contemplates values of 4.6 kWh/m2.
Table 1 shows the data compilation from 2010 describing the global-horizontal solar radiation for some locations in Mexico. The total irradiation provided by solar energy in the year 2018 was summarized in Reference [25]. Both reports agree that the Northwest of Mexico reached the maximum values, the states of Sonora, Baja California, Coahuila, and Chihuahua being the main producers/receivers.
Recently, diverse reports have studied the relations between some meteorological variables and the electric power generated in a photovoltaic system. In Reference [26], an analysis of the climatic factors and solar data from the Andes site was performed; nonetheless, all the data gathered for the study was averaged every 10 min and the power measurements were estimated. In [27], the effect of diverse meteorological variables (outdoor temperature, air pressure, humidity, wind speed, and solar radiation) on the generated energy are mentioned, although a representative model of the plant was not obtained and all the data was experimentally generated. In Reference [28], a statistical method was applied to forecast the energy generated by a solar plant, though a database of 30 plants from different locations is necessary and although it does not have real power measurements from the site and the computational load is heavy. In Reference [29], an artificial neural network (ANN) is used to obtain a model to forecast the photovoltaic energy; however, the solar radiation was the only meteorological variable analyzed.
Higher Education Institutions (HEI) currently play a major role in the generation of human capital and the associated impact on societal development; HEIs are ideal locations to focus the resources in terms of the deployment and experimentation of decarbonization technologies to demonstrate the best practice for a further replication within wider society [30].
Careful planning is required to manage the future electricity demand of PV systems due to its increasing potential demand in Mexico [31,32]. Therefore, it is vital to understand the influence of meteorological variables on energy consumption in which a better understanding of it can contribute to a more useful strategy in meeting the energy efficiency goal for the country.
According to the above and considering References [1,13,30], the aim of this work is to present a statistical analysis based on the gradient descent method that is easy to implement and has a low computational load to estimate the electric power generation from the meteorological data such as the solar radiation, outdoor temperature, wind speed, and daylight time collected from PV systems located in Mexico City and Sonora [33]. This is important because most of the current photovoltaic system deployments do not monitor these factors or employ them in adaptation and prediction tasks [34,35,36].
The proposed methodology achieves a satisfactory estimation of the PV power with a high determination coefficient and a fair error percentage value.
Table 1. The solar radiation on select places in Mexico (data in kWh/m2 per day) [37].
Table 1. The solar radiation on select places in Mexico (data in kWh/m2 per day) [37].
StateCityJanFebMarAprMayJunJulAugSepOctNovDecAvg
SonoraHermosillo4.04.65.46.68.38.56.96.66.76.04.73.96.0
SonoraGuaymas4.55.76.57.27.36.85.95.86.35.95.05.65.9
ChihuahuaChihuahua4.14.96.07.48.28.16.86.25.75.24.63.85.9
SLPSLP4.35.35.86.46.36.16.46.05.54.74.23.75.4
ZacatecasZacatecas4.95.76.67.57.86.26.25.95.44.84.84.15.8
GuanajuatoGuanajuato4.45.16.16.36.66.06.05.95.85.24.84.65.6
AguascalientesAguascalientes4.55.25.96.67.26.36.15.95.75.14.84.05.6
OaxacaSalina Cruz5.46.36.66.46.15.05.65.95.25.95.75.25.8
OaxacaOaxaca4.95.75.85.56.05.45.95.65.04.94.84.45.3
JaliscoColotlán4.65.76.57.58.26.65.85.65.85.34.94.15.9
JaliscoGuadalajara4.65.56.37.47.75.95.35.35.24.94.84.05.6
DurangoDurango4.45.46.57.07.56.865.65.75.14.83.95.7
Baja CaliforniaLa Paz4.45.56.06.66.56.66.36.25.95.84.94.25.7
Baja CaliforniaSan Javier4.24.65.36.26.57.16.46.36.45.14.73.75.5
Baja CaliforniaMexicali4.14.45.05.66.67.37.06.16.15.54.53.95.5
QuerétaroQuerétaro5.05.76.46.86.96.46.46.46.35.45.04.45.9
PueblaPuebla4.95.56.26.46.15.75.85.85.25.04.74.45.5
HidalgoPachuca4.65.15.66.86.05.75.95.85.34.94.64.25.4

2. Methodology

2.1. Statistical Method

The amount of data collected and the correlation among them are quite important to accomplish an estimation. If the total data is scarce or the correlation between any input and the output is low, the estimation will not be satisfactory.
The correlation analysis is one of the most used and reported statistical methods on scientist and medical researches; its visual representation is known as a dispersion graphic. It is used to prove or reject the existence of a relation between two different variables based on the Pearson correlation coefficient described by Equation (1).
ρ a , b = E a b σ a σ b
where E a b is the crossed correlation between a and b and where σ a 2 = E a 2 and σ b 2 = E b 2 are the variations of the variables a and b, respectively [38].
If “a” and “b” are the two considerate variables, then a dispersion graphic shows the location of each ordered pair in a coordinated system. If most points appear to be close to a straight line, this correlation is known as linear. If most points appear to be close to a curve, the correlation is known as nonlinear. On the other hand, if a clear pattern among the ordered pairs is inexistent, then there is no relation between both variables [39]. In order to represent the curve (curve or line regression) that best fits the behavior of ordered pairs, Equation (2) is used.
y = a 0 + a 1 x + a 2 x 2 + a 3 x 3 + + a m x m
The correlation coefficient or R requires that both a magnitude and a direction be either positive (from 0 to 1) or negative (from −1 to 0). If R gets close to ± 1 , the correlation will be stronger. The correlation does not depend on the direction or the sign: A correlation of 0.57 is equal to a −0.57 one.
It can be also stated that the greater the absolute value of R, the greater the correlation will be. The determination coefficient R 2 is defined as the percentage of the variation of the dependent variable values that can be explained as variations of the independent variable. In other words, a determination coefficient R 2 = 0.23 symbolizes that 23% of the dependent variables is attached to the changes of the independent variable. Therefore, if a correlation factor of R = 0.20 was found between two variables, the determination coefficient would be R 2 = 0.04 so that only 4% of the dependent variable is affected by the variations of the independent variable [40].
In order to create a statistical representative model of the generated power by solar energy, the concept of gradient descent optimization (GDO) was considered due its ability to minimize the model error by the LSR of a linear regression model, which is often used on estimation studies, but it is not as complex to implement as an intelligent technique [29,41,42].
This optimization method has multiple applications. In Reference [43], this method was used to optimize the geometry of the strut-and-tie truss to minimize the difference in the share of resisting actions with respect to the prediction of the multi-action shear model. In Reference [44], a model based on gradient descent is proposed that integrates several important parameters for ranking channels in order to select the best communication channel from a radio spectrum for transmission, enhanced by cognitive radio as an intelligent wireless solution.
The GDO is based on the linear regression method known as the relation among all the input variables and the output one. Depending on the number of input variables, the regression can be simple or multiple [45,46,47]. To organize all the gathered data, every input variable is considered as a column in a matrix called “X” and the output parameter is considered as a vector “y”, as seen in Equation (3).
X = x 11 x 12 x 13 x 1 k x 21 x 22 x 23 x 2 k x 31 x 32 x 33 x 3 k x n 1 x n 2 x n 3 x n k y = y 1 y 2 y n
The purpose of GDO is to find an estimation of the real output through an equation involving all the collected data as shown in Equation (4).
h θ x = θ 0 + θ 1 x 1 + θ 2 x 2 + + θ k x k + ε
where h θ x is the estimated output, x k is the kth input variable, θ k is the characteristic coefficient of every variable, and ϵ is the error between the model and the real data [48].
As can be seen from Equation (4), θ 0 is not a value of influence for any input; however, it does locate the resulting estimation graphic on the Y-axis. If θ 0 is greater than expected, the plot will have a wide gap against the real data in the lower values; nonetheless, if θ 0 is smaller than necessary, the gap between the real and estimated model will be large on the upper zone. Therefore, both cases have a greater error. In order to achieve the minimum error value and instead of a linear regression which uses the least squares method, the GDO finds the right θ 0 by recursive partial derivatives of the cost function described by Equation (5)
J θ = 1 2 m i = 1 m h θ x i y i 2
where m is the total amount of rows in the matrix, x i represents the ith row of “X”, and y i is the value of the ith row of “y”. The gradient descent, denoted by Equation (6), aims to converge to the cost function minimum through its partial derivative. The quickness of the convergence is given by α .
θ k : = θ k α θ k J θ
Equation (7) shows the substitution of Equation (5) into (6) which has to be repeated itself n-times until the convergence is done.
θ k : = θ k α 1 m i = 1 m h θ x i y i · x k i

2.2. Solar Data

The data gathered from the State of Sonora, specifically the city of Hermosillo, were obtained with the support of the University of Sonora (UNISON), while the ones from Mexico City were issued by the Centro de Investigación y de Estudios Avanzados (CINVESTAV) campus Zacatenco, with the solar radiation (global-horizontal), outside temperature, wind speed, and daylight hour (time) as the input meteorological variables and the electric power as the output for both PV systems. Considering that both the solar plants are from different locations and do not have the same arrangement, their respective output power magnitudes vary among them; Hermosillo site has an electric power maximum around 2500 W, while Mexico City site reaches values around 45 kW.
The matrix “X” and the vector “y” from Equation (3) are shown in Equation (8).
X = S o l a r r a d i a t i o n T e m p e r a t u r e W i n d s p e e d T i m e y = E l e c t r i c P o w e r
By using the same variables from both locations, the aim is to achieve an acceptable statistical estimation of the electric power even if those variables are gathered from distant geographic areas.

3. Results and Discussion

3.1. Hermosillo Site

Figure 3 shows the geographical zone where Hermosillo Site (HS) is located. This site covers a period data of 6 months. Figure 4, Figure 5, Figure 6 and Figure 7 depict the dispersion graphic and fitted curve for each input variable against the electric power. In Table 2 and Table 3, the monthly and total results of the characteristic equations for each fitted curve and their determination coefficients are described, respectively.
Equation (7) is applied in order to calculate the values of θ k for Hermosillo by using the Matlab® software. The optimal values found for θ are mentioned in Table 4 and substituted into Equation (4), thereby obtaining its mathematical representation given by Equation (9). Figure 8 is obtained by implementing six months of gathered data into Equation (9). It describes the behavior of the statistical estimation (red) against the real data (green) obtained by the PV system in Hermosillo. For a better appreciation of these results, Figure 9 shows a close-up from 24 November to 24 December 2018. A tracing practically coincident with the rise and fall times for each day as a satisfactory reach of the maximum and minimum values is observed.
h θ = 0.0044 + 2.2387 x 1 + 0.1156 x 2 + 0.0154 x 3 + 0.0022 x 4

3.2. Data from Mexico City

The geographical zone of the Mexico City Site (MCS) is presented in Figure 10. The data covers a period of 6 months. From Figure 11, Figure 12, Figure 13 and Figure 14 the dispersion graphic and fitted curve for each input variable against electric power are shown. In Table 5 and Table 6, the monthly and total results of the characteristic equations for each fitted curve and their determination coefficients are described, respectively.
The determination coefficient values for Mexico City are found applying the same procedure used in the case of Hermosillo, as shown in Table 6. Analogous to the case from Hermosillo, the coefficients θ for Mexico City data and its mathematical representation are displayed in Table 7 and Equation (10), respectively. Figure 15 shows a six months estimation result (red) against the real data behavior (green) from Mexico City. In order to a better appreciation, Figure 16 shows a close-up of Figure 15 covering the same period as that in the case of Hermosillo. The estimation result matches with the rise and fall times of the real power.
h θ = 0.1005 + 55.9815 x 1 + 2.2910 x 2 + 0.1290 x 3 + 0.0533 x 4

3.3. Error Analysis On-Site

As seen in Section 1 and according to Reference [51], the output power of a PV array greatly depends, among other parameters, on solar radiation; however, this variable has an intermittent nature and suffers from rapid fluctuations. In Reference [41,51], the above is considered and some parameters besides the solar radiation such as clear sky or weather data are added; nonetheless, for these cases, the solar radiation is either the estimated output or an estimated input, contrary to this paper where this variable is measured on-site.
Each coefficient θ k from Equations (9) and (10) represents the characteristic magnitude of its respective input variable x k . Its quantity values prove the influence in the mathematical model approximation, establishing an existing relationship between the meteorological variables involved and (responsible for solar fluctuations) the power output.
Comparing Figure 9 and Figure 16, an adequate behavior is observed between both results, achieving a valid statistical estimation. Regardless of the amount of stored data, similar outcomes were obtained from both locations.
This statement is mathematically proved through the correlation and error analysis.

3.3.1. Correlation Analysis

According to Section 2.1, the smaller an error between the estimation and real data, the greater the determination coefficient will be. Figure 17 and Figure 18 show the dispersion graphic of the estimation against the real electric power for Hermosillo and Mexico City respectively. Moreover, in Table 8, the monthly and total fitted curve characteristic equation and determination coefficients for both cases are displayed.
Both statistical models estimate favorably the real electric power for each tested day as we can observe through the values in Table 8. A great closeness can be confirmed between the estimation and the real data for each location, which shows that the implementation of a gradient descent optimization achieves a satisfactory result.

3.3.2. Error Analysis

The less the difference between the estimated and the real value, the better the estimation. Two kinds of errors were applied, where “ P m ” was the measured power, “ P e ” was the estimated power, “s” was the sample in consideration, and “N” was the total amount of samples. The first one was called Mean Absolute Error (MAE) defined by Equation (11) and was used as a standard statistical metric to measure the model performance in meteorology, air quality, and climate research studies [42,51,52]. The second one was the percentage of the MAE known as Mean Absolute Percentage Error (MAPE) defined by Equation (12).
Given that the results shown in Figure 9 and Figure 16 resemble a sinusoidal behavior, P m P e / P m does not present the real error between both signals. Considering Figure 19a, if the measured value is on a high level, then the rate between the P m P e and P m values will be low; however, if the measured value is on a low level, then the same rate will be high or close to 1, as seen in Figure 19b. According to this, the rate is changed regarding the full range of the measured signal (max-min) instead in order to obtain a reliable value. Equation (13) represents the modified MAPE.
M A E = s = 1 N P m P e N
M A P E % = s = 1 N P m P e P m N · 100
M A P E ( r a n g e ) % = s = 1 N P m P e max P m min P m N · 100
With respect to the data from HS and MCS, Table 9 shows the values of the errors mentioned in Equations (11) and (13).
Table 9 shows the errors for both locations, where the MAE result values for HS have a range from 138.1913 W to 299.7928 W and an overall of 200.6071 W. Considering MCS, its MAE goes from 590.6460 W to 1192.3 W, with an overall of 1130.4 W. From these results we can highlight that the values are low according to the total value range for each location, being up to 2500 W and 45 kW for HS and MCS, respectively. This statement is proven by column 4 from Table 9, achieving a MAPE of 6.9036% in general for HS; meanwhile, for MCS, its respective results is a MAPE of 1.9291%. The above demonstrates the estimation robustness of the statistical model.

4. Conclusions

A gradient descent method was used to calculate, through a characteristic equation, the relationship between several variables which can be easily computed using quantitative data.
The results showed that solar radiation and daylight time are relevant in estimating the electric power, with solar radiation being the one with the most influence over the photovoltaic power generation; nonetheless, temperature, wind speed, and daylight hour affect this process, with their inclusion being fundamental in the analysis to achieve a proper electric power estimation.
According to Figure 9 and Figure 16, an acceptable approximation between the statistical simulation results and real-time power generation has achieved. This is proven both by Table 8, where 97.63% of the estimation results matches with the real data for HS and 99.66% match for MCS, and by Table 9, achieving overall error values no greater than 7% and 2% for HS and MCS, respectively. An observable relationship between the determination coefficient values and error results is clear, concluding that the percentage error will be lower while R 2 is higher.
Several causes could perturb the correlation, resulting in a wider dispersion, with the weather conditions as the most important for this study. Contrary to Mexico City, Hermosillo has abrupt changes in climate throughout the year, reaching temperatures around 50 °C and 0 °C in summer and winter, respectively, as well as sudden weather changes from a sunny day to a cloudy or even a rainy day within hours.
However, according to the above and even if the results have a satisfactory behavior, by definition, a statistical estimation will always have an error value. Nevertheless, these error values can be reduced, increasing the gathered input data or by other alternative methods such as intelligent systems which have proven to be an efficient methodology [53,54,55,56]. The above will be discussed elsewhere.

Author Contributions

Research ideas and global design: N.P.-D., J.A.R.-H., and F.A.-V.; analysis of the results: Y.M., E.F.V.-C., and E.J.H.-L.; revision of the document: R.E.C.-L. and J.H.A.-G.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank the University of Sonora (UNISON), Centro de Investigación y de Estudios Avanzados (CINVESTAV) campus Zacatenco, and Consejo Nacional de Ciencia y Tecnología (CONACYT), with the mobility scholarship 291249.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ANNArtificial Neural Network
CSPConcentrated Solar Power Plants
GDOGradient Descent Optimization
HEIHigher Education Institutions
HSHermosillo Site
kWh/m2Quantity of kilowatts during an hour striking a squared meter area
LSRLeast Square Regression
MAEMean Absolute Error
MAPEMean Absolute Percentage Error
MCSMexico City Site
PVPhotovoltaic
SHCSolar Heating and Cooling
WWatt. It is a unit of power defined as a derived unit of 1 joule per second (J/s)

References

  1. D’Adamo, I. The Profitability of Residential Photovoltaic Systems. A New Scheme of Subsidies Based on the Price of CO2 in a Developed PV Market. Soc. Sci. 2018, 7, 148. [Google Scholar] [CrossRef]
  2. SEIA. Solar Energy Industries Association. 2018. Available online: https://www.seia.org/initiative-topics/solar-technologies (accessed on 12 April 2019).
  3. Järvelä, M.; Valkealahti, S. Ideal operation of a photovoltaic power plant equipped with an energy storage system on electricity market. Appl. Sci. 2017, 7, 749. [Google Scholar] [CrossRef]
  4. Awan, A.B.; Zubair, M.; Abokhalil, A.G. Solar energy resource analysis and evaluation of photovoltaic system performance in various regions of Saudi Arabia. Sustainability 2018, 10, 1129. [Google Scholar] [CrossRef]
  5. Zahedi, A. Solar photovoltaic (PV) energy; latest developments in the building integrated and hybrid PV systems. Renew. Energy 2006, 31, 711–718. [Google Scholar] [CrossRef]
  6. Perea-Moreno, A.J.; Hernandez-Escobedo, Q.; Garrido, J.; Verdugo-Diaz, J.D. Stand-Alone Photovoltaic System Assessment in Warmer Urban Areas in Mexico. Energies 2018, 11, 284. [Google Scholar] [CrossRef]
  7. Trujillo-Camacho, M.; García-Gómez, C.; Hinojosa-Palafox, J.; Castillón-Barraza, F. Evaluación de compositos TiO2/clinoptilolita en la fotodegradación del tinte MV-2B en un reactor-concentrador solar cpc. Revista Mexicana de Ingeniería química 2010, 9, 139–149. [Google Scholar]
  8. López-Sosa, L.; Hernández-Ramírez, L.; González-Avilés, M.; Servín-Campuzano, H.; Zárate-Medina, J. Desarrollo de un recubrimiento absorbente solar de bajo costo basado en hollín de biomasa forestal: Caracterización térmica y aplicación en un sistema de cocción solar. Revista Mexicana de Ingeniería Química 2018, 17, 651–668. [Google Scholar] [CrossRef]
  9. González-Avilés, M.; López-Sosa, L.; Servín-Campuzano, H.; González-Pérez, D. Adopción tecnológica sustentable de cocinas solares en comunidades indígenas y rurales de Michoacán. Revista Mexicana de Ingeniería Química 2017, 16, 273–282. [Google Scholar]
  10. Saga, T. Crystalline and Polycrystalline Silicon PV Technology. NPG Asia Mater. 2010, 2, 96–102. [Google Scholar] [CrossRef]
  11. Kazem, H.A.; Yousif, J.H.; Chaichan, M.T. Modeling of daily solar energy system prediction using support vector machine for Oman. Int. J. Appl. Eng. Res. 2016, 11, 10166–10172. [Google Scholar]
  12. Zhang, W.-T.; Wang, S.; Du, X.-H. Research of power prediction about photovoltaic power system: Based on BP neural network. J. Environ. Prot. Ecol. 2017, 18, 1614–1623. [Google Scholar]
  13. Shin, H.; Geem, Z.W. Optimal Design of a Residential Photovoltaic Renewable System in South Korea. Appl. Sci. 2019, 9, 1138. [Google Scholar] [CrossRef]
  14. Hogg, D.C.; Decker, M.; Guiraud, F.; Earnshaw, K.; Merritt, D.; Moran, K.; Sweezy, W.; Strauch, R.; Westwater, E.; Little, C. An automatic profiler of the temperature, wind and humidity in the troposphere. J. Clim. Appl. Meteorol. 1983, 22, 807–831. [Google Scholar] [CrossRef]
  15. Lawrence, M.G. The relationship between relative humidity and the dewpoint temperature in moist air: A simple conversion and applications. Bul. Am. Meteorol. Soc. 2005, 86, 225–234. [Google Scholar] [CrossRef]
  16. Dudley, B. BP Energy Outlook 2040. BP plc 2019. Available online: http://oilproduction.net/files/OilProduction-bp-energy-outlook-2019.pdf (accessed on 12 April 2019).
  17. IRENA. Renewable Power Generation Costs in 2017; International Renewable Energy Agency: Abu Dhabi, UAE, 2018; Available online: https://www.irena.org/-/media/Files/IRENA/Agency/Publication/2018/Jan/IRENA_2017_Power_Costs_2018.pdf (accessed on 12 April 2019).
  18. Rasero, C.M. Energía Solar Fotovoltaica. Energía Solar Fotovoltaica, Situación Actual 2011. Available online: https://static.eoi.es/savia/documents/componente75553.pdf (accessed on 12 April 2019).
  19. Philibert, C. Solar Energy Perspectives; OECD: Paris, France, 2011. [Google Scholar]
  20. Elshurafa, A.M.; Albardi, S.R.; Bigerna, S.; Bollino, C.A. Estimating the learning curve of solar PV balance–of–system for over 20 countries: Implications and policy recommendations. J. Clean. Prod. 2018, 196, 122–134. [Google Scholar] [CrossRef]
  21. Nemet, G.F.; O’Shaughnessy, E.; Wiser, R.; Darghouth, N.; Barbose, G.; Gillingham, K.; Rai, V. Characteristics of low-priced solar PV systems in the US. Appl. Energy 2017, 187, 501–513. [Google Scholar] [CrossRef]
  22. Acceso a Datos de Radiación Solar en España. 2018. Available online: http://www.adrase.com (accessed on 12 April 2019).
  23. The World Bank, Solar Resource Data: Solargis, World. 2017. Available online: https://globalsolaratlas.info/downloads/world (accessed on 12 April 2019).
  24. The World Bank, Solar Resource Data: Solargis, Mexico. 2017. Available online: https://globalsolaratlas.info/downloads/mexico (accessed on 12 April 2019).
  25. Yaneva, M.; Tisheva, P.; Tsanova, T. Mirecweek—The Big Mexico Renewable Energy Report. 2018. Available online: http://www.awex-export.be/files/library/Infos-sectorielles/Ameriques/2017/MEXIQUE/Mirec-Report-2018-The-BIG-Mexico-renewable-energy-report-ENG.pdf (accessed on 12 April 2019).
  26. Della Ceca, L.; Micheletti, M.; Freire, M.; Garcia, B.; Mancilla, A.; Salum, G.; Crinó, E.; Piacentini, R. Solar and Climatic High Performance Factors for the Placement of Solar Power Plants in Argentina Andes Sites—Comparison With African and Asian Sites. J. Sol. Energy Eng. 2019, 141, 041004. [Google Scholar] [CrossRef]
  27. Sağlam, Ş. Meteorological parameters effects on solar energy power generation. WSEAS Trans. Circuits Syst. 2010, 9, 637–649. [Google Scholar]
  28. Nam, S.; Hur, J. Probabilistic Forecasting Model of Solar Power Outputs Based on the Naïve Bayes Classifier and Kriging Models. Energies 2018, 11, 2982. [Google Scholar] [CrossRef]
  29. Ding, M.; Wang, L.; Bi, R. An ANN-based approach for forecasting the power output of photovoltaic system. Procedia Environ. Sci. 2011, 11, 1308–1315. [Google Scholar] [CrossRef]
  30. Horan, W.; Shawe, R.; O’Regan, B. Ireland’s Transition towards a Low Carbon Society: The Leadership Role of Higher Education Institutions in Solar Photovoltaic Niche Development. Sustainability 2019, 11, 558. [Google Scholar] [CrossRef]
  31. Viscidi, L. Mexico’s Renewable Energy Future. 2018. Available online: https://www.wilsoncenter.org/sites/default/files/mexico_renewable_energy_future_0.pdf (accessed on 12 April 2019).
  32. SENER. Prospectivas de Energías Renovables 2012–2026. 2012. Available online: https://www.gob.mx/cms/uploads/attachment/file/62954/Prospectiva_de_Energ_as_Renovables_2012-2026.pdf (accessed on 12 April 2019).
  33. Gordon, G.; Tibshirani, R. Gradient descent revisited. Optimization 2012, 10, 725. [Google Scholar]
  34. Hernández, L.; Baladrón, C.; Aguiar, J.M.; Calavia, L.; Carro, B.; Sánchez-Esguevillas, A.; Cook, D.J.; Chinarro, D.; Gómez, J. A study of the relationship between weather variables and electric power demand inside a smart grid/smart world framework. Sensors 2012, 12, 11571–11591. [Google Scholar] [CrossRef]
  35. Ooi, S.C.; Mardiana, A.; Yusup, Y. Analysing Meteorological Variables, Energy Consumption and Occupant Behaviour in an Office Building in Hot-Humid Climate Zone. Int. J. Sci. Res. IJSR 2015, 4, 88–93. [Google Scholar]
  36. Hor, C.L.; Watson, S.J.; Majithia, S. Analyzing the impact of weather variables on monthly electricity demand. IEEE Trans. Power Syst. 2005, 20, 2078–2085. [Google Scholar] [CrossRef]
  37. Mendoza, C. Viabilidad Técnica-Económica de una Central Solar Termoeléctrica de Colectores Cilíndricos Parabólicos para su Implementación de México. UNAM, Marzo. 2012. Available online: http://132.248.52.100:8080/xmlui/handle/132.248.52.100/277 (accessed on 12 April 2019).
  38. Benesty, J.; Chen, J.; Huang, Y.; Cohen, I. Pearson correlation coefficient. In Noise Reduction in Speech Processing; Springer: Berlin, Germany, 2009; pp. 1–4. [Google Scholar]
  39. Spiegel, M.; Stephens, L. Estadística, 4th ed.; Graw Hill: Contadero, Mexico, 2009; ISBN 978-6-071-51188-1. Available online: https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=2ahUKEwj12s7rtLDhAhWB7Z8KHQIiDwAQFjAAegQIAxAC&url=http%3A%2F%2Fensfep.edu.mx%2Fenlinea%2Fpluginfile.php%2F1531%2Fmod_folder%2Fcontent%2F0%2FEstad%25C3%25ADstica.%2520Serie%2520Schaum-%2520 (accessed on 12 April 2019).
  40. Taylor, R. Interpretation of the correlation coefficient: A basic review. J. Diagn. Med. Sonogr. 1990, 6, 35–39. [Google Scholar] [CrossRef]
  41. Clack, C.T. Modeling solar irradiance and solar PV power output to create a resource assessment using linear multiple multivariate regression. J. Appl. Meteorol. Climatol. 2017, 56, 109–125. [Google Scholar] [CrossRef]
  42. Aggarwal, S.; Saini, L. Solar energy prediction using linear and non-linear regularization models: A study on AMS (American Meteorological Society) 2013–14 Solar Energy Prediction Contest. Energy 2014, 78, 247–256. [Google Scholar] [CrossRef]
  43. Bairán, J.M.; Marí, A.; Cladera, A. Analysis of shear resisting actions by means of optimization of strut and tie models taking into account crack patterns. Hormigón y acero 2018, 69, 197–206. [Google Scholar] [CrossRef]
  44. Arjoune, Y.; Mrabet, Z.; Kaabouch, N. Multi-Attributes, Utility-Based, Channel Quality Ranking Mechanism for Cognitive Radio Networks. Appl. Sci. 2018, 8, 628. [Google Scholar] [CrossRef]
  45. Montgomery, D.C. Design and Analysis of Experiments; John Wiley & Sons: Hoboken, NJ, USA, 2017; ISBN 978-1-119-11347-8. Available online: https://www.wiley.com/en-mx/Design+and+Analysis+of+Experiments%2C+9th+Edition-p-9781119113478 (accessed on 12 April 2019).
  46. Gardiner, S.K.; Crabb, D.P. Examination of different pointwise linear regression methods for determining visual field progression. Investig. Ophthalmol. Vis. Sci. 2002, 43, 1400–1407. [Google Scholar]
  47. Verma, S.; Bartosova, A.; Markus, M.; Cooke, R.; Um, M.J.; Park, D. Quantifying the Role of Large Floods in Riverine Nutrient Loadings Using Linear Regression and Analysis of Covariance. Sustainability 2018, 10, 2876. [Google Scholar] [CrossRef]
  48. Ng, A. CS229 Lecture Notes. 2000, Volume 1, pp. 1–3. Available online: http://backspaces.net/temp/ML/CS229.pdf (accessed on 12 April 2019).
  49. Google Maps, UNISON. 2019. Available online: https://www.google.com/maps/place/Universidad+de+Sonora+-+UNISON/@23.9391074,-102.172878,2504267m/data=!3m1!1e3!4m5!3m4!1s0x86ce8447973925f7:0xcf527709b7555a3!8m2!3d29.0834761!4d-110.9603621 (accessed on 12 April 2019).
  50. Google Maps, CINVESTAV. 2019. Available online: https://www.google.com/maps/place/CINVESTAV+-+IPN/@23.9391074,-102.172878,2504267m/data=!3m1!1e3!4m5!3m4!1s0x85d1f77a076bd911:0x63af0ad86ca91f65!8m2!3d19.5099425!4d-99.129559 (accessed on 12 April 2019).
  51. Dev, S.; AlSkaif, T.; Hossari, M.; Godina, R.; Louwen, A.; van Sark, W. Solar Irradiance Forecasting Using Triple Exponential Smoothing. In Proceedings of the 2018 International Conference on Smart Energy Systems and Technologies (SEST), Sevilla, Spain, 10–12 September 2018; pp. 1–6. [Google Scholar]
  52. Chai, T.; Draxler, R.R. Root mean square error (RMSE) or mean absolute error (MAE)?—Arguments against avoiding RMSE in the literature. Geosci. Model Dev. 2014, 7, 1247–1250. [Google Scholar] [CrossRef]
  53. Peterson, M.G. Intelligent medical systems and the interface with statistics. In Proceedings of the 11th IEEE Symposium on Computer-Based Medical Systems, Lubbock, TX, USA, 12–14 June 1998; p. 300. [Google Scholar]
  54. Chen, S.H.; Jakeman, A.J.; Norton, J.P. Artificial intelligence techniques: An introduction to their use for modelling environmental systems. Math. Comput. Simul. 2008, 78, 379–400. [Google Scholar] [CrossRef]
  55. Mosavi, A.; Ozturk, P.; Chau, K.w. Flood prediction using machine learning models: Literature review. Water 2018, 10, 1536. [Google Scholar] [CrossRef]
  56. Inman, R.H.; Pedro, H.T.; Coimbra, C.F. Solar forecasting methods for renewable energy integration. Prog. Energy Combust. Sci. 2013, 39, 535–576. [Google Scholar] [CrossRef]
Figure 1. The worldwide horizontal solar radiation [23].
Figure 1. The worldwide horizontal solar radiation [23].
Applsci 09 01649 g001
Figure 2. The Mexican solar radiation atlas [24].
Figure 2. The Mexican solar radiation atlas [24].
Applsci 09 01649 g002
Figure 3. The geographical location of the Hermosillo Site (HS) [49] (29.0843293 N, −110.9583558 W).
Figure 3. The geographical location of the Hermosillo Site (HS) [49] (29.0843293 N, −110.9583558 W).
Applsci 09 01649 g003
Figure 4. The Power vs. Solar radiation of HS.
Figure 4. The Power vs. Solar radiation of HS.
Applsci 09 01649 g004
Figure 5. The Power vs. Temperature of HS.
Figure 5. The Power vs. Temperature of HS.
Applsci 09 01649 g005
Figure 6. The Power vs. Wind speed of HS.
Figure 6. The Power vs. Wind speed of HS.
Applsci 09 01649 g006
Figure 7. The Power vs. Time of HS.
Figure 7. The Power vs. Time of HS.
Applsci 09 01649 g007
Figure 8. The results obtained from HS by gradient descent.
Figure 8. The results obtained from HS by gradient descent.
Applsci 09 01649 g008
Figure 9. A close-up of the results from HS.
Figure 9. A close-up of the results from HS.
Applsci 09 01649 g009
Figure 10. The geographical location of Mexico City Site (MCS) [50] (19.5099425 N, −99.1317477 W).
Figure 10. The geographical location of Mexico City Site (MCS) [50] (19.5099425 N, −99.1317477 W).
Applsci 09 01649 g010
Figure 11. The Power vs. Solar radiation of MCS.
Figure 11. The Power vs. Solar radiation of MCS.
Applsci 09 01649 g011
Figure 12. The Power vs. Temperature of MCS.
Figure 12. The Power vs. Temperature of MCS.
Applsci 09 01649 g012
Figure 13. The Power vs. Wind speed of MCS.
Figure 13. The Power vs. Wind speed of MCS.
Applsci 09 01649 g013
Figure 14. The Power vs. Time of MCS.
Figure 14. The Power vs. Time of MCS.
Applsci 09 01649 g014
Figure 15. The results obtained from MCS by gradient descent.
Figure 15. The results obtained from MCS by gradient descent.
Applsci 09 01649 g015
Figure 16. A close-up of the results from MCS.
Figure 16. A close-up of the results from MCS.
Applsci 09 01649 g016
Figure 17. The estimated vs. real photovoltaic power of HS.
Figure 17. The estimated vs. real photovoltaic power of HS.
Applsci 09 01649 g017
Figure 18. The estimated vs. real photovoltaic power of MCS.
Figure 18. The estimated vs. real photovoltaic power of MCS.
Applsci 09 01649 g018
Figure 19. The difference errors between the high and low values: (a) Error value when measured data is high; (b) Error value when measured data is low.
Figure 19. The difference errors between the high and low values: (a) Error value when measured data is high; (b) Error value when measured data is low.
Applsci 09 01649 g019
Table 2. The fitted curve characteristic equation of every input from HS.
Table 2. The fitted curve characteristic equation of every input from HS.
MonthInput Variables
Solar radiationTemperature
Aug f x = 2.245 x + 92.59 f x = 17.98 x 2 + 1274 x 20,920
Sep f x = 2.308 x + 28.45 f x = 11.43 x 2 + 796.7 x 12,210
Oct f x = 2.276 x + 50.32 f x = 8.751 x 2 + 556.5 x 7217
Nov f x = 2.326 x 29.48 f x = 6.074 x 2 + 354.8 x 3505
Dec f x = 2.25 x 36.17 f x = 0.08391 x 2 + 45.57 x + 219.1
Jan f x = 2.398 x + 8.88 f x = 2.016 x 2 + 144.8 x 752.7
TOTAL f x = 2.296 x 1.747 f x = 1.751 x 2 + 124.2 x 581.9
Wind speedHour
Aug f x = 60.94 x 2 + 438.2 x + 753.3 f x = 43,250 x 2 + 44,070 x 8893
Sep f x = 21.89 x 2 + 111.6 x + 1352 f x = 42,970 x 2 + 43,290 x 8593
Oct f x = 20.63 x 2 + 118.4 x + 1292 f x = 46,940 x 2 + 46,740 x 9314
Nov f x = 4.208 x 2 + 12.31 x + 1348 f x = 50,610 x 2 + 50,160 x 10,180
Dec f x = 11.68 x 2 + 120.2 x + 1005 f x = 51,600 x 2 + 51,860 x 10,900
Jan f x = 14.18 x 2 + 206.8 x + 886.7 f x = 55,500 x 2 + 56,360 x 11,960
TOTAL f x = 14.47 x 2 + 138.7 x + 1096 f x = 49,950 x 2 + 50,110 x 10,360
Table 3. The determination coefficients (R2) of every input from HS.
Table 3. The determination coefficients (R2) of every input from HS.
MonthSolar RadiationTemperatureWind SpeedHour
Aug0.99170.16790.05410.9825
Sep0.99810.15020.02420.9968
Oct0.95230.16640.03800.9850
Nov0.97080.19550.01350.9962
Dec0.98510.09650.01360.8865
Jan0.96090.13900.05650.9760
TOTAL0.95320.10770.01400.8098
Table 4. The values of θ k from HS by gradient descent.
Table 4. The values of θ k from HS by gradient descent.
θ k HS
θ 0 0.0044
θ 1 2.2387
θ 2 0.1156
θ 3 0.0154
θ 4 0.0022
Table 5. The fitted curve characteristic equation of every input from MCS.
Table 5. The fitted curve characteristic equation of every input from MCS.
MonthInput Variables
Solar radiationTemperature
Jul f x = 53.58 x + 135 f x = 20.07 x 2 + 3392 x 44,530
Aug f x = 56.62 x + 241.2 f x = 43.78 x 2 + 4162 x 51,030
Sep f x = 59.48 x + 245.7 f x = 9.184 x 2 + 1875 x 25,960
Oct f x = 57.58 x + 189.8 f x = 6.143 x 2 + 2115 x 18,990
Nov f x = 57.64 x + 1.454 f x = 132 x 2 + 5978 x 36,790
Dec f x = 53.52 x 7.362 f x = 104.9 x 2 + 4312 x 18,320
TOTAL f x = 56.75 x + 65.85 f x = 31.9 x 2 + 203.1 x + 2471
Wind speedHour
Jul f x = 173.8 x 2 + 531.6 x + 20,420 f x = 584,100 x 2 + 632,300 x 137,600
Aug f x = 157.7 x 2 + 325.2 x + 21,850 f x = 599,500 x 2 + 649,400 x 141,400
Sep f x = 156.4 x 2 + 459.8 x + 19,960 f x = 632,500 x 2 + 686,200 x 153,700
Oct f x = 113 x 2 + 47.06 x + 23,490 f x = 693,400 x 2 + 735,600 x 159,300
Nov f x = 66.96 x 2 1894 x + 28,680 f x = 918,400 x 2 + 908,400 x 181,500
Dec f x = 44.08 x 2 706.9 x + 22,910 f x = 837,800 x 2 + 843,100 x 175,100
TOTAL f x = 130.7 x 2 + 61.23 x + 22,720 f x = 617,100 x 2 + 650,100 x + 136,600
Table 6. The determination coefficients (R2) of every input from MCS.
Table 6. The determination coefficients (R2) of every input from MCS.
MonthSolar RadiationTemperatureWind SpeedHour
Jul0.99970.49520.02950.6032
Aug0.99900.32610.02140.5652
Sep0.99500.36840.00930.5014
Oct0.99230.21480.01770.5803
Nov0.99950.16740.03200.8684
Dec0.99950.14400.00150.8865
TOTAL0.99890.20520.01780.5725
Table 7. The values of θ k from MCS by gradient descent.
Table 7. The values of θ k from MCS by gradient descent.
θ k MCS
θ 0 0.1005
θ 1 55.9815
θ 2 2.2910
θ 3 0.1290
θ 4 0.0533
Table 8. The monthly and total fitted curve characteristic equation and determination coefficient values of estimation for both locations.
Table 8. The monthly and total fitted curve characteristic equation and determination coefficient values of estimation for both locations.
LocationMonthFitted CurveR2
Aug f x = 1.005 x + 84.92 0.9916
Sep f x = 1.032 x + 20.92 0.9981
Oct f x = 1.031 x + 44.45 0.9877
HSNov f x = 1.022 x 29.74 0.9926
Dec f x = 1.05 x 37.57 0.9849
Jan f x = 1.009 x + 6.707 0.9982
TOTAL f x = 1.025 x 4.633 0.9763
Jul f x = 1.003 x + 91.51 0.9989
Aug f x = 1.01 x + 184.3 0.9961
Sep f x = 0.9843 x + 525.7 0.9946
MCSOct f x = 1.003 x + 140.6 0.9997
Nov f x = 1.007 x 40.05 0.9982
Dec f x = 1.004 x 44.28 0.9995
TOTAL f x = 1.014 x + 19.18 0.9966
Table 9. The error values between the estimated and real data.
Table 9. The error values between the estimated and real data.
LocationMonthMAEMAPE (Range)
Aug299.7928 W11.9436%
Sep203.7386 W7.5643%
Oct267.9841 W9.2382%
HSNov138.1913 W5.1418%
Dec178.9316 W7.0381%
Jan184.5001 W6.7978%
TOTAL200.0924 W6.8859%
Jul590.6460 W1.0480%
Aug1192.3 W2.0392%
Sep726.2990 W1.2646%
MCSOct741.8355 W1.2660%
Nov741.4300 W1.3534%
Dec812.8113 W1.7237%
TOTAL1130.4 W1.9291%

Share and Cite

MDPI and ACS Style

Ruz-Hernandez, J.A.; Matsumoto, Y.; Arellano-Valmaña, F.; Pitalúa-Díaz, N.; Cabanillas-López, R.E.; Abril-García, J.H.; Herrera-López, E.J.; Velázquez-Contreras, E.F. Meteorological Variables’ Influence on Electric Power Generation for Photovoltaic Systems Located at Different Geographical Zones in Mexico. Appl. Sci. 2019, 9, 1649. https://doi.org/10.3390/app9081649

AMA Style

Ruz-Hernandez JA, Matsumoto Y, Arellano-Valmaña F, Pitalúa-Díaz N, Cabanillas-López RE, Abril-García JH, Herrera-López EJ, Velázquez-Contreras EF. Meteorological Variables’ Influence on Electric Power Generation for Photovoltaic Systems Located at Different Geographical Zones in Mexico. Applied Sciences. 2019; 9(8):1649. https://doi.org/10.3390/app9081649

Chicago/Turabian Style

Ruz-Hernandez, Jose A., Yasuhiro Matsumoto, Fernando Arellano-Valmaña, Nun Pitalúa-Díaz, Rafael Enrique Cabanillas-López, José Humberto Abril-García, Enrique J. Herrera-López, and Enrique Fernando Velázquez-Contreras. 2019. "Meteorological Variables’ Influence on Electric Power Generation for Photovoltaic Systems Located at Different Geographical Zones in Mexico" Applied Sciences 9, no. 8: 1649. https://doi.org/10.3390/app9081649

APA Style

Ruz-Hernandez, J. A., Matsumoto, Y., Arellano-Valmaña, F., Pitalúa-Díaz, N., Cabanillas-López, R. E., Abril-García, J. H., Herrera-López, E. J., & Velázquez-Contreras, E. F. (2019). Meteorological Variables’ Influence on Electric Power Generation for Photovoltaic Systems Located at Different Geographical Zones in Mexico. Applied Sciences, 9(8), 1649. https://doi.org/10.3390/app9081649

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