UV and global irradiance measurements and analysis during the Marsaxlokk ( Malta ) campaign

A solar radiation measurement campaign was performed in the south-eastern village of Marsaxlokk (3550 N; 1433 E; 10 m a.s.l), Malta, between 15 May and 15 October 2012. Erythemal solar radiation data (from a UVB-1 pyranometer), and total horizontal solar radiation (global and diffuse components) from two CM21 pyranometer were recorded. A comparison of atmospheric compounds from ground measurements and satellites shows that TOC (total ozone column) data from the Ozone Monitoring Instrument OMI, TOMS and DOAS algorithms correlate well with ground-based recorded data. The water vapour column and the aerosol optical depth at 550 nm show a significant correlation at the confidence level of 99 %. Parametric models for evaluating the solar UV erythemal (UVER), global (G) and diffuse (D) horizontal irradiances are calibrated, from which aerosol effects on solar irradiance are evaluated using the Aerosol Modification Factor (AMF). The AMFUVER values are lower than AMFG, indicating a greater aerosol effect on UVER than on global solar irradiance. In this campaign, several dust event trajectories are identified by means of the HYbrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model and by synoptic conditions for characterizing desert dust events. Hence, changes in the UV index due to atmospheric aerosols are described.


Introduction and objectives
The UV erythemal (UVER) radiation affects human health (short-term effects include erythema or sunburn and longterm effects include photo-ageing/skin cancer), damages aquatic life, and affects plants as well as the conservation and durability of materials, in addition to impacting global energy balance and climate change (UNEP, 2010).Yet, not all its effects are harmful, with the synthesis of vitamin D being one of the beneficial effects of UV (Webb, 2006;Fioletov et al., 2009).Different action spectra are used to quantify these effects.McKinlay and Diffey (1987) established the erythemal action spectrum, which represents the spectral response of human skin to UV radiation being able to trigger erythema (or sunburn).Solar radiation weighted using this spectrum is called erythemal radiation (UVER; 280-400 nm) (Webb et al., 2011;CIE, 1998).
Various authors have proposed new models for estimating solar irradiance (e.g.Mateos et al. (2010) who assessed ultraviolet and global cloud modification factors in central Spain).The attenuation of UVER by different atmospheric components has been analysed by De Miguel et al. (2011); the results confirm that small aerosols scatter more solar UV wavelengths than larger ones.Román et al. (2013) described the effects of a desert dust episode at Granada (Spain) on global, direct and diffuse spectral UV irradiance, and reported that the attenuation of direct UV was about 50 %, while the diffuse irradiance increased by up to 40 %.Authors such as Pace et al. (2006) have studied the influence and identification of different aerosol types at Lampedusa Island (Central Mediterranean); their observations were combined with air mass trajectories in order to identify different types of particles and to determine their mean optical properties.
Also it is known that aerosols are atmospheric particles that can affect incoming solar radiation directly by absorption and scattering, and indirectly by acting as cloud condensation nuclei (modifying cloud microphysical properties).
Both effects contribute to cooling the Earth's surface and simultaneously warming the lower atmosphere (IPCC, 2007).
The aerosol radiative forcing in the Central Mediterranean has also been studied (Gómez-Amo et al., 2011); the results show the effects of different aerosols on the radiative budget during a dust event within the boundary layer.The sensitivity of shortwave radiative fluxes to changes in the vertical distribution of aerosols and a case of Sahara dust layer above urban aerosols have also been studied (Gómez-Amo et al., 2010;Kaskaoutis et al., 2010a, b).
The present work seeks to explore the influence of aerosols on solar irradiance in the area of Malta.In the study, parametric models for evaluating solar UVER, global (G) and diffuse (D) horizontal irradiances under cloudless conditions were calibrated, from which the aerosol Modification factor (AMF) was evaluated.The models are a function of the solar zenith angle (SZA), the total ozone column (TOC), and the aerosol optical depth at 550 nm (AOD 550 ).The models were validated through a graphical comparison of measured and estimated data and by different statistical indices (Bilbao et al., 2014).Atmospheric ground reference data were compared to those retrieved from satellites; the agreement was assessed.In this campaign, several dust event trajectories were identified by means of the HYbrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model and by the synoptic conditions (Bilbao et al., 2014).Changes in the UV index due to atmospheric aerosols were thus described.
In the following sections, the site description, instrument and data collection are detailed, the methodology is explained, and the observations and results are analyzed in relation with the influence of aerosols on UVER, G and D irradiances.

Measurement site
The Universities of Malta and Valladolid (Spain) conducted a solar radiation measurement campaign, which took place at the Institute for Sustainable Energy in the south-eastern village of Marsaxlokk (35 • 50 N; 14 • 33 E; 10 m a.s.l.), Malta.Measurements were recorded between May and October 2012.The campaign involved two European Institution Groups: the Institute for Sustainable Energy at the University of Malta, and the Atmosphere and Energy Laboratory at the University of Valladolid, Spain.Solar radiation instruments used for this study were located on the rooftop of the Institute for Sustainable Energy, which had an obstruction-free horizon.

Ground-based instruments
Global and diffuse horizontal irradiance were recorded using two CM21 (Kipp & Zonen) pyranometers, one of which was equipped with a shadow-band.Diffuse solar irradiance data were corrected following the method proposed by Batlles et al. (1995) and Perez et al. (1990) that take into account geometric and atmospheric (clearness index) corrections.The sampling rate of the global and diffuse irradiance pyranometers was 10 s although measurements were averaged and recorded at 1 and 10 min.As a result, 10 min solar irradiance measurements were obtained.The CM21 (Kipp & Zonen) instruments have a flat spectral response from 305 to 2800 nm and the cosine effect is below 3 % for a solar elevation above 10 • .The CM21 sensors were regularly calibrated by comparison with a reference sensor at the Kip & Zonen manufacturer; the differences obtained were below 2 % (Bilbao et al., 2014).
The UVER measurements were recorded using a Yankee Environmental Systems (YES) UVB-1 radiometer, which has a spectral range between 280 and 400 nm and a spectral sensitivity which resembles the erythemal action spectrum.The UVB-1 sensor is designed to operate continuously and autonomously in the field (Esteve et al., 2009).The cosine response is greater than ±5 % for solar zenith angles (SZA) below 60 • .This instrument is calibrated by the standard rules in the National Institute for Aerospace Technology (INTA) in Spain.The calibration consists of measuring the spectral response of the sensor indoors and comparing with a Brewer MKIII spectroradiometer outdoors.After this process, a double input matrix with the calibration factors depending on SZA and the total ozone column (TOC) was obtained using a radiative transfer model (Hülsen and Gröbner, 2007).The error given by the calibration matrix is below 9 % for a SZA of less than 70 • .The output-voltage signal of the sensor is converted into units of erythemal irradiance (W m −2 ) by multiplying the signal voltage by the calibration factor obtained from the matrix, taking into account the average SZA during the 10 min measure and the daily TOC.The experimental uncertainty of the sensor according to the results is in the 4.6-7 % range (De Miguel et al., 2012).
A Solar Light Microtops-II manual Sun photometer was used throughout the campaign to measure TOC and AOD at 1020 nm (AOD 1020 ), and water vapour column.The Microtops-II Sun photometer is equipped with five optical collimators, with a field of view of 2.5 • , to perform direct radiation measurements at the following nominal wavelengths: 305, 312, 320, 936, and 1020 nm.The instrument is calibrated every two years at the Mauna Loa Observatory (Hawaii) by Solar Light Company (Glenside, PA), further details are available in De Miguel et al. (2011).Furthermore, an inter-comparison of ozone data between the Microtops-II and the Brewer spectrometer was performed at Lampedusa island (in the Central Mediterranean) over five days.The differences observed between the TOC estimations of the two instruments were below 2 % (Mateos et al., 2014b).More details concerning solar irradiance sensors are given in De Miguel and Bilbao (2005).
Solar sensors were connected to a Campbell CR10X Data Logger, which was programmed to take measurements each 10 s from which 10 min average values were computed and stored.There are 6 readings for every 10 min rintervals in an hour which were added together providing a single hourly reading.Hourly and daily irradiances were then evaluated from the 10 min average values (Bilbao and de Miguel, 2013).Data were transmitted continuously from data-logger to the Atmosphere and Energy Laboratory via web.Measurements were collected in Greenwich Mean Time (GMT).

Satellite-borne instruments
TOC, aerosol optical depth measured at 550 nm wavelength (AOD 550 ), and precipitable water column (w), were retrieved from a satellite sensor (MODIS OMI).During the campaign, whenever the Microtops TOC measurements were not available because of cloudy conditions, data from the Ozone Monitoring Instrument (OMI; onboard Aura satellite) were used.If TOC from OMI was not available, data from the Global Ozone Monitoring Experiment-2 instrument (GOME-2; onboard MetOp-A satellite) were used.These satellite-based data were obtained from the Aura Validation Data Centre (AVDC: http://avdc.gsfc.nasa.gov)as overpass files.The TOC from OMI retrieved by the TOMS (OMIT-OMS) and DOAS (OMIDOAS) algorithms, and the TOC retrieved from GOME-2 (GOME2) were compared with Microtops measurements, and showed a high correlation between ground based and satellite-based TOC for the three databases.OMITOMS was seen to be the most similar TOC data series to ground measurements, followed by OMIDOAS and then GOME-2.The comparison can be found in Bilbao et al. (2014).The correlations obtained between the sun photometer and satellite TOC retrieved daily values are shown and compared in Bilbao et al. (2014), with RMSE values being below 2.5 %.In addition, the correlation coefficient between the OMI and Microtops-II TOC daily data of 0.95, the frequency distribution of the differences showed a standard deviation of 6 DU and also a RMSE of 2.2 % were obtained during a two month campaign in Trisaia campaign, Southern Italy (Mateos et al., 2014a).
Water-vapour column data were also obtained from noon measurements of the MODIS instrument (Moderate Resolution Imaging Spectroradiometer) on board the Terra satellite, using MODIS infrared channels.Comparisons of this MODIS product against ground-based water-vapour measurements in a Mediterranean region show that the combined uncertainty of MODIS water vapour increases from 0.38 to 0.52 cm when w rises from 0.5 to 3 cm.This may be found in Román et al. (2014).
The AOD 550 , was also obtained from the MODIS instrument on board the Terra satellite and the results were considered as constant daily values.The uncertainty in the AOD 550 daily values from MODIS was evaluated by the expression proposed by Kaufman et al. (1997).The Ångström exponent values were also observed from daily AERONET (http://aeronet.gsfc.nasa.gov/new_web/index.html)data at the nearest Italian station of Lampedusa (35.52 • N; 12.63 • E; 45 m a.s.l.) (Meloni et al., 2003) and from MODIS.AOD 1020 was calculated using the Ångström exponent and the AOD 550 from MODIS.Comparisons of Microtops II AOD 1020 values and calculated AOD 1020 show a difference of 0.06 % (Bilbao et al., 2014).
The MODIS data used in this work were downloaded from the GIOVANNI application (http://disc.sci.gsfc.nasa.gov/giovanni/overview/index.html) as a spatial the average value in a 0.2 • × 0.2 • square with centre at the Institute for Sustainable Energy, Marsaxlokk, Malta.

Modelling
Parametric models for evaluating solar UVER and horizontal global and diffuse irradiances under cloudless conditions were used; they are called the UVER, Global and Diffuse Model, respectively.The UVER model is based on an analytical expression for clear-sky conditions proposed by Madronich (2007).The Global model was derived from the analytical relation between the horizontal global solar irradiance, G, and the cosine of SZA.It has been observed (Bilbao et al., 2014) that G is a linear function of SZA, while D shows a power function with SZA.The inputs to the UVER model consist of SZA and TOC values.The input to the Global model is SZA and to the Diffuse model the SZA and the AOD.The algorithms of the models are detailed in Table 1.
Data from 15 May to 31 July 2012 were used to estimate the model coefficients and data from 1 August to 15 October 2012 were used for model validation.Statistical indices such as the mean bias, root mean square errors, square fits and correlation coefficients, were used for model validation.

Statistical estimators used
The accuracy of the different models was assessed by means of two widely used statistics: the root mean square error (RMSE) and mean bias error (MBE).The following expressions for RMSE and MBE, as a percentage of the measured average value, are used: where n is the number of data points, x isim is the ith simulated value, x imea is the ith measured value, and x ave is the mean of the measured values.The estimators were calculated for each models.
During the campaign, different atmospheric conditions were observed, particularly events with high aerosols due to   desert dust.The agreement between the measured solar irradiance data and the estimated values was also assessed by determining the percentages of the estimated values that agree with the ±10, ±20 and ±30 % of the measured reference data, which are denoted as W 10 , W 20 and W 30 (%) (Tanskanen et al., 2007).

AMF
In order to evaluate the aerosol effect, the AMF was defined as the relation between the measured solar irradiance and the cloudless sky one under the same atmospheric conditions.
where SR aerosol and SR aerosol−free represent solar irradiance under the corresponding conditions.AMF is the ratio of clear-sky solar irradiance with the atmosphere containing aerosols to that solar irradiance recorded under cloudless conditions and an atmosphere without aerosol, for the same SZA and TOC values.In this work, we compare two adjacent days that verify these conditions (26th with aerosols and 25th aerosol free).

Air mass trajectories
The HYSPLIT model http://ready.arl.noaa.gov/HYSPLIT.php was used to evaluate air mass backward trajectories in order to describe the influence of aerosols on solar irradiance.The input data for a backward trajectory evaluation are: measurement station, considered as the final position of the backward trajectory (35 • 50 N; 14 • 33 E); 72 h duration, concluding at 12:00 UTC; the day analyzed the trajectories at three heights above sea level, 500, 1500 and 3000 m.The dates analyzed are 25, 26, 27 August 2012.

Atmospheric component correlations
Satellite atmospheric component data were compared with ground-based data; frequency histograms were generated before analyzing their influence on solar irradiances.As observed by Koepke et al. (2002), AOD 550 depends on humidity.In this paper, the relation between AOD 1020 and waterapour column was evaluated, and the correlation and its significance was also calculated.In addition, the comparison of TOC measurements from the Microtops and from the OMI-TOMS, OMI-DOAS and GOME-2 sensors was carried out and their uncertainty frequency distributions were plotted so as to verify the agreement between the two data series.The linear fit of calculated AOD 1020 data from MODIS versus AOD 1020 data from Microtops was performed.The correlation coefficient value and the frequency distribution of the difference of its values were evaluated so as to study the agreement between MODIS and Microtops data.

UVI evaluation
The UV index (UVI) was also evaluated in the measurement campaign.To calculate UVI, it is necessary to calibrate the UVER sensor and to have the calibration matrix available, using the facilities at the INTA Laboratory, in our case.The raw signal of the UVB-1 pyranometer is multiplied by a cal-ibration coefficient in order to convert it to erythemal UV-B irradiance.The calibration coefficient depends on SZA and TOC, this coefficient being obtained each time by interpolation in a matrix calibration.Once the raw signal is converted into erythemal UV-B irradiance (W m −2 ), this amount is multiplied by 40 m 2 W −1 in order to convert it into UVI as indicated in WHO (2002).The influence of aerosol desert dust on UVI was also studied.This evaluation also serves as an important vehicle to raise public awareness and to alert people vis-à-vis the need to adopt protective measures when exposed to UV radiation.

Results
After examining the campaign data series, it was observed that on cloudless days, the diffuse horizontal irradiance increases with increasing atmospheric water vapour.The linear fit between AOD 550 from MODIS and w from MODIS shows a correlation coefficient of 0.51, with a confidence interval of (0.160-0.751), at a 99 % confidence level, and a p value of 0.00041.In addition, w uncertainty for the Microtops and MODIS data were evaluated and a standard deviation of the residuals of 0.5 cm was obtained.The linear fit between TOC correlation coefficients between 0.95 and 0.86, while the frequency distribution of the differences is in the range from 4.5 to 7.6 DU and a standard deviation value of 6 DU was obtained; this result can be used as uncertainty of satellite-derived TOC.A similar evaluation can be found in Román et al. (2014).

Estimated vs. ground-based measurements
In this study, three different models for UVER, global and diffuse irradiance under cloudless conditions were calibrated.The data set for the calibration models contain 10 min clear sky data for May, June and July 2012, a total of 684 10 min values.For the evaluation period (August, September and October), a total of 650 data points were used.The solar irradiance clear-sky model expressions and coefficient values are shown in Table 1.The UVER model fits as a SZA cosine and TOC power function.The solar global horizontal irradiance fits as a function of cos(SZA) cosine, and the solar diffuse horizontal irradiance is a function of cos(SZA) and AOD 550 .Similar results have been obtained by Bilbao et al. (2014) when studying the effects of atmospheric components on UVER and horizontal irradiances.A model relating clear sky UV-B irradiance with SZA and TOC was developed by Dubrovsky (2000) at two stations in the Czech Republic.Model validations have been evaluated using scatterplots.Figure 1 consists of three scattered graphs, and the diagonal line represents the ideal match between the estimated and measured values.Table 2 shows the linear regression equations and the correlation coefficient, r, between measured and estimated values.It is observed that UVER and the global solar models show the highest correlation coeffi-cients of 0.995.The diffuse model obtains a lower correlation coefficient of 0.921.
Table 3 shows the model performance, namely mean measured and mean estimated values, the RMSE, MBE deviations and the W i (%) variation coefficients between the measured and the estimated values.The W i coefficients for UVER indicate that 86.7 % of the data exhibit a variation within ±10 % in respect of the measured data, and 96.9 % of the data a variation within ±20 %.From Table 3, it can be deduced that 94.6 % of the global solar data show a variation within ±10 % compared to the measured data.A similar validation coefficient analysis was used by Mateos et al. (2013) and Tanskanen et al. (2007).

Case study: desert aerosols
Figure 2 shows the 10 min evolution and daily value of UVER and UVI 10 min evolution and UVI daily maximum values from the 22 to 28 August 2012.25 and 26 August were classified as cloudless days; however, recorded UVER and calculated UVI decreased on 26 August.In addition, TOC dropped from 296 to 286 DU starting with on 23 August, AOD 550 increased on 26 August only, and UVI increased on 25 August and decreased again on 26 August.While TOC does not show any change (< 1 DU), AOD shows a difference of two units, from 0.10 on 25 to 0.3 on 26 August.Although 27 August was not classified as a cloudless day, it can be observed that the mean daily UVER and AOD values are similar to those for 25 August.
As regards the UVI values, 25 August shows a maximum with 9.8, which drops to 8.6 on 26 August and increases to 9.7 on 27 August.It is also observed that cloudless days with high AOD and water vapour value of 3.8 ± 0.5 cm coincide with a decrease in UVI and UVER, and an increase in daily diffuse horizontal irradiance, which value was 5.05 ± 0.02 MJ m −2 on 26 August 2012.
In addition to this, it is known (Bilbao et al., 2014) that UVER irradiance does not depend on w because w has no absorption bands in the UV range.It was obtained that w reduces G irradiance between −2.44 and −4.53 % cm −1 , and D irradiance increasing with w between 42.15 % cm −1 and 20.30 % cm −1 from low to high SZA angles.It has been observed that D irradiance does not increase with w when AOD 550 is fixed (Bilbao et al., 2015).In addition, the increase in the D component is due to the dust episodes in Malta whose synoptic conditions also increase w and it makes possible the hygroscopic growth of aerosols and as consequence the increases of D irradiance depends on aerosols.Aerosols cause a UVER reduction between 28 and 52 % and a G irradiance reduction between 13 and 41 %.These results show that w impact on irradiances in Malta measurement station is a magnitude lower compared to aerosols, as it can be seen at Bilbao et al. (2014).Under this situation, and in order to study the aerosol impact, it was decided to study air mass trajectories.

Air-mass trajectories
Figure 3 shows the daily backward trajectories on 25, 26 and 27 August 2012.It can be observed that on 25 August, the air masses come from the eastern Mediterranean, that the situation changes on 26 August when the air masses come from the African continent.The air mass at 3000 m crosses Africa at very low altitude and might explain the high aerosol content.The situation changes on 27 August when the air masses come from the Iberian Peninsula.As a result, it is concluded that 26 August was characterized by a dust event from Africa.

AMF
AMF was calculated so as to determine the impact of AOD on horizontal solar radiation.For these evaluations, measurements on 25 August (cloudless sky) were used as reference measurements for an aerosol-free atmosphere.The AMFs were evaluated for erythemal UVER (AMF UVER ), solar global (AMF G ), and diffuse (AMF D ) solar components.An aerosol intrusion from Africa was described on 26 August 2012.In addition, we evaluate the AMF, defined previously, using 25 August measurements as a reference for an aerosol-free atmosphere.In this way, erythemal UVER, global horizontal and diffuse (AMF UVER ), (AMF G ) and (AMF D ) were calculated, respectively.First, from the data series it is observed that on 26 August, Fig. 2, the watervapour column increases to a value of 4.0 cm and the diffuse irradiation to a value of 5.3 MJ m −2 day −1 due to a desert dust intrusion.The AMFs were evaluated for 26 August.The influence of aerosol events on UVER, global and diffuse solar irradiances is compared in Fig. 4. The results show that AMF UVER (purple) values increase from 0.8 to 0.9 after four o'clock GMT.The AMF G does not show changes during the day; its value is around 0.9.On 26 August, the AMF D shows high values in the morning and lower in the afternoon, and AMF D values are higher than one throughout the whole day.The AMF UVER shows lower values than AMF G , indicating a more important aerosol effect on erythemal than on global solar irradiance, similar result having been obtained by Bilbao et al. (2014).In the afternoon, AMF D decreases, probably due to a fall in the aerosol concentration.This result may be confirmed taking into account the air mass origin in Fig. 3 (27 August).The desert dust event observed in Malta with MODIS data of AOD 550 and air-mass back trajectories is also corrobo-rated with the evolution (not shown here) of aerosol properties (AOD), Ångström exponent, among others at the nearby site of Lampedusa island (e.g., Meloni et al., 2003).The agreement of satellite-derived data with ground-based measured data was studied so as to validate satellite data.The TOC daily values from OMI were compared with Microtops measurements and the results showed a high correlation coefficient of 0.95.The frequency distribution of the differences showed a standard deviation of 6 DU.The daily water vapour column measured by Microtops and retrieved from MODIS was compared and the frequency distribution of the differences showed a standard deviation of 0.5 cm.A comparison between AOD 1020 from MODIS and from Microtops showed a correlation coefficient of 0.68 and an RMSE of 0.6.The frequency distribution of said values have a standard deviation of 0.06.
Three different clear-sky models for solar UVER, global and diffuse were calibrated taking into account that solar irradiance values depend on SZA, TOC and aerosols.The RMSE values obtained ranged between 6.6 and 4.5 % for UVER and global solar, while for the diffuse irradiance of 19 %.It was observed that the statistical estimators give better results for solar global irradiance followed by UVER and diffuse.
Aerosol effects have been considered as a case study of dessert aerosol influence.The AMFs were evaluated on 26 August.The influence of the aerosol events on UVER, global and diffuse solar irradiances was compared.The results showed that on 26 August 2012, the AMF UVER and AMF G changed during the day and their values were around 0.9.The AMF D showed high values in the morning and lower in the afternoon.The AMF UVER displays lower values than AMF G , pointing to a more important effect of aerosols on erythemal solar irradiance than on global.A similar result was obtained by Bilbao et al. (2014).The evolution of the AMF D values was due to an increase in AOD observed by daily backward air-mass trajectories.In the afternoon, both AMF D and AMF UVER decreased, probably due to a drop in aerosol concentration.This result may be confirmed taking into account the air-mass origin in Fig. 3 (27 August).
The Supplement related to this article is available online at doi:10.5194/asr-9-147-2015-supplement.

Figure 1 .
Figure 1.Comparison of 10 min measured and estimated UVER (left), global (centre) and diffuse (right) solar irradiances under cloudless conditions during the period May to October 2012 at Marsaxlokk.Red line is the regression fit.Blue line is the 1 : 1 line.

Figure 2 .
Figure 2. 10 min and mean daily values of UVER and UVI, and TOC, AOD 550 and w from satellites at Marsaxlokk in the period 22-28 August 2012.

Figure 3 .
Figure 3. Air-mass backward trajectories at 500, 1500 and 3000 m a.g.l.obtained from NCEP/NCAR reanalysis data and HYSPLIT model at Marsaxlokk in the period 25-27 August 2012.

Table 1 .
UVER, global and diffuse solar horizontal irradiance calibrated models under cloudless conditions, calibration coefficients and their expanded uncertainty, during the period May to October 2012 at Marsaxlokk, Malta.

Table 2 .
Estimation of solar UVER, global and diffuse irradiance values by means of linear regressions for clear-sky conditions in Marsaxlokk, Malta in the period August to October 2012.