Solar Radiation Components on a Horizontal Surface in a Tropical Coastal City of Salvador

: Renewable energy must be prioritized by humankind, mainly if there is an expected increase of 50% in energy consumption by 2030 and climate change scenarios are also conﬁrmed. Urban areas consume 70% of the available energy on the planet. Brazil, the largest country in South America, concentrates more than 85% of its population in urban areas, facing a challenge to increase the renewable power plants in its energy matrix. This work presents the solar radiation components behavior for the city of Salvador to contribute with initiatives for the use of solar energy resource. Firstly, a radiometric platform was implemented to obtain direct measurements of global ( E G ) and diffuse ( E DF ) components of incoming solar radiation at the surface. The knowledge of E DF is an important requirement to support photovoltaic system projects, and there is no information on direct measurements of this component in the State of Bahia. The diffuse radiation measurement device (DRMD) was designed and built for this purpose. The measurements of solar radiation components performed in this research were submitted to a speciﬁc quality control, statistically analyzed and used to evaluate the performance of different empirical models to represent the behavior of E DF in the tropical coastal city. The results demonstrate the potential to explore solar energy in the city of Salvador, with annual values of sunshine higher than 2200 h year − 1 and average daily intensities of E G and E DF equal to 18.7 MJm − 2 day − 1 and 7.35 MJm − 2 day − 1 , respectively. The analysis of the diurnal cycle shows that E G in summer reaches a maximum of 4.2 MJm − 2 h − 1 and in the rainy season it reaches a minimum of 3.7 MJm − 2 h − 1 , both at noon, and in summer the E hDF is 35% of the E hG and it is 46% in the rainy season.


Introduction
Energy is one of the most important resources available on Earth and humankind is clearly dependent on it to survive [1]. The urban areas are responsible for spending about 70% of the total energy generated in the world. These regions require large water reservoirs in order to produce electricity and stock fresh water for their own consumption [2]. The increase in energy demand and the constant fossil fuel exploration contribute to an unsustainable global energy production level [3,4].
Up to 2030, it is expected that the energy demand will increase more than 50% throughout the world [5]. At first, the COVID-19 pandemic changed the global energy consumption pattern, since billions of people on the planet changed their way of working and living, since most people were forced to slow down or even stop their work activities in order acquisition of E DF measures is not trivial, since they are reasonably expensive and with limited availability; therefore, many empirical correlations models have been developed for it as the most feasible alternative [32,33].
The evaluation of the solar radiation balance components has been performed over the last decades in Brazil. One of the first works to investigate the seasonal behavior of E G , E DF and direct (E DR ) solar radiation components for all Brazilian territory, was realized by Pereira et al. (1996) using a radiative transfer model derived from satellite measurements, finding out a maximum of E G daily value bigger than 20 MJm −2 day −1 over the NEB during the winter [31]. Martins et al. (2012) used an updated version of the model proposed by Pereira et al. (1996) and demonstrated that the NEB shows average annual values of E G close to 23.4 MJm −2 day −1 in drought areas and 18.7 MJm −2 day −1 near to oriental NEB [34]. The results obtained by Pereira et al. (2017) in the Brazilian Solar Atlas (BS-Atlas) confirm that NEB has a high solar radiation power potential available, and the city has an incredible potential to generate electricity by solar energy resource according to Brazilian Solar Atlas (BS-Atlas), which estimates a value of 20.0 MJm −2 day −1 for the annual average of daily global radiation [20]. These estimates were generated by the BRASIL-SR model, and they were validated against the collected observational measurements provided by surface stations from National Meteorology Institute (INMET).
The hourly and seasonal evolution of the components of the solar radiation balance based on radiometric measurements collected in São Paulo City were investigated by Oliveira et al. (2002) [35]. The results acquired for E G in summer (rainy season) were 18.41 MJm −2 day −1 , 55% over the winter (dry season) results. At the same time, the E DF values obtained reach a maximum of 9.52 MJm −2 day −1 in summer and a minimum of 4.1 MJm −2 day −1 in winter, matching with the estimates obtained by Pereira et al. (1996) [26]. Oliveira et al. (2002) suggested that cloudy weather and atmospheric pollution produce a strong influence on the scattering of solar radiation and this hypothesis was confirmed by Codato et al. (2008) [36].
In a similar work, conducted in Rio de Janeiro City, Marques Filho et al. (2016) got a daily maximum value of E G equal to 23.89 MJm −2 day −1 in summer and 10.78 MJm −2 day −1 in winter, presenting superior values in comparison with the same analysis performed in Sao Paulo City, although the cities are part of the same latitude range. The E DF component reaches a maximum of 7.34 MJm −2 day −1 in summer, which is twice more than the values obtained in winter. In both works, the cloudiness and marine aerosols influence the components and balance of solar radiation on the surface [37].
The seasonal behavior of E G for Alagoas (NEB State) was presented by De Souza et al. (2016) [38]. Based on radiometric measurements collected in three different experimental sites, a new coefficient for the Angström-Prescott empirical model was proposed [39] that uses information of S to predict a seasonal evolution of E G and can be applied to monthly periods with no temporal dependence for successive values of radiation [40]. These models reach better estimates, but the adjusting coefficients get a high variability.
The main objective of this work is to describe the seasonal variation of daily and hourly values of incoming solar radiation components, mainly E G , E DF and E DIR , at the surface in Salvador City, using in situ measurements and empirical models, because Salvador has a large potential for solar radiation, which has not been fully documented and confirmed in the literature using in situ observations. Due to the lack of surface measurements of E DF , for the first time, a specific equipment, known as diffuse radiation measurement device (DRMD), was designed and built for this research, to perform direct measurements of E DF in the City of Salvador, and continues in operation at the Federal University of Bahia (UFBA), contributing to provide measurements of the diffuse component to the solar energy research field.
The present research covers the construction of a diffuse radiation measurement device, experimental data collection on the surface, quality controls of radiation measurements, statistical analysis and predictions of solar radiation components using predict models. The results found in this work give information about the behavior of radiation components, sunshine duration during the year and the seasons and climatic characterization of the city that can be used in a infinitude of analyses related to solar radiation and photovoltaic generation; in addition, the empirical models developed here performed very well, and they can be used to estimate the diffuse component of solar radiation in other regions throughout Brazil with similar climate and characteristics to Salvador City.

Site, Instrumental and Database
In this work, solar radiation components and standard meteorological measurements collected by INMET and Micrometeorology and Modeling Laboratory (LabMiM) located in Salvador City were used. The LabMiM is located in a suburban area, with buildings up to four floors, surrounded by medium-sized vegetation, installed on top of a 26 m high building of the Center for Energy and Environment (CIEnAm), within the Ondina Campus of Federal University of Bahia (UFBA). The INMET meteorological station is located in the parking area, which is predominantly vegetated, and follows the standards recommended by the World Meteorological Organization (WMO). The distance between the INMET and LabMiM stations is only 0.6 km, with a vertical elevation of 46 and 48 m above mean sea level (AMSL), respectively. For instance, the altitude correction was performed in thermodynamic variables at LabMiM, based on extrapolation of dry adiabatic lapse rate down to the surface.
The climate characteristics for Salvador City were evaluated based on daily and hourly data, collected from 2010 to 2019. The climate normal (1981-2010) was provided by INMET and confronted with the seasonal estimates obtained from the measurements collected during the last decade. The location of each station, the variables collected and the observational periods are shown in Figure 1 and Table 1. The LabMiM instrumental set is systematically calibrated (every 6 months) and the maintenance is done weekly, with the objective to grant the proper functionality as well as to preserve the quality of the collected measurements.
In order to obtain directly measured data from the E DF component in Salvador City, we developed and manufactured a device named the diffuse radiation measuring device (DRMD) since 2018, based on Melo-Escobedo-Oliveira shadow-ring measuring method (MEO) presented by Dal Pai et al. (2016) [41]. Some design changes were implemented as follows: (i) latitude (φ) adjustable mechanism, with a range of 0 to 45 • ; (ii) completely assembled in structural aluminum, discarding the use of welds, with reduced weight and increased resistance to weathering corrosion, as well as the physical properties to provide mechanical strength; (iii) implementation of stainless steel threaded bar, bearings and shafts to support and perform horizontal movement of the pyranometer base; and (iv) use of a quick mechanism to turn the lever and lock the shadow ring. Figure 2 presents a 3D model and the manufactured device in operation.  In order to ensure the quality of the collected measurements, a visual inspection of each time series was performed to remove physical inconsistencies generated by the data acquisition system. Statistical criteria based on absolute values and confidence intervals were applied to remove spurious data [42].
The first step consisted in keeping the hourly radiation data. For daily data, months with less than 80% of measurements were discarded. After that, we kept 71.2% of the raw data. This step it is important to keep the analysis consistent, since sometimes the system can have a lack of communication between the sensor and the data logger and must be guaranteed that, for each measurement of E G related to E DF at the same time, because the analysis depends on the relationship between both solar radiation components.
Additionally, the quality control procedure that was used was based on some criteria proposed by Younes et al. (2005) and Journee et al. (2011), in order to exclude spurious data caused by some errors such as the following [43,44]: operation-related problems and errors; complete or partial shade-ring misalignment; dust, dew, water droplets, bird droppings, etc.; incorrect sensor levelling; shading caused by building structures; electric fields in the vicinity of cables; mechanical loading of cables; orientation and/or improper screening of the vertical sensors from ground-reflected radiation; and station shut-down.
After the first step the quality criteria of Equation (1) were applied. The criteria consist of keep the radiation components measurements where the elevation solar angle (α) > 2 • , in this case, the total of the measures kept from the initial raw data was 69.5%.
A comparison was realized between the hourly global radiation (E h G ) and theorical hourly extraterrestrial radiation (E h T ), because the first must have lower values than (E h T ). Equation (2) was used to perform these criteria, where K h T is the hourly clearness index. After this process, the total of the measures kept from the initial raw data was 68.3%.
This step is similar to the previous steps. The main difference is the comparison of the hourly diffuse radiation E h DF with (E h T ). In this case the quality control for E h DF is given by Equation (3). After the control the total of the measures kept from the initial raw data was 67.8%.
This last step of the quality control consists of comparing the hourly diffuse radiation where E h DIR is the hourly direct radiation. After the step, the total of the measures kept from the initial raw data was 67.3%.
The E T component was estimated analytically from the equations proposed by Iqbal (1983), considering the solar constant equal to 1.366 Wm −2 [45,46]. After quality control was applied, 67.3% of the radiometric measurements were considered valid.
To verify the consistency of the tests, the measurements of E d G collected by INMET were compared to those collected in the LabMiM (Figure 3), and the last one showed a greater control and data traceability than the first one due to the constant maintenance and verification of the measuring instruments. The behavior of E d G in both experimental sites is similar, with a linear fitting near to the identity function (0.98) and a high coefficient of determination (0.99). The smaller standard deviations of LabMiM (horizontal bars) show the representativeness of this dataset in comparison to INMET (vertical bars). The average block indicates that INMET underestimates E G only for values greater than 25 MJm −2 day −1 , which represent 13% of the measurements in the period. The scattering presented for INMET measurements when E d G assumes values between 6 MJm −2 day −1 and 10 MJm −2 day −1 , occurs for α < 5 • , and it should be occasioned by obstacles around the sensor.  (Table 1) and, therefore, this sample time is not enough to statistically characterize the climatic condition observed. In this case, the alternative found was to use empirical models [32,33] to represent the behavior of E DF in the related city. In the next section, some empirical models of the E h DF proposed by the literature will be discussed.

Diffuse Empirical Models
The accurate determination of the E DF values is an essential precondition in the construction of photovoltaic system design and to evaluate the performance of photovoltaic panels [47,48]. However, E DF measurements are difficult to obtain, and they are not part of the worldwide meteorological network. In this case, to estimate the hourly, daily and monthly E DF values in different locations on Earth, it has been proposed to use several empirical models [49]. Some existing empirical models represent the relationship between clearness index (K T ) and diffuse fraction (K D ) through the linear, polynomial, logarithmic or logistic sigmoidal functions that use one or more input parameters [32].
The pioneer Liu-Jordan (1960) work established a function based on values of K T to represent the behavior of K D , and to determine the intensity for E DF on a horizontal surface. The coefficient K T , defined as the ratio E G /E T , provides a measure of the effects of the local atmosphere on the incident solar radiation on the surface [50,51]. When K T → 0 , the sky condition is defined as a cloudy sky with a high portion of E DF . On the other hand, in a clear sky condition, K T → 1 and there is a large fraction of E T the surface, with a predominance of E DR .
The E DF component is an important fraction of E G and a key variable in solar energy evaluation [52]. The coefficient K D , defined as the ratio of E DF /E G , has great importance in the assessment of the effective solar radiation collected by a surface [53,54]. When K D → 0 , means that the fraction of E DR is predominant over the E DF . On the other hand, when K D → 1 , a large fraction E DR is being prevented from reaching the surface, with the predominance of E DF .
The polynomials derived by linear regression fit were able to model the daily and monthly values of diffuse radiation satisfactorily [35]. The main advantage is how easy it is to obtain the input parameters [55]. Otherwise, the models are not applicable for all intervals of K T , and they do not represent the variability observed in the Liu-Jordan diagram.
The sigmoid logistic function is statistically more significant in comparison with other correlation models to represent the K D as a function of K T , and it provides a physically reliable estimate for extreme values of the clearness index, though using fewer parameters than other tested models [56].
The Marques Filho et al. (2016) model (Marques Filho) proposed the following sigmoidal model to represent the behavior of the K T and K D radiation fraction observed in do Rio de Janeiro City [37]: The main disadvantage of Marques Filho model is that it does not reproduce the observed variability of K D , which used K T as the only input parameter [37].
Since the K D has a direct relationship with the sky condition, to improve the accuracy of empirical models, it is necessary to evaluate other parameters that include effects of cloud (cloudiness and cloud type), traditional meteorological variables (air temperature, relative humidity and atmospheric pressure observed at the surface) and air pollution (concentration of particulate matter observed at the surface) [57].
The Ridley et al. (2010) model (Ridley) is a multiple predictor empirical model developed based on a sigmoidal function that represents better (low) variability of K T in relation with K D , mainly for Southern Hemisphere locations, which has been used as a universal predictor model for E DF , according to Equation (6) [58].
where K h T is the hourly clearness index, K d T is the daily clearness index, α is the apparent solar time (AST) and the persistence of E G level (ψ) is used as predictors.
The Lemos et al. (2017) model (Lemos) was based on measurements of nine radiometric stations distributed throughout Brazilian territory, using the adjusted Ridley model in order to provide better estimates of E DF , according to Equation (7) [30].
In this work, the sigmoidal models proposed by Ridley [30] are evaluated to determine which model better represents the K D behavior in Salvador City.

Statistical Indicators
By evaluating the performance of the models when reproducing the values of E DF , we used a set of statistical indicators presented by Gueymard et al. (2014) and Funabashi (2016) [59,60].
The coefficient of determination (R 2 ) is a statistical measure of how close the data are to the fitted regression line. In this measure, 0 indicates that the model does not explain the variability of the response data at all, while a value around and at 1 means the regression actually fits the observation [61,62].
where p i and o i are the predicted and observed values, respectively, while o m and p m represent the average values of o i and p i , respectively. Additionally, n denotes the entire test data amount. Mean bias error/difference (MBE) captures the average bias in the prediction. The test of MBE provides information on the long-term performance of models studied. A positive MBE value gives the average amount of overestimation in the calculated values and vice versa. In general, a small MBE is desirable. It should be noted, however, that over-estimation of an individual observation will cancel out the underestimation in a separate observation [63].
The root mean square error (RMSE) denotes the root of the ratio of the square of the two norms of the error vector to the number of samples. RMSE provides information on the short-term performance of the models as it allows for a term-by-term comparison of the actual deviation between the calculated value and the measured value [64,65].
where N represents the number of observations. The t − statistic (TS) is considered a robust and great indicator used to rank the modeled results. When the (TS) value approaches zero, it means that the probability of significant difference between the analyzed values is lower [66,67]. To obtain the indicator there is a combination of MBE and RMSE. The t − statistic is defined by the following: In order to compare and evaluate the models, particularly to determine which models are more accurate, the index of agreement (d) was developed to measure errors from prediction models and to promote the standardization of measurements. The parameter can be easily interpreted by making cross-comparisons for a variety of models. When d → 1 , this indicates a great alignment between prediction and observational values; otherwise, when d → 0 , there is a complete misalignment [68,69]. The index of agreement is calculated by the expression:

Climatology
According to Köppen-Geiger's climate classification and climate normal estimates of the observations carried out for 30 years (1981-2010), the climate of Salvador City can be classified as Tropical Rainforest (Af) [70]. The annual amplitude of monthly average air temperature T 1981−2010 is small, with a minimum of 23.7 • C in August and a maximum of 27.1 • C in February and March (Figure 4a). The monthly average air relative humidity RH 1981−2010 is high throughout the year due to the vicinity of the South Atlantic (SA) ocean, with values higher than 78.5%. In Figure 4a, it can be observed that the last years have been drier with summers that are slightly warmer, in comparison to the regular average climate. In both periods evaluated, the estimates of T 2010−2019 and T 2017−2019 fit with the regular climate.
The estimates of RH present a high variability, remaining under the values registered on regular climate. In recent years, the RH 2017−2019 was wetter, but it kept the estimates within the deviation bars.
The seasonal evolution of monthly accumulated rainfall in the period from 2010 to 2019 presents a similar behavior of regular climate, with values 9.5% less in the average annual accumulated rainfall (Figure 4b). For the period between 2017 and 2019, the variability of accumulated rainfall was higher than the others periods during the year, principally in March, when the accumulated was 60% higher than observed in the period between 2010 and 2017. The previous result of rainfall is in accordance with what was discussed to RH in Figure 4a.
The seasonal evolution of monthly accumulated rainfall in the period from 2010 to 2019 presents a similar behavior of regular climate, with values 9.5% less in the average annual accumulated (Figure 4b). For the period between 2017 and 2019, the variability of accumulated rainfall is higher than the other periods during the year, principally in March, when the accumulated was 60% higher than observed in the period between 2010 and 2017. The previous result of rainfall is in accordance with what was discussed to RH in Figure 4a.
The monthly average values of sunshine hours (Figure 4b) are minimum, between the months of April and July, matching with the rainy season and reaching a maximum level in summer (the driest season). Even by having a high concentration of water vapor in the atmosphere, the region has a great availability of sunlight between the months of September and March with values greater than 200 h/month −1 .
The observation performed for both periods (2010-2019 and 2017-2019) is in accordance with the regular average climate. The differences found, however, may be related to the recorded period and the use of the land INMET (vegetated area) and Lab-MiM (suburban area). The results found are also in accordance with the climate change predictions for the NEB, which show an increase in air temperature and a decrease in accumulated rainfall [18,19].

Evaluation of Empirical Models
As previously informed in this study, the sigmoidal empirical models proposed by  [30] were applied to obtain E DF from K D for Salvador City. Ridley and Lemos are physically more consistent as they consider a total of five multiple predictors to represent K D behavior, which increases the accuracy of the estimated results [74,75]. On the other hand, Marques Filho developed a study for a coastal tropical region of Brazil with climate conditions similar to Salvador City. In order to verify what model better reproduces the behavior of K h D , we used the observational data of E h G and E h DF collected by LabMiM ( Table 1). The behavior of K h D estimates for Salvador reproduces the Liu-Jordan diagram, and all empirical models applied there are within the observational data range ( Figure 5). Considering the multi-parameter of sigmoidal models, Ridley and Lemos best represent the observational variability, something that does not occur for Marques Filho due to the simplicity of this last model that only used K h T as an input parameter. At the same time, when the average error bars of the experimental data are compared, Lemos underestimates the values of K h D when K h T > 0.4. In order to evaluate the performance of the empirical models applied, the statistical indicators described in Section 2.2 are shown in Table 2. All evaluated models present R 2 close to 0.90, with the best result (performance) acquired by Ridley. The bias evaluation evidence that Ridley and Marques Filho overestimates the measures, with MBE equal to 4.05% and 8.22%, respectively, while Lemos measures in 6.13%, as discussed above ( Figure 5). Although Lemos has been developed for Brazil, a large part of the experimental data used was collected in regions with distinct climate conditions and far from the coast, differently from Salvador City, which is located in a tropical latitude with a strong influence of marine sprays. These conditions are closer to those used by Ridley, which was developed for coastal cities of the Southern Hemisphere.

Solar Radiation Components
The seasonal evolution of the monthly average, collecting daily values of solar radiation components at the surface in Salvador City, is presented in Figure 6.  This result demonstrates a top potential season to explore solar energy in Salvador City since, in Brazil, the greater demand for energy occurs in summer [76,77]. In the rainy season, however, between April and June, the nebulosity increases (Figure 4b), and the values of E d G consequently decrease. In order to compare the seasonal solar radiation intensity in some locations of Brazil, the values of E d G to Salvador, Rio de Janeiro, São Paulo and Maceio City are presented in Table 3. The minimum values presented in Table 3 are related to the rainy season (April-June) for Salvador and Maceió [30], while for São Paulo [35] and Rio de Janeiro [37] they are related to winter. Since the maximum values are presented in Table 3 and occur during summertime, in this period the electricity consumption trends change, and it has been registered in recent years that the demand peak in urban areas happens in the hours of maximum air temperature, between 13 h and 15 h, because of the use of air-conditioning [78]. In the rainy season, there is a considerable drop of E d G 14.83 ± 2.1 MJM −2 day −1 , with estimates 65% lesser than the observed in summer and near those of wintertime 14.98 ± 1.8 MJm −2 day −1 .
In order to verify the E d G values' consistency, the model proposed by Angström-Prescott for solar radiation prediction was applied. In general, this model is used to obtain E d G a simple linear equation that uses S as the parameter, and coefficients a and b are estimated according to the ordinary regression approach [79]. The Angström-Prescott equation estimated for the city of Salvador is given by: where a and b are equal to 0.262 and 0.538, respectively, and where S and S 0 are the monthly average and maximum month average daily values of sunshine hours, respectively. The Angström-Prescott model adjusted for Salvador City reproduces the annual variation of the monthly average of E d G with RMSE equal to 3.24% and its correlation coefficient R 2 is equal to 99.6%. These results show that the estimates for E d G represent the climatic conditions observed before in the related city.
Additionally, to verify the representativeness of the solar radiation components to Salvador City, the estimates of the monthly average of E d G and E d DF were compared to those presented in the BS-Atlas dataset [20].
The newest edition of BS-Atlas, released in 2017, was used for this purpose. The solar radiation components published by BS-Atlas were estimated by a BRASIL-SR model that used images collected by Satellite GOES, and it was validated by ground stations' measurements from the SONDA Network project during the period defined between 1999 and 2015 [20]. In this work, the solar radiation components database is related to the coordinates (−13.00, −38.45) with a grid size of (10 × 10 km).
In comparison to BS-Atlas, E d G is overestimated, mainly between the spring and the beginning of summer, when RMSE is equal to 6.7% and R 2 is equal to 0.97. Nevertheless, E d DF is underestimated by BS-Atlas (Figure 7). The correlation coefficient for E d DF , between BS-Atlas and Ridley, is 0.92. The differences found of E d DF between these models may be related to the performance of radiative transfer model BRASIL-SR and conditions of maintenance of radiometric ground stations [80]. The diurnal evolution of monthly average E h T and E h G during summer and rainy season months are presented in Figure 6. Both variables show a well-defined diurnal cycle with a maximum peak at noon. In the summertime, the integrated value along the day of E h T and E h G is 27% and 46.0% higher in comparison to the values observed in the rainy season, respectively. During summer, the E h DF is 35% of the E h G , while this value in the rainy season is E h DF 46% higher than E h G . The vertical bars correspond to standard deviation. The monthly variation of diurnal evolution E h G shows a great availability of solar energy incidents on the surface throughout the year, between 10 h and 15 h, except for the rainy season ( Figure 6). In summer, we may observe E h G maximums up to 4.2 MJm −2 h −1 (February) and in rainy season a minimum of 3.7 MJm −2 h −1 , both at noon. The highest values of E h G in summer are a combination of astronomical factors and atmospheric conditions, since these are the driest months of the year. The longest and shortest days of sunshine hours occur in December (10.2 h) and July (9.1 h), respectively.

Salvador City Solar Energy Poential
This research provides substantial results regarding the behavior and information of the solar radiation components at Salvador City. In this section we present some comparisons between Salvador and other cities around the world.
The general pattern of E G for each place depends on geographic coordinates and is related with latitude, sky conditions, clouds, atmospheric aerosol and altitude, for example. Around the world, the solar radiation in general is higher in arid tropical and subtropical zones. Some countries such as South Africa, Australia, Mexico, Brazil and the United States are located in these zones [81]. Although these areas have greater potential for photovoltaic generation, some countries in Europe and Asia are included in the largest producers of energy through photovoltaic power plants [9]. Table 4 presents the solar PV global capacity from 10 different countries around the world and the theoretical potential for PV power [9,81].  Table 4 was adapted using information available from the Global Solar Atlas [81] and the REN21 Global Status Report [9]. Is important to remember that solar power generation depends on the type of cells and solar panel, time to solar exposition (sunshine hours), the seasons and other parameters.
In Table 4 is possible to see how the potential of solar generation in Brazil is poorly explored. In comparison with other countries, Brazil has the second largest theoretical potential and the third biggest evaluated area, using only 0.06% of PV equivalent area. If compare Brazil with Germany, which is one of the largest PV power generators of world, it is possible to see that the tropical country gets a theoretical potential that is 77.2% higher and evaluated area that is 24 times larger; however, the European country has seven times more PV capacity installed.
The collected measurements of solar radiation components performed by LabMiM sensors show that Salvador City has a measured potential for E d G equal to 18.65 MJm −2 day −1 , which is close to the theoretical average potential predicted to Brazil by the Global Solar atlas. The potential for E d DF for Salvador is equal to 7.35 MJm −2 day −1 , 28.7% less than the theoretical average potential predicted for the Netherlands by the Global Solar atlas, and still remembering that Salvador City has more than 2200 h of sunshine a year. In a hypothetical situation, photovoltaic panels could generate a significant amount of electricity in the city of Salvador, just considering the incidence of E DF , and it would be higher than in many countries where the PV panels are pointed towards the sun. Table 5 presents the values of E d G to some cities worldwide [82]. The cities chosen are spread around the world from five different continents, with different geographical coordinates and a ranking.

Conclusions
Nowadays, solar energy resource is the most promising resource available around the world. Brazil, especially where NEB is located, is a tropical zone with high effectiveness of solar radiation [20,31]. The radiation components balance was evaluated for Salvador City, which is the biggest and the most populous city in NEB, with a population of about three million people [12], who are using and consuming energy on a daily basis.
Salvador is climatologically classified as a Tropical Rainforest (Af) zone, where air temperature varies between 23.7 • C in winter (rainy) and 27.1 • C in summer, with relative humidity higher than 78.5%. The rainfall registers a minimum of 58.1 and a maximum of 1029 mm month −1 in July (winter/rainy) and December (summer), respectively. The sunshine minimum terms are between the months of April and July and the maximum are between the months of September and March (over 200 h/month −1 ). These observations are according to the normal climate behavior.
Observational data powered by INMET and experimental data collected from LabMiM were used in order to evaluate the main radiation solar components. Since the E DF data were not available, it was necessary to design and build a diffuse radiation measuring device (DRMD), based on MEO (Melo-Escobedo-Oliveira shadow-ring measuring method). In addition, to perform the analysis, some empirical models were used and compared with the observational measurements and BS-Atlas.
In this work, sigmoidal models, such as the Ridley; Marques Filho, and Lemos models, were evaluated to determine which model better represents the K D behavior of Salvador City. After the statistical analysis, all models obtained good results (R 2 close to 0.90). The Ridley and Marques Filho models overestimated the measures, with MBE equal to 4.05% and 8.22%, respectively. On the other hand, the Lemos MOD underestimated the measures by 6.13%. Ridley presents all statistical indicators better than the others (Table 2), with RMSE equal to 22.15%. Although the Lemos model was used in Brazilian territory, the Ridley model got better results for Salvador City, since it was built for coastal cities in the Southern Hemisphere.
The observational data values in comparison with BS-Atlas show that E d G is overestimated, mainly between the spring and the beginning of summer, with the RMSE equal to 6.7% and R 2 equal to 0.97. However, the E d DF is underestimated by BS-Atlas ( Figure 6). The correlation coefficient for E d DF , between BS-Atlas and Ridley, is 0.92. The differences found of E d DF between these models are related to the performance of radiative transfer model BRASIL-SR and maintenance conditions of radiometric ground stations.
Salvador City has a well-defined diurnal cycle, reaching a maximum value of radiation at noon. During the summer, the E h G is predominant and E h DF reaches 35%, while in the rainy season, the E h DF increases and reaches 46% of E h G . The longest and shortest sunshine hours occur in December and July, with 10.2 and 9.1 h, respectively. In the summertime, a maximum of E h G goes up to 4.2 MJm −2 h −1 and it occurs in February, while the minimum of 3.7 MJm −2 h −1 occurs in the rainy season.
The research described the seasonal variation of daily and hourly values of incoming solar radiation components, mainly E G , E DF and E DIR , at the surface in Salvador City, using in situ measurements and empirical models, and they results confirm the that the city has a large potential for solar radiation if compared with countries and cities around the world.
The measured potential for E d G was equal to 18.65 MJm −2 day −1 , which is next to that predicted for Brazil by the Global Solar atlas. The quality and potential of solar radiation found for Salvador is great, until the values of E d DF of 7.35 MJm −2 day −1 can reach values of E d G for other countries available on the Global Solar map. Another advantage of the city is its annual values of sunshine higher than 2200 h year −1 .
An objective comparison between Salvador and nine of the biggest cities in world shown that the Brazilian tropical city reached the third place on the list, regarding the E d G potential.
The design and construction of the diffuse radiation measurement device (DRMD) brought originality and novelty for the research, and it was fundamental to perform direct measurements of E DF . The device continues in operation at the Federal University of Bahia (UFBA), contributing to provide measurements of the diffuse component to the solar energy research field.
The present research covers the construction of a diffuse radiation measurement device, experimental data collection on the surface, quality controls of radiation measurements, statistical analysis and predictions of solar radiation components using predict models. The results found in this work give information about the behavior of radiation components, sunshine duration during the year and the seasons and climatic characterization of the city that can be used in a infinitude of analyses related to solar radiation and photovoltaic generation, in addition the empirical models developed here, which performed very well and t can be used to estimate the diffuse component of solar radiation in other regions throughout Brazil with similar climate and characteristics as Salvador City.