Analysis of Long-Term Aerosol Optical Properties Combining AERONET Sunphotometer and Satellite-Based Observations in Hong Kong

: This study analyzes seasonal characteristics and long-term variations in aerosol optical parameters in Hong Kong from 2006 to 2021 using AERONET data and satellite-based observations based on the extreme-point symmetric mode decomposition (ESMD) model. The dominant aerosol types in Hong Kong are mixed aerosols and urban/industrial aerosols with ﬁne-mode sizes, and slightly absorbing or non-absorbing properties. Aerosol optical depth (AOD), Angstrom exponent (AE) and single scattering albedo (SSA) varied seasonally with a lower AOD but higher AE and SSA in summer, and elevated AOD but lower AE and SSA in spring and winter. The long-term variations show the year 2012 to be a turning point, with an upward trend in AOD and AE before 2012 and then downwards after 2012. However, for SSA, a rising trend was exhibited in both pre- and post-2012 periods, but with a larger gradient in the ﬁrst period. The ESMD analysis shows shorter-term, non-linear ﬂuctuations in aerosol optical parameters, with alternating increasing and declining trends. The examination of the relationships between AOD and meteorological factors based on the extreme gradient boosting (XGBoost) method shows that the effects of weather conditions on AOD are complex and non-monotonic. A lower relative humidity, higher wind speed in southwest directions and lower temperature are beneﬁcial to the abatement of aerosol loads in Hong Kong. In conclusion, the ﬁndings of this study enhance the understanding of aerosol properties and the interactions between aerosol loading and meteorological factors.


Introduction
Aerosols are liquid or solid particles of different sizes suspended in the atmosphere, and are emitted from anthropogenic and natural sources. Extensive studies have been conducted to investigate the aerosol distribution and analyze the related effects on climate, the atmospheric environment and human health [1][2][3][4][5][6]. They can influence the Earth-atmosphere system and climate both indirectly (changing lifetime, amount and radiative properties of cloud by modifying cloud physical properties) and directly (altering the radiative forcing via scattering and absorbing radiation) [7]. Moreover, aerosol leads to considerable uncertainty in Earth-atmosphere interactions and climate change [8]. Thus, quantifying the long-term trends of aerosol parameters and the underlying causes are beneficial to understanding regional air quality conditions and climate variability.
The optical parameters of aerosols, aerosol optical depth (AOD), Angstrom exponent (AE) and single scattering albedo (SSA) indicate some important characteristics of aerosols, Therefore, to fill the aforementioned gaps, and to analyze the impacts of meteorological conditions on AOD, the aerosol parameters measured by the PolyU AERONET station in Hong Kong and the satellite-based observations were integrated and investigated based on ESMD and a data-driven method (XGBoost). This study was conducted to (1) identify the seasonal characteristics in aerosol optical parameters and the dominant aerosol type in Hong Kong; (2) investigate the long-term variations in AOD, AE and SSA from annual and seasonal perspectives between 2006 and 2021; and (3) reveal how AOD responds to the changes in weather conditions.

AERONET Data and Sampling Site
AERONET (https://aeronet.gsfc.nasa.gov/ accessed on 18 October 2022) is a global aerosol monitoring network that provides accurate measurements of aerosol optical parameters at multiple wavelengths [13]. There are three quality levels: Levels 1.0, 1.5 and 2.0, which represent raw observations, cloud-screened data with quality controls but not final calibration, and quality-assured data after pre-and post-field calibrations, respectively. Two AERONET stations are deployed in Hong Kong, namely Sheung and PolyU stations. The Sheung station is a suburban site, while the PolyU station shows subtropical urban aerosol characteristics (22.3033°N, 114.1797°E, 30 m) with a longer service period (starting from November 2005), located on the roof of the Hong Kong Polytechnic University ( Figure 1). In this study, daily and monthly Version 3 Level 2.0 AE at 440-675 nm, AOD at 500 nm and Level 1.5 SSA at 440 nm from 2006 to 2021 were adopted. The uncertainty in the obtained AOD is up to 0.02, while SSA is an inversion product with an uncertainty of about 0.03-0.05 [25].

Satellite-Based Observations
Due to the maintenance and calibration of sunphotometers, or weather conditions, missing values of AERONET measurements are inevitable. Therefore, satellite-based observations were considered to ensure continuous time-series sequences. MCD19A2 (https://lpdaac.usgs.gov/ accessed on 18 October 2022), an AOD product at 470 nm and 550 nm wavelengths provided by MODIS with daily intervals and a 1 km spatial resolution via the multi-angle inversion atmospheric correction (MAIAC) approach, was adopted, which is highly consistent with AERONET ground observations [26,27]. A comparison between AERONET AOD measurements and the MODIS MAIAC AOD dataset in different seasons in Hong Kong is displayed in Figure 2. Daily mean satellite-based AOD at 470 nm and 550 nm was extracted by defining with the PolyU AERONET station in the center of an area of 5×5 pixels to interpolate the AE at 440-675 nm and AOD at 500 nm according to the equations presented in this section [28]. The results show that the

Satellite-Based Observations
Due to the maintenance and calibration of sunphotometers, or weather conditions, missing values of AERONET measurements are inevitable. Therefore, satellite-based observations were considered to ensure continuous time-series sequences. MCD19A2 (https://lpdaac.usgs.gov/ accessed on 22 August 2022), an AOD product at 470 nm and 550 nm wavelengths provided by MODIS with daily intervals and a 1 km spatial resolution via the multi-angle inversion atmospheric correction (MAIAC) approach, was adopted, which is highly consistent with AERONET ground observations [26,27]. A comparison between AERONET AOD measurements and the MODIS MAIAC AOD dataset in different seasons in Hong Kong is displayed in Figure 2. Daily mean satellite-based AOD at 470 nm and 550 nm was extracted by defining with the PolyU AERONET station in the center of an area of 5 × 5 pixels to interpolate the AE at 440-675 nm and AOD at 500 nm according to the equations presented in this section [28]. The results show that the satellite-based observations are in good agreement with the ground-based measurements, with a correlation coefficient from 0.72 to 0.77.
where τ(λ i ) represents AOD at λ i wavelength, α denotes the Angstrom exponent and β is the turbidity coefficient. Subsequently, monthly mean AOD and AE were calculated based on valid daily measurements to complement the missing AERONET data. The percentage of the satellite-based observations used in this study to substitute the AERONET missing data accounts for about 38%. satellite-based observations are in good agreement with the ground-based measurements, with a correlation coefficient from 0.72 to 0.77.
where ( ) represents AOD at wavelength, α denotes the Angstrom exponent and β is the turbidity coefficient. Subsequently, monthly mean AOD and AE were calculated based on valid daily measurements to complement the missing AERONET data. The percentage of the satellite-based observations used in this study to substitute the AERONET missing data accounts for about 38%. For the missing values of SSA, the Multi-angle Imaging SpectroRadiometer (MISR) measurements were used, which can provide SSA inversion results with a spatial resolution of 4.4 km and a repeat coverage of about seven to nine days for the sampling site (https://misr.jpl.nasa.gov/ accessed on 18 October 2022). Daily mean SSA at 446 nm was extracted by defining the PolyU AERONET station in the center of an area of 3 × 3 pixels. Some studies have stated that the retrieval SSA of MISR was overestimated as compared to the AERONET, indicating that the effects of some highly absorbing particles were ignored [29][30][31]. Thus, the retrieved SSA of MISR was corrected by using the SSA inversion For the missing values of SSA, the Multi-angle Imaging SpectroRadiometer (MISR) measurements were used, which can provide SSA inversion results with a spatial resolution of 4.4 km and a repeat coverage of about seven to nine days for the sampling site (https:// misr.jpl.nasa.gov/ accessed on 22 August 2022). Daily mean SSA at 446 nm was extracted by defining the PolyU AERONET station in the center of an area of 3 × 3 pixels. Some studies have stated that the retrieval SSA of MISR was overestimated as compared to the AERONET, indicating that the effects of some highly absorbing particles were ignored [29][30][31]. Thus, the retrieved SSA of MISR was corrected by using the SSA inversion results of AERONET using the linear regression. After this, the monthly average SSA was calculated to replace the missing values of AERONET.

Meteorological Data
In this study, the measurements of the King's Park weather station (22. Figure 1b). The daily ground-level mean measurements can be acquired from the Hong Kong Observatory, including wind direction, wind speed, relative humidity, temperature, pressure and precipitation (https://www.hko.gov.hk/en/index.html accessed on 22 August 2022). Precipitation was not involved since the ground-level measurements of aerosol optical parameters are not available on rainy days. Additionally, to ensure accuracy and reduce uncertainties, the weather station records consistent with valid measurements of the PolyU AERONET site were extracted and used to explore the impacts of meteorological variables on AOD.

Extreme-Point Symmetric Mode Decomposition (ESMD)
ESMD was developed based on ensemble empirical mode decomposition (EEMD) [32], which can be used to detect the time-dependent non-linear tendency without antecedent assumption. This is a newly popular approach to probing non-linear varying trends of climate-related issues [33,34]. The original signals can be decomposed into several intrinsic mode functions (IMFs) and a residual curve (R).
where M i (m) is the i th IMF, R(m) is the tendency curve and n denotes the total number of IMFs. More detailed descriptions can be found in the Supplementary Materials. In our study, R denotes the variation tendency of aerosol optical parameters from 2006 to 2021.

Extreme Gradient Boosting (XGBoost)
XGBoost is an ensemble decision tree model of machine learning with high accuracy for classification and prediction issues. It has been widely applied to explore the relationships between dependent and independent variables. More detailed descriptions can be found in the literature [35]. In this study, XGBoost was adopted to identify the relative importance of meteorological factors and investigate their interactions with AOD, which can be generally demonstrated as follows: where y AOD means the daily averaged measurements of the PolyU AERONET station, and x i represents meteorological factors, including pressure, temperature, wind direction, wind speed and relative humidity. The relationships between AOD and meteorological factors were estimated by the partial dependence plot results based on the XGBoost model, which can visualize how AOD responds to changes in a specific variable by marginalizing other impact factors [36]. In other words, the variable sensitivity can be illustrated by how the dependent variable reacts when a given predictor changes while other predictors remain constant. To avoid overfitting and to examine reliability, the 10-fold cross-validation (CV) was adopted, whereby the dataset was divided into ten equal-sized subgroups. One subgroup was used as the test dataset and the other samples are for training. The training process is repeated ten times to ensure that every subgroup is tested once. Mean absolute error (MAE), root-mean-square error (RMSE) and the coefficient of determination R-squared (R 2 ) were used to assess the performance of XGBoost. in summer [37]. A much higher AOD was observed in spring, with about 58% of AOD levels above 0.5. Summer, autumn and winter have lower AOD values, especially in summer, when almost 50% of AOD is below 0.3. Thus, spring has the largest mean AOD (0.58), followed by autumn and winter. This may be attributed to long-distance dust and biomass-burning sources, as well as rapid production mechanisms by secondary aerosols [38]. Summer has the lowest mean AOD, as the prevailing wind is southerly in summer, bringing clean air from the South China Sea [39]. Additionally, there is a large amount of rainfall in summer, which can scavenge the pollutants and cleanse the atmosphere [7].

Aerosol Types and Absorbing Aerosol Types
Two different classification criteria suggested by Salinas et al. [41] and Zheng et al. [42] were used to define the aerosol types and absorbing aerosol types in Hong Kong during the study period (see Table 1). For the first one, a daily averaged AOD at 500 nm and AE at 440-675 nm were used to classify the aerosol types, as suggested by Salinas et al. [41]. As depicted in Figure 4a, mixed aerosol (MIXA1 and MIXA2) is the representative type in Hong Kong, accounting for nearly 50%. Instead of pure aerosol, it is composed of different types of aerosol (i.e., dust, urban and biomass burning aerosols). The second dominant type is urban/industrial aerosol (UIA), with a percentage of 33%, followed by biomass-burning aerosol (BBA). Marine aerosol (MA) and dust aerosol (DA) only make up a small fraction (6%). For the second classification criteria, to further explore the dominant aerosol size and properties of absorbing and scattering aerosols, a daily mean Level 2.0 AE at 440-675 nm, fine-mode fraction (FMF) at 500nm and Level 1.5 SSA at 440 nm were used as the aerosol optical parameters, as proposed by Zheng et al. [42]. FMF represents the proportion of fine-mode aerosols in the atmosphere as compared to total aerosols [43,44]. The results are displayed in Figure 4b. Consistent with the above statements, anthropogenic fine-mode aerosols with non-absorbing (FNA) or slight absorbing (FSA) characteristics are the prevailing type in Hong Kong, accounting for over 66% overall. Meanwhile, there are some fine-mode aerosols exhibiting medium (FMA) and highly absorbing (FHA), with percentages of 16.3% and 4.9%, respectively. Marine absorbing (M-A) and non-absorbing aerosols (M-NA) have comparable proportions, 6.5% and 4.8%. The percentage of coarse-mode aerosols is only small (0.7%), and most of them show absorption properties. Analogously, AE also exhibited a unimodal structure centered at 1.3-1.4 in spring, summer and winter, while centered at 1.4-1.5 in autumn, suggesting the dominance of fine-mode aerosol ( Figure 3b). Significant seasonal disparities of AE can be found in summer, with a wide range from 0.17 to 2.47 and a high SD of 0.32. This indicates variations in aerosol particle types in the atmosphere (i.e., fine-mode pollutants and coarse-mode dust). However, Figure 3b suggests that anthropogenic fine-mode particles make more contributions than coarse-mode pollutants over the study period in Hong Kong. Figure 3c indicates that SSA has a similar distribution frequency histogram in all seasons except summer. The wide range of SSA, 0.721 to 0.996, indicates different aerosolabsorbing characteristics from strongly scattering to strongly absorbing aerosols. A large percentage of high SSA above 0.9 can be observed in all seasons, suggesting that strongly scattering aerosols dominated in Hong Kong during the study period. The highest mean SSA was found in summer, indicating strong scattering properties. This is partly ascribed to the aerosol particles tending to capture atmospheric water vapor under high relative humidity in summer, making a larger particle radius and further increasing the aerosol scattering [40]. Due to the effects of the winter monsoon, the transport of pollutants from China's mainland was facilitated by the prevailing northerly wind [39], and hence absorbing particles (i.e., black carbon) are increased in winter in Hong Kong.

Aerosol Types and Absorbing Aerosol Types
Two different classification criteria suggested by Salinas et al. [41] and Zheng et al. [42] were used to define the aerosol types and absorbing aerosol types in Hong Kong during the study period (see Table 1). For the first one, a daily averaged AOD at 500 nm and AE at 440-675 nm were used to classify the aerosol types, as suggested by Salinas et al. [41]. As depicted in Figure 4a, mixed aerosol (MIXA1 and MIXA2) is the representative type in Hong Kong, accounting for nearly 50%. Instead of pure aerosol, it is composed of different types of aerosol (i.e., dust, urban and biomass burning aerosols). The second dominant type is urban/industrial aerosol (UIA), with a percentage of 33%, followed by biomass-burning aerosol (BBA). Marine aerosol (MA) and dust aerosol (DA) only make up a small fraction (6%). For the second classification criteria, to further explore the dominant aerosol size and properties of absorbing and scattering aerosols, a daily mean Level 2.0 AE at 440-675 nm, fine-mode fraction (FMF) at 500nm and Level 1.5 SSA at 440 nm were used as the aerosol optical parameters, as proposed by Zheng et al. [42]. FMF represents the proportion of fine-mode aerosols in the atmosphere as compared to total aerosols [43,44]. The results are displayed in Figure 4b. Consistent with the above statements, anthropogenic fine-mode aerosols with non-absorbing (FNA) or slight absorbing (FSA) characteristics are the prevailing type in Hong Kong, accounting for over 66% overall. Meanwhile, there are some fine-mode aerosols exhibiting medium (FMA) and highly absorbing (FHA), with percentages of 16.3% and 4.9%, respectively. Marine absorbing (M-A) and non-absorbing aerosols (M-NA) have comparable proportions, 6.5% and 4.8%. The percentage of coarsemode aerosols is only small (0.7%), and most of them show absorption properties. Table 1. Classification criteria based on aerosol optical parameters. Additionally, seasonal aerosol types and absorbing aerosol types were also investigated to identify the discrepancies between aerosol types in different seasons. As shown in Figure  S1, MIXA and UIA are the predominant types in all seasons, and the percentage of UIA in winter is the highest, reaching up to 42.6%. This may be because of the long-range transport of increased urban aerosols during the cold winter heating season in the Chinese mainland [37]. For BBA, it has a higher proportion in spring, possibly caused by frequent forest fires and the increasing burden of biomass burning in spring. The percentages of MA and DA remain relatively constant in different seasons, except in spring, during which MA particles were not found. The seasonal absorbing aerosol types are defined and displayed in Figure S2. Discernible differences occur between summer and winter, with larger quantities of FNA in summer and decreasing frequencies in winter, implying strong scattering properties in summer but more absorption in winter. This is consistent with the statistics of SSA given in Figure 3c Figure S3); additionally, the annual mean AOD decreased gradually and a significant decline trend can be seen for 2012 to 2021, which is consistent with AOD trends in the PRD region suggested by previous studies [45]. Thus, the study period was divided into two subperiods, namely pre-2012 and post-2012. The linear regression shows an increasing tendency in AOD and AE during the pre-2012 period with a significant decreasing trend after 2012 (Figure 5a,b). Contrary to the decrease in AOD and AE, a slight increase in SSA was found, at a rate of 0.0001 during the post-2012 period (Figure 5c).

Types Aerosol Optical Parameters Thresholds
A time series of aerosol optical parameters were decomposed into five IMFs and one variation trend curve based on the ESMD approach. The variance contribution rate (VCR) and the significance test at 95% confidence level of each IMF and the residual curve were calculated and are displayed in Table 2. It can be seen that the IMFs and trend curves are above the 95% confidence level and significantly correlated with the original time series, except the IMF5 of AOD. As suggested by Zheng et al. [21], the IMFs contain periodic oscillations and irregular non-periodic changes. Take the decomposed results of AOD time series as an example: for the irregular fluctuations (i.e., IMF1 and IMF2 in Figure S4a), the unstable amplitudes of IMFs may be caused by various unexpected factors, including weather conditions, short-term pollution control policies and human activity intensity. Periodic oscillations generally include seasonal (IMF3 in Figure S4a) and inter-annual periodic changes (IMF4 and IMF5 in Figure S4a) related to geographical phenomena. The residual curve can be used to denote the variation trend, which is the main focus of our study. The varying trends analyzed by ESMD are not completely consistent with the linear trends, presenting more fluctuating and detailed variations in aerosol optical parameters. Using AOD as an illustration, the ESMD-based AOD trend (blue dashed line in Figure 4a) shows that there is a decline from 2006 to 2008 with the value changing from above 0.5 to around 0.4, and the values increased subsequently in 2009. The descending trend seen up to 2008 may be caused by emissions control [46] and the reduction in industrial air pollution during the global economic crisis, and an increasing trend can be seen after 2009 since the emissions from economic production and transport returned to pre-recession levels, as suggested by Buchholz et al. [47]. After this, a remarkable descending trend can be found, with the value decreasing to around 0.3 in 2021, which has been attributed to the encouragement of renewable sources of energy and relative air pollution control policies [48,49]. For AE, an increasing tendency similar to the linear trend can be observed before 2012, whereas noticeable fluctuating variations (decreasing-increasing-decreasing) are present in the second period. Especially, there is a steep declining trend after 2020, possibly caused by the reduction in human activities due to the effects of COVID-19 lockdown policies [50]. Nonetheless, most of the monthly averaged AE values are larger than 1.0, suggesting the dominance of anthropogenic aerosols. Figure 5c shows more oscillations in SSA, alternating between upward and downward tendencies. SSA increased substantially after 2017, implying a reduction in absorbing aerosols (long-range transport biomass burning and local vehicles) and more aerosols with strong scattering properties.
The underlying cause appears to be the formulation of air pollution control policies and the promotion of new energy sources. well as the proportion of absorbing aerosols in Hong Kong, were reduced during the study period, suggesting that the aforementioned air quality control countermeasures were effective.

Seasonal Variations of Aerosol Optical Parameters
To reveal a detailed pattern of seasonal variations in aerosol optical parameters, linear and ESMD-based non-linear trends were derived using monthly averaged measurements directly. Figure 6a shows that the overall linear seasonal trends in AOD are con-  In conclusion, the linear trend is useful for summarizing the overall variable trends during the study period, whereas the ESMD-based non-linear trends can accurately capture more detailed fluctuations in data, and provide more detailed and plausible variations of aerosol optical properties. The overall results indicate that AOD and AE levels, as well as the proportion of absorbing aerosols in Hong Kong, were reduced during the study period, suggesting that the aforementioned air quality control countermeasures were effective.

Seasonal Variations of Aerosol Optical Parameters
To reveal a detailed pattern of seasonal variations in aerosol optical parameters, linear and ESMD-based non-linear trends were derived using monthly averaged measurements directly. Figure 6a shows that the overall linear seasonal trends in AOD are consistent with the long-term variation directions, but the magnitudes of trends vary in different seasons. For AE, the linear seasonal trends are similar to the long-term variations, except in winter. There is an increase from 2013 to 2021 in winter, attenuating the substantial decline presented in other three seasons and leading to a moderate decrease in the annual trend in the post-2012 period. Similarly, descending variations in linear seasonal trends of SSA are observed in winter, which is different from the increasing pattern in other seasons.

Quantitative Impacts of Meteorological Factors
Meteorological conditions are vital influencing factors for AOD, which have nonnegligible impacts on the composition and transmission path of aerosols [9]. Pressure (hpa), wind speed (km/h), wind direction (degrees), relative humidity (%), and temperature (°C) are used to quantify the relative contributions of each factor and complex rela- A time series of seasonal aerosol optical parameters was also decomposed into several IMFs and the trend curve was based on the ESMD approach. The variance contribution rate (VCR) and the significance test at a 95% confidence level of the varying trends were calculated and are displayed in Table 3, and the estimation results of each IMF are shown in Table S1. The results show that the varying trends are above the 95% confidence level and have a stronger correlation with the original time series, except for the varying trend of summer AOD. For the ESMD-based non-linear trends, these are not entirely aligned with the long-term tendencies shown by the linear analysis. Specifically, in the pre-2012 period, a decreasing-increasing trend of AOD was exhibited in summer and autumn, while a steady upward trend was seen in spring and winter. In terms of the post-2012 period, significant declining trends dominated in all seasons. Strikingly, during the COVID-19 lockdown period (2020-2021), AOD variations leveled off in spring and summer but pronounced decreasing trends occurred in autumn and winter. Apart from the effects of restricted human activities during COVID-19, the reduction in background pollution in Guangdong province is beneficial as well, the upwind area of Hong Kong under winter monsoon [51,52]. Fan et al. [53] also suggested that a considerable mitigation of air pollution can be found across the Pearl River Delta (PRD) area due to the COVID-19 lockdown. As depicted in Figure 6b, spring and autumn have similar changing patterns, with fluctuating variations (increasing-decreasing-increasing) in the pre-2012 period and a continuous downward trend from 2013 to 2021. Although a slightly rising trend in winter AE was presented after 2017, it was offset by the rapid decline in summer AE and a sustained descent in autumn and spring, leading to an overall continuous downward trend in AE after 2017. With regard to the seasonal variations of SSA, the downward trend in spring SSA occurred more frequently, presumably caused by the considerable biomass-burning activities across Indochina Peninsula in spring [54,55]. In addition, during the pre-2012 period, the noticeable increase in winter compensated for the decreasing trend in other seasons, making a similar varying trend to the long-term changes. Although there is a slight descending linear trend in winter SSA during the post-2012 period, a significant rising trend dominated after 2017, indicating a reduction in absorbing aerosols (i.e., black carbon) in recent years. The striking increase in SSA can be found in all seasons from 2020 to 2021, which may be ascribed to the reduction in absorbing aerosols emanating from long-range biomass burning, as well as the restrictions on vehicle transport and human activities during COVID-19 lockdown period [56,57]. Overall, both linear and ESMD-based non-linear seasonal variations of aerosol optical parameters are not completely consistent with the long-term trends. This is because a higher AOD can be observed in spring than in other seasons and spring AOD has been exhibiting a stable trend in recent years. The results imply that attempting to mitigate spring AOD is beneficial to achieving a remarkable drop in overall AOD in Hong Kong. Table 3. Variance contribution rate (VCR) and correlation coefficient of the residual curve.

Quantitative Impacts of Meteorological Factors
Meteorological conditions are vital influencing factors for AOD, which have nonnegligible impacts on the composition and transmission path of aerosols [9]. Pressure (hpa), wind speed (km/h), wind direction (degrees), relative humidity (%), and temperature ( • C) are used to quantify the relative contributions of each factor and complex relationships between explanatory variables and AOD based on the Xgboost model. In our study, after the 10-fold CV test, the XGBoost method achieves an R 2 of 0.80, MAE of 0.08 and RMSE of 0.11. The estimated relative importance of meteorological factors is displayed in Table 4, which is calculated by evaluating the increase in the model prediction error when randomly shuffling a single feature, with the range from 0 to 1 [58]. The results suggest that temperature exerts the greatest influence on AOD, while relative humidity is the second primary influencing factor, followed by wind direction, with an importance rate of about 0.2. Wind speed and pressure make lower contributions to AOD, accounting for 0.18 and 0.16, respectively. Overall, the meteorological variables used in this study have comparable influences in estimating AOD, with a relative importance of around 0.2. Instead of Pearson correlation analysis and multi-linear regression, partial dependence derived by XGBoost was employed to investigate the non-monotonic relationship between a given variable and AOD (Figure 7). The results indicate that AOD increases significantly with relative humidity, increasing from 50% to 70% due to the hygroscopic growth of aerosols [59], whereas there is a noticeable decline in AOD when the relative humidity exceeds 70%. This may be because a higher relative humidity tends to induce frequent precipitation events and the aerosol particles may be scoured, resulting in lower AOD levels [22]. Notably, AOD remains almost unchanged when relative humidity is lower than 50%, which may be responsible for the sparse records of lower relative humidity (see Figure S5), making the correlation insignificant [60]. In contrast, a strong negative correlation can be found between wind speed and AOD, as a higher wind speed promotes aerosol diffusion, dilution and circulation, leading to better air quality [61]. Nonetheless, there are some slight positive trends between wind speed and AOD. This can be attributed to the surface soil and dust raised by high wind speeds, thus increasing the aerosol particles [62]. As depicted in Figure 7c, there is a strong relationship between wind direction and AOD. Specifically, significant upward trends in AOD occur when wind directions change from north to northeast (0-45 degrees) and southwest to west (225-270 degrees). Higher AOD levels can be observed in the northeast (45 degrees) and west (270 degrees) directions. This is more likely caused by the contributions of the long-distance transport of industrial emissions from the mainland (i.e., Henan, Shandong and Anhui provinces) [9], as well as dust and agricultural pollution from Southeast Asia [53]. With the wind directions changing from northeast to east (45-90 degrees) and from west to northwest (270-315 degrees), the neighboring regions (i.e., Guangdong, Fujian and Guangxi provinces) of Hong Kong have cleaner background pollution, leading to a decrease in AOD. The southern part of Hong Kong is surrounded by sea (see Figure 1), which is most likely the reason why AOD remains stable at a wind direction of 90-225 degrees. As can be seen in Figure 7d, AOD is sensitive to variations in temperature, and AOD increases with increased temperature when the temperature is below 25 • C. However, when the temperature is in the range of 28-29 • C, a negative association with AOD is presented. According to the partial dependence results of pressure, AOD is less susceptible to changes in pressure (Figure 7e). AOD decreases gradually with an increasing pressure from 1005 to 1010 hpa, whereas a positive correlation can be observed when the pressure is in the range of 1010-1015 hpa. Subsequently, AOD stabilizes between 0.45 and 0.5 when the pressure exceeds 1015 hpa. In summary, unlike for typical linear correlations, partial dependence analysis implemented by the XGBoost can derive complex and non-linear relationships between AOD and meteorological factors. How AOD responds to the changes in driving factors can be estimated and visualized, thus recommending the use of black-box models for facilitating the understanding of interactions between AOD and meteorological conditions. to the partial dependence results of pressure, AOD is less susceptible to changes in pressure (Figure 7e). AOD decreases gradually with an increasing pressure from 1005 to 1010 hpa, whereas a positive correlation can be observed when the pressure is in the range of 1010-1015 hpa. Subsequently, AOD stabilizes between 0.45 and 0.5 when the pressure exceeds 1015 hpa. In summary, unlike for typical linear correlations, partial dependence analysis implemented by the XGBoost can derive complex and non-linear relationships between AOD and meteorological factors. How AOD responds to the changes in driving factors can be estimated and visualized, thus recommending the use of black-box models for facilitating the understanding of interactions between AOD and meteorological conditions.

Discussion
Using ground-based AERONET measurements taken over 16 years and combined with satellite-based observations, our study explores the dominant aerosol types, long-term varying trends and the interactions with meteorological factors. The main contributions of our study can be summarized as follows: First, in addition to classifying the aerosol types based on AOD and AE values in previous studies [63][64][65][66], the aerosol-absorbing properties in Hong Kong were also categorized using AE, SSA and FMF. Moreover, the multi-seasonal AOD, AE and SSA distribution patterns in Hong Kong were discussed, which provides a basis for understanding the seasonal aerosol optical characteristics in subtropical urban areas such as Hong Kong. Notably, the thresholds of the aerosol optical parameters are generally used for dominant aerosol type classification, while the proportions of specific aerosol types depend on the selected aerosol parameters and the thresholds used for the classification [63]. Additionally, using aerosol optical parameters for aerosol type identification is suitable for areas with a stable aerosol composition, but may not be well-suited for locations with complex topography, Arctic regions and remote marine areas, as aerosols in these regions are difficult to represent using median optical parameters [67]. Second, variations in AOD were mainly analyzed in most previous studies in the literature [16][17][18][19], whereas changes in AE and SSA are rarely included, leading to an incomplete view of aerosol variations. Apart from the linear trend, the non-linear fluctuations of aerosol optical parameters were also discussed based on the ESMD method, which can reflect more detailed varying trends. Finally, different from previous studies [21,22], non-monotonic interactions between AOD and meteorological factors were analyzed using a data-driven model, XGBoost, which contributes to the knowledge of how AOD responds to the changes in weather conditions. However, there are still some limitations in our study. First, the minimum threshold for the number of AERONET AOD observations per day was not considered. To enrich the sample size of the AOD observations, the daily mean AOD was not filtered, although there was only one observation per day. In the future, multi-source AOD datasets (i.e., reanalysis data and satellite-based observations) or more accurately retrieved AOD should be included for the substitute of AERONET AOD to generate accurate and continuous AOD time series. Second, due to the availability of meteorological factors, the interactions between AOD and other meteorological factors (i.e., planetary boundary layer, PBL) were ignored, which has been suggested to have a significant impact on aerosol vertical distributions and the dispersion of air pollutants in previous studies [68,69].

Conclusions
In this study, combining aerosol optical parameters measured by the PolyU AERONET station with satellite observations (i.e., MODIS and MISR aerosol products) from 2006 to 2021 was adopted to analyze the seasonal distribution patterns, varying trends and the association with meteorological factors in Hong Kong. Four primary conclusions can be summarized from the results. Firstly, aerosol optical parameters present significant seasonal characteristics. The seasonal mean AOD is at its highest in spring, while summer has the lowest mean AOD but the largest averaged AE and SSA. Secondly, mixed aerosol and urban/industrial aerosol, which have a fine-mode size with slight or non-absorbing properties, are predominant types in Hong Kong during the study period. Thirdly, the year 2012 can be deemed as the turning point. Both AOD and AE exhibit an increasing trend in the pre-2012 period and a downward tendency after 2012, whereas a consistent rising trend of SSA can be observed over both periods. The use of ESMD to investigate the non-linear annual and seasonal varying trends, was able to provide more detailed insights than linear trends. AOD, AE and SSA present "decreasing The results suggest that the sharp AOD decrease in autumn and winter in recent years contributed to the decline in overall AOD in Hong Kong. Finally, instead of a linear correlation, XGBoost was employed to quantify the non-monotonic relationships between meteorological factors and AOD, as well as rank the relative importance of each variable. It can be concluded that meteorological factors have comparable contributions to AOD and that interactions are complex. The results imply that meteorological conditions with lower relative humidity, higher wind speed, winds from east and southwest directions, lower temperature and higher pressure tend to lower AOD levels. In conclusion, our study can enhance the understanding of aerosol characteristics and long-term variations in Hong Kong. Furthermore, the interactions between meteorological factors and AOD can serve as a basis for developing air quality promotion strategies.