Acessibilidade / Reportar erro

Analysis of Measurements of Direct Normal and Global Solar Radiation for Maceió, Alagoas State, Northeastern Brazil

Análise das Medidas de Radiação Solar Direta Normal e Global para Maceió, Estado de Alagoas, Nordeste do Brasil

Abstract

Studying solar radiation is essential for human knowledge, since it is present in practically all its activities. Thus, the aim of this work was to analyze the climatic and seasonal variation of direct normal and global solar radiation in the region of Maceió, Alagoas State, Northeastern Brazil with sky conditions characterized by clearness index (Kt). The Kt was determined by the ratio between global solar irradiance and solar irradiance at the top of the atmosphere. The highest occurrences of daily direct normal solar irradiance under conditions of Kt ≥ 0.6 were recorded between 400 W m−2 and 700 W m−2 for all seasons. Under conditions of 0.4 ≤ Kt < 0.6, the daily direct normal solar irradiance occurred between 200 W m−2 and 500 W m−2 and for conditions of Kt < 0.4, its maximum value was 200 W m−2. It was observed that the levels of solar incidence in the study region depend on cloud cover conditions, with little influence of seasonality.

Keywords
solar irradiance; climatic and seasonal analysis; sky conditions

Resumo

Estudar a radiação solar é essencial para o conhecimento humano, visto que ela está presente em praticamente todas suas atividades. Assim, o objetivo desse trabalho é analisar a variação climática e sazonal da radiação solar direta normal e global na região de Maceió, estado de Alagoas, Nordeste do Brasil com condições do céu caracterizada pelo índice de claridade (Kt). O Kt foi determinado pela razão entre irradiância solar global e irradiância solar no topo da atmosfera. As maiores ocorrências de irradiância solar direta normal diária em condições de Kt ≥ 0,6 foram registradas entre 400 W m−2 e 700 W m−2 em todas as estações do ano. Sob condições de 0,4 ≤ Kt < 0,6, a irradiância solar direta normal diária ocorreu entre 200 W m−2 e 500 W m−2 e para condições de Kt < 0,4, seu valor máximo foi de 200 W m−2. Observou-se que os níveis de incidência solar na região de estudo dependem das condições de cobertura de nuvens, com pouca influência da sazonalidade.

Palavras-chave:
irradiância solar; análise climática e sazonal; condição do céu

1. Introduction

The radiation emitted by the Sun that arrives at Earth (top of the atmosphere and terrestrial surface) is the main source of energy that influences practically all human activities (Liou, 2002Liou, K.N. An Introduction to Atmospheric Radiation. 2. Ed. San Diego: Academic Press, 2002, 583p.; Muneer, 2004Muneer, T. Solar Irradiation and Daylight Models. 2nd ed. Amsterdam: Elsevier, 2004.; Santos et al., 2014Santos, C.M.; Souza, J.L.; Ferreira Junior, R.A.; Tiba, C.; Melo, R.O. et al. On modeling global solar irradiation using air temperature for Alagoas State, Northeastern Brazil. Energy, v. 71, p. 388-398, 2014.). The knowledge of the global solar irradiance (IG) which is measured on a horizontal surface, direct normal solar irradiance (ID) and diffuse solar irradiance (Id) has great importance in several scientific and technological areas, e.g. meteorology, climatology, energy generation, agrometeorology, and others (Vick et al., 2012Vick, B.D.; Daryl, R.; Myers, D.R.; Boyson, W.E. Using direct normal irradiance models and utility electrical loading to assess benefit of a concentrating solar power plant. Solar Energy, v. 86, n. 12, p. 3519-3530, 2012.; Tabari and Talaee, 2014Tabari, H.; Talaee, P.H. Sensitivity of evapotranspiration to climatic change in different climates. Glob. Planet. Chang., v. 115, p. 16-23, 2014.; Evseev and Kudish, 2015Evseev, E.G.; Kudish, A.I. Analysis of solar irradiation measurements at Beer Sheva, Israel from 1985 through 2013. Energy Conversion and Management, v. 97, p. 307-314, 2015.; Ferreira Junior et al., 2015Ferreira Junior, R.A.; Souza, J.L.; Lyra, G.B.; Escobedo, J.F.; Santos, M.V.C. C. Energy conversion efficiency in sugarcane under two row spacings in northeast of Brazil. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 19, n. 8, p. 741-747, 2015.; Santos et al., 2016Santos, C.M.; Escobedo, J.F.; Teramoto, E.T.; Silva, S.H.M.G. Assessment of ANN and SVM models for estimating normal direct irradiation (Hb). Energy Conversion and Management, v. 126, p. 826-836, 2016.).

Thus, solar radiation studies require reliable and readily available data series for applications at the site of interest. The instrumentation to obtain IG and ID time series is given by the installation of a network of solar monitoring stations with radiometers and data acquisition system. The amount and spectral distribution of IG and ID depend on astronomical, atmospheric and local parameters such as: solar zenith angle, aerosol concentration and properties, amount of ozone, clouds, altitude of the site and surface albedo (Kaskaoutis and Kambezidis, 2009Kaskaoutis, D.G.; Kambezidis, H.D. The diffuse-to-global and diffuse-to-direct-beam spectral irradiance ratios as turbidity indexes in an urban environment. Journal of Atmospheric and Solar-Terrestrial Physics, v. 71, n. 2, p. 246-256, 2009.). ID is very sensitive to the optical depth of the aerosol and exhibits high frequency variation to the dynamic effects of clouds, showing greater spatial and temporal variability than IG, making it difficult to characterize daily and seasonal variability of solar radiation and its components (Kotti et al., 2014Kotti, M.C.; Argiriou, A.A.; Kazantzidis, A. Estimation of direct normal irradiance from measured global and corrected diffuse horizontal irradiance. Energy, v. 70, n. 1, p. 382-392, 2014.).

The climatic characterization of the solar radiation requires a series of long-term measurements. However, the availability of solar databases, especially long time and continuous ID is scarce. With the availability of measurements, different researchers studied the climatic characteristics of solar irradiance, especially the work of Liu and Jordan (1960)Liu, B.Y.H.; Jordan, R.C. The interrelationship and characteristic distribution of direct, diffuse and total solar radiation. Solar Energy, v. 4, n. 3, p. 1-19, 1960.. In Beer Sheva - Israel, Evseev and Kudish (2015)Evseev, E.G.; Kudish, A.I. Analysis of solar irradiation measurements at Beer Sheva, Israel from 1985 through 2013. Energy Conversion and Management, v. 97, p. 307-314, 2015. analyzed solar irradiance from 1985 to 2013 and observed that the site has high frequency of clear sky days. In Seville - Spain, Moreno-Tejera et al. (2016)Moreno-Tejera, S.; Silva-Perez, M.A.; Lillo-Bravo, I.; Ramırez-Santigosa, L. L. Solar resource assessment in Seville, Spain. Statistical characterisation of solar radiation at different time resolutions. Solar Energy, v. 132, p. 430-441, 2016. showed an analysis of solar irradiance variation, applying a statistical characterization in different time partitions.

Despite the importance and interest of the scientific and governmental community, the knowledge of solar irradiance components on the surface is still in its beginning in the Brazilian territory, mainly because it is a country of great territorial extent. In Brazil, there are few solarimetric stations that simultaneously perform high quality measurements of IG and ID components. Studies analyzing databases of solar components are not so often found in specialized literature; however global solar radiation for several sites in Northeastern Brazil (NEB) has been evaluated (Tiba, 2001Tiba, C. Solar radiation in the Brazilian Northeast. Renewable Energy, v. 22, p. 565-578, 2001.). The seasonal variations of the diurnal evolution of solar irradiance components based on measurements collected on the surface of the city of São Paulo - Brazil (Oliveira et al., 2002Oliveira, A.P.; Machado, J. A.; Escobedo, J. F.; Soares, J. Diurnal evolution of solar irradiation at the surface in the city of Sao Paulo: Seasonal variation and modeling, Theor. Appl. Climatol., v. 71, n. 3-4, p. 231-249, 2002.). The global solar irradiation was evaluated in Maceió (NEB) (Souza et al., 2005Souza, J.L.; Nicácio, R.M.; Moura, M.A.L. Global solar radiation measurements in Maceió, Brazil. Renewable Energy, v. 30, n. 8, p. 1203-1220, 2005.), showing a lower monthly average daily value in July and higher in November, with dry season (October-February) average hourly values 12.26% higher than values in the rainy season (April-August). Analyzes and comparisons of solar components were made for the city of São Paulo and rural area of Botucatu-SP (Codato et al., 2008Codato, G.; Oliveira, A.P.; Soares, J.; Escobedo, J.F.; Gomes, E.N.; Pai, A.D. D. Global and diffuse solar irradiances in urban and rural areas in southeast Brazil, Theor. Appl. Climatol., v. 93, p. 57-73, 2008.), in which, although located in similar latitude and altitude, the seasonal evolution of the daily values indicates that São Paulo receives, during the days of open sky, 7.8% less in August and 5.1% less in June than Botucatu. An observational characterization of the diurnal and seasonal variation of the solar irradiance components for Rio de Janeiro has been Marques Filho et al. (2016)Marques Filho, E.P.; Oliveira, A.P.; Vita, W.A. et al. Global, diffuse and direct solar radiation at the surface in the city of Rio de Janeiro: Observational characterization and empirical modeling. Renewable Energy, v. 91, p. 64-74, 2016..

The research site is located in the metropolitan region of Maceió, with sugarcane cultivation in the vicinity, bathed by the Atlantic Ocean and large lagoons; thus, the local atmosphere is polluted by particulates originating from industries, burning and local fires, and elevation of the water vapor concentration that attenuates IG and, mainly, ID. These very small solid or liquid particles, called aerosols, are determinants of the spread of ID. On the other hand, water vapor acts in the earth's atmosphere as an ID absorber. The combined action of aerosol concentration and water vapor content significantly affects the magnitude and variability of ID, influencing the transparency level of the atmosphere. The importance of solar radiation characterization for the region where there is a growing demand for solarimetric information for solar energy, fundamental for establishing public policies in the area of renewable energy is also highlighted. Thus, this paper presents a climatic and seasonal analysis of the effects of sky cover on IG and ID values in the metropolitan region of Maceió between years of 2011 and 2014.

2. Material and Methods

2.1. Study area

The measurements used in this study were obtained from the solarimetric station located at the Federal University of Alagoas (9°33' S; 35°46' W and 100 m a.s.l.), in Maceió, capital of the State of Alagoas, located in eastern NEB (Fig. 1). The microregion of Maceió is located on the coast, with variation of altitude between (0 and 20) m in the lower region (coastal plain), between (20 and 180) m on the slopes and coastal boards, reaching 300 m at the northern end of the municipality (Fig. 1). This variation can cause changes in air temperature and winds. The solarimetric station is installed at 15 km from the Atlantic Ocean, with herbaceous vegetation and Atlantic Rainforest and hot and humid climate.

Figure 1
Geographic location of northeastern Brazil (in dark gray), from the State of Alagoas with its mesoregions, Coastal/Humid Zones (1), Hinterland (2) and Interior (3) and municipality of Maceió (highlighted).

Alagoas is subdivided into three geographic mesoregions (Coastal/Humid Zones, Hinterland and Interior) with different types of climate and further details are described in Santos et al. (2014Santos, C.M.; Souza, J.L.; Ferreira Junior, R.A.; Tiba, C.; Melo, R.O. et al. On modeling global solar irradiation using air temperature for Alagoas State, Northeastern Brazil. Energy, v. 71, p. 388-398, 2014.), Santos et al. (2016)Santos, C.M.; Escobedo, J.F.; Teramoto, E.T.; Silva, S.H.M.G. Assessment of ANN and SVM models for estimating normal direct irradiation (Hb). Energy Conversion and Management, v. 126, p. 826-836, 2016. and Barros et al. (2012Barros, A.H.C.; Araujo Filho, J.C.; Silva, A.B.; Santiago, G.A.C.F. Climatologia do Estado de Alagoas. Recife: EMBRAPA SOLOS, 32p., 2012.). The microregion of Maceió is part of mesoregions of Coastal/Humid Zones. This mesoregion have hot and humid climate, with annual rainfall varying from 600 to 3000 mm, the rainy season is concentrated between autumn and winter (from May to August), with 55.2% of the annual rainfall and the dry period occurs between spring and summer (November to February) with only 12.5% of the total rainfall (Lima, 1991Lima, M.C. Variabilidade da Precipitação no Litoral Leste da Região Nordeste do Brasil. Dissertação de Mestrado em Meteorologia, INPE - 5283 - TDI/453, 1991.; Souza and Lima, 1995Souza, J.L.; Lima, F.Z. Clima - Estudo de Impactos Ambiental, Levantamento Ambiental e Relatório, no Meio Ambiente de área de Proteção Ambiental de Piaçabuçu - AL. In: Relatório Técnico, Maceió: UFAL, p. 180-238, 1995.; Molion and Bernardo, 2002Molion, L.C.B.; Bernardo, S.O. Uma revisão da dinâmica das chuvas no nordeste brasileiro. Revista Brasileira de Meteorologia, v. 17, n. 1, p. 1-10, 2002.). The rainfall regime of the State of Alagoas is associated with several meteorological systems common to much NEB. Among the main systems, the coastal instability lines, seasonal variation of trade winds, position and intensity of the Intertropical Convergence Zone (ITCZ), Upper-Level Cyclonic Vortexes (ULCV), Eastern Waves (EW), are highlighted as the main ones (Souza et al., 2004Souza, J.L.; Moura Filho, G.; Fonseca, R.F.; Teodoro, I.; Santos, E.A. et al. Análise da precipitação pluvial e temperatura do ar na região do Tabuleiro Costeiro de Maceió, AL, período 1972 - 2001. Revista Brasileira de Agrometeorologia, v. 11, n. 2, p.131-141, 2004.; Lyra et al., 2014Lyra, G.B.; Oliveira Júnior, J.F.; Zeri, M. Cluster analysis applied to the spatial and temporal variability of monthly rainfall in Alagoas state, Northeast of Brazil. Internat Journ Climat, v. 34, n. 13, p. 3546-3558, 2014.). In general, Alagoas has rainy season between May and July, when there is the maximum convergence of trade winds with the terrestrial breeze, which should be stronger during the autumn and winter seasons when the temperature contrast between land and sea is higher (Kousky, 1980Kousky, V.E. Diurnal raifall variation on Northeast Brazil. Monthly Weather Review, v. 108, n. 4, p. 488-498, 1980.), and a dry season between November and January, a period in which a shift in the direction of trade winds takes place, which impairs the formation of clouds in the region.

2.2. Obtaining data

In this study measurements of global solar irradiance (IG) and direct normal solar irradiance (ID) were obtained from 2011 to 2014. Due to technical/operational problems, the ID database used is not continuous over the measurement period. The lack of electricity and non-linearization of the sensor with the solar disk are some factors responsible for significantly reducing the database. The IG component was measured with a piranometer model CM5 Kipp & Zonen (Delft, The Netherlands) (Fig. 2). The ID component was measured by sNIP model secondary standard pyrheliometer assembled on an automatic solar tracker (Fig. 2) manufactured by Eppley. All sensors used followed the maintenance and calibration procedure according to manufacturer's recommendations (periodic calibration). Periodic checks throughout the measurements were made to ensure that the radiometer domes were clean, obtaining accurate and reliable measurements.

Figure 2
Solarimetric station: Kipp & Zonen CM5 pyranometer and pyrheliometer assembled on the Eppley solar tracker.

The radiometers used were connected to a Micrologger 21XL (Campbell Scientific Inc), with frequency of 1 Hz and programmed to take readings every 10 s and store the average every 10 min. Hourly and daily data are calculated from the average values of 10 min. The collected IG and ID data were integrated to produce global solar irradiation (HG) and direct normal solar irradiation (HD) on hourly (HGh and HDh) and daily (HGd and HDd) basis, respectively. A quality control process is essential before using any database of the place of interest. The quality of individual measures was verified according to recommendations of the World Meteorological Organization (WMO, 2008WMO. Guide to Meteorological Instruments and Methods of Observation. WMO-n°8, Seventh Edition, 157-197, Geneva, Switzerland, 2008.). http://www.soda-pro.com/home
http://www.soda-pro.com/home...

2.3. Statistical and climatic analysis of measurements

The measures of IG (HG) and ID (HD) in the hourly and daily time partitions were analyzed considering seasonality and local cloud cover. For statistical and climatic analysis of measurements, solar irradiance components were separated into periods and characterized in the different seasons of the year according to the region climatology, which is strongly dependent on the seasonality of precipitation and air temperature (Souza et al., 2004Souza, J.L.; Moura Filho, G.; Fonseca, R.F.; Teodoro, I.; Santos, E.A. et al. Análise da precipitação pluvial e temperatura do ar na região do Tabuleiro Costeiro de Maceió, AL, período 1972 - 2001. Revista Brasileira de Agrometeorologia, v. 11, n. 2, p.131-141, 2004.). In the Maceió region, the rainy season begins in April and ends in August, while the dry season begins in October and extends through February. Based on this seasonality and for further study, the climatic periods were separated into four: 1 - the dry (less rainy) period was separated into Dry Season (DS) (November to February); 2 - the transition period between dry and rainy season (March and April) was separated and defined as Pre-Rainy Season (PRS); 3 - the rainy season (May to August), will be called here as the Rainy Season (RS); 4 - the transition period between the rainy and dry season (September and October) was denominated as Pre-Dry Season (PDS).

The cloud cover was characterized from clearness index (Kt) values in the hourly (Kth) and daily (Ktd) partitions, determined by the Kt= IG/I0 ratio. Thus, based on Iqbal (1983)Iqbal, M. An Introduction to Solar Radiation. New York: Academic Press, 390 p., 1983., some adaptations were made regarding Kt intervals to the atmosphere of Maceió, when (Kt ≥ 0.60), the conditions are of Clear Sky (CS); when (0.40 ≤ Kt < 0.60) characterizes Partially Cloudy Sky (PCS) condition and for (Kt < 0.40) it characterizes Overcast Sky (OS). For the representation of irradiance curves (irradiations) and estimates, days representing the different cloud cover conditions in each of the stations were randomly chosen. The IG and ID variations were elaborated with average measures of 10 min, and hourly curves of HGh and HDh were plotted from the integration of the irradiance curves for each hour. Irradiance and irradiation analyses were performed between the interval between 06:00 am and 05:00 pm, comprising practically the entire daytime period of solar incidence.

The average value of solar irradiance (global and direct) were separated into three different cloud cover conditions and four-time intervals. The average hourly of global solar irradiance (I¯Gh) with the respective standard deviations (δIGh) and the average hourly of direct normal solar irradiance (I¯Dh) with the respective standard deviations (δIDh) were calculated at each time interval. Histograms of direct normal daily solar irradiance (IDd) frequencies of each season of the year under study and under different cloud cover conditions were performed with the purpose of analyzing seasonal tendencies.

2.4. Calculation of solar irradiance at the top of the atmosphere

The solar irradiance at the top of the atmosphere (I0) was calculated by Eq. (1):

(1) I 0 = I S C E 0 c o s θ Z

in which ISC is the solar constant (1367 W m−2), E0 is the correction factor of the eccentricity of the, further details are described in Souza et al. (2005Souza, J.L.; Nicácio, R.M.; Moura, M.A.L. Global solar radiation measurements in Maceió, Brazil. Renewable Energy, v. 30, n. 8, p. 1203-1220, 2005.) and Iqbal (1983)Iqbal, M. An Introduction to Solar Radiation. New York: Academic Press, 390 p., 1983.. I0 data were integrated to obtain the solar irradiation values at the top of the atmosphere (H0), hourly (H0h) and daily (H0d).

3. Results and Discussion

3.1. Climate and Weather system

The climate of Maceió is “As” type, characterized as rainy tropical, with rain in the winter and autumn and dry in the summer (Köppen classification). The air temperature is high (> 20.0 °C), with monthly thermal amplitude equal to 2.9 °C, average air temperature of 23.9 °C in the coldest month (July) and 26.6 °C in the hottest (March), the average is 25.4°C (Fig. 3a). In spring and summer, it is common to have days with maximum air temperature above 30.0 °C, whereas in autumn and winter, daily air temperatures rarely reduce values below 20.0 °C. The relative air humidity presents low annual fluctuation, with monthly average values greater than 75% (Fig. 3a), which is due to its proximity to the Atlantic Ocean. Rainfall in the wettest month is 341.9 mm (June) and in the least rainy month 30.2 mm (December); summer precipitation is approximately 293.4 mm (Fig. 3b). Spring is the season with the lowest precipitation (168.8 mm).

Figura 3
Monthly average variation: a) Air temperature (T, °C) and Relative humidity (RH,%); b) Precipitation (P, mm) and sunshine duration (SD, hours).

From November to March, precipitation is usually attributed to the convergence between the land breeze and the southeastern trade winds; however, in this time of year, this does not occur often, since a shift in the direction of trade winds decreases the intensity of the encounter with the terrestrial breeze, with lower contrast between the surface and the ocean temperatures. Another meteorological phenomenon contributing to the precipitation during this time of year is the ULCV, characterized as closed cyclone circulation with subsiding air in the center, and in these systems, the center has Clear Sky; however, its peripheries are submitted to large amounts of precipitation (Molion and Bernardo, 2002Molion, L.C.B.; Bernardo, S.O. Uma revisão da dinâmica das chuvas no nordeste brasileiro. Revista Brasileira de Meteorologia, v. 17, n. 1, p. 1-10, 2002.; Reboita et al., 2012Reboita, M.S.; Krusche, N.; Ambrizzi, T.; Porfírio, R. Entendendo o tempo e o clima na América do Sul. Terrae Didatica, v. 8, n. 1, p.34-50, 2012.). The precipitation reaches 820.2 mm and 527.2 mm in the autumn and winter, respectively. In the rainy season (April to August) precipitation is caused due to the action of some weather systems dependent on the effect of some surface and atmosphere configurations. We can highlight the occurrence of the penetration of frontal systems that may in some cases reach equatorial latitudes (Kousky, 1979Kousky, V.E. Frontal influences on Northeast Brazil. Monthly Weather Review, v. 107, n. 9, p. 1140-1153, 1979.), the convergence of trade winds with the terrestrial breeze that reaches its peak during the autumn and winter seasons with the increasing contrast of temperature between earth and ocean (Kousky, 1980Kousky, V.E. Diurnal raifall variation on Northeast Brazil. Monthly Weather Review, v. 108, n. 4, p. 488-498, 1980.), the eastern undulating disturbances known as eastern waves (EW) (Chou, 1990Chou, S.C. Análise de Distúrbios Ondulatórios de Leste Sobre o Oceano Atlântico Equatorial Sul. Dissertação de Mestrado em Meteorologia, Instituto Nacional de Pesquisas Espaciais, São José dos Campos, 134 p., 1990.; Espinoza, 1996Espinoza, E.S. Distúrbios nos Ventos de Leste no Atlântico Tropical. São José dos Campos: INPE, 149p., 1996.), and finally the action of the South Atlantic Convergence Zone (SACZ) when moved to lower latitudes (near Equator). The SACZ configuration is related to the ITCZ positioning, causing good part of the precipitation in Maceió (Molion and Bernardo, 2000Molion, L.C.B.; Bernardo, S.O. Dinâmica das chuvas sobre o Nordeste do Brasil. In: Congresso Brasileiro de Meteorologia, 11, Rio de Janeiro, p. 1334-1342, 2000.). The month of December has the longest sunshine duration (SD) (260.4 hours), followed by November (259.2 hours) (Fig. 3b). The months of June (164.7 hours) and July (165.4 hours) are those with the lowest SD.

In addition to its proximity to the Atlantic Ocean, Maceió is also influenced by the presence of two large lagoons: Mundaú and Manguaba. The Mundaú lagoon has an approximate reservoir of 23 km2 and Manguaba with 42 km2. Both are located within a radius of 15 km south of Maceió. Their dimensions represent a local microclimate, causing interference in the atmosphere and in the microphysics of clouds. With the convection of the ocean and lagoons, the region of Maceió is strongly influenced by water vapor for most of the year (Fig. 4), helping in the formation and predominance of cumulus clouds in the local atmosphere. Water vapor accompanies the annual trend of air temperature and SD, increasing values in spring - summer and decreasing in autumn - winter.

Figura 4
Monthly average of water vapor in the atmosphere of Maceió.

3.2. Analysis of solar irradiance under clear sky conditions (CS)

Regardless of astronomical effect, on days whose sky has high atmospheric transmittance, the daytime solar irradiance variation is attributed only to atmospheric constituents (aerosols and water vapor, for example), mainly in the ID magnitude (Gueymard, 2012Gueymard, C.A. Temporal variability in direct and global irradiance at various time scales as affected by aerosols. Solar Energy, v. 86, p. 3544-3553, 2012.). Aerosols can reach levels that attenuate ID by approximately 70%, or more in extreme conditions (Ruiz-Arias et al., 2016Ruiz-Arias, J.A.; Gueymard, C.A.; Santos-Alamillos, F.J.; Pozo-Vázquez, D. D. Worldwide impact of aerosol's time scale on the predicted long-term concentrating solar power potential. Scientific Reports, 6, 30546, 2016.). Under CS, the daytime radiation components did not have large oscillations for DS, with total: HGd = 27.52 MJ m−2 and HDd = 31.33 MJ m−2, morning values (HG = 16.34 MJ m−2 and HD = 18.41 MJ m−2) and afternoon values (HG = 11.18 MJ m−2 and HD = 12.92 MJ m−2) (Fig. 5a). HD values in the morning (from sunrise to noon) and afternoon (from noon to sunset) are higher than HG values for the same period. This can be caused by the atmospheric optical mass and the method of measuring the solar components that are different. While ID (HD) is measured directly (with the pyrheliometer) pointed at the solar disk (following diurnal solar trajectory), IG (HG) is measured with fixed sensor (pyranometer) on the horizontal surface. The daily average IG and ID were (687.2 ± 292.0 W m−2) and (787.0 ± 148.6 W m−2), respectively. The maximum IG value is close to that observed in Maceió (Porfírio et al., 2012Porfírio, A.C.S.; Souza, J.L.; Lyra, G.B.; Lemes, M.A.M. An assessment of the global UV solar radiation under various sky conditions in Maceió-Northeastern Brazil. Energy, v. 4 4, n. 1, p. 584-592, 2012.) and higher than that found by Podstawczynska (2010)Podstawczynska, A. UV and global solar radiation in Lodz, central Poland. International Journal Climatology, v. 30, n. 1, p. 1-10, 2010.. For six climatic zones in Nigeria and considering all types of sky cover condition, Chukwujindu and Julie (2017)Chukwujindu, N.S.; Julie, O. Estimation of direct normal irradiance under various sky conditions in data sparse tropical ecological zones in Nigeria. Internacional Journal of Advanced Astronomy, v. 5, n. 2, p. 90-105, 2017. found average values ranging from 6.30 to 28.8 MJm−2.

Figura 5
Variation of global (IG) and direct normal (ID) solar irradiances under clear sky condition: (a) (12/12/2013) in dry season; (b) (03/17/2011) in pre-rainy season; (c) (06/24/2011) in rainy season; (d) (10/09/2012) in pre-dry season.

In PRS, under CS conditions, daytime variation has a marked fluctuation caused by the passage of dispersed clouds in the atmosphere (Fig. 5b): with daily totals for HGd = 26.99 MJ m−2 and HDd = 28.07 MJ m−2. Although it is characterized as clear sky day, the daily formation of transient clouds of the cumulus type is common at this time of year due to the proximity of the region to the Atlantic Ocean, which in some cases may develop into other types of clouds, both of horizontal or vertical development. The maximum IG and ID values were 1107.0 W m−2 and 905.8 W m−2, respectively. The daily average was IG = 672.1 ± 329.5 W m−2 and ID = 702.4 ± 183.6 W m−2. During this period of year, negative variation may occur due to the decrease of the solar radiation available at the top of the atmosphere and the increase of the water vapor content in the atmosphere, which is absorbed with greater intensity. Similar values for IG and ID irradiance was found for different sites in Portugal (Cavaco et al., 2016Cavaco, A.; Canhoto, P.; Costa, M.J.; Collares-Pereira, M. DNI measurements in the South of Portugal: Long term results through direct comparison with global and diffuse radiation measurements and existing time series. AIP Conference Proceedings 1734, 150001, 2016.).

Also under CS conditions, for RS, the existence of greater oscillations and reduction of diurnal IG and ID values can be observed (Fig. 5c): with daily totals equal 20.90 MJ m−2 to HG and 24.02 MJ m−2 to HDd The maximum IG and ID values were 870.3 W m−2 and 824.0 W m−2, with daily average of 529.8 ± 279.4 W m−2 and 632.0 ± 224.7 W m−2, respectively. During this period, there is a decrease in the available solar radiation due to the greater solar declination in the Northern Hemisphere, increase in cloudiness and reduction in the water vapor content in the local atmosphere. In addition, it is also verified that between 10:00 am and 01:00 pm, irradiance values are approximately constant, varying from (812.6 to 870.3 W m−2) for IG and from (815.2 to 816.3 W m−2) for ID.

Greater variations in values occurred in PDS during the daytime period, in which a sharp drop in the irradiance values was observed when IG decreased from (627.3 to 385.1 W m−2), and ID from (738.6 to 175.4 W m−2) (Fig. 5d). This oscillation is the result of the dynamics of cumulus clouds in the local atmosphere, because in this period (PDS), the water vapor in the atmosphere remains approximately constant. The maximum values observed for the day were IG = 1149.1 W m−2 and ID = 923.2 W m−2, with daily totals HGd = 26.50 MJ m−2 and HDd = 24.52 MJ m−2. It was observed that ID increases with the elevation of IG, revealing that ID is controlled by the variation of IG. Results similar to those from Maceió for ID were obtained by Padova (northern Italy) and Trisaia (Southern Italy) (Dugaria et al., 2015Dugaria, S.; Padovan, A.; Sabatelli, V.; Del Col, D. Assessment of estimation methods of DNI resource in solar concentrating systems. Solar Energy, v. 121, p. 103-115, 2015.).

3.3. Analysis of solar irradiance under partially cloudy sky conditions (PCS)

Under PCS condition, the irradiance curves have larger fluctuations throughout the day, caused mainly by the dynamics of cumulus clouds in the atmosphere (Fig. 6a-d). The maximum peaks IG and ID values were higher under PCS condition than on CS days. This is explained by the multi-reflections between the surface and the cloud bases, being phenomena typical of partially cloudy days, causing an increase of the descending reflection that increases the irradiance in comparison with cloud-free days (Porfírio et al., 2012Porfírio, A.C.S.; Souza, J.L.; Lyra, G.B.; Lemes, M.A.M. An assessment of the global UV solar radiation under various sky conditions in Maceió-Northeastern Brazil. Energy, v. 4 4, n. 1, p. 584-592, 2012.). Souza et al. (2005Souza, J.L.; Nicácio, R.M.; Moura, M.A.L. Global solar radiation measurements in Maceió, Brazil. Renewable Energy, v. 30, n. 8, p. 1203-1220, 2005.) in a study on the solar irradiance for Maceió observed the occurrence of multiple reflections, confirmed by IG values very close to I0. In different places in the State of Alagoas / Brazil, it has been reported that this phenomenon is not uncommon and occurs on at least one third of the days in a given month (Andrade and Tiba, 2016Andrade, R.C.; Tiba, C. Extreme global solar irradiance due to cloud enhancement in northeastern Brazil. Renewable Energy, v. 86, p. 1433-1441, 2016.).

Figura 6
Variation of global (IG) and direct normal (ID) solar irradiances under partially cloudy sky condition: (a) in dry season (01/09/2014), (b) in pre-rainy season (03/02/2013), (c) in rainy season (06/07/2012) and (d) in pre-dry season (09/17/2012).

As expected, the HGd (20.52 MJ m−2) and HDd values (18.19 MJ m−2) under PCS condition for DS were lower than those obtained under CS condition (Fig. 6a). The maximum IG and ID values were 1123.0 W m−2 and 929.8 W m−2, respectively, with average daily of IG = 512.2 ± 292.6 W m−2 and ID = 458.3 ± 314.3 W m−2. The HGd values under these conditions were close to those observed by Porfírio et al. (2012)Porfírio, A.C.S.; Souza, J.L.; Lyra, G.B.; Lemes, M.A.M. An assessment of the global UV solar radiation under various sky conditions in Maceió-Northeastern Brazil. Energy, v. 4 4, n. 1, p. 584-592, 2012.. In the PRS, the irradiance that reached the surface was lower than in PDS, with daily average IG and ID values (maximum) of 507.6 ± 278.3 W m−2 (1099.6 W m−2) and 306.4 ± 241.6 W m −2 (859.0 W m−2), respectively (Fig. 6b). The irradiation values were 20.05 and 11.58 MJ m−2 for HGd and HDd, respectively. Similar values were observed by Polo et al. (2017Polo, J.; Fernández-Peruchena, C.; Gastón, M. Analysis on the long-term relationship between DNI and CSP yield production for different technologies. Solar Energy, v. 155, p. 1121-1129, 2017.).

The solar irradiance values in RS were smaller when compared to the other seasons, with a sharp decrease (between 09:00 am and 12:00 am) in IG from (673.7 to 302.2 W m−2), while ID decreases from (500.7 to 18.8 W m−2) (Fig. 6c). The irradiance remained low and increased to the maximum peak of IG (1015.5 W m−2) and ID (827.5 W m−2). The average daily irradiance values were IG = 378.2 ± 244.5 W m−2 and ID = 250.24 ± 257.60 W m−2. During this period, cloudiness is the main responsible for the irradiance variations throughout the day. Irradiation values were 14.62 MJ m−2 and 9.49 MJ m−2 for HGd and HDd, respectively. The average PRS values [IG = 480.5 ± 288.0 W m−2 (HGd = 19.29 MJ m−2) and ID = 272.2 ± 238.2 W m−2 (HDd = 10.94 MJ m−2)] are higher than those observed in RS and smaller than those observed in PRS and DS (Fig. 6d). In four sites in Australia, Grantham et al. (2017Grantham, A.P; Pudney, P.J.; Ward, L.A.; Belusko, M.; Boland, J.W. Generating synthetic five-minute solar irradiance values from hourly observations. Solar Energy, v. 147, p. 209-221, 2017.) found diurnal curves for IG and ID under CS and OS conditions with the same behavior of results of this research. On day of OS, evidences of the transient effects of clouds causing fluctuations in ID are observed. On day of CS, IG and ID do not exceed the value of 1000 W m−2. In Maceió, values above this value are observed.

3.4. Analysis of solar irradiance under overcast sky conditions (OS)

Under OS condition, low IG and ID values are expected, since much of the solar irradiance cannot cross the Earth's atmosphere, being absorbed, reflected or scattered by clouds. In this condition, the clouds result in a reduction in ID values. Therefore, the absence of clouds causes the elevation of IG and ID in CS days, while in OS days, ID tends to zero and the IG is significantly reduced to smaller values (IGId). Thus, under OS condition, the daily IG values are greater than the daily ID values. Studies have shown that ID is much more sensitive than IG for changes in the amounts of clouds and aerosols (Ruiz-Arias et al., 2016Ruiz-Arias, J.A.; Gueymard, C.A.; Santos-Alamillos, F.J.; Pozo-Vázquez, D. D. Worldwide impact of aerosol's time scale on the predicted long-term concentrating solar power potential. Scientific Reports, 6, 30546, 2016.). Moreover, in this sky cover, it is expected that the irradiance values for DS [IG = 362.4 ± 161.3 W m−2 and ID = 123.7 ± 188.2 W m−2] are lower than those observed in the other conditions (Fig. 7a). Peak IG and ID values were 732.7 W m−2 and 625.3 W m−2, respectively. It was observed that the morning, evening and daily HG totals were one of the lowest values found, especially during the morning, whose ID was almost non-existent, caused by cloud attenuation.

Figura 7
Variation of global (IG) and direct normal (ID) solar irradiances under overcast sky condition: (a) (12/30/2012) in dry season, (b) (04/07/2014) in pre-rainy season, (c) (08/05/2013) in rainy season and (d) (10/01/2012) in pre-dry season.

In PRS, the average daily IG and ID values were equal to 335.5 ± 200.3 W m−2 and 208.3 ± 226.6 W m−2, respectively (Fig. 7b). The daily HG and HD totals were equal to 13.46 and 8.24 MJ m−2, respectively. In RS, ID had values close to zero during the daytime period, the daily averages were IG = 265.0 ± 166.6 W m−2 (HGd= 10.47 MJ m−2) and ID = 49.6 ± 107.5 W m−2 (HDd = 1.94 MJ m−2), with daily maximum values of IG = 801.9 W m−2 and ID = 603.0 W m−2 (Fig. 7c). The daily totals were the lowest, since the selected day belongs to the rainy season and has low Kt value. In relation to OS days in RS (Souza et al., 2005Souza, J.L.; Nicácio, R.M.; Moura, M.A.L. Global solar radiation measurements in Maceió, Brazil. Renewable Energy, v. 30, n. 8, p. 1203-1220, 2005.; Porifírio et al., 2012Porfírio, A.C.S.; Souza, J.L.; Lyra, G.B.; Lemes, M.A.M. An assessment of the global UV solar radiation under various sky conditions in Maceió-Northeastern Brazil. Energy, v. 4 4, n. 1, p. 584-592, 2012.) found similar values with measurements in the same location from different years. The average daily IG and ID values were 328.6 ± 243.3 W m−2 (HGd = 13.18 MJ m−2) and 61.2 ± 126.3 W m−2 (HDd = 2.38 MJ m−2) and daily peaks reached 1089.5 W m−2 and 723.7 W m−2, respectively (Fig. 7d). Study IG and IDSouza et al. (2005)Souza, J.L.; Nicácio, R.M.; Moura, M.A.L. Global solar radiation measurements in Maceió, Brazil. Renewable Energy, v. 30, n. 8, p. 1203-1220, 2005., Vilani et al. (2006Vilani, M.T.; Sanches, L.; Nogueiram, J.S.; Priante Filho, N. Sazonalidade da radiação, temperatura, e umidade em uma floresta de transição Amazônia Cerrado. Revista Brasileira de Meteorologia.v. 21, n. 3b, p. 331-343, 2006.) and Porifírio et al. (2012)Porfírio, A.C.S.; Souza, J.L.; Lyra, G.B.; Lemes, M.A.M. An assessment of the global UV solar radiation under various sky conditions in Maceió-Northeastern Brazil. Energy, v. 4 4, n. 1, p. 584-592, 2012. found values with the same quantitative characteristics regarding the periods of greatest and smallest oscillations. The influence of cloudiness on IG and ID levels is remarkable. In addition, atmospheric water vapor has strong absorption in the region of the infrared spectrum, being able to reduce significantly the values of irradiances incident on the surface. The small differences observed along the analyzed periods show that energy totals are more affected by cloudiness than by seasonal variation, when compared with days of CS.

3.5. Hourly average of the global (I¯Gh) and direct normal (I¯Dh) solar irradiance

In DS under CS condition, the maximum I¯Gh value was 765.6 ± 100.0 W m−2 occurred between 09:00 am and 12:00 pm and was 728.2 ± 29.9 W m−2 (12:00 pm - 03:00 pm), which were the highest values found for CS (Table 1). Under PCS condition, the maximum I¯Gh and I¯Dh values were 608.6 ± 91.6 W m−2 and 420.5 ± 26.8 W m−2, respectively. Under OS condition, I¯Gh had maximum value of 354.9 ± 49.2 W m−2 and I¯Dh with maximum value of 197.3 ± 16.7 W m−2. The lowest values are found for OS condition and the highest values for the CS condition, evidencing the effect of clouds and atmospheric constituents on the attenuation of irradiances.

Table 1
Hourly average of the global (I¯Gh, W m−2) and direct normal (I¯Dh, W m−2) solar irradiances under different cloud cover for the dry season (November-February) with the respective standard deviations (δIGh and δIDh).

In PRS under CS condition, the maximum I¯Gh and I¯Dh values were 829.3 ± 105.0 W m−2 and 739.9 ± 31.8 W m−2, which represents an increase of 8.31 and 1.59%, when compared to DS values (Table 2). For the PCS condition, the highest global and direct irradiance values were 674.2 ± 97.8 W m−2 and 425.2 ± 18.9 W m−2, respectively, indicating an increase of 10.76% in I¯Gh and 1.10% the I¯Dh in relation to ES. Under OS condition, the maximum I¯Gh and I¯Dh values were 367.2 ± 40.2 W m−2 and 196.1 ± 3.6 W m−2, respectively, an increase of 6.47% and decrease of 0.60% in comparison to DS. The lowest I¯Dh values are observed in PRS under OS condition. As RS is prolonged, there is a relatively higher presence of clouds, increased humidity and absorption of near infrared (NIR) radiation of the solar spectrum, which produces low I¯Dh magnitude.

Table 2
Hourly average of the global (I¯Gh, W m−2) and direct normal (I¯Dh, W m−2) solar irradiances under different cloud cover for the pre-rainy season (March and April) with the respective standard deviations (δIGh and δIDh).

Under CS condition for RS, irradiances had maximum I¯Gh values = 712.2 ± 103.4 W m−2 and I¯Dh values = 649.0 ± 37.7 W m−2, reduction of ≈ 16.42% and ≈14% in relation to PRS, respectively (Table 3). Under PCS condition, the maximum I¯Gh value was 572.6 ± 65.7 W m−2 (decrease of ≈17.73%) and I¯Dh was 347.7 ± 10.5 W m−2 (reduction of ≈ 22.28%) compared to PRS. Under OS condition, I¯Gh has maximum value of 354.0 ± 45.6 W m−2, representing a reduction of ≈ 3.74% when compared to PRS; I¯Dh reached maximum value of 169.1 ± 59.7 W m−2, with a decrease of ≈ 15.98% compared to PRS.

Table 3
Hourly average of the global (I¯Gh, W m−2) and direct normal (I¯Dh, W m−2) solar irradiances under different cloud cover for the rainy season (May-August) with the respective standard deviations (δIGh and δIDh).

The maximum I¯Gh and I¯Dh values under CS condition for PDS were 887.8 ± 103.7 W m−2 and 657.3 ± 46.7 W m−2, respectively, increasing ≈ 24.64% and ≈ 1.27% when compared to RS (Table 4). Under PCS conditions, the maximum I¯Gh and I¯Dh values were 690.9 ± 103.2 W m−2 and 406.9 ± 21.6 W m−2, increase of 20.66% and 17.01% with respect to RS. Under OS condition, he maximum I¯Gh and I¯Dh values were 442.9 ± 70.2 W m−2 and 532.0 ± 35.8 W m−2, increase of ≈ 25.12% for I¯Gh and ≈ 314.66% for I¯Dh.

Table 4
Hourly average of the global (I¯Gh, W m−2) and direct normal (I¯Dh, W m−2) solar irradiances under different cloud cover for the pre-dry season (September and October) with the respective standard deviations (δIGh and δIDh).

3.6. Hourly solar irradiation

Figure 8 shows the average hourly global solar irradiation (H¯Gh) observed for all four seasons (obtained through the average irradiation of each hour interval). The diurnal variation among the different seasons is the same, increasing from morning until reaching maximum near noon and then falling until the end of the afternoon. In relation to this behavior, the hourly values tend to spread symmetrically to the solar half day. The largest curve amplitude is in DS, followed by PDS, PRS and RS. This behavior evidences the solar declination effect among seasons, the amplitude decreases as the day of the year distances from the summer solstice. The diurnal evolution is characterized by the effect of the zenith angle, since H¯Gh is measured in the horizontal plane. In DS, the total HG value in the day is 22.47 MJ m−2, and the morning and afternoon values were 12.01 and 10.47 MJ m −2, respectively. The diurnal evolution in PRS has a decrease in values when compared to DS, HGd was 20.45 MJ m−2, with maximum hourly value (H¯Gh = 2.90 MJ m−2). In RS, the hourly averages and the daily total (HGd = 15.45 MJ m−2) had the greatest reduction, and the maximum value at this season was H¯Gh = 2.24 MJ m−2. The diurnal evolution in PDS for HGd = 20.45 MJ m−2 and maximum hourly H¯Gh = 2.91 MJ m−2. Approximate H¯Gh values were also found by researchers Souza et al. (2005)Souza, J.L.; Nicácio, R.M.; Moura, M.A.L. Global solar radiation measurements in Maceió, Brazil. Renewable Energy, v. 30, n. 8, p. 1203-1220, 2005., Marques Filho et al. (2016)Marques Filho, E.P.; Oliveira, A.P.; Vita, W.A. et al. Global, diffuse and direct solar radiation at the surface in the city of Rio de Janeiro: Observational characterization and empirical modeling. Renewable Energy, v. 91, p. 64-74, 2016. and Tahir and Asim (2018)Tahir, Z.R.; Asim, M. Surface measured solar radiation data and solar energy resource assessment of Pakistan: A review. Renewable and Sustainable Energy Reviews, v. 81, Parte 2, p. 2839-2861, 2018.. Although they are not located at the same latitude, the results of Maceió are similar to those obtained at sites in Nigeria by Osinowo et al. (2015)Osinowo, A.A.; Okogbue, E.C.; Ogungbenro, S.B.; Fashanu, O. Analysis of Global Solar Irradiance over Climatic Zones in Nigeria for Solar Energy Applications. Journal of Solar Energy, Volume 2015, Article ID 819307, 9 p., which found a uniform pattern in the distribution of H¯Gh values, with peaks ranging from (15.01 to 12.01) MJ m−2 and with results observed in Saudi Arabia (AlYahya and Irfan, 2016AlYahya, S.; Irfan, M.A. Analysis from the new solar radiation Atlas for Saudi Arabia. Solar Energy, v. 130, p. 116-127, 2016.).

The daytime variation of H¯Gh has a decrease after 08:00 am returning to increase soon after. This is due to the greater cloudiness in the morning, caused by the effect of local wind circulation (sea breezes), which are responsible for creating a zone of convergence near the coast (front breeze), contributing to the formation of cumulus clouds until near 01:00 pm, when it is verified the formation of an area of suppression of cumulus clouds related to the displacement of the zones of breeze subsidence and convergence (National Research Council, 1992National Research Council. 1992. Coastal Meteorology: A Review of the State of the Science. Washington: The National Academy Press, 1992. 112 p.; Silva, 2003Silva, G.R. Características de Vento na Região Nordeste: Análise, Modelagem, e Aplicações Para Projetos de Centrais Eólicas. Dissertação de Mestrado em Engenharia Mecânica, Universidade Federal de Pernambuco, Recife, 141 p., 2003.). This effect is more intense in HD because this component is measured with the sensor directly pointed at the solar disk. In DS, HD has value of 17.20 MJ m−2 and maximum hourly average direct solar normal (H¯Dh) value = 2.00 MJ m−2 (Fig. 8b). In PRS, the HD values were higher than in DS, with HD = 18.38 MJ m−2 and maximum H¯Dh = 2.23 MJ m−2. H¯Dh has the lowest values in RS, with = 12.48 MJ m−2 and maximum peak H¯Dh = 1.46 MJ m−2. In PDS, HD was 14.88 MJ m−2 and maximum H¯Dh = 1.90 MJ m−2. The seasonal variation in the H¯Dh transmission is associated with sky cover conditions: large (small) Kt values are related to large (small) values of the H¯Dh transmission. The local results resemble those obtained for two sites in Cyprus by Pashiardis et al. (2017)Pashiardis, S.; Kalogirou, S.A.; Pelengaris, A. Statistical analysis for the characterization of solar energy utilization and inter-comparison of solar radiation at two sites in Cyprus. Applied Energy, v. 190, p. 1138-1158, 2017., who performed a statistical analysis of HD values, results found by Marques Filho et al. (2016)Marques Filho, E.P.; Oliveira, A.P.; Vita, W.A. et al. Global, diffuse and direct solar radiation at the surface in the city of Rio de Janeiro: Observational characterization and empirical modeling. Renewable Energy, v. 91, p. 64-74, 2016. for Rio de Janeiro and for locations in Kenya (Kariuki and Sato, 2018Kariuki, B.W.; Sato, T. Interannual and spatial variability of solar radiation energy potential in Kenya using Meteosat satellite. Renewable Energy, v. 116, Parte A, p. 88-96, 2018.), with values from satellite and two places of Morocco (Gairaa et al., 2017Gairaa, K.; Khellaf, A.; Benkacilai, S.; Guermoui, M. Solar radiation measurements in Algeria: case of Ghardaïa station as member of the enerMENA meteorological network. Energy Procedia, v. 141, p. 50-54, 2017.).

Figura 8
The hourly average solar irradiation: a) global (H¯Gh); b) direct solar normal (H¯Dh). dry season (DS), pre-rainy season (PRS), rainy season (RS) and pre-dry season (PDS).

3.7. Frequency distribution of daily direct normal solar irradiance (IDd) in the different sky cover conditions

In the frequency distribution of IDd for CS condition in the four seasons, it was observed that the IDd values in DS had greater occurrence between 400 and 600 W m−2, comprising a relative occurrence rate of 68.8% of days, with higher occurrence of IDd for 440 W m−2 (Fig. 9a). Between 160 and 400 W m−2, IDd maintained at a relative occurrence rate of 13.6% and from 600 W m−2 up to the maximum value (740 W m−2), the relative occurrence rate corresponded to 17.6%. IDd in PRS had relative occurrence rate of 65.91% between 400 and 600 W m−2, between 257.5 and 400 W m−2, the occurrence of IDd remained with relative occurrence rate of 15.69%, but from 600 to 700 W m−2, the relative occurrence rate corresponded 18.19% of days (Fig. 9b). For RS, IDd had higher relative occurrence rate between 400 and 600 W m−2, which corresponds to about 60.92% of days of this season and from 252.5 to 400 W m−2, the occurrence reached 39.08% of cases (Fig. 9c). In PDS, IDd values between 385 and 635 W m−2 corresponded to half the relative occurrence rate of all values, the other half was between 235 and 385 W m−2 and the interval with the highest relative occurrence rate was between 385 and 435 W m −2, which totaled 27.77% (Fig. 9d). For the city of Daejeon, in South Korea, Lee et al. (2017Lee, H-J.; Kim, S-Y.; Yun, C-Y. Comparison of Solar Radiation Models to Estimate Direct Normal Irradiance for Korea. Energies, v. 10, n. 5, p. 594, 2017.) showed higher relative occurrence rate for 50 W m−2, followed by 850 W m−2 and with lower relative occurrence rate for 950 W m−2.

Figura 9
Frequency distribution (F) of daily direct normal solar irradiance (IDd) under clear sky conditions: a) dry season; b) pre-rainy season; c) rainy season and d) pre-dry season. N = number of days.

Under PCS condition, IDd for DS showed higher occurrence between 250 and 450 W m−2, totaling a relative occurrence rate of 55.43%, from 10 to 250 W m−2, the relative occurrence rate was 41.30%, in the interval from 450 W m−2 up to the maximum value (IDd = 590 W m−2), the relative occurrence rate was 3.27% (Fig. 10a). In PRS, IDd had relative occurrence rate of 60.47% for the interval from 200 to 400 W m−2 (Fig. 10b). From 60 to 200 W m−2, the percentage of IDd was 19.76% and from 400 to 600 W m−2, corresponding to relative occurrence rate of 19.77%. The highest occurrence value of IDd was for the 300 W m−2. In the relative occurrence rate histograms of IDd for RS (Fig. 10c), higher occurrence was observed close to 230 W m−2, with higher relative occurrence rate (80%) between 38 and 308 W m−2. In PDS, most occurrences of IDd were between 25 and 275 W m−2, totaling a relative occurrence rate of 70.83%, between 275 and 675 W m−2 the relative occurrence rate of was 29.17% (Fig. 10d). In Seville - Spain, Moreno-Tejera et al. (2016)Moreno-Tejera, S.; Silva-Perez, M.A.; Lillo-Bravo, I.; Ramırez-Santigosa, L. L. Solar resource assessment in Seville, Spain. Statistical characterisation of solar radiation at different time resolutions. Solar Energy, v. 132, p. 430-441, 2016. found that 42.4% of ID values F are higher than 600 W m−2, with higher relative occurrence rate of high values close to the summer solstice of the Southern Hemisphere.

Figura 10
Frequency distribution (F) of daily direct normal solar irradiance (IDd) under partially cloudy sky conditions: a) dry season; b) pre-rainy season; c) rainy season and d) pre-dry season. N = number of days.

Under OS condition, the behavior of IDd is different from the other sky covers, since the higher relative occurrence rate for values smaller than 100 W m−2 as observed in DS (Fig. 11a), whose greatest occurrences were up to 188 W m−2, corresponding to a relative occurrence rate of 71.42% and in the interval from 188 to 563 W m−2, the corresponding relative occurrence rate is 28.58%. For PRS, all IDd measurements were not higher than 200 W m−2, with 55% of the relative occurrence rate up to 100 W m−2 (Fig. 11b). In RS, the maximum value was 155 W m−2, and up to 105 W m−2, the relative occurrence rate was 80.64% and from 105 to 155 W m−2, the relative occurrence rate was 19.36% (Fig. 11c). In PDS, the IDd values increased and the accumulated relative occurrence rate reached 100% up to a maximum of 488 W m−2 and 50% up to 338 W m−2 (Fig. 11d).

Figura 11
frequency distribution (F) of daily direct normal solar irradiance (IDd) under overcast sky conditions (Kt < 0.40): a) dry season; b) pre-rainy season; c) rainy season and d) pre-dry season. N = number of days.

4. Final considerations

Cumulus clouds were mainly responsible for seasonal oscillation and irradiance attenuation, significantly reducing surface values. It was observed the phenomenon of multireflections that resulted in peaks and maximum values in the condition of partly cloudy days.

The instantaneous global and direct normal solar irradiance for seasonal periods showed maximum values close to noon, especially for clear sky conditions (global atmospheric transmittance, Kt ≥ 0.60).

Higher oscillations were found in the intermediate seasonal periods (pre-rainy season - PRS and pre-dry season - PDS).

The daily totals of global and direct normal solar irradiation have greater variations due to the cloudiness than astronomical factors (seasonality), being the greatest values found, in the majority, in the morning during the dry and pre-rainy season.

The average daily totals of the global hourly solar irradiation were higher during the dry season, different from the direct totals that were higher in the pre-rainy season.

The lowest cloudiness levels were observed in the dry season, followed by pre-rainy, pre-dry and rainy seasons.

The relative occurrence rate of direct normal irradiance in the daily partition under conditions of Kt ≥ 0.60 was greater between 400 and 700 W m−2. Under conditions of 0.4 ≤ Kt <0.6, IDd occurred more between 200 and 500 W m−2 and under conditions of Kt < 0.4, it was below 200 W m−2.

Therefore, according to the obtained results, it could be inferred that the patterns of solar incidence in the region of study are dependent on the cloudiness conditions. On the other hand, seasonality is a secondary factor that has little influence on the configuration of energy patterns. Thus, it was possible to characterize the distribution of direct normal solar irradiance for the region of Maceió, but there is still a need to intensify these studies for better understanding and application in several human activities, notably engineering, energy, agriculture, health and environment.

Acknowledgments

The authors thank the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Centrais Elétricas Brasileiras S. A (Eletrobrás).

References

  • AlYahya, S.; Irfan, M.A. Analysis from the new solar radiation Atlas for Saudi Arabia. Solar Energy, v. 130, p. 116-127, 2016.
  • Andrade, R.C.; Tiba, C. Extreme global solar irradiance due to cloud enhancement in northeastern Brazil. Renewable Energy, v. 86, p. 1433-1441, 2016.
  • Barros, A.H.C.; Araujo Filho, J.C.; Silva, A.B.; Santiago, G.A.C.F. Climatologia do Estado de Alagoas Recife: EMBRAPA SOLOS, 32p., 2012.
  • Cavaco, A.; Canhoto, P.; Costa, M.J.; Collares-Pereira, M. DNI measurements in the South of Portugal: Long term results through direct comparison with global and diffuse radiation measurements and existing time series. AIP Conference Proceedings 1734, 150001, 2016.
  • Chou, S.C. Análise de Distúrbios Ondulatórios de Leste Sobre o Oceano Atlântico Equatorial Sul Dissertação de Mestrado em Meteorologia, Instituto Nacional de Pesquisas Espaciais, São José dos Campos, 134 p., 1990.
  • Chukwujindu, N.S.; Julie, O. Estimation of direct normal irradiance under various sky conditions in data sparse tropical ecological zones in Nigeria. Internacional Journal of Advanced Astronomy, v. 5, n. 2, p. 90-105, 2017.
  • Codato, G.; Oliveira, A.P.; Soares, J.; Escobedo, J.F.; Gomes, E.N.; Pai, A.D. D. Global and diffuse solar irradiances in urban and rural areas in southeast Brazil, Theor. Appl. Climatol., v. 93, p. 57-73, 2008.
  • Dugaria, S.; Padovan, A.; Sabatelli, V.; Del Col, D. Assessment of estimation methods of DNI resource in solar concentrating systems. Solar Energy, v. 121, p. 103-115, 2015.
  • Espinoza, E.S. Distúrbios nos Ventos de Leste no Atlântico Tropical São José dos Campos: INPE, 149p., 1996.
  • Evseev, E.G.; Kudish, A.I. Analysis of solar irradiation measurements at Beer Sheva, Israel from 1985 through 2013. Energy Conversion and Management, v. 97, p. 307-314, 2015.
  • Ferreira Junior, R.A.; Souza, J.L.; Lyra, G.B.; Escobedo, J.F.; Santos, M.V.C. C. Energy conversion efficiency in sugarcane under two row spacings in northeast of Brazil. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 19, n. 8, p. 741-747, 2015.
  • Gairaa, K.; Khellaf, A.; Benkacilai, S.; Guermoui, M. Solar radiation measurements in Algeria: case of Ghardaïa station as member of the enerMENA meteorological network. Energy Procedia, v. 141, p. 50-54, 2017.
  • Grantham, A.P; Pudney, P.J.; Ward, L.A.; Belusko, M.; Boland, J.W. Generating synthetic five-minute solar irradiance values from hourly observations. Solar Energy, v. 147, p. 209-221, 2017.
  • Gueymard, C.A. Temporal variability in direct and global irradiance at various time scales as affected by aerosols. Solar Energy, v. 86, p. 3544-3553, 2012.
  • Iqbal, M. An Introduction to Solar Radiation New York: Academic Press, 390 p., 1983.
  • Kariuki, B.W.; Sato, T. Interannual and spatial variability of solar radiation energy potential in Kenya using Meteosat satellite. Renewable Energy, v. 116, Parte A, p. 88-96, 2018.
  • Kaskaoutis, D.G.; Kambezidis, H.D. The diffuse-to-global and diffuse-to-direct-beam spectral irradiance ratios as turbidity indexes in an urban environment. Journal of Atmospheric and Solar-Terrestrial Physics, v. 71, n. 2, p. 246-256, 2009.
  • Kotti, M.C.; Argiriou, A.A.; Kazantzidis, A. Estimation of direct normal irradiance from measured global and corrected diffuse horizontal irradiance. Energy, v. 70, n. 1, p. 382-392, 2014.
  • Kousky, V.E. Diurnal raifall variation on Northeast Brazil. Monthly Weather Review, v. 108, n. 4, p. 488-498, 1980.
  • Kousky, V.E. Frontal influences on Northeast Brazil. Monthly Weather Review, v. 107, n. 9, p. 1140-1153, 1979.
  • Lee, H-J.; Kim, S-Y.; Yun, C-Y. Comparison of Solar Radiation Models to Estimate Direct Normal Irradiance for Korea. Energies, v. 10, n. 5, p. 594, 2017.
  • Lima, M.C. Variabilidade da Precipitação no Litoral Leste da Região Nordeste do Brasil Dissertação de Mestrado em Meteorologia, INPE - 5283 - TDI/453, 1991.
  • Liou, K.N. An Introduction to Atmospheric Radiation. 2. Ed. San Diego: Academic Press, 2002, 583p.
  • Liu, B.Y.H.; Jordan, R.C. The interrelationship and characteristic distribution of direct, diffuse and total solar radiation. Solar Energy, v. 4, n. 3, p. 1-19, 1960.
  • Lyra, G.B.; Oliveira Júnior, J.F.; Zeri, M. Cluster analysis applied to the spatial and temporal variability of monthly rainfall in Alagoas state, Northeast of Brazil. Internat Journ Climat, v. 34, n. 13, p. 3546-3558, 2014.
  • Marques Filho, E.P.; Oliveira, A.P.; Vita, W.A. et al. Global, diffuse and direct solar radiation at the surface in the city of Rio de Janeiro: Observational characterization and empirical modeling. Renewable Energy, v. 91, p. 64-74, 2016.
  • Molion, L.C.B.; Bernardo, S.O. Dinâmica das chuvas sobre o Nordeste do Brasil. In: Congresso Brasileiro de Meteorologia, 11, Rio de Janeiro, p. 1334-1342, 2000.
  • Molion, L.C.B.; Bernardo, S.O. Uma revisão da dinâmica das chuvas no nordeste brasileiro. Revista Brasileira de Meteorologia, v. 17, n. 1, p. 1-10, 2002.
  • Moreno-Tejera, S.; Silva-Perez, M.A.; Lillo-Bravo, I.; Ramırez-Santigosa, L. L. Solar resource assessment in Seville, Spain. Statistical characterisation of solar radiation at different time resolutions. Solar Energy, v. 132, p. 430-441, 2016.
  • Muneer, T. Solar Irradiation and Daylight Models 2nd ed. Amsterdam: Elsevier, 2004.
  • National Research Council. 1992. Coastal Meteorology: A Review of the State of the Science. Washington: The National Academy Press, 1992. 112 p.
  • Oliveira, A.P.; Machado, J. A.; Escobedo, J. F.; Soares, J. Diurnal evolution of solar irradiation at the surface in the city of Sao Paulo: Seasonal variation and modeling, Theor. Appl. Climatol., v. 71, n. 3-4, p. 231-249, 2002.
  • Osinowo, A.A.; Okogbue, E.C.; Ogungbenro, S.B.; Fashanu, O. Analysis of Global Solar Irradiance over Climatic Zones in Nigeria for Solar Energy Applications. Journal of Solar Energy, Volume 2015, Article ID 819307, 9 p.
  • Pashiardis, S.; Kalogirou, S.A.; Pelengaris, A. Statistical analysis for the characterization of solar energy utilization and inter-comparison of solar radiation at two sites in Cyprus. Applied Energy, v. 190, p. 1138-1158, 2017.
  • Podstawczynska, A. UV and global solar radiation in Lodz, central Poland. International Journal Climatology, v. 30, n. 1, p. 1-10, 2010.
  • Polo, J.; Fernández-Peruchena, C.; Gastón, M. Analysis on the long-term relationship between DNI and CSP yield production for different technologies. Solar Energy, v. 155, p. 1121-1129, 2017.
  • Porfírio, A.C.S.; Souza, J.L.; Lyra, G.B.; Lemes, M.A.M. An assessment of the global UV solar radiation under various sky conditions in Maceió-Northeastern Brazil. Energy, v. 4 4, n. 1, p. 584-592, 2012.
  • Reboita, M.S.; Krusche, N.; Ambrizzi, T.; Porfírio, R. Entendendo o tempo e o clima na América do Sul. Terrae Didatica, v. 8, n. 1, p.34-50, 2012.
  • Ruiz-Arias, J.A.; Gueymard, C.A.; Santos-Alamillos, F.J.; Pozo-Vázquez, D. D. Worldwide impact of aerosol's time scale on the predicted long-term concentrating solar power potential. Scientific Reports, 6, 30546, 2016.
  • Santos, C.M.; Escobedo, J.F.; Teramoto, E.T.; Silva, S.H.M.G. Assessment of ANN and SVM models for estimating normal direct irradiation (Hb). Energy Conversion and Management, v. 126, p. 826-836, 2016.
  • Santos, C.M.; Souza, J.L.; Ferreira Junior, R.A.; Tiba, C.; Melo, R.O. et al. On modeling global solar irradiation using air temperature for Alagoas State, Northeastern Brazil. Energy, v. 71, p. 388-398, 2014.
  • Silva, G.R. Características de Vento na Região Nordeste: Análise, Modelagem, e Aplicações Para Projetos de Centrais Eólicas Dissertação de Mestrado em Engenharia Mecânica, Universidade Federal de Pernambuco, Recife, 141 p., 2003.
  • Souza, J.L.; Moura Filho, G.; Fonseca, R.F.; Teodoro, I.; Santos, E.A. et al. Análise da precipitação pluvial e temperatura do ar na região do Tabuleiro Costeiro de Maceió, AL, período 1972 - 2001. Revista Brasileira de Agrometeorologia, v. 11, n. 2, p.131-141, 2004.
  • Souza, J.L.; Lima, F.Z. Clima - Estudo de Impactos Ambiental, Levantamento Ambiental e Relatório, no Meio Ambiente de área de Proteção Ambiental de Piaçabuçu - AL In: Relatório Técnico, Maceió: UFAL, p. 180-238, 1995.
  • Souza, J.L.; Nicácio, R.M.; Moura, M.A.L. Global solar radiation measurements in Maceió, Brazil. Renewable Energy, v. 30, n. 8, p. 1203-1220, 2005.
  • Tabari, H.; Talaee, P.H. Sensitivity of evapotranspiration to climatic change in different climates. Glob. Planet. Chang, v. 115, p. 16-23, 2014.
  • Tahir, Z.R.; Asim, M. Surface measured solar radiation data and solar energy resource assessment of Pakistan: A review. Renewable and Sustainable Energy Reviews, v. 81, Parte 2, p. 2839-2861, 2018.
  • Tiba, C. Solar radiation in the Brazilian Northeast. Renewable Energy, v. 22, p. 565-578, 2001.
  • Vick, B.D.; Daryl, R.; Myers, D.R.; Boyson, W.E. Using direct normal irradiance models and utility electrical loading to assess benefit of a concentrating solar power plant. Solar Energy, v. 86, n. 12, p. 3519-3530, 2012.
  • Vilani, M.T.; Sanches, L.; Nogueiram, J.S.; Priante Filho, N. Sazonalidade da radiação, temperatura, e umidade em uma floresta de transição Amazônia Cerrado. Revista Brasileira de Meteorologiav. 21, n. 3b, p. 331-343, 2006.
  • WMO. Guide to Meteorological Instruments and Methods of Observation WMO-n°8, Seventh Edition, 157-197, Geneva, Switzerland, 2008.

Internet Resources

Nomenclature

    Publication Dates

    • Publication in this collection
      27 Nov 2020
    • Date of issue
      Oct-Dec 2020

    History

    • Received
      19 May 2019
    • Accepted
      19 Mar 2020
    Sociedade Brasileira de Meteorologia Rua. Do México - Centro - Rio de Janeiro - RJ - Brasil, +55(83)981340757 - São Paulo - SP - Brazil
    E-mail: sbmet@sbmet.org.br