Nothing Special   »   [go: up one dir, main page]

Next Issue
Volume 12, January
Previous Issue
Volume 11, November
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 

Atmosphere, Volume 11, Issue 12 (December 2020) – 116 articles

Cover Story (view full-size image): The Antarctic Peninsula is one of the Earth’s regions with the strongest warming since the mid-20th century. However, its northwest region cooled from 2000 to 2015. Increased snow accumulation was recorded in Livingston Island between 2009 and 2014, with sites that were snow-free in the summer becoming snow covered. The ground thermal regimes from two boreholes show the effects of the increasing snow insulation, resulting in the disappearance of the active layer and on permafrost aggradation. This regime shift may be used as an analogue for the transition from subaerial periglacial to subglacial ground thermal conditions. 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:
24 pages, 11539 KiB  
Article
Characteristics of Clouds and Raindrop Size Distribution in Xinjiang, Using Cloud Radar Datasets and a Disdrometer
by Yong Zeng, Lianmei Yang, Zuyi Zhang, Zepeng Tong, Jiangang Li, Fan Liu, Jinru Zhang and Yufei Jiang
Atmosphere 2020, 11(12), 1382; https://doi.org/10.3390/atmos11121382 - 21 Dec 2020
Cited by 22 | Viewed by 4019
Abstract
Observation data from March to May 2020 of the Ka-band millimeter-wave cloud radar and disdrometer, located in Xinjiang, a typical arid region of China, were used to study the diurnal variation of clouds and precipitation, raindrop size distribution (DSD), and the physical parameters [...] Read more.
Observation data from March to May 2020 of the Ka-band millimeter-wave cloud radar and disdrometer, located in Xinjiang, a typical arid region of China, were used to study the diurnal variation of clouds and precipitation, raindrop size distribution (DSD), and the physical parameters of raindrops. The results showed that there are conspicuous diurnal changes in clouds and precipitation. There is a decreasing trend of the cloud base height (CBH) from 05:00 to 19:00 CST (China Standard Time, UTC +8) and a rising trend of CBHs from 20:00 to 04:00 CST. The cloud top height (CTH) and the cloud thickness show a rising trend from 03:00 to 05:00 CST, 12:00 to 14:00 CST, and 20:00 to 01:00 CST. The diurnal variation of clouds is mainly driven by wind and temperature closely related to the topography of the study area. There are three apparent precipitation periods during the day, namely, 02:00–09:00 CST, 12:00 CST, and 17:00–21:00 CST. The changes in the physical parameters of raindrops are more drastic and evident with a lower CBH, lower CTH, and higher number of cloud layers from 12:00 to 21:00 CST than other times, which are closely related to day-to-day variations of systems moving through, and incoming solar radiation and the mountain–valley wind circulation caused by the trumpet-shaped topography that opens to the west played a secondary role. The DSD is in agreement with a normalized gamma distribution, and the value of the shape factor μ is significantly different from the fixed μ value in the Weather Research and Forecasting (WRF) Model. The rain in arid Xinjiang had a higher concentration of raindrops and a smaller average raindrop diameter than the rain in other humid regions of the Central and Southeast Asian continent. In the ZR (radar reflectivity–rain rate) relationship, Z=249R1.20 is derived for stratiform rain, and it is significantly different from humid regions. Using Z/Dm (mass–weighted mean diameter) and R, a new empirical relationship Z/Dm=214R1.20 is established, and improvement is obtained in rain retrieval by using the Z/DmR relation relative to the conventional ZR relation. Additionally, the NtR, DmR, NwR, and NtNw relationships with larger differences from humid regions are established by fitting the power-law equations. These results are useful for improving the data parameters of microphysical processes of WRF and the accuracy of quantitative precipitation estimation in arid regions. Full article
(This article belongs to the Section Meteorology)
Show Figures

Figure 1

Figure 1
<p>Location of the Xinyuan Meteorological Station (XY, 83.25° E, 43.45° N, 929.7 m above sea level) and the photos of the millimeter-wave cloud radar (MMCR) and disdrometer at the Xinyuan Meteorological Observatory.</p>
Full article ">Figure 2
<p>Flow chart of classifying stratiform, convective, stratocumulus rain, and non-rainfall.</p>
Full article ">Figure 3
<p>Diurnal variations in (<b>a</b>) the cloud base height (CBH) (m), (<b>b</b>) the cloud top height (CTH) (m), (<b>c</b>) the cloud thickness (CTK) (m), (<b>d</b>) the frequency distribution of the CTK, and (<b>e</b>) the frequency distribution of the cloud layer number (CLN) during the observation period. On each box, the blue square represents the average, the central crossbar indicates the median, and the bottom and top edges of the box indicate the 25th and 75th percentiles, respectively. The black whiskers extend to the maximum and minimum of the data points, respectively. The positions of the purple and green vertical lines are midday and midnight in <a href="#atmosphere-11-01382-f003" class="html-fig">Figure 3</a>a, respectively.</p>
Full article ">Figure 3 Cont.
<p>Diurnal variations in (<b>a</b>) the cloud base height (CBH) (m), (<b>b</b>) the cloud top height (CTH) (m), (<b>c</b>) the cloud thickness (CTK) (m), (<b>d</b>) the frequency distribution of the CTK, and (<b>e</b>) the frequency distribution of the cloud layer number (CLN) during the observation period. On each box, the blue square represents the average, the central crossbar indicates the median, and the bottom and top edges of the box indicate the 25th and 75th percentiles, respectively. The black whiskers extend to the maximum and minimum of the data points, respectively. The positions of the purple and green vertical lines are midday and midnight in <a href="#atmosphere-11-01382-f003" class="html-fig">Figure 3</a>a, respectively.</p>
Full article ">Figure 4
<p>Diurnal variations in (<b>a</b>) accumulated rain amount, (<b>b</b>) accumulated rainy minute number, (<b>c</b>) precipitation intensity (<span class="html-italic">R</span>) (mm · h<sup>−1</sup>), (<b>d</b>) number concentration of total raindrops (<span class="html-italic">N<sub>t</sub></span>) (m<sup>−3</sup>), (<b>e</b>) mass-weighted mean diameter (<span class="html-italic">D<sub>m</sub></span>) (mm), (<b>f</b>) normalized intercept parameter (<span class="html-italic">N<sub>w</sub></span>) (mm<sup>−1</sup> · m<sup>−3</sup>), (<b>g</b>) radar reflectivity (<span class="html-italic">Z</span>) (dBZ), and (<b>h</b>) liquid water content (<span class="html-italic">LWC</span>) (g · m<sup>−</sup><sup>3</sup>) during the observation period. The green line represents the average, on each box, the central crossbar indicates the median, and the bottom and top edges of the box indicate the 25th and 75th percentiles, respectively. The black whiskers extend to the maximum and minimum of the data points, respectively.</p>
Full article ">Figure 5
<p>Raindrop size distributions (DSDs) and Gamma distribution fitting curves for stratiform and stratocumulus precipitation.</p>
Full article ">Figure 6
<p>(<b>a</b>) Scatterplots of the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>log</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <msub> <mi>N</mi> <mi>t</mi> </msub> <mo>−</mo> <mi>R</mi> </mrow> </semantics></math> relationship, (<b>b</b>) scatterplots of the <math display="inline"><semantics> <mrow> <msub> <mi>D</mi> <mi>m</mi> </msub> <mo>−</mo> <mi>R</mi> </mrow> </semantics></math> relationship, (<b>c</b>) scatterplots of the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>log</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <msub> <mi>N</mi> <mi>w</mi> </msub> <mo>−</mo> <mi>R</mi> </mrow> </semantics></math> relationship, and (<b>d</b>) scatterplots of the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>log</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <msub> <mi>N</mi> <mi>w</mi> </msub> <mo>−</mo> <msub> <mrow> <mi>log</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <msub> <mi>N</mi> <mi>t</mi> </msub> </mrow> </semantics></math> relationship. The red dashed curve is the fitted curve for stratiform precipitation in Xinjiang.</p>
Full article ">Figure 7
<p>Scatterplots of <span class="html-italic">Z</span> versus <span class="html-italic">R</span> and the fitted relationship for stratiform rainfall. The red curve is the fitted curve for stratiform precipitation in Xinjiang, and the orange circles represent the scattered points of the <span class="html-italic">Z–R</span> relationship obtained from the disdrometer data.</p>
Full article ">Figure 8
<p>Scatterplots of <math display="inline"><semantics> <mrow> <mi>Z</mi> <mo>/</mo> <msub> <mi>D</mi> <mi>m</mi> </msub> </mrow> </semantics></math> versus <span class="html-italic">R</span> and the fitted relationship for stratiform rainfall. The red dashed curve is the fitted curve for stratiform precipitation in Xinjiang.</p>
Full article ">Figure 9
<p>(<b>a</b>) Temporal variation of estimated rain intensity (mm · h<sup>−</sup><sup>1</sup>) by using <math display="inline"><semantics> <mrow> <mi>Z</mi> <mo>−</mo> <mi>R</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>Z</mi> <mo>/</mo> <msub> <mi>D</mi> <mi>m</mi> </msub> <mo>−</mo> <mi>R</mi> </mrow> </semantics></math> relations and measured by the disdrometer and (<b>b</b>) <span class="html-italic">Z</span> (dBZ), <span class="html-italic">Dm</span> (mm), <span class="html-italic">LWC</span> (g · m<sup>−</sup><sup>3</sup>), and <span class="html-italic">N<sub>t</sub></span> (m<sup>−</sup><sup>3</sup>) on 20–21 May 2020.</p>
Full article ">Figure 10
<p>Time–height profile of the (<b>a</b>) radar reflectivity factor (dBZ), (<b>b</b>) Doppler velocity (m · s<sup>−</sup><sup>1</sup>), and (<b>c</b>) velocity spectrum width (m · s<sup>−</sup><sup>1</sup>).</p>
Full article ">Figure 10 Cont.
<p>Time–height profile of the (<b>a</b>) radar reflectivity factor (dBZ), (<b>b</b>) Doppler velocity (m · s<sup>−</sup><sup>1</sup>), and (<b>c</b>) velocity spectrum width (m · s<sup>−</sup><sup>1</sup>).</p>
Full article ">Figure 11
<p>Diurnal variations in (<b>a</b>) temperature at 2 m (°C), and (<b>b</b>) wind direction at 10 m (°) during the observation period at Xinyuan Meteorological Station. On each box, the central crossbar indicates the median, and the bottom and top edges of the box indicate the 25th and 75th percentiles, respectively. The black whiskers extend to the maximum and minimum of the data points, respectively. The green dashed, blue dashed and purple dashed lines are the positions of east wind, south wind and west wind, respectively.</p>
Full article ">Figure 11 Cont.
<p>Diurnal variations in (<b>a</b>) temperature at 2 m (°C), and (<b>b</b>) wind direction at 10 m (°) during the observation period at Xinyuan Meteorological Station. On each box, the central crossbar indicates the median, and the bottom and top edges of the box indicate the 25th and 75th percentiles, respectively. The black whiskers extend to the maximum and minimum of the data points, respectively. The green dashed, blue dashed and purple dashed lines are the positions of east wind, south wind and west wind, respectively.</p>
Full article ">Figure 12
<p>Schematic diagram of the diurnal variation of clouds (the green five-pointed star represents the position of the new source).</p>
Full article ">Figure 13
<p>The start, duration and end time of each rain event.</p>
Full article ">Figure 14
<p>The TBB (black-body temperature) from FY-2G (Fengyun 2-G satellite) at the beginning of each rainfall event, rainfall event 1 to rainfall event 17 are (<b>a</b>–<b>p</b>), respectively (except rainfall event 16, because rainfall event 16 and rainfall event 15 are affected by the same system).</p>
Full article ">
16 pages, 1701 KiB  
Article
Increases in Biogenic Volatile Organic Compound Concentrations Observed after Rains at Six Forest Sites in Non-Summer Periods
by Takafumi Miyama, Tomoaki Morishita, Yuji Kominami, Hironori Noguchi, Yukio Yasuda, Natsuko Yoshifuji, Michiaki Okano, Katsumi Yamanoi, Yasuko Mizoguchi, Satoru Takanashi, Kenzo Kitamura and Kazuho Matsumoto
Atmosphere 2020, 11(12), 1381; https://doi.org/10.3390/atmos11121381 - 21 Dec 2020
Cited by 8 | Viewed by 3221
Abstract
Since biogenic volatile organic compounds (BVOCs) are important precursors of ozone, the monitoring of the BVOC concentration distributions is needed. In general, forest BVOC concentrations increase in summer as well as in other seasons. This study aims to detect temporally sporadic increases in [...] Read more.
Since biogenic volatile organic compounds (BVOCs) are important precursors of ozone, the monitoring of the BVOC concentration distributions is needed. In general, forest BVOC concentrations increase in summer as well as in other seasons. This study aims to detect temporally sporadic increases in BVOC concentrations in the non-summer months and to analyze the occurring climatic conditions. Using a uniform sampling system and shared gas chromatography–mass spectrometry, the concentrations of isoprene and monoterpenes in six Japanese forests were observed approximately once a month for 3 years. Using the observed data, the relations between the BVOC concentration increases and meteorological factors were evaluated. Twenty instances of temporal increases in BVOC concentrations were observed. These mainly occurred in spring for isoprene and in autumn for monoterpenes. Most of the increases in the non-summer months were observed after a rainfall event, when the daily temperature range was large, suggesting that wind, rain, and a rapid diurnal temperature rise could be factors in the non-summer months. Thus, the network monitoring of BVOC concentrations might be effective for understanding the effects of factors other than temperature, and the mechanisms and frequency of the temporal increases, on the BVOC concentrations in various forests. Full article
(This article belongs to the Special Issue Tropospheric Ozone Observations)
Show Figures

Figure 1

Figure 1
<p>Locations of the six measurement sites in Japan.</p>
Full article ">Figure 2
<p>Relation between the air temperature (Ta) and (<b>a</b>) isoprene (ISO, above) and (<b>b</b>) monoterpene (MT, below) concentrations.</p>
Full article ">Figure 3
<p>Seasonal variation in the (<b>a</b>) isoprene (ISO, above) or (<b>b</b>) monoterpene (MT, below) concentration.</p>
Full article ">Figure 3 Cont.
<p>Seasonal variation in the (<b>a</b>) isoprene (ISO, above) or (<b>b</b>) monoterpene (MT, below) concentration.</p>
Full article ">Figure 4
<p>Box plot of the (<b>a</b>) isoprene (ISO, left) and (<b>b</b>) monoterpene (MT, right) concentrations. The blank small squares, symbols, whiskers, and boxes indicate the mean, outlier, minimum, maximum, and 25, 50, and 75 percentiles, respectively. The means, identified by letters, significantly differ (<span class="html-italic">p</span> &lt; 0.05) based on the Bonferroni post-hoc tests.</p>
Full article ">
16 pages, 231 KiB  
Commentary
Lessons Learned from Coupled Fire-Atmosphere Research and Implications for Operational Fire Prediction and Meteorological Products Provided by the Bureau of Meteorology to Australian Fire Agencies
by Mika Peace, Joseph Charney and John Bally
Atmosphere 2020, 11(12), 1380; https://doi.org/10.3390/atmos11121380 - 21 Dec 2020
Cited by 6 | Viewed by 2736
Abstract
Coupled fire-atmosphere models are simulators that integrate a fire component and an atmospheric component, with the objective of capturing interactions between the fire and atmosphere. As a fire releases energy in the combustion process, the surrounding atmosphere adjusts in response to the energy [...] Read more.
Coupled fire-atmosphere models are simulators that integrate a fire component and an atmospheric component, with the objective of capturing interactions between the fire and atmosphere. As a fire releases energy in the combustion process, the surrounding atmosphere adjusts in response to the energy fluxes; coupled fire-atmosphere (CFA) models aim to resolve the processes through which these adjustments occur. Several CFA models have been developed internationally, mostly by meteorological institutions and primarily for use as a research tool. Research studies have provided valuable insights into some of the atmospheric processes surrounding a fire. The potential to run CFA models in real time is currently limited due to the intensive computational requirements. In addition, there is a need for systematic verification to establish their accuracy and the appropriate circumstances for their use. The Bureau of Meteorology (the Bureau) is responsible for providing relevant and accurate meteorological information to Australian fire agencies to inform decisions for the protection of life and property and to support hazard management activities. The inclusion of temporally and spatially detailed meteorological fields that adjust in response to the energy released by a fire is seen as a component in developing fire prediction systems that capture some of the most impactful fire and weather behavior. The Bureau’s ten-year research and development plan includes a commitment to developing CFA models, with the objective of providing enhanced services to Australian fire agencies. This paper discusses the operational use of fire predictions and simulators, learnings from CFA models and potential future directions for the Bureau in using CFA models to support fire prediction activities. Full article
(This article belongs to the Special Issue Coupled Fire-Atmosphere Simulation)
10 pages, 3509 KiB  
Article
Atmospheric Mercury Deposition in Macedonia from 2002 to 2015 Determined Using the Moss Biomonitoring Technique
by Trajče Stafilov, Lambe Barandovski, Robert Šajn and Katerina Bačeva Andonovska
Atmosphere 2020, 11(12), 1379; https://doi.org/10.3390/atmos11121379 - 21 Dec 2020
Cited by 13 | Viewed by 2708
Abstract
The moss biomonitoring technique was used in 2002, 2005, 2010 and 2015 in a potentially toxic elements study (PTEs) in Macedonia. For that purpose, more than 70 moss samples from two dominant species (Hypnum cupressiforme and Homalothecium lutescens) were collected during the [...] Read more.
The moss biomonitoring technique was used in 2002, 2005, 2010 and 2015 in a potentially toxic elements study (PTEs) in Macedonia. For that purpose, more than 70 moss samples from two dominant species (Hypnum cupressiforme and Homalothecium lutescens) were collected during the summers of the mentioned years. Total digestion of the samples was done using a microwave digestion system, whilst mercury was analyzed by cold vapour atomic absorption spectrometry (CV–AAS). The content of mercury ranged from 0.018 mg/kg to 0.26 mg/kg in 2002, from 0.010 mg/kg to 0.42 mg/kg in 2005, from 0.010 mg/kg to 0.60 mg/kg in 2010 and from 0.020 mg/kg to 0.25 mg/kg in 2015. Analysis of the median values shows the increase of the content in the period 2002–2010 and a slight reduction of the air pollution with Hg in the period 2010–2015. Mercury distribution maps show that sites with increased concentrations of mercury in moss are likely impacted by anthropogenic pollution. The results were compared to similar studies done during the same years in neighboring countries and in Norway—which is a pristine area and serves as a reference, and it was concluded that mercury air pollution in Macedonia is significant primarily in industrialized regions. Full article
Show Figures

Figure 1

Figure 1
<p>Map of the Republic of Macedonia.</p>
Full article ">Figure 2
<p>Location of sampling points.</p>
Full article ">Figure 3
<p>Box plots of Hg according to the year of sampling.</p>
Full article ">Figure 4
<p>Distribution of Hg in moss samples collected in 2002, 2005, 2010 and 2015.</p>
Full article ">
20 pages, 8133 KiB  
Article
Circulation Specific Precipitation Patterns over Svalbard and Projected Future Changes
by Andreas Dobler, Julia Lutz, Oskar Landgren and Jan Erik Haugen
Atmosphere 2020, 11(12), 1378; https://doi.org/10.3390/atmos11121378 - 21 Dec 2020
Cited by 6 | Viewed by 3393
Abstract
Precipitation on Svalbard can generally be linked to the atmospheric circulation in the Northern Atlantic. Using an automated circulation type classification, we show that weather type statistics are well represented in the Max Planck Institute Earth System Model at base resolution (MPI-ESM-LR). For [...] Read more.
Precipitation on Svalbard can generally be linked to the atmospheric circulation in the Northern Atlantic. Using an automated circulation type classification, we show that weather type statistics are well represented in the Max Planck Institute Earth System Model at base resolution (MPI-ESM-LR). For a future climate projection following the Representative Concentration Pathway scenario RCP8.5, we find only small changes in the statistics. However, convection permitting simulations with the regional climate model from the Consortium for Small-scale Modeling in climate mode (COSMO-CLM) covering Svalbard at 2.5 km demonstrate an increase in precipitation for all seasons. About 74% of the increase are coming from changes under cyclonic weather situations. The precipitation changes are strongly related to differences in atmospheric conditions, while the contribution from the frequencies of weather types is small. Observations on Svalbard suggest that the general weather situation favouring heavy precipitation events is a strong south-southwesterly flow with advection of water vapour from warmer areas. This is reproduced by the COSMO-CLM simulations. In the future projections, the maximum daily precipitation amounts are further increasing. At the same time, weather types with less moisture advection towards Svalbard are becoming more important. Full article
Show Figures

Figure 1

Figure 1
<p>Orography in the intermediate (<b>left</b>) and high-resolution (<b>right</b>) regional climate model domain. The location of Svalbard Airport is indicated by a circle in the high-resolution domain.</p>
Full article ">Figure 2
<p>Annual and seasonal mean sea level pressure from ERA-Interim (1979–2008, as contour lines) and the current climate MPI-ESM-LR simulation (1971–2000, in colours). The numbered points and dashed box denote the grid points and bounding box used for the weather type classification, respectively. The labels in the top left corners refer to the annual mean (ANN) or the months of the specific season of the year.</p>
Full article ">Figure 3
<p>Mean monthly occurrence of weather types in ERA-Interim (1979–2008, <b>left</b>) and the MPI-ESM-LR current climate simulation (1971–2000, <b>right</b>), respectively. The dots mark statistically significant differences.</p>
Full article ">Figure 4
<p>Annual mean sea level pressure from MPI-ESM-LR (1971–2000, in colours) and ERA-Interim (1979–2008, as contour lines) for the current climate (top left) and the single weather types. The bold grey line shows the 1010 hPa level in the ERA-Interim fields and the dashed box the area used for the weather type classification. The labels in the top left corners refer to the weather type.</p>
Full article ">Figure 5
<p>Annual mean precipitation (top left) from the MPI-ESM-LR driven COSMO-CLM current climate (1971–2000) simulation and its decomposition into single weather types over Svalbard. For the annual mean, the average daily precipitation over the whole land area is shown in the bottom right corner. For the different weather types, the contribution to the total precipitation is given in the bottom right corner, while the number in the top left corresponds to the number of occurrences per year.</p>
Full article ">Figure 6
<p>(<b>Left</b>): Difference between the annual and seasonal mean precipitation from the MPI-ESM-LR and the ERA-Interim driven COSMO-CLM simulation for the current climate (1971–2000 and 2004–2017, respectively). Statistically significant differences are outlined. (<b>Center</b>): The part resulting from differences between the weather type frequencies in MPI-ESM-LR and ERA-Interim. (<b>Right</b>): The part resulting from large-scale condition differences. The labels in the top left corners refer to the annual mean (ANN) or the months of the specific season of the year.</p>
Full article ">Figure 7
<p>Projected changes from 1971–2000 to 2071–2100 in the mean annual and seasonal sea level pressure in the MPI-ESM-LR climate simulation. Statistically significant changes are marked with dots. The labels in the top left corners refer to the annual mean (ANN) or the months of the specific season of the year.</p>
Full article ">Figure 8
<p>Mean monthly changes in the occurrence of weather types in the MPI-ESM-LR RCP8.5 projection (2071–2100) compared to the historical simulation (1971–2000). The dots mark statistically significant changes.</p>
Full article ">Figure 9
<p>Annual mean sea level pressure from the future (2071–2100, in colours) and current (1971–2000, as contour lines) MPI-ESM-LR simulation (top left) and the decomposition for single weather types. The bold grey line shows the 1010 hPa level in the current climate and the dashed box the area used for the weather type classification. In the top left corner, the number of events per year in the future and current (in parentheses) simulation is given.</p>
Full article ">Figure 10
<p>Annual mean precipitation change from 1971–2000 to 2071–2100 (top left) and its decomposition into single weather types. Statistically significant changes are outlined. For the annual mean, the number in the bottom right corner gives the mean increase. For the weather types, the numbers give the contributions to the change and the contribution to the future mean precipitation (in parentheses), respectively. In the top left corner, the frequencies (number of days per year) for the future and current climate (in parentheses) are given.</p>
Full article ">Figure 11
<p>Mean precipitation changes (from 1971–2000 to 2071–2100) for the weather type Nc in winter (<b>left</b>) and SWc in summer (<b>right</b>). Statistically significant changes are outlined.</p>
Full article ">Figure 12
<p>Changes in precipitation (from 1971–2000 to 2071–2100) decomposed into the contributions from frequency and large-scale condition changes. (<b>Left</b>): Changes in the annual (top) and seasonal mean precipitation climate projections over Svalbard. Statistically significant changes are outlined. (<b>Center</b>): Precipitation changes resulting from frequency changes only. (<b>Right</b>): Differences resulting from large-scale condition changes only.</p>
Full article ">Figure 13
<p>Distributions for current (1971–2000, orange) and future (2071–2100, brown) mean daily precipitation values over Svalbard associated with the different weather types. The crosses denote the contribution to the total precipitation (in %) and the circles the mean precipitation (in mm/day) for each weather type. The horizontal line in the box indicates the median, while bottom and top edges indicate the first and third quartile, respectively. Whiskers extend to the most extreme value which is no more than 1.5 times the interquartile range from the box. Data points outside the whisker range are plotted as outliers. The width of the box indicates the number of events.</p>
Full article ">Figure 14
<p>Precipitable water (colours) and sea level pressure (contours) for the top three precipitation events in the 1971–2000 (<b>left</b>) and 2071–2100 (<b>right</b>) climate simulations. The labels in the top left corner denote the corresponding weather type.</p>
Full article ">
15 pages, 2614 KiB  
Article
Quantifying the Independent Influences of Land Cover and Humidity on Microscale Urban Air Temperature Variation in Hot Summer: Methods of Path Analysis and Genetic SVR
by Weifang Shi, Nan Wang, Aixuan Xin, Linglan Liu, Jiaqi Hou and Yirui Zhang
Atmosphere 2020, 11(12), 1377; https://doi.org/10.3390/atmos11121377 - 21 Dec 2020
Cited by 2 | Viewed by 2165
Abstract
Mitigating high air temperatures and heat waves is vital for decreasing air pollution and protecting public health. To improve understanding of microscale urban air temperature variation, this paper performed measurements of air temperature and relative humidity in a field of Wuhan City in [...] Read more.
Mitigating high air temperatures and heat waves is vital for decreasing air pollution and protecting public health. To improve understanding of microscale urban air temperature variation, this paper performed measurements of air temperature and relative humidity in a field of Wuhan City in the afternoon of hot summer days, and used path analysis and genetic support vector regression (SVR) to quantify the independent influences of land cover and humidity on air temperature variation. The path analysis shows that most effect of the land cover is mediated through relative humidity difference, more than four times as much as the direct effect, and that the direct effect of relative humidity difference is nearly six times that of land cover, even larger than the total effect of the land cover. The SVR simulation illustrates that land cover and relative humidity independently contribute 16.3% and 83.7%, on average, to the rise of the air temperature over the land without vegetation in the study site. An alternative strategy of increasing the humidity artificially is proposed to reduce high air temperatures in urban areas. The study would provide scientific support for the regulation of the microclimate and the mitigation of the high air temperature in urban areas. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic of the field measurement.</p>
Full article ">Figure 2
<p>Two priori models for path analysis. (<b>a</b>) Priori model one; (<b>b</b>) Priori model two.</p>
Full article ">Figure 3
<p>Framework of the genetic support vector regression.</p>
Full article ">Figure 4
<p>Over flowchart of genetic support vector regression (SVR) training to optimize parameters.</p>
Full article ">Figure 5
<p>Box plots of Air Temperature Difference (∆<span class="html-italic">Ta</span>) (<b>a</b>) and Relative Humidity Difference (∆<span class="html-italic">RH</span>) (<b>b</b>).</p>
Full article ">Figure 6
<p>Air Temperature Difference (∆<span class="html-italic">Ta</span>) varying with Distance (∆<span class="html-italic">X</span>) at ∆<span class="html-italic">RH</span> = 0.</p>
Full article ">Figure 7
<p>Air Temperature Difference (∆<span class="html-italic">Ta</span>) varying with Relative Humidity Difference (∆<span class="html-italic">RH</span>) at ∆<span class="html-italic">X</span> = 0.</p>
Full article ">Figure 8
<p>Air Temperature Difference (∆<span class="html-italic">Ta</span>) varying with Relative Humidity Difference (∆<span class="html-italic">RH</span>) at ∆<span class="html-italic">X</span> = 30 m.</p>
Full article ">
13 pages, 6844 KiB  
Article
A 5000-Year Sedimentary Record of East Asian Winter Monsoon from the Northern Muddy Area of the East China Sea
by Yanping Chen, Yan Li, Wenzhe Lyu, Dong Xu, Xibin Han, Tengfei Fu and Liang Yi
Atmosphere 2020, 11(12), 1376; https://doi.org/10.3390/atmos11121376 - 20 Dec 2020
Cited by 8 | Viewed by 3111
Abstract
The variability of the winter monsoon is one of the key components of the Asian monsoon, significantly influencing paleoenvironmental evolution in East Asia. However, whether the winter or the summer monsoon is the dominated factor controlling sedimentary dynamics of the muddy areas of [...] Read more.
The variability of the winter monsoon is one of the key components of the Asian monsoon, significantly influencing paleoenvironmental evolution in East Asia. However, whether the winter or the summer monsoon is the dominated factor controlling sedimentary dynamics of the muddy areas of the continental shelves of the East China Sea is debated, due to lack of consistency between various winter monsoon proxies in previous studies. In this work, the sediments of the upper part of core ECS-DZ1 with several marine surface samples were studied in terms of sediment grain size and radiocarbon dating, and changes in sedimentary dynamics of the northern muddy area of the ECS over the past 5000 years were documented. The main findings are as follows: (1) regional sedimentary dynamics were low and did not significantly change since the middle Holocene; (2) coarse particles are the dominated component in the sediments; (3) a proxy can be derived to indicate changes in winter monsoon. Based on this reconstructed winter monsoon record, we found that this record was generally negatively correlated to the stalagmite-based summer monsoon variability over the past 3500 years, but positively correlated before that. Moreover, this record can be well correlated to changes in the Kuroshio Current and the Bond ice-rafting debris events in the North Atlantic on millennial timescales, inferring large-scale and common atmospheric dynamics across the Asian continent over the past 5000 years. Therefore, we concluded that the winter monsoon is the predominant factor controlling sedimentary dynamics in the northern part of the ECS and proposed that the contribution of coarse particles may be one of potential indices to identify the role of the winter and the summer monsoons in sedimentary evolution. Full article
(This article belongs to the Section Climatology)
Show Figures

Figure 1

Figure 1
<p>Location of core ECS-DZ1 with 6 sites of surface samples (H1-H6) and regional climatic and oceanographic systems. The solid arrows denote warm currents (YSWC, Yellow Sea warm current; KC, Kuroshio current), and the dashed arrows indicate coastal cold waters (YSCC, Yellow Sea coast cold water). The base map data was generated using the open and free software DIVA-GIS 7.5 (<a href="http://www.diva-gis.org/" target="_blank">http://www.diva-gis.org/</a>).</p>
Full article ">Figure 2
<p>Core photos and profile of core ECS-DZ1 about component percentages of clay (&lt;4 μm), silt (4–63 μm), and sand (&gt;63 μm), median size (Md), mean size (Mz), C value (the one percentile of grain size distribution), and AMS 14C dating results with an age-depth model.</p>
Full article ">Figure 3
<p>Characteristics of sediment grain size of core ECS-DZ1 and six surface samples. (<b>a</b>) Ternary diagrams; (<b>b</b>) C-M diagrams; (<b>c</b>) grain size distribution; (<b>d</b>) principal component analysis. Lcomp1 and Lcomp2 are the two characterized grain-size components through mathematical partitioning of sediment grain-size spectrum; and Fcomp1, Fcomp2, and Fcomp3 are the three components by VPCA.</p>
Full article ">Figure 4
<p>Changes in sediment grain-size parameters of core ECS-DZ1.</p>
Full article ">Figure 5
<p>Comparison between various paleoenvironmental proxies over the past 5000 years. (<b>a</b>) 5-point moving average of five grain-size parameters of core ECS-DZ1; (<b>b</b>) Fcom2 record of core ECS-DZ1 indicating changes in the winter monsoon; (<b>c</b>) DET YS-GS index, indicating winter dust storms in Central Asia [<a href="#B40-atmosphere-11-01376" class="html-bibr">40</a>]; (<b>d</b>) ice-rafting debris (IRD), stack ice-rafting debris data from Bond, et al. [<a href="#B41-atmosphere-11-01376" class="html-bibr">41</a>]; (<b>e</b>) Stack KC, integrated indicator from four time series [<a href="#B42-atmosphere-11-01376" class="html-bibr">42</a>]; (<b>f</b>) Stalagmite δ18O series of Sanbao Cave [<a href="#B43-atmosphere-11-01376" class="html-bibr">43</a>], indicating East Asian summer monsoon (EASM) intensity.</p>
Full article ">
26 pages, 25251 KiB  
Article
Graz Lagrangian Model (GRAL) for Pollutants Tracking and Estimating Sources Partial Contributions to Atmospheric Pollution in Highly Urbanized Areas
by Aleksey A. Romanov, Boris A. Gusev, Egor V. Leonenko, Anastasia N. Tamarovskaya, Alexander S. Vasiliev, Nikolai E. Zaytcev and Ilia K. Philippov
Atmosphere 2020, 11(12), 1375; https://doi.org/10.3390/atmos11121375 - 19 Dec 2020
Cited by 19 | Viewed by 9525
Abstract
Computational modeling allows studying the air quality problems in depth and provides the best solution reducing the population risks. This research demonstrates the Graz Lagrangian model effectiveness for assessing emission sources contributions to the air pollution: particles tracking and accumulation estimate. The article [...] Read more.
Computational modeling allows studying the air quality problems in depth and provides the best solution reducing the population risks. This research demonstrates the Graz Lagrangian model effectiveness for assessing emission sources contributions to the air pollution: particles tracking and accumulation estimate. The article describes model setting up parameters and datasets preparation for the analysis. The experiment simulated the dispersion from the main groups of emission sources for real weather conditions during 96 h of December 2018, when significant excess of NO2, CO, SO2, PM10, and benzo(a)pyrene concentrations were observed in the Krasnoyarsk surface atmospheric layer. The computational domain was a parallelepiped of 40 × 30 × 2.5 km, which was located deep inside the Eurasian continent on a heterogeneous landscape exaggerated by high-rise buildings, with various pollutions sources and the ice-free Yenisei River. The results demonstrated an excellent applicability of the Lagrange model for hourly tracking of particle trajectories, taking into account the urban landscape. For values <1 MPC (maximum permissible concentration) of peak pollutants concentrations, the coincidences were 93 cases, and for values < 0.1 shares of MPC, there were 36 cases out of the total number of 97. The same was found for the average daily concentration for values <1 MPC—31, and for values <0.1 MPC—5 matches out of 44. Wind speeds COR—65.3%, wind directions COR—68.6%. The Graz Lagrangian model showed the ability to simulate air quality problems in the Krasnoyarsk greater area conditions. Full article
(This article belongs to the Special Issue Atmospheric Trace Gas Source Detection and Quantification)
Show Figures

Figure 1

Figure 1
<p>Krasnoyarsk city view from a high elevation point: (<b>A</b>)—2 November 2020 (wind direction—SW, wind speed 2 m/s); (<b>B</b>)—28 November 2020 (calm, air temperature −14 Celsius, humidity—85%, atmospheric pressure 1040 mb, AQI from 175–638). Credit: A Romanov.</p>
Full article ">Figure 2
<p>A heterogeneous landscape of the Krasnoyarsk greater area. Spatial distribution of the main emission sources, meteorological stations, and observation points (air quality measuring stations, AQM).</p>
Full article ">Figure 3
<p>The average monthly wind speeds for the period from 1966 to 2019.</p>
Full article ">Figure 4
<p>Calm frequency for the period from 1966 to 2019.</p>
Full article ">Figure 5
<p>Averaging concentration values in a cell with AQM station for further model verification.</p>
Full article ">Figure 6
<p>Wind speeds and directions by the data from meteorological station 29572 (Krasnoyarsk—downtown).</p>
Full article ">Figure 7
<p>Emission sources time profiles designed to account for the intensity from different groups of sources.</p>
Full article ">Figure 8
<p>The Krasnoyarsk city atmosphere vertical profile is characterizing the Yenisei river beds influence on a calm air temperature of minus 25 degrees Celsius and water temperature plus 4 degrees Celsius.</p>
Full article ">Figure 9
<p>Simulated (<b>A</b>) and measured (<b>B</b>) wind flows and speeds distribution for meteorological station code 29570.</p>
Full article ">Figure 10
<p>Simulated and observed wind speeds (meteorological station code 29570).</p>
Full article ">Figure 11
<p>Simulated and observed wind directions (meteorological on code 29570).</p>
Full article ">Figure 12
<p>Spatial distribution of benzo(a)pyrene (BaP) in the surface layer of the Krasnoyarsk atmosphere (2 m above the ground) at the beginning (<b>A</b>) and the end (<b>B</b>) of adverse weather conditions (AWC).</p>
Full article ">Figure 13
<p>Spatial distribution of PM10 in the surface layer of Krasnoyarsk atmosphere (2 m above the ground) at the beginning (<b>A</b>) and the end (<b>B</b>) of AWC.</p>
Full article ">Figure 14
<p>Simulated contributions of the main groups of emission sources in air pollution of BaP at 2 m height above the ground in the cells corresponding to the Severny observation point.</p>
Full article ">Figure 15
<p>Simulated contributions of the main groups of emission sources in air pollution of PM at 2 m height above the ground in the cells corresponding to the Pokrovka observation point.</p>
Full article ">Figure 16
<p>Simulated contributions of the main groups of emission sources in air pollution of NO<sub>2</sub> at 2 m height above the ground in the cells corresponding to the Vetluzhanka observation point.</p>
Full article ">
25 pages, 6870 KiB  
Article
Comparing Approaches for Reconstructing Groundwater Levels in the Mountainous Regions of Interior British Columbia, Canada, Using Tree Ring Widths
by Stephanie C. Hunter, Diana M. Allen and Karen E. Kohfeld
Atmosphere 2020, 11(12), 1374; https://doi.org/10.3390/atmos11121374 - 19 Dec 2020
Cited by 6 | Viewed by 2903
Abstract
Observed groundwater level records are relatively short (<100 years), limiting long-term studies of groundwater variability that could provide valuable insight into climate change effects. This study uses tree ring data from the International Tree Ring Database (ITRDB) and groundwater level data from 22 [...] Read more.
Observed groundwater level records are relatively short (<100 years), limiting long-term studies of groundwater variability that could provide valuable insight into climate change effects. This study uses tree ring data from the International Tree Ring Database (ITRDB) and groundwater level data from 22 provincial observation wells to evaluate different approaches for reconstructing groundwater levels from tree ring widths in the mountainous southern interior of British Columbia, Canada. The twenty-eight reconstruction models consider the selection of observation wells (e.g., regional average groundwater level vs. wells classified by recharge mechanism) and the search area for potential tree ring records (climate footprint vs. North American Ecoregions). Results show that if the climate footprint is used, reconstructions are statistically valid if the wells are grouped according to recharge mechanism, with streamflow-driven and high-elevation recharge-driven wells (both snowmelt-dominated) producing valid models. Of all the ecoregions considered, only the Coast Mountain Ecoregion models are statistically valid for both the regional average groundwater level and high-elevation recharge-driven systems. No model is statistically valid for low-elevation recharge-driven systems (rainfall-dominated). The longest models extend the groundwater level record to the year 1500, with the highest confidence in the later portions of the reconstructions going back to the year 1800. Full article
(This article belongs to the Special Issue Past Climate Reconstructed from Tree Rings)
Show Figures

Figure 1

Figure 1
<p>Location map of observation wells within the southern Interior Plateau Region of British Columbia and tree ring records available from the International Tree Ring Database (ITRDB). Tree ring records used for the Explorative Data Analysis (EDA) are circled in red. The inset map shows the location of the study region within western Canada and the northwestern United States. Observation wells are colour-coded by aquifer–stream system type and classified using hysteresis plots as discussed in <a href="#sec2dot1-atmosphere-11-01374" class="html-sec">Section 2.1</a>. Source for basemap digital elevation model (DEM): HydroSHEDS Database, available from <a href="http://hydrosheds.cr.usgs.gov" target="_blank">http://hydrosheds.cr.usgs.gov</a>.</p>
Full article ">Figure 2
<p>Representative hysteresis plots for (<b>a</b>) a streamflow-driven system characterized by a negative or counter-clockwise loop, and (<b>b</b>) a recharge-driven system characterized by a positive or clockwise hysteresis loop.</p>
Full article ">Figure 3
<p>Climate footprint for the BC Interior showing the spatial extent of significant correlation (<span class="html-italic">p</span> &lt; 0.10) between average depth to groundwater levels and gridded 12-month Standardized Precipitation Evaporation Index (SPEI) (Koninklijk Nederlands Meteorologisch Instituut (KNMI) Climate Explorer).</p>
Full article ">Figure 4
<p>Example groundwater (black solid line) and streamflow (shaded grey area) hydrographs for (<b>a</b>) a streamflow-driven system and b) a recharge-driven system. Temperature (grey dashed line) and precipitation (blue/light blue bars distinguishing rain from snow) are also shown in a) and (<b>b</b>). Corresponding tree ring width/depth to groundwater level correlation graphs are shown for (<b>c</b>) a streamflow-driven system and (<b>d</b>) a recharge-driven system. The standard chronologies produced using dplR are shown for the two sample chronologies which were located near these observations wells: (<b>e</b>) Cana 229 and (<b>f</b>) Cana 234, with the grey solid line indicating annual ring width indices, grey shading indicating the sample depth through time, and the solid red line showing a 20-year smoothing spline.</p>
Full article ">Figure 5
<p>Tree ring records used in the (<b>a</b>) low-elevation recharge-driven, (<b>b</b>) streamflow-driven, and (<b>c</b>) high-elevation recharge-driven groundwater level models, using the climate footprint to select potential tree ring records.</p>
Full article ">Figure 6
<p>Final groundwater level reconstructions and actual groundwater levels (GWL, shown in dark blue) for the average of wells in the (<b>a</b>) streamflow-driven and (<b>b</b>) high-elevation recharge-driven models using the climate footprint, and the (<b>c</b>) all-wells and (<b>d</b>) high-elevation recharge-driven models created using the Coast Mountain Ecoregions. The number of tree ring records used in each model throughout time is shown with the dashed red line. The black arrow indicates the time at which the reconstruction models might be truncated due to increasing uncertainty moving farther back into time (see <a href="#sec4-atmosphere-11-01374" class="html-sec">Section 4</a>).</p>
Full article ">
19 pages, 1528 KiB  
Review
Monitoring of Selected CBRN Threats in the Air in Industrial Areas with the Use of Unmanned Aerial Vehicles
by Anna Rabajczyk, Jacek Zboina, Maria Zielecka and Radosław Fellner
Atmosphere 2020, 11(12), 1373; https://doi.org/10.3390/atmos11121373 - 19 Dec 2020
Cited by 11 | Viewed by 4926
Abstract
Unmanned aerial vehicles (UAVs) play an increasingly important role in various areas of life, including in terms of protection and security. As a result of fires, volcanic eruptions, or other emergencies, huge amounts of toxic gases, dust, and other substances are emitted into [...] Read more.
Unmanned aerial vehicles (UAVs) play an increasingly important role in various areas of life, including in terms of protection and security. As a result of fires, volcanic eruptions, or other emergencies, huge amounts of toxic gases, dust, and other substances are emitted into the environment, which, together with high temperature, often leads to serious environmental contamination. Based on the available literature and patent databases, an analysis of the available UAVs models was carried out in terms of their applicability in air contaminated conditions in industrial areas, in the event of emergencies, such as fire, chemical contamination. The possibilities of using the devices were analyzed in terms of weather conditions, construction, and used materials in CBRN (chemical, biological, radiological, nuclear) threat situations. It was found that, thanks to the use of appropriate sensors, cameras, and software of UAVs integrated with a given system, it is possible to obtain information on air quality at a given moment, which is very important for the safety of people and the environment. However, several elements, including the possibility of use in acidification conditions, requires refinement to changing crisis conditions. Full article
(This article belongs to the Section Air Quality)
Show Figures

Figure 1

Figure 1
<p>The place of the incident during the exercises in Dortmund during the e-Notice project. The place of the simulated leakage is marked with a red circle (21.09.2019 Radoslaw Fellner).</p>
Full article ">Figure 2
<p>Comparison of images from RGB and thermal imaging cameras—a railway wagon tank (21.09.2019 Radoslaw Fellner).</p>
Full article ">
13 pages, 1991 KiB  
Article
Improving Air Quality for Operators of Mobile Machines in Underground Mines
by Andrzej Szczurek, Monika Maciejewska, Marcin Przybyła and Wacław Szetelnicki
Atmosphere 2020, 11(12), 1372; https://doi.org/10.3390/atmos11121372 - 18 Dec 2020
Cited by 1 | Viewed by 2941
Abstract
In underground mines, mobile mining equipment is critical for the production system. The microenvironment inside the mobile machine may cause exposure to strongly polluted mine air, which adversely affects the health and working performance of the operator. Harmful pollutants may access the cabin [...] Read more.
In underground mines, mobile mining equipment is critical for the production system. The microenvironment inside the mobile machine may cause exposure to strongly polluted mine air, which adversely affects the health and working performance of the operator. Harmful pollutants may access the cabin together with the ventilation air delivered from the machine’s surroundings. This work proposes a solution that is able to ensure that the air for the machine operator is of proper quality. The proposal emerged from an analysis of the compliance of cabins of mobile machines working underground in mines with occupational health and safety (H&S) standards. An analytical model of air quality in a well-mixed zone was utilized for this purpose. The cabin atmosphere was investigated with regard to the concentration of gaseous species in the surrounding air, the cabin ventilation rate, and human breathing parameters. The analysis showed that if currently available ventilation approaches are used, compliance with multiple H&S standards cannot be attained inside the cabin if standards are exceeded in the surroundings of the machine. The proposed solution overcomes this problem by combining elements that are already in place, i.e., ventilation, air-conditioning, and filtration with a personal supply of clean air. The concept is generic and may be adapted to various specific requirements. Full article
(This article belongs to the Special Issue Exposure Assessment of Air Pollution)
Show Figures

Figure 1

Figure 1
<p>Cabin of a loader—autonomous mobile machine operating in underground mines. Internal dimensions: A—1275 mm, C—1335 mm, D—800 mm.</p>
Full article ">Figure 2
<p>Time when the concentration of a substance in cabin air increases from the clean air level to the predefined threshold value. The time is displayed as a function of the air exchange rate between the cabin and its surroundings, for various concentrations of a substance in the mine air. The concentration in the mine air refers to the threshold level multiple, using <math display="inline"><semantics> <mi>b</mi> </semantics></math>. Solid lines refer to substances emitted only outside the cabin, such as NO, NO<sub>2</sub>, and CO. Dotted lines refer to CO<sub>2</sub>.</p>
Full article ">Figure 3
<p>H&amp;S standard compliance in cabin air with respect to substances emitted exclusively outside cabin, such as NO, NO<sub>2</sub>, and CO.</p>
Full article ">Figure 4
<p>H&amp;S standard compliance in cabin air with respect to CO<sub>2</sub>.</p>
Full article ">Figure 5
<p>Isolines display the fraction of time when the threshold concentration may be exceeded in the cabin air (no more than threshold multiple), while in the remaining time the concentration remains below the threshold (no more than the threshold fraction). This ensures H&amp;S standard compliance.</p>
Full article ">Figure 6
<p>Isolines display the fraction of the clean air, <math display="inline"><semantics> <mrow> <msub> <mi>f</mi> <mrow> <mi>P</mi> <mi>U</mi> <mi>R</mi> </mrow> </msub> <mo>,</mo> </mrow> </semantics></math> that has to be delivered to the breathing zone of the mobile machine operator, to ensure the required quality of the breathing air (threshold fraction) for the particular composition of the cabin air (threshold multiple). Solid lines refer to the substances emitted only outside the cabin such as NO, NO<sub>2</sub>, and CO. Dotted lines refer to CO<sub>2</sub>.</p>
Full article ">
12 pages, 2079 KiB  
Article
Vertical Profiles of Atmospheric Species Concentrations and Nighttime Boundary Layer Structure in the Dry Season over an Urban Environment in Central Amazon Collected by an Unmanned Aerial Vehicle
by Patrícia Guimarães, Jianhuai Ye, Carla Batista, Rafael Barbosa, Igor Ribeiro, Adan Medeiros, Tianning Zhao, Wei-Chun Hwang, Hui-Ming Hung, Rodrigo Souza and Scot T. Martin
Atmosphere 2020, 11(12), 1371; https://doi.org/10.3390/atmos11121371 - 18 Dec 2020
Cited by 22 | Viewed by 3704
Abstract
Nighttime vertical profiles of ozone, PM2.5 and PM10 particulate matter, carbon monoxide, temperature, and humidity were collected by a copter-type unmanned aerial vehicle (UAV) over the city of Manaus, Brazil, in central Amazon during the dry season of 2018. The vertical [...] Read more.
Nighttime vertical profiles of ozone, PM2.5 and PM10 particulate matter, carbon monoxide, temperature, and humidity were collected by a copter-type unmanned aerial vehicle (UAV) over the city of Manaus, Brazil, in central Amazon during the dry season of 2018. The vertical profiles were analyzed to understand the structure of the urban nighttime boundary layer (NBL) and pollution within it. The ozone concentration, temperature, and humidity had an inflection between 225 and 350 m on most nights, representing the top of the urban NBL. The profile of carbon monoxide concentration correlated well with the local evening vehicular congestion of a modern transportation fleet, providing insight into the surface-atmosphere dynamics. In contrast, events of elevated PM2.5 and PM10 concentrations were not explained well by local urban emissions, but rather by back trajectories that intersected regional biomass burning. These results highlight the potential of the emerging technologies of sensor payloads on UAVs to provide new constraints and insights for understanding the pollution dynamics in nighttime boundary layers in urban regions. Full article
(This article belongs to the Special Issue Atmospheric Measurements Using Unmanned Systems)
Show Figures

Figure 1

Figure 1
<p>Representative image of the study area in the city of Manaus, Brazil, in central Amazon. The urban region is crossed by rivers and interspersed with forests. Nighttime flights to collect vertical profiles took place on the campus of the Amazonas State University (marked as a red disc and labeled “UEA”) in the middle of the city. The green region in the inset shows the state of Amazonas, Brazil, in South America.</p>
Full article ">Figure 2
<p>Box-whisker plots of the statistics of ozone concentration at different altitudes. Results are presented for the (<b>a</b>) wet and (<b>b</b>) dry seasons of 2018 (57 and 84 flights, respectively). The red line represents the median, the box edges show the quartiles, and the horizontal black lines indicate the minimum and maximum values, excluding outliers. Flights took place between 20:00 and 00:00 (local time (LT)). LT was 4 h earlier than UTC. The data of panel (<b>a</b>) are from Guimarães et al. [<a href="#B10-atmosphere-11-01371" class="html-bibr">10</a>].</p>
Full article ">Figure 3
<p>Box-whisker statistics plots of the height of the nighttime boundary layer during the dry season of 2018 based on (<b>a</b>) ozone concentration, (<b>b</b>) potential temperature, and (<b>c</b>) specific humidity (84 flights). The median (red line), quartiles (blue box edges), and the minimum and maximum values (black lines) are represented, excluding outliers. Flights took place between 20:00 and 00:00 (LT). Local time (LT) was 4 h earlier than UTC.</p>
Full article ">Figure 4
<p>Vertical profiles of the concentrations of ozone, carbon monoxide, PM<sub>2.5</sub>, and PM<sub>10</sub> from land surface to an altitude of 500 m for 8 flights on 29 September 2018 between 20:00 and 00:00 (LT). The red line represents the median of the data at each altitude. The dashed line represents the instrumental limit of detection for ozone (3 ppbv). Individual profiles are plotted in <a href="#app1-atmosphere-11-01371" class="html-app">Figure S2</a>.</p>
Full article ">Figure 5
<p>Evolution of vertical profiles of the concentrations of ozone, carbon monoxide, PM<sub>2.5</sub>, and PM<sub>10</sub> from land surface to an altitude of 500 m on 27 September 2018 from 20:00 to 22:00 (LT).</p>
Full article ">
16 pages, 3752 KiB  
Article
Local Weather Conditions Create Structural Differences between Shallow Firn Columns at Summit, Greenland and WAIS Divide, Antarctica
by Ian E. McDowell, Mary R. Albert, Stephanie A. Lieblappen and Kaitlin M. Keegan
Atmosphere 2020, 11(12), 1370; https://doi.org/10.3390/atmos11121370 - 17 Dec 2020
Cited by 2 | Viewed by 3376
Abstract
Understanding how physical characteristics of polar firn vary with depth assists in interpreting paleoclimate records and predicting meltwater infiltration and storage in the firn column. Spatial heterogeneities in firn structure arise from variable surface climate conditions that create differences in firn grain growth [...] Read more.
Understanding how physical characteristics of polar firn vary with depth assists in interpreting paleoclimate records and predicting meltwater infiltration and storage in the firn column. Spatial heterogeneities in firn structure arise from variable surface climate conditions that create differences in firn grain growth and packing arrangements. Commonly, estimates of how these properties change with depth are made by modeling profiles using long-term estimates of air temperature and accumulation rate. Here, we compare surface meteorology and firn density and permeability in the depth range of 3.5–11 m of the firn column from cores collected at Summit, Greenland and WAIS Divide, Antarctica, two sites with the same average accumulation rate and mean annual air temperature. We show that firn at WAIS Divide is consistently denser than firn at Summit. However, the difference in bulk permeability of the two profiles is less statistically significant. We argue that differences in local weather conditions, such as mean summer temperatures, daily temperature variations, and yearly wind speeds, create the density discrepancies. Our results are consistent with previous results showing density is not a good indicator of firn permeability within the shallow firn column. Future modeling efforts should account for these weather variables when estimating firn structure with depth. Full article
(This article belongs to the Special Issue Modeling and Measuring Snow Processes across Scales)
Show Figures

Figure 1

Figure 1
<p>Location of firn core collection and automated weather station sites: (<b>a</b>) Summit, Greenland (GIS Summit) and (<b>b</b>) West Antarctic Ice Sheet (WAIS) Divide, Antarctica (WAIS Divide). Coring locations are plotted on the Greenland 5 km digital elevation model (DEM) [<a href="#B18-atmosphere-11-01370" class="html-bibr">18</a>] and the Antarctic 1 km DEM [<a href="#B19-atmosphere-11-01370" class="html-bibr">19</a>]. Greenland and Antarctica are displayed using the WGS 1984 Arctic and Antarctic Polar Stereographic Coordinate systems, respectively. The Summit firn core was collected from ~3200 m above sea level and the elevation at WAIS Divide was ~1800 m.</p>
Full article ">Figure 2
<p>Firn density profiles (<b>a</b>,<b>b</b>) and permeability profiles (<b>c</b>,<b>d</b>) between 3.5 m and 11 m depth at GIS Summit (<b>a</b>,<b>c</b>) and WAIS Divide (<b>b</b>,<b>d</b>). The solid black line in panels (<b>a</b>,<b>b</b>) is the modeled firn density profile generated from Herron and Langway (1980) [<a href="#B12-atmosphere-11-01370" class="html-bibr">12</a>]. In panels (<b>c</b>,<b>d</b>), the solid black line represents the modeled permeability profile from Adolph and Albert (2014) [<a href="#B14-atmosphere-11-01370" class="html-bibr">14</a>] and the dashed black line shows modeled permeability profiles using the model of Freitag et al. (2002) [<a href="#B15-atmosphere-11-01370" class="html-bibr">15</a>]. All modeled density and permeability profiles shown were generated using an accumulation rate of 0.22 m yr<sup>−1</sup>, mean air temperature of −28 °C, and surface snow density of 0.34 kg m<sup>−3</sup>.</p>
Full article ">Figure 3
<p>Distribution of density (<b>a</b>,<b>c</b>) and permeability (<b>b</b>,<b>d</b>) in the depth range of 3.5–11 m in the firn column at GIS Summit (<b>a</b>,<b>b</b>) and WAIS Divide (<b>c</b>,<b>d</b>). The number of measurements, <span class="html-italic">n</span>, used to create each histogram are as follows: GIS Summit density <span class="html-italic">n</span> = 137, GIS Summit permeability <span class="html-italic">n</span> = 137, WAIS Divide density <span class="html-italic">n</span> = 86, WAIS Divide permeability <span class="html-italic">n</span> = 71. Bin widths were set to 0.2 kg m<sup>−3</sup> for density histograms and 7.0 × 10<sup>−10</sup> m<sup>2</sup> for permeability histograms. Box plots above each distribution show interquartile range and whiskers extend to encompass 95% of the data. The median value is indicated by a line within the interquartile range.</p>
Full article ">Figure 4
<p>Meteorological data from both Summit (panels (<b>a</b>,<b>c</b>,<b>e</b>)) and WAIS Divide (panels (<b>b</b>,<b>d</b>,<b>f</b>). Air temperatures at the two sites are shown in panels (<b>a</b>,<b>b</b>), relative humidity values are displayed in panels (<b>c</b>,<b>d</b>), and wind speeds are plotted in panels (<b>e</b>,<b>f</b>). Hourly data are plotted from the years 1997–2007 for Summit and 2009–2019 for WAIS Divide vs. the day of the collection year with grey transparent lines, so that progressively darker greys show overlapping data. Black lines show a mean meteorological value for the given day of the year, averaged over the collection period.</p>
Full article ">Figure 5
<p>Variability of daily air temperatures at Summit (<b>a</b>) and WAIS Divide (<b>b</b>). Averaged daily minimum and maximum temperatures are plotted, with the shaded area representing the daily temperature variability. The difference between the highest and lowest daily temperatures are plotted in panel (<b>c</b>).</p>
Full article ">Figure A1
<p>NCEP-NCAR reanalysis temperature and relative humidity (grey) plotted with 30-day running means of AWS observed variables (black) for Summit, Greenland (<b>a</b>,<b>c</b>) and WAIS Divide, Antarctica (<b>b</b>,<b>d</b>). Dashed grey lines show least-squares regression fits to the reanalysis data, showing that meteorological conditions remained essentially constant over the previous 25 years. Discrepancies between observed data and reanalysis data are attributed to local variability not captured in the regional nature of the 2.5° × 2.5° gridded reanalysis output.</p>
Full article ">
22 pages, 4312 KiB  
Article
The Use of Gaussian Mixture Models with Atmospheric Lagrangian Particle Dispersion Models for Density Estimation and Feature Identification
by Alice Crawford
Atmosphere 2020, 11(12), 1369; https://doi.org/10.3390/atmos11121369 - 17 Dec 2020
Cited by 10 | Viewed by 5001
Abstract
Atmospheric Lagrangian particle dispersion models, LPDM, simulate the dispersion of passive tracers in the atmosphere. At the most basic level, model output consists of the position of computational particles and the amount of mass they represent. In order to obtain concentration values, this [...] Read more.
Atmospheric Lagrangian particle dispersion models, LPDM, simulate the dispersion of passive tracers in the atmosphere. At the most basic level, model output consists of the position of computational particles and the amount of mass they represent. In order to obtain concentration values, this information is then converted to a mass distribution via density estimation. To date, density estimation is performed with a nonparametric method so that output consists of gridded concentration data. Here we introduce the use of Gaussian mixture models, GMM, for density estimation. We compare to the histogram or bin counting method for a tracer experiment and simulation of a large volcanic ash cloud. We also demonstrate the use of the mixture model for automatic identification of features in a complex plume such as is produced by a large volcanic eruption. We conclude that use of a mixture model for density estimation and feature identification has potential to be very useful. Full article
(This article belongs to the Special Issue Forecasting the Transport of Volcanic Ash in the Atmosphere)
Show Figures

Figure 1

Figure 1
<p>The top panels, (<b>a</b>–<b>c</b>) show particle positions at (<b>a</b>,<b>b</b>) 09/18/1983 21:00 UTC and (<b>c</b>) 09/18/1983 from 21 to 24 UTC with particle positions output every 5 min. The Gaussian fits are shown in light blue and the colors of the particles indicate which cluster the particles were assigned to by the algorithm. The bottom panels, (<b>d</b>–<b>f</b>) show 3 h averaged concentrations calculated from the fits on a <math display="inline"> <semantics> <mrow> <msup> <mn>0.05</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.05</mn> <mo>∘</mo> </msup> </mrow> </semantics> </math> 25 m grid from the ground to 25 m. A threshold of <math display="inline"> <semantics> <mrow> <mn>1</mn> <mspace width="3.33333pt"/> <mi>pg</mi> <mspace width="3.33333pt"/> <msup> <mi mathvariant="normal">m</mi> <mn>3</mn> </msup> </mrow> </semantics> </math> has been applied and station measurements for the CAPTEX1 experiment are shown by the circles. (<b>d</b>,<b>e</b>) were calculated from separate fits to each time period and the panels above show the first such fit. (<b>f</b>) was calculated from one fit to all time periods which is shown in (<b>c</b>).</p>
Full article ">Figure 2
<p>Comparison of density estimations. Color shows concentrations in pg m<math display="inline"> <semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> </semantics> </math>. Shown is a cross section of the simulated plume from the CAPTEX 1 experiment with station measurements shown by the circles. Simulated concentrations are 3 h averages on 09/18/1983 from 21 to 24 UTC, 3 h after the release. all figures have horizontal resolution of <math display="inline"> <semantics> <mrow> <msup> <mn>0.05</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.05</mn> <mo>∘</mo> </msup> <mspace width="3.33333pt"/> <mo>×</mo> </mrow> </semantics> </math> 25 m except for the control runE, (<b>d</b>) has resolution of <math display="inline"> <semantics> <mrow> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> </mrow> </semantics> </math>100 m. Number of particles released per hour is shown at the top of each column. (<b>a</b>–<b>e</b>) are calculated using histogram method. (<b>f</b>–<b>k</b>) are calculated using a GMM with number of components annotated on the graph. Particles from the ground to 500 m were used in the fit. A threshold of 10 pg m<math display="inline"> <semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> </semantics> </math> was applied to (<b>f</b>–<b>k</b>). (<b>i</b>,<b>j</b>) were calculated from fits to each time step in the averaging period, while (<b>f</b>,<b>g</b>,<b>h</b>,<b>k</b>) were calculated from a fit to all particles in the averaging period. (<b>h</b>) and (<b>k</b>) were calculated from Run C.</p>
Full article ">Figure 3
<p>Statitics for CAPTEX experiments using a variety of density reconstructions. The rank, described by Equation (<a href="#FD1-atmosphere-11-01369" class="html-disp-formula">1</a>) is shown in panel (<b>a</b>). The four statistics which go into calculating rank are shown in panels (<b>b</b>–<b>e</b>) and the root mean square error, RMSE, is shown in (<b>f</b>). Values for the standard simulation, runE, with standard concentration grid are shown by the black square. Values for concentrations calculated with the histogram and a high resolution grid are plotted slightly to the left of the black square. Values for concentrations calculated with the mixture model are plotted to the right of the black square. The high resolution grid of <math display="inline"> <semantics> <mrow> <msup> <mn>0.05</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.05</mn> <mo>∘</mo> </msup> <mo>×</mo> </mrow> </semantics> </math> 25 m was used unless otherwise noted in the legend. For the GMM, the time averaging was performed with a separate fit to each output time unless tmave = aggregate is noted in the legend.</p>
Full article ">Figure 4
<p>Mass loading from simulations of the 2008 eruption of Kasatochi. The first and third rows (<b>a</b>–<b>c</b>,<b>f</b>,<b>h</b>,<b>i</b>) display plots of mass loading with the simulation used and method of density reconstruction noted in the bottom right corner. The GMM used 50 Gaussians and output displayed at the same resolution as the histogram plots (<math display="inline"> <semantics> <mrow> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> </mrow> </semantics> </math> grid). Directly below each mass loading plot is a histogram of the mass loading values with a threshold of 0.01 g m<math display="inline"> <semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>2</mn> </mrow> </msup> </semantics> </math> applied. The histogram for (<b>a</b>) is displayed as the red shaded region in each histogram plot for easy comparison.</p>
Full article ">Figure 5
<p>Concentrations at longitude-151. Top row (<b>a</b>–<b>c</b>) concentrations calculated using histogram method resolution <math display="inline"> <semantics> <mrow> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> <mn>1000</mn> <mspace width="3.33333pt"/> <mi mathvariant="normal">m</mi> </mrow> </semantics> </math>. Second row (<b>d</b>–<b>f</b>) and third (<b>g</b>–<b>i</b>) row concentrations calculated using GMM with 50 Gaussians and plotted at resolution of <math display="inline"> <semantics> <mrow> <msup> <mn>0.1</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.1</mn> <mo>∘</mo> </msup> <mo>×</mo> <mn>200</mn> <mspace width="0.277778em"/> <mi mathvariant="normal">m</mi> </mrow> </semantics> </math>. For the second row concentrations below below 0.01 g m<math display="inline"> <semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> </semantics> </math> are not shown. The third row shows contour levels of output shown in second row and the black dots indicate position of computational particles within <math display="inline"> <semantics> <mrow> <msup> <mn>0.125</mn> <mo>∘</mo> </msup> </mrow> </semantics> </math> of the −151 line of longitude. The dots are shown smaller in (<b>g</b>) for clarity. First column (<b>a</b>,<b>d</b>,<b>g</b>) Run KA, second column (<b>b</b>,<b>e</b>,<b>h</b>) Run KB and third column (<b>c</b>,<b>f</b>,<b>i</b>) run KD.</p>
Full article ">Figure 6
<p>Concentrations at longitude −145. Top row (<b>a</b>–<b>c</b>) concentrations calculated using histogram method. Second row (<b>d</b>–<b>f</b>) and third (<b>g</b>–<b>i</b>) row concentrations calculated using GMM with 50 Gaussians and plotted at the same resolution as the histogram method in the second row and the contours in the third row are for with a higher vertical resolution (<math display="inline"> <semantics> <mrow> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> </mrow> </semantics> </math> 100 m). First column (<b>a</b>,<b>d</b>,<b>g</b>) Run KA, second column (<b>b</b>,<b>e</b>,<b>h</b>) Run KB and third column (<b>c</b>,<b>f</b>,<b>i</b>) run KD. The black dots in the third row indicate position of particles within <math display="inline"> <semantics> <mrow> <msup> <mn>0.125</mn> <mo>∘</mo> </msup> </mrow> </semantics> </math> of the −145 line of longitude. The red dots indicate position of particles within <math display="inline"> <semantics> <mrow> <msup> <mn>0.5</mn> <mo>∘</mo> </msup> </mrow> </semantics> </math> of the −145 line of longitude.</p>
Full article ">Figure 7
<p>The score, S<math display="inline"> <semantics> <msub> <mrow/> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </semantics> </math>, as a function of time. The score is described in the text.</p>
Full article ">Figure 8
<p>log probability of each particle position. The left column shows a projection onto the latitude, longitude plane and the right column shows a projection onto the height, longitude plane. The date is noted at the top of each figure. <math display="inline"> <semantics> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </semantics> </math> indicates that the probabilities were calculated using the particle positions <span class="html-italic">j</span> and the fit to particle positions <span class="html-italic">i</span>. For example, <math display="inline"> <semantics> <msub> <mi>S</mi> <mn>41</mn> </msub> </semantics> </math> shows positions of the smallest particle size with how probable it is that they could belong to a fit to the largest particle size. Note the differing color scales. Low values indicate poor fit and regions where the overlap between the location of the two particle sizes is low.</p>
Full article ">Figure 9
<p>The top two rows show features identified by a BGMM with 10 Gaussians at three times. The black stars show the centers of the Gaussians. The colored dots indicate the position of the computational particles. The different colors belong to different groups. The width of each Gaussian is indicated by the blue shaded regions. The top row is a projection onto the height-longitude plane while the second row is a projection onto the latitude-longitude plane. The bottom row shows the position of the centers of the 10 Gaussians every hour from 10 August 12 UTC to 11 August 11 UTC. The plot in the left corner is a projection onto the latitude-longitude plane while the plot on the right is a three dimensional rendering of the positions.</p>
Full article ">Figure A1
<p>Examples of shot noise in the histogram method. The top row (<b>a</b>–<b>c</b>) shows normalized histograms (blue bars) of the number of points found in a volume. The black line shows the Poisson distribution with the same mean. The position of computational particles in the volume are shown in the second and third rows. Positions from three of the simulations are plotted with the different colors indicating the computational particles from different simulations. The first and second columns (<b>a</b>,<b>b</b>,<b>d</b>,<b>e</b>,<b>g</b>,<b>h</b>) are from 50 simulations of the runKB with output every 10 min from 9 August, 12:10 UTC to 13:00 UTC. The first column (<b>a</b>,<b>d</b>,<b>g</b>) looks at a smaller volume (<math display="inline"> <semantics> <mrow> <msup> <mn>0.1</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.1</mn> <mo>∘</mo> </msup> <mo>×</mo> </mrow> </semantics> </math> 500 m) and only one particle size while the second column uses a volume of <math display="inline"> <semantics> <mrow> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> </mrow> </semantics> </math> 1 km and all the particle sizes. The last column (<b>c</b>,<b>f</b>,<b>i</b>) is from 100 simulations of runB with output every 5 min from 25 September 18:00 UTC to 21:00 UTC and a volume of <math display="inline"> <semantics> <mrow> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> <msup> <mn>0.25</mn> <mo>∘</mo> </msup> <mo>×</mo> </mrow> </semantics> </math> 100 m. Although the histogram in (<b>c</b>) looks somewhat flatter than the Poisson, this may be simply because 100 points is not enough to represent the distribution well.</p>
Full article ">Figure A2
<p>Score, <math display="inline"> <semantics> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </semantics> </math> as a function of number of Gaussians used in the GMM fit. Black and green lines indicate score for <math display="inline"> <semantics> <mrow> <mi>i</mi> <mo>=</mo> <mi>j</mi> </mrow> </semantics> </math> and red and blue lines indicate scores for <math display="inline"> <semantics> <mrow> <mi>i</mi> <mo>≠</mo> <mi>j</mi> </mrow> </semantics> </math>. (<b>a</b>) RunC and RunD for CAPTEX1 at 09/20/1983 0 UTC. (<b>b</b>) RunC and RunD for CAPTEX1 at 09/20/1983 03 UTC. (<b>c</b>) RunKC and runKD at 08/09/2008 at 04 UTC. (<b>d</b>) RunKC and RunKD at 08/10/2008 12 UTC.</p>
Full article ">
16 pages, 4819 KiB  
Article
Carbonaceous Aerosols in PM1, PM2.5, and PM10 Size Fractions over the Lanzhou City, Northwest China
by Xin Zhang, Zhongqin Li, Feiteng Wang, Mengyuan Song, Xi Zhou and Jing Ming
Atmosphere 2020, 11(12), 1368; https://doi.org/10.3390/atmos11121368 - 17 Dec 2020
Cited by 30 | Viewed by 3552
Abstract
Carbonaceous particles have been confirmed as major components of ambient aerosols in urban environments and are related to climate impacts and environmental and health effects. In this study, we collected different-size particulate matter (PM) samples (PM1, PM2.5, and PM [...] Read more.
Carbonaceous particles have been confirmed as major components of ambient aerosols in urban environments and are related to climate impacts and environmental and health effects. In this study, we collected different-size particulate matter (PM) samples (PM1, PM2.5, and PM10) at an urban site in Lanzhou, northwest China, during three discontinuous one-month periods (January, April, and July) of 2019. We measured the concentrations and potential transport pathways of carbonaceous aerosols in PM1, PM2.5, and PM10 size fractions. The average concentrations of OC (organic carbon) and EC (elemental carbon) in PM1, PM2.5, and PM10 were 6.98 ± 3.71 and 2.11 ± 1.34 μg/m3, 8.6 ± 5.09 and 2.55 ± 1.44 μg/m3, and 11.6 ± 5.72 and 4.01 ± 1.72 μg/m3. The OC and EC concentrations in PM1, PM2.5, and PM10 had similar seasonal trends, with higher values in winter due to the favorable meteorology for accumulating pollutants and urban-increased emissions from heating. Precipitation played a key role in scavenge pollutants, resulting in lower OC and EC concentrations in summer. The OC/EC ratios and principal component analysis (PCA) showed that the dominant pollution sources of carbon components in the PMs in Lanzhou were biomass burning, coal combustion, and diesel and gasoline vehicle emissions; and the backward trajectory and concentration weight trajectory (CWT) analysis further suggested that the primary pollution source of EC in Lanzhou was local fossil fuel combustion. Full article
(This article belongs to the Special Issue Chemistry of Aqueous Surfaces in the Atmospheric Context)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Geographical location of Lanzhou, China and (<b>b</b>) locations of the sampling site on the topographic map. Background map was retrieved from Google Earth.</p>
Full article ">Figure 2
<p>Time series of the EC and OC concentrations in PM<sub>1</sub>, PM<sub>2.5</sub>, and PM<sub>10</sub>, respectively. The temperature, relative humidity (RH), precipitation, and wind speed (WS) and direction (WD) were recorded by the automatic weather station.</p>
Full article ">Figure 3
<p>Seasonal variations of OC and EC in the PMs, where each box indicates the 99%, 75%, 50%, 25%, and 1% quartiles of the data from the top to the bottom.</p>
Full article ">Figure 4
<p>Seasonal variation of POC, SOC, OC/EC, and POC/SOC in different-sized PMs.</p>
Full article ">Figure 5
<p>Regressions between OC and EC in PM<sub>1</sub>, PM<sub>2.5</sub>, and PM<sub>10.</sub> Different lines (colour) represent the Regressions between in different seasons.</p>
Full article ">Figure 6
<p>Seasonal variation and size distribution of carbon components.</p>
Full article ">Figure 7
<p>Five-day backward trajectories of Lanzhou and the ratio between OC and EC in each cluster in 500-m layer in (<b>a</b>) winter, (<b>b</b>) spring, and (<b>c</b>) summer. Different lines (colour) represent different clusters.</p>
Full article ">Figure 8
<p>Concentration-weighted trajectory (CWT) of OC in (<b>a</b>) winter, (<b>b</b>) spring, and (<b>c</b>) summer and of EC in (<b>d</b>) winter, (<b>e</b>) spring, and (<b>f</b>) summer in PM<sub>10</sub> in Lanzhou.</p>
Full article ">
29 pages, 10275 KiB  
Article
Rainfall Threshold for Shallow Landslides Initiation and Analysis of Long-Term Rainfall Trends in a Mediterranean Area
by Anna Roccati, Guido Paliaga, Fabio Luino, Francesco Faccini and Laura Turconi
Atmosphere 2020, 11(12), 1367; https://doi.org/10.3390/atmos11121367 - 17 Dec 2020
Cited by 30 | Viewed by 4949
Abstract
The effects of climate change on landslide activity may have important environmental, socio-economic, and political consequences. In the last decades, several short-term extreme rainfall events affected Mediterranean regions, resulted in damaging geo-hydrological processes and casualties. It is unequivocal that the impact of landslides [...] Read more.
The effects of climate change on landslide activity may have important environmental, socio-economic, and political consequences. In the last decades, several short-term extreme rainfall events affected Mediterranean regions, resulted in damaging geo-hydrological processes and casualties. It is unequivocal that the impact of landslides in several Mediterranean countries is increasing with time, but until now, there has been little or no quantitative data to support these increases. In this paper, both rainfall conditions for the occurrence of shallow landslides and rainfall trends were investigated in the Portofino promontory, which extends in the Ligurian Sea, where heavy rainfall and related ground effects often occur. Adopting a frequentist approach, the empirical intensity-duration threshold was estimated. Our findings highlight that the rainfall intensity required to trigger landslides is lower for the same duration than those expected in other similar environments, suggesting a high susceptibility to rainfall-induced landslides in the Portofino territory. Further, the Mann-Kendall test and Hurst exponent were used for detecting potential trends. Analysis of long-term rainfall time series showed statistically significant increasing trends in short duration precipitation occurrence and rainfall rates, suggesting a possible future scenario with a more frequent exceedance of the threshold triggering value and an increase of landslide risk. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the Portofino promontory. The red line shows the study area; white lines represent the boundaries of the three municipalities.</p>
Full article ">Figure 2
<p>Geographical setting of the Portofino promontory.</p>
Full article ">Figure 3
<p>Geological and geomorphological sketch map of the Portofino promontory: 1. Anthropic deposits; 2. Alluvial deposits; 3. Debris covers; 4. Conglomerate of Portofino; 5. Flysch of Mt. Antola; 6. Fault and presumed fault; 7. Fold axis; 8. Areas affected by rockfalls and widespread shallow landslides; 9. Deep-seated gravitational slope deformation; 10. Downcutting talwegs; 11. Cliff; 12. Active landslide or degradation scarp; 13. Inactive landslide or degradation scarp.</p>
Full article ">Figure 4
<p>Rain gauges available in the Portofino promontory and surrounding area classified by networks and type of rainfall data. Colored symbols represent different types of rainfall data: red, daily; blue, hourly; green, daily (before 2003) and hourly (after 2003). For rain gauge list, see <a href="#atmosphere-11-01367-t002" class="html-table">Table 2</a>.</p>
Full article ">Figure 5
<p>(<b>A</b>) Spatial distribution of the rainfall-induced shallow landslides and mud-debris flows that affected the Portofino promontory over the 1910–2019 period. (<b>B</b>) Monthly distribution of slope instabilities triggered by precipitations over the 1910–2019 period. (<b>C</b>) Distribution of annual landslide occurrences over the 1910–2019 period.</p>
Full article ">Figure 6
<p>(<b>A</b>) Debris flow affected slope downstream San Rocco along the western coast on 28 October 1961 (photo R. Terranova); (<b>B</b>) New debris flow occurred on 24 March 1964, in the same coastal sector previously involved in instability processes (photo R. Terranova); (<b>C</b>) Erosional processes and widespread shallow landslides along the slope at Le Gave, after the events in 1987, 1995 and 1996 (photo R. Bovolenta); (<b>D</b>) Warning notice for possible rockfall along the road crossing the ancient coastal landslide body at Le Gave (photo R. Bovolenta). See location in <a href="#atmosphere-11-01367-f002" class="html-fig">Figure 2</a>.</p>
Full article ">Figure 7
<p>(<b>A</b>) San Fruttuoso bay and its hamlet in 1910 when the sea lapped against rocks and the ancient arcade (photo from the historical archive of Agenzia Bozzo, Camogli); (<b>B</b>) The new small beach in front of the abbey in 1920 generated by debris and materials carried downstream the Fontanini creek by running water during the damaging rainfall event on 25 September 1915 (photo from the historical archive of Agenzia Bozzo, Camogli); (<b>C</b>) Suspended sediment load inside the bay transported downstream the San Fruttuoso creek by intense and short-lived rainfall on 26 July 2014; (<b>D</b>) Damage to tourist structures and facilities in San Fruttuoso hamlet caused by the mud-debris flow occurred on 26 July 2014.</p>
Full article ">Figure 8
<p>Rainfall thresholds obtained for the possible initiation of shallow landslides and mud/debris flows on the Portofino promontory corresponding to a 5% exceedance probability level (red curve) and associated uncertainty (grey pattern).</p>
Full article ">Figure 9
<p>Comparison between the ID threshold established for the Portofino promontory and threshold curves available in the literature for the possible initiation of shallow landslides in Italy and worldwide (<a href="#atmosphere-11-01367-t004" class="html-table">Table 4</a>). The red curve is the new threshold obtained in the present work. Black curves represent global thresholds, dotted and broken black curves portray local and regional thresholds, respectively, grey curves symbolize thresholds defined for single catchments.</p>
Full article ">Figure 10
<p>(<b>A</b>) Rainfall conditions (gray dots) that resulted in shallow landslides on Portofino promontory: colored dots show rainfall conditions for the 7 September 2020 (red) and the 1–2 October 2020 (yellow) events; the black line represents the 5% rainfall threshold. (<b>B</b>) Hourly rainfall (blue bars) and cumulated event rainfall (red line) for the 7 September 2020 recorded at Santa Margherita rain gauge. (<b>C</b>) Hourly rainfall (blue bars) and cumulated event rainfall (red line) for the 1–2 October 2020 recorded at Rapallo rain gauge. Red arrows identify the time of landslides occurrence (true or assumed).</p>
Full article ">Figure 11
<p>Rainy days (<b>A</b>) and rainfall rate (<b>B</b>) series scatterplot for Genoa University rain gauge (1833–2019); MM30: mean mobile at 30 years; S: Sen’s trend line; Sl: Sen’s lower confidence (95%) slope; Su: Sen’s upper confidence (95%) slope.</p>
Full article ">Figure 11 Cont.
<p>Rainy days (<b>A</b>) and rainfall rate (<b>B</b>) series scatterplot for Genoa University rain gauge (1833–2019); MM30: mean mobile at 30 years; S: Sen’s trend line; Sl: Sen’s lower confidence (95%) slope; Su: Sen’s upper confidence (95%) slope.</p>
Full article ">Figure 12
<p>Rainy days (<b>A</b>) and rainfall rate (<b>B</b>) series scatterplot for Chiavari rain gauge (1877–2019); MM30: mean mobile at 30 years; S: Sen’s trend line; Sl: Sen’s lower confidence (95%) slope; Su: Sen’s upper confidence (95%) slope.</p>
Full article ">Figure 13
<p>Maximum rainfall in 1 (<b>A</b>), 3 (<b>B</b>), 6 (<b>C</b>), 12 (<b>D</b>), and 24 (<b>E</b>) h in the yearly series of Chiavari (1932–2019) and Genova Ponte Carrega (1939–2019).</p>
Full article ">Figure 14
<p>Comparison between the rainfall threshold established for different periods, from the past to the present: 1910–1999 (yellow curve), 2000–2019 (green curve), and 1910–2019 (red curve).</p>
Full article ">
16 pages, 2176 KiB  
Article
Pollution Caused by Potentially Toxic Elements Present in Road Dust from Industrial Areas in Korea
by Hyeryeong Jeong, Jin Young Choi, Jaesoo Lim and Kongtae Ra
Atmosphere 2020, 11(12), 1366; https://doi.org/10.3390/atmos11121366 - 17 Dec 2020
Cited by 16 | Viewed by 3509
Abstract
We examined the pollution characteristics of potentially toxic elements (PTEs) in road dust (RD) from nine industrial areas in South Korea to assess PTE pollution levels and their environmental risks for devising better strategies for managing RD. The median concentrations (mg/kg) were in [...] Read more.
We examined the pollution characteristics of potentially toxic elements (PTEs) in road dust (RD) from nine industrial areas in South Korea to assess PTE pollution levels and their environmental risks for devising better strategies for managing RD. The median concentrations (mg/kg) were in the order Zn (1407) > Cr (380) > Cu (276) > Pb (260) > Ni (112) > As (15) > Cd (2) > Hg (0.1). The concentration of PTEs was the highest at the Onsan Industrial Complex, where many smelting facilities are located. Our results show that Onsan, Noksan, Changwon, Ulsan, Pohang, and Shihwa industrial areas are heavily polluted with Cu, Zn, Cd, and Pb. The presence of these toxic elements in RD from the impervious layer in industrial areas may have a moderate to severe effect on the health of the biota present in these areas. The potential ecological risk index (Eri) for PTEs was in the decreasing order of Cd > Pb > Hg > Cu > As > Zn > Ni > Cr, indicating that the dominant PTE causing ecological hazards is Cd owing to its high toxicity. Our research suggests the necessity for the urgent introduction of an efficient management strategy to reduce RD, which adds to coastal pollution and affects human health. Full article
Show Figures

Figure 1

Figure 1
<p>Map of the study area showing road dust sampling locations from 9 different industrial areas of Korea.</p>
Full article ">Figure 2
<p>Comparison of PTEs concentrations in RD from different Industrial areas of Korea. The bar height and error bar represents the median values and standard deviation (SD) of the data (note the log scale on the <span class="html-italic">y</span>–axis).</p>
Full article ">Figure 3
<p>Ecological risk factors (<math display="inline"><semantics> <mrow> <msubsup> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">r</mi> <mi mathvariant="normal">i</mi> </msubsup> </mrow> </semantics></math>) of individual toxic elements in RD samples in this study.</p>
Full article ">Figure 4
<p>Comparison of median potentially toxic elements (PTEs) load (mg/m<sup>2</sup>) in RD samples from different industrial areas in this study. Log-scale was used for arithmetic scaling of the y-axis. Values are expressed as median ± standard deviation (SD) of the measured metal data.</p>
Full article ">Figure 5
<p>Cumulative curves (<b>a</b>) and particle size distributions (<b>b</b>) in the RD samples from different national industrial areas of this study.</p>
Full article ">
17 pages, 2392 KiB  
Article
Spatial Distribution, Source Apportionment, Ozone Formation Potential, and Health Risks of Volatile Organic Compounds over a Typical Central Plain City in China
by Kun He, Zhenxing Shen, Jian Sun, Yali Lei, Yue Zhang and Xin Wang
Atmosphere 2020, 11(12), 1365; https://doi.org/10.3390/atmos11121365 - 16 Dec 2020
Cited by 9 | Viewed by 3125
Abstract
The profiles, contributions to ozone formation, and associated health risks of 56 volatile organic compounds (VOCs) species were investigated using high time resolution observations from photochemical assessment monitoring stations (PAMs) in Luoyang, China. The daily averaged concentration of total VOCs (TVOCs) was 21.66 [...] Read more.
The profiles, contributions to ozone formation, and associated health risks of 56 volatile organic compounds (VOCs) species were investigated using high time resolution observations from photochemical assessment monitoring stations (PAMs) in Luoyang, China. The daily averaged concentration of total VOCs (TVOCs) was 21.66 ± 10.34 ppbv in urban areas, 14.45 ± 7.40 ppbv in suburbs, and 37.58 ± 13.99 ppbv in an industrial zone. Overall, the VOCs levels in these nine sites followed a decreasing sequence of alkanes > aromatics > alkenes > alkyne. Diurnal variations in VOCs exhibited two peaks at 8:00–9:00 and 19:00–20:00, with one valley at 23:00–24:00. Source apportionment indicated that vehicle and industrial emissions were the dominant sources of VOCs in urban and suburban sites. The industrial site displayed extreme levels, with contributions from petrochemical-related sources of up to 38.3%. Alkenes and aromatics displayed the highest ozone formation potentials because of their high photochemical reactivity. Cancer and noncancer risks in the industrial site were higher than those in the urban and suburban areas, and USEPA possible risk thresholds were reached in the industrial site, indicating PAMs VOC–related health problems cannot be ignored. Therefore, vehicle and industrial emissions should be prioritized when considering VOCs and O3 control strategies in Luoyang. Full article
Show Figures

Figure 1

Figure 1
<p>Locations of nine sampling sites in Luoyang, China.</p>
Full article ">Figure 2
<p>Daily variations of TVOCs concentrations and the main compositions in urban, suburban and industrial areas (1,2,3,4 in x-axis denotes morning period 8:00–9:00, afternoon period 15:00–16:00, evening period 19:00–20:00, midnight period 23:00–24:00, respectively).</p>
Full article ">Figure 3
<p>Source contributions to VOCs determined by PMF in (<b>a</b>) urban, (<b>b</b>) suburban and (<b>c</b>) industrial sites.</p>
Full article ">Figure 4
<p>The ratios of concentration and OFP for four types of VOCs to TVOCs in urban, suburban and industrial areas.</p>
Full article ">
19 pages, 1983 KiB  
Article
Impact of Environmental Conditions on Grass Phenology in the Regional Climate Model COSMO-CLM
by Eva Hartmann, Jan-Peter Schulz, Ruben Seibert, Marius Schmidt, Mingyue Zhang, Jürg Luterbacher and Merja H. Tölle
Atmosphere 2020, 11(12), 1364; https://doi.org/10.3390/atmos11121364 - 16 Dec 2020
Cited by 3 | Viewed by 3359
Abstract
Feedbacks of plant phenology to the regional climate system affect fluxes of energy, water, CO2, biogenic volatile organic compounds as well as canopy conductance, surface roughness length, and are influencing the seasonality of albedo. We performed simulations with the regional climate model COSMO-CLM [...] Read more.
Feedbacks of plant phenology to the regional climate system affect fluxes of energy, water, CO2, biogenic volatile organic compounds as well as canopy conductance, surface roughness length, and are influencing the seasonality of albedo. We performed simulations with the regional climate model COSMO-CLM (CCLM) at three locations in Germany covering the period 1999 to 2015 in order to study the sensitivity of grass phenology to different environmental conditions by implementing a new phenology module. We provide new evidence that the annually-recurring standard phenology of CCLM is improved by the new calculation of leaf area index (LAI) dependent upon surface temperature, day length, and water availability. Results with the new phenology implemented in the model show a significantly higher correlation with observations than simulations with the standard phenology. The interannual variability of LAI improves the representation of vegetation in years with extremely warm winter/spring (e.g., 2007) or extremely dry summer (e.g., 2003) and shows a more realistic growth period. The effect of the newly implemented phenology on atmospheric variables is small but tends to be positive. It should be used in future applications with an extension on more plant functional types. Full article
Show Figures

Figure 1

Figure 1
<p>Map (<b>center</b>) with the three experimental locations (Lindenberg, Linden, and Selhausen) surrounded by their (Selhausen—<b>topleft</b>, Linden—<b>bottomleft</b>, Lindenberg—<b>right</b>) climate diagrams (data from Reference [<a href="#B42-atmosphere-11-01364" class="html-bibr">42</a>], 1982–2012), where the monthly precipitation sums are marked with blue bars and the annual cycle of monthly mean temperatures is shown in red.</p>
Full article ">Figure 2
<p>Leaf area index (LAI) satellite (dotted) and in-situ (solid lines) observations at Linden and Selhausen for the years shown in the legend on the left in different colors. In-situ measurements are only available for the given years and dates. At Linden, the shown simulated years are (except 1998) the same as the in-situ observations. At Selhausen, the six years of simulations before the in-situ observations are shown. The mean annual cycle of the satellite LAI for the given years is shown in black (<span style="color: #000000">–</span>).</p>
Full article ">Figure 3
<p>Mean (1999–2015) annual cycle of LAI. Results with the standard phenology(_old, <span style="color:purple">–</span>), with only the dependence on temperature implemented (_T, <span style="color: #0000FF">–</span>), with the dependence on day length added (_TD, <span style="color: #00FF00">–</span>), with the fully implemented new phenology (_TDW, <span style="color: #FF0000">–</span>), and satellite observations (_Obs, <span style="color: #000000">–</span>) are shown at the three experimental domains Lindenberg, Linden, and Selhausen.</p>
Full article ">Figure 4
<p>Start of the growing season (SGS) in number of days for each year from 1999 to 2015 and each domain (Lindenberg, Linden, and Selhausen) for satellite observations (_Obs, <span style="color: #000000">–</span>), the standard phenology simulations (_old, <span style="color:purple">–</span>), and the new phenology simulations (_TDW, <span style="color: #FF0000">–</span>).</p>
Full article ">Figure 5
<p>Annual cycle of LAI of the year 2003 with extremely dry summer (<b>top</b>) and the year 2007 with extremely warm spring (<b>bottom</b>) at Lindenberg (<b>left</b>), Linden (<b>middle</b>), and Selhausen (<b>right</b>). In black (<span style="color: #000000">–</span>) are the satellite observations (_Obs) and in different colors the simulations with the standard phenology (_old, <span style="color:purple">–</span>), with only the dependence on temperature (_T, <span style="color: #0000FF">–</span>), the dependence on temperature and day length (_TD, <span style="color: #00FF00">–</span>), and with the new phenology (_TDW, <span style="color: #FF0000">–</span>).</p>
Full article ">Figure 6
<p>Mean diurnal cycle of latent heat flux during March 2007 and during July 2003 at Lindenberg for the simulations with the standard phenology (_old, <span style="color:purple">–</span>), the simulations with the new phenology (_TDW, <span style="color: #FF0000">–</span>) and the observations (_Obs, <span style="color: #000000">–</span>).</p>
Full article ">Figure 7
<p>Heavy precipitation events with more than 20 mm per day (<b>top</b>) and very warm days within the 90 th Percentile of the observed maximum temperatures (<b>bottom</b>) in each year of the period 1999 to 2015 at Lindenberg, Linden and Selhausen for the HYRAS observations (_Obs, <span style="color: #000000">–</span>), the simulations with the standard phenology (_old, <span style="color:purple">–</span>) and the simulations with the new phenology (_TDW, <span style="color: #FF0000">–</span>).</p>
Full article ">
27 pages, 8504 KiB  
Article
A Long-Term, 1-km Resolution Daily Meteorological Dataset for Modeling and Mapping Permafrost in Canada
by Yu Zhang, Budong Qian and Gang Hong
Atmosphere 2020, 11(12), 1363; https://doi.org/10.3390/atmos11121363 - 16 Dec 2020
Cited by 5 | Viewed by 4568
Abstract
Climate warming is causing permafrost thaw and there is an urgent need to understand the spatial distribution of permafrost and its potential changes with climate. This study developed a long-term (1901–2100), 1-km resolution daily meteorological dataset (Met1km) for modeling and mapping permafrost at [...] Read more.
Climate warming is causing permafrost thaw and there is an urgent need to understand the spatial distribution of permafrost and its potential changes with climate. This study developed a long-term (1901–2100), 1-km resolution daily meteorological dataset (Met1km) for modeling and mapping permafrost at high spatial resolutions in Canada. Met1km includes eight climate variables (daily minimum, maximum, and mean air temperatures, precipitation, vapor pressure, wind speed, solar radiation, and downward longwave radiation) and is suitable to drive process-based permafrost and other land-surface models. Met1km was developed based on four coarser gridded meteorological datasets for the historical period. Future values were developed using the output of a new Canadian regional climate model under medium-low and high emission scenarios. These datasets were downscaled to 1-km resolution using the re-baselining method based on the WorldClim2 dataset as spatial templates. We assessed Met1km by comparing it to climate station observations across Canada and a gridded monthly anomaly time-series dataset. The accuracy of Met1km is similar to or better than the four coarser gridded datasets. The errors in long-term averages and average seasonal patterns are small. The error occurs mainly in day-to-day fluctuations, thus the error decreases significantly when averaged over 5 to 10 days. Met1km, as a data generating system, is relatively small in data volume, flexible to use, and easy to update when new or improved source datasets are available. The method can also be used to generate similar datasets for other regions, even for the entire global landmass. Full article
(This article belongs to the Special Issue Climate Data for Agricultural Applications: Downscaling and Scenarios)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Air temperature and (<b>b</b>) precipitation generated by Met1km for a grid at Yellowknife, Northwest Territories (62.4540° N, 114.3718° W). Black curves are for the historical period (1901–2017), blue and red curves are for future scenarios under RCP 4.5 and 8.5, respectively. The dash curves are annual values, and the bold curves are 10-year moving averages.</p>
Full article ">Figure 2
<p>Spatial distributions of mean air temperature and changes with time generated by Met1km. (<b>a</b>) Mean air temperature in the 2000s (averaged from 2001 to 2010), (<b>b</b>) changes from the 1900s (averaged from 1901 to 1910) to the 2000s, and changes from the 2000s to the 2090s (averaged from 2091 to 2100) under scenarios of (<b>c</b>) RCP 4.5 and (<b>d</b>) RCP 8.5.</p>
Full article ">Figure 3
<p>Correlation coefficients (R) between two climate stations for (<b>a</b>) daily mean air temperature and (<b>b</b>) daily precipitation and mean absolute differences between two climate stations for (<b>c</b>) daily mean air temperature and (<b>d</b>) daily precipitation. They were calculated using the deviations (differences for air temperature and ratios for precipitation) from the long-term seasonal patterns, which were linearly interpolated from the long-term monthly averages.</p>
Full article ">Figure 4
<p>Changes of correlation coefficients (R) (<b>a</b>,<b>b</b>) and mean absolute differences (<b>c</b>,<b>d</b>) between climate stations with duration of running average for air temperature (<b>a</b>,<b>c</b>) and precipitation (<b>b</b>,<b>d</b>). They were calculated using the deviations (differences for air temperature and ratios for precipitation) from the long-term daily averages, which were interpolated from the long-term monthly averages. Each curve represents the median values calculated from all the stations within a certain distance range.</p>
Full article ">Figure 5
<p>Comparisons between Met1km dataset and observations at the Yellowknife airport climate station in 1953 for difference climate variable. (<b>a</b>) Daily minimum air temperature; (<b>b</b>) daily maximum air temperature; (<b>c</b>) precipitation; (<b>d</b>) vapor pressure; (<b>e</b>) solar radiation; and (<b>f</b>) wind speed.</p>
Full article ">Figure 6
<p>(<b>a</b>) Correlation coefficients (R) between Met1km and measurements at the Yellowknife airport climate station, and (<b>b</b>) the total error and the three error components of Met1km calculated by comparing to the observations at the climate station. The daily observations for precipitation are from 1948 to 2016, and other observations are from 1953 to 2005.</p>
Full article ">Figure 7
<p>Changes of mean absolute error (MAE) with the duration of running average (<b>a</b>) for the original data and (<b>b</b>) for the error in deviations from the long-term average seasonal patterns. The daily observations for precipitation are from 1948 to 2016, and other observations are from 1953 to 2005.</p>
Full article ">Figure 8
<p>Comparisons of the probability distributions of Met1km data with that of the observations at the Yellowknife airport climate station. The x-axis is expressed as (<b>a</b>–<b>f</b>) deviations from or (<b>g</b>) ratio to the long-term average seasonal pattern, which is interpolated from long-term monthly averages. The inset in panel (<b>g</b>) shows the same data but in different scales. The daily data are from 1948 to 2016 for precipitation and from 1953 to 2005 for other climate variables.</p>
Full article ">Figure 9
<p>The errors of Met1km comparing with the homogenized daily air temperature and precipitation station data for the period from 1945 to 2014 (1948 to 2014 in eastern Canada) and from 1901 to 1944 (1901 to 1947 for eastern Canada). Panels (<b>a</b>–<b>c</b>) are for the MAE, and panels (<b>b</b>,<b>d</b>) are for the three error components (error in long-term averages (absolute of mean error), MAE in average seasonal patterns, and MAE for the deviations from long-term average seasonal patterns). The large bars and the error bars are for the median and the standard deviations of the errors of all the available climate stations.</p>
Full article ">Figure 10
<p>Correlation coefficients between the monthly anomalies of Met1km and the CANGRD dataset. The daily data from Met1km were converted to monthly anomalies corresponding to the grids of the CANGRD dataset. The large bars and the error bars are the averages and standard deviations of the 4890 grids of the CANGRD dataset.</p>
Full article ">Figure 11
<p>Correlation coefficients between the CANGRD and the anomalies of Met1km datasets for (<b>a</b>) annual mean air temperature and (<b>b</b>) annual total precipitation. The daily data from Met1km were converted to anomalies corresponding to the CANGRD dataset. The data were from 1901 in southern Canada (south of 60° N) and from 1948 in northern Canada.</p>
Full article ">Figure 12
<p>Comparisons of the errors of the source datasets, the spatially downscaled source datasets and Met1km (the three bars in each group respectively) (<b>a</b>–<b>p</b>). The errors were calculated by comparing with the homogenized daily climate station data [<a href="#B43-atmosphere-11-01363" class="html-bibr">43</a>]. The four rows are for different climate variables, and the four columns are for the errors in long-term averages and the mean absolute errors for daily values, 10-day averages and 30-day averages. The large bars and the error bars are the averages and the standard deviations of the errors of all the climate stations. The errors of Met1km (the red hatched bars) are different in each panel due to differences in periods of the datasets and spatial domains of the NRCANmet and PNWNAmet datasets.</p>
Full article ">
19 pages, 591 KiB  
Article
Sensitivity of Surface Fluxes in the ECMWF Land Surface Model to the Remotely Sensed Leaf Area Index and Root Distribution: Evaluation with Tower Flux Data
by David Stevens, Pedro M. A. Miranda, René Orth, Souhail Boussetta, Gianpaolo Balsamo and Emanuel Dutra
Atmosphere 2020, 11(12), 1362; https://doi.org/10.3390/atmos11121362 - 16 Dec 2020
Cited by 11 | Viewed by 3983
Abstract
The surface-atmosphere turbulent exchanges couple the water, energy and carbon budgets in the Earth system. The biosphere plays an important role in the evaporation process, and vegetation related parameters such as the leaf area index (LAI), vertical root distribution and stomatal resistance are [...] Read more.
The surface-atmosphere turbulent exchanges couple the water, energy and carbon budgets in the Earth system. The biosphere plays an important role in the evaporation process, and vegetation related parameters such as the leaf area index (LAI), vertical root distribution and stomatal resistance are poorly constrained due to sparse observations at the spatio-temporal scales at which land surface models (LSMs) operate. In this study, we use the Carbon Hydrology Tiled European Center for Medium-Range Weather Forecasts (ECMWF) Scheme for Surface Exchanges over Land (CHTESSEL) model and investigate the sensitivity of the simulated turbulent fluxes to these vegetation related parameters. Observed data from 17 FLUXNET towers were used to force and evaluate model simulations with different vegetation parameter configurations. The replacement of the current LAI climatology used by CHTESSEL, by a new high-resolution climatology, representative of the station’s location, has a small impact on the simulated fluxes. Instead, a revision of the root profile considering a uniform root distribution reduces the underestimation of evaporation during water stress conditions. Despite the limitations of using only one model and a limited number of stations, our results highlight the relevance of root distribution in controlling soil moisture stress, which is likely to be applicable to other LSMs. Full article
Show Figures

Figure 1

Figure 1
<p>Latent (<b>a</b>–<b>d</b>, top panels) and sensible (<b>e</b>–<b>h</b>, bottom panels) heat flux evaluation in terms of mean bias error (MBE), SD, normalized mean error (NME) and r for CTR, LAI, LAI_RSMIN and LAI_ROOT. The boxplots represent each metric’s distribution for the 17 towers showing the percentiles of 25, 50 and 75. Symbols denote outliers for values greater than 1.5 times the interquartile range from the nearest 25th or 75th percentile.</p>
Full article ">Figure 2
<p>Latent heat flux (Qle) and soil moisture in Amplero (<b>a</b>), Blodgett (<b>b</b>) and Espirra (<b>c</b>). For the Qle (left axis), the observations are in grey and the simulations from MLAI in blue and from MLAI_NOSM in orange. The soil moisture (right axis) for MLAI simulation is shown for the top 3 layers (red solid) and deep layers (red dashed). The light dotted red and blue lines represent soil moisture at field capacity and the wilting point, respectively.</p>
Full article ">Figure 3
<p>Simulated and observed latent heat flux during one example year with a 14 day running mean smoothing. Observations are in grey, CTR in black, MLAI in blue, MLAI_ROOT in green and MLAI_RSMIN in red.</p>
Full article ">
2 pages, 128 KiB  
Editorial
Special Issue Air Quality and Smoke Management
by Scott L. Goodrick
Atmosphere 2020, 11(12), 1361; https://doi.org/10.3390/atmos11121361 - 15 Dec 2020
Viewed by 1749
Abstract
The Atmosphere Special Issue “Special Issue Air Quality and Smoke Management” explores our ability to simulate wildland fire smoke events and the potential to link such modeling to future studies of human health impacts [...] Full article
(This article belongs to the Special Issue Air Quality and Smoke Management)
14 pages, 2106 KiB  
Article
Evaluating Satellite Sounding Temperature Observations for Cold Air Aloft Detection
by Rebekah Esmaili, Nadia Smith, Mark Schoeberl and Chris Barnet
Atmosphere 2020, 11(12), 1360; https://doi.org/10.3390/atmos11121360 - 15 Dec 2020
Cited by 3 | Viewed by 2334
Abstract
Cold Air Aloft (CAA) can impact commercial flights when cold air descends below 12,192 m (40,000 ft) and temperatures drop dramatically. A CAA event is identified when air temperature falls below −65 °C, which decreases fuel efficiency and poses a safety hazard. This [...] Read more.
Cold Air Aloft (CAA) can impact commercial flights when cold air descends below 12,192 m (40,000 ft) and temperatures drop dramatically. A CAA event is identified when air temperature falls below −65 °C, which decreases fuel efficiency and poses a safety hazard. This manuscript assesses the performance of the National Oceanic and Atmospheric Administration Unique Combined Atmospheric Processing System (NUCAPS) in detecting CAA events using sounders on polar-orbiting satellites. We compare NUCAPS air temperature profiles with those from Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) for January–March 2018. Of 1311 collocated profiles, 236 detected CAA. Our results showed that NUCAPS correctly detects CAA in 48.1% of profiles, while 17.2% are false positives and 34.7% are false negatives. To identify the reason for these detection states, we used a logistic regression trained on NUCAPS diagnostic parameters. We found that cloud cover can impact the skill even at higher vertical levels. This work indicates that a CAA-specific quality flag is feasible and may be useful to help forecasters to diagnose NUCAPS in real-time. Furthermore, the inclusion of an additional sounder data source (e.g., NOAA-20) may increase CAA forecast accuracy. Cloud scenes change rapidly, so additional observations provide more opportunities for correct detection. Full article
(This article belongs to the Special Issue Weather and Aviation Safety)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Illustration of how the tropopause height is lower in the poles. (<b>b</b>) Flights are often diverted when widespread Cold Air Aloft (CAA) is present in the upper troposphere and lower stratosphere, which coincides with commercial flight zones in high latitudes. Air temperatures below −65 °C can cause water within the jet fuel to freeze and common fuels begins to form wax crystals, which can reduce fuel efficiency or pose a safety hazard.</p>
Full article ">Figure 2
<p>Radiosonde temperature profiles from 14 Alaskan sites for February through March 2018. A total of 1180 profiles are shown. Radiosondes are released twice daily at 00 UTC and 12 UTC. Red points indicate where CAA is detected. The vertical layer between the two dashed lines in each plot shows the typical transpolar flight cruising altitudes, assuming a standard atmosphere.</p>
Full article ">Figure 3
<p>A National Oceanic and Atmospheric Administration Unique Combined Atmospheric Processing System (NUCAPS) horizontal slice along the 200 hPa isobar on 26 February 2018 over the study domain, 40–80° N, 169° E–131° W. (<b>a</b>) Several combined overpasses between 1100 and 1700 UTC and (<b>b</b>) matchups between with Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) profiles within 150 km for the same time period. Blue coloring indicates where the air temperature is below 208 K (−65 °C) to distinguish the presence of a CAA event. Profiles annotated in (<b>b</b>) are shown in <a href="#atmosphere-11-01360-f004" class="html-fig">Figure 4</a>.</p>
Full article ">Figure 4
<p>(<b>a</b>) NUCAPS profiles showing the matchups between COSMIC (black dashed line) and NUCAPS (blue solid lines), where the shaded blue box indicates the temperature and pressure criteria for CAA. The profile locations are annotated in <a href="#atmosphere-11-01360-f003" class="html-fig">Figure 3</a>. (<b>b</b>) The mean (solid line) and standard deviation (error bars) of the averaging kernel matrix (AKM) diagonal vectors are a function of pressure in the study domain on 26 February 2018.</p>
Full article ">Figure 5
<p>Percentage of matchup profiles between January and March 2018 for the true positive, false positive, and false negative CAA detection states that are summarized in <a href="#atmosphere-11-01360-t001" class="html-table">Table 1</a>. Profiles categorized as true negative are the most numerous (N = 1026) but for clarity, are not shown.</p>
Full article ">Figure 6
<p>Results for a logistic regression on diagnostic parameters for the four CAA detection states described in <a href="#atmosphere-11-01360-t001" class="html-table">Table 1</a>. X indicates diagnostic parameters that do not meet statistical significance requirements.</p>
Full article ">
2 pages, 147 KiB  
Editorial
Special Issue Editorial: Study of Brake Wear Particle Emissions
by Jens Wahlström
Atmosphere 2020, 11(12), 1359; https://doi.org/10.3390/atmos11121359 - 15 Dec 2020
Cited by 2 | Viewed by 1808
(This article belongs to the Special Issue Study of Brake Wear Particle Emissions)
31 pages, 2404 KiB  
Review
Radiation in the Atmosphere—A Hazard to Aviation Safety?
by Matthias M. Meier, Kyle Copeland, Klara E. J. Klöble, Daniel Matthiä, Mona C. Plettenberg, Kai Schennetten, Michael Wirtz and Christine E. Hellweg
Atmosphere 2020, 11(12), 1358; https://doi.org/10.3390/atmos11121358 - 14 Dec 2020
Cited by 26 | Viewed by 9452
Abstract
Exposure of aircrew to cosmic radiation has been recognized as an occupational health risk for several decades. Based on the recommendations by the International Commission on Radiological Protection (ICRP), many countries and their aviation authorities, respectively have either stipulated legal radiation protection regulations, [...] Read more.
Exposure of aircrew to cosmic radiation has been recognized as an occupational health risk for several decades. Based on the recommendations by the International Commission on Radiological Protection (ICRP), many countries and their aviation authorities, respectively have either stipulated legal radiation protection regulations, e.g., in the European Union or issued corresponding advisory circulars, e.g., in the United States of America. Additional sources of ionizing and non-ionizing radiation, e.g., due to weather phenomena have been identified and discussed in the scientific literature in recent years. This article gives an overview of the different generally recognized sources due to weather as well as space weather phenomena that contribute to radiation exposure in the atmosphere and the associated radiation effects that might pose a risk to aviation safety at large, including effects on human health and avionics. Furthermore, potential mitigation measures for several radiation sources and the prerequisites for their use are discussed. Full article
(This article belongs to the Special Issue Weather and Aviation Safety)
Show Figures

Figure 1

Figure 1
<p>The radiation field at aviation altitudes is created in repeated interactions of the primary and secondary radiation with the constituents of the atmosphere.</p>
Full article ">Figure 2
<p>Variations of the effective dose rate over time (<b>a</b>) and the dependency on the magnetic shielding expressed by the cut-off rigidity <span class="html-italic">R<sub>C</sub></span> (<b>b</b>) at 29,000 ft. and 41,000 ft. during solar minimum (dose rates are calculated with the PANDOCA model [<a href="#B29-atmosphere-11-01358" class="html-bibr">29</a>]).</p>
Full article ">Figure 3
<p>Model predicted GCR total neutron fluence rate dependence on the altitude from PANDOCA [<a href="#B22-atmosphere-11-01358" class="html-bibr">22</a>] (black solid line) and CARI-7 [<a href="#B25-atmosphere-11-01358" class="html-bibr">25</a>] (black dashed line) and for neutron energies greater than 10 MeV from PANDOCA (grey line) without geomagnetic shielding (<span class="html-italic">R<sub>C</sub></span> = 0 GV) during solar minimum.</p>
Full article ">Figure 4
<p>Dose rates and accumulated doses were calculated for the altitude of 41,000 ft. (FL410) and locations at 90° N 0° for GLE42 and 70° N 50° E for GLE70, respectively.</p>
Full article ">Figure 5
<p>D-Scale for aviation with Space Weather D-Indices from 0 to 5.</p>
Full article ">
21 pages, 2674 KiB  
Article
Air Pollution Measurements and Land-Use Regression in Urban Sub-Saharan Africa Using Low-Cost Sensors—Possibilities and Pitfalls
by Asmamaw Abera, Kristoffer Mattisson, Axel Eriksson, Erik Ahlberg, Geremew Sahilu, Bezatu Mengistie, Abebe Genetu Bayih, Abraham Aseffaa, Ebba Malmqvist and Christina Isaxon
Atmosphere 2020, 11(12), 1357; https://doi.org/10.3390/atmos11121357 - 14 Dec 2020
Cited by 18 | Viewed by 5020
Abstract
Air pollution is recognized as the most important environmental factor that adversely affects human and societal wellbeing. Due to rapid urbanization, air pollution levels are increasing in the Sub-Saharan region, but there is a shortage of air pollution monitoring. Hence, exposure data to [...] Read more.
Air pollution is recognized as the most important environmental factor that adversely affects human and societal wellbeing. Due to rapid urbanization, air pollution levels are increasing in the Sub-Saharan region, but there is a shortage of air pollution monitoring. Hence, exposure data to use as a base for exposure modelling and health effect assessments is also lacking. In this study, low-cost sensors were used to assess PM2.5 (particulate matter) levels in the city of Adama, Ethiopia. The measurements were conducted during two separate 1-week periods. The measurements were used to develop a land-use regression (LUR) model. The developed LUR model explained 33.4% of the variance in the concentrations of PM2.5. Two predictor variables were included in the final model, of which both were related to emissions from traffic sources. Some concern regarding influential observations remained in the final model. Long-term PM2.5 and wind direction data were obtained from the city’s meteorological station, which should be used to validate the representativeness of our sensor measurements. The PM2.5 long-term data were however not reliable. Means of obtaining good reference data combined with longer sensor measurements would be a good way forward to develop a stronger LUR model which, together with improved knowledge, can be applied towards improving the quality of health. A health impact assessment, based on the mean level of PM2.5 (23 µg/m3), presented the attributable burden of disease and showed the importance of addressing causes of these high ambient levels in the area. Full article
(This article belongs to the Special Issue Challenges in Measuring and Assessing Environmental Health)
Show Figures

Figure 1

Figure 1
<p>The location of Adama on the East African continent and the selected classes of variables for the land-use regression (LUR) model.</p>
Full article ">Figure 2
<p>Time series of particulate matter (PM2.5) were obtained by changing the positions of 5 sensors such that 20 sites were sampled during the campaign (see <a href="#sec2dot2-atmosphere-11-01357" class="html-sec">Section 2.2</a> and <a href="#sec4dot2-atmosphere-11-01357" class="html-sec">Section 4.2</a> for details). For comparison see <a href="#app1-atmosphere-11-01357" class="html-app">Figure S3</a> in which sensors were not moved.</p>
Full article ">Figure 3
<p>Original classes for road type in open street maps were reclassified to seven classes.</p>
Full article ">Figure 4
<p>Road distribution for all road types included in Adama and measurement sites locations.</p>
Full article ">Figure 5
<p>A 4-day time series measured by Purple Air illustrating the typical diurnal pattern, with elevated levels during early morning and late afternoon rush hours, seen throughout the Adama measurement sites.</p>
Full article ">Figure 6
<p>Frequency of wind direction for 2018–2019 during (<b>a</b>) the dry season, and (<b>b</b>) the wet season.</p>
Full article ">
12 pages, 3756 KiB  
Article
The Lagged Effect of Anthropogenic Aerosol on East Asian Precipitation during the Summer Monsoon Season
by Su-Jung Lee, Yong-Cheol Jeong and Sang-Wook Yeh
Atmosphere 2020, 11(12), 1356; https://doi.org/10.3390/atmos11121356 - 14 Dec 2020
Cited by 2 | Viewed by 2899
Abstract
The authors investigated the lagged effect of anthropogenic aerosols (AAs) during the premonsoon season (April–May–June) on the East Asian precipitation during the postmonsoon season (July–August) using the aerosol optical depth (AOD) from a satellite dataset and reanalysis datasets. When the AOD is high [...] Read more.
The authors investigated the lagged effect of anthropogenic aerosols (AAs) during the premonsoon season (April–May–June) on the East Asian precipitation during the postmonsoon season (July–August) using the aerosol optical depth (AOD) from a satellite dataset and reanalysis datasets. When the AOD is high in Eastern China during the premonsoon season, the amount of precipitation increases in the western North Pacific, including the Korean Peninsula and Japan, during the postmonsoon season. The amount of cloud in the western-to-central North Pacific in the premonsoon season increases during the high-AOD period. Subsequently, it cools the sea surface temperature until the postmonsoon season, which strengthens the North Pacific High. The strengthened North Pacific High in the postmonsoon season expands to the western North Pacific, which leads to the enhancement of the moisture flows from the ocean. This results in the increase in precipitation in the western North Pacific, including the Korean Peninsula and Japan, during the postmonsoon season. Full article
(This article belongs to the Special Issue Aerosol-Climate Interaction)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Distribution of climatological (2000–2017) aerosol optical depth (AOD). The green box indicates the East Asian region (110.5–122.5° E, 22.5–40.5° N), where AOD is the highest. (<b>b</b>) The black line denotes time series data of the AOD index, averaged over the green box area shown in (<b>a</b>). The blue line denotes normalized time series data with the linear trend removed from the black line. (<b>c</b>) The composite difference of AOD during April–May–June (AMJ; high AOD minus low AOD). Regions with black cross-hatching indicate a statistical significance at a 95% confidence level.</p>
Full article ">Figure 2
<p>The composite difference between high AOD and low AOD of (<b>a</b>) AMJ cloud fraction (CF), (<b>b</b>) AMJ surface net solar radiation (SSRD; 10<sup>−7</sup> J/m<sup>2</sup>), and (<b>c</b>) AMJ surface air temperature (SAT; ℃). Regions with black cross-hatching indicate a statistical significance at a 95% confidence level. Solid black lines denote the climatological (2000–2017) mean.</p>
Full article ">Figure 3
<p>The difference (high AOD minus low AOD) of composited precipitation from the Global Precipitation Climatology Project (GPCP) between the high AOD and the low AOD in July–August (JA). Regions with a red dot indicate a statistical significance at a 95% confidence level. The unit is mm/day.</p>
Full article ">Figure 4
<p>The composite difference between high AOD and low AOD of (<b>a</b>) JA surface air temperature (SAT; ℃), (<b>b</b>) JA sea level pressure (SLP; hPa), (<b>c</b>) JA wind speed (shading) (m/s) and moisture transport (vector) (kg*m/s), and (<b>d</b>) JA latent heat flux (10<sup>−6</sup> J/m<sup>2</sup>). Regions with black cross-hatching indicate a statistical significance at a 95% confidence level. Solid black lines denote the climatological (2000–2017) mean.</p>
Full article ">Figure 5
<p>The composite difference of AOD during JA (high AOD minus low AOD). Regions with black cross-hatching indicate a statistical significance at a 95% confidence level.</p>
Full article ">Figure 6
<p>Schematic diagram showing the process through which AA variability in AMJ affects precipitation in JA.</p>
Full article ">
12 pages, 10892 KiB  
Article
Air Pollution and Long Term Mental Health
by Younoh Kim, James Manley and Vlad Radoias
Atmosphere 2020, 11(12), 1355; https://doi.org/10.3390/atmos11121355 - 14 Dec 2020
Cited by 28 | Viewed by 5136
Abstract
We study the long-term consequences of air pollution on mental health, using a natural experiment in Indonesia. We find that exposure to severe air pollution has significant and persistent consequences on mental health. An extra standard deviation in the pollution index raises the [...] Read more.
We study the long-term consequences of air pollution on mental health, using a natural experiment in Indonesia. We find that exposure to severe air pollution has significant and persistent consequences on mental health. An extra standard deviation in the pollution index raises the probability of clinical depression measured 10 years past exposure by almost 1%. Women in particular seem to be more affected, but some effects persist for men as well. Pollution exposure increases the likelihood of clinical depression for women and also the severity of depressive symptoms for both sexes. It is not clear if men are more resistant to pollution or they simply recover faster from it. Education, perceived economic status, and marriage seem to be the best mitigators for these negative effects. Full article
(This article belongs to the Special Issue Contributions of Aerosol Sources to Health Impacts)
24 pages, 6405 KiB  
Article
Effect of Various Types of ENSO Events on Moisture Conditions in the Humid and Subhumid Tropics
by Daria Gushchina, Irina Zheleznova, Alexander Osipov and Alexander Olchev
Atmosphere 2020, 11(12), 1354; https://doi.org/10.3390/atmos11121354 - 13 Dec 2020
Cited by 17 | Viewed by 3848
Abstract
Moisture anomaly conditions within humid and subhumid tropics that are associated with different types of El Niño and La Niña phenomena are described and analyzed with a focus on their spatial distribution and seasonal variability. Five dryness indices (Keetch–Byram Drought Index, Weighted Anomaly [...] Read more.
Moisture anomaly conditions within humid and subhumid tropics that are associated with different types of El Niño and La Niña phenomena are described and analyzed with a focus on their spatial distribution and seasonal variability. Five dryness indices (Keetch–Byram Drought Index, Weighted Anomaly Standardized Precipitation Index, Standardized Precipitation Index, Palmer Drought Severity Index, and Percent of Normal Precipitation) were derived from ECMWF (European Centre for Medium-Range Weather Forecasts) fifth generation reanalysis (ERA5) reanalysis and University Corporation for Atmospheric Research (UCAR) datasets for the period from 1979 to 2019. Cross-correlation analysis was used to evaluate the relationships between the El Niño Southern Oscillation (ENSO) and selected dryness indices. To describe the seasonal variability of the ENSO–surface moisture relationships, the composite maps of dryness indices in different seasons were analyzed. The results showed a significant heterogeneity of the ENSO-induced moisture anomaly conditions both within and across various geographical regions. Four main areas in humid and subhumid tropics with the maximum effects of El Niño/La Niña events on the surface moisture conditions were found: Southeast Asia and Australia, Eastern and South Africa, Northeastern and Eastern South America, and Central America. It was shown that the effects of La Niña were usually opposite to those of El Niño, while the responses to the two types of El Niño differed mostly in the moisture anomaly intensity and its spatial patterns. Full article
Show Figures

Figure 1

Figure 1
<p>Correlation coefficients between the Eastern Pacific (EP) El Niño and the dryness indices: Keetch–Byram Drought Index (KBDI) (<b>a</b>), Weighted Anomaly Standardized Precipitation Index (WASP) (<b>b</b>), Palmer Drought Severity Index (PDSI) (<b>c</b>), Standardized Precipitation Index using 3-month averages (SPI-3) (<b>d</b>), Percent of Normal Precipitation (PERCENT) (<b>e</b>), and precipitation (<b>f</b>). The hatched areas correspond to a significant correlation at the 95% confidence level.</p>
Full article ">Figure 2
<p>Correlation coefficients between the Central Pacific (CP) El Niño and the dryness indices: KBDI (<b>a</b>), WASP (<b>b</b>), PDSI (<b>c</b>), SPI-3 (<b>d</b>), PERCENT (<b>e</b>), and precipitation (<b>f</b>). The hatched areas correspond to a significant correlation at the 95% confidence level.</p>
Full article ">Figure 3
<p>Correlation coefficients between La Niña and the dryness indices: KBDI (<b>a</b>), WASP (<b>b</b>), PDSI (<b>c</b>), SPI-3 (<b>d</b>), PERCENT (<b>e</b>), and PRECIPITATION (<b>f</b>). The hatched areas correspond to a significant correlation at the 95% confidence level.</p>
Full article ">Figure 4
<p>The Walker circulation cell (averaged between 5° S–5° N) in DJF (December, January, February). The composites for (<b>a</b>) EP El Niños, (<b>b</b>) CP El Niños, and (<b>c</b>) La Niñas (see <a href="#atmosphere-11-01354-t001" class="html-table">Table 1</a> for the selected events). The arrows are the vectors of the wind velocity: divergent zonal wind (m/s) along the horizontal axes, pressure vertical velocity (×10<sup>−2</sup> Pa/s) along the vertical axes; the color field is the value of the pressure vertical velocity (×10<sup>−2</sup> Pa/s).</p>
Full article ">Figure 5
<p>The Hadley circulation cell over (<b>a</b>,<b>d</b>,<b>g</b>) Indonesia and Australia (averaged between 90–160° E); (<b>b</b>,<b>e</b>,<b>h</b>) Africa (averaged between 20° W–50° E); (<b>c</b>,<b>f</b>,<b>i</b>) South America (averaged between 30–80° W) in DJF. The composites are for EP El Niños (left column), CP El Niños (middle column), and La Niñas (right column) (see <a href="#atmosphere-11-01354-t001" class="html-table">Table 1</a> for the selected events). The arrows are the vectors of wind velocity: divergent meridional wind (m/s) along the horizontal, pressure vertical velocity (×10<sup>−2</sup> Pa/s) along the vertical; the color field is the value of the pressure vertical velocity (×10<sup>−2</sup> Pa/s).</p>
Full article ">Figure 6
<p>The composites of dryness indices during DJF (left column) and MAM (March, April, May) (right column) seasons of the EP El Niño years (see <a href="#atmosphere-11-01354-t001" class="html-table">Table 1</a> for selected events).</p>
Full article ">Figure 7
<p>The same as <a href="#atmosphere-11-01354-f006" class="html-fig">Figure 6</a> but for JJA (June, July, August) (left column) and SON (September, October, November) (right column) seasons.</p>
Full article ">Figure 8
<p>The composites of the linear (left column: <b>a</b>,<b>c</b>,<b>e</b>,<b>g</b>,<b>i</b>) and non-linear (right column: <b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>,<b>j</b>) responses of the dryness indices to the ENSO during DJF. The linear response was estimated as El Niño − La Niña; the non-linear response as 0.5 × (El Niño + La Niña).</p>
Full article ">
16 pages, 3167 KiB  
Article
Farming Practices for Reducing Ammonia Emissions in Polish Agriculture
by Arkadiusz Piwowar
Atmosphere 2020, 11(12), 1353; https://doi.org/10.3390/atmos11121353 - 13 Dec 2020
Cited by 11 | Viewed by 3736
Abstract
The main source of ammonia emissions in Poland is agriculture. In 2017, approximately 94% of the total ammonia emissions in Poland came from agriculture, of which the largest part (78%) was related to livestock manure and 22% to nitrogen fertilization. This study presents [...] Read more.
The main source of ammonia emissions in Poland is agriculture. In 2017, approximately 94% of the total ammonia emissions in Poland came from agriculture, of which the largest part (78%) was related to livestock manure and 22% to nitrogen fertilization. This study presents the results of representative research on the implementation of technologies and techniques that reduce ammonia emissions on farms in Poland. The research methodology, including statistical data analysis (multiple correspondence analysis), allowed comparisons to be made of the applied low-carbon practices, taking into account farmers’ characteristics (e.g., age and education) and farm attributes (area size, location, etc.). According to the research, both in the case of mineral fertilization and animal production, farmers in Poland relatively rarely undertake pro-ecological practices aimed at reducing ammonia emissions. The most frequently undertaken activities include dividing the doses of nitrogen fertilizers (in terms of plant production) and the use of feed additives (in terms of livestock production). Empirical studies, supported by correspondence analyses, confirmed a significant differentiation of coexistence and strength of the relationship between the studied variables. The use of correspondence analysis made it possible to precisely recognize the differentiation and co-occurrence of variable categories. In the course of analytical work, a relatively strong correlation was found between the use of divided doses of nitrogen fertilizers and the economic size of farms (φ2 = 0.11571). In turn, the use of feed additives was most strongly determined by the economic size of farms (φ2 = 0.072614) and the location of farms (φ2 = 0.072223). Full article
(This article belongs to the Section Biosphere/Hydrosphere/Land–Atmosphere Interactions)
Show Figures

Figure 1

Figure 1
<p>Main food producing areas in Poland and spatial extent of field studies [<a href="#B37-atmosphere-11-01353" class="html-bibr">37</a>].</p>
Full article ">Figure 2
<p>Declarations of respondents engaged in plant production on the methods of reducing ammonia emissions. Source: own study based on questionnaire surveys (<span class="html-italic">n</span> = 1034)</p>
Full article ">Figure 3
<p>Declarations of respondents engaged in livestock production on the methods of reducing ammonia emissions. Source: own study based on questionnaire surveys (<span class="html-italic">n</span> = 520).</p>
Full article ">Figure 4
<p>Graphical presentation of the results of the correspondence analysis between the declarations of sharing nitrogen fertilizer doses with the researched characteristics of farmers and farms. Source: own study based on questionnaire surveys (<span class="html-italic">n</span> = 1034).</p>
Full article ">Figure 5
<p>Graphical presentation of the results of the correspondence analysis between the declarations of using appropriately selected feed additives with the tested characteristics of farmers and farms. Source: own study based on questionnaire surveys (<span class="html-italic">n</span> = 520).</p>
Full article ">Figure 5 Cont.
<p>Graphical presentation of the results of the correspondence analysis between the declarations of using appropriately selected feed additives with the tested characteristics of farmers and farms. Source: own study based on questionnaire surveys (<span class="html-italic">n</span> = 520).</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop