Variation of Aerosol Optical Properties and Trends over Cluj-Napoca, Romania, Based on 10 Years of AERONET Data and Satellite-Based Observations

1. Introduction

According to the Intergovernmental Panel on Climate Change (IPCC), aerosols play an important role in Earth’s climate system, due to the absorption and scattering of the solar radiation, and interactions with the clouds [1]. These processes are still sources of uncertainty in climate modelling and the magnitude of aerosol radiative forcing [2]. In addition, the atmospheric aerosols may affect human health [3], [4], and pose socio economic risks [5]. Regional studies are essential to assess the aerosol trends, as the optical parameters vary from one region to another.

In present, some ground networks are monitoring the aerosols properties, the most developed one being the AERONET network, with hundreds of stations at global scale [6]. Other major networks are SKYNET (MEXT Sky Radiometer Network) [7] and the Global Atmospheric Watch/Precision Filter Radiometer (GAW/PFR) [8].

AERONET data are accurate and considered a benchmark for validating other aerosol measurement data. AERONET products are accessible, making them a valuable resource for researchers and atmospheric scientists to analyze and interpret data for various applications. In Asia, Tu et al. (2021) used AERONET data to determine aerosols characteristics over four observation sites in East China, highlighting the frequency of urban pollution from fossil fuel burning during a period of 20 years [9]. Yu et al. (2022), determined the dominant aerosol types in Hong Kong and their interactions with meteorological factors by using both AERONET data and satellite-based observations from 2006 to 2021 [10]. Natural, large-scale events of aerosols intrusions have been researched by Sun et al. (2022), who studied two extreme dust events over North China in March, 2021, by using AERONET in situ observations and CALIPSO satellite data to obtain the distribution characteristics of dust aerosol optical properties along their transport pathways [11]. Additionally, Dementeva et al. (2022) analyzed AERONET data collected over a 10-year period and found that there was a significant increase in smoke aerosols during the summer months, which was caused by large-scale wildfires in the boreal forests surrounding Lake Baikal [12]. In Southeastern Europe, long-term aerosol trends have been researched using AERONET datasets. Carstea et al. (2019), analyzed the climatology of aerosol optical and microphysical properties over Romania, based on 9 years of AERONET data to highlight the efficiency of EU regulations on particulate matter emissions in Bucharest [13]. Evgenieva et al. (2022), used a two-year AERONET dataset to highlight the main characteristics and transport models of aerosol loads over Sofia, Bulgaria, quantifying the high content of urban aerosols and rare occurrences of desert dust and biomass-burning aerosols [14]. In Greece, Raptis et al. (2020), studied the aerosols seasonality and trends over Athens, using a decade-worth of AERONET and satellite data to underline an increase in aerosol loads during spring and summer months [15]. Measuring campaigns in Thessaloniki and data from AERONET observations, carried out by Voudouri et al. (2022), were used to study the intrusions of biomass-burning aerosols over Greece [16]. AERONET long-term measurements were also used to determine the climate impacts of aerosols. Markowicz et al. (2021), studied the climate interaction trends of aerosols distributed over Poland and their effects on incoming radiation fluxes, using a 10-year AERONET dataset [17]. Damiano et al. (2022), characterized the columnar aerosol optical and microphysical properties to determine the prevailing aerosol type in the Naples Mediterranean area, using an AERONET dataset from a 5-year period [18]. AERONET and satellite datasets for a period of 16 years were also used by researchers to determine the annual variability of aerosol intrusions episodes in Morocco, and the location’s susceptibility to desert dust transport on a seasonal basis [19]. Timpu et al (2020) analyzed the tropospheric dust and associated atmospheric circulations over the Mediterranean region using modeling and AERONET data [20].

The Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm [21], applied to the Moderate Resolution Imaging Spectroradiometer (MODIS) level 1B land-surface radiances generated for a global 1 km resolution sinusoidal grid, is widely used in air quality and epidemiological studies due to its high 1 km resolution and lack of "urban" bias. The MAIAC AOD has been evaluated over a few geographical regions such as North America [22,23] and South America [24], South Asia [25], and Europe ([26]–[28], the Black Sea, arid areas of the Dead sea [29], the alpine region [30] and Australia [31]. Some studies suggest that the algorithm may be systematically underestimated in desert areas of Western China and exhibit a positive bias at low-moderate aerosol loading in eastern China [32,33]. Additionally, the MAIAC AOD product has been found to perform well in South Asia, outperforming other AOD products [34]. The algorithm has also been used to derive surface particulate matter concentration over various regions such as the USA [35], Mexico City [36], Italy [28] and Israel [37]. These studies have highlighted the ability of the MAIAC AOD product to capture spatial variations of PM2.5 with higher accuracy than other methods, as well as its improved correspondence with ground-based measurements in certain regions. Despite the extensive use of MAIAC AOD, few studies have been conducted to verify its accuracy and robustness on a global scale. The uncertainty of the MAIAC AOD retrievals was found to be heavily dependent on satellite and solar geometries, aerosol types, particle size and aerosol loading. The algorithm performs particularly well over densely vegetated areas, bright surfaces and when retrieving smoke AOD [34], [38], [39].

The main objective of the study is to comprehensively investigate the variability and trends of aerosol optical properties over Cluj-Napoca, Romania, using AERONET data. This involves characterizing the temporal distribution of aerosol optical properties, assessing data availability, comparing the results obtained with other data sources, and identifying the factors contributing to the observed variability and trends. The study aims to provide valuable insights into the dynamics of atmospheric aerosols over the study area, as well as fill the knowledge gap regarding the distribution of aerosols in the region (only site on a 200 km radius), thereby making a novel contribution to the field of aerosol research in Eastern Europe. In order to achieve the main goal of this study, several complementary objectives were pursued. These included conducting a detailed investigation of the seasonal patterns of aerosol optical properties, examining the relationship between aerosol optical properties and meteorological parameters, and assessing the accuracy and precision of the MODIS MAIAC data.

The MODIS MAIAC AOD product is known to improve the retrieval over urban areas while subsequently reducing the scale down to 1 km [21]. These improvements over similar algorithms such as Deep Blue (DB) and Dark Target (DT) make it a viable candidate for analyzing urban AOD and aerosol climatology studies. Since validation efforts in Eastern Europe have not been reported, this study covers this knowledge gap while also comparing the robustness of the retrieval algorithm to other urban locations in different parts of the globe.

2. Materials and Methods

2.1. Site description

Cluj-Napoca is a city of aprox. 400 000 inhabitants located in the North-Western part of Romania, on Someș river valley surrounded by hills. The measurement site is located at the Faculty of Environmental Science and Engineering (Lat: 46.76833° N, Long: 23.55139° E, elevation: 405 m). The measurement station is part of ACTRIS (Aerosols, Clouds and Trace gases Research Infrastructure, https://www.actris.eu/) network and conducts aerosol measurements since 2010. The Cluj-Napoca AERONET measurement station can be considered representative for the whole Transylvanian territory, and also for the Pannonian Basin, given the limited coverage in Hungary, Austria, and South-Western Ukraine [40]. The main sources of local emissions are represented by heavy traffic and heating facilities during winter, with few emissions related to industrial facilities. The positioning of the city in a valley with neighboring high topography may cause frequent temperature inversions during the winter season and contributes to the accumulation of air pollution in the city. In addition, long-range transported aerosols such as Saharan dust [40] and biomass burning [41] may also influence the area. A detailed view of Cluj-Napoca and the surrounding area are shown in Figure 1.

2.2. AERONET data

The CLUJ_UBB AERONET station (46.768 ° N; 23.551° E) is located in Cluj-Napoca city, on the rooftop (20 m above ground level) of the Faculty of Environmental Science and Engineering, Babeș Bolyai University at an elevation of 405 m. The instrument is a Cimel Electronics (Paris, France) Automatic Sun Tracking Photometer CE 318, which performs direct Sun and sky radiance measurements at nine band-pass filters between 340 and 1640 nm. We analyzed data from the start of the measurements in July 2010 until October 2020. Gaps in the time series are expected due to the annual calibrations demanded by the AERONET network protocol. The calibrations were performed at Laboratoire d’Optique Appliquée (LOA) at Université de Lille, France. The measurements are processed centrally and are available on Aerosol Robotic Network (AERONET) [42] website. Aerosol optical depth (AOD) data are produced for three data quality levels: Level 1.0 (unscreened), Level 1.5 (cloud screened), and Level 2.0 (cloud-screened and quality-assured). Besides the gaps from the calibration of the instruments (2 or 3 months for every calibration), other gaps in time series may be present because of the cloudy conditions and/or because of the instrument malfunction.

In this study, AERONET version 3 retrievals have been used at level 2.0 for direct Sun (AOD and Ångström exponent) in order to have the highest quality of data. For the Single Scattering Albedo (SSA) inversion retrievals, level 1.5 was used. Due to the very strict AERONET level 2.0 criteria, these type of data is limited in many areas, including Cluj-Napoca. Nevertheless, this approach should be performed taking into account the higher uncertainty of this product which could be up to 0.04 higher than level 2, depending on the AOD and Solar Zenith Angle. A similar approach was carried out in [43] as this practice is relatively common in climatology studies.

Further on, an aerosol classification was conducted based on the thresholds described by Dubovik et al. [44] using both direct measurements and inversion products. Raptis et al. employed a comparable methodology in their previous investigation [45]. It should be noted that this classification has the role of highlighting the dominant aerosol types and their prevalent proportions, as the columnar atmospheric mixture over Cluj-Napoca is highly unlikely to have only one type of aerosols [41].

Following the classification process, a statistical analysis was performed to determine the distribution of monthly observations for each aerosol type, with the aim of ascertaining any potential seasonal patterns. A Pearson correlation coefficient (noted "λr" not to be confused with the correlation coefficient presented in the MODIS data comparison section) was calculated in order to measure the linear correlation between the monthly distribution for each class of aerosols and the monthly distribution for the total number of observations.

2.3. MODIS data

The Moderate Resolution Imaging Spectroradiometer (MODIS) sensors on board the Terra and Aqua satellites have been providing columnar aerosol properties since the year 2000 and 2002 respectively [46] [47]. The two satellites perform an equatorial crossing at 10:30 (Terra) and 13:30 (Aqua) local time. MODIS sensor makes use of 36 spectral bands at spatial resolutions of 250, 500 and 1000 meters, while scanning a 2330 km wide swath to provide near-daily global coverage of the Earth. The MODIS aerosol retrieval algorithms have undergone several updates resulting in four operational algorithms: Dark Target (DT) over land [47], DT over ocean [48] Deep Blue (DB) for bright surfaces [49] and land surface [46] and Multiple Angle Implication of Atmospheric Correction (MAIAC) algorithm. The standard MODIS aerosol products (DT and DB) are provided with two spatial resolutions; one of 10 km at the nadir and another of 3 km [50], while MAIAC retrieves high resolution data at 1 km [51,21]. This finer spatial scale provides an excellent opportunity for aerosol research at city levels. In this analysis, we used collection 6 (C6) data products of the MAIAC AOD retrieval algorithm.

MODIS MAIAC aerosol retrieval algorithm

The Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm utilizes a physical atmosphere-surface model that is obtained from measurements with minimal assumptions [21]. The algorithm is applied to the MODIS level 1B land-surface radiances, which are generated for a global 1 km resolution sinusoidal grid. To simultaneously retrieve aerosol and surface parameters, the algorithm utilizes multi-angle information from time series of MODIS observations for up to 16 days for a given pixel at a resolution of 1 km [51]. The surface characterization in MAIAC is achieved through the derivation of spectral regression coefficients (SRC) that relate the surface BRF in the blue (470 nm), green (550 nm), and shortwave infrared (2130 nm) bands of MODIS. The MAIAC cloud mask is combined with the detection of absorbing aerosols (smoke or dust) allowing MAIAC to characterize most of the aerosol emission sources at 1 km resolution, including high intensity plumes, without masking them as clouds. By characterizing and storing unique surface signatures for each 1 km grid cell, MAIAC is able to detect clear-sky, clouds, or snow cover pixels [52] and enhance the accuracy of atmospheric correction and AOD retrieval in partially cloudy days [21]. The C6 version of the MAIAC algorithm has also improved the estimation of surface reflectance at a pixel scale, reducing the "Blockiness" of AOD retrievals caused by random SRC bias [21]. MAIAC selects among three aerosol types including two absorbing coarse and fine mode aerosols (dust and smoke) and one less absorbing background aerosol type, which is regionally derived using one of eight prescribed aerosol models based on long-term climatology of aerosol properties obtained from AERONET [21].

Collocation and validation approach

To validate the MODIS MAIAC AOD product we used the level 2.0 direct-sun AOD AERONET version 3 as reference [53]. The MAIAC AOD provides data at 550 nm and 470 nm while the closest AERONET Direct Sun AOD is at 440 nm and 500 nm. In this study, we estimated the AERONET AOD at 550 nm by linear interpolation, using the AERONET Angstrom Exponent, from the 440 nm and 675 nm wavelength pair [54]. The uncertainty in the AERONET AOD is dominated by sensor calibration and is approximated at 0.01 – 0.02 [54]. The sampling interval is dependent on site specifications and is achieved between 5 and 15 minutes, in cloud-free conditions.

To accurately compare satellite AOD retrievals with AERONET AOD measurements, it is crucial to account for atmospheric variability and temporal offset in spatio-temporal collocation. Averaging two or more AERONET measurements taken around the satellite overpass time is recommended to ensure reliable results [55], [56]. Thus, AOD validation typically averages AERONET data within a ± 30–60 min temporal window matching satellite retrievals over a ± 25 x 25 km area; however, this criterion should be adjusted for the fine spatial resolution of the MODIS MAIAC product. Several studies have reported validation efforts of the MAIAC AOD product, utilizing various temporal collocations intervals. A temporal window of ± 60 minutes of mean AERONET AOD was reported in [34], [57]. The most common interval reported in the literature is ± 30 minutes as seen in [23], [32], [38], [39]. [58] suggest the optimal temporal interval of ±15 min which is also reported in [23], [38]. However, this comes at the expense of lesser matches with minimal improvements in overall statistics. Since the validations in this paper were conducted over a 10-year period, we selected a more stringent temporal interval of ±15 min with two or more AERONET retrievals. Regarding the spatial collocation, values range from single pixel overlap of 1 x 1 km boxes ([38], [39], [59]), up to 25 x 25 km boxes [57] or 25 km radius [32]. We selected a 5 km radius in this study, similar to [23], [38] as the best representation of the urban area around the AERONET location in Cluj-Napoca. Matchups were discarded if less than < 20 % of pixels were retrieved in the selected 5 km radius.

For the C6 MODIS MAIAC retrievals, the highest quality data was selected for validation. The MAIAC AOD uncertainty is estimated at ± (0.05 + 0.1⋅AOD) [21] however, regional influences were reported ranging from ± (0.05 + 0.15AOD) over South Asia [34] and China [60], to ± (0.05 + 0.05AOD) over South America [57]. We decided to evaluate the reported uncertainty estimate of ± (0.05 + 0.01⋅AOD) by [21]. To achieve a rigorous validation effort we opted for the method described in [58].

The total expected discrepancy (ED) between a single matchup MODIS MAIAC and AERONET, accounts for both datasets uncertainties

ED = sqrt (ε²_S + ε²_A)

(1)

where ε2S is the MODIS MAIAC uncertainty while ε2A is the AERONET uncertainty. Based on the expected discrepancy we can then define a normalized error ΔN as the ratio of the actual error to the ED,

ΔN = (AOD_MODIS – AOD_AERONET)/ED

(2)

In the ideal case, EA << ES, the shape of ΔN is dominated by the uncertainty and errors on the satellite retrieved AOD. If the uncertainties on satellite and reference AOD have been calculated appropriately and the sample is sufficiently large, then the normalized error ΔN should approximate a Gaussian distribution with mean 0 and variance 1. If ΔN is normally distributed, 68.3 % of values should fall within the range [-1, +1]. If the fraction is smaller, then uncertainties are underestimated, if it is larger, then uncertainties are overestimated. Additionally, the average uncertainties, correlation coefficient (r), mean bias (Bias), root mean square error (RMSE).

3. Results

3.1. Meteorological parameters

Meteorological parameters such as temperature, relative humidity, and wind speed and direction, have important effects on the evolution of aerosols [61]. In order to analyze the meteorological conditions in the study period, daily data received from the in- situ DAVIS Vantage Pro2 (Hayward, USA) weather station was used. For wind direction and wind speed, we used only the daytime data during the sun-photometer measurements.

Due to the temperate continental climate, the monthly average air temperature in Cluj-Napoca can vary greatly, ranging between – 2.17 °C in January and 21.1 °C in July, with a mean annual air temperature of 10.3 °C (2010-2020) (Figure 2). The official national meteorological station annual mean temperature in Cluj-Napoca for the same period was 10.1°C [62]. Comparing it with the value recorded over a longer period (between 1961-2019), of 8.8°C [63], a significant increase in the local temperature was observed in the investigated period.

The relative humidity, which is generally inversely proportional with air temperature shows the highest values in the cold season, mainly between November and January (around 90 % in average), while the lowest values are specific to the March-April period (< 70%). Monthly, the mean values of the relative humidity vary between 92 % in December and 68 % in April, while the annual mean is 79 % (Figure 2).

The wind direction in Cluj-Napoca is largely influenced by the local topography, being rather consistent with the Someșul Mic Valley orientation (generally from Southwest to Northeast). In addition, the wind direction pattern can have significant variations along the seasons. During the cold season (October - March) when the regional influence of the East-European anticyclone is increasing, the daytime winds tend to blow from the East-North-East (36.8 % of all observations in the range of 45-105°), bringing polar-continental air masses over Romania. In the warm season (April - September) the winds are more likely to come from North-West (23.8 % of all observations in the range of 285-345°), although the East-North-East directions have maintained their high prevalence. High percentages are also specific to winds that come from the South-West (Figure 3). The average daytime wind speed is low in the study area, being higher in the warm season (1.94 m/s) and lower in the cold season (1.86 m/s). Analyzing the wind speed (m/s) against wind direction it can be observed that the highest average wind speeds are clearly related to the northwesterly directions that usually bring unstable polar-maritime air masses (Figure 3). Due to the sheltering conditions induced by Western Carpathians and Somesul Mic Valley, the atmospheric calm has a high frequency (47.2%), having minimum values in the warm season and on daytime conditions, thus producing more often pollution episodes.

3.2. Seasonality

As one of the objectives of this study is to determine the climatology of aerosols in the Cluj-Napoca region, the analysis encompassed 44,070 observations collected between 2010 and 2020. Figure 4 illustrates the monthly distribution of the data points during the entire measuring period. It can be observed that the quantity of data obtained is correlated with the number of sunny days since measurements rely on direct solar irradiance. Thus, most observations were made in the summer months (mainly July and August), December being the month with the lowest number of observations (1169). This proportionality is apparent if we consider the simulated historical climate and weather data for Cluj-Napoca available on the meteoblue® website [64]. This data shows that July and August are the months with the sunniest and partly cloudy days (over twenty-four days per month).

A further aspect that influenced the availability of data during the analyzed period is associated with periodic calibration procedures that imply data gaps in the respective period.

The monthly mean variation of AOD at 500 nm is presented in Figure 5. Retrievals at this wavelength show peak monthly AOD values in July-August (0.23 and 0.21 ± 0.11, respectively), while the minimum values are recorded in December and January (0.15 ± 0.09). In the absence of precipitation, the low scavenging of atmospheric aerosols is linked to the warm season maximum. In addition, aerosols resulted from wildfires caused by severe drought events in Eastern Europe have influenced the aerosol load in the atmosphere in the summer months. On the opposite, higher precipitation values in the colder months played a dominant role in below-cloud aerosol wet deposition. Also, the behavior of the population in Eastern Europe plays a role in contributing to the total aerosol load in the atmosphere in the cold season by generating controlled fires (biomass burnings) in agricultural regions [65]. The majority of these fires occur during the pre-spring (February) and late autumn (November) periods.

Figure 6 shows the monthly mean variation of Ångström exponent (440-870nm). Upon initial observation, it becomes apparent that fine particles prevail throughout the year as indicated by AE values surpassing 1. The overall mean value of the AE for the entire measurement period is 1.5 ± 0.29, with a maximum of 1.62 ± 0.26 in August, and a minimum of 1.33 ± 0.39 in April. The April minimum can be mainly attributed to the presence of larger aerosol particles such as mineral dust. This might be a consequence of the Saharan dust intrusion events [65]–[68].

Single Scattering Albedo at 440 nm is presented in Figure 7. Despite utilizing level 1.5 data for this parameter, just over 5000 observation points were generated during the analyzed period. The overall mean SSA value at 440 nm is 0.92 ± 0.05, with a maximum of 0.96 ± 0.03 in June, and a minimum of 0.87 ± 0.07 in December. The lowest values are typically observed during the cold season, particularly in October, November and December, when the atmosphere contains higher concentrations of absorbing aerosols due to atmospheric conditions and increased use of fossil fuels for heating and transportation. This can result in elevated levels of black carbon in the atmosphere.

The aerosol classification, which was used in similar studies [44], [45], is presented visually in Figure 8. As mentioned in the Methodology section, this classification is only meant to highlight the prevailing type(s) of aerosols and their proportions based on previously defined threshold values for each class. According to this classification, we observed higher proportions of continental (19.1%), polluted (29.4%), and mixed aerosols (32%). In contrast, the least frequently detected as dominant aerosol classes are those resulting from biomass burning (7.9%) and marine sources (4.5%). Although marine and continental aerosols are prevalent in cases of low AOD, the scarcity of observed marine aerosols can be attributed to the considerable distance from the sea. Additionally, the physical properties of marine aerosols are not conducive to their long-range transport. Mineral dust comprises 7.1% of the overall observations, this proportion being influenced by the seasonal occurrence of mineral dust intrusions.

The graphics presented in Figure 9 depict the monthly distribution of observations for each aerosol class to highlight a potential seasonal character for different types of aerosols. While mixed and polluted aerosols appear to be more prevalent in the summer months, their high correlation coefficient with the total number of observations (λr=0.86 for mixed and λr=0.91 for polluted) suggests that they are present throughout the year. It is possible that their apparent prevalence in the summer months is due to proportional representation.

The distribution of monthly observations for continental aerosols reveals a shift towards the months of August through October. Although the high number of observations in August can be attributed to a relatively elevated correlation coefficient (λr=0.71) the prevalence of such observations in early autumn can be a result of increased agricultural activities, like harvest and plowing in addition to natural processes related to vegetation decay. In addition, weather patterns (relatively calm weather conditions) facilitate the buildup of aerosols from previously mentioned sources.

For biomass burning aerosols, the high frequency of observations during summer months is attributed to the increased occurrence of wildfires over the last decade, which have been exacerbated by episodes of drought. Some examples of such events are the 2010 Russian wildfires [69], the 2015 Ukrainian wildfires [70], [71], and the 2019 Siberian wildfires [72]. Biomass burning aerosols from the specific Eastern European agricultural practices of burning the farmland are also observed in smaller fractions during the months of March-April and November. During the spring months, the atmosphere contains a significant fraction of coarse aerosols (such as mineral dust), which can contribute to a lower number of observations where only the fine fraction is detected.

Similar to specific events that were highlighted in previous studies on dust intrusions over Cluj-Napoca [65], [67], [68], [73], this statistical analysis confirm that dust aerosols are notably prevalent during the months of April, May, and June. This seasonality is mostly determined by atmospheric circulation patterns.

The atmospheric circulation patterns previously mentioned facilitate also the transport of marine aerosols, particularly in May, as air masses move from the southwest (Adriatic Sea) towards the continental region of Eastern Europe. The small value of the correlation coefficient (λr= 0.05) clearly shows the seasonal character of these aerosols in the atmospheric column above Cluj-Napoca.

3.3. MODIS MAIAC validation against AERONET

Figure 10 shows the linear correlation between satellite (MODIS MAIAC) and ground based observations (AERONET) over Cluj-Napoca, Romania. The slope and intercept of the linear fit are known as good indicators of the robustness of aerosol model assumptions and surface reflectance estimates as part of the retrieval algorithm [32], [47]. The slope (0.73) and negative bias (-0.03) indicate a slight underestimation of MAIAC AOD in overall low AOD conditions averaging 0.17 observed by AERONET and 0.14 retrieved by MAIAC. The correlation coefficient r = 0.77 shows fairly good agreement between the two datasets while the RMSE indicate a low value of 0.068.

The distribution of normalized errors is also presented in Figure 10. The expected error (EE) envelope contains the absolute (0.05) and relative (0.10) uncertainty estimates as described in [21] and should account for the majority of error sources. At the Cluj-Napoca site, 71% of errors fall within EE while 23% fall below and 6 % fall above this threshold. This indicates that the MAIAC retrieval seems to overestimate the uncertainties in low AOD conditions. Since both AOD and error values were relatively small, we assessed the contribution of the AERONET uncertainty to the overall distribution of normalized errors. As such there was only a 3% increase in the below threshold suggesting that the MAIAC error budget is somewhat pessimistic. Thus, at least for this location, a slight adjustment to the first coefficient of the EE expression is warranted.

The cumulative distribution function (CDF) in Figure 11 suggests that the normalized errors are randomly distributed with no apparent deviation from an ideal curve. The bottom row shows 12 equally divided bins of the expected discrepancy vs the 1 σ retrieval errors and twice the expected discrepancy vs the 2 σ retrieval errors as suggested in [58]. The MAIAC uncertainties seem to be overestimated on the low end and slightly underestimated on the high end with the majority of values being clustered close to the 1:1 line. When assessing the 2 σ retrieval errors, the overestimation of uncertainties at Cluj-Napoca is slightly more apparent although bin values seem to cluster more together suggesting that the retrieval errors are roughly Gaussian. Overall, these findings indicate that the uncertainties of the MAIAC AOD retrievals over Cluj-Napoca are fairly well estimated.

4. Discussion

In this study, we used ground-based AERONET measurements taken between 2010 and 2020 to identify the dominant aerosol types, trends, and interactions with meteorological factors. First, we characterized the meteorological parameters considering the monthly values of temperature, relative humidity and wind speed and direction. Secondly, we analyzed the optical parameters such as the AOD at 500 nm and the Ångström exponent in order to determine the climatology of aerosols in the Cluj-Napoca region. The maximum average value of AOD was reached in July and August and the minimum in December and January. This is in agreement with AERONET average values of AOD at stations from South East Europe [74], [75]. The Ångström exponent showed minimum values in April and May and maximum values in August. The minimum is influenced by relatively frequent intrusions of mineral dust, as shown in other studies [20], [40]. The maximum values from August are mainly due to the presence of smoke from wildfires as other studies suggest [74].

Regarding the proportions of the primary aerosol classes, it is noteworthy to consider the similarities and differences with previous studies that utilized a similar methodology for classification. Similar to these study findings, the statistical analysis of 10 years of AERONET data from the station in Athens [45], showed that the most prevalent aerosol types were continental (19%), mixed (23%), and polluted aerosols (27%). Higher proportions of marine (11%) and dust aerosols (16%) were measured in Athens, a fact determined by the geographical positioning in a coastal area much closer to the African continent. A slightly smaller fraction of biomass burning aerosols was identified over Athens (5%). Stefan et al. in a similar study of the aerosol optical characteristics over the Romanian Black Sea Coast using AERONET data [75] showed the predominant presence of a fine fraction of aerosols of anthropogenic origin (mixed, polluted, biomass). In a study examining the typology of long-range transported aerosols over Europe from 2008 to 2018, [74] reported that smoke, continental, continental polluted, dust, and marine aerosols were present in Southeast Europe, with smoke being the most prevalent (43%), followed by continental (28%), continental polluted (12%), marine aerosols (8%), and dust (4%). [76] classified the aerosol types at the AERONET sites in the eastern Mediterranean and Black Sea using a method based on the sensitivity of microphysical aerosol properties and identified a higher proportion of marine aerosols, specific to coastal regions. Generally, in Eastern Europe studies focused on the proportion of primary aerosol classes showed similar results, with variations that can be attributed to differences in geographic locations, meteorological conditions, and analytical techniques.

In low AOD conditions, such as the one presented in this study, the MAIAC retrieval algorithm may exhibit a slight underestimation over an urban area (possible overestimation of the surface reflectance) dominated by absorbing aerosols (possible overestimation of the SSA) [32], [34], [57]. Figure 1 shows the land cover types present in the collocation area. The vast majority of land use is classified as urban, hence predominantly bright surfaces. A negative bias in such conditions has been observed by [34], [77], with similar RMSE values. [32] reported similar underestimations in the western part of China for low to moderate AOD. [38] showed that negative bias for 0.1 <AOD> 0.3 is common regardless of location while lower correlations may result from retrievals over brighter surfaces. [57] showed slight underestimations in urban and mixed areas over South America for low AOD. [23] reported better overall correlations and lower RMSE in urban areas of North America with similar slight negative bias. Complex urban areas are known to impose challenges for satellite aerosol retrievals where even slight errors in surface reflectance estimations may propagate and enhance overall bias [32], [57]. In cases of low AOD, the static aerosol models used by the MAIAC retrieval algorithm may not account for on-site conditions with respect to spatiotemporal variations of environmental attributes and aerosol properties [21].

Similar findings regarding the distribution of uncertainties as reported in our study are consistent with [34] who reported 70 % within EE for low AOD over bright surfaces. [38] reported between 62% and 75 % within EE for low AOD conditions while [39] reported 74% within EE over urban land. [57] reported between 45% and 57% within EE for urban areas and 64% to 68% for mixed areas although these values were indicative of a more stringent EE envelope of ± (0.05 + 0.05 × AOD). Neither of the authors reported a significant percentage below EE for a low AOD urban site. As such these values may be indicative of complex absorbing aerosols mixtures as seen in Figure 8 and the seasonal variation of SSA observed in Figure 7.

5. Conclusions

In this study, we analyzed a decade of AERONET measurements for the city of Cluj-Napoca, Romania, from 2010 to 2020. Version 3 Level 2.0 data have been used, except for SSA where Level 1.5 data were used. Maximum AOD values are observed in July and August (0.23 and 0.21 ± 0.11, respectively), while the minimum values are recorded in December and January (0.15 ± 0.09). The overall mean value of the AE for the entire measurement period is 1.5 ± 0.29, with a maximum of 1.62 ± 0.26 in August, and a minimum of 1.33 ± 0.39 in April.

Based on the classification proposed by Dubovik et al. [44], the dominant types of aerosols in Cluj-Napoca area are represented as mixed (32%), polluted (29.4%), and continental aerosols (19.1%).

The MAIAC algorithm was fairly well correlated with AERONET over the urban environment in Cluj-Napoca, Romania. The slightly negative bias in low AOD conditions seems to be a common feature reported in scientific literature and may be indicative of complex absorbing aerosol mixtures. The uncertainties roughly followed a Gaussian distribution with small overestimations on the low end and some underestimated on the high end. Since the MAIAC retrieval algorithm utilizes static aerosol models, the spatiotemporal variations of aerosol properties in urban environments may not be sufficiently represented, thus inducing additional errors.

This study comes to complete and improve other similar studies from other regions. Further analysis of different sources of data could more clearly improve the signatures of aerosol types in the region.