1. Introduction
The accelerated urban expansion and intensifying urbanization have imposed substantial environmental pressures, especially air quality, leading to a marked increase in PM2.5 concentrations (Zhou et al., 2023). These fine matter particulates, predominantly derived from transportation, industrial activities, and biomass combustion, are also transported across regions. PM2.5 denotes particles with an aerodynamic diameter of less than 2.5 μm, which are harmful to human health and ecosystems (Wang et al., 2021; Gao et al., 2017; Chen & Chen, 2021). Exposure to fine particulate matter is associated with reduced lung function and an increase in respiratory symptoms, including airway irritation, coughing, and breathing difficulties, and it may exacerbate COVID-19 symptoms (Carretero-Peña et al., 2019). According to a 2012 report by the United States Environmental Protection Agency (EPA), air pollution is responsible for 1.8% to 6.4% of child mortality (ages 0 to 4) across Europe, resulting in an estimated 100,000 annual deaths in urban areas (EPA, 2012). Additionally, in 2015, exposure to ambient PM2.5 was estimated to contribute to approximately 4.2 million premature deaths worldwide (Cohen et al., 2017).
The swift expansion of urban areas and the construction of towering skyscrapers, serving as both residential and commercial spaces, have resulted in human exposure to air pollution at various heights. Recently, numerous studies have focused on the vertical distribution of air pollution, with particular attention to PM2.5 levels in urban settings. These studies have employed diverse methodologies, including Lidar technology (Tao et al., 2016; Lyu et al., 2018; Liu et al., 2019; Wang et al., 2016), UAV-based measurements (Zhao et al., 2021; Qu et al., 2022; Wu et al., 2021), high-rise building observations (Yang et al., 2009; Liu et al., 2018; Ding et al., 2005), and computational modeling techniques (Liu et al., 2021; Lee et al., 2017; Yang et al., 2020).
Hanoi, the capital and most populous city of Vietnam had an estimated population of 8.6 million in 2023 (Statistical Publishing House, 2023). The city’s rapid socio-economic growth and urbanization have exacerbated severe air pollution, particularly concerning PM2.5 levels. From 2015 to 2021, the annual average PM2.5 concentration in Hanoi surpassed the QCVN 05:2013/BTNMT standard by over 1.3 times (MONRE, 2023), highlighting persistent pollution challenges. This has led to extensive research on PM2.5, including emission sources (Bang, 2022; Amann et al., 2019; Le et al., 2022), composition (Vo et al., 2022; Hai & Oanh, 2013; Makkonen et al., 2023), spatial distribution (Vuong et al., 2023), meteorological impacts (Luong et al., 2021; Ly et al., 2021; Hien et al., 2002; Tham et al., 2018; Ngo et al., 2012) and health consequences (Nguyen et al., 2022; Nhung et al., 2022). However, studies addressing the vertical distribution of PM2.5 in Vietnam remain relatively scarce.
Therefore, it is necessary to conduct a study evaluating the vertical distribution of air pollution in Hanoi. In this research, PM2.5 concentrations were measured at different heights across the city, adhering to EPA reference standards. The findings from the assessment of PM2.5 distribution will offer a more detailed understanding of its spatial patterns and the magnitude of its impact at various heights within the urban environment in Hanoi.
2. Materials and Methods
2.1. Methods
The primary methodology employed in this study involved measuring meteorological parameters and PM2.5 concentrations on the rooftops of high-rise buildings. Due to the absence of meteorological or television towers of approximately 200 meters in height in Hanoi, the research team opted for tall building rooftops to conduct simultaneous measurements. This approach enabled the study to effectively capture PM2.5 data at various elevations within urban areas.
In this study, sampling was conducted at two urban locations within the inner city of Hanoi, with each site having simultaneous measurements taken at five distinct points situated on the rooftops of buildings at varying elevations as shown in Figure 1. Vertical distribution data for the examined parameters were obtained by sampling at the highest floors of buildings, with elevations ranging from 40 meters to over 300 meters (Table 1). Wind speed, humidity, and temperature were directly measured at different heights using a microclimate measuring device. At the same time, PM2.5 concentrations were sampled using the Instrumex dust sampler and analyzed with an AD-421D-32 analytical balance.
Five Instrumex samplers, designed to concurrently collect PM2.5 and PM10 particles in compliance with U.S. EPA reference standards, were employed simultaneously for sampling. Each sample was collected over 6 hours, with two samples taken daily between 8 AM and 8 PM.
Figure 1. Collecting sample images. (a) KengNam 49; (b) CEO Tower.
Table 1. Parameters considered for correlation assessment.
Independent variables |
Elements |
X0 |
Height |
X1 |
Temperature |
X2 |
Wind speed |
X3 |
Humidity |
X4 |
PM10 Concertration |
2.2. Analysis of the Relationship between PM2.5 and Other Factors
There are multiple approaches for analyzing the relationship between PM2.5 concentrations and meteorological variables, with Multivariable Linear Regression (MLR) and Multivariable Nonlinear Regression (MNLR) being the most commonly employed techniques (Liu et al., 2023; Li et al., 2017; Hao et al., 2022). In this study, both MLR and MNLR models were utilized to assess the relationship between PM2.5 concentrations and key variables, including elevation, PM10 concentrations, and various meteorological conditions.
Multivariable Linear Regression (MLR), an extension of simple linear regression, is utilized to predict a response variable based on multiple explanatory variables. This study employed a MLR to assess the extent to which fine dust concentration can be estimated using meteorological factors. The regression function was accordingly adjusted to incorporate several predictor variables as follows:
(1)
In this statistical model, 𝑌 is designated as the dependent variable, with 𝑋0, X1 and additional terms representing the explanatory variables, also referred to as independent regressors. The termnotes the stochastic error component. It is observed that as the quantity of predictor variables escalates, the corresponding constants 𝛽 exhibit a proportional increase.
Multivariable Nonlinear Regression (MNLR) is a sophisticated statistical technique that enables the modeling of relationships between a dependent variable and one or more independent variables when linear representations cannot adequately capture such relationships. In this study, the regression model is formulated as a nonlinear function of the parameters in conjunction with one or more independent variables, allowing for greater flexibility and accuracy in representing complex data patterns.
Informed by the comprehensive analysis of MNLR in the studied (Yin et al., 2016), the model proposed in this study is delineated as follows:
(2)
where Y represents the dependent variable (PM2.5 dust), a is the constant, and
,
,
are the coefficients associated with the independent variables
, which include PM10 dust, temperature, wind speed, rainfall, humidity, and other relevant factors.
It is crucial to note that the considerations for interpreting the results of the nonlinear regression model are analogous to those applicable in the context of multivariate linear regression models.
2.3. Data
This study collected observational data during two distinct periods across two separate areas, with five concurrent samples obtained at each location. Table 2 describes the geographic coordinates and sampling heights for each sampling point. The collected data encompassed wind speed, humidity, temperature, and concentrations of PM10 and PM2.5. Additionally, surface data were gathered from the Environmental Protection Agency’s air quality monitoring station located at 105.799˚E, 21.015˚N.
Table 2. Information about PM2.5 sampling locations.
|
Name |
Height (m) |
Coordinates |
Time |
Area 1 |
CT 3-1 Me Tri Ha |
40 |
21˚00'56"N 105˚46'55"E |
Period 1: 25-31/08/2023. Period 2: 09-15/10/2023 |
CEO Tower |
102 |
21˚00'56"N 105˚46'58"E |
Vinhome Skylake |
151 |
21˚01'11"N 105˚46'53"E |
Keangnam |
212 |
21˚01'06"N 105˚47'04"E |
Keangnam |
336 |
21˚01'00"N 105˚47'03"E |
Surface |
The Sub-department of Environment Protection Hà Nội |
|
21˚00'54"N 105˚47'59"E |
Automatic Station |
Area 2 |
Sunny Hotel |
40 |
21˚01'40"N 105˚48'46"E |
Period 1: 04-10/09/2023. Period 2: 16-22/10/2023. |
VIT Tower |
80 |
21˚01'46"N 105˚48'42"E |
Lieu Giai Tower |
122 |
21˚02'07"N 105˚48'50"E |
Vinhome Metropolis |
160 |
21˚01'52"N 105˚48'55"E |
Lotte Center |
272 |
21˚01'56"N 105˚48'45"E |
3. Results
3.1. Trends in PM2.5 Concentrations
The trends in PM2.5 concentrations were analyzed for two investigation periods: from August 25, 2023, to September 10, 2023, and from October 4, 2023, to October 22, 2023. These periods fall within the autumn season in Northern Vietnam and align with the pre-harvest and post-harvest phases of the summer-fall rice crop in the surrounding areas of Hanoi. PM2.5 concentrations during these periods ranged from 9 to 109 μg/m3. The average concentration for the first period was 22.33 μg/m3 in Area 1 and 31.07 μg/m3 in Area 2, while for the second period, it increased to 33.39 μg/m3 and 31.93 μg/m3 for Area 1 and Area 2, respectively. These findings indicate higher fine particulate matter concentrations during the second period, likely due to the impact of agricultural residue burning, consistent with previous studies in Hanoi (Le et al., 2022; Vuong et al., 2023; Pham et al., 2024; Pham et al., 2021).
In this study, the majority of PM2.5 concentrations were within the lower range of the daily limit of 50 μg/m3 specified by QCVN 05:2013/BTNMT. However, these values exceeded the recommended 15 μg/m3 threshold established by the World Health Organization (WHO).
3.2. Distribution of PM2.5 Concentrations by Heights
Figure 2 shows the variations in meteorological factors, along with PM10 and PM2.5 concentrations by Height during two measurement periods in Area 1. In the first measurement, PM2.5 concentrations were observed to increase at a height of 40 m before decreasing sharply with further increases in height. This trend is attributed to the significant rise in wind speed at higher heights. At 40 m, the PM2.5 concentration was approximately 34.76 ± 2.38 μg/m3, which declined to 13.95 ± 1.7 μg/m3 at 336 m, representing a reduction by a factor of two.
A similar pattern was observed in Area 2 during the first measurement period, where PM2.5 concentrations also peaked at 40 m and then gradually decreased with increasing height (Figure 3). At 40 m, the concentration was approximately 33.92 ± 2.98 μg/m3, decreasing to 23.99 ± 4.53 μg/m3 at 272 m, a reduction of about 1.4 times. Differences in the ground-level measurement location and the PM2.5 monitoring methods can explain the concentration increase at 40 m. This trend of decreasing PM2.5 concentration with height aligns with findings from several other studies worldwide (Liu et al., 2018; Chan et al., 2005).
Figure 2. Measurement data and analysis in area 1.
PM10 and PM2.5 concentrations decreased with increasing height during the first period; however, the second period did not show significant variations in these values. Despite these differences, fine particulate matter concentrations in both periods consistently remained within the permissible environmental limits.
Similarly, Figure 3 presents the results for meteorological factors and fine dust concentrations by Height in Area 2. As observed in Area 1, discrepancies were noted between ground-level measurements and those at 40 m. Temperature, wind speed, and humidity varied at different heights, showing consistent trends across both periods. PM2.5 and PM10 concentrations in Area 2 were comparable to those in Area 1, with no significant differences in concentration levels.
Figure 3. Measurement data and analysis in area 2.
Meteorological factors are vital in influencing the distribution of PM10 and PM2.5 concentrations. Lower wind speeds contribute to poor dust dispersion in the air, and activities such as traffic and construction at ground level lead to higher dust concentrations. At greater heights, the influence of topographical obstructions diminishes, resulting in a gradual decrease in dust concentrations.
The measurements were conducted during a transitional period affected by light cold air masses, which emit higher dust values in the second period compared to the first. This pattern aligns with Hanoi’s climatic characteristics, as cold air masses hinder the dispersion of air layers, causing dust particles to remain stagnant.
A comparison of PM2.5 values with those from published studies shows variations across regions. For instance, in the research by Zauli-Sajani et al. (2018), PM2.5 concentrations at height during warm conditions were measured at 15 μg/m3, while cold conditions yielded 22 μg/m3. Similarly, the research by C.Y. Chan (Chan et al., 2005) in Beijing, which measured air quality at heights ranging from 8 to 325 m, found an average PM2.5 concentration of 65 μg/m3 and PM10 concentration of 150 μg/m3. These studies also demonstrate a similar trend of minimal concentration changes with height, consistent with the results observed in Hanoi.
3.3. Correlation between PM2.5 Concentrations and Factors
MLR and MNLR analyses were employed to assess the influence of variables such as height, temperature, wind speed, humidity, and PM10 concentration on PM2.5 levels. The detailed results of the MLR model are presented in Table 3, while the MNLR model results are shown in Table 4. The MLR analysis reveals that the selected factors account for approximately 50% to 80% of the variance in PM2.5 concentrations. In comparison, the MNLR model explains between 80% and 94% of the variance in PM2.5 concentrations. These findings demonstrate a clear influence of meteorological variables and PM10 concentrations on PM2.5 levels. Wind speed emerged as the most significant factor affecting PM2.5 concentrations, aligning with the findings of several previous studies (Vuong et al., 2023; Luong et al., 2021).
Table 3. The results of MLR model.
Coefficients of the MLR |
Area 1 |
Area 2 |
Period 1 |
Period 2 |
Period 1 |
Period 2 |
|
−17 |
−34 |
−23 |
−0.60 |
|
−0.02 |
0.02 |
0.02 |
0.02 |
|
0.12 |
1.23 |
0.56 |
0.08 |
|
0.64 |
−1.43 |
−0.49 |
−1.47 |
|
0.12 |
0.06 |
0.20 |
0.07 |
|
0.72 |
0.54 |
0.49 |
0.61 |
Table 4. The results of MNLR model.
Coefficients of the MNLR |
Area 1 |
Area 2 |
Period 1 |
Period 2 |
Period 1 |
Period 2 |
a |
−128 |
−32 |
3.47 |
−123 |
b0 |
0.23 |
0.06 |
1.58 |
−0.07 |
b1 |
11 |
1.45 |
−7.12 |
4.16 |
b2 |
12 |
4 |
−49 |
16 |
b3 |
−0.29 |
−0.04 |
3 |
1.76 |
b4 |
−3 |
0.25 |
−1.34 |
0.84 |
b10 |
0.00 |
0.00 |
−0.02 |
0.00 |
b20 |
0.00 |
0.00 |
−0.01 |
0.02 |
b30 |
0.00 |
0.00 |
−0.01 |
0.00 |
b40 |
0.00 |
0.00 |
0.00 |
0.00 |
b21 |
0.07 |
−0.11 |
1.16 |
−0.31 |
b31 |
−0.01 |
0.00 |
−0.02 |
−0.02 |
b41 |
0.06 |
0.00 |
0.03 |
−0.01 |
b32 |
−0.05 |
−0.01 |
0.13 |
0.02 |
b42 |
−0.03 |
0.00 |
0.02 |
−0.17 |
b43 |
0.02 |
0.01 |
0.01 |
0.00 |
b00 |
0.00 |
0.00 |
0.00 |
0.00 |
b11 |
−0.18 |
−0.02 |
0.14 |
−0.03 |
b22 |
−1.80 |
−0.21 |
1.75 |
−1.04 |
b33 |
0.00 |
0.00 |
−0.02 |
−0.01 |
b44 |
0.01 |
0.00 |
0.00 |
0.00 |
Figure 4. MAPE of the models. a) Area 1—Period 1; b) Area 1—Period 2; c) Area 2—Period 1; d) Area 2—Period 2.
Table 5. Comparison of the results from the linear regression model (MLR) and the nonlinear regression model (MNLR).
Area |
Phase |
R2 |
MAPE (error) |
MLR |
MNLR |
MLR |
MNLR |
1 |
1 |
0.795 |
0.892 |
21% |
26% |
2 |
0.787 |
0.923 |
24% |
21% |
2 |
1 |
0.517 |
0.815 |
53% |
16% |
2 |
0.883 |
0.942 |
28% |
9% |
Since PM2.5 is a subset of PM10, numerous researchers in previous studies have posited that PM2.5 concentration has a linear relationship with PM10 concentration (Makkonen et al., 2023). As a result, they have employed Multivariable Linear Regression (MLR) to estimate PM2.5 concentrations based solely on PM10 levels. However, as observed in Table 3 and Table 4, meteorological factors such as temperature, wind speed, and humidity also significantly influence PM2.5 concentrations with a proportional relationship with PM10. Given that meteorological conditions vary by season, their impact on PM concentrations fluctuates accordingly. Multivariable Nonlinear Regression (MNLR) was applied for each survey period. Table 5 presents a comparison between the MLR and MNLR models. The results indicate that through all periods, the adjusted R² values for each area and survey period using MNLR are consistently higher than those obtained with MLR. Except for the first period in Area 1, the MNLR model also yields substantially better Mean Absolute Percentage Error (MAPE) values compared to the MLR model. Thus, MNLR proves to be a more reliable approach for estimating PM2.5 concentrations. Figure 4 shows a comparison between the estimated PM2.5 concentrations from both models and the observed values at different heights, under corresponding meteorological conditions.
4. Conclusion
This study provided an overview of PM2.5 distribution by height across Hanoi, measured at multiple high-rise buildings following EPA standards. Most PM2.5 concentrations ranged between 30 and 38 μg/m3, lower than QCVN 05:2023 but more than double the WHO limit of 15 μg/m3.
PM2.5 distribution by height was uneven, with notable differences from ground level to over 300 m. In the first monitoring period, concentrations peaked at 40 m and decreased twofold from 34.76 μg/m3 at 40 m to 13.95 μg/m3 at 336 m. During the second period, there was less variation, likely due to cold air masses and reduced wind speeds.
MLR and MNLR models identified wind speed as the most significant factor influencing PM2.5 concentrations, explaining the minimal height variation in the second period.
The synchronization of meteorological conditions, along with the adjustment of corresponding height values across experimental phases, will more accurately elucidate the distribution of PM2.5 and PM10 concentrations as influenced by height variations.
Due to logistical constraints, measurements were taken simultaneously at five high-rise locations, which may have been affected by local meteorological conditions and building proximity. Future studies will aim to use technologies such as UAVs and Lidar for more accurate height-based assessments.