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

Atmosphere, Volume 12, Issue 5 (May 2021) – 130 articles

Cover Story (view full-size image): Aerosols and clouds play critical roles in Earth’s climate system, air quality, and hydrological cycle. Lidar measurements provide essential vertical profiles of aerosols and clouds. This study is the first ever to train and employ a Convolutional Neural Network (CNN) for the detection of aerosols and clouds using space-borne lidar data, improving the horizontal resolution needed to detect atmospheric features from the Cloud-Aerosol Transport System (CATS) by a factor of 12 and enabling more accurate cloud-aerosol discrimination.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:
19 pages, 2383 KiB  
Article
Potential Human and Plant Pathogenic Species in Airborne PM10 Samples and Relationships with Chemical Components and Meteorological Parameters
by Salvatore Romano, Mattia Fragola, Pietro Alifano, Maria Rita Perrone and Adelfia Talà
Atmosphere 2021, 12(5), 654; https://doi.org/10.3390/atmos12050654 - 20 May 2021
Cited by 6 | Viewed by 3597
Abstract
A preliminary local database of potential (opportunistic) airborne human and plant pathogenic and non-pathogenic species detected in PM10 samples collected in winter and spring is provided, in addition to their seasonal dependence and relationships with meteorological parameters and PM10 chemical species. The PM10 [...] Read more.
A preliminary local database of potential (opportunistic) airborne human and plant pathogenic and non-pathogenic species detected in PM10 samples collected in winter and spring is provided, in addition to their seasonal dependence and relationships with meteorological parameters and PM10 chemical species. The PM10 samples, collected at a Central Mediterranean coastal site, were analyzed by the 16S rRNA gene metabarcoding approach, and Spearman correlation coefficients and redundancy discriminant analysis tri-plots were used to investigate the main relationships. The screening of 1187 detected species allowed for the detection of 76 and 27 potential (opportunistic) human and plant pathogens, respectively. The bacterial structure of both pathogenic and non-pathogenic species varied from winter to spring and, consequently, the inter-species relationships among potential human pathogens, plant pathogens, and non-pathogenic species varied from winter to spring. Few non-pathogenic species and even fewer potential human pathogens were significantly correlated with meteorological parameters, according to the Spearman correlation coefficients. Conversely, several potential plant pathogens were strongly and positively correlated with temperature and wind speed and direction both in winter and in spring. The number of strong relationships between presumptive (human and plant) pathogens and non-pathogens, and meteorological parameters slightly increased from winter to spring. The sample chemical composition also varied from winter to spring. Some potential human and plant pathogens were correlated with chemicals mainly associated with marine aerosol and/or with soil dust, likely because terrestrial and aquatic environments were the main habitats of the detected bacterial species. The carrier role on the species seasonal variability was also investigated. Full article
(This article belongs to the Special Issue Bioaerosols: Composition, Meteorological Impact, and Transport)
Show Figures

Figure 1

Figure 1
<p>Relative percentage contribution of the overall amount of the 27 potential plant pathogenic species, 76 potential human pathogenic species, 1084 potential non-pathogenic species, and unclassified species in (<b>a</b>) winter (S1–S10) and (<b>b</b>) spring (S11–S20) samples. Bray–Curtis dissimilarity dendrograms highlighting the relatedness of the species-level bacterial communities in (<b>a</b>) winter and (<b>b</b>) spring samples have also been reported.</p>
Full article ">Figure 2
<p>Mean percentage contribution (on a logarithmic scale) of the 10 most abundant potential and pervasive human pathogenic species in (<b>a</b>) winter and (<b>b</b>) spring: <span class="html-italic">Sphingobacterium multivorum</span> (<span class="html-italic">S. multivorum)</span>, <span class="html-italic">Johnsonella ignava</span> (<span class="html-italic">J. ignava</span>), <span class="html-italic">Streptococcus bovis</span> (<span class="html-italic">S. bovis</span>), <span class="html-italic">Staphylococcus aureus (S. aureus)</span>, <span class="html-italic">Clostridium cadaveris</span> (<span class="html-italic">C. cadaveris</span>), <span class="html-italic">Peptococcus niger</span> (<span class="html-italic">P. niger</span>), <span class="html-italic">Propionibacterium avidum</span> (<span class="html-italic">P. avidum</span>)<span class="html-italic">, Propionibacterium acnes</span> (<span class="html-italic">P. acnes</span>), <span class="html-italic">Providencia rettgeri</span> (<span class="html-italic">P. rettgeri</span>)<span class="html-italic">, Acinetobacter lwoffii</span> (<span class="html-italic">A. lwoffii</span>)<span class="html-italic">, Acinetobacter ursingii</span> (<span class="html-italic">A. ursingii</span>)<span class="html-italic">, Acinetobacter johnsonii</span> (<span class="html-italic">A. johnsonii</span>), <span class="html-italic">Enterobacter aerogenes</span> (<span class="html-italic">E. aerogenes</span>)<span class="html-italic">, Enterobacter amnigenus</span> (<span class="html-italic">E. amnigenus</span>), and <span class="html-italic">Enterobacter hormaechei</span> (<span class="html-italic">E. hormaechei</span>). Error bars represent the standard error of the mean. Phyla related to each species are also reported on the left (Prot.: Proteobacteria, Actin.: Actinobacteria, Firm.: Firmicutes, Bact.: Bacteroidetes).</p>
Full article ">Figure 3
<p>Mean percentage contribution (on a logarithmic scale) of the most abundant potential and pervasive plant pathogenic species in (<b>a</b>) winter and (<b>b</b>) spring: <span class="html-italic">Bacillus megaterium</span> (<span class="html-italic">B. megaterium</span>)<span class="html-italic">, Rathayibacter tritici</span> (<span class="html-italic">R. tritici</span>)<span class="html-italic">, Clavibacter michiganensis</span> (<span class="html-italic">C. michiganensis</span>)<span class="html-italic">, Curtobacterium flaccumfaciens</span> (<span class="html-italic">C. flaccumfaciens</span>)<span class="html-italic">, Agrobacterium larrymoorei</span> (<span class="html-italic">A. larrymoorei</span>)<span class="html-italic">, Erwinia mallotivora</span> (<span class="html-italic">E. mallotivora</span>)<span class="html-italic">, Janthinobacterium agaricidamnosum</span> (<span class="html-italic">J. agaricidamnosum</span>)<span class="html-italic">, Sphingomonas melonis</span> (<span class="html-italic">S. melonis</span>)<span class="html-italic">, Pseudomonas viridiflava</span> (<span class="html-italic">P. viridiflava</span>), and <span class="html-italic">Enterobacter cloacae</span> (<span class="html-italic">E. cloacae</span>). Error bars represent the standard error of the mean. Phyla related to each species are also reported on the left (Prot.: Proteobacteria, Actin.: Actinobacteria, Firm.: Firmicutes).</p>
Full article ">
5 pages, 203 KiB  
Editorial
Climate Change, Aquatic Ecosystems and Human Infectious Diseases in a Globalised World
by Arturo Sousa, Mónica Aguilar-Alba and Leoncio García-Barrón
Atmosphere 2021, 12(5), 653; https://doi.org/10.3390/atmos12050653 - 20 May 2021
Cited by 1 | Viewed by 2263
Abstract
One of the greatest challenges that human society currently faces is the hazard of climate change with respect to human health [...] Full article
20 pages, 6499 KiB  
Article
Dry Spells in Croatia: Observed Climate Change and Climate Projections
by Ivana Marinović, Ksenija Cindrić Kalin, Ivan Güttler and Zoran Pasarić
Atmosphere 2021, 12(5), 652; https://doi.org/10.3390/atmos12050652 - 20 May 2021
Cited by 19 | Viewed by 4003
Abstract
This study performs a systematic analysis of the recent and future changes of dry spells (DS) in Croatia. DS are defined as consecutive sequences of days with daily precipitation less than 5 mm of the precipitation-per-day threshold (DS5). Daily precipitation data come from [...] Read more.
This study performs a systematic analysis of the recent and future changes of dry spells (DS) in Croatia. DS are defined as consecutive sequences of days with daily precipitation less than 5 mm of the precipitation-per-day threshold (DS5). Daily precipitation data come from a dense national rain gauge network (covering seven regions) and span the period 1961–2015. The spatial and temporal changes of the observed mean (MDS5) and maximum (MxDS5) seasonal and annual dry spells were analysed by means of the Kendall tau method and the partial trend method. Future changes of DS5 were assessed by employing the three regional climate models (RegCM4, CCLM4 and RCA4) covering the EURO-CORDEX domain with a 12.5-km horizontal resolution, resulting in a realistic orography and land–sea border over Croatia. The models were forced at their boundaries by the four CMIP5 global climate models. For the reference period 1971–2000, the observed, as well as modelled, DS5 were analysed, and the systematic model errors were assessed. Finally, the projections and future changes of the DS5 statistics based on simulations under the high and medium greenhouse gases concentration scenarios (i.e., RCP8.5 and RCP4.5) with a focus on the climate change signal between 1971–2000 and two future periods, 2011–2040 and 2041–2070, were examined. A prevailing increasing trend of MDS5 was found in the warm part of the year, being significant in the mountainous littoral and North Adriatic coastal region. An increasing trend of MxDS5 was also found in the warm part of the year (both the spring and summer), and it was particularly pronounced along the Adriatic coast, while a coherent negative trend pattern was found in the autumn. By applying the partial trend methodology, an increase was found in the very long DS5 (above the 90th percentile) in the recent half of the analysed 55-year period in all seasons, except in the autumn when shortening in the DS5 was detected. The climate change signal during the two analysed future periods was positive for the summer in all regions, weakly negative for the winter and not conclusive for the spring, autumn and year. It was found that no RCM-GCM combination is the best in all cases, since the most successful model combinations depend on the season and location. Full article
Show Figures

Figure 1

Figure 1
<p>Locations of the meteorological stations used in this study and the regions listed in <a href="#atmosphere-12-00652-t001" class="html-table">Table 1</a>. The nine stations selected for the regional climate model validations and climate change projection analysis are highlighted with black diamonds.</p>
Full article ">Figure 2
<p>Seasonal (upper four panels) and annual (bottom panel) mean dry spell durations in days (brown bars) for the 5-mm precipitation-per-day thresholds. Red dots indicate the corresponding median values, while the interquartile range (25th and 75th percentiles) for each station is depicted in grey. The seven regions are distinguished by bold dashed vertical lines.</p>
Full article ">Figure 3
<p>Seasonal (upper four panels) and annual (bottom panel) maximum dry spell durations (MxDS5, in days) for the 5-mm precipitation-per-day thresholds. The seven regions are distinguished by bold dashed vertical lines.</p>
Full article ">Figure 4
<p>The trends in the seasonal and annual MDS5 (in days/10 years). Dots refer to trend values (Sen’s slope) for each station in the associated region. Blue dots present nonsignificant trends, whereas significant trends are highlighted by red dots.</p>
Full article ">Figure 5
<p>The same as <a href="#atmosphere-12-00652-f004" class="html-fig">Figure 4</a> but for MxDS5.</p>
Full article ">Figure 6
<p>The partial trend scatter diagram for DS5 for nine selected stations in Croatia. Seasons are distinguished by different colours. The 95th percentile values from each period (1961–1987 and 1988–2014) are denoted by circles.</p>
Full article ">Figure 7
<p>Average MDS5 from observations (O) and the median of 30-year average MDS5 from the models (M) for seven stations on the annual and seasonal time scale. Seasons are distinguished by different colours.</p>
Full article ">Figure 8
<p>The median difference in MDS5 duration between two future periods and historical period for RCP4.5 and RCP8.5 scenarios, for seven stations. The circles refer to the period P1 (2011–2040) vs. P0 (1971–2000) and triangles refer to the second period P2 (2041–2070) vs. P0. Seasons are distinguished by different colours. For each station, left symbols present results for RCP4.5, while right symbols present results for RCP8.5 scenario. The corresponding differences in MDS5 (in days) are given on the vertical axis.</p>
Full article ">
17 pages, 4597 KiB  
Article
A Short-Term Wind Speed Forecasting Model Based on a Multi-Variable Long Short-Term Memory Network
by Anqi Xie, Hao Yang, Jing Chen, Li Sheng and Qian Zhang
Atmosphere 2021, 12(5), 651; https://doi.org/10.3390/atmos12050651 - 19 May 2021
Cited by 37 | Viewed by 4595
Abstract
Accurately forecasting wind speed on a short-term scale has become essential in the field of wind power energy. In this paper, a multi-variable long short-term memory network model (MV-LSTM) based on Pearson correlation coefficient feature selection is proposed to predict the short-term wind [...] Read more.
Accurately forecasting wind speed on a short-term scale has become essential in the field of wind power energy. In this paper, a multi-variable long short-term memory network model (MV-LSTM) based on Pearson correlation coefficient feature selection is proposed to predict the short-term wind speed. The proposed method utilizes multiple historical meteorological variables, such as wind speed, temperature, humidity, and air pressure, to predict the wind speed in the next hour. Hourly data collected from two ground observation stations in Yanqing and Zhaitang in Beijing were divided into training and test sets. The training sets were used to train the model, and the test sets were used to evaluate the model with the root-mean-square error (RMSE), mean absolute error (MAE), mean bias error (MBE), and mean absolute percentage error (MAPE) metrics. The proposed method is compared with two other forecasting methods (the autoregressive moving average model (ARMA) method and the single-variable long short-term memory network (LSTM) method, which inputs only historical wind speed data) based on the same dataset. The experimental results prove the feasibility of the MV-LSTM method for short-term wind speed forecasting and its superiority to the ARMA method and the single-variable LSTM method. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location of the stations in Beijing, Northeast China.</p>
Full article ">Figure 2
<p>Real meteorological data from Yanqing station.</p>
Full article ">Figure 3
<p>Real meteorological data from Zhaitang station.</p>
Full article ">Figure 4
<p>LSTM cell structure.</p>
Full article ">Figure 5
<p>The framework of the proposed multi-variable LSTM network method for wind speed forecasting.</p>
Full article ">Figure 6
<p>The wind speed forecasting results obtained by the different models on the dataset of Yanqing Station. In the top sub-figure, the x-axis is the wind speed value (unit: m/s), and the y-axis is the time of the data points (time interval: 1 h); in the bottom sub-figure, the x-axis is the residual errors of the MV-LSTM method (unit: m/s), and the y-axis is the same as it in the top sub-figure.</p>
Full article ">Figure 7
<p>The wind speed forecasting results obtained by the different models on the dataset of Zhaitang Station. In the top sub-figure, the x-axis is the wind speed value (unit: m/s), and the y-axis is the time of the data points (time interval: 1 h); in the bottom sub-figure, the x-axis is the residual errors of the MV-LSTM method (unit: m/s), and the y-axis is the same as it in the top sub-figure.</p>
Full article ">Figure 8
<p>Boxplot of the residual errors of each method with different ranges of the wind speed on the test set of Yanqing Station.</p>
Full article ">Figure 9
<p>Boxplot of the residual errors of each method with different ranges of the wind speed on the test set of Zhaitang Station.</p>
Full article ">
36 pages, 57857 KiB  
Review
Recent Advances in Our Understanding of Tropical Cyclone Intensity Change Processes from Airborne Observations
by Robert F. Rogers
Atmosphere 2021, 12(5), 650; https://doi.org/10.3390/atmos12050650 - 19 May 2021
Cited by 19 | Viewed by 6294
Abstract
Recent (past ~15 years) advances in our understanding of tropical cyclone (TC) intensity change processes using aircraft data are summarized here. The focus covers a variety of spatiotemporal scales, regions of the TC inner core, and stages of the TC lifecycle, from preformation [...] Read more.
Recent (past ~15 years) advances in our understanding of tropical cyclone (TC) intensity change processes using aircraft data are summarized here. The focus covers a variety of spatiotemporal scales, regions of the TC inner core, and stages of the TC lifecycle, from preformation to major hurricane status. Topics covered include (1) characterizing TC structure and its relationship to intensity change; (2) TC intensification in vertical shear; (3) planetary boundary layer (PBL) processes and air–sea interaction; (4) upper-level warm core structure and evolution; (5) genesis and development of weak TCs; and (6) secondary eyewall formation/eyewall replacement cycles (SEF/ERC). Gaps in our airborne observational capabilities are discussed, as are new observing technologies to address these gaps and future directions for airborne TC intensity change research. Full article
(This article belongs to the Special Issue Rapid Intensity Changes of Tropical Cyclones)
Show Figures

Figure 1

Figure 1
<p>Intensity forecast errors (kt) from consensus of two statistical/dynamical models (Decav-SHIPS, Logistical Growth Equation Model) and two regional deterministic models (intelpolated HWRE, COAMPS-TC) for TCs between 2015 and 2017. Black line shows intensity forecast errors for all TCs not undergoing RI; red line shows errors for all TC’s undergoing RI.</p>
Full article ">Figure 2
<p>Composite azimuthally-averaged fields of three-dimensional analyses from cases in Rogers et al. (2012). (<b>a</b>) tangential wind (m s<sup>−1</sup>); (<b>b</b>) Reflectivity (dBZ); (<b>c</b>) Radial wind (m s<sup>−1</sup>). (<b>d</b>) Relative vorticity (×10<sup>−4</sup> s<sup>−1</sup>); (<b>e</b>) Vertical velocity (m s<sup>−1</sup>); (<b>f</b>) Horizontal divergence (×10<sup>−4</sup> s<sup>−1</sup>). Data from a minimum of 20 analyses are required for plotting. All composites plotted as a function of normalized radius r* and height above ground level. The dashed line denotes the axis of peak axisymmetric tangential wind from 0.5- to 10-km altitude calculated from the composite in (<b>a</b>). Adapted from [<a href="#B21-atmosphere-12-00650" class="html-bibr">21</a>].</p>
Full article ">Figure 3
<p>Quadrant-average cross sections of shear-relative composite reflectivity (dBZ, shaded), vertical velocity (m s<sup>−1</sup>, black contours), and radial velocity (m s<sup>−1</sup>, gray contours) for cases with peak winds &gt;31 m s<sup>−1</sup> and 850–200 hPa vertical shear &gt;7 m s<sup>−1</sup>. The quadrants are arranged such that the shear vector points to the right of the page. Regions of negative vertical motion are highlighted by the 0 and −0.25 m s<sup>−1</sup> dashed black contours. Contours of positive vertical motion (solid black) are drawn at 0.25, 0.5, 1, 1.5, 2, and 2.5 m s<sup>−1</sup>. The contour interval for radial inflow (dashed gray) and outflow (solid gray) is 1 m s<sup>−1</sup> (zero contour omitted). The radial coordinate r* represents the radius normalized bythe RMWat 2 km altitude. Adapted from [<a href="#B30-atmosphere-12-00650" class="html-bibr">30</a>].</p>
Full article ">Figure 4
<p>(<b>a</b>) Three-dimensional schematic depicting the vertical motion distribution in each quadrant. The environmental shear vector is denoted by an arrow pointing toward the top of the figure, and quadrants are labeled according to their direction relative to the shear vector (DR, DL, UL, and UR). Concentric circles below the clouds show the locations of the eyewall normalized radius values between 0.75 and 1.25), which connect to gray dashed boxes encompassing the eyewall up to 10-km altitude. Vertical arrows denote the vertical motion distribution, where size is proportional to magnitude. (<b>b</b>) As in (<b>a</b>), except illustrating the mean circulation that occurs when intense updrafts are present. Arrow width corresponds to composite velocity magnitude, as noted in the legend at the bottom. Dashed arrows refer to features that were present in the composite analysis but were shown to be not significant. (<b>c</b>) As in (<b>b</b>), except illustrating the mean circulation that occurs when intense downdrafts are present. The question mark denotes the structure existing UR is weak, disorganized, and not robust. Adapted from [<a href="#B31-atmosphere-12-00650" class="html-bibr">31</a>].</p>
Full article ">Figure 5
<p>Frequency of convective burst points as a function of normalized radial location (bar chart, %) overlain on normalized radius-height plot of composite-mean axisymmetric vertical vorticity (color shaded, ×10<sup>−4</sup> s<sup>−1</sup>) for the (<b>a</b>) intensifying (IN) and (<b>b</b>) steady-state (SS) TC composites. Adapted from Rogers et al. (2013).</p>
Full article ">Figure 6
<p>The 24-h intensity change (m s<sup>−1</sup>) of North Atlantic TCs with an inner-core lightning burst (ICLB) relative to the RMW (r/RMW; logarithmic scale) with aircraft reconnaissance within 1 h of the observed ICLB. The error bars indicate one standard deviation of 10,000 random errors (+0.2°) added to the NHC best-track position. Adapted from [<a href="#B47-atmosphere-12-00650" class="html-bibr">47</a>].</p>
Full article ">Figure 7
<p>Three-dimensional schematics summarizing the hypothesized hindrances to precipitation symmetry intropical cyclones. The top schematic depicts the more asymmetric case (Cristobal), while the bottomschematic depicts the more symmetric case (Bertha). Adapted from [<a href="#B60-atmosphere-12-00650" class="html-bibr">60</a>].</p>
Full article ">Figure 8
<p>(<b>a</b>) Locations of all points in Doppler analysis in a radial coordinate system normalized by the 2-km RMW (x*, y*) where peak upward motion in the 8–16-km layer is greater than 3 m s<sup>−1</sup> (green ×s) and 5 m s<sup>−1</sup> (black dots) for all center passes during the 14 Sep (intensifying) mission. (<b>b</b>) As in (<b>a</b>), but forthe 16 Sep (steady-state) mission. Arrow denotes shear direction; shear-relative quadrants are labeled. Adapted from Rogers et al. (2016) [<a href="#B67-atmosphere-12-00650" class="html-bibr">67</a>].</p>
Full article ">Figure 9
<p>Average and 95% confidence intervals of the M (solid) and 20-dBZ (dashed) slopes for (<b>a</b>) intensifying and (<b>b</b>) weakening/steady cases. Higher (lower) values of slope denote a more sloped (upright) configuration. Inset in each panel shows schematic illustrating relationship between momentum and reflectivitym surfaces for each intensity change category. Adapted from Hazelton et al. [<a href="#B26-atmosphere-12-00650" class="html-bibr">26</a>].</p>
Full article ">Figure 10
<p>(<b>a</b>) Boxplot comparison of the maximum updraft magnitude between the shear-relative quadrants and TC intensity change—intensifying (IN) and steady-state (SS); (<b>b</b>) as in (<b>a</b>), but for the height of the maximumupdraft; (<b>c</b>) as in (<b>a</b>), but for the maximum height of the 15-dBZ echo top. In the boxplots, the solid line in the box denotes the median value, while the upper and lower edges of the box represent 75th and 25th quartiles, respectively. The difference between the 75th and 25th percentiles represents the interquartile range (IOR). Whiskers extending above and below the box represent either 1.5 × IOR above the 75th and below the 25th percentiles respectively, or to the maximum/minimum values. Beyond the whiskers, values are statisticaloutliers and are represented as plus signs. Adapted from Wadler et al. (2018) [<a href="#B68-atmosphere-12-00650" class="html-bibr">68</a>].</p>
Full article ">Figure 11
<p>Conceptual diagram highlighting the measurements and analysis from the HIWRAP, HAMSR, and WP-3D TDR during the Hurricane Karl (2010) sampling. The arrows at lower levels represent the mesoscale flow and the arrows in the clouds represent the convective-scale flow. Red arrows indicate anomalously warm, buoyant air with blue arrows the opposite. Adapted from [<a href="#B72-atmosphere-12-00650" class="html-bibr">72</a>].</p>
Full article ">Figure 12
<p>Composite analysis of (<b>a</b>) the relative tangential wind (m s<sup>−1</sup>: contour interval 2 m s<sup>−1</sup>), (<b>b</b>) radial wind (m s<sup>−1</sup>; contour interval 2 m s<sup>−1</sup>), and (<b>c</b>) virtual potential temperature (K; contour interval 0.5 K) as a function of altitude and the normalized radius to the storm center. The thick black lines in (<b>a</b>) are the 40 and 50 m s<sup>−1</sup> contours and the black dashed line depicts the height of the maximum wind speed varying with radius. The thick black line in (<b>b</b>) depicts the inflow laver height defined as the height where the radial wind speed is 10% of the peak inflow. In (<b>b</b>), negative values are contoured in dashed lines. The white dashed line in (<b>b</b>) represents the zero contour, and the black × represents the location of the maximum tangential wind speed. The thick black lines in (<b>c</b>) are the 305- and 310-K contours. Adapted from Zhang et al. (2011) [<a href="#B78-atmosphere-12-00650" class="html-bibr">78</a>].</p>
Full article ">Figure 13
<p>(<b>a</b>) Plots of the wavenumber-1 component of equivalent potential temperature (shaded; K) at 50 m from dropsondes and contours of vertical velocity measured by TDR (contours, 0.2 m s<sup>−l</sup> contour interval) at 1.5 km ltitude. Gray contours represent updrafts and black contours represent downdrafts. Shear direction is shown by the black arrow. The thick dashed black line represents the radius of maximum wind speed. Adapted from [<a href="#B82-atmosphere-12-00650" class="html-bibr">82</a>].</p>
Full article ">Figure 14
<p>Distributions of surface (<b>a</b>) latent heat flux and (<b>b</b>) sensible heat flux in non-intensifying (NI; blue), slowly intensifying (SI; green), and rapidly intensifying (RI; red) cases in each shear-relative quadrant. Boxes encompass the 25th, 50th, and 75th percentiles, while whiskers extend to the 10th and 90th percentile values. Brackets denote the differences in the means of the distributions that were significant at the 95% confidence level as determined by a bootstrap test. Adapted from [<a href="#B84-atmosphere-12-00650" class="html-bibr">84</a>].</p>
Full article ">Figure 15
<p>Plots of (<b>a</b>) P-3 aircraft (N43) horizontal track varving with altitude in color and dropsonde locations (squares) released by a second P-3 aircraft (N42) flying above N43, (<b>b</b>) vertical profiles of wind speed from dropsondes with the mean wind profile (dashed line) and smoothed mean profile (black line), (<b>c</b>) vertical profile of the vertical wind shear with low-level aircraft altitudes (circles), and (<b>d</b>) vertical eddy diffusivity from a multitude of dropsonde and flight-level gmeasurement locations as a function of wind speed. Reproduced and adapted from [<a href="#B99-atmosphere-12-00650" class="html-bibr">99</a>,<a href="#B100-atmosphere-12-00650" class="html-bibr">100</a>].</p>
Full article ">Figure 16
<p>(<b>a</b>) Vertical cross section of potential temperature (°C; filled contours) and the cold-point tropopause height (green lines) along the high-altitude WB-57 transects over the center of Hurricane Patricia (2015) while it was rapidly intensifying on 22 Oct 2015. Dropsonde locations and compass directions indicated. Dashed vertica lines mark the storm center and hatching indicates regions of missing values, where linear interpolation is performed in the radial direction. (<b>b</b>) As in (<b>a</b>), but vertical cross sections of the potential temperature anomaly (°C) for same transect. (<b>c</b>) As in (<b>a</b>), but storm-relative radial flow (shaded, m s<sup>−1</sup>) for same transect. Adapted from [<a href="#B112-atmosphere-12-00650" class="html-bibr">112</a>].</p>
Full article ">Figure 17
<p>Meso-α-scale (600 km × 600 km) circulation tendencies with height at (<b>a</b>) 0000 and (<b>b</b>) 1200 UTC 13 Sep. Budget terms are mean stretching tendency (solid blue), eddy flux tendency (dash-dotted red), tilting tendency (long-dashed gray), friction tendency (short-dashed purple), and net tendency calculated by summation of thecomponent tendencies (thick solid black). Adapted from [<a href="#B124-atmosphere-12-00650" class="html-bibr">124</a>].</p>
Full article ">Figure 18
<p>(<b>a</b>) Plot of 20 dBZ echo top heights (shaded, km) and storm-relative winds (vectors, m s<sup>−1</sup>) at 2-km altitude forindividual center passes centered at 1742 UTC August 31 in Tropical Depression #9 (future Hermine); (<b>b</b>) As in (<b>a</b>), but for winds at 5-km altitude; (<b>c</b>) As in (<b>a</b>), but for 2034 UTC August 31; (<b>d</b>) As in (<b>b</b>), but for 2034UTC August 31; (<b>e</b>) As in (<b>a</b>), but for 2310 UTC August 31; (<b>f</b>) As in (<b>b</b>), but for 2310 UTC August 31. “L”and “M” denote locations of subjectively-determined circulation centers at 2- and 5-km altitudes, respectively. Lighter, smaller letters denote locations from previous center passes. Inset in lower right corner in (<b>a</b>,<b>c</b>,<b>e</b>) denotes SHIPS-derived 850–200 hPa shear. Adapted from [<a href="#B135-atmosphere-12-00650" class="html-bibr">135</a>].</p>
Full article ">Figure 19
<p>Azimuthally averaged tangential wind (m s<sup>−1</sup>; shaded) and vertical velocity (contours) obtained from (<b>a</b>) the 2358 UTC 3 Sep center pass. Positive vertical velocity values are shown in the solid contours every 0.5 m s<sup>−1</sup>, while negative vertical velocity values are shown in the dashed contours every 0.2 m s<sup>−1</sup>. The inset of the panel depicts the magnitude of the maximum azimuthally averaged tangential wind (V<sub>T3</sub>; m s<sup>−1</sup>) and RMW (km) at a height of 3 km. (<b>b</b>) As in (<b>a</b>), but for the 0904 UTC 4 Sep center pass. (<b>c</b>) As in (<b>a</b>), but for the 1027 UTC 4 Sep center pass. (<b>d</b>) As in (<b>a</b>), but for the 1142 UTC 4 Sep center pass. (<b>e</b>) As in (<b>a</b>), but for the 1313 UTC 4 Sep center pass. (<b>f</b>) As in (<b>a</b>), but for the 2140 UTC 4 Sep center pass. For each flight, only radii withat least 33% azimuthal coverage are shown. Adapted from [<a href="#B140-atmosphere-12-00650" class="html-bibr">140</a>].</p>
Full article ">Figure 20
<p>Number of publications from 2006–2020 using NOAA IFEX aircraft data that include NOAA/AOML/HRD and affiliated scientists, stratified by primary TC physical process topic addressed. Adapted from [<a href="#B152-atmosphere-12-00650" class="html-bibr">152</a>].</p>
Full article ">
27 pages, 20853 KiB  
Article
High NO2 Concentrations Measured by Passive Samplers in Czech Cities: Unresolved Aftermath of Dieselgate?
by Michal Vojtisek-Lom, Miroslav Suta, Jitka Sikorova and Radim J. Sram
Atmosphere 2021, 12(5), 649; https://doi.org/10.3390/atmos12050649 - 19 May 2021
Cited by 5 | Viewed by 4371
Abstract
This work examines the effects of two problematic trends in diesel passenger car emissions—increasing NO2/NOx ratio by conversion of NO into NO2 in catalysts and a disparity between the emission limit and the actual emissions in everyday driving—on ambient [...] Read more.
This work examines the effects of two problematic trends in diesel passenger car emissions—increasing NO2/NOx ratio by conversion of NO into NO2 in catalysts and a disparity between the emission limit and the actual emissions in everyday driving—on ambient air quality in Prague. NO2 concentrations were measured by 104 membrane-closed Palmes passive samplers at 65 locations in Prague in March–April and September–October of 2019. NO2 concentrations measured by city stations during those periods were comparable with the average values during 2016–2019. The average measured NO2 concentrations at the selected locations, after correcting for the 18.5% positive bias of samplers co-located with a monitoring station, were 36 µg/m3 (range 16–69 µg/m3, median 35 µg/m3), with the EU annual limit of 40 µg/m3 exceeded at 32% of locations. The NO2 concentrations have correlated well (R2 = 0.76) with the 2019 average daily vehicle counts, corrected for additional emissions due to uphill travel and intersections. In addition to expected “hot-spots” at busy intersections in the city center, new ones were identified, i.e., along a six-lane road V Holešovičkách. Comparison of data from six monitoring stations during 15 March–30 April 2020 travel restrictions with the same period in 2016–2019 revealed an overall reduction of NO2 and even a larger reduction of NO. The spatial analysis of data from passive samplers and time analysis of data during the travel restrictions both demonstrate a consistent positive correlation between traffic intensity and NO2 concentrations along/near the travel path. The slow pace of NO2 reductions in Prague suggests that stricter vehicle NOx emission limits, introduced in the last decade or two, have so far failed to sufficiently reduce the ambient NO2 concentrations, and there is no clear sign of remedy of Dieselgate NOx excess emissions. Full article
(This article belongs to the Special Issue Ambient Air Quality in the Czech Republic)
Show Figures

Figure 1

Figure 1
<p>Locations of the passive samplers and air quality monitoring stations used for comparison in this study. Photo of a sampler is shown in the upper right corner. (Map source: <a href="http://www.mapy.cz" target="_blank">www.mapy.cz</a> (accessed on 18 May 2021), © Seznam.cz, a.s., used with permission).</p>
Full article ">Figure 2
<p>Comparison of passive sampler reported NO<sub>2</sub> concentrations to the corresponding average values from corresponding monitoring stations. Larger points circled in red denote the colocation of the sampler at the monitoring station.</p>
Full article ">Figure 3
<p>Comparison of monitoring station NO<sub>2</sub> averages during sampling periods with four-year average.</p>
Full article ">Figure 4
<p>Comparison of spring and fall NO<sub>2</sub> concentrations.</p>
Full article ">Figure 5
<p>Relationship between traffic intensity and NO<sub>2</sub> concentrations measured by passive samplers in spring and fall of 2019 and by the national monitoring network (average of 2016–2019).</p>
Full article ">Figure 6
<p>Relationship between traffic intensity and NO<sub>2</sub> concentrations measured by passive samplers (average of all measurement periods) and by the national monitoring network (average of 2016–2019).</p>
Full article ">Figure 7
<p>Relationship between adjusted traffic intensity (traffic count × (1 + fraction of vehicles travelling uphill + 3 × fraction of vehicles stopping at an intersection)) and NO<sub>2</sub> concentrations measured by passive samplers (average of all measurement periods) and by the national monitoring network (average of 2016–2019).</p>
Full article ">Figure 8
<p>Relationship between adjusted traffic intensity (traffic count × (1 + fraction of vehicles travelling uphill + 3 × fraction of vehicles stopping at an intersection)) and NO<sub>2</sub> concentrations measured by passive samplers (average of all measurement periods) and by the national monitoring network (average of 2016–2019). EU annual limit of 40 µg/m<sup>3</sup> NO<sub>2</sub> shown as a red line, Swiss federal limit of 30 µg/m<sup>3</sup> NO<sub>2</sub> shown as a dotted green line.</p>
Full article ">
18 pages, 4334 KiB  
Article
Improving the Retrieval of Cloudy Atmospheric Profiles from Brightness Temperatures Observed with a Ground-Based Microwave Radiometer
by Qing Li, Ming Wei, Zhenhui Wang, Sulin Jiang and Yanli Chu
Atmosphere 2021, 12(5), 648; https://doi.org/10.3390/atmos12050648 - 19 May 2021
Cited by 5 | Viewed by 2656
Abstract
Atmospheric temperature and humidity retrievals from ground-based microwave remote sensing are useful in a variety of meteorological and environmental applications. Though the influence of clouds is usually considered in current retrieval algorithms, the resulting temperature and humidity estimates are still biased high in [...] Read more.
Atmospheric temperature and humidity retrievals from ground-based microwave remote sensing are useful in a variety of meteorological and environmental applications. Though the influence of clouds is usually considered in current retrieval algorithms, the resulting temperature and humidity estimates are still biased high in overcast conditions compared to radiosonde observations. Therefore, there is a need to improve the quality of retrievals in cloudy conditions. This paper presents an approach to make brightness temperature (TB) correction for cloud influence before the data can be used in the inversion of vertical profiles of atmospheric temperature and humidity. A three-channel method is proposed to make cloud parameter estimation, i.e., of the total 22 channels of the ground-based radiometer, three are adopted to set up a relationship between cloud parameters and brightness temperatures, so that the observations from the three channels can be used to estimate cloud thickness and water content and complete the cloud correction for the rest of the channels used in the retrieval. Based on two years of data from the atmosphere in Beijing, a comparison of the retrievals with radiosonde observations (RAOB) shows: (1) the temperature retrievals from this study have a higher correlation with RAOB and are notably better than in the vendor-provided LV2. The bias of the temperature retrievals from this study is close to zero at all heights, and the RMSE is greatly reduced from >5 °C to <2 °C in the layer, from about 1.5 km up to 5 km. The temperature retrievals from this study have higher correlation with RAOB data compared to the vendor-provided LV2, especially at and above a 2 km height. (2) The bias of the water vapor density profile from this study is near to zero, while the LV2 has a positive bias as large as 4 g/m3. The RMSE of the water vapor density profile from this study is <2 g/m3, while the RMSE for LV2 is as large as 10 g/m3. That is, both the bias and RMSE from this study are evidently less than the LV2, with a greater improvement in the lower troposphere below 5 km. Correlation with RAOB is improved even more for the water vapor density. The correlation of the retrievals from this study increases to one within the boundary layer, but the correlation of LV2 with RAOB is only 0.8 at 0.5 km height, 0.7 at 1 km, and even less than 0.5 at 2 km. (3) A parameter named the Cloud Impact index, determined by cloud water concentration and cloud thickness, together with the cloud base height, has been defined to show that both BIAS and RMSE of “high-CI subsample” are larger than those of the “low-CI subsample”, indicating that high-CI cloud has a higher impact on the retrievals and the correction for cloud influence is more necessary. Full article
(This article belongs to the Section Atmospheric Techniques, Instruments, and Modeling)
Show Figures

Figure 1

Figure 1
<p>Scatter plots showing the dependence of cloud radiation contribution (K) on (<b>left</b> panels) cloud water content (g/m<sup>3</sup>) and (<b>right</b> panels) cloud thickness (km) for channels 2, 7 and 10 (sequentially from top to bottom) of the microwave radiometer.</p>
Full article ">Figure 1 Cont.
<p>Scatter plots showing the dependence of cloud radiation contribution (K) on (<b>left</b> panels) cloud water content (g/m<sup>3</sup>) and (<b>right</b> panels) cloud thickness (km) for channels 2, 7 and 10 (sequentially from top to bottom) of the microwave radiometer.</p>
Full article ">Figure 2
<p>The technical flow chart for this study. (<b>a</b>) Flow-chart for the “3-channel cloud parameter estimation” and cloud impact identification. (<b>b</b>) Flow-chart for quantitative estimation of the cloud brightness temperature contribution and its correction in front of the inversion procedure for atmospheric temperature and water vapor density profiles. <span class="html-italic">CI</span> is equal or close to zero in the case when cloud impact can be neglected and TB correction for cloud is not necessary, as shown by the dotted box in <a href="#atmosphere-12-00648-f002" class="html-fig">Figure 2</a>a, which leads a “goto” circled A shown in <a href="#atmosphere-12-00648-f002" class="html-fig">Figure 2</a>b.</p>
Full article ">Figure 3
<p>Comparison of bias and RMSE of (<b>a</b>) temperature and (<b>b</b>) humidity profiles constructed from <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>f</mi> <mo>⇀</mo> </mover> <mo>−</mo> </msub> <msub> <mi>T</mi> <mrow> <mi>B</mi> <mi>O</mi> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> retrievals (red lines) with those from LV2 products (green lines). (<b>c</b>) Comparison of correlation coefficients between temperature and water vapor density profiles and radiosonde. Sample size = 109.</p>
Full article ">Figure 4
<p>Comparison of atmospheric temperature (<b>left</b> panels) and humidity (<b>right</b> panels) profiles from this study (red curves) and LV2 products (green lines) with RAOB (blue lines) for the four typical cases. (<b>a</b>) Case 1 at 00:00 UTC on 17 November 2010; (<b>b</b>) Case 2 at 00:00 UTC on 18 November 2010; (<b>c</b>) Case 3 at 00:00 UTC on 30 July 2010; (<b>d</b>) Case 4 at 00:00 UTC on 17 November 2011.</p>
Full article ">Figure 4 Cont.
<p>Comparison of atmospheric temperature (<b>left</b> panels) and humidity (<b>right</b> panels) profiles from this study (red curves) and LV2 products (green lines) with RAOB (blue lines) for the four typical cases. (<b>a</b>) Case 1 at 00:00 UTC on 17 November 2010; (<b>b</b>) Case 2 at 00:00 UTC on 18 November 2010; (<b>c</b>) Case 3 at 00:00 UTC on 30 July 2010; (<b>d</b>) Case 4 at 00:00 UTC on 17 November 2011.</p>
Full article ">Figure 5
<p>The performance of the retrieved atmospheric temperature and humidity profiles for the two kinds of cloud types in terms of Cloud Impact index. (<b>a</b>,<b>b</b>) BIAS and RMS of the vendor-provided LV2 products, (<b>c</b>,<b>d</b>) BIAS and RMS of the retrievals from <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mrow> <mi>B</mi> <mi>O</mi> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> by this study.</p>
Full article ">Figure 5 Cont.
<p>The performance of the retrieved atmospheric temperature and humidity profiles for the two kinds of cloud types in terms of Cloud Impact index. (<b>a</b>,<b>b</b>) BIAS and RMS of the vendor-provided LV2 products, (<b>c</b>,<b>d</b>) BIAS and RMS of the retrievals from <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mrow> <mi>B</mi> <mi>O</mi> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> by this study.</p>
Full article ">
22 pages, 13192 KiB  
Article
Field Study of Metal Oxide Semiconductor Gas Sensors in Temperature Cycled Operation for Selective VOC Monitoring in Indoor Air
by Tobias Baur, Johannes Amann, Caroline Schultealbert and Andreas Schütze
Atmosphere 2021, 12(5), 647; https://doi.org/10.3390/atmos12050647 - 19 May 2021
Cited by 35 | Viewed by 5641
Abstract
More and more metal oxide semiconductor (MOS) gas sensors with digital interfaces are entering the market for indoor air quality (IAQ) monitoring. These sensors are intended to measure volatile organic compounds (VOCs) in indoor air, an important air quality factor. However, their standard [...] Read more.
More and more metal oxide semiconductor (MOS) gas sensors with digital interfaces are entering the market for indoor air quality (IAQ) monitoring. These sensors are intended to measure volatile organic compounds (VOCs) in indoor air, an important air quality factor. However, their standard operating mode often does not make full use of their true capabilities. More sophisticated operation modes, extensive calibration and advanced data evaluation can significantly improve VOC measurements and, furthermore, achieve selective measurements of single gases or at least types of VOCs. This study provides an overview of the potential and limits of MOS gas sensors for IAQ monitoring using temperature cycled operation (TCO), calibration with randomized exposure and data-based models trained with advanced machine learning. After lab calibration, a commercial digital gas sensor with four different gas-sensitive layers was tested in the field over several weeks. In addition to monitoring normal ambient air, release tests were performed with compounds that were included in the lab calibration, but also with additional VOCs. The tests were accompanied by different analytical systems (GC-MS with Tenax sampling, mobile GC-PID and GC-RCP). The results show quantitative agreement between analytical systems and the MOS gas sensor system. The study shows that MOS sensors are highly suitable for determining the overall VOC concentrations with high temporal resolution and, with some restrictions, also for selective measurements of individual components. Full article
(This article belongs to the Special Issue The Future of Air Quality Monitoring)
Show Figures

Figure 1

Figure 1
<p>Temperature cycle of the SGP30.</p>
Full article ">Figure 2
<p>Schematic overview of the calibration setup with the gas mixing apparatus and the sensor systems.</p>
Full article ">Figure 3
<p>Schematic top view of the field test room, a standard office in our building. The locations of the trolly containing sensor systems and analytical measurements, release test and a fan are indicated (modified from [<a href="#B30-atmosphere-12-00647" class="html-bibr">30</a>]).</p>
Full article ">Figure 4
<p>Flowchart of the data evaluation for an initial and drift-compensated regression model.</p>
Full article ">Figure 5
<p>PLSR model for quantification of VOC<sub>sum</sub> for (<b>a</b>) training and testing with data from the initial calibration (initial regression model), (<b>b</b>) training with initial calibration, testing with 1st recalibration, (<b>c</b>) the drift compensated regression model (training with initial calibration plus background only of the 1st recalibration), testing with extended range data from 1st recalibration as well as 2nd recalibration. Dashed lines indicate the root-mean-square-error of validation (RMSEV) based on the training data set and of testing (RMSET), respectively.</p>
Full article ">Figure 6
<p>RMSE of the models for different target VOCs for the initial (RMSEV<sub>IRM</sub>, RMSET<sub>IRM</sub>) and the drift compensated regression model (RMSEV<sub>DCRM</sub>, RMSET<sub>DCRM</sub>). For each model RMSEV for 10-fold validation during training (error bars indicate the standard deviation of the RMSEV for the different folds) and testing (20% holdout).</p>
Full article ">Figure 7
<p>Results recorded during four release tests with toluene (<b>a</b>), acetone (<b>b</b>), ethanol (<b>c</b>), and the simultaneous release of toluene, acetone, and ethanol (<b>d</b>). The upper graphs show the PLSR prediction of the MOS sensor model for the released gases and, if available, reference data (dots); all other signals are shown in the lower graphs. The signals are smoothed over five points (10 min). Numbers in parenthesis behind the released substances refer to the release tests, cf. <a href="#atmosphere-12-00647-t004" class="html-table">Table 4</a>.</p>
Full article ">Figure 8
<p>Results recorded during three release tests with hydrogen (<b>a</b>,<b>b</b>) and the simultaneous release of toluene and acetone (<b>c</b>). The upper graphs show the PLSR prediction of the MOS sensor model for the released gases and, if available, reference data (dots indicate distinct sampling times); all other signals are shown in the lower graphs. The signals are smoothed over five points (10 min). Numbers in parenthesis behind the released substances refer to the release tests, cf. <a href="#atmosphere-12-00647-t004" class="html-table">Table 4</a>.</p>
Full article ">Figure 9
<p>Results recorded during four release tests with toluene (<b>a</b>), m/p-xylene (<b>b</b>), simultaneous release of toluene and m/p-xylene (<b>c</b>), and isopropyl alcohol (<b>d</b>). The upper graphs show the PLSR prediction of the MOS sensor model for the same types of VOC and, if available, reference data (dots); all other signals are shown in the lower graphs. The signals are smoothed over five points (10 min). Numbers in parenthesis behind the released substances refer to the release tests, cf. <a href="#atmosphere-12-00647-t004" class="html-table">Table 4</a>.</p>
Full article ">Figure A1
<p>Flowchart of different data evaluation steps. Indices i: number of features, j: number of PLSR components and k: current fold during 10-fold cross validation.</p>
Full article ">
18 pages, 1297 KiB  
Article
Electric Field Multifractal Features in the High-Latitude Ionosphere: CSES-01 Observations
by Giuseppe Consolini, Virgilio Quattrociocchi, Giulia D’Angelo, Tommaso Alberti, Mirko Piersanti, Maria Federica Marcucci and Paola De Michelis
Atmosphere 2021, 12(5), 646; https://doi.org/10.3390/atmos12050646 - 19 May 2021
Cited by 36 | Viewed by 2662
Abstract
In the polar ionosphere, the electric field is characterized by broadband and power law spectral densities at small/short spatio-temporal scales, which support a possible turbulent nature of the electric field fluctuations. Here, we investigate the multifractal character of the full three-dimensional electric field [...] Read more.
In the polar ionosphere, the electric field is characterized by broadband and power law spectral densities at small/short spatio-temporal scales, which support a possible turbulent nature of the electric field fluctuations. Here, we investigate the multifractal character of the full three-dimensional electric field in the polar ionosphere as recorded on board the first Chinese Seismo-Electromagnetic Satellite (CSES-01). The results of our analysis prove a clear different degree of multifractality of the electric field fluctuations approaching either the polar cap trailing edge or the auroral region. The observed differences in the multifractal character are interpreted in terms of the different natures of the particle precipitation in the polar cap and in the auroral region. A possible link between the multifractal nature of electric field fluctuations, parallel to the geomagnetic field, and filamentation of field aligned currents (FACs) is established. Full article
Show Figures

Figure 1

Figure 1
<p>The CSES-01 trajectory for the selected time interval from 21:41 UT to 21:45 UT on 11 August 2018. The auroral display refers to observation by SSUSI instrument on board DMSP-F17 satellite at 21:45 UT. The red dashed curves refer to the upper and lower boundary of the auroral oval. The red curve refers to CSES-01 orbit. Reference system is AACGM (Lat, MLT) coordinate system.</p>
Full article ">Figure 2
<p>The ground measurements of H and Z magnetic field components at the MAW observatory.</p>
Full article ">Figure 3
<p>Electric field measurements, corrected for the <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="bold">E</mi> <mi>s</mi> </msub> <mo>=</mo> <msub> <mi mathvariant="bold">v</mi> <mi>s</mi> </msub> <mo>×</mo> <mi mathvariant="bold">B</mi> </mrow> </semantics></math>, collected by CSES-01 EFD instrument. The three components are in the geographical (GEO) reference system. The vertical dashed lines enclose the two selected time intervals considered in this study (see text for more details).</p>
Full article ">Figure 4
<p>The CSES-01 trajectory for the two selected time intervals superimposed on the instantaneous convection cells as reconstructed by SuperDARN observations in Antarctica. The green and cyan parts of CSES-01 trajectory (red line) refer to interval #1 and interval #2, respectively. Reference system is AACGM MLat-MLT.</p>
Full article ">Figure 5
<p>The electric field components in the minimum variance reference frame for the two selected time intervals. From top to bottom the components are along the minimum, the middle and the maximum variance directions, respectively.</p>
Full article ">Figure 6
<p>Trace of the power spectral density (PSD) of the electric field in the two selected time intervals. The solid and dashed lines refer to power laws <math display="inline"><semantics> <msup> <mi>f</mi> <mrow> <mo>−</mo> <mi>β</mi> </mrow> </msup> </semantics></math>, characterized by exponents <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>5</mn> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>, respectively. The secondary bottom axis is computed considering the satellite speed <math display="inline"><semantics> <msub> <mi>v</mi> <mi>s</mi> </msub> </semantics></math> as <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mi>f</mi> <mo>/</mo> <msub> <mi>v</mi> <mi>s</mi> </msub> </mrow> </semantics></math>, being <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>s</mi> </msub> <mo>≃</mo> <mn>8</mn> </mrow> </semantics></math> km/s.</p>
Full article ">Figure 7
<p>The PSD of electric field fluctuations in the two selected intervals in the Min/Med variance plane (eft panel) and in the maximum variance direction (right panel), respectively. The dashed and solid lines are power laws <math display="inline"><semantics> <msup> <mi>f</mi> <mrow> <mo>−</mo> <mi>α</mi> </mrow> </msup> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math> for dashed lines, and <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.5</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.8</mn> </mrow> </semantics></math> for solid lines in left and right panels, respectively.</p>
Full article ">Figure 8
<p><math display="inline"><semantics> <mrow> <msub> <mo>∑</mo> <mi>i</mi> </msub> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>,</mo> <mi>τ</mi> <mo>)</mo> </mrow> <mo form="prefix">ln</mo> <msub> <mi>p</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>τ</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> (upper panel) and <math display="inline"><semantics> <mrow> <msub> <mo>∑</mo> <mi>i</mi> </msub> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>,</mo> <mi>τ</mi> <mo>)</mo> </mrow> <mo form="prefix">ln</mo> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>,</mo> <mi>τ</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> (lower panel) as a function of <math display="inline"><semantics> <mrow> <mo form="prefix">ln</mo> <mi>τ</mi> </mrow> </semantics></math> at different <span class="html-italic">q</span> values in the case of the first selected time interval and for the the maximum variance direction component. Solid and dashed lines are linear fits. The vertical dashed line refers to the separation scale <math display="inline"><semantics> <mrow> <msup> <mi>τ</mi> <mo>*</mo> </msup> <mo>=</mo> <mn>0.075</mn> </mrow> </semantics></math> s between the two different dynamical ranges.</p>
Full article ">Figure 9
<p>The multifractal singularity spectra, <math display="inline"><semantics> <mrow> <mi>f</mi> <mo>(</mo> <mi>α</mi> <mo>)</mo> </mrow> </semantics></math>, relative to the two range of scales for the two selected time intervals. The results refer to the maximum variance direction.</p>
Full article ">Figure 10
<p>The multifractal singularity spectra, <math display="inline"><semantics> <mrow> <mi>f</mi> <mo>(</mo> <mi>α</mi> <mo>)</mo> </mrow> </semantics></math>, relative to the two range of scales for the two selected time intervals. The results refer to the Min/Med variance plane.</p>
Full article ">
23 pages, 2565 KiB  
Article
Improving Air Pollutant Metal Oxide Sensor Quantification Practices through: An Exploration of Sensor Signal Normalization, Multi-Sensor and Universal Calibration Model Generation, and Physical Factors Such as Co-Location Duration and Sensor Age
by Kristen Okorn and Michael Hannigan
Atmosphere 2021, 12(5), 645; https://doi.org/10.3390/atmos12050645 - 19 May 2021
Cited by 11 | Viewed by 3574
Abstract
As low-cost sensors have become ubiquitous in air quality measurements, there is a need for more efficient calibration and quantification practices. Here, we deploy stationary low-cost monitors in Colorado and Southern California near oil and gas facilities, focusing our analysis on methane and [...] Read more.
As low-cost sensors have become ubiquitous in air quality measurements, there is a need for more efficient calibration and quantification practices. Here, we deploy stationary low-cost monitors in Colorado and Southern California near oil and gas facilities, focusing our analysis on methane and ozone concentration measurement using metal oxide sensors. In comparing different sensor signal normalization techniques, we propose a z-scoring standardization approach to normalize all sensor signals, making our calibration results more easily transferable among sensor packages. We also attempt several different physical co-location schemes, and explore several calibration models in which only one sensor system needs to be co-located with a reference instrument, and can be used to calibrate the rest of the fleet of sensor systems. This approach greatly reduces the time and effort involved in field normalization without compromising goodness of fit of the calibration model to a significant extent. We also explore other factors affecting the performance of the sensor system quantification method, including the use of different reference instruments, duration of co-location, time averaging, transferability between different physical environments, and the age of metal oxide sensors. Our focus on methane and stationary monitors, in addition to the z-scoring standardization approach, has broad applications in low-cost sensor calibration and utility. Full article
(This article belongs to the Special Issue Atmospheric Trace Gas Source Detection and Quantification)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Map and timeline of co-locations for methane.</p>
Full article ">Figure 2
<p>(<b>a</b>) Traditional co-location and deployment setup, used for the individual, z-scored individual, median, and sensor-specific normalization calibration approaches; (<b>b</b>) One-Cal method—generating one pod calibration model to apply to all pods; (<b>c</b>) One-Hop method—using one co-located pod as a reference to generate individual pod calibrations.</p>
Full article ">Figure 3
<p>Light VOC metal oxide co-location sensor readings in their raw format (<b>left</b>) and after being z-scored (<b>right</b>).</p>
Full article ">Figure 4
<p>R<sup>2</sup> (<b>left</b>) and MBE vs. CRMSD target plot (<b>right</b>) for all sensor signal normalization and multi-pod calibration approaches studied for methane against: a reference-grade instrument (<b>top</b>); a pod as a secondary standard (<b>bottom</b>). The box-and-whiskers are each in the following order from left to right: Greeley calibration, Los Angeles calibration, Los Angeles validation, Wiggins validation* (* bottom only).</p>
Full article ">Figure 5
<p>R<sup>2</sup> (<b>left</b>) and MBE vs. CRMSD target plot (<b>right</b>) for all sensor signal normalization and multi-pod calibration approaches studied for ozone against a reference-grade instrument. The box-and-whiskers are each in the following order from left to right: Boulder validation (Boulder Campus), Boulder calibration (South Boulder Creek), Shafter calibration, Shafter validation.</p>
Full article ">Figure 6
<p>Methane R<sup>2</sup> boxplot (<b>left</b>) and MBE vs. CRMSD target plot (<b>right</b>) for calibration models applied to sensors of different ages (light) and deployed in different environments (dark shading). Within each boxplot column, calibration data is leftmost, different deployment environment data is center, and different sensor age data is rightmost.</p>
Full article ">Figure 7
<p>R<sup>2</sup> (<b>left</b>) and MBE vs. CRMSD target plot (<b>right</b>) for individually calibrated pods as compared with two different reference instruments, using different durations of co-location. In each boxplot column, the leftmost box-and-whisker represents the Picarro data, while the rightmost is the Thermo data.</p>
Full article ">Figure 8
<p>Timeseries of Picarro and Thermo reference instruments, with the period of data used for analysis boxed in.</p>
Full article ">Figure 9
<p>Methane R<sup>2</sup> boxplot (<b>left</b>) and MBE vs. CRMSD target plot (<b>right</b>) for duration of co-located “hop” using the one-hop approach as tabulated against a regulatory-grade reference instrument and an individually calibrated pod utilized as a secondary standard. Within each box-and whisker column, the data is as follows from left to right: Greeley calibration, Los Angeles calibration, Wiggins validation (bottom plot only), and Los Angeles validation.</p>
Full article ">Figure 10
<p>Methane R<sup>2</sup> boxplot (<b>left</b>) and MBE vs. CRMSD target plot (<b>right</b>) for time averaging of co-located “hop” using the one-hop method as tabulated against a regulatory-grade reference instrument and an individually calibrated pod utilized as a secondary standard. Within each box-and whisker column, the data is as follows from left to right: Greeley calibration, Los Angeles calibration, Wiggins validation (bottom plot only), and Los Angeles validation.</p>
Full article ">
13 pages, 1580 KiB  
Article
Wind Shear Prediction from Light Detection and Ranging Data Using Machine Learning Methods
by Jingyan Huang, Michael Kwok Po Ng and Pak Wai Chan
Atmosphere 2021, 12(5), 644; https://doi.org/10.3390/atmos12050644 - 18 May 2021
Cited by 11 | Viewed by 3504
Abstract
The main aim of this paper is to propose a statistical indicator for wind shear prediction from Light Detection and Ranging (LIDAR) observational data. Accurate warning signal of wind shear is particularly important for aviation safety. The main challenges are that wind shear [...] Read more.
The main aim of this paper is to propose a statistical indicator for wind shear prediction from Light Detection and Ranging (LIDAR) observational data. Accurate warning signal of wind shear is particularly important for aviation safety. The main challenges are that wind shear may result from a sustained change of the headwind and the possible velocity of wind shear may have a wide range. Traditionally, aviation models based on terrain-induced setting are used to detect wind shear phenomena. Different from traditional methods, we study a statistical indicator which is used to measure the variation of headwinds from multiple headwind profiles. Because the indicator value is nonnegative, a decision rule based on one-side normal distribution is employed to distinguish wind shear cases and non-wind shear cases. Experimental results based on real data sets obtained at Hong Kong International Airport runway are presented to demonstrate that the proposed indicator is quite effective. The prediction performance of the proposed method is better than that by the supervised learning methods (LDA, KNN, SVM, and logistic regression). This model would also provide more accurate warnings of wind shear for pilots and improve the performance of Wind shear and Turbulence Warning System. Full article
(This article belongs to the Special Issue Low Level Windshear and Turbulence for Aviation Safety)
Show Figures

Figure 1

Figure 1
<p>Map of Hong Kong International Airport and Lantau Island (height contours: 100 m) [<a href="#B1-atmosphere-12-00644" class="html-bibr">1</a>], with the location of the LIDAR (red square). Runway corridors are shown as pink arrows with the names marked alongside.</p>
Full article ">Figure 2
<p>Diagram illustrating a LIDAR glide path scan along the 3-degree glide path for the western approach towards the north runway.</p>
Full article ">Figure 3
<p>Histograms of <span class="html-italic">k<sub>i,</sub></span><sub>1</sub> and <span class="html-italic">k<sub>i,</sub></span><sub>2</sub> for wind shear and non-wind shear cases in March 2015. (<b>left</b>) <span class="html-italic">k<sub>i,</sub></span><sub>1</sub> for wind shear and non-wind shear cases; (<b>right</b>) <span class="html-italic">k<sub>i,</sub></span><sub>2</sub> for wind shear and non-wind shear cases.</p>
Full article ">Figure 4
<p>QQ Plot of <span class="html-italic">k<sub>i,</sub></span><sub>1</sub> for wind shear and null wind shear cases in March 2015. (<b>a</b>) <span class="html-italic">k<sub>i,</sub></span><sub>1</sub> for wind shear; (<b>b</b>) <span class="html-italic">k<sub>i,</sub></span><sub>1</sub> for non-wind shear cases.</p>
Full article ">
17 pages, 5206 KiB  
Article
Using Copernicus Atmosphere Monitoring Service (CAMS) Products to Assess Illuminances at Ground Level under Cloudless Conditions
by William Wandji Nyamsi, Philippe Blanc, Dominique Dumortier, Ruben Mouangue, Antti Arola and Lucien Wald
Atmosphere 2021, 12(5), 643; https://doi.org/10.3390/atmos12050643 - 18 May 2021
Cited by 5 | Viewed by 2671
Abstract
Natural daylight is recognized as an important variable in the energy performance of buildings. A method that estimates the global illuminance received on a horizontal surface at ground level and its direct component at normal incidence under cloudless conditions is presented. The method [...] Read more.
Natural daylight is recognized as an important variable in the energy performance of buildings. A method that estimates the global illuminance received on a horizontal surface at ground level and its direct component at normal incidence under cloudless conditions is presented. The method uses the k-distribution method and the correlated-k approximation to compute a set of clearness indices integrated over 13 spectral bands covering the range 380–780 nm. A spectral resampling technique, including a spectral disaggregation and a spectral linear interpolation, is applied to these indices for providing a detailed set of solar irradiances at 1 nm in spectral resolution over the whole range. Then, these are weighted by the standardized CIE action spectrum for human eye for assessing the illuminance. Inputs to the method include the total column contents of ozone and water vapor as well as aerosol optical properties produced by the Copernicus Atmosphere Monitoring Service. Estimates of illuminance were compared to high-quality 1 min measurements of illuminance that were collected from two experimental sites located in two different climatic zones. A slight overestimation is observed for the global illuminance: the bias is between +1 klx and +3 klx, i.e., between +1% and +4% in relative value. The root mean square error varies between 5 klx (8%) and 6 klx (9%). The squared correlation coefficient ranges between 0.95 and 0.97. At the site providing the direct illuminance at normal incidence, the performance of the method is lower compared to global illuminance with a lower squared correlation coefficient of 0.53. The bias, relative bias, RMSE, and rRMSE are +7 klx, +9%, 12 klx, and 15%, respectively. The uncertainty of the method is of the order of the uncertainty of the measurements. The method offers accurate estimates of illuminance in cloudless conditions at high spatial and temporal resolutions useful for construction industries and operators as well as thermal simulation tools for optimal building design strategies. Full article
(This article belongs to the Section Meteorology)
Show Figures

Figure 1

Figure 1
<p>Illustration of the spectral resampling technique. The clearness indices obtained between 363 and 791 nm by two runs of libRadtran are shown in green and brown respectively for detailed calculations at 1 nm resolution and for the Kato et al. [<a href="#B12-atmosphere-12-00643" class="html-bibr">12</a>] scheme. The selected fine bands FB are shown by crosses in magenta. The linear interpolation provides the clearness indices every 1 nm drawn in blue.</p>
Full article ">Figure 2
<p>Two-dimensional (2D) histogram between measurements and estimates of global illuminance on horizontal surface at Golden. The color bar depicts the number of points in the area of 1.5 klx × 1.5 klx. The two thin magenta lines delimit the area of relative errors of ±8% from the measurements.</p>
Full article ">Figure 3
<p>Box plots of the ratio (<b>top</b>) of the estimated (Esti) to the measured (Meas) global illuminance on horizontal surface and difference (<b>bottom</b>) between estimated and measured global illuminance on horizontal surface with CAMS input data at Golden. The mean of the box is marked with a red dot. The 1st, 2nd, and 3rd quartiles are marked with a blue line. The notch in the box shows the 95% confidence interval around the median. The number of data for a given range are reported in pink number.</p>
Full article ">Figure 3 Cont.
<p>Box plots of the ratio (<b>top</b>) of the estimated (Esti) to the measured (Meas) global illuminance on horizontal surface and difference (<b>bottom</b>) between estimated and measured global illuminance on horizontal surface with CAMS input data at Golden. The mean of the box is marked with a red dot. The 1st, 2nd, and 3rd quartiles are marked with a blue line. The notch in the box shows the 95% confidence interval around the median. The number of data for a given range are reported in pink number.</p>
Full article ">Figure 4
<p>Diurnal variability of AOD at 550 nm from CAMS and measurements at Golden on 15 April 2017. SCI stands for selected cloudless instant.</p>
Full article ">Figure 5
<p>Two-dimensional (2D) histogram between measurements and estimates of global illuminance at Vaulx-en-Velin. The color bar depicts the number of points in the area of 1.5 klx × 1.5 klx. The two thin magenta lines delimit the area of relative errors of ±5 % from the measurements.</p>
Full article ">Figure 6
<p>Two-dimensional (2D) histogram between measurements and estimates of direct illuminance at normal incidence at Vaulx-en-Velin. The colorbar depicts the number of points in the area of 1.5 klx × 1.5 klx. The two thin magenta lines delimit the area of relative errors of ±5% from the measurements.</p>
Full article ">
14 pages, 5059 KiB  
Article
Aerosol Optical Properties around the East China Seas Based on AERONET Measurements
by Qianguang Tu, Zengzhou Hao, Yunwei Yan, Bangyi Tao, Chuyong Chung and Sumin Kim
Atmosphere 2021, 12(5), 642; https://doi.org/10.3390/atmos12050642 - 18 May 2021
Cited by 6 | Viewed by 2737
Abstract
Understanding aerosols optical properties over the oceans is vital for enhancing our knowledge of aerosol effects on climate and pollutant transport between continents. In this study, the characteristics of aerosol optical thickness (AOT) at 500 nm (τ500nm), Ångström exponent [...] Read more.
Understanding aerosols optical properties over the oceans is vital for enhancing our knowledge of aerosol effects on climate and pollutant transport between continents. In this study, the characteristics of aerosol optical thickness (AOT) at 500 nm (τ500nm), Ångström exponent for the wavelength pair 440–870 nm (α) and volume size distribution (VSD), are presented and analyzed over the East China seas based on the observations at four AERONET sites during 1999–2019. The main results are: (1) the mean τ500nm (α) value ranged from 0.31 to 0.36 (1.17–1.31); (2) the distribution of τ500nm (α) is similar to a log-normal distribution with a right-skewed long tail larger than 0.5 (closer to the normal distribution); (3) large AOT (τ500nm>0.6) was frequently observed in summer (June and July) and spring (March to May), followed by autumn and winter; (4) all aerosol types were observed, and urban/industrial aerosols and mixed types were dominant throughout the period. The atmospheric column aerosol was characterized by a bimodal lognormal size distribution with a fine mode at effective radius, Reff = 0.16 ± 0.01 μm, and coarse mode at Reff = 2.05 ± 0.1 μm. Full article
(This article belongs to the Special Issue Atmospheric Aerosol Optical Properties)
Show Figures

Figure 1

Figure 1
<p>The locations of the four AERONET sites.</p>
Full article ">Figure 2
<p>The monthly Level 2.0 data volume plot of aerosol optical thickness (<b>a</b>) and aerosol inversions (<b>b</b>) for each site used in this study.</p>
Full article ">Figure 3
<p>Frequency distribution of AOT (500 nm) for (<b>a</b>) Baengnyeong, (<b>b</b>) Anmyon, (<b>c</b>) Gosan_SNU, and (<b>d</b>) Cape_Fuguei_Station.</p>
Full article ">Figure 4
<p>Frequency distribution of Ångström exponent for (<b>a</b>) Baengnyeong, (<b>b</b>) Anmyon, (<b>c</b>) Gosan_SNU, and (<b>d</b>) Cape_Fuguei_Station.</p>
Full article ">Figure 5
<p>(<b>a</b>) Mean daily values of <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mrow> <mn>500</mn> <mi>nm</mi> </mrow> </msub> </mrow> </semantics></math>, (<b>b</b>) daily standard deviations of <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mrow> <mn>500</mn> <mi>nm</mi> </mrow> </msub> </mrow> </semantics></math>, (<b>c</b>) mean daily values of Ångström exponent, and (<b>d</b>) precipitable water vapor in the total atmospheric column at Baengnyeong, the Yellow Sea.</p>
Full article ">Figure 6
<p>(<b>a</b>) Mean daily values of <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mrow> <mn>500</mn> <mi>nm</mi> </mrow> </msub> </mrow> </semantics></math>, (<b>b</b>) daily standard deviations of <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mrow> <mn>500</mn> <mi>nm</mi> </mrow> </msub> </mrow> </semantics></math>, (<b>c</b>) mean daily values of Ångström exponent, and (<b>d</b>) precipitable water vapor in the total atmospheric column at Anmyon, the Yellow Sea.</p>
Full article ">Figure 7
<p>(<b>a</b>) Mean daily values of <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mrow> <mn>500</mn> <mi>nm</mi> </mrow> </msub> </mrow> </semantics></math>, (<b>b</b>) daily standard deviations of <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mrow> <mn>500</mn> <mi>nm</mi> </mrow> </msub> </mrow> </semantics></math>, (<b>c</b>) mean daily values of Ångström exponent, and (<b>d</b>) precipitable water vapor in the total atmospheric column at Gosan SNU, the East China Sea.</p>
Full article ">Figure 8
<p>(<b>a</b>) Mean daily values of <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mrow> <mn>500</mn> <mi>nm</mi> </mrow> </msub> </mrow> </semantics></math>, (<b>b</b>) daily standard deviations of <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mrow> <mn>500</mn> <mi>nm</mi> </mrow> </msub> </mrow> </semantics></math>, (<b>c</b>) mean daily values of Ångström exponent, and (<b>d</b>) precipitable water vapor in the total atmospheric column at Cape_Fuguei_Station, the East China Sea.</p>
Full article ">Figure 9
<p>Density plots of the Ångström exponent against the AOT (500 nm) for (<b>a</b>) Baengnyeong, (<b>b</b>) Anmyon, (<b>c</b>) Gosan_SNU, and (<b>d</b>) Cape_Fuguei_Station. Legend: MA = marine aerosol, DA = dust aerosol, UIA = urban/industrial aerosol, BBA = biomass burning aerosol, MIX = mixed aerosol.</p>
Full article ">Figure 10
<p>Multi-year average aerosol particle volume size distribution for Baengnyeong, Anmyon, Gosan_SNU, and Cape_Fuguei_Station.</p>
Full article ">
14 pages, 1996 KiB  
Article
Geochemical Characterization and Heavy Metal Sources in PM10 in Arequipa, Peru
by Jianghanyang Li, Greg Michalski, Elizabeth Joy Olson, Lisa R. Welp, Adriana E. Larrea Valdivia, Juan Reyes Larico, Francisco Alejo Zapata and Lino Morales Paredes
Atmosphere 2021, 12(5), 641; https://doi.org/10.3390/atmos12050641 - 18 May 2021
Cited by 5 | Viewed by 3424
Abstract
Particulate matter smaller than 10 μm (PM10) is an important air pollutant that adversely affects human health by increasing the risk of respiratory and cardiovascular diseases. Recent studies reported multiple extreme PM10 levels at high altitude Peruvian cities, which resulted [...] Read more.
Particulate matter smaller than 10 μm (PM10) is an important air pollutant that adversely affects human health by increasing the risk of respiratory and cardiovascular diseases. Recent studies reported multiple extreme PM10 levels at high altitude Peruvian cities, which resulted from a combination of high emissions and limited atmospheric circulation at high altitude. However, the emission sources of the PM10 still remain unclear. In this study, we collected PM10 samples from four sites (one industrial site, one urban site, and two rural sites) at the city of Arequipa, Peru, during the period of February 2018 to December 2018. To identify the origins of PM10 at each site and the spatial distribution of PM10 emission sources, we analyzed major and trace element concentrations of the PM10. Of the observed daily PM10 concentrations at Arequipa during our sampling period, 91% exceeded the World Health Organization (WHO) 24-h mean PM10 guideline value, suggesting the elevated PM10 strongly affected the air quality at Arequipa. The concentrations of major elements, Na, K, Mg, Ca, Fe, and Al, were high and showed little variation, suggesting that mineral dust was a major component of the PM10 at all the sites. Some trace elements, such as Mn and Mo, originated from the mineral dust, while other trace elements, including Pb, Sr, Cu, Ba, Ni, As and V, were from additional anthropogenic sources. The industrial activities at Rio Seco, the industrial site, contributed to significant Pb, Cu, and possibly Sr emissions. At two rural sites, Tingo Grande and Yarabamba, strong Cu emissions were observed, which were likely associated with mining activities. Ni, V, and As were attributed to fossil fuel combustion emissions, which were strongest at the Avenida Independencia urban site. Elevated Ba and Cu concentrations were also observed at the urban site, which were likely caused by heavy traffic in the city and vehicle brake wear emissions. Full article
Show Figures

Figure 1

Figure 1
<p>Map of the sampling sites. Industrial site (Rio Seco), urban site (Av. Independencia), and suburban sites (Tingo Grande and Yarabamba). Urban land cover appears grey in this image.</p>
Full article ">Figure 2
<p>PM<sub>10</sub> concentrations during the sampling period at all the sites. (<b>A</b>) Time series of observed PM<sub>10</sub> concentrations in this study; (<b>B</b>) box plot of PM<sub>10</sub> concentrations at each site.</p>
Full article ">Figure 3
<p>Box plot of major element concentrations (<b>A</b>: Fe, <b>B</b>: Na, <b>C</b>: K, <b>D</b>: Mg, in %) at each site.</p>
Full article ">Figure 4
<p>Box plot of the air concentrations (in ng/m<sup>3</sup>) of trace elements at each site.</p>
Full article ">Figure 5
<p>Box plot of mass-weighted concentrations (in ppm) of trace elements in the PM<sub>10</sub> at each site.</p>
Full article ">
15 pages, 7600 KiB  
Article
Accelerated Time and High-Resolution 3D Modeling of the Flow and Dispersion of Noxious Substances over a Gigantic Urban Area—The EMERGENCIES Project
by Olivier Oldrini, Patrick Armand, Christophe Duchenne, Sylvie Perdriel and Maxime Nibart
Atmosphere 2021, 12(5), 640; https://doi.org/10.3390/atmos12050640 - 18 May 2021
Cited by 5 | Viewed by 2724
Abstract
Accidental or malicious releases in the atmosphere are more likely to occur in built-up areas, where flow and dispersion are complex. The EMERGENCIES project aims to demonstrate the operational feasibility of three-dimensional simulation as a support tool for emergency teams and first responders. [...] Read more.
Accidental or malicious releases in the atmosphere are more likely to occur in built-up areas, where flow and dispersion are complex. The EMERGENCIES project aims to demonstrate the operational feasibility of three-dimensional simulation as a support tool for emergency teams and first responders. The simulation domain covers a gigantic urban area around Paris, France, and uses high-resolution metric grids. It relies on the PMSS modeling system to model the flow and dispersion over this gigantic domain and on the Code_Saturne model to simulate both the close vicinity and the inside of several buildings of interest. The accelerated time is achieved through the parallel algorithms of the models. Calculations rely on a two-step approach: the flow is computed in advance using meteorological forecasts, and then on-demand release scenarios are performed. Results obtained with actual meteorological mesoscale data and realistic releases occurring both inside and outside of buildings are presented and discussed. They prove the feasibility of operational use by emergency teams in cases of atmospheric release of hazardous materials. Full article
Show Figures

Figure 1

Figure 1
<p>Example of parallel settings: (<b>a</b>) domain decomposition and time frame parallel treatment for the PSWIFT model using 17 cores on a domain decomposed into 8 tiles; (<b>b</b>) domain decomposition and particle distribution for the PSPRAY model.</p>
Full article ">Figure 2
<p>View of the domain boundaries (green line) and the department boundaries (blue line) (© OpenStreetMap Contributors).</p>
Full article ">Figure 3
<p>3D view of the topography and buildings.</p>
Full article ">Figure 4
<p>Domain decomposition of the domain in 1088 tiles.</p>
Full article ">Figure 5
<p>Aerial view zoom of the city of Paris with the specific buildings of interest (red circles).</p>
Full article ">Figure 6
<p>Closeup view of the flow streamlines computed in the streets south and east of the train station. The wind speed ranges from 0 m/s in cyan to 2 m/s in magenta.</p>
Full article ">Figure 7
<p>3D view of the flow streamlines for the nested domain both in the vicinity of and inside the administration building. The wind speed ranges from 0 m/s in deep blue to 3.5 m/s in orange.</p>
Full article ">Figure 8
<p>View of the plume inside and in the vicinity of the museum: (<b>a</b>) 1 min after the release; (<b>b</b>) 3 min; (<b>c</b>) 4 min; (<b>d</b>) 6 min (<b>e</b>) 8 min; (<b>f</b>) 10 min. The change from the 1 m resolution of the inner domain to the 3 m resolution of the large domain can be noticed. The threshold for 500 µg/m<sup>3</sup> of concentration is in red, while that for 1 µg/m<sup>3</sup> is in orange and that for 0.001 µg/m<sup>3</sup> is in green.</p>
Full article ">Figure 9
<p>View of the three release locations and of the plumes, at the ground level, using tiled navigation. The threshold for 0.01 µg/m<sup>3</sup> of concentration is in green, while that for 1 µg/m<sup>3</sup> is in red (background satellite image © TerraMetrics, LLC—<a href="http://www.terrametrics.com" target="_blank">www.terrametrics.com</a> © Google 2021, accessed on 30 April 2021).</p>
Full article ">
15 pages, 1337 KiB  
Article
Impact of Industrial Air Pollution on Agricultural Production
by Wei Wei and Zanxin Wang
Atmosphere 2021, 12(5), 639; https://doi.org/10.3390/atmos12050639 - 18 May 2021
Cited by 20 | Viewed by 6336
Abstract
This paper aimed to study how industrial air pollution impacts crop yield by investigating the relationship between output and changes in factors. A translog production function was estimated in the context of stochastic frontier analysis using data collected from a field survey in [...] Read more.
This paper aimed to study how industrial air pollution impacts crop yield by investigating the relationship between output and changes in factors. A translog production function was estimated in the context of stochastic frontier analysis using data collected from a field survey in the case of corn. The interaction between the factors as well as the impact of industrial air pollution on the relationship between factors was analyzed using numerical simulation, followed by the estimation of economic losses of corn yield in the polluted area. Results show that industrial air pollution causes a decrease in crop yield for two reasons. First, industrial air pollution changes the output elasticities of production factors and reduces its absolute amount. Second, industrial air pollution causes the relationship between labor and capital, labor and chemicals, capital and seeds to change from substitutable to complementary; it also resulted in an opposite result for the relationship between capital and chemicals. The paper presents a new explanation of how industrial air pollution affects agricultural production from an economic perspective. Full article
Show Figures

Figure 1

Figure 1
<p>The output elasticity of K along with the changes of D in three models.</p>
Full article ">Figure 2
<p>Elasticity of substitution between factors in two groups.</p>
Full article ">Figure 3
<p>The economic loss of crop production.</p>
Full article ">
15 pages, 7732 KiB  
Article
Evaluation of a Flexible Single Ice Microphysics and a Gaussian Probability-Density-Function Macrophysics Scheme in a Single Column Model
by Jiabo Li, Xindong Peng, Xiaohan Li, Yanluan Lin and Wenchao Chu
Atmosphere 2021, 12(5), 638; https://doi.org/10.3390/atmos12050638 - 17 May 2021
Cited by 1 | Viewed by 2369
Abstract
Scale-aware parameterizations of subgrid scale physics are essentials for multiscale atmospheric modeling. A single-ice (SI) microphysics scheme and Gaussian probability-density-function (Gauss-PDF) macrophysics scheme were implemented in the single-column Global-to-Regional Integrated forecast System model (SGRIST) and they were tested using the Tropical Warm Pool-International [...] Read more.
Scale-aware parameterizations of subgrid scale physics are essentials for multiscale atmospheric modeling. A single-ice (SI) microphysics scheme and Gaussian probability-density-function (Gauss-PDF) macrophysics scheme were implemented in the single-column Global-to-Regional Integrated forecast System model (SGRIST) and they were tested using the Tropical Warm Pool-International Cloud Experiment (TWP-ICE) and the Atmospheric Radiation Measurement Southern Great Plains Experiment in 1997 (ARM97). Their performance was evaluated against observations and other reference schemes. The new schemes simulated reasonable precipitation with proper fluctuations and peaks, ice, and liquid water contents, especially in lower levels below 650 hPa during the wet period in the TWP-ICE. The root mean square error (RMSE) of the simulated cloud fraction was below 200 hPa was 0.10/0.08 in the wet/dry period, which showed an obvious improvement when compared to that, i.e., 0.11/0.11 of original scheme. Accumulated ice water content below the melting level decreased by 21.57% in the SI. The well-matched average liquid water content displayed between the new scheme and observations, which was two times larger than those with the referencing scheme. In the ARM97 simulations, the SI scheme produced considerable ice water content, especially when convection was active. Low-level cloud fraction and precipitation extremes were improved using the Gauss-PDF scheme, which displayed the RMSE of cloud fraction of 0.02, being only half of the original schemes. The study indicates that the SI and Gauss-PDF schemes are promising approaches to simplify the microphysics process and improve the low-level cloud modeling. Full article
(This article belongs to the Section Atmospheric Techniques, Instruments, and Modeling)
Show Figures

Figure 1

Figure 1
<p>Temporal-vertical plot of cloud fraction (left), ice water content (middle column) and liquid water content (right) in addition to the atmospheric temperature (solid, °C) in the MG-PBR (<b>d</b>–<b>f</b>), SI-PBR (<b>g</b>–<b>i</b>), MG-GPDF (<b>j</b>–<b>l</b>), and SI-GPDF (<b>m</b>–<b>o</b>) simulation of the TWP-ICE in comparison with the observations and retrievals (<b>a</b>–<b>c</b>). The melting level is marked with bold lines in the middle panels.</p>
Full article ">Figure 2
<p>Vertical profiles of cloud fraction and the deviation from the observation (<b>a</b>,<b>d</b>), ice water content (<b>b</b>,<b>e</b>) and liquid water content (<b>c</b>,<b>f</b>) averaged in the wet (18–25 January, top) and dry periods (26–30 January, bottom) in the four tests, MG-PBR (red), SI-PBR (cyan), MG-GPDF (orange), and SI-GPDF (brown) in comparison with the observations (black).</p>
Full article ">Figure 3
<p>Vertical profiles of microphysics processes for ice in the MG–PBR and SI–PBR experiments, sublimation (red), melting (cyan), deposition (orange), and auto-conversion of cloud ice to snow (brown), respectively, averaged in the wet period.</p>
Full article ">Figure 4
<p>(<b>a</b>) Total precipitation rate, (<b>b</b>) large-scale precipitation rate, (<b>c</b>) convective precipitation rate and (<b>d</b>) liquid cloud water path in the simulations MG-PBR (red), SI-PBR (cyan), MG-GPDF (orange), and SI-GPDF (brown) compared to the observation (black).</p>
Full article ">Figure 5
<p>Temporal-vertical plot of cloud fraction (left), ice water content (middle column) and liquid water content (right) in addition to the atmospheric temperature (solid, °C) in the MG-PBR (<b>d</b>–<b>f</b>), SI-PBR (<b>g</b>–<b>i</b>), MG-GPDF (<b>j</b>–<b>l</b>), and SI-GPDF (<b>m</b>,<b>n</b>,<b>o</b>) simulation of the ARM97 in comparison with the observations and retrievals (<b>a</b>–<b>c</b>). The melting level is marked with bold lines in the middle panels.</p>
Full article ">Figure 6
<p>Vertical profiles of cloud fraction and the deviation from the observation (<b>a</b>), ice water content (<b>b</b>) and liquid water content (<b>c</b>) averaged during 19 June and 17 July in the four tests, MG-PBR (red), SI-PBR (cyan), MG-GPDF (orange), and SI-GPDF (brown) in comparison with the observations (black).</p>
Full article ">Figure 7
<p>Temporal-vertical plot of mass-weighted ice falling speed in (<b>a</b>) the SI-PBR and (<b>b</b>) the SI-GPDF.</p>
Full article ">Figure 8
<p>(<b>a</b>) Total precipitation rate, (<b>b</b>) large-scale precipitation rate, (<b>c</b>) convective precipitation rate and (<b>d</b>) liquid cloud water path in the simulations MG-PBR (red), SI-PBR (cyan), MG-GPDF (orange), and SI-GPDF (brown) compared to the observation (black).</p>
Full article ">
11 pages, 2533 KiB  
Article
Emissions of Volatile Organic Compounds (VOCs) from an Open-Circuit Dry Cleaning Machine Using a Petroleum-Based Organic Solvent: Implications for Impacts on Air Quality
by Hyeonji Lee, Kyunghoon Kim, Yelim Choi and Daekeun Kim
Atmosphere 2021, 12(5), 637; https://doi.org/10.3390/atmos12050637 - 17 May 2021
Cited by 10 | Viewed by 3985
Abstract
Volatile organic compounds (VOCs) are known to play an important role in tropospheric chemistry, contributing to ozone and secondary organic aerosol (SOA) generation. Laundry facilities, using petroleum-based organic solvents, are one of the sources of VOCs emissions. However, little is known about the [...] Read more.
Volatile organic compounds (VOCs) are known to play an important role in tropospheric chemistry, contributing to ozone and secondary organic aerosol (SOA) generation. Laundry facilities, using petroleum-based organic solvents, are one of the sources of VOCs emissions. However, little is known about the significance of VOCs, emitted from laundry facilities, in the ozone and SOA generation. In this study, we characterized VOCs emission from a dry-cleaning process using petroleum-based organic solvents. We also assessed the impact of the VOCs on air quality by using photochemical ozone creation potential and secondary organic aerosol potential. Among 94 targeted compounds including toxic organic air pollutants and ozone precursors, 36 compounds were identified in the exhaust gas from a drying machine. The mass emitted from one cycle of drying operation (40 min) was the highest in decane (2.04 g/dry cleaning). Decane, nonane, and n-undecane were the three main contributors to ozone generation (more than 70% of the total generation). N-undecane, decane, and n-dodecane were the three main contributors to the SOA generation (more than 80% of the total generation). These results help to understand VOCs emission from laundry facilities and impacts on air quality. Full article
(This article belongs to the Section Air Quality)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Heat map of relative VOCs concentrations over the first 18 min of drying process. Principal component analysis was performed using R studio (version 1.3.1073) for aggregating analytes into different groups depending on the time point at which maximum concentration appeared. The relative concentrations are on the natural logarithmic scale.</p>
Full article ">Figure 2
<p>The concentrations of TVOCs in exhaust gas during drying processes with different operational conditions: (<b>a</b>) operating temperature in dry unit; (<b>b</b>) operating time in dry unit; (<b>c</b>) weight of laundry; (<b>d</b>); type of laundry cloth. For the operating condition in the dry-cleaning process, defaults were applied, except for the variable parameters. Default conditions were 3 kg of laundry weight (cotton); 40 min of drying time; 40 °C of operating temperature. Abbreviation: total volatile organic compound (TVOC).</p>
Full article ">Figure 3
<p>POCP<sub>weighted emission</sub> contributions (%) of dry-cleaning VOCs emitted from small-scale dry-cleaning operation.</p>
Full article ">Figure 4
<p>SOAP<sub>weighted emission</sub> contributions (%) of dry-cleaning VOCs emitted from small-scale dry-cleaning operation.</p>
Full article ">
24 pages, 6660 KiB  
Article
“Military Parade Blue Skies” in Beijing: Decisive Influence of Meteorological Factors on Transport Channel and Atmospheric Pollutant Concentration Level
by Shujian Yang, Yang Zhang, Jing Shang, Zhengqiang Li, Benjamin de Foy, James Jay Schauer and Yuanxun Zhang
Atmosphere 2021, 12(5), 636; https://doi.org/10.3390/atmos12050636 - 17 May 2021
Cited by 1 | Viewed by 2493
Abstract
The severity of high atmospheric pollution has been a major social problem in northern China. To improve the air quality in the Beijing–Tianjin–Hebei (BTH) region and guarantee a suitable environment during the military parade and other celebrating activities for the 70th anniversary of [...] Read more.
The severity of high atmospheric pollution has been a major social problem in northern China. To improve the air quality in the Beijing–Tianjin–Hebei (BTH) region and guarantee a suitable environment during the military parade and other celebrating activities for the 70th anniversary of the victory for anti-Fascist Warcraft in the year 2015, a series of strict air quality control policies were carried out. To analyze the reduction extents of PM2.5 and organic matter components during the control period and to examine the meteorological conditions in this region and their decisive influence on the air quality, PM2.5 samples were collected and the Lagrangian particle dispersion model FLEXPART was performed to calculate potential source locations within the BTH region. PM2.5, organic carbon (OC), elementary carbon (EC), and three species types were specifically analyzed. Although the results showed that PM2.5, OC, and EC reduced by 64.55%, 48.74%, and 60.75% during the control period, the air mass transport patterns showed great difference at certain periods, which altered the dominant transport direction of air mass and the potential source region of pollutants and organic matters. This alteration completely changed major atmospheric pollutants sources contribution and caused huge concentration changes. Parallel cases also showed that meteorological conditions could avoid massive atmospheric transported from a major emission source region to a receptor site. The meteorological conditions changed the main contribution source region in control and non-control periods and proved the air quality control measures were less necessary in some southern Hebei cities during special events periods. Full article
(This article belongs to the Special Issue Aerosol Pollution in Asia)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) shows the concentration of PM<sub>2.5</sub>, OC, and EC concentration from off-line samples and online monitoring system from RADI. The dashed lines are the time boundaries of three periods. (<b>b</b>) shows the average concentration of PM<sub>2.5</sub>, OC, and EC in 3 periods.</p>
Full article ">Figure 2
<p>(<b>a</b>) shows the correlation of PM<sub>2.5</sub> with OC and with EC. (<b>b</b>) is the correlation of OC and EC.</p>
Full article ">Figure 3
<p>The column plots for other atmospheric chemical compound concentration in three periods (CO, O<sub>3</sub>, NO<sub>2</sub>, SO<sub>2</sub>).</p>
Full article ">Figure 4
<p>Plots (<b>a</b>–<b>d</b>), respectively, show the correlation graph array between PM<sub>2.5</sub> and SO<sub>2</sub>, CO, NO<sub>2</sub>, and O<sub>3</sub>. The correlation coefficient R values are given in the subplots, respectively.</p>
Full article ">Figure 5
<p>Concentration of SOC and POC in 3 periods. The average concentrations are marked inside the columns.</p>
Full article ">Figure 6
<p>Average concentrations for aldehydes, alkanols, and sterols categories in three time periods. The concentration data are marked on the column bars.</p>
Full article ">Figure 7
<p>RTAs of the 3 time periods and the whole sampling period. Subplots (<b>a</b>–<b>c</b>) are the RTA fields for the “before”, “during”, and “after” periods. Subplot (<b>d</b>) is the whole sampling period. The red dot is our sampling site and the particle release point in FLEXPART model. The color bar represents the relative contribution percentage in each grid to the sampling site. From bottom to top side of the color bar, the cell contribution gets higher.</p>
Full article ">Figure 8
<p>Wind rose plots of the three separate periods and the total sampling period. The (<b>a</b>–<b>c</b>) subplots are for the before, during, and after periods. The total wind rose plot is the subplot (<b>d</b>). The same scale is used for wind direction frequencies in the three periods. The data used in this plot are hourly instantaneous data from the China National Air Pollution Control Monitoring System. Subplots (<b>e</b>,<b>f</b>) are the wind rose plots for two separate PM2.5 high concentration periods: 16–17 August, 16–17 September.</p>
Full article ">Figure 9
<p>CFA fields plots of PM<sub>2.5</sub>, OC, and EC in 3 time periods. Subplots (<b>a</b>–<b>d</b>) are CFAs of PM<sub>2.5</sub> in 3 time periods and the whole sampling period, respectively. The red dot is the sampling site and particle release point in FLEXPART. The color bar only represents the relative contribution percentage in each grid. Beijing city, Tianjin city, and Hebei Province are marked in the maps. From bottom to top side of the color bar, the cell contribution gets higher.</p>
Full article ">Figure 10
<p>RTAs of four time periods (before, during, after, total) in 2014 and 2016. Subplots (<b>a</b>,<b>c</b>,<b>e</b>,<b>g</b>) are for the year of 2014 and (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) are for the year of 2016. From bottom to top side of the color bar, the cell contribution gets higher.</p>
Full article ">Figure 11
<p>PM<sub>2.5</sub> CFAs comparison of 2014 and 2016 during the same sampling periods. Subplots (<b>a</b>,<b>c</b>,<b>e</b>,<b>g</b>) are the CFAs of the four periods in 2014. Subplots (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) are the CFAs of four periods in 2016. Other configurations and set-ups are the same in the MATLAB calculation scripts. From bottom to top side of the color bar, the cell contribution gets higher.</p>
Full article ">Figure 11 Cont.
<p>PM<sub>2.5</sub> CFAs comparison of 2014 and 2016 during the same sampling periods. Subplots (<b>a</b>,<b>c</b>,<b>e</b>,<b>g</b>) are the CFAs of the four periods in 2014. Subplots (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) are the CFAs of four periods in 2016. Other configurations and set-ups are the same in the MATLAB calculation scripts. From bottom to top side of the color bar, the cell contribution gets higher.</p>
Full article ">Figure 12
<p>PM<sub>2.5</sub> average concentration in three periods in 2014, 2015, and 2016. The averaged time periods in 2014 and 2016 are both in consistent with 2015 as demonstrated in previous sessions.</p>
Full article ">
16 pages, 1993 KiB  
Review
The Research Progress of the Influence of Agricultural Activities on Atmospheric Environment in Recent Ten Years: A Review
by Pengxiang Ge, Mindong Chen, Yan Cui and Dongyang Nie
Atmosphere 2021, 12(5), 635; https://doi.org/10.3390/atmos12050635 - 17 May 2021
Cited by 11 | Viewed by 4932
Abstract
In recent years, the industrial emission of air pollution has been reduced via a series of measures. However, with the rapid development of modern agriculture, air pollution caused by agricultural activities is becoming more and more serious. Agricultural activities can generate a large [...] Read more.
In recent years, the industrial emission of air pollution has been reduced via a series of measures. However, with the rapid development of modern agriculture, air pollution caused by agricultural activities is becoming more and more serious. Agricultural activities can generate a large amount of air pollutants, such as ammonia, methane, nitrogen oxides, volatile organic compounds, and persistent organic pollutants, the sources of which mainly include farmland fertilization, livestock breeding, pesticide use, agricultural residue burning, agricultural machinery, and agricultural irrigation. Greenhouse gases emitted by agricultural activities can affect regional climate change, while atmospheric particulates and persistent organic pollutants can even seriously harm the health of surrounding residents. With the increasing threat of agricultural air pollution, more and more relevant studies have been carried out, as well as some recommendations for reducing emissions. The emissions of ammonia and greenhouse gases can be significantly reduced by adopting reasonable fertilization methods, scientific soil management, and advanced manure treatment systems. Regarding pesticide use and agricultural residues burning, emission reduction are more dependent on the restriction and support of government regulations, such as banning certain pesticides, prohibiting open burning of straw, and supporting the recycling and reuse of residues. This review, summarizing the relevant research in the past decade, discusses the current situation, health effects, and emission reduction measures of agricultural air pollutants from different sources, in order to provide some help for follow-up research. Full article
(This article belongs to the Special Issue Agricultural Pollutants in the Atmosphere)
Show Figures

Figure 1

Figure 1
<p>The spatial distribution of NH<sub>3</sub> emissions in China [<a href="#B29-atmosphere-12-00635" class="html-bibr">29</a>] and the spatial distribution of agricultural and synthetic fertilizer application sources in the Beijing-Tianjin-Hebei region [<a href="#B72-atmosphere-12-00635" class="html-bibr">72</a>].</p>
Full article ">Figure 2
<p>GHG emissions during beef production in Saskatchewan, Canada [<a href="#B32-atmosphere-12-00635" class="html-bibr">32</a>].</p>
Full article ">Figure 3
<p>Pesticides can influence the environment and human health in different ways.</p>
Full article ">Figure 4
<p>(<b>a</b>) Annual PM<sub>2.5</sub> emissions from total anthropogenic sources and biomass burning in some years in China and (<b>b</b>) annual PM<sub>2.5</sub> emissions from biomass burning in different regions of China [<a href="#B41-atmosphere-12-00635" class="html-bibr">41</a>].</p>
Full article ">Figure 5
<p>Ratio of emissions from agricultural machinery and road vehicles during busy season (<b>a</b>) Spring and (<b>b</b>) Autumn in the Yangtze River Delta of China [<a href="#B45-atmosphere-12-00635" class="html-bibr">45</a>].</p>
Full article ">
26 pages, 38633 KiB  
Article
The 3D Neural Network for Improving Radar-Rainfall Estimation in Monsoon Climate
by Nurulhani Roslan, Mohd Nadzri Md Reba, Syarawi M. H. Sharoni and Mohammad Shawkat Hossain
Atmosphere 2021, 12(5), 634; https://doi.org/10.3390/atmos12050634 - 17 May 2021
Cited by 4 | Viewed by 3056
Abstract
The reflectivity (Z)—rain rate (R) model has not been tested on single polarization radar for estimating monsoon rainfall in Southeast Asia, despite its widespread use for estimating heterogeneous rainfall. The artificial neural network (ANN) regression has been applied to the radar [...] Read more.
The reflectivity (Z)—rain rate (R) model has not been tested on single polarization radar for estimating monsoon rainfall in Southeast Asia, despite its widespread use for estimating heterogeneous rainfall. The artificial neural network (ANN) regression has been applied to the radar reflectivity data to estimate monsoon rainfall using parametric Z-R models. The 10-min reflectivity data recorded in Kota Bahru radar station (in Malaysia) and hourly rain record in nearby 58 gauge stations during 2013–2015 were used. The three-dimensional nearest neighbor interpolation with altitude correction was applied for pixel matching. The non-linear Levenberg Marquardt (LM) regression, integrated with ANN regression minimized the spatiotemporal variability of the proposed Z-R model. Results showed an improvement in the statistical indicator, when LM and ANN overestimated (6.6%) and underestimated (4.4%), respectively, the mean total rainfall. For all rainfall categories, the ANN model has a positive efficiency ratio of >0.2. Full article
(This article belongs to the Special Issue Weather Radar in Rainfall Estimation)
Show Figures

Figure 1

Figure 1
<p>Weather radar map, showing locations of rain gauges in Kota Bahru, Peninsular Malaysia.</p>
Full article ">Figure 2
<p>Illustration of single radar beam observing one target (A) at the atmosphere is a function of the range (<span class="html-italic">r</span>), elevation angles (<span class="html-italic">θ</span>) and azimuth angles (<span class="html-italic">ϕ</span>) to be converted into Cartesian position in a function of the latitude (<span class="html-italic">x</span>), longitude (<span class="html-italic">y</span>), and height (<span class="html-italic">z</span>).</p>
Full article ">Figure 3
<p>The error-bar plots illustrating reflectivity measured in different gauge locations interpolated by the nearest distance (<b>a</b>) and the 3D interpolation (<b>b</b>), and the performance of (<b>c</b>) measured in terms of residual mean square error (RMSE). Horizontal bars indicate the standard error of the means.</p>
Full article ">Figure 4
<p>Time series plot for the corresponding (<b>a</b>) reflectivity (dBZ) and (<b>b</b>) rainfall intensity (mm/h) from 2013 to 2015. Higher reflectivity found up to 60 dBZ and high rainfall intensity was recorded at 150 mm/h. Yellow line represents the time average of each parameter.</p>
Full article ">Figure 5
<p>Trend plot of α (solid line with circle mark) and β (dash line with cross mark) at each month from January 2013 to March 2015.</p>
Full article ">Figure 6
<p>The Taylor diagram illustrating a statistical comparison of four different hourly radar rainfall estimates with rain gauge (RG) observation under different rainfall intensities: (<b>a</b>) low, (<b>b</b>) moderate, (<b>c</b>) high, and (<b>d</b>) all intensity. The azimuthal angle denotes correlation; the radial distance is standard deviation; and the semicircles centered at the RG represent residual mean square error (RMSE).</p>
Full article ">Figure 7
<p>Bar plot of bias and NSE for the estimated radar rainfall model.</p>
Full article ">Figure 8
<p>Plot of R<sup>2</sup> and G/R of radar rainfall estimates on a monthly basis for all rainfall intensity for (<b>a</b>) LM, (<b>b</b>) MP, (<b>c</b>) ROS, and (<b>d</b>) ANN models.</p>
Full article ">Figure 9
<p>Plot of RMSE for monthly radar rainfall at (<b>a</b>) all, (<b>b</b>) low, (<b>c</b>) medium, and (<b>d</b>) high rain intensities.</p>
Full article ">Figure 10
<p>Scatter plots of radar rainfall estimates based on ANN versus rainfall measurements from the gauges. (<b>a</b>) training, (<b>b</b>) testing and (<b>c</b>) all regressions of the reflectivity for all rainfall intensity. The color lines are fitted line and the dashed lines are 95% confidence bounds of the fitted line.</p>
Full article ">Figure 11
<p>Scatter plots of radar rainfall estimates based on ANN versus rainfall measurements from the gauges at (<b>a</b>) low, (<b>b</b>) medium, and (<b>c</b>) high intensity rainfall. The straight line is the regression line.</p>
Full article ">Figure 12
<p>Plot of the rain gauge measurement and estimated rainfall from the radar rainfall on hourly basis for (<b>a</b>) 16 to 19 December 2014 (96 h) and (<b>b</b>) 20 to 24 December 2014 (102 h).</p>
Full article ">Figure 13
<p>CAPPI map on daily basis generated from the LM model for (<b>a</b>) 16 December 2014, (<b>b</b>) 17 December 2014, (<b>c</b>) 18 December 2014, and (<b>d</b>) 19 December 2014 for KB weather radar station.</p>
Full article ">Figure 14
<p>CAPPI map on daily basis generated from the ANN model from the KB weather radar for similar date as <a href="#atmosphere-12-00634-f013" class="html-fig">Figure 13</a>. (<b>a</b>) 16 December 2014, (<b>b</b>) 17 December 2014, (<b>c</b>) 18 December 2014, and (<b>d</b>) 19 December 2014.</p>
Full article ">
23 pages, 93677 KiB  
Article
Exposure Assessment of Climate Extremes over the Europe–Mediterranean Region
by Mehmet Barış Kelebek, Fulden Batibeniz and Barış Önol
Atmosphere 2021, 12(5), 633; https://doi.org/10.3390/atmos12050633 - 17 May 2021
Cited by 19 | Viewed by 4914
Abstract
The use of a compact set of climate change indexes enhances our understanding of the combined impacts of extreme climatic conditions. In this study, we developed the modified Climate Extremes Index (mCEI) to obtain unified information about different types of extremes. For this [...] Read more.
The use of a compact set of climate change indexes enhances our understanding of the combined impacts of extreme climatic conditions. In this study, we developed the modified Climate Extremes Index (mCEI) to obtain unified information about different types of extremes. For this purpose, we calculated 10 different climate change indexes considering the temperature extremes, extreme precipitation, and moisture surplus and drought over the Europe–Mediterranean (EURO–MED) region for the 1979–2016 period. As a holistic approach, mCEI provides spatiotemporal information, and the high-resolution grid-based data allow us to accomplish detailed country-based and city-based analyses. The analyses indicate that warm temperature extremes rise significantly over the EURO–MED region at a rate of 1.9% decade−1, whereas the cold temperature extremes decrease. Extreme drought has a significant increasing trend of 3.8% decade−1. Although there are regional differences, extreme precipitation indexes have a significant increasing tendency. According to the mCEI, the major hotspots for the combined extremes are the Mediterranean coasts, the Balkan countries, Eastern Europe, Iceland, western Russia, western Turkey, and western Iraq. The decadal changes of mCEI for these regions are in the range of 3–5% decade−1. The city-scale analysis based on urbanized locations reveals that Fes (Morocco), Izmir (Turkey), Marseille and Aix-en-Provence (France), and Tel Aviv (Israel) have the highest increasing trend of mCEI, which is greater than 3.5% decade−1. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>mCEI is the arithmetic average of five indicators that measure the annual temporal intensities of extreme climatic conditions.</p>
Full article ">Figure 2
<p>Spatial distribution of the thresholds for the extreme indexes.</p>
Full article ">Figure 3
<p>The climatological means (1979–2016) of the extreme indexes.</p>
Full article ">Figure 4
<p>Spatial distribution of the trends in the extreme indexes.</p>
Full article ">Figure 5
<p>Spatial averages of the extreme indexes over the EURO–MED region.</p>
Full article ">Figure 6
<p>Zonal average of the extreme indexes over the EURO–MED region.</p>
Full article ">Figure 7
<p>(<b>a</b>) Reference period (1981–2010) climatology, and (<b>b</b>) decadal anomalies of mCEI.</p>
Full article ">Figure 8
<p>(<b>a</b>) Zonally averaged mCEI anomaly, (<b>b</b>) decadal trend map of mCEI (% decade<sup>−1</sup>), and (<b>c</b>) confidence level of the trend.</p>
Full article ">Figure 9
<p>Spatially-average of mCEI (<b>upper</b>) and annual anomalies (<b>bottom</b>) over the EURO–MED region.</p>
Full article ">Figure 10
<p>Extreme indexes, mCEI, and mCEI anomaly for 2003.</p>
Full article ">Figure 11
<p>Extreme indexes, mCEI, and mCEI anomaly for 2010.</p>
Full article ">Figure 12
<p>Country-based annual and climatological means, and decadal trends of mCEI. Boldface country name indicates a statistically significant trend at 95% confidence level.</p>
Full article ">Figure 13
<p>mCEI climatology over the urban settlements of the EURO–MED region.</p>
Full article ">Figure 14
<p>Spatial averages of mCEI over the urban settlements of the EURO–MED region.</p>
Full article ">Figure 15
<p>City-based decadal trends and climatological means of mCEI over the major urban settlements of the EURO–MED region. Boldface city name indicates a statistically significant trend at 95% confidence level.</p>
Full article ">
14 pages, 4195 KiB  
Article
Characteristics of Summer Hourly Extreme Precipitation Events and Its Local Environmental Influencing Factors in Beijing under Urbanization Background
by Zuofang Zheng, Guirong Xu and Hua Gao
Atmosphere 2021, 12(5), 632; https://doi.org/10.3390/atmos12050632 - 16 May 2021
Cited by 10 | Viewed by 2634
Abstract
Studies on urban extreme precipitation and its influencing factors are significant for prevention and reduction of meteorological disasters; however, few studies focus on hourly extreme precipitation (HEP) events due to the limitation of observation. By using the summer hourly precipitation data in Beijing [...] Read more.
Studies on urban extreme precipitation and its influencing factors are significant for prevention and reduction of meteorological disasters; however, few studies focus on hourly extreme precipitation (HEP) events due to the limitation of observation. By using the summer hourly precipitation data in Beijing from 1980 to 2020, the spatial distribution and temporal variation of HEP as well as its local environmental influencing factors are investigated. It is found that both summer precipitation amount and frequency of HEP are affected by topography, with high values in windward slope area. The summer precipitation amount of HEP is 160–200 mm, accounting for 42–47% of the annual summer precipitation amount, while the frequency proportion of HEP is only 5.5–6.0%. Although the summer precipitation amount and frequency in Beijing both decrease in the past 41 years, those for HEP present an opposite trend mainly due to the increasing HEP since 2003 and this is a phenomenon worthy of attention. A similar bimodal pattern in diurnal variation is found for the summer precipitation amount and frequency of HEP, with two peaks in 19–22 LT and 01–05 LT, respectively, indicating that HEPs are more concentrated in the evening and early morning especially in urban area. Moreover, the urbanization process of Beijing is consistent with the change trend of HEP, implying that the stronger the urban heat island intensity (UHI), the higher the probability of HEP. Furthermore, the convergence lines of terrain are also conducive to local heavy rainfall, and lower tropospheric stability (LTS) and convective available potential energy (CAPE) as well as aerosols may also affect HEP in Beijing. Full article
(This article belongs to the Special Issue Hazards, Urbanization, and Climate Change)
Show Figures

Figure 1

Figure 1
<p>Terrain height (unit: m, gray filled) and urban impervious surface (red area) in Beijing as well as HEP thresholds (unit: mm, the numbers) of 20 stations.</p>
Full article ">Figure 2
<p>Spatial distribution of summer precipitation in Beijing from 1980 to 2020: (<b>a</b>) mean annual summer precipita-tion amount (unit: mm, color filled) and frequency (unit: h, isoline); (<b>b</b>) mean summer precipitation amount of HEP (unit: mm, color filled) and its proportion in mean annual summer precipitation amount (unit: %, isoline); (<b>c</b>) mean summer precipitation frequency of HEP (unit: h, color filled) and its proportion in mean annual summer precipitation frequency (unit: %, isoline); (<b>d</b>) mean summer precipitation intensity (unit: mm/h, color filled) and that for HEP (unit: mm/h, isoline). The "A", "B" and "C" in the figure represent the windward slope of northeast mountain area, central urban area and windward slope of southwest mountain area in Beijing, respectively.</p>
Full article ">Figure 2 Cont.
<p>Spatial distribution of summer precipitation in Beijing from 1980 to 2020: (<b>a</b>) mean annual summer precipita-tion amount (unit: mm, color filled) and frequency (unit: h, isoline); (<b>b</b>) mean summer precipitation amount of HEP (unit: mm, color filled) and its proportion in mean annual summer precipitation amount (unit: %, isoline); (<b>c</b>) mean summer precipitation frequency of HEP (unit: h, color filled) and its proportion in mean annual summer precipitation frequency (unit: %, isoline); (<b>d</b>) mean summer precipitation intensity (unit: mm/h, color filled) and that for HEP (unit: mm/h, isoline). The "A", "B" and "C" in the figure represent the windward slope of northeast mountain area, central urban area and windward slope of southwest mountain area in Beijing, respectively.</p>
Full article ">Figure 3
<p>Interannual variations of the summer maximum hourly precipitation in urban area (red) and its difference (blue) between urban and suburban areas in Beijing during 1980–2020. The dash line is the 5-point moving average series.</p>
Full article ">Figure 4
<p>Interannual variations of summer precipitation in urban area in Beijing: (<b>a</b>) precipitation amount (left) and proportion of HEP precipitation amount in annual summer precipitation amount (right); (<b>b</b>) precipitation frequency (left) and proportion of HEP precipitation frequency in annual summer precipitation frequency (right); (<b>c</b>) precipitation intensity anomaly; (<b>d</b>) Mann–Kendall test on precipitation amount of HEP, and the thick straight lines indicate the values at significance level of 0.05.The anomaly values in the figure are the differences between the meteorological elements of the current year and the climatic mean values during 1981–2010, and the dash line is the 5-point moving average series.</p>
Full article ">Figure 5
<p>Diurnal variations of summer precipitation in Beijing: (<b>a</b>) precipitation amount (left) and proportion of HEP precipitation amount in annual summer precipitation amount (right); (<b>b</b>) precipitation frequency (left) and proportion of HEP precipitation frequency in annual summer precipitation frequency (right); (<b>c</b>) change trends for precipitation frequency (left) and amount (right) of HEP. The blue curve is a polynomial fitting data series.</p>
Full article ">Figure 6
<p>Annual changes of built-up area and UHI anomaly in Beijing during 1980–2020.</p>
Full article ">Figure 7
<p>The spatial distributions of (<b>a</b>) UHI and (<b>b</b>) haze days in summer in Beijing during 1980–2020 and their correlation coefficients with HEP frequency proportion. The UHI (unit: °C) and haze days are color filled, while blue isoline indicates the correlation coefficient. The wind shafts at the stations on the left panel indicate the prevailing surface wind field in summer in Beijing, and the thick dotted line indicates the wind convergence line caused by terrain.</p>
Full article ">Figure 8
<p>Diurnal variations of (<b>a</b>) sensible heat flux (QE) and latent heat flux (QH) in summer in Beijing during 2015–2016, and (<b>b</b>) those of lower tropospheric stability (LTS) and convective available potential energy (CAPE) in summer in Beijing during 2008–2018.</p>
Full article ">
14 pages, 1177 KiB  
Article
Pesticide Inhalation Exposure of Applicators and Bystanders Using Conventional and Innovative Cropping Systems in the Valencian Region, Spain
by Esther Fuentes, Antonio López, María Ibáñez, Vicent Yusà, Amalia Muñoz, Teresa Vera, Esther Borrás, Héctor Calvete-Sogo and Clara Coscollà
Atmosphere 2021, 12(5), 631; https://doi.org/10.3390/atmos12050631 - 15 May 2021
Cited by 4 | Viewed by 2934
Abstract
This paper provides scientific results from a European LIFE project carried out in the Valencian region of Spain during the 2017 to 2018 time frame. In 2018, more than 60,000 tons of pesticides were commercialized in Spain, with approximately 15% destined for Valencian [...] Read more.
This paper provides scientific results from a European LIFE project carried out in the Valencian region of Spain during the 2017 to 2018 time frame. In 2018, more than 60,000 tons of pesticides were commercialized in Spain, with approximately 15% destined for Valencian crops. In order to improve the air quality in the agricultural areas of this region, an innovative cropping system based on irrigation was developed and compared to conventional treatments based on hand-spray and turbo application. After applying conventional treatments to five types of crops (citrus, persimmon, nectarine, watermelon, and other stone fruits), a total of 13 active substances were detected in the air. The same active substances were applied to crops using the novel irrigation system, and no pesticide was detected in the air. Moreover, applicator and bystander populations in the region were assessed for their risk of inhalation exposure to pesticides, and no risk was found when either of the techniques, the innovative and the conventional agricultural one, were applied. Full article
(This article belongs to the Special Issue Agricultural Pollutants in the Atmosphere)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Perfect Life Project and (<b>b</b>) Irrilife project.</p>
Full article ">Figure 2
<p>Main injection point (circled in purple) and secondary injection points (red) in the IRRILIFE drip irrigation network. The drip irrigation network’s total surface area was &gt;100 has. The surface of crops treated during the Irrilife project (green square) was &gt;5.7 has. Location: L’Alcudia, Valencia, Spain.</p>
Full article ">Figure 3
<p>PUF+XAD2+PUF sandwich for the gas phase collection.</p>
Full article ">
18 pages, 7310 KiB  
Article
Longer Time-Scale Variability of Atmospheric Vertical Motion over the Tibetan Plateau and North Pacific and the Climate in East Asia
by Rongxiang Tian, Yaoming Ma, Weiqiang Ma, Xiuyi Zhao and Duo Zha
Atmosphere 2021, 12(5), 630; https://doi.org/10.3390/atmos12050630 - 15 May 2021
Cited by 2 | Viewed by 2589
Abstract
The vertical motion of air is closely related to the amount of precipitation that falls in a particular region. The Tibetan Plateau and the North Pacific are important determinants of the East Asian climate. We use climate diagnosis and statistical analysis to study [...] Read more.
The vertical motion of air is closely related to the amount of precipitation that falls in a particular region. The Tibetan Plateau and the North Pacific are important determinants of the East Asian climate. We use climate diagnosis and statistical analysis to study the vertical motion of the air over the North Pacific and Tibetan Plateau and the relationship between the vertical motion of air over them and the climate in East Asia. Here we show that there is a downward movement of air over the Tibetan Plateau during the winter, with a maximum velocity of downward movement at 500 hPa, whereas there is an upward movement of air with a maximum velocity of upward movement at 600 hPa during the summer. Precipitation in East Asia has a significant negative correlation (The correlation coefficient exceeds −0.463 and confidence level is greater than 99%) with the vertical motion of air over the Tibetan Plateau and the North Pacific during both the winter and summer. There is also a negative correlation of precipitation in the region south of the Yangtze River with the vertical motion of air over the Tibetan Plateau in winter, whereas the area of negative correlation to the vertical motion of air over the North Pacific in winter is located to the east of the Tibetan Plateau and the Yangtze–Huaihe river basin. The research results provide a climatic framework for the vertical motion of air over both the Tibetan Plateau and the North Pacific. Full article
(This article belongs to the Section Climatology)
Show Figures

Figure 1

Figure 1
<p>Distribution of the vertical velocity of air over the Tibetan Plateau. (<b>a</b>) Vertical velocity at 500 hPa in spring (March–May), summer (June–August), autumn (September–November) and winter (December–February). (<b>b</b>) Profiles of vertical velocity centers in summer (33° N, 97° E) and winter (100° E, 30° N). (<b>c</b>) Spatial distribution of the primary EOF-analyzed mode for vertical velocities at 200 and 500 hPa in summer and winter. (<b>d</b>) Temporal variation of the primary mode of the vertical velocities at 200 and 500 hPa in summer and winter. <span class="html-italic">PC</span> denotes principal component and gray areas denote large fluctuations. (<b>e</b>) Wavelet power spectrum of temporal coefficients of the primary mode of the vertical velocity at 500 hPa in both summer and winter. The red line delineates the cone of influence and the yellow areas show confidence levels &gt;95%.</p>
Full article ">Figure 2
<p>Monthly means of vertical velocity over the Tibetan Plateau at 500 hPa from 1981 to 2010.</p>
Full article ">Figure 3
<p>Distribution of the vertical velocity of air over the North Pacific. (<b>a</b>) Vertical velocity at 500 hPa in spring (March–May), summer (June–August), autumn (September–November) and winter (December–February). (<b>b</b>) Primary EOF-analyzed mode of vertical velocities at 500 and 850 hPa in summer and winter. (<b>c</b>) Temporal variation in the primary mode of the vertical velocities at 500 and 850 hPa in summer and winter. PC denotes principal component and gray areas denote large fluctuations. (<b>d</b>) Wavelet power spectrum of the temporal coefficients of the primary mode of vertical velocity at 500 hPa in both summer and winter. The red line delineates the cone of influence and the yellow areas show confidence levels &gt;95%.</p>
Full article ">Figure 4
<p>Correlation analysis between the surface meteorological variables and the vertical motion of air for (<b>a</b>) the Tibetan Plateau and (<b>b</b>) the North Pacific. Red shades denote a positive correlation and blue shades denote a negative correlation. Note: the South China Sea has not been marked due to layout reasons.</p>
Full article ">Figure 5
<p>Wind velocity fields. (<b>a</b>) Mean meridional circulation at 90° E in January. (<b>b</b>) Vertical velocity field at 500 hPa in January. (<b>c</b>) Mean zonal circulation at 30° N in January. The grey area in part (<b>a</b>) represents the Tibetan Plateau area, and the white rectangle (<b>b</b>) represents the downdraft area and the red triangle (<b>c</b>) represents the intersection of the updraft and the downdraft.</p>
Full article ">
22 pages, 65908 KiB  
Article
Novel Integrated and Optimal Control of Indoor Environmental Devices for Thermal Comfort Using Double Deep Q-Network
by Sun-Ho Kim, Young-Ran Yoon, Jeong-Won Kim and Hyeun-Jun Moon
Atmosphere 2021, 12(5), 629; https://doi.org/10.3390/atmos12050629 - 14 May 2021
Cited by 7 | Viewed by 2994
Abstract
Maintaining a pleasant indoor environment with low energy consumption is important for healthy and comfortable living in buildings. In previous studies, we proposed the integrated comfort control (ICC) algorithm, which integrates several indoor environmental control devices, including an air conditioner, a ventilation system, [...] Read more.
Maintaining a pleasant indoor environment with low energy consumption is important for healthy and comfortable living in buildings. In previous studies, we proposed the integrated comfort control (ICC) algorithm, which integrates several indoor environmental control devices, including an air conditioner, a ventilation system, and a humidifier. The ICC algorithm is operated by simple on/off control to maintain indoor temperature and relative humidity within a defined comfort range. This simple control method can cause inefficient building operation because it does not reflect the changes in indoor–outdoor environmental conditions and the status of the control devices. To overcome this limitation, we suggest the artificial intelligence integrated comfort control (AI2CC) algorithm using a double deep Q-network(DDQN), which uses a data-driven approach to find the optimal control of several environmental control devices to maintain thermal comfort with low energy consumption. The suggested AI2CC showed a good ability to learn how to operate devices optimally to improve indoor thermal comfort while reducing energy consumption. Compared to the previous approach (ICC), the AI2CC reduced energy consumption by 14.8%, increased the comfort ratio by 6.4%, and decreased the time to reach the comfort zone by 54.1 min. Full article
(This article belongs to the Special Issue Zero Energy Building and Indoor Thermal)
Show Figures

Figure 1

Figure 1
<p>Comfort range in a cooling season on a psychrometric chart.</p>
Full article ">Figure 2
<p>Control modes based on indoor air initial state.</p>
Full article ">Figure 3
<p>Flowchart of ICC in the cooling season.</p>
Full article ">Figure 4
<p>Reinforcement learning model.</p>
Full article ">Figure 5
<p>Floor plan of the BICT and indoor environmental control devices.</p>
Full article ">Figure 6
<p>Co-simulation for AI2CC.</p>
Full article ">Figure 7
<p>Convergence of AI2CC.</p>
Full article ">Figure 8
<p>Temperature and relative humidity during the uncontrolled period (8 June).</p>
Full article ">Figure 9
<p>Energy consumption of AI2CC per episode.</p>
Full article ">Figure 10
<p>Comfort ratio of AI2CC per episode.</p>
Full article ">Figure 11
<p>Indoor status profile of AI2CC per episode.</p>
Full article ">Figure 12
<p>Indoor temperature and device states using ICC.</p>
Full article ">Figure 13
<p>Indoor relative humidity and humidifier state using ICC.</p>
Full article ">Figure 14
<p>Indoor temperature and device states using AI2CC (episode 1984).</p>
Full article ">Figure 15
<p>Indoor relative humidity and humidifier state using AI2CC (episode 1984).</p>
Full article ">Figure 16
<p>Cooling load of an air conditioner (<b>upper</b>: ICC, <b>lower</b>: AI2CC).</p>
Full article ">
19 pages, 3506 KiB  
Article
Overview of Low-Level Wind Shear Characteristics over Chinese Mainland
by Caiyan Lin, Kaijun Zhang, Xintao Chen, Sheng Liang, Junjie Wu and Wei Zhang
Atmosphere 2021, 12(5), 628; https://doi.org/10.3390/atmos12050628 - 14 May 2021
Cited by 14 | Viewed by 3649
Abstract
The characteristics of low-level wind shear (LLWS) over the Chinese mainland were investigated using reports from pilots, air traffic controllers and the number of arriving/departing flights from 2016 to 2020. A preliminary analysis of the impact of the flights on the LLWS reports [...] Read more.
The characteristics of low-level wind shear (LLWS) over the Chinese mainland were investigated using reports from pilots, air traffic controllers and the number of arriving/departing flights from 2016 to 2020. A preliminary analysis of the impact of the flights on the LLWS reports was carried out, and the cause of LLWS was also investigated. LLWS reports from most airports indicate that LLWS is most likely to occur within 600 m AGL with a higher density below 300 m, but with some exceptions, as wind shear is reported at higher altitudes at some airports. Airports with a high frequency of LLWS reports are almost all located in or around regions with complex topography and in regions with prevailing weather conditions favorable to LLWS. The variation in overall LLWS reports displays a steady increase from 2016 to 2019 and a decrease from 2019 to 2020, consistent with the trend in the number of flights, but with no evidently similar trends for individual airports. Seasonal variations in LLWS reports are observed and demonstrate no notable impact caused by the number of flights, implying that the main cause is the monthly variation of weather conditions. Diurnal variation is also evident and largely associated with the variation in number of flights during the busy period in addition to weather conditions, such as common strong winds, in the afternoon. Full article
(This article belongs to the Special Issue Weather and Aviation Safety)
Show Figures

Figure 1

Figure 1
<p>Location of the 115 airports with LLWS reports in Chinese mainland. The size and color of the dots indicate the number of LLWS reports at each airport from 2016 to 2020.</p>
Full article ">Figure 2
<p>The ratio of number of LLWS reports to flights handled at the 20 airports.</p>
Full article ">Figure 3
<p>Overall vertical distribution of wind shear reports during 2016–2020 at or close to airports in Chinese mainland.</p>
Full article ">Figure 4
<p>Overall cumulative frequency distribution of LLWS at different heights.</p>
Full article ">Figure 5
<p>Cumulative frequency distribution of total wind shear events (overall) and that at the 20 airports within different heights.</p>
Full article ">Figure 6
<p>Variation of total LLWS reports and flights in Chinese mainland from 2016 to 2020.</p>
Full article ">Figure 7
<p>As in <a href="#atmosphere-12-00628-f006" class="html-fig">Figure 6</a>, but for the 20 airports.</p>
Full article ">Figure 8
<p>As in <a href="#atmosphere-12-00628-f006" class="html-fig">Figure 6</a>, but for monthly frequency distribution.</p>
Full article ">Figure 9
<p>Cumulative monthly variation in LLWS reports at the 20 airports.</p>
Full article ">Figure 10
<p>Cumulative seasonal variation of LLWS reports at the 20 airports.</p>
Full article ">Figure 11
<p>Overall diurnal variation of total LLWS reports.</p>
Full article ">Figure 12
<p>As in <a href="#atmosphere-12-00628-f011" class="html-fig">Figure 11</a>, but for the 20 airports.</p>
Full article ">
19 pages, 4125 KiB  
Article
Spatial Characteristics of Precipitation in the Greater Sydney Metropolitan Area as Revealed by the Daily Precipitation Concentration Index
by Kevin K. W. Cheung, Aliakbar A. Rasuly, Fei Ji and Lisa T.-C. Chang
Atmosphere 2021, 12(5), 627; https://doi.org/10.3390/atmos12050627 - 13 May 2021
Cited by 5 | Viewed by 3461
Abstract
In this study; the spatial distribution of the Daily Precipitation Concentration Index (DPCI) has been analyzed inside the Greater Sydney Metropolitan Area (GSMA). Accordingly, the rainfall database from the Australian Bureau of Meteorology archive was utilized after comprehensive quality control. The compiled data [...] Read more.
In this study; the spatial distribution of the Daily Precipitation Concentration Index (DPCI) has been analyzed inside the Greater Sydney Metropolitan Area (GSMA). Accordingly, the rainfall database from the Australian Bureau of Meteorology archive was utilized after comprehensive quality control. The compiled data contains a set of 41 rainfall stations indicating consistent daily precipitation series from 1950 to 2015. In the analysis of the DPCI across GSMA the techniques of Moran’s Spatial Autocorrelation has been applied. In addition, a cross-covariance method was applied to assess the spatial interdependency between vector-based datasets after performing an Ordinary Kriging interpolation. The results identify four well-recognized intense rainfall development zones: the south coast and topographic areas of the Illawarra district characterized by Tasman Sea coastal regions with DPCI values ranging from 0.61 to 0.63, the western highlands of the Blue Mountains, with values between 0.60 and 0.62, the inland regions, with lowest rainfall concentrations between 0.55 and 0.59, and lastly the districts located inside the GSMA with DPCI ranging 0.60 to 0.61. Such spatial distribution has revealed the rainstorm and severe thunderstorm activity in the area. This study applies the present models to identify the nature and mechanisms underlying the distribution of torrential rains over space within the metropolis of Sydney, and to monitor any changes in the spatial pattern under the warming climate. Full article
Show Figures

Figure 1

Figure 1
<p>The location map of the GSMA within New South Wales of Australia (upper). The lower panel shows the boundaries of the local government areas within the GSMA, with the names of some major local cities (dots).</p>
Full article ">Figure 2
<p>Spatial distribution of the rainfall stations with abbreviated codes as in <a href="#atmosphere-12-00627-t0A1" class="html-table">Table A1</a> (<a href="#app1-atmosphere-12-00627" class="html-app">Appendix A</a>). Topography in the area can refer to the DEM in Figure 4.</p>
Full article ">Figure 3
<p>The empirical concentration curves for Albion Park and Wombeyan rain stations with the dash straight line the reference equidistribution line. The area <span class="html-italic">S</span>′ is that bounded by the diagonal equidistribution line and the concentration curve, while area <span class="html-italic">S</span> is the remaining area underneath the equidistribution line.</p>
Full article ">Figure 4
<p>Spatial distribution of DPCI (contours) overlaid on a DEM inside the GSMA.</p>
Full article ">Figure 5
<p>Spatial distribution of the constant “<span class="html-italic">b</span>” (shaded) with classification into four groups of intensity and occurrence locations of the thundery flash flooding days during 1989–2015 (lightning symbol) based on the BoM storm archive.</p>
Full article ">Figure 6
<p>The result of cross−covariance surfaces for all pairs of variables: (CI,AP); (CI,CV); (CI,TN); (CI,MxR); (CI,a); (CI,b). CI is the same as the DPCI. Six bins have been set for the categories of covariance values. The arrows (with the blue and red lines) are examples of directional searches of high covariance values over the surfaces. See <a href="#app3-atmosphere-12-00627" class="html-app">Appendix C</a> for details.</p>
Full article ">Figure A1
<p>Illustration of the Global Moran’s statistic and association with the dispersed, random and clustered precipitation patterns.</p>
Full article ">Figure A2
<p>Cross-covariance surface or cloud between the DPCI and the constant “<span class="html-italic">b</span>” in the concentration curve. The upper panel has all the covariance values according to the distance between the two points for computing the covariance. The lower panel has the cross-covariance map and an example of the directional search arrow. The search arrow is set based on the parameters on the right hand side (such as the direction, lag size and number of lags).</p>
Full article ">
13 pages, 4031 KiB  
Article
Impact of Black Carbon on Surface Ozone in the Yangtze River Delta from 2015 to 2018
by Yue Tan, Delong Zhao, Honglei Wang, Bin Zhu, Dongping Bai, Ankang Liu, Shuangshuang Shi and Qihang Dai
Atmosphere 2021, 12(5), 626; https://doi.org/10.3390/atmos12050626 - 13 May 2021
Cited by 5 | Viewed by 2696
Abstract
Despite the yearly decline in PM2.5 in China, surface ozone has been rapidly increasing recently, which makes it imperative to coordinate and control both PM2.5 and ozone in the atmosphere. This study utilized the data of pollutant concentrations and meteorological elements [...] Read more.
Despite the yearly decline in PM2.5 in China, surface ozone has been rapidly increasing recently, which makes it imperative to coordinate and control both PM2.5 and ozone in the atmosphere. This study utilized the data of pollutant concentrations and meteorological elements during 2015 to 2018 in Nanjing, China to analyze the daily correlation between black carbon and ozone (CBO), and the distribution of the pollutant concentrations under different levels of CBO. Besides, the diurnal variations of pollutant concentrations and meteorological elements under high positive and negative CBO were discussed and compared. The results show that the percentage of positive CBO had been increasing at the average rate of 7.1%/year, and it was 38.7% in summer on average, nearly twice of that in other seasons (19.2%). The average black carbon (BC), PM2.5 and NO2 under positive CBO was lower than those under negative CBO. It is noticeable that the surface ozone began to ascend when CBO was up to 0.2, with PM2.5 and NO2 decreasing and BC remaining steady. Under negative CBO, pollutant concentrations and meteorological elements showed obvious diurnal variations: BC showed a double-peak pattern and surface ozone, PM2.5, SO2 and CO showed single-peak patterns and NO2 showed a trough from 10:00 to 19:00. Wind speed and visibility showed a single-peak pattern with little seasonal difference. Relative humidity rose first, then it lowered and then it rose. Under positive CBO, the patterns of diurnal variations became less obvious, and some of them even showed no patterns, but just fluctuated at a certain level. Full article
Show Figures

Figure 1

Figure 1
<p>The scatter plot and the linear fit of O<sub>3</sub> and BC concentrations.</p>
Full article ">Figure 1 Cont.
<p>The scatter plot and the linear fit of O<sub>3</sub> and BC concentrations.</p>
Full article ">Figure 2
<p>The annual variations of the frequency of negative CBO and positive CBO.</p>
Full article ">Figure 3
<p>The monthly variations of the frequency of negative CBO and positive CBO.</p>
Full article ">Figure 4
<p>The seasonal variations of the frequency of negative CBO and positive CBO.</p>
Full article ">Figure 5
<p>Mass concentrations of pollutants (BC, PM<sub>2.5</sub>, O<sub>3</sub>, NO<sub>2</sub>) under different levels of CBO.</p>
Full article ">Figure 6
<p>Diurnal variations of BC in two situations: (<b>left</b>): −1.0 &lt; CBO &lt; −0.5 and (<b>right</b>): 0.5 &lt; CBO &lt; 1.0.</p>
Full article ">Figure 7
<p>Diurnal variations of O<sub>3</sub> in two situations: (<b>left</b>): −1.0 &lt; CBO &lt; −0.5 and (<b>right</b>): 0.5 &lt; CBO &lt; 1.0.</p>
Full article ">Figure 8
<p>Diurnal variations of PM<sub>2.5</sub> in two situations: (<b>left</b>): −1.0 &lt; CBO &lt; −0.5 and (<b>right</b>): 0.5 &lt; CBO &lt; 1.0.</p>
Full article ">Figure 9
<p>Diurnal variations of NO<sub>2</sub> in two situations: (<b>left</b>): −1.0 &lt; CBO &lt; −0.5 and (<b>right</b>): 0.5 &lt; CBO &lt; 1.0.</p>
Full article ">Figure 10
<p>Diurnal variations of SO<sub>2</sub> in two situations: (<b>left</b>): −1.0 &lt; CBO &lt; −0.5 and (<b>right</b>): 0.5 &lt; CBO &lt; 1.0.</p>
Full article ">Figure 11
<p>Diurnal variations of CO in two situations: (<b>left</b>): −1.0 &lt; CBO &lt; −0.5 and (<b>right</b>): 0.5 &lt; CBO &lt; 1.0.</p>
Full article ">Figure 12
<p>Diurnal variations of wind speed in two situations: (<b>left</b>): −1.0 &lt; CBO &lt; −0.5 and (<b>right</b>): 0.5 &lt; CBO &lt; 1.0.</p>
Full article ">Figure 13
<p>Diurnal variations of relative humidity in two situations: (<b>left</b>): −1.0 &lt; CBO &lt; −0.5 and (<b>right</b>): 0.5 &lt; CBO &lt; 1.0.</p>
Full article ">Figure 14
<p>Diurnal variations of visibility in two situations: (<b>left</b>): −1.0 &lt; CBO &lt; −0.5 and (<b>right</b>): 0.5 &lt; CBO &lt; 1.0.</p>
Full article ">
22 pages, 1751 KiB  
Article
The Diurnal Variation in Stratospheric Ozone from MACC Reanalysis, ERA-Interim, WACCM, and Earth Observation Data: Characteristics and Intercomparison
by Ansgar Schanz, Klemens Hocke, Niklaus Kämpfer, Simon Chabrillat, Antje Inness, Mathias Palm, Justus Notholt, Ian Boyd, Alan Parrish and Yasuko Kasai
Atmosphere 2021, 12(5), 625; https://doi.org/10.3390/atmos12050625 - 13 May 2021
Cited by 5 | Viewed by 3637
Abstract
In this study, we compare the diurnal variation in stratospheric ozone of the MACC (Monitoring Atmospheric Composition and Climate) reanalysis, ECMWF Reanalysis Interim (ERA-Interim), and the free-running WACCM (Whole Atmosphere Community Climate Model). The diurnal variation of stratospheric ozone results from photochemical and [...] Read more.
In this study, we compare the diurnal variation in stratospheric ozone of the MACC (Monitoring Atmospheric Composition and Climate) reanalysis, ECMWF Reanalysis Interim (ERA-Interim), and the free-running WACCM (Whole Atmosphere Community Climate Model). The diurnal variation of stratospheric ozone results from photochemical and dynamical processes depending on altitude, latitude, and season. MACC reanalysis and WACCM use similar chemistry modules and calculate a similar diurnal cycle in ozone when it is caused by a photochemical variation. The results of the two model systems are confirmed by observations of the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) experiment and three selected sites of the Network for Detection of Atmospheric Composition Change (NDACC) at Mauna Loa, Hawaii (tropics), Bern, Switzerland (midlatitudes), and Ny-Ålesund, Svalbard (high latitudes). On the other hand, the ozone product of ERA-Interim shows considerably less diurnal variation due to photochemical variations. The global maxima of diurnal variation occur at high latitudes in summer, e.g., near the Arctic NDACC site at Ny-Ålesund, Svalbard. The local OZORAM radiometer observes this effect in good agreement with MACC reanalysis and WACCM. The sensed diurnal variation at Ny-Ålesund is up to 8% (0.4 ppmv) due to photochemical variations in summer and negligible during the dynamically dominated winter. However, when dynamics play a major role for the diurnal ozone variation as in the lower stratosphere (100–20 hPa), the reanalysis models ERA-Interim and MACC which assimilate data from radiosondes and satellites outperform the free-running WACCM. Such a domain is the Antarctic polar winter where a surprising novel feature of diurnal variation is indicated by MACC reanalysis and ERA-Interim at the edge of the polar vortex. This effect accounts for up to 8% (0.4 ppmv) in both model systems. In summary, MACC reanalysis provides a global description of the diurnal variation of stratospheric ozone caused by dynamics and photochemical variations. This is of high interest for ozone trend analysis and other research which is based on merged satellite data or measurements at different local time. Full article
(This article belongs to the Section Air Quality)
Show Figures

Figure 1

Figure 1
<p>Relative diurnal variation in ozone from WACCM for March 2000, MACC reanalysis and ERA-Interim for March 2012 at 5 hPa over (<b>a</b>) Bern, Switzerland (46.9<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> N, 7.4<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> E) and (<b>b</b>) Mauna Loa, Hawaii (19.5<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> N, 204.5<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> E). The figures show the relative diurnal variation over local time according to Equation (1). The SMILES curves (orange) are for a related period (1 March–21 April 2010) and averaged over belts from 20–50<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> N for Bern and 20<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> S–20<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> N for Hawaii.</p>
Full article ">Figure 2
<p>Zonal-mean <math display="inline"> <semantics> <msub> <mi>D</mi> <msub> <mi mathvariant="normal">O</mi> <mn>3</mn> </msub> </msub> </semantics> </math>/O3 as function of latitude and pressure derived from MACC reanalysis. The figure shows monthly means in the middle and upper stratosphere for (<b>a</b>) March, (<b>b</b>) June, (<b>c</b>) September, and (<b>d</b>) December of 2012 (according to Equation (<a href="#FD7-atmosphere-12-00625" class="html-disp-formula">7</a>)). The dashed, magenta lines refer to the polar circles.</p>
Full article ">Figure 3
<p>Same as <a href="#atmosphere-12-00625-f002" class="html-fig">Figure 2</a> but derived from WACCM for the year 2000.</p>
Full article ">Figure 4
<p>Same as <a href="#atmosphere-12-00625-f002" class="html-fig">Figure 2</a> but derived from ERA-Interim.</p>
Full article ">Figure 5
<p>Seasonal behavior of <math display="inline"> <semantics> <msub> <mi>D</mi> <msub> <mi mathvariant="normal">O</mi> <mn>3</mn> </msub> </msub> </semantics> </math>/O3 (see Equation (<a href="#FD7-atmosphere-12-00625" class="html-disp-formula">7</a>)) at 3 hPa derived from (<b>a</b>) WACCM, (<b>b</b>) MACC reanalysis, and (<b>c</b>) ERA-Interim (the latter two from 2012). The presented continuous monthly means are achieved by a sliding time window of <math display="inline"> <semantics> <mrow> <mo>±</mo> <mn>15</mn> </mrow> </semantics> </math> days. The solid contour lines refer to the solar zenith angle at noon. Dashed contour lines show the sunshine duration given in hours.</p>
Full article ">Figure 6
<p>Relative diurnal variation in ozone as function of local time (LT) for MACC reanalysis (red markers), ERA-Interim (black markers), WACCM (blue line), and OZORAM (grey line) at 5 hPa (38 km) over Ny-Ålesund, Svalbard (78.9<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> N, 11.9<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> E). The figure shows the relative diurnal variation according to Equation (1) for a summer (<b>a</b>) and winter (<b>b</b>) month. The SMILES climatology does not cover the high latitude of Ny-Ålesund. The summer period is taken from 2011 due to technical problems of the instrument in summer 2012.</p>
Full article ">Figure 7
<p>Ozone time series from MACC reanalysis and ERA-Interim (<b>a</b>) and OZORAM (<b>b</b>) during June 2011 at 3 hPa (40 km) over Ny-Ålesund, Svalbard (78.9<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> N, 11.9<math display="inline"> <semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics> </math> E). The uncertainty range stands for the combined random and systematic standard deviation.</p>
Full article ">Figure 8
<p>Similar to <a href="#atmosphere-12-00625-f007" class="html-fig">Figure 7</a> but for December 2012. The label 01/01 corresponds to 2013/01/01.</p>
Full article ">Figure 9
<p><math display="inline"> <semantics> <msub> <mi>D</mi> <msub> <mi mathvariant="normal">O</mi> <mn>3</mn> </msub> </msub> </semantics> </math>/O3 of WACCM (<b>a,d</b>), MACC reanalysis (<b>b,e</b>), and ERA-Interim (<b>c,f</b>) for June 2012 at 5 hPa. The upper panel is for Arctic summer in the Northern Hemisphere, while the lower panel is for Antarctic winter in the Southern Hemisphere. In the upper panel, the magenta markers refer to the geographical position of OZORAM at Ny-Ålesund, Svalbard. In the lower panel, the marker indicate Troll, Antarctica, where the BAS-MRT has been operating from February 2008 and throughout January 2010. The geographical location of the BAS-MRT is not where the strong effects are indicated by MACC reanalysis and ERA-Interim. The magenta dashed line is the Arctic and Antarctic polar circle.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop