Long-Term Tropospheric Ozone Data Analysis 1997–2019 at Giordan Lighthouse, Gozo, Malta

: Long-term data analysis of the hourly ozone volume fractions in the middle of the Mediterranean Seawas carried out covering a period of 22 years. It was noticed that the amount of ozone during this period very rarely exceeded the recommended upper limit value of 80 ppb and that the amount of tropospheric ozone in the area is rather low. Fourier data analysis shows the presence of only a seasonal cycle in ozone concentrations. Statistical analysis of the data is showing a slightly negative trend in ozone concentrations of − 0.46 ± 0.08 ppb/year for average values and a slightly higher negative trend of − 0.54 ± 0.11 ppb/year for the 95 th percentile values. These results obtained through simple linear regression were conﬁrmed using the more appropriate Mann–Kendall test. The possible quadratic trend was not observed for the whole series of data. Air mass trajectories were calculated for those days in the year with the highest pollution, indicating that during those days horizontal air transfer, in most cases, brings the air mass from the North and from Sicily in Southern Italy.


Introduction
For almost a century and especially after World War II, due to a large global population increase, air pollution problems have become of significant importance for the future of humankind.The everlasting negative impacts of various components in the air, especially those emitted by human activity, are of great concern.Therefore, continuous monitoring and analysis of data are of the greatest importance to improving air quality and the quality of life in general.As for the components of the air, the most widely recognised pollutants in the troposphere are various oxides (for example NO x ), ground-level ozone, O 3 and particulate matter (PM) [1][2][3][4][5][6].
Of those, tropospheric ozone is one of the most important.Unlike its stratospheric counterpart, which is of utmost importance for preserving terrestrial life forms [7], in the troposphere it is one of the most dangerous pollutants which causes tissue damage in both plants and animals, including humans [8][9][10].Although not as important as carbon dioxide, ozone is also a greenhouse gas and its influence on the greenhouse effect cannot be omitted [11].Its presence in the troposphere is caused by formation in photochemical reactions (about 90%) and by stratospheric intrusions or diffusion (about 10%) [12].It can also be formed from lightning discharges during thunderstorms.Photochemical reactions that are essential for ozone production include primary pollutants as volatile organic compounds VOC and nitrogen oxides NO x [13,14].Since diesel and heavy fuel oil engines are commonly used as a form of power for ships, and since Malta is located in the middle of the Mediterranean Sea on the major maritime transport routes, this area is extremely prone to pollution by nitrogen oxides [15] despite the low efficiency of NO x emitted by ship engines as the ozone precursors in the marine boundary layer [16,17].
That is in addition to typical Mediterranean meteorological conditions that include high insolation, anthropogenic, and biogenic ozone precursors concentrations, which all favour photochemical ozone production [18][19][20].
Meteorological conditions also have a large impact on ozone concentrations [21,22].Wind speed and direction can influence a horizontal transfer; higher temperatures usually accompanied by higher insolation and humidity can also influence its deposition (i.e., loss of concentration in air) [23][24][25].
Giordan Lighthouse is located on one of the main maritime trade routes in Europe and the world.It also faces Sicily and its volcano, Mount Etna, which is the main natural source of particulate matter and gaseous atmospheric pollutants [26].For a quality mathematical analysis of photochemical pollution, it is very useful to have a long-term data set, which is a fairly common way of analysing ozone trends [27].Together with forecasting methods [28], long-term data analyses gives the best insight into the real situation.By obtaining data in the two-decade time span available, it is now possible to obtain a real insight into the situation in this area.It is known that global tropospheric ozone concentrations in the background increase by an annual rate of around 2% [29].In the Mediterranean area it has already been noticed that trends can have even higher positive rates [30,31], which is another good incentive for this research.On the contrary, in some cases in the wider Mediterranean area [32] or near it [33] during the observed period, even negative rates were observed.
Since vegetation in Malta is very scarce, there are no significant natural sources of ozone in the area or even in the whole state.The aim of this paper is, therefore, to give an insight into the photochemical pollution of the area which is crucial for both the local population as well as the wider Mediterranean and European area, which, given the circumstances, can be assessed from the ozone data.

Experimental
The Giordan Lighthouse GAW station is located on the northern side of Gozo, Malta (36 • 4 N, 14 • 13 E, 167 m above sea level) and started its operations in 1997 [34,35].The station is equipped to monitor meteorological parameters, trace gases and aerosols.For the analysis of ozone, a reactive gas detector has been used.The ozone monitor is calibrated by means of a long path length cell from Landesanstalt für Umwelt Baden-Württemberg (LUBW), Karlsruhe, Germany [15].
A detailed description of the instruments and funding acquired is described by [26].
The trace gases and meteorological data at the Giordan Lighthouse station are recorded at one-minute intervals.On a daily basis, the data are flagged for any effects originating from local anthropogenic emissions.On an annual basis, the data are processed to hourly averages and then uploaded on the EBAS data.Therefore, the data set complied for this publication is available for download from the EBAS database: https://ebas-data.nilu.no(accessed on 16 September 2023.).Meteorological data are available for download from the EBAS database.
Over the span of the 25 years of ground level ozone measurements, the ozone background concentration levels were measured by two commercial ozone analysers-both using the UV absorption technique.In 1997, the department acquired an ozone analyser (Dasibi, 1006 AH).This analyser was donated from an institute in Germany.On an annual basis, the ozone analyser was calibrated against a portable ozone transfer standard (Model O3, UMEG, Karlsruhe, Germany), which itself was checked against an absolute long-path optical cell (Model R-UV-PH, UMEG, Karlsruhe, Germany).In 2010, when the department acquired the ERDF funds, the ozone analyser was replaced with a new analyser (Thermo, 49i, Waltham, MA, USA).On an annual basis, the instrument is calibrated against the ozone laboratory standard (UV-Photometer, Karlsruhe, Germany).
On average, the ozone background level at the Giordan Lighthouse station is around 50 ppb.

Results and Discussion
The monitoring station at Giordan Lighthouse on Gozo Island was established as a Global Atmosphere Watch station in 2001 and upgraded to a fully-fledged station in 2010, using funding from European Regional Development Funds allocated to Malta.This station is a background station located on the isolated part of the island, relatively distant from traffic and other anthropogenic influences (Figure 1).This location for air quality measurements is highly influenced by the maritime traffic and Etna's volcanic particle and gas emissions from the northerly direction.The wind rose plot in Figure 2 was processed with SPSS.A simple technique was used.First, all 25 years of hourly wind data were put in an Excel file.The data were transferred into SPSS.Then, a recoded technique was used in order to assign the wind averages into the respective 36 bins of 10-degree intervals.Then, the frequency of each bin was recoded and converted into a percentage fraction for the whole data set.From this analysis, it was determined that the prevailing wind direction for the Maltese islands is from the NW sectors.The nearby vegetation is rather scarce so the influence of volatile organic compounds (VOC) as precursors of ozone is negligible.

Results and Discussion
The monitoring station at Giordan Lighthouse on Gozo Island was established as a Global Atmosphere Watch station in 2001 and upgraded to a fully-fledged station in 2010, using funding from European Regional Development Funds allocated to Malta.This station is a background station located on the isolated part of the island, relatively distant from traffic and other anthropogenic influences (Figure 1).This location for air quality measurements is highly influenced by the maritime traffic and Etna's volcanic particle and gas emissions from the northerly direction.The wind rose plot in Figure 2 was processed with SPSS.A simple technique was used.First, all 25 years of hourly wind data were put in an Excel file.The data were transferred into SPSS.Then, a recoded technique was used in order to assign the wind averages into the respective 36 bins of 10-degree intervals.Then, the frequency of each bin was recoded and converted into a percentage fraction for the whole data set.From this analysis, it was determined that the prevailing wind direction for the Maltese islands is from the NW sectors.The nearby vegetation is rather scarce so the influence of volatile organic compounds (VOC) as precursors of ozone is negligible.As can be seen from Figure 3 and Figure S1, the majority of hourly ozone volume fractions fall below 50 ppb, and over 80% below 59 ppb, with some years having even lower averages.Distribution also indicates that most of the values are between 35 and 59 ppb indicating rather low daily changes in hourly ozone volume ratios.These low changes, which are even more visible on Figure 4, indicates relatively low ozone pollution in the area in general, despite the influence of heavy traffic, but average values of ozone volume fractions in the area with almost no natural precursors still indicates that artificial precursors are present in significant quantities.Since the reactions for ozone production are photochemical and the Maltese islands are an area with a high sunlight intensity, together with enough precursors some at-site production can be observed.Low daily changes indicate lower production and lower depletion during the average day.Occasionally, there are some episodes of high ozone volume ratios.Those can be attributed mainly to horizontal transport, as discussed later.Still, it has to be said that those values are rather low given the fact that values over 80 ppb are the ones that start causing health problems.As can be seen from Figures 3 and S1, the majority of hourly ozone volume fractions fall below 50 ppb, and over 80% below 59 ppb, with some years having even lower averages.Distribution also indicates that most of the values are between 35 and 59 ppb indicating rather low daily changes in hourly ozone volume ratios.These low changes, which are even more visible on Figure 4, indicates relatively low ozone pollution in the area in general, despite the influence of heavy traffic, but average values of ozone volume fractions in the area with almost no natural precursors still indicates that artificial precursors are present in significant quantities.Since the reactions for ozone production are photochemical and the Maltese islands are an area with a high sunlight intensity, together with enough precursors some at-site production can be observed.Low daily changes indicate lower production and lower depletion during the average day.Occasionally, there are some episodes of high ozone volume ratios.Those can be attributed mainly to horizontal transport, as discussed later.Still, it has to be said that those values are rather low given the fact that values over 80 ppb are the ones that start causing health problems.From Figure 5 we can observe the very common annual periodicity in ozone volume fractions.Even though Malta does not have any long cold periods or "true" winter, it can be seen that during the "winter" months average levels of ozone are lower.This can be attributed to fewer sunlight hours than during summer months; that is, a sunlight intensity seasonality.Still, absolute minimum values of hourly ozone volume ratios are usually lower during summer due to faster photochemical loss and larger vapor concentrations.Similar observations have already been confirmed by [36] albeit for a period that preceded across the Mediterranean Sea show positive contribution to the summertime average ozone volume ratios of up to 12 ppb and up to 40 ppb near busy lines, such as should be the situation in Malta [41].The seasonality of ozone concentrations in Mediterranean area was observed in Greece too, with lower concentrations in the winter and higher in summer [42,43].In comparison with other background stations like Mace Head in Ireland, ozone levels in the Mediterranean area are consistently higher, as can be seen from this analysis and the previously mentioned one [44].From Figure 5 we can observe the very common annual periodicity in ozone volume fractions.Even though Malta does not have any long cold periods or "true" winter, it can be seen that during the "winter" months average levels of ozone are lower.This can be attributed to fewer sunlight hours than during summer months; that is, a sunlight intensity seasonality.Still, absolute minimum values of hourly ozone volume ratios are usually lower during summer due to faster photochemical loss and larger vapor concentrations.Similar observations have already been confirmed by [36] albeit for a period that preceded the one described here.Fourier data analysis was conducted to confirm this further (Figure 6).One year periodicity is indicated by this analysis with the occurrence of a high value first harmonic.However, no higher harmonics can be seen from Figure 6.Sine fit was applied to the dataset prepared for the previously shown Fourier analysis to see where the maximum ozone values occurred.From Figure 7, in addition to Figure 5, maximal ozone concentration occurs around June with the minimum being in December.The corresponding equation that describes the plot is: y = 46.6 ± 0.6 ppb + 6.6 ± 0.5 ppb × sin(x − 2.3 ± 0.3 rad).In this equation x shows a one-year period represented by 0 to 2π rad.The second term represents the first harmonic.In comparison with previous works [37], we can see the lack of the second harmonic which caused our sine function to be reduced by the last term.Also, if this result is compared with some other marine boundary layer stations [37], the annual average is among the highest and amplitude is relatively low which again shows a low daily change.Regarding the offset of x, it is caused by the climate conditions in the region, i.e., early summer condition occurrences.Figure 4  ---- ----  As was indicated earlier, hourly ozone volume fractions are decreasing in the period observed.To further prove that indication, linear regression analysis with the 95% confidence interval has been used.From Figure 8 it can be seen that average values are decreasing with a slope of −0.46 ± 0.08 ppb/year, which is equivalent to an about 1% decrease rate.Similar values are obtained for the summer (April to September) values where a decrease equals slightly less −0.43 ± 0.10 ppb/year.For a better insight, 95th percentile values are also taken into account.The rationale behind this is that high values are of greater concern regarding the environmental and health safety issues.Still, since the maxima can be isolated events which could impact the conclusion significantly, it is better to use slightly lower values, the 95th percentile being a fairly good choice.From Figure 8 it can be seen that higher values are decreasing with an even higher rate, which is −0.54 ± 0.11 ppb/year in both cases-based on the whole year data or just during the summer period.This leads to the conclusion that there is an obvious decrease in hourly ozone volume fractions or, in relation to this, the concentration of ozone in the maritime boundary layer in this area based on near surface measurements.This analysis also shows the importance of longer periods of measurement.In the previous analysis [34], only the first 4 years were considered, during which period a positive trend was observed.But this analysis shows that in the much longer period, and what is more important, closer to the present, such a conclusion is no longer correct.As for the situation in the Mediterranean area, various relevant investigations have been conducted.Photochemical formation of ozone has been reported with mostly slow increase in the Western Mediterranean area in background stations, following the low precursor emissions of NO X and VOC, while the increase in ozone volume ratios was much higher at urban stations [38].On account of the influence of VOC on ozone formation in the maritime boundary layer over the sea, research is conducted directly on the sea with no significant impact on the ozone concentration found [39,40].Ship emissions across the Mediterranean Sea show positive contribution to the summertime average ozone volume ratios of up to 12 ppb and up to 40 ppb near busy lines, such as should be the situation in Malta [41].The seasonality of ozone concentrations in Mediterranean area was observed in Greece too, with lower concentrations in the winter and higher in summer [42,43].In comparison with other background stations like Mace Head in Ireland, ozone levels in the Mediterranean area are consistently higher, as can be seen from this analysis and the previously mentioned one [44].
As was indicated earlier, hourly ozone volume fractions are decreasing in the period observed.To further prove that indication, linear regression analysis with the 95% confidence interval has been used.From Figure 8 it can be seen that average values are decreasing with a slope of −0.46 ± 0.08 ppb/year, which is equivalent to an about 1% decrease rate.Similar values are obtained for the summer (April to September) values where a decrease equals slightly less −0.43 ± 0.10 ppb/year.For a better insight, 95th percentile values are also taken into account.The rationale behind this is that high values are of greater concern regarding the environmental and health safety issues.Still, since the maxima can be isolated events which could impact the conclusion significantly, it is better to use slightly lower values, the 95th percentile being a fairly good choice.From Figure 8 it can be seen that higher values are decreasing with an even higher rate, which is −0.54 ± 0.11 ppb/year in both cases-based on the whole year data or just during the summer period.This leads to the conclusion that there is an obvious decrease in hourly ozone volume fractions or, in relation to this, the concentration of ozone in the maritime boundary layer in this area based on near surface measurements.This analysis also shows the importance of longer periods of measurement.In the previous analysis [34], only the first 4 years were considered, during which period a positive trend was observed.But this analysis shows that in the much longer period, and what is more important, closer to the present, such a conclusion is no longer correct.Values on the y-axis are φ(O3), which stands for ozone volume fractions.Values on the x-axis are the year for which an annual average was calculated.A confidence interval of 95% is applied.
However, linear regression is just one way of analysing trends.To obtain a different insight, a well-known Mann-Kendall statistical test was used.In comparison with the previously obtained results, the decrease rate for the whole period whole year average data obtained by Mann-Kendall's test was almost identical −0.44 ppb/year.In this case, the 95% confidence interval for the slope lies between −0.22 ppb/year and −0.59 ppb/year.An even more rigorous 99% confidence interval shows only negative values and lies between −0.14 ppb/year and −0.65 ppb/year.Some other analyses were carried out.For the 95th percentile whole period whole year data are considered, then the slightly different result of −0.50 ppb/year is obtained here with the 95% confidence interval lying between −0.27 ppb/year and −0.74 ppb/year, while the 99% confidence interval is between −0.16 ppb/year and −0.77 ppb/year.If only seasonal data are taken into account, average values of the slope rate for the whole period are a slightly lower −0.41 ppb/year, with a 95% confidence interval between −0.12 ppb/year and −0.62 ppb/year.The 99% confidence interval lies in the range between 0.00 ppb/year (no change at all) and -0.68 ppb/year.For the more interesting 95th percentile data, the slope rate obtained by Mann-Kendall's test is also a fairly close −0.51 ppb/year.The 95% confidence interval is in the range of −0.20 ppb/year and −0.75 ppb/year, while the 99% confidence interval lies between 0.00 and −0.80 ppb/year.From all these results, a negative trend in ozone concentrations can clearly be seen.Mann-Kendall's test confirmed linear regression analysis for the whole dataset.For other datasets, i.e., seasonal data and the 95th percentile values instead of averages which were chosen as representatives of the highest pollution, a rather high negative trend was also observed, again in accordance with the results of linear regression analyses.All of this indicates a relatively favourable situation for ozone in the marine boundary layer of the troposphere around Malta.Graphical depiction of Mann-Kendall's test can be found in the supplemental data, Figure S3.
In addition to linear regression analyses, sometimes higher order polynomial regression analyses are conducted [37,44,45].Fourier analysis of our data (Figure 6) indicates Values on the y-axis are ϕ(O 3 ), which stands for ozone volume fractions.Values on the x-axis are the year for which an annual average was calculated.A confidence interval of 95% is applied.
However, linear regression is just one way of analysing trends.To obtain a different insight, a well-known Mann-Kendall statistical test was used.In comparison with the previously obtained results, the decrease rate for the whole period whole year average data obtained by Mann-Kendall's test was almost identical −0.44 ppb/year.In this case, the 95% confidence interval for the slope lies between −0.22 ppb/year and −0.59 ppb/year.An even more rigorous 99% confidence interval shows only negative values and lies between −0.14 ppb/year and −0.65 ppb/year.Some other analyses were carried out.For the 95th percentile whole period whole year data are considered, then the slightly different result of −0.50 ppb/year is obtained here with the 95% confidence interval lying between −0.27 ppb/year and −0.74 ppb/year, while the 99% confidence interval is between −0.16 ppb/year and −0.77 ppb/year.If only seasonal data are taken into account, average values of the slope rate for the whole period are a slightly lower −0.41 ppb/year, with a 95% confidence interval between −0.12 ppb/year and −0.62 ppb/year.The 99% confidence interval lies in the range between 0.00 ppb/year (no change at all) and -0.68 ppb/year.For the more interesting 95th percentile data, the slope rate obtained by Mann-Kendall's test is also a fairly close −0.51 ppb/year.The 95% confidence interval is in the range of −0.20 ppb/year and −0.75 ppb/year, while the 99% confidence interval lies between 0.00 and −0.80 ppb/year.From all these results, a negative trend in ozone concentrations can clearly be seen.Mann-Kendall's test confirmed linear regression analysis for the whole dataset.For other datasets, i.e., seasonal data and the 95th percentile values instead of averages which were chosen as representatives of the highest pollution, a rather high negative trend was also observed, again in accordance with the results of linear regression analyses.All of this indicates a relatively favourable situation for ozone in the marine boundary layer of the troposphere around Malta.Graphical depiction of Mann-Kendall's test can be found in the Supplemental Data, Figure S3.
In addition to linear regression analyses, sometimes higher order polynomial regression analyses are conducted [37,44,45].Fourier analysis of our data (Figure 6) indicates and sine fit confirms follows (Figure 7) a strong 1 year periodicity.However, no higher harmonics can be observed.Therefore, linear regression should be the most appropriate.Nevertheless, we made a further calculation with the second order polynomic function (Figure 5), but this has very low coefficients in front of the second order of ppb/month.From this, we can conclude with even more certainty that the trend is linear with a negative slope and that there is no slowing in the decrease.In comparison with the Mace Head station, quadratic regression is not obtained but it must be noticed that even at that station, the maximum of the quadratic function is around the year 2000, which is at the beginning of our data set.Also, quadratic regression might be obtained if only the first half of the set was taken into account since at the beginning a rise in the O 3 concentrations is observed (also in [34]).If calculated, the equation of such regression is given as y = −0.0014± 0.0005 ppb/month x 2 + 0.16 ± 0.02 ppb/month x + 47.6 ± 0.8 ppb/month.A negative trend after that year is observed here too (Figure 9).and sine fit confirms follows (Figure 7) a strong 1 year periodicity.However, no higher harmonics can be observed.Therefore, linear regression should be the most appropriate.Nevertheless, we made a further calculation with the second order polynomic function (Figure 5), but this has very low coefficients in front of the second order of ppb/month.From this, we can conclude with even more certainty that the trend is linear with a negative slope and that there is no slowing in the decrease.In comparison with the Mace Head station, quadratic regression is not obtained but it must be noticed that even at that station, the maximum of the quadratic function is around the year 2000, which is at the beginning of our data set.Also, quadratic regression might be obtained if only the first half of the set was taken into account since at the beginning a rise in the O3 concentrations is observed (also in [34]).If calculated, the equation of such regression is given as y = −0.0014± 0.0005 ppb/month x 2 + 0.16 ± 0.02 ppb/month x + 47.6 ± 0.8 ppb/month.A negative trend after that year is observed here too (Figure 9).Apart from the in situ production of tropospheric ozone in the area, there is always the possibility of horizontal or even vertical air transfer.Since hourly ozone volume fraction values higher than 100 ppb have not been observed during the whole period, it can be safely concluded that vertical transfer did not happen during this time.Vertical transfer of fluids is not very common in general and since the ozone volume fractions are much higher in the lower stratosphere than is indicated with strong spikes of ozone volume fractions measured in the troposphere.Such spikes are absent during the whole period of observation.The observed period is very long, so it can be presumed with a high probability that vertical transfers, like intrusion from the stratosphere, are not very likely here.
Horizontal transfers, which are much more common, can influence the air quality in the area.To acquire a better insight, in every year the days with the highest values of hourly ozone volume fractions were chosen for the trajectories of horizontal transfer analyses.As can be seen from Figures 10 and S4, in almost all cases the air masses came from the north, most commonly over Sicily.Also, since most of the maximum ozone concentrations were measured during the summer, the weather was sunny and warm.This is in accordance with the wind rose, which also shows that most winds are from the northwesterly direction.Since those areas are often highly polluted with anthropogenic and natural sources (Mount Etna), it can be safely concluded that horizontal air transfers do indeed have a negative impact on the air quality in the area.It is worth mentioning that in some Apart from the in situ production of tropospheric ozone in the area, there is always the possibility of horizontal or even vertical air transfer.Since hourly ozone volume fraction values higher than 100 ppb have not been observed during the whole period, it can be safely concluded that vertical transfer did not happen during this time.Vertical transfer of fluids is not very common in general and since the ozone volume fractions are much higher in the lower stratosphere than is indicated with strong spikes of ozone volume fractions measured in the troposphere.Such spikes are absent during the whole period of observation.The observed period is very long, so it can be presumed with a high probability that vertical transfers, like intrusion from the stratosphere, are not very likely here.
Horizontal transfers, which are much more common, can influence the air quality in the area.To acquire a better insight, in every year the days with the highest values of hourly ozone volume fractions were chosen for the trajectories of horizontal transfer analyses.As can be seen from Figure 10 and Figure S4, in almost all cases the air masses came from the north, most commonly over Sicily.Also, since most of the maximum ozone concentrations were measured during the summer, the weather was sunny and warm.This is in accordance with the wind rose, which also shows that most winds are from the northwesterly direction.Since those areas are often highly polluted with anthropogenic and natural sources (Mount Etna), it can be safely concluded that horizontal air transfers do indeed have a negative impact on the air quality in the area.It is worth mentioning that in some cases there is an observable horizontal transfer from north Africa as well.It is not uncommon that such air masses bring Saharan sand and warmth when they occur.Such occurrences influence particulate matter concentrations much more than ozone concentrations.
sphere 2023, 14, x FOR PEER REVIEW cases there is an observable horizontal transfer from north Africa as we mon that such air masses bring Saharan sand and warmth when they rences influence particulate matter concentrations much more than ozo Figure 10.Air trajectories calculated for the Giordan Lighthouse observation sit ware and kinematic model [46] on the most polluted day during the observati 2003.Trajectories were calculated for 5 days backwards from the date mentione for the most polluted days per year can be found in the supplemental material Wind direction and ozone concentration were compared to give shore and onshore wind influences (Figure 11).The station is located in with northern winds coming from the sea and southern winds coming Gozo itself and further south from the island of Malta.It can be see ozone volume fractions are connected with the northern winds (betw and lower ozone volume fraction values with southern winds (90° to proves the conclusion made from the wind rose and air trajectories t transferred to the site from the north.In addition, ozone is less depleted over the sea due to a lack of photochemical reactions.The average diff high but is in accordance with the relatively low diurnal production an Wind direction and ozone concentration were compared to give an insight into offshore and onshore wind influences (Figure 11).The station is located in the north of Gozo with northern winds coming from the sea and southern winds coming from the island of Gozo itself and further south from the island of Malta.It can be seen here that higher ozone volume fractions are connected with the northern winds (between 270 • and 90 • ) and lower ozone volume fraction values with southern winds (90 • to 270 • ).This further proves the conclusion made from the wind rose and air trajectories that most ozone is transferred to the site from the north.In addition, ozone is less depleted when transferred over the sea due to a lack of photochemical reactions.The average difference is not very high but is in accordance with the relatively low diurnal production and depletion cycle.

Figure 1 .
Figure 1.Location of the monitoring station at Giordan Lighthouse in Gozo, Malta (marked with a drop point).Figure 1. Location of the monitoring station at Giordan Lighthouse in Gozo, Malta (marked with a drop point).

Figure 1 .
Figure 1.Location of the monitoring station at Giordan Lighthouse in Gozo, Malta (marked with a drop point).Figure 1. Location of the monitoring station at Giordan Lighthouse in Gozo, Malta (marked with a drop point).

Figure 2 .
Figure 2. Wind rose based on the whole year data for the Giordan Lighthouse station showing wind direction frequencies.

Figure 2 . 15 Figure 3 .
Figure 2. Wind rose based on the whole year data for the Giordan Lighthouse station showing wind direction frequencies.Atmosphere 2023, 14, x FOR PEER REVIEW 5 of 15

Figure 3 .
Figure 3. Distribution of the hourly average ozone volume fractions from 1997 to 2019.Hourly averages of the ozone volume fraction are distributed in sets with the range of 5 ppb and shown with vertical columns.The red line shows the percentage of hourly averages of ozone volume fractions considered until the given set.

Figure 4 .
Figure 4. Box and whiskers plot of the hourly ozone volume fractions with statistical values from the whole dataset.Maxima and minima (upper and lower extremes) are shown with a dot, the upper whisker shows the 90th percentile, the upper box line shows the 75th percentile, the inner box line shows the median, the lower box line shows the 25the percentile and the lower whisker line shows the 10the percentile.φ(O3) stands for hourly ozone volume fractions.Values on the x-axis are the hours in the day when the data were obtained.

Figure 4 .
Figure 4. Box and whiskers plot of the hourly ozone volume fractions with statistical values from the whole dataset.Maxima and minima (upper and lower extremes) are shown with a dot, the upper whisker shows the 90th percentile, the upper box line shows the 75th percentile, the inner box line shows the median, the lower box line shows the 25the percentile and the lower whisker line shows the 10the percentile.ϕ(O 3 ) stands for hourly ozone volume fractions.Values on the x-axis are the hours in the day when the data were obtained.

15 Figure 5 .
Figure 5. Monthly averages of ozone volume fractions φ(O3) for the whole period of observation.φ(O3) stands for ozone volume fractions.The black line represents the linear regression trend in ozone volume fractions.The dotted line represents quadratic polynomial regression with the corresponding equation.Units given in the equation are ppb/month on the appropriate polynomial order.Intercept points represent the value for the year 1997 as the first year of measurements.

Figure 7 .
Figure 7. Sine fit of O3 during the whole observed period.A curve shows the least square regression of the first harmonic shown in Figure 6.

Figure 7 .
Figure 7. Sine fit of O 3 during the whole observed period.A curve shows the least square regression of the first harmonic shown in Figure 6.

Atmosphere 2023 , 15 Figure 8 .
Figure 8. Linear regression based on the annual average of ozone volume fractions, marked with dots, for the whole year (blue) and for the summer season (grey) and based on the 95th percentile of the hourly ozone volume fractions for the whole year (red) and for the summer season (orange).Values on the y-axis are φ(O3), which stands for ozone volume fractions.Values on the x-axis are the year for which an annual average was calculated.A confidence interval of 95% is applied.

Figure 8 .
Figure 8. Linear regression based on the annual average of ozone volume fractions, marked with dots, for the whole year (blue) and for the summer season (grey) and based on the 95th percentile of the hourly ozone volume fractions for the whole year (red) and for the summer season (orange).Values on the y-axis are ϕ(O 3 ), which stands for ozone volume fractions.Values on the x-axis are the year for which an annual average was calculated.A confidence interval of 95% is applied.

Figure 9 .
Figure 9. Monthly averages of ozone volume fractions φ(O3) for the period from 1997 to 2008 in which maximum values of ozone volume fractions occurred.The dotted line represents quadratic polynomial regression with the corresponding equation.Units given in the equation are ppb/month on the appropriate polynomial order.Intercept points represent the value for the year 1997 as the first year of measurements.

6 φFigure 9 .
Figure 9. Monthly averages of ozone volume fractions ϕ(O 3 ) for the period from 1997 to 2008 in which maximum values of ozone volume fractions occurred.The dotted line represents quadratic polynomial regression with the corresponding equation.Units given in the equation are ppb/month on the appropriate polynomial order.Intercept points represent the value for the year 1997 as the first year of measurements.

Figure 10 .
Figure10.Air trajectories calculated for the Giordan Lighthouse observation site using METEX software and kinematic model[46] on the most polluted day during the observations-14th of August 2003.Trajectories were calculated for 5 days backwards from the date mentioned above.Trajectories for the most polluted days per year can be found in the Supplemental Material.