[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Issue
Volume 14, April
Previous Issue
Volume 14, February
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 

Atmosphere, Volume 14, Issue 3 (March 2023) – 189 articles

Cover Story (view full-size image): A numerical model (WRF-FDDA) has been implemented in the perspective of operational use for the nowcasting of the 10 m wind in the Venice Lagoon, which will be part of a real-time system aimed at the forecast of sea level height in the area. Two case studies were selected considering that the real-time forecasts missed their evolution, so that the sea level height was significantly underpredicted. They refer, respectively, to the passage of a small-scale cyclone across the lagoon and to the crossing of a frontal system from west to east. The comparison of the simulated wind with the observations shows that WRF-FDDA system represents a promising tool and a valuable support to the decisio makers for nowcasting in case of high tide in the Venice Lagoon. View this paper
  • Issues are regarded as officially published after their release is announced to the table of contents alert mailing list.
  • You may sign up for e-mail alerts to receive table of contents of newly released issues.
  • PDF is the official format for papers published in both, html and pdf forms. To view the papers in pdf format, click on the "PDF Full-text" link, and use the free Adobe Reader to open them.
Order results
Result details
Section
Select all
Export citation of selected articles as:
26 pages, 19193 KiB  
Article
Evaluating CMIP6 Historical Mean Precipitation over Africa and the Arabian Peninsula against Satellite-Based Observation
by Isaac Kwesi Nooni, Faustin Katchele Ogou, Abdoul Aziz Saidou Chaibou, Francis Mawuli Nakoty, Gnim Tchalim Gnitou and Jiao Lu
Atmosphere 2023, 14(3), 607; https://doi.org/10.3390/atmos14030607 - 22 Mar 2023
Cited by 11 | Viewed by 4151
Abstract
This study evaluated the historical precipitation simulations of 49 global climate models (GCMs) of the Coupled Model Intercomparison Project Phase 6 (CMIP6) in reproducing annual and seasonal precipitation climatology, linear trends, and their spatial correlation with global SST across Africa and the Arabian [...] Read more.
This study evaluated the historical precipitation simulations of 49 global climate models (GCMs) of the Coupled Model Intercomparison Project Phase 6 (CMIP6) in reproducing annual and seasonal precipitation climatology, linear trends, and their spatial correlation with global SST across Africa and the Arabian Peninsula during the period of 1980–2014, using Global Precipitation Climatology Centre (GPCP) data as a reference. Taylor’s diagram was used to quantify the strengths and weaknesses of the models in simulating precipitation. The CMIP6 multi-mean ensemble (MME) and the majority of the GCMs replicated the dominant features of the spatial and temporal variations reasonably well. The CMIP6 MME outperformed the majority of the individual models. The spatial variation of the CMIP6 MME closely matched the observation. The results showed that at annual and seasonal scales, the GPCP and CMIP6 MME reproduced a coherent spatial pattern in terms of the magnitude of precipitation. The humid region received >300 mm and the arid region received <50 mm across Africa and the Arabian Peninsula. The models from the same modeling centers replicated the precipitation levels across different seasons and regions. The CMIP6 MME and the majority of the individual models overestimate (underestimate) in humid (arid and semi-arid)-climate zones. The annual and pre-monsoon seasons (i.e., DJFMA) were better replicated in the CMIP6 GCMs than in the monsoon-precipitation model (MJJASON). The CMIP6 MME (GPCP) showed stronger wetting (drying) trends in the northern hemisphere. In contrast, a strong drying trend in the CMIP6 MME and a weak wetting trend in the GPCP were shown in the Southern Hemisphere. The CMIP6 MME captures the spatial pattern of linear trends better than individual models across different climate zones and regions. The relationship between precipitation and sea-surface temperature (SST) exhibited a high spatial correlation (−0.80 and 0.80) with large variability across different regions and climate zones. The GPCP (CMIP6 MME) exhibited a heterogenous (homogeneous) spatial pattern, with higher correlation coefficients recorded in the CMIP6 MME than in the GPCP in all cases. Individual models from the same modeling centers showed spatial homogeneity in correlation values. The differences exhibited by the individual GCMs highlight the significance of each model’s unique dynamics and physics; however, model selection should be considered for specific applications. Full article
(This article belongs to the Special Issue Precipitation in Africa)
Show Figures

Figure 1

Figure 1
<p>Digital elevation model (DEM) of the African continent and Arabian Peninsula.</p>
Full article ">Figure 2
<p>Taylor diagram comparing PREC observation (GPCP) with MME and 49 models (CMIP6) during 1980–2014. (<b>a</b>) Annual, (<b>b</b>) MJJASON, (<b>c</b>) December, (<b>d</b>) JFMA.</p>
Full article ">Figure 3
<p>Spatial distribution of multi-year mean precipitation in GPCP (<b>a1</b>–<b>a4</b>), and MME (<b>b1</b>–<b>b4</b>) for annual (<b>a1</b>,<b>b1</b>), JFMA (<b>a2</b>,<b>b2</b>), MJJASON (<b>a3</b>,<b>b3</b>), and December (<b>a4</b>,<b>b4</b>) periods during 1980–2014. The unit is mm year<sup>−1</sup>.</p>
Full article ">Figure 4
<p>Spatial distribution of multi-year mean precipitation in 49 selected individual models. The unit is mm year<sup>−1</sup>.</p>
Full article ">Figure 5
<p>Interannual variations of observation (i.e., GPCP, black line) in CMIP6 MME- and 49 CMIP6-simulated precipitation anomalies during 1980–2014 across Africa and the Arabian Peninsula. Anomalies were calculated using the mean over the period.</p>
Full article ">Figure 6
<p>Annual cycle of precipitation (mm) in GPCP (dotted black line) and CMIP6 models (colored lines) averaged over the period from 1980 to 2014.</p>
Full article ">Figure 7
<p>Linear trends of GPCP (<b>a1</b>–<b>a4</b>) (left), and MME (<b>b1</b>–<b>b4</b>) (right) at annual scale (<b>a1</b>,<b>b1</b>), JFMA (<b>a2</b>,<b>b2</b>), MJJASON (<b>a3</b>,<b>b3</b>), and December (<b>a4</b>,<b>b4</b>). The unit is mm decade<sup>−1</sup>. The black dots indicate that the trend passes the 0.05 significance test.</p>
Full article ">Figure 8
<p>Heatmap plots comparing correlation coefficients of mean precipitation between GPCP and 49 GCMs and their (<b>a</b>) annual, (<b>b</b>), MJJASON, (<b>c</b>) December, and (<b>d</b>) JFMA MME during 1980–2014.</p>
Full article ">Figure 9
<p>Correlation coefficient between observation GPCP and GCMs and the corresponding sea-surface temperature (SST) at annual scale (<b>a1</b>,<b>b1</b>), MJJASON (<b>a2</b>,<b>b2</b>), December (<b>a3</b>,<b>b3</b>), and JFMA (<b>a4</b>,<b>b4</b>). from 1980 to 2014. Hatched area indicates a 95% confidence level.</p>
Full article ">Figure 10
<p>Spatial distributions of correlation coefficients of GPCP and CMIP6 MME PREC and SST (i.e., annual scale (<b>a1</b>,<b>b1</b>), JFMA (<b>a2</b>,<b>b2</b>), MJJASON (<b>a3</b>,<b>b3</b>), and December (<b>a4</b>,<b>b4</b>)). during 1980–2014. Hatched area indicates a 95% confidence level.</p>
Full article ">Figure 11
<p>Spatial distributions of annual correlation coefficients of 49 selected CMIP6 models and SST during 1980–2014. Hatched area indicates a 95% confidence level.</p>
Full article ">
19 pages, 2649 KiB  
Article
Tree Traits and Microclimatic Conditions Determine Cooling Benefits of Urban Trees
by Mahmuda Sharmin, Mark G. Tjoelker, Sebastian Pfautsch, Manuel Esperón-Rodriguez, Paul D. Rymer and Sally A. Power
Atmosphere 2023, 14(3), 606; https://doi.org/10.3390/atmos14030606 - 22 Mar 2023
Cited by 10 | Viewed by 5044
Abstract
Trees play a key role in mitigating urban heat by cooling the local environment. This study evaluated the extent to which street trees can reduce sub-canopy air temperature relative to ambient conditions (ΔT), and how ΔT relates to tree traits and microclimatic variables. [...] Read more.
Trees play a key role in mitigating urban heat by cooling the local environment. This study evaluated the extent to which street trees can reduce sub-canopy air temperature relative to ambient conditions (ΔT), and how ΔT relates to tree traits and microclimatic variables. Air temperature under the canopies of 10 species was recorded within residential areas in Western Sydney, Australia, during summer 2019–2020. Tree and canopy traits, namely tree height, specific leaf area, leaf dry matter content, leaf area index, crown width and the Huber value (the ratio of sapwood area to leaf area) were then measured for all species. Species differed significantly in their ΔT values, with peak cooling (maximum ΔT −3.9 °C) observed between 9–10 am and sub-canopy warming (i.e., positive ΔT values) typically occurring during afternoon and overnight. Trees with high LAI and wider canopies were associated with the greatest daytime cooling benefits and lower levels of nighttime warming. ΔT was also negatively related to windspeed and vapor pressure deficit, and positively to solar irradiance. This study provides valuable information on how tree characteristics and microclimate influence potential cooling benefits that may aid planning decisions on the use of trees to mitigate heat in urban landscapes. Full article
(This article belongs to the Special Issue Strategies for Mitigation and Adaptation to Urban Heat)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Example of the urban morphology of one of the studied suburbs, Cranebrook, (<b>b</b>) street trees (<span class="html-italic">Eucalyptus microcorys</span> located on Arafura Avenue, Cranebrook), and (<b>c</b>) sensor and shielding affixed to a tree. The street map of Cranebrook shows that trees are growing on nature strips (verges) and at varying distances from houses. Images of (<b>a</b>,<b>b</b>) were sourced from google maps.</p>
Full article ">Figure 2
<p>Correlation between air temperatures recorded by three weather stations (EucFace, Richmond and Cranebrook); (<b>a</b>,<b>b</b>) maximum and minimum temperatures of Richmond BOM vs maximum and minimum temperatures from the EucFace meteorological station; (<b>c</b>,<b>d</b>) maximum and minimum temperatures of Richmond BOM vs maximum and minimum temperatures of Cranebrook BOM; (<b>e</b>,<b>f</b>) maximum and minimum temperatures of Cranebrook BOM vs maximum and minimum temperatures of the EucFace meteorological station.</p>
Full article ">Figure 3
<p>Diurnal variations in ΔT averaged across 37 summer days for individual trees (n = 96) (<b>a</b>) and species averages (9 to 10 replicates for each species) (<b>b</b>). The dashed red lines indicate no difference between the sub-canopy and ambient air temperatures. Circle colours correspond to tree species. For abbreviations of species names and number of trees for each species see <a href="#atmosphere-14-00606-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 4
<p>Species differences in ΔT (<b>a</b>) in the morning (07:00 to 13:00), (<b>b</b>) in the afternoon (15:00 to 19:00), and (<b>c</b>) at night (01:00 to 05:00). Different letters indicate significant differences at <span class="html-italic">p</span> &lt; 0.05. The black points indicate mean values (±SE). For abbreviations of species names and number of replicates see <a href="#atmosphere-14-00606-t001" class="html-table">Table 1</a>. Chi-square (X<sup>2</sup>) and <span class="html-italic">p</span> values for each of the <span class="html-italic">lmer</span> models are given in the figures. Analysis of variance (ANOVA) results of <span class="html-italic">lmer</span> models for each plot are in <a href="#app1-atmosphere-14-00606" class="html-app">Supplementary Table S1</a>.</p>
Full article ">Figure 5
<p>Relationship between ΔT and (<b>a</b>) vapor pressure deficit (VPD), (<b>b</b>) solar irradiance, (<b>c</b>) wind speed, (<b>d</b>) canopy width, (<b>e</b>) Huber value and (<b>f</b>) leaf area index (LAI) in the morning (07:00 to 13:00 h), using all recorded data from 96 individual trees over a three-month period. Total number of observations was 127,822. The dashed red lines indicate no difference between the canopy and ambient air temperatures and solid blue lines indicate best-fit trend predictions. Only significant relationships are shown. For detailed results of the <span class="html-italic">lmer</span> model see <a href="#app1-atmosphere-14-00606" class="html-app">Supplementary Table S2</a>.</p>
Full article ">Figure 6
<p>Relationship between ΔT and (<b>a</b>) vapor pressure deficit (VPD), (<b>b</b>) wind speed and (<b>c</b>) leaf area index (LAI) in the afternoon, using all recorded data from 96 individual trees over a three-month period. The dashed red lines indicate no difference between the sub-canopy and ambient air temperatures and soild blue lines correspond to regression slopes. There was a total of 127,822 data points. For detailed results of the <span class="html-italic">lmer</span> model see <a href="#app1-atmosphere-14-00606" class="html-app">Supplementary Table S3</a>.</p>
Full article ">Figure 7
<p>Relationships between ΔT and (<b>a</b>) vapor pressure deficit (VPD), (<b>b</b>) wind speed and (<b>c</b>) leaf area index at nighttime, using all recorded data from 96 individual trees over a three-month period. Only the first three strongest relationships with high chi-square (<span class="html-italic">X</span><sup>2</sup>) values are shown. The dashed red lines indicate no difference between the sub-canopy and ambient air temperatures and the solid blue lines represent regression slopes. There was a total of 106,560 data points. Only significant relationships are shown. For detailed results of the <span class="html-italic">lmer</span> model see <a href="#app1-atmosphere-14-00606" class="html-app">Supplementary Table S4</a>.</p>
Full article ">
11 pages, 277 KiB  
Article
Effects of Feeding a Commercial Starch Binding Agent during Heat Stress on Enteric Methane Emission, Rumen Volatile Fatty Acid Contents, and Diet Digestibility of Merino Lambs
by Pragna Prathap, Surinder S. Chauhan, Brian J. Leury, Jeremy J. Cottrell, Aleena Joy, Minghao Zhang and Frank R. Dunshea
Atmosphere 2023, 14(3), 605; https://doi.org/10.3390/atmos14030605 - 22 Mar 2023
Cited by 2 | Viewed by 1663
Abstract
Twenty-four Merino lambs were allocated to three dietary treatment groups to determine the effects of a dietary starch and protein binding agent and heat stress on methane (CH4) emissions and rumen parameters. The diets were a wheat-based diet (WD), a 2% [...] Read more.
Twenty-four Merino lambs were allocated to three dietary treatment groups to determine the effects of a dietary starch and protein binding agent and heat stress on methane (CH4) emissions and rumen parameters. The diets were a wheat-based diet (WD), a 2% Bioprotect™ (Bioprotect™, RealisticAgri, Rutland, UK) treated wheat-based diet (BD), and a maize-based diet (MD) for 3 periods of 1-week duration. During Period 1 (P1) the lambs were maintained under thermoneutral conditions and at a 1.7 × Maintenance (M) level. During P2 and P3, the lambs were maintained under cyclic heat stress conditions and fed at 1.7 × M and 2.0 × M, respectively. Total CH4 production was lower for the BD diet than the WD diet, which in turn was lower than the MD diet (p < 0.001). Total CH4 production was lower during P2 than P1 with P3 intermediate (p = 0.04). Rumen total volatile fatty acid (TVFA) concentrations were higher for the WD diet than the MD diet with the BD diet intermediate (p = 0.01). Rumen TVFA concentrations were lower during heat stress than under thermoneutral conditions (p < 0.001). Whole tract starch digestibility was higher for the BD and WD diets than the MD diet (p < 0.001). In conclusion, feeding Merino lambs the BD diet reduces CH4 emissions without reducing starch digestibility. Full article
(This article belongs to the Special Issue Observations and Management of Livestock Production Emissions)
16 pages, 2353 KiB  
Article
Understanding Temporal Patterns and Determinants of Ground-Level Ozone
by Junshun Wang, Jin Dong, Jingxian Guo, Panli Cai, Runkui Li, Xiaoping Zhang, Qun Xu and Xianfeng Song
Atmosphere 2023, 14(3), 604; https://doi.org/10.3390/atmos14030604 - 22 Mar 2023
Viewed by 2512
Abstract
Ground-level ozone pollution causes adverse health effects, and the detailed influences of meteorological factors and precursors on ozone at an hourly scale need to be further understood. We conducted an in-depth analysis of the phase relationships and periods of ground-level ozone in Shunyi [...] Read more.
Ground-level ozone pollution causes adverse health effects, and the detailed influences of meteorological factors and precursors on ozone at an hourly scale need to be further understood. We conducted an in-depth analysis of the phase relationships and periods of ground-level ozone in Shunyi station, Beijing, and contributing factors using wavelet analysis and geographic detectors in 2019. The combined effects of different factors on ozone were also calculated. We found that temperature had the strongest influence on ozone, and they were in phase over time. NO2 had the greatest explanatory power for the temporal variations in ozone among precursors. The wavelet power spectrum indicated that ozone had a periodic effect on multiple time scales, the most significant being the 22–26 h period. The wavelet coherence spectrum showed that in January–March and October–December, NO2 and ozone had an antiphase relationship, largely complementary to the in-phase relationship of temperature and ozone. Thus, the main influencing factors varied during the year. The interactions of temperature with NO2 significantly affected the temporal variations in ozone, and explanatory power surpassed 70%. The findings can deepen understanding of the effects of meteorological factors and precursors on ozone and provide suggestions for mitigating ozone pollution. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Ozone concentration determinants and their proxies.</p>
Full article ">Figure 2
<p>Z-score normalization time series and wavelet power spectrum for (<b>a</b>,<b>b</b>) O<sub>3</sub>; (<b>c</b>,<b>d</b>) NO<sub>2</sub>; (<b>e</b>,<b>f</b>) surface pressure; (<b>g</b>,<b>h</b>) 2 m temperature; (<b>i</b>,<b>j</b>) 10 m V wind; and (<b>k</b>,<b>l</b>) net surface net solar radiation.</p>
Full article ">Figure 3
<p>(<b>a</b>) Cross wavelet spectrum (XWT); (<b>b</b>) wavelet coherence spectrum (WTC) of NO<sub>2</sub> and O<sub>3</sub>.</p>
Full article ">Figure 4
<p>(<b>a</b>) Cross wavelet spectrum (XWT); (<b>b</b>) wavelet coherence spectrum (WTC) of 2 m temperature and O<sub>3</sub>.</p>
Full article ">
18 pages, 5519 KiB  
Article
Simulation and Analysis of the Influence of Sounding Rocket Outgassing on In-Situ Atmospheric Detection
by Zhiliang Zhang, Yueqiang Sun, Yongping Li, Jiangzhao Ai, Xiaoliang Zheng and Wei Wang
Atmosphere 2023, 14(3), 603; https://doi.org/10.3390/atmos14030603 - 22 Mar 2023
Cited by 1 | Viewed by 1844
Abstract
The Meridian Project’s sounding rocket mission uses a mass spectrometer to conduct in-situ atmospheric detection. In order to assess the influence of surface material outgassing and the attitude control jet on the spectrometer’s detection, a sounding rocket platform was modeled and simulated. Using [...] Read more.
The Meridian Project’s sounding rocket mission uses a mass spectrometer to conduct in-situ atmospheric detection. In order to assess the influence of surface material outgassing and the attitude control jet on the spectrometer’s detection, a sounding rocket platform was modeled and simulated. Using the physical field simulation software COMSOL and the Monte Carlo method, this study investigated whether the gas molecules from the two cases could enter the in-situ atmospheric mass spectrometer’s sensor sampling port after colliding with the background atmosphere. The simulation results show that the influence of surface material outgassing on the in-situ atmospheric detection is very small, even under the conditions of medium solar activity and medium geomagnetic activity, while the influence of the attitude control jet on the in-situ atmospheric detection is large but can be reduced by reducing the low-altitude attitude control operation and decreasing the transmission probability. Through simulation optimization and according to engineering needs, increasing the nozzle outlet cross-sectional area, increasing the temperature of the gas used for attitude control, increasing the nozzle rotation angle, increasing the nozzle outlet angle, or increasing the nozzle center height can reduce the transmission probability. This model can simulate and analyze the influence of both surface material outgassing and attitude control jets on in-situ atmospheric detection, optimize relevant parameters, and provide new ideas for relevant work. Full article
(This article belongs to the Section Upper Atmosphere)
Show Figures

Figure 1

Figure 1
<p>Diagram of the incident molecular motion.</p>
Full article ">Figure 2
<p>Structure of the de Laval nozzle. The three straight arrows are the inlet, throat and outlet from left to right, where <span class="html-italic">A</span><sup>*</sup> is the cross-sectional area of throat and the speed is 1 <span class="html-italic">Ma</span> at this point. The only curved arrow indicates the outlet angle <span class="html-italic">θ</span>.</p>
Full article ">Figure 3
<p>Relationship between <span class="html-italic">Ma</span> and <span class="html-italic">A</span><sup>*</sup>/<span class="html-italic">A</span>(<span class="html-italic">x</span>). The x-axis represents the Mach number, which is the ratio of the jet speed to the local speed of sound (the blue line intersects the x-axis at <math display="inline"><semantics> <mrow> <msqrt> <mn>6</mn> </msqrt> </mrow> </semantics></math>, and the green line intersects the x-axis at <math display="inline"><semantics> <mrow> <mrow> <mrow> <msqrt> <mrow> <mn>4389</mn> </mrow> </msqrt> </mrow> <mo>/</mo> <mrow> <mn>33</mn> </mrow> </mrow> </mrow> </semantics></math>). The y-axis represents the ratio of the cross-sectional area of the throat to any cross-sectional area of the nozzle.</p>
Full article ">Figure 4
<p>Flowchart of building model.</p>
Full article ">Figure 5
<p>Structure of the model: (<b>a</b>) the sounding rocket platform; (<b>b</b>) the sensor sampling port of the mass spectrometer for in-situ atmospheric detection; (<b>c</b>) the outlet of the de Laval nozzle for the attitude control jet, which is the blue area; and (<b>d</b>) the collision domain used to simulate the collision process.</p>
Full article ">Figure 6
<p>Simulation results of surface material outgassing in the MM case: (<b>a</b>) 120 km; (<b>b</b>) 320 km.</p>
Full article ">Figure 7
<p>Simulation results of the attitude control jet under different environments. The x-axis represents the gradual intensification of solar activity and geomagnetic activity, where F10.7 represents the degree of solar activity and Ap represents the degree of geomagnetic activity. The y-axis represents the transmission probability, which is the ratio of the number of particles received at the sensor sampling port to the total number of particles.</p>
Full article ">Figure 8
<p>Simulation results of the attitude control jet in the LL case: (<b>a</b>) 120 km; (<b>b</b>) 320 km.</p>
Full article ">Figure 9
<p>Relationship between the speed of jetting N<sub>2</sub> and transmission probability. The simulation results are indicated by circles in the line graph. The primary x-axis at the bottom indicates the speed of jetting N<sub>2</sub>, and the secondary x-axis at the top indicates the nozzle outlet radius at the corresponding speed. The y-axis indicates the transmission probability, which is the ratio of the number of particles received at the sensor sampling port to the total number of particles.</p>
Full article ">Figure 10
<p>Relationship between the speed of jetting He and transmission probability. The simulation results are indicated by stars in the line graph. The primary x-axis at the bottom indicates the speed of jetting He, and the secondary x-axis at the top indicates the nozzle outlet radius at the corresponding speed. The y-axis indicates the transmission probability, which is the ratio of the number of particles received at the sensor sampling port to the total number of particles.</p>
Full article ">Figure 11
<p>Relationship between speed of jetting N<sub>2</sub> and transmission probability. The simulation results are indicated by triangles in the line graph. The primary x-axis at the bottom indicates the speed of jetting N<sub>2</sub>, and the secondary x-axis at the top indicates the temperature at the corresponding speed. The y-axis indicates the transmission probability, which is the ratio of the number of particles received at the sensor sampling port to the total number of particles.</p>
Full article ">Figure 12
<p>Relationship between nozzle rotation angle and transmission probability. The simulation results are indicated by squares in the line graph. The x-axis represents the nozzle rotation angle, which is the angle between the outer normal vector of the nozzle outlet and the x-axis in <a href="#atmosphere-14-00603-f005" class="html-fig">Figure 5</a>. The y-axis represents the transmission probability, which is the ratio of the number of particles received at the sensor sampling port to the total number of particles.</p>
Full article ">Figure 13
<p>Relationship between nozzle outlet angle and transmission probability. The simulation results are indicated by pentagrams in the line graph. The x-axis represents the nozzle outlet angle, as shown in <a href="#atmosphere-14-00603-f002" class="html-fig">Figure 2</a> and <a href="#atmosphere-14-00603-t001" class="html-table">Table 1</a>. The y-axis represents the transmission probability, which is the ratio of the number of particles received at the sensor sampling port to the total number of particles.</p>
Full article ">Figure 14
<p>Relationship between nozzle center height and transmission probability. The simulation results are indicated by x-marks in the line graph. The x-axis represents the nozzle center height, as shown in <a href="#atmosphere-14-00603-f005" class="html-fig">Figure 5</a>c, which is the location of the nozzle. The y-axis represents the transmission probability, which is the ratio of the number of particles received at the sensor sampling port to the total number of particles.</p>
Full article ">
14 pages, 2112 KiB  
Article
Screening of Absorbents for Viscose Fiber CS2 Waste Air and Absorption–Desorption Process
by Ruixue Xiao, Kefan Chao, Ju Liu, Muhua Chen, Xinbao Zhu and Bo Fu
Atmosphere 2023, 14(3), 602; https://doi.org/10.3390/atmos14030602 - 22 Mar 2023
Cited by 2 | Viewed by 2084
Abstract
Screening of absorbents is essential for improving the removal rate of carbon disulfide (CS2) waste air by absorption. In this work, the UNIFAC model in Aspen Plus was utilized to calculate the excess Gibbs function and absorption potential of the binary [...] Read more.
Screening of absorbents is essential for improving the removal rate of carbon disulfide (CS2) waste air by absorption. In this work, the UNIFAC model in Aspen Plus was utilized to calculate the excess Gibbs function and absorption potential of the binary system of CS2 with various alcohols, ethers, esters, amines, and aromatic hydrocarbons. The results were used to quantitatively compare the efficiency of each solvent for CS2 absorption. The theoretical predictions were then verified by absorption experiments in a packed tower. The results showed that the performance of various solvents to CS2 roughly followed the order of esters < alcohols < amines < heavy aromatics < glycol ethers. Meanwhile, N-methyl-2-pyrrolidone (NMP) is the optimal absorbent for CS2 waste air treatment. Additionally, the process parameters of absorption and desorption of NMP were optimized. The results illustrated that the average mass removal efficiency of CS2 by NMP is 95.2% under following conditions: liquid–gas ratio of 3.75 L·m−3, a temperature of 20 °C, and inlet concentration lower than 10,000 mg·m−3. Under the conditions of 115 °C, 10 kPa, and a desorption time of 45 min, the average desorption rate of CS2 is 99.6%, and the average water content after desorption is 0.39%. Furthermore, the recycled lean liquid can maintain an excellent CS2 purification effect during the recycling process. Full article
(This article belongs to the Section Air Pollution Control)
Show Figures

Figure 1

Figure 1
<p>CS<sub>2</sub> standard curve.</p>
Full article ">Figure 2
<p>Absorption experiment setup: 1, 2, 12—rotameter; 3, 4, 10, 13—valve; 5, 9, 16, 17—stop valve; 6—CS<sub>2</sub> storage tank; 7—peristaltic pump; 8—mixing tank; 11, 15, 18—sampling test point; 14—absorption tower.</p>
Full article ">Figure 3
<p>Average mass removal efficiency of CS<sub>2</sub> waste air by six different absorbents.</p>
Full article ">Figure 4
<p>Absorption process research: (<b>a</b>) The effect of liquid-gas ratio on absorption, (<b>b</b>) The effect of temperature on absorption, (<b>c</b>) The effect of inlet concentration on absorption, (<b>d</b>) Repeated absorption experiments.</p>
Full article ">Figure 5
<p>Desorption process research: (<b>a</b>) The effect of water content on absorption (CS<sub>2</sub>-free), (<b>b</b>) The effect of CS<sub>2</sub> content on absorption (water-free), (<b>c</b>) The effect of water and CS<sub>2</sub> content on absorption (0.5‰ CS<sub>2</sub>), (<b>d</b>) The effect of temperature on desorption, (<b>e</b>) The effect of pressure on desorption, (<b>f</b>) The effect of time on desorption.</p>
Full article ">Figure 6
<p>Absorbent regeneration experiment: (<b>a</b>) Repeated desorption experiments, (<b>b</b>) Regenerated lean liquid circulation absorption experiment.</p>
Full article ">
16 pages, 4882 KiB  
Article
Seismo Ionospheric Anomalies around and over the Epicenters of Pakistan Earthquakes
by Munawar Shah, Rasim Shahzad, Muhsan Ehsan, Bushra Ghaffar, Irfan Ullah, Punyawi Jamjareegulgarn and Ahmed M. Hassan
Atmosphere 2023, 14(3), 601; https://doi.org/10.3390/atmos14030601 - 22 Mar 2023
Cited by 11 | Viewed by 3618
Abstract
Global Navigation Satellite System (GNSS)-based ionospheric anomalies are nowadays used to identify a possible earthquake (EQ) precursor and hence a new research topic in seismic studies. The current study also aims to provide an investigation of ionospheric anomalies associated to EQs. In order [...] Read more.
Global Navigation Satellite System (GNSS)-based ionospheric anomalies are nowadays used to identify a possible earthquake (EQ) precursor and hence a new research topic in seismic studies. The current study also aims to provide an investigation of ionospheric anomalies associated to EQs. In order to study possible pre-and post-seismic perturbations during the preparation phase of large-magnitude EQs, statistical and machine learning algorithms are applied to Total Electron Content (TEC) from the Global Positioning System (GPS) and Global Ionosphere Maps (GIMs). We observed TEC perturbation from the Sukkur (27.8° N, 68.9° E) GNSS station near the epicenter of Mw 5.4 Mirpur EQ within 5–10 days before the main shock day by implementing machine learning and statistical analysis. However, no TEC anomaly occurred in GIM-TEC over the Mirpur EQ epicenter. Furthermore, machine learning and statistical techniques are also implemented on GIM TEC data before and after the Mw 7.7 Awaran, where TEC anomalies can be clearly seen within 5–10 days before the seismic day and the subsequent rise in TEC during the 2 days after the main shock. These variations are also evident in GIM maps over the Awaran EQ epicenter. The findings point towards a large emission of EQ energy before and after the main shock during quiet storm days, which aid in the development of lithosphere ionosphere coupling. However, the entire analysis can be expanded to more satellite and ground-based measurements in Pakistan and other countries to reveal the pattern of air ionization from the epicenter through the atmosphere to the ionosphere. Full article
(This article belongs to the Special Issue Structure and Dynamics of Mesosphere and Lower Thermosphere)
Show Figures

Figure 1

Figure 1
<p>Study area map with red stars showing the location of the Awaran EQ (24 September 2013) epicenter and the Mirpur EQ (24 September 2013), while orange lines represent fault lines of Pakistan and green lines represent Pakistan regional boundaries.</p>
Full article ">Figure 2
<p>Temporal TEC variations with upper and lower confidence bounds of Mirpur EQ.</p>
Full article ">Figure 3
<p>Spatial variations of TEC from GIMs over the Mirpur EQ epicenter on 22 September 2019 (2 days before the main shock).</p>
Full article ">Figure 4
<p>(<b>a</b>) Comparison between GPS VTEC and NARX-predicted VTEC of Mirpur EQ. (<b>b</b>) Values exceeding the bounds depicting the anomalous behavior. The ‘*’ is for multiplication.</p>
Full article ">Figure 5
<p>(<b>a</b>) Comparison between GPS VTEC and MLP-predicted VTEC of the Mirpur EQ. (<b>b</b>) Values exceeding the bounds depicting the anomalous behavior. The ‘*’ is for multiplication.</p>
Full article ">Figure 6
<p>Temporal TEC variations with upper and lower confidence bounds of Awaran EQ.</p>
Full article ">Figure 7
<p>Spatial variations of TEC from GIMs over the Awaran EQ epicenter on 24 September 2013 (1 day after the main shock).</p>
Full article ">Figure 8
<p>(<b>a</b>) Comparison between GPS VTEC and NARX-predicted VTEC of Awaran EQ. (<b>b</b>) Values exceeding the bounds depicting the anomalous behavior. The ‘*’ is for multiplication.</p>
Full article ">Figure 9
<p>(<b>a</b>) Comparison between GPS VTEC and MLP-predicted VTEC of Awaran EQ. (<b>b</b>) Values exceeding the bounds depicting the anomalous behavior.</p>
Full article ">Figure 10
<p>Geomagnetic indices including Kp, Dst F10.7, and AE around the main shock day, where (<b>a</b>–<b>d</b>) represent Mirpur EQ and (<b>a’</b>–<b>d’</b>) for Awaran EQ.</p>
Full article ">
18 pages, 3450 KiB  
Article
Improvement of Maximum Air Temperature Forecasts Using a Stacking Ensemble Technique
by Linna Zhao, Shu Lu and Dan Qi
Atmosphere 2023, 14(3), 600; https://doi.org/10.3390/atmos14030600 - 21 Mar 2023
Cited by 3 | Viewed by 2092
Abstract
Due to the influence of complex factors such as atmospheric dynamic processes, physical processes and local topography and geomorphology, the prediction of near-surface meteorological elements in the numerical weather model often has deviation. The deep learning neural networks are more flexible but with [...] Read more.
Due to the influence of complex factors such as atmospheric dynamic processes, physical processes and local topography and geomorphology, the prediction of near-surface meteorological elements in the numerical weather model often has deviation. The deep learning neural networks are more flexible but with high variance. Here, we proposed a stacking ensemble model named FLT, which consists of a fully connected neural network with embedded layers (ED-FCNN), a long short-term memory (LSTM) network and a temporal convolutional network (TCN) to overcome the high variance of a single neural network and to improve prediction of maximum air temperature. The case study of daily maximum temperature forecast evaluated with observation of almost 2400 weather stations shows substantial improvement over that of single neural network model, ECMWF-IFS and statistical post-processing model. The FLT model can more effectively improve the forecast bias of the ECMWF-IFS model than that of any of the above single neural network model, with the RMSE reduced by 52.36% and the accuracy of temperature forecast increased by 43.12% compared with the ECMWF-IFS model. The average RMSEs of the FLT model decreases by 8.39%, 1.50%, 2.96% and 16.03%, respectively, compared with ED-FCNN, LSTM, TCN and the decaying average method. Full article
(This article belongs to the Section Meteorology)
Show Figures

Figure 1

Figure 1
<p>The illustration of the research area and the location of weather stations (black dots).</p>
Full article ">Figure 2
<p>Architectures of the ensemble neural network model (FLT) constructed in this study.</p>
Full article ">Figure 3
<p>Flow chart of the proposed FLT base on three neural network models by stacking generalization technique.</p>
Full article ">Figure 4
<p>Architectures of the sub-neural network model in an ensemble neural network (<b>a</b>) long short-term memory (LSTM) and (<b>b</b>) temporal convolutional network (TCN).</p>
Full article ">Figure 5
<p>Schematic design of time sliding window in the construction of time-series data set.</p>
Full article ">Figure 6
<p>Spatial distributions of average RMSE of daily maximum air temperature forecasts at all stations using different neural work models and the ECMWF-IFS model (unit: °C; the ring chart in the lower left corner represents the proportion of the number of stations (less than 6% are not shown) with different RMSE interval values to the total number of stations). (<b>a</b>) ECMWF-IFS, (<b>b</b>) DA, (<b>c</b>) ED-FCNN, (<b>d</b>) LSTM, (<b>e</b>) TCN and (<b>f</b>) FLT.</p>
Full article ">Figure 7
<p>The average RMSE at all stations over the Chinese mainland on each calendar day during the validation period using five post-processing models and the ECMWF-IFS model (the average RMSE of each model is also shown in the top; unit: °C).</p>
Full article ">Figure 8
<p>Bivariate of temperature forecast accuracy for (<b>a</b>) ED-FCNN, (<b>b</b>) LSTM, (<b>c</b>) TCN, (<b>d</b>) FLT and (<b>e</b>) DA model relative to the ECMWF-IFS model (horizontal coordinates are the ATF of the daily maximum temperature by the ECMWF-IFS model on the test set, vertical coordinates represent that of neural network models and DA model, red solid lines represent the contour of temperature forecast accuracy, blue solid lines are diagonal lines, the histogram at the top of the main figure is the distribution of temperature forecast accuracy of the ECMWF-IFS model, and the histogram at the right of the main figure is the frequency of temperature forecast accuracy of each neural network model).</p>
Full article ">
15 pages, 5201 KiB  
Article
Improving Intra-Urban Prediction of Atmospheric Fine Particles Using a Hybrid Deep Learning Approach
by Zhengyu Zhang, Jiuchun Ren and Yunhua Chang
Atmosphere 2023, 14(3), 599; https://doi.org/10.3390/atmos14030599 - 21 Mar 2023
Cited by 3 | Viewed by 1823
Abstract
Growing evidence links intra-urban gradients in atmospheric fine particles (PM2.5), a complex and variable cocktail of toxic chemicals, to adverse health outcomes. Here, we propose an improved hierarchical deep learning model framework to estimate the hourly variation of PM2.5 mass [...] Read more.
Growing evidence links intra-urban gradients in atmospheric fine particles (PM2.5), a complex and variable cocktail of toxic chemicals, to adverse health outcomes. Here, we propose an improved hierarchical deep learning model framework to estimate the hourly variation of PM2.5 mass concentration at the street level. By using a full-year monitoring data (including meteorological parameters, hourly concentrations of PM2.5, and gaseous precursors) from multiple stations in Shanghai, the largest city in China, as a training dataset, we first apply a convolutional neural network to obtain cross-domain and time-series features so that the inherent features of air quality and meteorological data associated with PM2.5 can be effectively extracted. Next, a Gaussian weight calculation layer is used to determine the potential interaction effects between different regions and neighboring stations. Finally, a long and short-term memory model layer is used to efficiently extract the temporal evolution characteristics of PM2.5 concentrations from the previous output layer. Further comparative analysis reveals that our proposed model framework significantly outperforms previous benchmark methods in terms of the stability and accuracy of PM2.5 prediction, which has important implications for the intra-urban health assessment of PM2.5-related pollution exposures. Full article
(This article belongs to the Special Issue Investigate Secondary Aerosol Formation and Source by Stable Isotopes)
Show Figures

Figure 1

Figure 1
<p>TOC art.</p>
Full article ">Figure 2
<p>Architecture overview of the HDL-learning framework.</p>
Full article ">Figure 3
<p>CNN structure of the model.</p>
Full article ">Figure 4
<p>Spatial effect of neighbors.</p>
Full article ">Figure 5
<p>Vector merge operation.</p>
Full article ">Figure 6
<p>Fitting trends of the different models. (<b>a</b>) LSTM; (<b>b</b>) ARIMA; (<b>c</b>) SVR; (<b>d</b>) HDL-learning.</p>
Full article ">Figure 7
<p>Fitting trends of different stations: (<b>a</b>) Chaunsha; (<b>b</b>) Zhangjiang; (<b>c</b>) Huinan; (<b>d</b>) Pudong; (<b>e</b>) Zhoupu; (<b>f</b>) Lingang.</p>
Full article ">Figure 7 Cont.
<p>Fitting trends of different stations: (<b>a</b>) Chaunsha; (<b>b</b>) Zhangjiang; (<b>c</b>) Huinan; (<b>d</b>) Pudong; (<b>e</b>) Zhoupu; (<b>f</b>) Lingang.</p>
Full article ">Figure 8
<p>The distance matrix of six counties.</p>
Full article ">Figure 9
<p>The weight matrix of six counties.</p>
Full article ">Figure 10
<p>Scatter plot showing the model performance over several locations in Shanghai. (<b>a</b>) Chuansha corr = 0.873; (<b>b</b>) Zhangjiang corr = 0.891; (<b>c</b>) Huinan corr = 0.954; (<b>d</b>) Pudong corr = 0.965; (<b>e</b>) Zhoupu corr = 0.948; (<b>f</b>) Lingang corr = 0.932.</p>
Full article ">Figure 10 Cont.
<p>Scatter plot showing the model performance over several locations in Shanghai. (<b>a</b>) Chuansha corr = 0.873; (<b>b</b>) Zhangjiang corr = 0.891; (<b>c</b>) Huinan corr = 0.954; (<b>d</b>) Pudong corr = 0.965; (<b>e</b>) Zhoupu corr = 0.948; (<b>f</b>) Lingang corr = 0.932.</p>
Full article ">
13 pages, 4503 KiB  
Article
High Wet Deposition of Black Carbon over the Sichuan Basin of China
by Yu Zhou, Xiaolin Zhang and Yuanzhi Wang
Atmosphere 2023, 14(3), 598; https://doi.org/10.3390/atmos14030598 - 21 Mar 2023
Cited by 1 | Viewed by 1776
Abstract
The wet deposition flux of black carbon (BC) over the Sichuan Basin is studied on the basis of the MERRA-2 data from 1981 to 2020, aiming at investigating high BC wet deposition flux in China in terms of long-term spatial-temporal trends and influences [...] Read more.
The wet deposition flux of black carbon (BC) over the Sichuan Basin is studied on the basis of the MERRA-2 data from 1981 to 2020, aiming at investigating high BC wet deposition flux in China in terms of long-term spatial-temporal trends and influences of BC column mass density and precipitation. In China, the largest BC wet deposition flux with a regionally-averaged value of 1.00 × 10−2 μg m−2 s−1 over the Sichuan Basin is observed, especially in the western and southern regions of the Basin with values as high as 2.20 × 10−2 μg m−2 s−1. The seasonality of BC wet deposition flux over the Sichuan Basin depicts maximum levels in autumn, moderate levels in spring and winter, and minimum levels in summer. The monthly mean BC wet deposition flux varies almost twofold, ranging from the lowest average value of 8.05 × 10−3 μg m−2 s−1 in July to the highest 1.28 × 10−2 μg m−2 s−1 in October. This study suggests that BC column mass density and precipitation are two significant factors affecting high BC wet deposition flux, whereas BC wet deposition flux is more related to BC column mass density than to precipitation over the Sichuan Basin. Full article
(This article belongs to the Special Issue Carbonaceous Aerosols)
Show Figures

Figure 1

Figure 1
<p>Location of the site selected in this study (<b>a</b>). Black box denotes the study area of the Sichuan basin (28° N–32° N, 103° E–108° E) in China (<b>b</b>). Red boxes denote the study areas with high BC wet deposition values over the western (region A) and southern (region B) parts of Sichuan Basin. Unit: μg m<sup>−2</sup> s<sup>−1</sup>.</p>
Full article ">Figure 2
<p>Spatial distribution of decade mean black carbon wet deposition flux over the study area in the Sichuan Basin during 1981–1990 (<b>a</b>), 1991–2000 (<b>b</b>), 2001–2010 (<b>c</b>), and 2011–2020 (<b>d</b>). Unit: μg m<sup>−2</sup> s<sup>−1</sup>.</p>
Full article ">Figure 3
<p>Annual variations of area-averaged black carbon wet deposition flux over the Sichuan basin from 1981 to 2020 (<b>a</b>) and in four periods (<b>b</b>) during 1981–1990, 1991–2000, 2001–2010 and 2011–2020.</p>
Full article ">Figure 4
<p>Mean spatial distribution of black carbon wet deposition flux in spring (<b>a</b>), summer (<b>b</b>), autumn (<b>c</b>), and winter (<b>d</b>) over the Sichuan Basin during 1981–2020. Unit: μg m<sup>−2</sup> s<sup>−1</sup>.</p>
Full article ">Figure 5
<p>Seasonal variations of area-averaged black carbon wet deposition flux over the Sichuan Basin from 1981 to 2020 (MAM, JJA, SON and DJF represent spring, summer, autumn, and winter, respectively) (<b>a</b>) and in four periods (<b>b</b>) during 1981–1990, 1991–2000, 2001–2010, and 2011–2020.</p>
Full article ">Figure 6
<p>Spatial distribution of monthly variations of black carbon wet deposition flux over the Sichuan Basin during 1981–2020. Unit: μg m<sup>−2</sup> s<sup>−1</sup>.</p>
Full article ">Figure 7
<p>Monthly variations of area-averaged black carbon wet deposition flux over the Sichuan basin from 1981 to 2020 (<b>a</b>) and in four periods (<b>b</b>) during 1981–1990, 1991–2000, 2001–2010, and 2011–2020.</p>
Full article ">Figure 8
<p>Relations between BC wet deposition flux and BC column mass density (<b>a</b>), and precipitation (<b>b</b>).</p>
Full article ">Figure 9
<p>The spatial distribution of BC wet deposition flux (unit: μg m<sup>−2</sup> s<sup>−1</sup>) (<b>a</b>), and spatial distributions of correlation coefficients of BC wet deposition flux with BC column mass density (<b>b</b>) and precipitation (<b>c</b>) over the Sichuan basin.</p>
Full article ">Figure 10
<p>The wet deposition-concentration-precipitation (WETD-C-P) relationships. The relationships of monthly BC wet deposition flux with precipitation and BC column mass density in region A (28° N–31° N, 103° E–104° E) (<b>a</b>) and region B (28° N–29° N, 104° E–106° E) (<b>b</b>).</p>
Full article ">
14 pages, 716 KiB  
Article
The Relationship between Ambient Fine Particulate Matter (PM2.5) Pollution and Depression: An Analysis of Data from 185 Countries
by Ravi Philip Rajkumar
Atmosphere 2023, 14(3), 597; https://doi.org/10.3390/atmos14030597 - 21 Mar 2023
Cited by 5 | Viewed by 2668
Abstract
Several studies have identified a relationship between air pollution and depression, particularly in relation to fine particulate matter (PM2.5) exposure. However, the strength of this association appears to be moderated by variables such as age, gender, genetic vulnerability, physical activity, and [...] Read more.
Several studies have identified a relationship between air pollution and depression, particularly in relation to fine particulate matter (PM2.5) exposure. However, the strength of this association appears to be moderated by variables such as age, gender, genetic vulnerability, physical activity, and climatic conditions, and has not been assessed at a cross-national level to date. Moreover, certain studies in this field have yielded negative results, and there are discrepancies between the results obtained in high-income countries and those from low- and middle-income countries. The current study examines cross-sectional and longitudinal associations between the incidence of depression in each country, based on Global Burden of Disease Study data, and the average national level of PM2.5 based on the World Health Organization’s database, over the past decade (2010–2019). The observed associations were adjusted for age, gender, level of physical activity, income, education, population density, climate, and type of depression. It was observed that while PM2.5 levels showed significant cross-sectional associations with the incidence of depression, longitudinal analyses were not suggestive of a direct causal relationship. These findings are discussed in the light of recent contradictory results in this field, and the need to consider the intermediate roles of a number of individual and environmental factors. Full article
Show Figures

Figure 1

Figure 1
<p>Scatter plots of the relationship between average annual PM<sub>2.5</sub> and the incidence of depression in 2010 and 2019.</p>
Full article ">Figure 2
<p>The complex nature of the relationship between PM<sub>2.5</sub> and depression.</p>
Full article ">
26 pages, 6859 KiB  
Article
The Effect of Abrupt Changes to Sources of PM10 and PM2.5 Concentrations in Three Major Agglomerations in Mexico
by Jorge Méndez-Astudillo and Ernesto Caetano
Atmosphere 2023, 14(3), 596; https://doi.org/10.3390/atmos14030596 - 21 Mar 2023
Viewed by 2307
Abstract
In the three major urban agglomerations in Mexico (Mexico City, Monterrey, and Guadalajara), a significant change to anthropogenic sources of air pollution happened in March–May 2020, when policies implemented to stop the spread of the COVID-19 virus in Mexico caused the reduction of [...] Read more.
In the three major urban agglomerations in Mexico (Mexico City, Monterrey, and Guadalajara), a significant change to anthropogenic sources of air pollution happened in March–May 2020, when policies implemented to stop the spread of the COVID-19 virus in Mexico caused the reduction of some anthropogenic sources of air pollution. We study the effect of these significant changes to air pollution sources using satellite-retrieved aerosol optical depth (AOD) and particulate matter (PM10 and PM2.5) concentrations from ground stations. The Chow test was applied to study trend changes in PM concentrations from 1 January to 30 May 2020. The Mann–Whitney non-parametric test was then used to compare average PM concentrations in April and May pre-lockdown, during lockdown in 2020, and post-lockdown in 2021. The assessment was further performed by evaluating the exceedance of national air quality standard maxima. The trend analysis showed that PM10 concentrations were reduced during lockdown in Mexico City and Monterrey, whereas no change was found for PM10 in Guadalajara and PM2.5 in the three cities. Further analysis showed that in Mexico City and Guadalajara, average PM10 and PM2.5 concentrations decreased by 12% in April and May 2020. However, in Monterrey, average PM10 and PM2.5 concentrations increased by 2.76% and 11.07%, respectively, in April 2021 due to a severe drought that caused dry soils and dust around the city. The results of this research can be used to implement policies for reducing anthropogenic sources to improve the air quality in urban areas. Full article
(This article belongs to the Section Air Pollution Control)
Show Figures

Figure 1

Figure 1
<p>Locations of Guadalajara, Monterrey, and Mexico City, and locations of the stations used in each city. The light grey areas are urban areas, whereas dark green areas are forest and/or mountains. All stations report meteorological and PM data. See <a href="#atmosphere-14-00596-t0A1" class="html-table">Table A1</a> in the <a href="#app1-atmosphere-14-00596" class="html-app">Appendix A</a> for more information about the stations.</p>
Full article ">Figure 2
<p>Flow diagram of the methodology of this study.</p>
Full article ">Figure 3
<p>AOD (dimensionless) in GMC, MAG, and MTY in April and May 2015–2022.</p>
Full article ">Figure 4
<p>PM concentrations at (<b>a</b>) GMC, (<b>b</b>) MTY, and (<b>c</b>) MAG from 1 January to 31 May 2020. The lockdown period is marked in yellow (from 23 March 2020, DOY 83).</p>
Full article ">Figure 5
<p>Representative meteograms showing the relations between wind speed and PM at stations in (<b>a</b>) Mexico City, (<b>b</b>) Monterrey, and (<b>c</b>) Guadalajara. Red bars show days that the PM<sub>2.5</sub> limit was exceeded, and blue bars show exceedance of the PM<sub>10</sub> limit. In the figure, PM<sub>10</sub> is represented by a blue line and PM<sub>2.5</sub> by a red line. The wind speed and direction are shown in the bottom of each panel.</p>
Full article ">Figure 6
<p>Number of days that exceeded the PM<sub>10</sub> and PM<sub>2.5</sub> limits in Mexico City (<b>a</b>,<b>b</b>). Number of days that exceeded the PM<sub>10</sub> and PM<sub>2.5</sub> limits in Monterrey (<b>c</b>,<b>d</b>). Number of days that exceeded the PM<sub>10</sub> limits in Guadalajara (<b>e</b>). In all cases, the period is April and May 2015–2021.</p>
Full article ">Figure A1
<p>AOD distribution in GMC, MTY, and MAG in May 2019–2021.</p>
Full article ">Figure A2
<p>(<b>a</b>) Surface wind maps in Mexico City (red point) at daytime and nighttime. (<b>b</b>) Surface wind maps in Guadalajara (red point) at daytime and nighttime. (<b>c</b>) Surface wind maps in Monterrey (red point) at daytime and nighttime.</p>
Full article ">Figure A2 Cont.
<p>(<b>a</b>) Surface wind maps in Mexico City (red point) at daytime and nighttime. (<b>b</b>) Surface wind maps in Guadalajara (red point) at daytime and nighttime. (<b>c</b>) Surface wind maps in Monterrey (red point) at daytime and nighttime.</p>
Full article ">Figure A3
<p>Meteograms for stations in Monterrey with a significant relation between wind speed and PM concentrations. Blue bars represent days over the PM<sub>10</sub> maximum established by the national standard for air quality. The yellow bars indicate exceeding both PM<sub>10</sub> and PM<sub>2.5</sub> maxima.</p>
Full article ">Figure A4
<p>Meteograms for stations in Guadalajara with a significant relation between wind speed and PM concentrations. Blue bars represent days over the PM<sub>10</sub> maximum established by the national standard for air quality.</p>
Full article ">Figure A5
<p>Meteograms for stations in Mexico City with a significant relation between wind speed and PM concentrations. Blue bars represent days over the PM<sub>10</sub> maximum established by the national standard for air quality. The yellow bars indicate exceeding both PM<sub>10</sub> and PM<sub>2.5</sub> maxima, and the red bars show days over the limit for PM<sub>2.5</sub>.</p>
Full article ">Figure A5 Cont.
<p>Meteograms for stations in Mexico City with a significant relation between wind speed and PM concentrations. Blue bars represent days over the PM<sub>10</sub> maximum established by the national standard for air quality. The yellow bars indicate exceeding both PM<sub>10</sub> and PM<sub>2.5</sub> maxima, and the red bars show days over the limit for PM<sub>2.5</sub>.</p>
Full article ">Figure A5 Cont.
<p>Meteograms for stations in Mexico City with a significant relation between wind speed and PM concentrations. Blue bars represent days over the PM<sub>10</sub> maximum established by the national standard for air quality. The yellow bars indicate exceeding both PM<sub>10</sub> and PM<sub>2.5</sub> maxima, and the red bars show days over the limit for PM<sub>2.5</sub>.</p>
Full article ">
14 pages, 5684 KiB  
Article
Decadal Prediction of the Summer Extreme Precipitation over Southern China
by Huijie Wang, Yanyan Huang, Dapeng Zhang and Huijun Wang
Atmosphere 2023, 14(3), 595; https://doi.org/10.3390/atmos14030595 - 21 Mar 2023
Cited by 4 | Viewed by 1925
Abstract
The decadal variability of the summer extreme precipitation over southern China (EPSC) is remarkable, especially for the significant decadal enhancement after the 1990s. The study documented that the summer sea surface temperature (SST) over the North Atlantic and spring sea ice concentration (SIC) [...] Read more.
The decadal variability of the summer extreme precipitation over southern China (EPSC) is remarkable, especially for the significant decadal enhancement after the 1990s. The study documented that the summer sea surface temperature (SST) over the North Atlantic and spring sea ice concentration (SIC) over the East Siberian Sea can significantly affect the EPSC. The summer SST over the North Atlantic influences the low-pressure cyclone in the western Pacific by modulating the SST over the tropical Pacific, thus affecting EPSC. A decrease in the SIC of the East Siberian Sea induces a negative Arctic Oscillation, which induces the increased SST over northwest Pacific and the anomalous cyclone over there, in turn, affecting EPSC. Both predictors have a quasi-period of 10–14 years, which provides useful predictive signals for EPSC. The leading 7-year SST and the leading 5-year SIC are chosen to establish the prediction model based on the decadal increment method, which can well predict the EPSC, especially for the shift in the early 1990s. These results provide a clue to the limited predictability of decadal-scale extreme climate events. Full article
(This article belongs to the Special Issue Long-Term Variability of Atmospheric Precipitation)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) The spatial pattern of the leading EOF mode and (<b>b</b>) the corresponding normalized time series of the 5-year smoothing mean of the total summer extreme precipitation in East China during 1963–2018 (PC1, red line), as well as the time series of the EPSC (R95p, blue line) and DI_EPSC during 1966–2018. EPSC is defined as the 5-year running mean of regional averaged summer extreme precipitation over southern China. DI_EPSC is calculated as the 3-year decadal increment of the EPSC.</p>
Full article ">Figure 2
<p>Spatial pattern of correlation coefficients between the predictand (DI_EPSC) during the time period 1966–2018 and the predictors: (<b>a</b>) the 3-year increment of the summer SST 7 years ahead of the DI_EPSC during 1959–2011, (<b>b</b>) the 3-year increment of the spring SIC 5 years ahead of DI_EPSC during 1961–2013. The dotted regions indicate significant variability at the 95% confidence level based on the Student’s <span class="html-italic">t</span>-test. The rectangles indicate the area-weighted averaged regions of the predictors, including the DI_SST (26°–39° N, 43°–60° W minus 49°–57° N, 14°–49° W), and the DI_SIC (71°–77° N, 149°–157° E).</p>
Full article ">Figure 3
<p>(<b>a</b>) Wavelet analysis and (<b>b</b>) the autocorrelation coefficients of the DI_SST predictor. In Panel (<b>a</b>), the dotted regions indicate significant variability at the 95% confidence level. In Panel (<b>b</b>), the horizontal coordinates indicate the years lagged by the predictors and the dashed lines indicate significant variability at the 90% confidence level. Additionally, the effective sample size is 21.</p>
Full article ">Figure 4
<p>Correlation coefficients between the DI_SST with (<b>a</b>) the summer surface winds, (<b>b</b>) the SST, (<b>d</b>) the 850 hPa winds, and (<b>e</b>) the 500 hPa winds during 1959–2011. Correlation coefficients between the eastern Pacific SST (red rectangle region in Panel (<b>b</b>), SST is multiplied by −1) with (<b>c</b>) the mean SLP during 1959–2011. Correlation coefficients between the DI_SST during 1959–2011 with (<b>f</b>) the 500 hPa vertical velocity during 1966–2018. Thus, the negative value in (<b>f</b>) should indicate the descent motion due to the negative autocorrelation with a lag of 7 years. The dotted regions and shaded regions indicate significant variability at the 90% confidence level, and the red rectangle region in Panels (<b>d</b>–<b>f</b>) is the study area (southern China). The variables including SST, surface winds, SLP, 850 hPa winds, 500 hPa winds, and 500 hPa vertical velocity are in the form of a 3-year decadal increment.</p>
Full article ">Figure 5
<p>(<b>a</b>) Wavelet analysis and (<b>b</b>) the autocorrelation coefficients of the DI_SIC predictor. In Panel (<b>a</b>), the dotted regions indicate significant variability at the 95% confidence level. In Panel (<b>b</b>), the horizontal coordinates indicate the years lagged by the predictors and the dashed lines indicate significant variability at the 90% confidence level. Additionally, the effective sample size is 18.</p>
Full article ">Figure 6
<p>Correlation coefficients between the DI_SIC with (<b>a</b>) the spring SLP, (<b>b</b>) the SST, and (<b>d</b>) 500 hPa vertical velocity during 1961–2013. Correlation coefficients between the SST over North Pacific (red rectangle region of Panel (<b>b</b>)), with the summer (<b>c</b>) 850 hPa winds and (<b>e</b>) vertical integral of the divergence of the water vapor flux during 1961–2013. The dotted and shaded regions indicate significant variability at the 90% confidence level, and the red rectangle region of Panels (<b>c</b>,<b>d</b>) represent the study area (southern China). The variables including SST, SLP, 850 hPa winds, 500 hPa vertical velocity, and vertical integral of the divergence of the water vapor flux are in the form of a 3-year decadal increment.</p>
Full article ">Figure 7
<p>Time series of the DI_EPSC during 1966–2019 (bars), the leading 7-year DI_SST during 1959–2015 (red line), and the leading 5-year DI_SIC during 1961–2017 (blue line). The results in this figure were standardized.</p>
Full article ">Figure 8
<p>The results of the cross-validation for the period of 1966–2019 (<b>a</b>,<b>c</b>) and the independent hindcast for the period of 2009–2019 (<b>b</b>,<b>d</b>). The panels show the decadal increment of the EPSC (DI_EPSC) predicted by the statistical model (<b>a</b>,<b>b</b>); the final predicted EPSC, achieved by adding the predicted DI_EPSC to the observed EPSC at 3 years ago (<b>c</b>,<b>d</b>), where the light pink regions indicate significant variability at the 95% prediction interval. The results in this figure were standardized. Cor indicates the correlation coefficient between the observations and the prediction results, and the MSSS is the prediction skill of the statistical model.</p>
Full article ">Figure 9
<p>Spatial pattern of the anomaly extreme precipitation (unit: mm) over eastern China after a 5-year running mean in 2019 related to the climatology (1991–2020). The black rectangle region represents the study area (southern China).</p>
Full article ">
16 pages, 12761 KiB  
Article
Temporal Variation of NO2 and O3 in Rome (Italy) from Pandora and In Situ Measurements
by Annalisa Di Bernardino, Gabriele Mevi, Anna Maria Iannarelli, Serena Falasca, Alexander Cede, Martin Tiefengraber and Stefano Casadio
Atmosphere 2023, 14(3), 594; https://doi.org/10.3390/atmos14030594 - 21 Mar 2023
Cited by 7 | Viewed by 3081
Abstract
To assess the best measures for the improvement of air quality, it is crucial to investigate in situ and columnar pollution levels. In this study, ground-based measurements of nitrogen dioxide (NO2) and ozone (O3) collected in Rome (Italy) between [...] Read more.
To assess the best measures for the improvement of air quality, it is crucial to investigate in situ and columnar pollution levels. In this study, ground-based measurements of nitrogen dioxide (NO2) and ozone (O3) collected in Rome (Italy) between 2017 and 2022 are analyzed. Pandora sun-spectrometers provided the time series of the NO2 vertical column density (VC-NO2), tropospheric column density (TC-NO2), near-surface concentration (SC-NO2), and the O3 vertical column density (VC-O3). In situ concentrations of NO2 and O3 are provided by an urban background air quality station. The results show a clear reduction of NO2 over the years, thanks to the recent ecological transition policies, with marked seasonal variability, observable both by columnar and in situ data. Otherwise, O3 does not show inter-annual variations, although a clear seasonal cycle is detectable. The results suggest that the variation of in situ O3 is mainly imputable to photochemical reactions while, in the VC-O3, it is triggered by the predominant contribution of stratospheric O3. The outcomes highlight the importance of co-located in situ and columnar measurements in urban environments to investigate physical and chemical processes driving air pollution and to design tailored climate change adaptation strategies. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Geographical map of the area under investigation. The black dashed line depicts the urban center of Rome, and cyan and magenta markers show the positions of the Sapienza (SAP) and Villa Ada (VA) stations. (<b>b</b>) Aerial photograph of the area surrounding Sapienza (SAP). (<b>c</b>) Aerial photograph of the area surrounding Villa Ada (VA).</p>
Full article ">Figure 2
<p>Pan#138 vs. Pan#117 scatter plots for the period 24 July—11 September 2020 (<b>a</b>) VC-NO<sub>2</sub>, (<b>b</b>) TC-NO<sub>2</sub>, (<b>c</b>) SC-NO<sub>2</sub>, and (<b>d</b>) VC-O<sub>3</sub>. The reported statistical parameters are root-mean-square of differences (RMSE), linear Pearson correlation coefficient (R), mean bias (MB), and number of valid retrievals (n). Orange lines indicate the linear fit between the two instruments, while green lines indicate the ideal fit <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mi>x</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Weekly average time series of VC-NO<sub>2</sub> and TC-NO<sub>2</sub> for (<b>a</b>) weekdays and (<b>b</b>) weekends and SC-NO<sub>2</sub> and in situ NO<sub>2</sub> concentration for (<b>c</b>) weekdays and (<b>d</b>) weekends. VC-NO<sub>2</sub>, TC-NO<sub>2,</sub> and SC-NO<sub>2</sub> are measured at Sapienza (SAP) from Pandora. In situ NO<sub>2</sub> concentrations are measured at the Villa Ada (VA) station.</p>
Full article ">Figure 4
<p>Hourly and monthly variation of VC-NO<sub>2</sub> for (<b>a</b>) weekdays and (<b>b</b>) weekends over the period 2017–2022.</p>
Full article ">Figure 5
<p>Hourly and monthly variation of TC-NO<sub>2</sub> for (<b>a</b>) weekdays and (<b>b</b>) weekends over the period 2017–2022.</p>
Full article ">Figure 6
<p>Hourly and monthly variation of SC-NO<sub>2</sub> for (<b>a</b>) weekdays and (<b>b</b>) weekends over the period 2017–2022.</p>
Full article ">Figure 7
<p>Hourly and monthly variation of in situ NO<sub>2</sub> for (<b>a</b>) weekdays and (<b>b</b>) weekends over the period 2017–2022.</p>
Full article ">Figure 8
<p>Weekly average time series of VC-O<sub>3</sub> and in situ O<sub>3</sub> concentration for (<b>a</b>) weekdays and (<b>b</b>) weekends. VC-O<sub>3</sub> is measured at Sapienza (SAP) from Pandora. In situ O<sub>3</sub> concentration is measured at the Villa Ada (VA) station.</p>
Full article ">
11 pages, 1748 KiB  
Communication
Seasonal Characteristics of Fine Particulate Carbonaceous Species in Taiyuan, North China
by Wenbo Li, Xinya Zhao, Fengbo Guo and Kankan Liu
Atmosphere 2023, 14(3), 593; https://doi.org/10.3390/atmos14030593 - 21 Mar 2023
Cited by 2 | Viewed by 1626
Abstract
To characterize seasonal carbonaceous aerosol pollution in Taiyuan, a typical city in North China that mainly relies heavily on coal, a total of 124 PM2.5 samples were collected from August 2018 to the next May. The annual mean PM2.5 concentration was [...] Read more.
To characterize seasonal carbonaceous aerosol pollution in Taiyuan, a typical city in North China that mainly relies heavily on coal, a total of 124 PM2.5 samples were collected from August 2018 to the next May. The annual mean PM2.5 concentration was 83.8 ± 48.5 μg m−3, with a seasonal rank of winter (117.4 ± 47.6 μg m−3) > spring (79.2 ± 34.3 μg m−3) > fall (67.3 ± 34.7 μg m−3) > summer (31.8 ± 6.5 μg m−3), suggesting that fine particulate pollution was still serious in cold seasons. Organic carbon (OC) and elemental carbon (EC) showed similar seasonal patterns with PM2.5. The mean concentration values of OC in summer, fall, winter, and spring were 5.1 ± 0.9, 11.8 ± 6.4, 22.1 ± 14.9, and 12.2 ± 6.7 μg m−3, respectively. The mean concentration values of EC in summer, fall, winter, and spring were 1.5 ± 0.3, 2.5 ± 1.6, 4.4 ± 2.8, and 2.4 ± 1.5 μg m−3, respectively. The proportion of total carbon aerosol (TCA) was about 31.7%, 33.8%, 30.0%, and 27.0% in PM2.5 in summer, fall, winter, and spring, respectively. The good correlation between OC vs. EC and the high value of OC/EC suggests that coal and biomass combustion were the main emissions in cold seasons, aggravated by adverse meteorological conditions and the dustpan-shaped terrain. The mean annual secondary organic carbon (SOC) concentration was 6.1 ± 7.1μg m−3, representing 38.7% of the OC content. The present results presented the serious carbonaceous particulate pollution, which might affect haze pollution in cold seasons. Full article
(This article belongs to the Section Aerosols)
Show Figures

Figure 1

Figure 1
<p>The sampling site (red star) in downtown Taiyuan city with the Wusu international airport (WIA, blue plane) and the agrometeorological station (AS, green pin).</p>
Full article ">Figure 2
<p>Daily PM<sub>2.5</sub>, OC, EC, and percentage of TCA in PM<sub>2.5</sub> during the sampling period.</p>
Full article ">Figure 3
<p>Seasonal wind roses with the intensity of atmospheric temperature inversion.</p>
Full article ">Figure 4
<p>Seasonal correlations between OC and EC in PM<sub>2.5</sub>.</p>
Full article ">Figure 5
<p>Seasonal SOC/OC percentage as a function of daily hours of sunlight and average solar irradiance during daylight hours.</p>
Full article ">Figure 6
<p>Scatter plots and linear regression equations of SOC/OC percentage vs. solar irradiance.</p>
Full article ">
22 pages, 4551 KiB  
Article
Apportionment of PM2.5 Sources across Sites and Time Periods: An Application and Update for Detroit, Michigan
by Zhiyi Yang, Md Kamrul Islam, Tian Xia and Stuart Batterman
Atmosphere 2023, 14(3), 592; https://doi.org/10.3390/atmos14030592 - 20 Mar 2023
Cited by 5 | Viewed by 3564
Abstract
Identifying sources of air pollutants is essential for informing actions to reduce emissions, exposures, and adverse health impacts. This study updates and extends apportionments of particulate matter (PM2.5) in Detroit, MI, USA, an area with extensive industrial, vehicular, and construction activity [...] Read more.
Identifying sources of air pollutants is essential for informing actions to reduce emissions, exposures, and adverse health impacts. This study updates and extends apportionments of particulate matter (PM2.5) in Detroit, MI, USA, an area with extensive industrial, vehicular, and construction activity interspersed among vulnerable communities. We demonstrate an approach that uses positive matrix factorization models with combined spatially and temporally diverse datasets to assess source contributions, trend seasonal levels, and examine pandemic-related effects. The approach consolidates measurements from 2016 to 2021 collected at three sites. Most PM2.5 was due to mobile sources, secondary sulfate, and secondary nitrate; smaller contributions arose from soil/dust, ferrous and non-ferrous metals, and road salt sources. Several sources varied significantly by season and site. Pandemic-related changes were generally modest. Results of the consolidated models were more consistent with respect to trends and known sources, and the larger sample size should improve representativeness and stability. Compared to earlier apportionments, contributions of secondary sulfate and nitrate were lower, and mobile sources now represent the dominant PM2.5 contributor. We show the growing contribution of mobile sources, the need to update apportionments performed just 5–10 years ago, and that apportionments at a single site may not apply elsewhere in the same urban area, especially for local sources. Full article
(This article belongs to the Special Issue The Michigan-Ontario Ozone Source Experiment (MOOSE))
Show Figures

Figure 1

Figure 1
<p>Map showing study area, monitoring sites (red diamonds), and major point sources (numbered black circles). Apportionments obtained at the three speciation sites.</p>
Full article ">Figure 2
<p>Trends of monthly average concentrations for selected pollutants at the three sites. Plots use 3-month running averages. “All” is the three-site average.</p>
Full article ">Figure 3
<p>Factor profiles in final models. Four sets of profiles are shown: colored bars are for three individual sites using approach 1; hatched bar is profile for combined dataset using approach 2.</p>
Full article ">Figure 4
<p>Trends of monthly factor contributions for the two PMF approaches.</p>
Full article ">
17 pages, 2614 KiB  
Article
The Source Apportionment of Heavy Metals in Surface Dust in the Main District Bus Stops of Tianshui City Based on the Positive Matrix Factorization Model and Geo-Statistics
by Chunyan Li, Xinmin Wang, Shun Xiao and Hai Wang
Atmosphere 2023, 14(3), 591; https://doi.org/10.3390/atmos14030591 - 20 Mar 2023
Cited by 4 | Viewed by 2509
Abstract
For the pollution assessment and quantitative source apportionment of heavy metals in surface dust, a total of 52 surface dust samples were collected from bus stops in Tianshui City. The geoaccumulation index (Igeo) and potential ecological risk index (RI [...] Read more.
For the pollution assessment and quantitative source apportionment of heavy metals in surface dust, a total of 52 surface dust samples were collected from bus stops in Tianshui City. The geoaccumulation index (Igeo) and potential ecological risk index (RI) were used to analyze the pollution levels caused by heavy metals. The Positive Matrix Factorization (PMF) of the receptor modeling and geo-statistics were employed to analyze the source of the heavy metals. The results were as follows. ① Except for Mn, Co and V, the mean concentrations of other heavy metals have exceeded the local background value of Gansu. The percentage of excessive concentrations of Cu, Zn, Sr, Ba and Pb in the samples was 100%, and that of Cr, Ni and As were 96.15%, 94.23%, and 96.15%, respectively. ② Semivariogram model fitting showed that the block-based coefficients of Cu, Zn, Sr, Ba, Pb, Cr, Ni, and As were between 0.25 and 0.75, indicating that they were mainly affected by human factors. The high values of Pb, Zn, Ni and As were mainly distributed in the eastern part of the study area, and the high values of Cu, Sr, Ba and Cr were distributed in a spot-like pattern in the study area. ③ The Igeo results showed that As, Cu, Zn, and Pb were the main contamination factors, and the optimized RI showed that the heavy metals were the overall ecological risk of intensity, among which Pb, As and Cu were the main ecological factors and should be taken as the priority control objects. ④ Based on the PMF, there are four main sources of eleven heavy metals. V, Mn, and Co were attributed to natural sources, accounting for 18.33%; Cu, Sr, and Ba were from mixed sources of pollution from transportation and industrial alloy manufacturing, accounting for 26.99%; Cr and Ni were from sources of construction waste pollution, accounting for 17.17%, As, Zn and Pb were mainly produced by coal-traffic mixed pollution emissions, accounting for 37.52%. Overall, the study area was dominated by coal-traffic emissions. Full article
Show Figures

Figure 1

Figure 1
<p>A schematic map of the study area and sampling points location.</p>
Full article ">Figure 2
<p>The spatial distribution of the concentrations of heavy metals.</p>
Full article ">Figure 3
<p>The ecological risk index of heavy metals at different sampling points.</p>
Full article ">Figure 4
<p>The contribution rates of different sources by the PMF model.</p>
Full article ">Figure 5
<p>The factor contribution rates of heavy metals are calculated by the PMF.</p>
Full article ">
5 pages, 193 KiB  
Editorial
Bioaerosols: Composition, Meteorological Impact, and Transport
by Salvatore Romano
Atmosphere 2023, 14(3), 590; https://doi.org/10.3390/atmos14030590 - 20 Mar 2023
Cited by 5 | Viewed by 3245
Abstract
The characterization and the main properties of biogenic airborne particles (or bioaerosols) comprising both living and dead microorganisms (such as bacteria, fungi, viruses, pollen, and microbial fragments) are currently of increasing interest in the scientific community [...] Full article
(This article belongs to the Special Issue Bioaerosols: Composition, Meteorological Impact, and Transport)
15 pages, 3771 KiB  
Article
Statistical PM2.5 Prediction in an Urban Area Using Vertical Meteorological Factors
by Jutapas Saiohai, Surat Bualert, Thunyapat Thongyen, Kittichai Duangmal, Parkpoom Choomanee and Wladyslaw W. Szymanski
Atmosphere 2023, 14(3), 589; https://doi.org/10.3390/atmos14030589 - 19 Mar 2023
Cited by 15 | Viewed by 4113
Abstract
A key concern related to particulate air pollution is the development of an early warning system that can predict local PM2.5 levels and excessive PM2.5 concentration episodes using vertical meteorological factors. Machine learning (ML) algorithms, particularly those with recognition tasks, show [...] Read more.
A key concern related to particulate air pollution is the development of an early warning system that can predict local PM2.5 levels and excessive PM2.5 concentration episodes using vertical meteorological factors. Machine learning (ML) algorithms, particularly those with recognition tasks, show great potential for this purpose. The objective of this study was to compare the performance of multiple linear regression (MLR) and multilayer perceptron (MLP) in predicting PM2.5 levels. The software was trained to predict PM2.5 levels up to 7 days in advance using data from long-term measurements of vertical meteorological factors taken at five heights above ground level (AGL)—10, 30, 50, 75, and 110 m—and PM2.5 concentrations measured 30 m AGL. The data used were collected between 2015 and 2020 at the Microclimate and Air Pollutants Monitoring Tower station at Kasetsart University, Bangkok, Thailand. The results showed that the correlation coefficients of PM2.5 predicted and observed using MLR and MLP were in the range of 0.69–0.86 and 0.64–0.82, respectively, for 1–3 days ahead. Both models showed satisfactory agreement with the measured data, and MLR performed better than MLP at PM2.5 prediction. In conclusion, this study demonstrates that the proposed approach can be used as a component of an early warning system in cities, contributing to sustainable air quality management in urban areas. Full article
(This article belongs to the Special Issue Atmospheric Particulate Matter Hazard Mapping)
Show Figures

Figure 1

Figure 1
<p>Location of data acquisition instruments at Kasetsart University, Bangkok, Thailand.</p>
Full article ">Figure 2
<p>Monthly averages (vertical bars) and annual averages (horizontal bars) of PM<sub>2.5</sub> concentrations from 2015 to 2020 measured at Kasetsart University using TEOM.</p>
Full article ">Figure 3
<p>Daily averages (over 24 h) of PM<sub>2.5</sub> concentrations from 2015 to 2020 obtained at Kasetsart University using TEOM. Some days are missing data due to equipment maintenance.</p>
Full article ">Figure 4
<p>Polar plot of actual PM<sub>2.5</sub> concentrations measured at 30 m (<b>a</b>), averaged over the years 2015–2020, and the corresponding wind direction and wind speed distribution measured at the KU tower at the indicated levels (<b>b</b>).</p>
Full article ">Figure 5
<p>Quality of predictions (R) of PM<sub>2.5</sub> concentrations for up to 7 days using multiple linear regression compared with observed PM<sub>2.5</sub> data. Forward prediction for: (<b>a</b>) 3 h, (<b>b</b>) 12 h, (<b>c</b>) 1 day, (<b>d</b>) 2 days, (<b>e</b>) 2 days, (<b>f</b>) 7 days.</p>
Full article ">Figure 6
<p>Measured and predicted PM<sub>2.5</sub> concentrations obtained using multiple linear regression modeling. Forward prediction for: (<b>a</b>) 3 h, (<b>b</b>) 12 h, (<b>c</b>) 1 day, (<b>d</b>) 2 days, (<b>e</b>) 2 days, (<b>f</b>) 7 days.</p>
Full article ">Figure 7
<p>Quality of prediction (R) of PM<sub>2.5</sub> concentrations for up to 7 days using the multilayer perceptron approach compared with observed PM<sub>2.5</sub> data. Forward prediction for: (<b>a</b>) 3 h, (<b>b</b>) 12 h, (<b>c</b>) 1 day, (<b>d</b>) 2 days, (<b>e</b>) 2 days, (<b>f</b>) 7 days.</p>
Full article ">Figure 8
<p>Variations between observed PM<sub>2.5</sub> concentrations and PM<sub>2.5</sub> concentrations predicted using the multilayer perceptron method. Forward prediction for: (<b>a</b>) 3 h, (<b>b</b>) 12 h, (<b>c</b>) 1 day, (<b>d</b>) 2 days, (<b>e</b>) 2 days, (<b>f</b>) 7 days.</p>
Full article ">Figure 9
<p>Relative error between predicted PM<sub>2.5</sub> concentrations and those observed on 8 January and 9 January 2020, showing the accuracy of the developed method as a function of time.</p>
Full article ">
11 pages, 1667 KiB  
Article
Assessment of Natural Radioactivity and Radon Exhalation in Peruvian Gold Mine Tailings to Produce a Geopolymer Cement
by Rafael Liza, Patrizia Pereyra, Jose Rau, Maribel Guzman, Laszlo Sajo-Bohus and Daniel Palacios
Atmosphere 2023, 14(3), 588; https://doi.org/10.3390/atmos14030588 - 19 Mar 2023
Cited by 8 | Viewed by 3464
Abstract
Mining generates significant amounts of waste that can represent a source of contamination for areas close to the extraction area, generating a negative impact both on the environment and the health of people. This study aims to evaluate the radiological risk derived from [...] Read more.
Mining generates significant amounts of waste that can represent a source of contamination for areas close to the extraction area, generating a negative impact both on the environment and the health of people. This study aims to evaluate the radiological risk derived from exposure to natural radionuclides contained in tailings from Peruvian gold mines and to establish whether the tailings can be used as raw materials in building materials. The mine tailings come from a mining project in the northern highlands of Peru. Radon exhalation was measured using Rad7 in a closed chamber and activity concentration of 226Ra, 232Th, and 40K radioisotopes by gamma spectrometry using NaI 3” × 3” detector. Maximum activity concentrations measured for 226Ra and 232Th were 15.38 Bq kg−1 and 11.9 Bq kg−1, respectively; meanwhile, activity concentration for 40K ranged from 182.7 Bq kg−1 to 770.8 Bq kg−1. All activity concentrations were below the worldwide average except for 40K. The radon exhalation rate varied from 2.8 to 7.2 mBq kg−1 h−1. The gamma index (Iγ), and radiological parameters, including the Radium equivalent activity (Raeq), and the external hazard index (Hex), being below the recommended levels by UNSCEAR, ensure the safe use of these mines tailing to produce a geopolymer cement. Full article
(This article belongs to the Special Issue Radon and NORM: Impact on Air Quality)
Show Figures

Figure 1

Figure 1
<p>Gold mine location map in the northern highlands of Peru.</p>
Full article ">Figure 2
<p>Experimental setup to assess <sup>222</sup>Rn exhalation rates from mine tailings samples using the Rad7.</p>
Full article ">Figure 3
<p>Experimental <sup>222</sup>Rn concentration growth and exponential fit of measured radon activity concentration inside the accumulation chamber.</p>
Full article ">
24 pages, 3526 KiB  
Article
A Multipurpose Simulation Approach for Hybrid Electric Vehicles to Support the European CO2 Emissions Framework
by Alessandro Tansini, Georgios Fontaras and Federico Millo
Atmosphere 2023, 14(3), 587; https://doi.org/10.3390/atmos14030587 - 18 Mar 2023
Cited by 3 | Viewed by 3051
Abstract
Hybrid Electric Vehicles (HEVs) are a prominent solution for reducing CO2 emissions from transport in Europe. They are equipped with at least two propulsion energy converters, an Internal Combustion Engine (ICE) and one or more Electric Machines (EMs), operated in a way [...] Read more.
Hybrid Electric Vehicles (HEVs) are a prominent solution for reducing CO2 emissions from transport in Europe. They are equipped with at least two propulsion energy converters, an Internal Combustion Engine (ICE) and one or more Electric Machines (EMs), operated in a way to exploit synergies and achieve fuel efficiency. Because of the variety in configurations and strategies, the use of simulation is essential for vehicle development and characterisation of energy consumption. This paper introduces a novel simulation approach to estimate the CO2 emissions from different hybrid architectures (series, parallel, power-split) and electrification degrees (mild, full, plug-in and range extender) that is relatively simple, flexible and accurate. The approach identifies the optimal power split between the energy converters for any given time in a driving cycle according to three evaluation levels: supervisor, ICE manager and optimiser. The latter relies on the Equivalent Consumption Minimisation Strategy (ECMS) and the limitations imposed by the other two layers. Six light-duty HEVs with different hybrid architectures were tested to support the development of the approach. The results show an indicative accuracy of ±5%, enabling to run assessments of hybrid powertrain solutions and supporting regulatory and consumer information initiatives. Full article
(This article belongs to the Special Issue Vehicle Emissions: New Challenges and Potential Solutions)
Show Figures

Figure 1

Figure 1
<p>Generic hybrid powertrain model [<a href="#B22-atmosphere-14-00587" class="html-bibr">22</a>].</p>
Full article ">Figure 2
<p>EM efficiency and generic full load curve (<b>left</b>) and maximum efficiency as function of rated power (<b>right</b>). Right figure based on Larsson [<a href="#B24-atmosphere-14-00587" class="html-bibr">24</a>] and authors’ further elaboration.</p>
Full article ">Figure 3
<p>Generic Electric Power System model [<a href="#B22-atmosphere-14-00587" class="html-bibr">22</a>].</p>
Full article ">Figure 4
<p>Equivalent Circuit Battery Model—Battery cell modelling [<a href="#B22-atmosphere-14-00587" class="html-bibr">22</a>].</p>
Full article ">Figure 5
<p>SoC–OCV curves for different battery chemistries (<b>left</b>) [<a href="#B28-atmosphere-14-00587" class="html-bibr">28</a>] and battery cell internal resistance as a function of SoC and operation (<b>right</b>) [<a href="#B22-atmosphere-14-00587" class="html-bibr">22</a>]. Left figure based on Bharathraj et al. [<a href="#B28-atmosphere-14-00587" class="html-bibr">28</a>] and authors’ further elaboration.</p>
Full article ">Figure 6
<p>Supervisor logic [<a href="#B22-atmosphere-14-00587" class="html-bibr">22</a>].</p>
Full article ">Figure 7
<p>ICE Manager logic [<a href="#B22-atmosphere-14-00587" class="html-bibr">22</a>].</p>
Full article ">Figure 8
<p>Optimiser logic [<a href="#B22-atmosphere-14-00587" class="html-bibr">22</a>].</p>
Full article ">Figure 9
<p>Cost of electric energy as a function of SoC. Based on [<a href="#B31-atmosphere-14-00587" class="html-bibr">31</a>] and authors’ further elaboration.</p>
Full article ">Figure 10
<p>Comparison of absolute balanced CO<sub>2</sub> emissions between experimental and simulation results.</p>
Full article ">Figure 11
<p>Simulation results for (<b>a</b>): Vehicle 1 (series REx), (<b>b</b>): Vehicle 2 (parallel full hybrid), (<b>c</b>): Vehicle 3 (series-parallel plug-in hybrid), (<b>d</b>): Vehicle 4 (parallel plug-in hybrid), (<b>e</b>): Vehicle 5 (parallel mild hybrid), (<b>f</b>): Simulation results for Vehicle 6 (power-split full hybrid).</p>
Full article ">Figure 12
<p>Vehicle 4: RDE trip characteristics (<b>top</b>) and CO<sub>2</sub> emissions results (<b>bottom</b>).</p>
Full article ">Figure A1
<p>Powertrain configurations of the tested vehicles [<a href="#B22-atmosphere-14-00587" class="html-bibr">22</a>].</p>
Full article ">
12 pages, 7034 KiB  
Article
What Causes the Arabian Gulf Significant Summer Sea Surface Temperature Warming Trend?
by Kamal A. Alawad, Abdullah M. Al-Subhi, Mohammed A. Alsaafani and Turki M. Alraddadi
Atmosphere 2023, 14(3), 586; https://doi.org/10.3390/atmos14030586 - 18 Mar 2023
Cited by 5 | Viewed by 3297
Abstract
The present study investigated the significant sea surface temperature (SST) warming trend during the summer season over the Arabian Gulf (AG) and its links with the large-scale atmospheric driver, namely, the Atlantic multidecadal oscillation (AMO), from 1900 to 2021. The link between the [...] Read more.
The present study investigated the significant sea surface temperature (SST) warming trend during the summer season over the Arabian Gulf (AG) and its links with the large-scale atmospheric driver, namely, the Atlantic multidecadal oscillation (AMO), from 1900 to 2021. The link between the AMO and the AGs oceanic circulations has received little scientific attention. It has been found that there is a significant spatial positive trend, with a maximum of up to 0.6 °C per decade over the far northern end, while the time series trend shows a significant shift after 1995, with an average value of about 0.36 °C per decade. The spatial trend in the AG is eight times and four times higher than the global value from 1980 to 2005 using HadISST and OISST, respectively. The AMOs significant role in the AGs SST significant warming trend has been confirmed by the spatial and temporal correlation coefficient, which is above 0.50 and 0.48, respectively, with statistical significance at the 99% level. The underlying mechanisms that explained the AMO-related AGs SST decadal variability can be explained as follows: when the AMO is in a positive phase, the surface northwesterly wind weakens, leading to (1) less advection of the relatively cold air masses from Southern Europe to the AG and surrounding areas, (2) less evaporation, and thus less surface cooling (3); thus, this enhances the water masses stratification and decreases the upwelling process, and vice versa occurs for the negative phase. For the air temperature, the positive AMO phase coincides with the occurrence of warm air masses covering all of the Arabian Peninsula, North Africa, and Southern Europe. These processes prove that the AMO is a possible candidate for the AGs SST decadal variability, hence enabling a better evaluation of future climate scenarios for this important region. Our results provide initial insights into the AMO-driven spatiotemporal variability in the SST over the AG and prove that the relation is nonstationary over time. Further analyses are required to explore whether the impacts of the AMO are extended to other oceanic variables such as evaporation rate, heat transport, etc. Full article
(This article belongs to the Section Meteorology)
Show Figures

Figure 1

Figure 1
<p>Arabian Gulf (AG) location and bathymetry in (m).</p>
Full article ">Figure 2
<p>Time series of annual SST (dot red line), summer (solid red line), Atlantic multidecadal oscillation (AMO) (solid blue line) during 1900–2020, and OISST (dot black line) during 1982–2020, and their associated correlation analysis.</p>
Full article ">Figure 3
<p>Summer climatology of the OISST and HadISST (1982–2018).</p>
Full article ">Figure 4
<p>The SST difference between pre- and post-1996.</p>
Full article ">Figure 5
<p>The spatial HadISST trend in °C decade<sup>−1</sup> during the summer season. The dot area indicates the significant area at 95%.</p>
Full article ">Figure 6
<p>The spatial correlation between the AMO time series and SST. The dot area indicates the significant area at 95%.</p>
Full article ">Figure 7
<p>The difference in the U-wind component between pre- and post-1996 (<b>a</b>), correlation analysis between the AMO time series and U-wind (<b>b</b>), and correlation analysis between the AGs SST time series and U-wind (<b>c</b>). All the data are from 1971–2021.</p>
Full article ">Figure 8
<p>Same as <a href="#atmosphere-14-00586-f007" class="html-fig">Figure 7</a>, but for the surface temperate.</p>
Full article ">Figure 9
<p>Same as <a href="#atmosphere-14-00586-f007" class="html-fig">Figure 7</a>, but for upper-air temperate at 850 hpa.</p>
Full article ">
23 pages, 4694 KiB  
Article
N2O Temporal Variability from the Middle Troposphere to the Middle Stratosphere Based on Airborne and Balloon-Borne Observations during the Period 1987–2018
by Gisèle Krysztofiak, Valéry Catoire, Thierry Dudok de Wit, Douglas E. Kinnison, A. R. Ravishankara, Vanessa Brocchi, Elliot Atlas, Heiko Bozem, Róisín Commane, Francesco D’Amato, Bruce Daube, Glenn S. Diskin, Andreas Engel, Felix Friedl-Vallon, Eric Hintsa, Dale F. Hurst, Peter Hoor, Fabrice Jegou, Kenneth W. Jucks, Armin Kleinböhl, Harry Küllmann, Eric A. Kort, Kathryn McKain, Fred L. Moore, Florian Obersteiner, Yenny Gonzalez Ramos, Tanja Schuck, Geoffrey C. Toon, Silvia Viciani, Gerald Wetzel, Jonathan Williams and Steven C. Wofsyadd Show full author list remove Hide full author list
Atmosphere 2023, 14(3), 585; https://doi.org/10.3390/atmos14030585 - 18 Mar 2023
Viewed by 3470
Abstract
Nitrous oxide (N2O) is the fourth most important greenhouse gas in the atmosphere and is considered the most important current source gas emission for global stratospheric ozone depletion (O3). It has natural and anthropogenic sources, mainly as an unintended [...] Read more.
Nitrous oxide (N2O) is the fourth most important greenhouse gas in the atmosphere and is considered the most important current source gas emission for global stratospheric ozone depletion (O3). It has natural and anthropogenic sources, mainly as an unintended by-product of food production activities. This work examines the identification and quantification of trends in the N2O concentration from the middle troposphere to the middle stratosphere (MTMS) by in situ and remote sensing observations. The temporal variability of N2O is addressed using a comprehensive dataset of in situ and remote sensing N2O concentrations based on aircraft and balloon measurements in the MTMS from 1987 to 2018. We determine N2O trends in the MTMS, based on observations. This consistent dataset was also used to study the N2O seasonal cycle to investigate the relationship between abundances and its emission sources through zonal means. The results show a long-term increase in global N2O concentration in the MTMS with an average of 0.89 ± 0.07 ppb/yr in the troposphere and 0.96 ± 0.15 ppb/yr in the stratosphere, consistent with 0.80 ppb/yr derived from ground-based measurements and 0.799 ± 0.024 ppb/yr ACE-FTS (Atmospheric Chemistry Experiment Fourier Transform Spectrometer) satellite measurements. Full article
(This article belongs to the Section Air Quality)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Spatial distribution of all aircraft and balloon data for pressures between 600 and 5 hPa. Seasons are represented by color codes: December-January-February (deep blue), March-April-May (light blue), June-July-August (green), and September-October-November (orange). (<b>b</b>) Pressure levels vs. latitude.</p>
Full article ">Figure 2
<p>N<sub>2</sub>O mole fraction (in volume mixing ratios with units of parts per billion) versus pressure for all the observations after rebinning all the concentrations onto the same pressure grid. Also shown are the two modes obtained by blind source separation with arbitrary amplitudes.</p>
Full article ">Figure 3
<p>Latitudinal dependence of the amplitude A of modes 1 and 2. The latitude is the median latitude during the profile.</p>
Full article ">Figure 4
<p>Number of measurements per 10° latitude grid box and 60 levels pressure (600–5 hPa) for all years and seasons: December-January-February (DJF), March-April-May (MAM), June-July-August (JJA), and September-October-November (SON). The dashed black line represents the tropopause pressure according to ECMWF ERA-interim [<a href="#B48-atmosphere-14-00585" class="html-bibr">48</a>].</p>
Full article ">Figure 5
<p>N<sub>2</sub>O trend in the upper troposphere (from +20 to +50 hPa compared to the tropopause pressure) versus time for the South Polar region (deep blue), southern midlatitudes (turquoise), tropics (red), northern midlatitudes (orange) and North Polar region (yellow). The error bars represent the standard deviation for each box (0.25-year range and pressure level). For each latitude band, the slope in ppb/yr is calculated by using linear regression, taking into account the standard deviation of each box; the standard deviation of the slope is also calculated (σ) along with the associated F value (Fisher’s test) and the critical t-value at 95% confidence. The more the F value exceeds the critical t-value, the more the trend can be considered linear.</p>
Full article ">Figure 6
<p>N<sub>2</sub>O trend in the UTLS from aircraft and balloon measurements versus pressure relative to the tropopause (tropopause altitude from ERA-interim climatology [<a href="#B48-atmosphere-14-00585" class="html-bibr">48</a>]). The color symbols represent the latitude band (deep blue circle: South Pole: 90° S–60° S; turquoise square: southern midlatitudes: 60° S–30° S; red diamond: tropical latitudes: 20° S–20° N; orange triangle: northern midlatitudes: 30° N–60° N; and yellow circle: North Pole: 60° N–90° N). The error bars represent one standard deviation on the slope for the plot concentration versus time for each pressure level. The tropospheric and stratospheric averages are represented by the black line and the grey areas represented the +/−1 sigma mean values. The trend according to WMO [<a href="#B9-atmosphere-14-00585" class="html-bibr">9</a>] mean surface data (0.80 ppb/yr) is represented by a black circle and for convenient visualization at +750 hPa pressure relative to the tropopause (WMO ground measurements: AGAGE in situ NOAA, flask and in situ, CSIRO, flask, WMO/GAW) [<a href="#B9-atmosphere-14-00585" class="html-bibr">9</a>]. ACE v4 mean (0.799 ± 0.024 ppb/yr) for latitude ranges ±60° and altitudes between 5.5 and 10.5 km is represented by a purple circle at +170 hPa relative to the tropopause [<a href="#B49-atmosphere-14-00585" class="html-bibr">49</a>].</p>
Full article ">Figure 7
<p>N<sub>2</sub>O mixing ratio zonal mean in 2018 (the trend observed in <a href="#atmosphere-14-00585-f005" class="html-fig">Figure 5</a> applied to the data) for December-January-February (DJF), March-April-May (MAM), June-July-August (JJA), and September-October-November (SON). Zonal means are calculated from error-weighted measurements (see <a href="#sec2dot3-atmosphere-14-00585" class="html-sec">Section 2.3</a>) with 10° latitude bins and 60 levels of pressure (from 600 hPa to 5 hPa). The dashed line represents the tropopause pressure according to ECMWF ERA-interim [<a href="#B48-atmosphere-14-00585" class="html-bibr">48</a>].</p>
Full article ">Figure 8
<p>Global tropospheric N<sub>2</sub>O VMR (600–400 hPa) versus years from measurements (weighted average of aircraft and balloon measurements) with error bars representing standard deviation. The dashed line represents the N<sub>2</sub>O VMR with an increase of 0.89 ppb/yr. The colored lines represent the N<sub>2</sub>O VMR according to RCP scenarios [<a href="#B60-atmosphere-14-00585" class="html-bibr">60</a>].</p>
Full article ">
29 pages, 2175 KiB  
Review
A Prompt Decarbonization Pathway for Shipping: Green Hydrogen, Ammonia, and Methanol Production and Utilization in Marine Engines
by Jie Shi, Yuanqing Zhu, Yongming Feng, Jun Yang and Chong Xia
Atmosphere 2023, 14(3), 584; https://doi.org/10.3390/atmos14030584 - 17 Mar 2023
Cited by 72 | Viewed by 17467
Abstract
The shipping industry has reached a higher level of maturity in terms of its knowledge and awareness of decarbonization challenges. Carbon-free or carbon-neutralized green fuel, such as green hydrogen, green ammonia, and green methanol, are being widely discussed. However, little attention has paid [...] Read more.
The shipping industry has reached a higher level of maturity in terms of its knowledge and awareness of decarbonization challenges. Carbon-free or carbon-neutralized green fuel, such as green hydrogen, green ammonia, and green methanol, are being widely discussed. However, little attention has paid to the green fuel pathway from renewable energy to shipping. This paper, therefore, provides a review of the production methods for green power (green hydrogen, green ammonia, and green methanol) and analyzes the potential of green fuel for application to shipping. The review shows that the potential production methods for green hydrogen, green ammonia, and green methanol for the shipping industry are (1) hydrogen production from seawater electrolysis using green power; (2) ammonia production from green hydrogen + Haber–Bosch process; and (3) methanol production from CO2 using green power. While the future of green fuel is bright, in the short term, the costs are expected to be higher than conventional fuel. Our recommendations are therefore as follows: improve green power production technology to reduce the production cost; develop electrochemical fuel production technology to increase the efficiency of green fuel production; and explore new technology. Strengthening the research and development of renewable energy and green fuel production technology and expanding fuel production capacity to ensure an adequate supply of low- and zero-emission marine fuel are important factors to achieve carbon reduction in shipping. Full article
(This article belongs to the Special Issue Shipping Emissions and Air Pollution)
Show Figures

Figure 1

Figure 1
<p>Planning system of hydropower [<a href="#B55-atmosphere-14-00584" class="html-bibr">55</a>].</p>
Full article ">Figure 2
<p>Potential onboard hydrogen storage technology.</p>
Full article ">Figure 3
<p>Offshore liquefied hydrogen production and ship refueling plant [<a href="#B145-atmosphere-14-00584" class="html-bibr">145</a>].</p>
Full article ">Figure 4
<p>Green ammonia production process using green hydrogen + HB process.</p>
Full article ">Figure 5
<p>HyMethShip project diagram.</p>
Full article ">Figure 6
<p>Production methods of green power, green hydrogen, green ammonia, and green methanol.</p>
Full article ">
15 pages, 3016 KiB  
Article
A Study on the Wide Range of Relative Humidity in Cirrus Clouds Using Large-Ensemble Parcel Model Simulations
by Miao Zhao and Xiangjun Shi
Atmosphere 2023, 14(3), 583; https://doi.org/10.3390/atmos14030583 - 17 Mar 2023
Cited by 2 | Viewed by 1862
Abstract
This study investigates the possible mechanisms related to the wide range of relative humidity in cirrus clouds (RHi). Under the closed adiabatic assumption, the impacts of vertical motion and ice crystal deposition/sublimation on RHi are investigated through in situ observations [...] Read more.
This study investigates the possible mechanisms related to the wide range of relative humidity in cirrus clouds (RHi). Under the closed adiabatic assumption, the impacts of vertical motion and ice crystal deposition/sublimation on RHi are investigated through in situ observations and parcel model simulations. Vertical motion is an active external force that changes the RHi, and ice crystal deposition/sublimation plays a role in mitigating the change in the RHi. They are the two most important mechanisms involved in controlling the RHi fluctuation during cirrus evolution and could well explain the wide range of RHi in wave-related cirrus clouds. Furthermore, a comparison of statistical cloud characteristics from both observations and simulations shows that a very low value (e.g., 0.001) for the water vapor ice deposition coefficient is highly unlikely. Full article
(This article belongs to the Special Issue Properties of Cirrus Cloud by Lidars: Observation and Theory)
Show Figures

Figure 1

Figure 1
<p>The statistical characteristics of the 50 prescribed W time series used to drive the parcel model (first row). They are the occurrence frequencies of W and ΔT (<b>left</b>) and the spectrum analysis of the W time series (<b>right</b>). The second and third rows present two cases (Wave1 and Wave2, respectively). The left columns show the time series of W (gray solid line), the ΔH (black solid line), and ΔT (red dotted line). The right columns show the spectrum analyses of the corresponding W time series. The ΔT and ΔH were calculated using the parcel model, with the initial conditions of T = 225 K and P = 250 hPa.</p>
Full article ">Figure 2
<p>Cirrus evolutions from parcel model simulations using the W time series shown in <a href="#atmosphere-14-00583-f001" class="html-fig">Figure 1</a> (Wave1 or Wave2). Shown are RH<sub>i</sub> (blue thick solid line), W (gray thin solid line), N<sub>i</sub> (orange thin solid line), q<sub>i</sub> (green thin solid line). Black dotted line represents the line of W = 0, also the line of RH<sub>i</sub> = 100%. The shaded figures indicate the N<sub>i</sub> contribution from each radius bin. A total of 50 bins were used in this study. Experiment names and the corresponding W time series are marked in the upper left corners. Note that the ranges of the N<sub>i</sub>-axis and RH<sub>i</sub>-axis from the ICadL experiment are different from the others.</p>
Full article ">Figure 3
<p>The scatter plots of W, RH<sub>i</sub>, and N<sub>i</sub>R<sub>i</sub> from observations (<b>upper panel</b>) and the REF experiment (<b>lower panel</b>). The dotted lines indicate RH<sub>i</sub> = 100% or W = 0 m s<sup>−1</sup>. Note that the 1000 samples selected for parcel model simulations are marked by red dots.</p>
Full article ">Figure 4
<p>The occurrence frequencies of N<sub>i</sub>R<sub>i</sub> and RH<sub>i</sub> from observations (black lines) and the REF experiment (red lines). The dotted line indicates RH<sub>i</sub> = 100%.</p>
Full article ">Figure 5
<p>The occurrence frequency of RH<sub>i</sub> from the REF, Wamp, and Wfre experiments (<b>upper panel</b>). The lower panel shows the spectrum analysis of the W time series (<b>left</b>), the distributions of W and ΔT (<b>middle</b>), and the scatter plots of RH<sub>i</sub> vs. W (<b>right</b>) from the Wamp (<b>second row</b>) and Wfre (<b>third row</b>) experiments.</p>
Full article ">Figure 6
<p>The occurrence frequency of RH<sub>i</sub> from the Ref, ICnosub, ICadH, and ICadL experiments (<b>upper panel</b>). The lower panel shows scatter plots similar to those in <a href="#atmosphere-14-00583-f004" class="html-fig">Figure 4</a> but for the ICnosub (<b>second row</b>), ICadH (<b>third row</b>), and ICadL (<b>fourth row</b>) experiments.</p>
Full article ">
15 pages, 8129 KiB  
Article
Application of Radar Radial Velocity Data Assimilation in the Forecasts of Typhoon Linfa Based on Different Horizontal Length Scale Factors
by Huimin Bian, Jinzhong Min and Feifei Shen
Atmosphere 2023, 14(3), 582; https://doi.org/10.3390/atmos14030582 - 17 Mar 2023
Cited by 1 | Viewed by 2004
Abstract
In order to explore the improvement of radar radial velocity data assimilation on the initial and forecast fields of typhoons, this study assimilates the quality-controlled radial velocity data in the case of Typhoon Linfa (2015) using the three-dimensional variational data assimilation system of [...] Read more.
In order to explore the improvement of radar radial velocity data assimilation on the initial and forecast fields of typhoons, this study assimilates the quality-controlled radial velocity data in the case of Typhoon Linfa (2015) using the three-dimensional variational data assimilation system of the weather research and forecasting model (WRF-3DVAR), and then conducts several sensitivity experiments with different horizontal length scale factors. The results show that reducing the horizontal length scale factor of the background error covariance can effectively assimilate the micro- and meso-scale information from radar data and improve the forecasting effect of Linfa. Following the optimization of the horizontal length scale factor, the radial velocity data assimilation can improve the typhoon wind field structure, produce reasonable cyclonic wind field increments, and further improve the dynamic and thermal structure of the inner core area of the typhoon. Then, we can obtain a better initial field of model forecasting, and thus typhoon track and intensity forecasting are improved. Full article
(This article belongs to the Special Issue Data Assimilation for Predicting Hurricane, Typhoon and Storm)
Show Figures

Figure 1

Figure 1
<p>Comparison of the radial velocity (unit: m s<sup>−1</sup>) from the Shantou Doppler radar (STRD) (<b>a</b>) before and (<b>b</b>) after the quality control at 1800 UTC on 8 July 2015.</p>
Full article ">Figure 2
<p>The weather research and forecasting (WRF) model domains with the China Meteorological Administration’s (CMA) best tracking of Typhoon Linfa from 0000 UTC on July 5 to 1800 UTC on 9 July 2015. The location of the STRD is indicated by asterisks. The maximum detection range of the STRD is circled at 230 km. Tropical cyclones are distinguished by their color: tropical depression (TD), tropical storm (TS), severe tropical storm (STS), typhoon (TY), and severe typhoon (STY). The shading indicates the terrain height (unit: m).</p>
Full article ">Figure 3
<p>Flow chart of the experiments.</p>
Full article ">Figure 4
<p>500 hPa geopotential height (solid blue line, unit: dagpm), wind field (wind bar, long line = 4 m s<sup>−1</sup>, short line = 2 m s<sup>−1</sup>, wind triangle = 20 m s<sup>−1</sup>) and isotherm (red line, unit: °C) at (<b>a</b>) 1200 UTC on 8 July and (<b>b</b>) 0000 UTC on 9 July 2015.</p>
Full article ">Figure 5
<p>700 hPa wind speed (shaded, unit: m s<sup>−1</sup>) and wind vector (arrow, reference vector = 20 m s<sup>−1</sup>) of the (<b>a</b>) CTRL, (<b>b</b>) DA_len1.0 and (<b>c</b>) DA_len0.25 experiments at 1800 UTC on 8 July 2015.</p>
Full article ">Figure 6
<p>The analysis increments of the wind speed (shaded, unit: m s<sup>−1</sup>) and wind vector (arrow, reference vector = 20 m s<sup>−1</sup>) at 700 hPa of (<b>a</b>) DA_len1.0 and (<b>b</b>) DA_len0.25 at 1800 UTC on 8 July 2015. The observed typhoon center (red asterisk) is also indicated.</p>
Full article ">Figure 7
<p>Root mean square errors (unit: m s<sup>−1</sup>) of the radial wind speed between DA_len1.0 (red), DA_len0.25 (blue) and the observation during the cyclic assimilation.</p>
Full article ">Figure 8
<p>Sea-level pressure field (black contour), 10-m wind speed (shaded, unit: m s<sup>−1</sup>) and wind vector (arrow, reference vector = 20 m s<sup>−1</sup>) of the (<b>a</b>) CTRL, (<b>b</b>) DA_len1.0 and (<b>c</b>) DA_len0.25 experiments at 0000 UTC on 9 July 2015.</p>
Full article ">Figure 9
<p>Vertical profiles of the horizontal wind speed (shaded, unit: m s<sup>−1</sup>) and potential temperature (black contours, unit: K) in the (<b>a</b>) CTRL, (<b>b</b>) DA_len1.0 and (<b>c</b>) DA_len0.25 experiments at 0000 UTC on 9 July 2015.</p>
Full article ">Figure 10
<p>Azimuthally averaged tangential wind speed (unit: m s<sup>−1</sup>) of the (<b>a</b>) CTRL, (<b>b</b>) DA_len1.0 and (<b>c</b>) DA_len0.25 experiments at 0000 UTC on 9 July 2015.</p>
Full article ">Figure 11
<p>(<b>a</b>) Forecast tracks in the three experiments and the CMA’s best track; (<b>b</b>) track errors (unit: km) of the three experiments compared to the CMA’s best track during the forecast period.</p>
Full article ">Figure 12
<p>(<b>a</b>) Minimum sea-level pressure (unit: mb) and (<b>b</b>) maximum surface wind speed (unit: m s<sup>−1</sup>) given by the three sets of experiments and the CMA observations.</p>
Full article ">
7 pages, 1503 KiB  
Communication
Potential Distribution and Priority Conservation Areas of Pseudotsuga sinensis Forests under Climate Change in Guizhou Province, Southwesten China
by Wangjun Li, Yingqian Yu, Tu Feng, Bin He, Xiaolong Bai and Shun Zou
Atmosphere 2023, 14(3), 581; https://doi.org/10.3390/atmos14030581 - 17 Mar 2023
Cited by 1 | Viewed by 1451
Abstract
Priority conservation areas are the key areas of biodiversity maintenance and ecosystem conservation. Based on a Maxent model, this study predicted the potential distribution of Pseudotsuga sinensis under the current climate and future climate change scenarios in Guizhou province, and then, assessed three [...] Read more.
Priority conservation areas are the key areas of biodiversity maintenance and ecosystem conservation. Based on a Maxent model, this study predicted the potential distribution of Pseudotsuga sinensis under the current climate and future climate change scenarios in Guizhou province, and then, assessed three kinds of priority conservation area under climate change. The results were as follows: (1) The AUC (Area Under the Curve) values showed excellent prediction accuracy of the model. (2) The areas of the potential habitats of P. sinensis forests under the current climate and future climate change scenarios were 22,062.85 km2 and 18,395.92 km2, respectively. As for their spatial distribution, the potential habitats of P. sinensis forests were distributed in the Bijie, Zunyi, Tongren, Liupanshui and Xingyi regions under the current climate, and in the Kaili region, in addition to the above-mentioned cities, under future climate change scenarios. (3) The total area of priority conservation areas under climate change was 25,350.26 km2. The area of the predicted sustainable potential habitats was 15,075.96 km2, of the vulnerable potential habitats was 7256.59 km2 and of the derivative potential habitats was 3017.71 km2. Full article
(This article belongs to the Special Issue Vegetation and Climate Relationships (2nd Edition))
Show Figures

Figure 1

Figure 1
<p>Location and administrative division of Guizhou province.</p>
Full article ">Figure 2
<p>(<b>A</b>) Potential habitats under current climate; (<b>B</b>) Potential habitats under future climate change scenarios.</p>
Full article ">Figure 3
<p>Priority conservation areas under climate change.</p>
Full article ">
19 pages, 10748 KiB  
Article
Impact of Hyperspectral Infrared Sounding Observation and Principal-Component-Score Assimilation on the Accuracy of High-Impact Weather Prediction
by Qi Zhang and Min Shao
Atmosphere 2023, 14(3), 580; https://doi.org/10.3390/atmos14030580 - 17 Mar 2023
Cited by 1 | Viewed by 1969
Abstract
Observations from a hyperspectral infrared (IR) sounding interferometer such as the Infrared Atmospheric Sounding Interferometer (IASI) and the Cross-Track Infrared Sounder (CrIS) are crucial to numerical weather prediction (NWP). By measuring radiance at the top of the atmosphere using thousands of channels, these [...] Read more.
Observations from a hyperspectral infrared (IR) sounding interferometer such as the Infrared Atmospheric Sounding Interferometer (IASI) and the Cross-Track Infrared Sounder (CrIS) are crucial to numerical weather prediction (NWP). By measuring radiance at the top of the atmosphere using thousands of channels, these observations convey accurate atmospheric information to the initial condition through data assimilation (DA) schemes. The massive data volume has pushed the community to develop novel approaches to reduce the number of assimilated channels while retaining as much information content as possible. Thus, channel-selection schemes have become widely accepted in every NWP center. Two significant limitations of channel-selection schemes are (1) the deficiency in retaining the observational information content and (2) the higher cross-channel correlation in the observational error (R) matrix. This paper introduces a hyperspectral IR observation DA scheme in the principal component (PC) space. Four-month performance comparison case studies using the Weather Research and Forecasting model (WRF) as a forecast module between PC-score assimilation and the selected-channel assimilation experiment show that the PC-score assimilation scheme can reduce the initial condition’s root-mean-squared error for temperature and water vapor compared to the channel-selection scheme and thus improve the forecasting of precipitation and high-impact weather. Case studies using the Unified Forecast System Short-Range Weather (UFS-SRW) application as forecast module also indicate that the positive impact can be retained among different NWP models. Full article
Show Figures

Figure 1

Figure 1
<p>Observational error covariance in PC space (<b>a</b>,<b>c</b>) and channel space (<b>b</b>,<b>d</b>) for CrIS (<b>a</b>,<b>b</b>) and IASI (<b>c</b>,<b>d</b>).</p>
Full article ">Figure 2
<p>Spatial coverage of the forecast domain.</p>
Full article ">Figure 3
<p>Forecast system workflow.</p>
Full article ">Figure 4
<p>Temperature (<b>a</b>) and specific humidity (<b>b</b>) RMSE departure of PC-score assimilation (red line) and selected-channel radiance assimilation (green line).</p>
Full article ">Figure 5
<p>CSI departure for precipitation within 2.5 and 7.5 mm/h (<b>a</b>–<b>c</b>) and above 7.5 mm/h (<b>d</b>–<b>f</b>) from PC-score assimilation (red line) and selected-channel radiance assimilation (green line).</p>
Full article ">Figure 6
<p>PC-score assimilation (red line) and selected-channel radiance assimilation (green line) STP departure for the 2019Mar03 (<b>a</b>), 2020Mar03 (<b>b</b>), and 2020Apr12 (<b>c</b>) cases.</p>
Full article ">Figure 7
<p>The temperature (<b>a</b>,<b>e</b>), specific humidity (<b>b</b>,<b>f</b>), u-component wind (<b>c</b>,<b>g</b>), and v-component (<b>d</b>,<b>h</b>) wind RMSE profiles from the PC-score assimilation experiment and CTL analysis (<b>a</b>–<b>d</b>) and 12 h lead-time forecast (<b>e</b>–<b>h</b>).</p>
Full article ">Figure 8
<p>The ACC profiles of temperature (<b>a</b>), specific humidity (<b>b</b>), u-component wind (<b>c</b>), and v-component wind (<b>d</b>) from the PC-score assimilation system (red) and CTL (black).</p>
Full article ">Figure 9
<p>Hanssen and Kuipers discriminant (<b>a</b>), multi-category Hanssen and Kuipers discriminant (<b>b</b>), and Kling–Gupta efficiency (<b>c</b>) time series from the PC-score assimilation system (red) and CTL (black).</p>
Full article ">Figure 10
<p>Time series of the fixed-layer significant tornado parameter (<b>a</b>) and its POD (<b>b</b>) derived from the PC-score assimilation system forecast result (red) and CTL (black).</p>
Full article ">Figure 11
<p>The contribution of each variable in the calculation of the significant tornado parameter, with the red (black) line representing the PC-score experiment forecast result (CTL).</p>
Full article ">Figure 12
<p>Time series of the EHI (<b>a</b>) and its POD (<b>b</b>) derived from the PC-score assimilation experiment forecast result (red) and CTL (black).</p>
Full article ">Figure 13
<p>Comparison between the HRRR and pseudo-operational forecasts initialized at 00:00 UTC 7 September (<b>a</b>,<b>b</b>) and 00:00 UTC 7 September (<b>c</b>,<b>d</b>). The black and gray marks (with different shapes representing the outbreak time) are the hail outbreak locations.</p>
Full article ">Figure 14
<p>Comparison between the HRRR and pseudo-operational forecasts initialized at 12:00 UTC 15 December (<b>a</b>,<b>b</b>) and 00:00 UTC 16 December 00:00 UTC (<b>c</b>,<b>d</b>). The black and gray marks (with the different shapes representing the outbreak time) are the hail outbreak locations.</p>
Full article ">Figure A1
<p>Correlation coefficient matrix in principle component (<b>a</b>,<b>c</b>) and radiance (<b>b</b>,<b>d</b>) space for CrIS (<b>a</b>,<b>b</b>) and IASI (<b>c</b>,<b>d</b>).</p>
Full article ">Figure A2
<p>Weighting function for CrIS and IASI in PC and radiance space. Profile dataset comes from European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) 60-level sample profile dataset from the Monitoring Atmospheric Composition and Climate (MACC) project (available at <a href="https://nwp-saf.eumetsat.int/site/download/profile_datasets/60l_macc.dat.tar.bz2" target="_blank">https://nwp-saf.eumetsat.int/site/download/profile_datasets/60l_macc.dat.tar.bz2</a> (accessed on 3 February 2023)).</p>
Full article ">
14 pages, 2136 KiB  
Article
Emission Source Areas of Fine Particulate Matter (PM2.5) in Ho Chi Minh City, Vietnam
by Tuyet Nam Thi Nguyen, Nguyen Xuan Du and Nguyen Thi Hoa
Atmosphere 2023, 14(3), 579; https://doi.org/10.3390/atmos14030579 - 17 Mar 2023
Cited by 8 | Viewed by 4663
Abstract
This study aims to determine emission source areas of fine particulate matter (PM2.5) in Ho Chi Minh (HCM) City, Vietnam, using a conditional bivariate probability function (CBPF) and hybrid receptor models, including three-dimensional potential source contribution function (3D-PSCF) and concentration-weighted trajectory [...] Read more.
This study aims to determine emission source areas of fine particulate matter (PM2.5) in Ho Chi Minh (HCM) City, Vietnam, using a conditional bivariate probability function (CBPF) and hybrid receptor models, including three-dimensional potential source contribution function (3D-PSCF) and concentration-weighted trajectory (3D-CWT), considering latitudes, longitudes, and height of trajectory segments. Uncertainties of the CBPF and 3D-PSCF/3D-CWT were evaluated based on the 95th confidence intervals and 95% confidence levels, respectively. For the local scale, PM2.5 in HCM City was primarily emitted from shallow or common ground sources (e.g., vehicle emissions) throughout the year. Regarding non-local source areas, PM2.5 in HCM City is contributed by those originated from the East Sea (e.g., shipping emissions) and southeastern Vietnam (e.g., Binh Duong and Dong Nai provinces) having several industrial zones with PM2.5 emission sources, especially in the dry season (December to April of the following year). In the rainy season (May–November), PM2.5 derived from emission sources in the Mekong Delta (e.g., biomass burning) might be transported to HCM City. However, contribution of the non-local sources to PM2.5 pollution in HCM City during the rainy season is less important because of PM2.5 deposition stemmed from the high rainfall amount in this season. Full article
(This article belongs to the Special Issue Air Quality in Metropolitan Areas and Megacities)
Show Figures

Figure 1

Figure 1
<p>Location of the PM<sub>2.5</sub> monitoring site and the meteorological stations in Ho Chi Minh City of Vietnam.</p>
Full article ">Figure 2
<p>Mean of PM<sub>2.5</sub> mass concentrations in HCM City shown in monthly and seasonal variation, diurnal variation in the dry season, and diurnal variation in the rainy season. The whiskers represent standard deviations of the mean.</p>
Full article ">Figure 3
<p>CBPF plots of PM<sub>2.5</sub> in the dry and rainy seasons of HCM City during the study period.</p>
Full article ">Figure 4
<p>Bootstrapped 3D-PSCF and 3D-CWT results for PM<sub>2.5</sub> at the 95% upper confidence levels in the dry and rainy seasons. The yellow circles represent the receptor site in HCM City.</p>
Full article ">Figure 5
<p>Trajectory cluster analysis plots of PM<sub>2.5</sub> in the dry and rainy seasons of HCM City. The numbers in each circle indicate cluster number. The percentages represent the contribution of each trajectory cluster.</p>
Full article ">
22 pages, 3364 KiB  
Article
Assessment of Air Pollution in Different Areas (Urban, Suburban, and Rural) in Slovenia from 2017 to 2021
by Maja Ivanovski, Kris Alatič, Danijela Urbancl, Marjana Simonič, Darko Goričanec and Rudi Vončina
Atmosphere 2023, 14(3), 578; https://doi.org/10.3390/atmos14030578 - 17 Mar 2023
Cited by 5 | Viewed by 4016
Abstract
Air pollution can have a significant effect on human health. The present work is focused on the investigation of daily, monthly, and annual concentration levels of five typical air pollutants (SO2, NO2, NOX, PM10, and [...] Read more.
Air pollution can have a significant effect on human health. The present work is focused on the investigation of daily, monthly, and annual concentration levels of five typical air pollutants (SO2, NO2, NOX, PM10, and PM2.5) in the Republic of Slovenia (RS) from January 2017 to December 2021. The study was conducted at five different monitoring stations of the following kind: traffic (A), industrial (D), and background (B, C, E). The obtained results showed a decline in the average concentrations for all the studied air pollutants through the years, respectively. The daily average SO2 concentrations were the lowest in the year 2021 at location B, which is classified as background location, while the highest were detected in the year 2018 at location E, which is also classified as background location. The average daily concentrations of NO2 and NOX were the highest at location A in the year 2017, whereas the lowest were detected in the year 2010 and 2021. It is believed that those results are a consequence of measures set by the Slovenian government during the COVID-19 pandemic. The PM10 and PM2.5 daily average concentrations were the highest at location A in 2017, while the lowest were observed in the year 2019 at location C. Meteorological parameters (temperature, wind speed, and relative humidity) were studied in addition. In general, the high temperatures in ambient air are responsible for the intense concentrations of air pollutants. It was found in the study results for temperature, wind speed, and relative humidity that no significant difference was shown between studied years. Full article
Show Figures

Figure 1

Figure 1
<p>The geographical location of the study area with marked monitoring stations (A—traffic, B—background, C—background, D—industrial, and E—background).</p>
Full article ">Figure 2
<p>Average monthly SO<sub>2</sub> concentrations at locations A, B, D, and E between 2017 and 2021 (from left to right, from top to bottom).</p>
Full article ">Figure 2 Cont.
<p>Average monthly SO<sub>2</sub> concentrations at locations A, B, D, and E between 2017 and 2021 (from left to right, from top to bottom).</p>
Full article ">Figure 3
<p>Average yearly SO<sub>2</sub> concentrations at locations A, B, D, and E between 2017 and 2021.</p>
Full article ">Figure 4
<p>Average monthly NO<sub>2</sub> concentrations at locations A, B, D, and E between 2017 and 2021 (from left to right, from top to bottom).</p>
Full article ">Figure 5
<p>Average yearly NO<sub>2</sub> concentrations at locations A, B, D, and E between 2017 and 2021.</p>
Full article ">Figure 6
<p>Average monthly NO<sub>X</sub> concentrations at locations A, B, D, and E between 2017 and 2021 (from left to right, from top to bottom).</p>
Full article ">Figure 7
<p>Average annual NO<sub>X</sub> concentrations at locations A, B, D, and E between 2017 and 2021.</p>
Full article ">Figure 8
<p>Average monthly PM<sub>10</sub> concentrations at locations A, B, C, and D between 2017 and 2021 (from left to right, from top to bottom).</p>
Full article ">Figure 9
<p>Average annual PM<sub>10</sub> concentrations at locations A, B, D, and E between 2017 and 2021.</p>
Full article ">Figure 10
<p>Average monthly PM<sub>2.5</sub> concentrations at locations A and D between 2017 and 2021 (from left to right).</p>
Full article ">Figure 11
<p>Average annual PM<sub>2.5</sub> concentrations at locations A and D between 2017 and 2021.</p>
Full article ">Figure 12
<p>Wind roses for each studied location and year.</p>
Full article ">Figure 12 Cont.
<p>Wind roses for each studied location and year.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop