ALBERT

Description: Due to some harmful effects on humans and the environment, particulate matter (PM) has recently become among the most studied atmospheric pollutants. Given the growing sensitivity to the problem and, since production and accumulation phenomena involving both primary and secondary P M 10 fractions are complex and non-linear, environmental authorities need tools to assess their plans designed to improve the air quality as requested from environmental laws. Multi-criteria decision analysis (MCDA) can be applied to support decision makers, by processing quantitative opinions provided by pools of experts, especially when different views on social aspects should be considered. The results obtained through this approach, however, can be highly dependent on the subjectivity of experts. To partially overcome these challenges, this paper suggests a two-step methodology in which an MCDA is fed with the solution of a multi-objective analysis (MOA). The methodology has been applied to a test case in northern Italy and the results show that this approach is a viable solution for the inclusion of subjective criteria in decision making, while reducing the impact of uncertain expert opinions for data that can be computed through the MOA.

Description: Both the global circulation model (GCM) and regional climate model (RCM) simulations suffer from model biases that eventually result in significant errors in regional forecasts. This model bias issue is addressed using the bias correction approach. This study examines the influence of bias correction on the performance of downscaling simulations of the East Asian winter climate using the Global/Regional Integrated Model system (GRIMs). To assess the bias correction approach, we conducted three sets of simulations for 25 winters (December to February) from 1982 to 2006 over East Asia. The GRIMs were forced by the (1) National Centers for Environmental Prediction (NCEP) Department of Energy (DOE) reanalysis data, (2) original NCEP Climate Forecast System (CFS) data, and (3) bias-corrected CFS data. The GCM climatological means were adjusted based on the NCEP–DOE reanalysis data. The bias correction method was applied to zonal and meridional wind, temperature, geopotential height, specific humidity, and sea surface temperature of the CFS data. The GCM-driven experiments with/without bias correction were compared with the reanalysis-driven simulation. The results of this comparison suggest that the application of bias correction improves the downscaled climate in terms of the climatological mean, inter-annual variability, and extreme events owing to the elimination of errors in large-scale circulations. The effect of bias correction on the simulated extreme event is not as significant as those on the climatological mean and inter-annual variability, but the increased skill appears to be a clue for potential use for predicting extreme events.

Description: Greenhouse gases emitted from soil play a crucial role in the atmospheric environment and global climate change. The theory and technique of detecting stable isotopes in the atmosphere has been widely used to an investigate greenhouse gases from soil. In this paper, we review the current literature on greenhouse gases emitted from soil, including carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O). We attempt to synthesize recent advances in the theory and application of stable isotopes in greenhouse gases from soil and discuss future research needs and directions.

Description: The vertical profiles and trends of temperature and humidity at the South Pole up to 10 km above mean sea level (amsl) were investigated by using radiosonde data collected from March 2005 to February 2018. During an average year between 2005 and 2018, the highest (lowest) temperature in the lower troposphere was approximately −25 °C (−60 °C) in December (July) at a height of about 500 m above the surface (at the surface). A temperature inversion layer above the surface was found during the whole year but was weaker during the summer, while the inversion layers at the tropopause (about 8 km amsl) mostly disappeared during spring and winter. General warming trends were found at all heights and months, but in a few heights and months cooling trends still occurred (e.g., in September below 7 km amsl). Nevertheless, seasonal and yearly averaged temperatures all presented warming trends: 1.1, 1.3, 0.6, 1.5 and 1.1 °C/decade at the surface, and 0.7, 1.0, 0.3, 0.3 and 0.6 °C/decade for the layer average from the surface to 10 km amsl, for spring, summer, autumn, winter, and yearly average, respectively. Most of the water vapor was confined in the lowermost 3 km of the atmosphere with a maximum of 0.35 g kg−1 in December at a 200 m height above surface, and the specific humidity had the similar characteristic of annual cycle and inversion layers as the temperature. At heights below 5 km amsl, increasing trends of specific humidity larger than 0.02 g kg−1/decade occurred during summer months, including the late spring and early autumn, and the annual mean showed an increasing trend of about 0.01–0.02 g kg−1/decade. Meanwhile, above 5 km amsl, the trends became small and generally less than 0.02 g kg−1/decade in all the months, and beyond 7 km amsl the specific humidity remained almost invariant due to its small moisture content as compared with lower levels. From the surface to 10 km amsl, the specific humidity averaged trends of 0.0062, 0.019, 0.0013, 0.002 and 0.007 g kg−1/decade for spring, summer, autumn, winter and yearly average, respectively.

Description: In several studies conducted recently, it was shown that equations pertinent to the electric and magnetic fields produced by electrical charges in motion can be used to calculate the electromagnetic fields produced by current pulses propagating along linearly restricted paths. An example includes the case of current pulses propagating along conductors and conducting channels such as lightning. In this paper, it is shown how the technique can be applied to estimate the electromagnetic fields generated by current and charge distributions moving in arbitrary directions in space. The analysis shows that, depending on the way the problem is formulated using the field equations pertinent to accelerating charges, one procedure leads to the generalized dipole equations, which are independent of the velocity of propagation of the current, and the other procedure leads to a set of equations that depend on the velocity. Using the well-tested transmission line model of lightning return strokes as an example, it is shown that both sets of field equations give rise to the same total electromagnetic field.

Description: Urban-induced thermal stress can threaten human health, especially during heat waves (HWs). The growth of cities further exacerbates this effect. Here, weather research and forecasting (WRF) with an urban canopy model (UCM) is used to assess the effects of megacities and their growth on the thermal regime of proximal cities during heat waves. Analysis of the heat fluxes shows that advection impacts cities downwind. Results indicate that as urban areas change size (50%–100% and 100–150% of their current size), the local 2 m temperature increases by 2.7 and 1.7 °C, and the 2 m specific humidity decreases by 2.1 and 1.4 g kg−1, respectively. A small city downwind is impacted with a 0.3–0.4 °C increase in 2 m temperature. Green roof is a potential mitigation strategy for these regions (i.e., beyond the megacity). With 50% green roofs in an urban area, a 0.5 °C decrease in 2 m temperature and 0.6 g kg−1 increase in specific humidity is simulated. Urbanization upwind of a megacity will contribute to regional climate change.

Description: Runoff in snowy alpine regions is sensitive to climate change in the context of global warming. Exploring the impact of climate change on the runoff in these regions is critical to understand the dynamics of the water cycle and for the improvement of water resources management. In this study, we analyzed the long-term variations in annual runoff in the headwaters region of the Yellow River (HRYR) (a typical snowy mountain region) during the period of 1956–2012. The Soil and Water Assessment Tool (SWAT) with different elevation bands was employed to assess the performance of monthly runoff simulations, and then to evaluate the impacts of climate change on runoff. The results show that the observed runoff for the hydrological stations at lower relative elevations (i.e., Maqu and Tangnaihai stations) had a downward trend, with rates of 1.91 and 1.55 mm/10 years, while a slight upward trend with a rate of 0.26 mm/10 years was observed for the hydrological station at higher elevation (i.e., Huangheyan station). We also found that the inclusion of five elevation bands could lead to more accurate runoff estimates as compared to simulation without elevation bands at monthly time steps. In addition, the dominant cause of the runoff decline across the whole HRYR was precipitation (which explained 64.2% of the decrease), rather than temperature (25.93%).

Description: In order to conduct real-time quantitative monitoring of dust storms, Ka-band millimeter wave radar (MMWR) was utilized for the consecutive detection of dust storms over the Taklimakan Desert from April to June 2018. The retrievals of the reflectivity factor, dust spectrum distribution and dust mass concentration were carried out with the power spectrum data detected by MMWR for three dust storm processes. The analysis shows that: The probability density distribution of dust conforms to the lognormal distribution. During the dust storm processes, the effective detection height of the reflectivity factor was within 2000 m and the range of the reflectivity factors was between-25 dBZ and 25 dBZ. During the floating dust period, the effective height of the dust spectrum distribution was lower than 300 m and the values of dust mass concentration were less than 31.62 μg·m-3,at a height of 200 m. Furthermore, during the blowing sand stage, the effective height of the dust spectrum distribution was normally lower than 600 m and the values of dust mass concentration were mainly less than 316.23 μg·m-3,at a height of 200 m. During the dust storm period, the effective height of the dust spectrum distributionexceeded1000 m; when the height was 100 m, the values of dust mass concentration were between 1220 μg∙m-3 and 42,146 μg∙m-3 and the average mass concentration was 9287 μg·m-3; whereas, the values of dust mass concentration were between 2 μg∙m-3 and 820 μg∙m-3 when the height was 1200m and the average mass concentration was 24 μg∙m-3. The relationship between the reflectivity factor Z and the dust mass concentration M isdefined as Z=651.6M0.796. Compared with the observational data from Grimm180 particle detector, the data of the retrieved dust mass concentration are basically accurate and this retrieved method proves to be feasible. Thus, the MMWR cans be used as a new device for quantitative monitoring of dust storms.

Description: Agricultural emissions are crucial to regional air quality in the autumn and spring due to the intense agricultural activities in Northeast China. However, information on rural ambient particulate matter (PM) in Northeast China is rare, limiting the accurate estimation of agricultural atmospheric particulate matter emissions. In this study, we monitored hourly ambient PM2.5 (PM with a diameter of less than 2.5 μm) concentrations and analyzed daily chemical components (i.e., water-soluble ions, trace elements, organic carbon, and element carbon) at a rural site in Northeast China during the autumn and spring and assessed the impact of agricultural activities on atmospheric PM2.5 concentrations. The results showed that the daily average concentrations of PM2.5 were 143 ± 109 (range: 39–539) μg m−3 from 19 October to 23 November 2017 (i.e., typical harvesting month) and 241 ± 189 (range: 97–976) μg m−3 from 1 April to 13 May 2018 (i.e., typical tilling month). In autumn, the ambient PM2.5 concentrations were high with a Southwest wind, while a Southeast wind caused high PM2.5 concentrations during spring in the rural site. The concentrations of selected water-soluble ions, trace elements, and carbonaceous fractions accounted for 33%, 4%, and 26% of PM2.5 mass concentrations, respectively, in autumn and for 10%, 5%, and 3% of PM2.5 mass concentrations, respectively, in spring. On the basis of the component analysis, straw burning, agricultural machinery, and soil dust driven by wind and tilling were the main contributors to high rural PM2.5 concentrations. In addition, the increasing coal combustion around the rural site was another important source of PM2.5.

Description: Using composite, regular, and partial regression analyses in the six consecutive seasons from spring of the El Niño–Southern Oscillation (ENSO)-/Indian Ocean Dipole (IOD)-developing year through summer following the ENSO/IOD mature phase, the individual and combined impacts of El Niño and positive Indian Ocean Dipole (pIOD) on the evolution of precipitation in China are diagnosed for the period 1950–2013. It is shown that the seasonal responses of precipitation in China to El Niño and pIOD events, and their relationship with the large-scale atmospheric circulations, differ from one season to another. For the pure El Niño years, there is a seasonal reversal of precipitation over southeastern and northwestern China, with deficient precipitation occurring in these two regions before the onset of anomalous wet conditions in the developing autumn. Meanwhile, North China tends to be drier than normal in the developing seasons, but wetter than normal in the decaying seasons. For the pure pIOD events, southern China suffers a precipitation deficit (surplus) in the developing spring (summer and autumn). Furthermore, both North China and northwestern China experience excessive precipitation in the developing autumn and decaying summer. In addition, there is reduced precipitation in northeastern China during both the developing and decaying summers, whereas increased precipitation occurs in the developing autumn and decaying winter. For the combined years, southern China experiences enhanced moisture supply and suffers from increased precipitation from the developing summer through the subsequent spring, but reduced precipitation in the developing spring and decaying summer. Similar to the pure El Niño, northwestern (North) China becomes wetter than normal after the developing summer (autumn) in the combined years. In general, the ENSO/IOD-related precipitation variability could be explained by the associated anomaly circulations.

Description: We analyse Swarm satellite magnetic field and electron density data one month before and one month after 12 strong earthquakes that have occurred in the first 2.5 years of Swarm satellite mission lifetime in the Mediterranean region (magnitude M6.1+) or in the rest of the world (M6.7+). The search for anomalies was limited to the area centred at each earthquake epicentre and bounded by a circle that scales with magnitude according to the Dobrovolsky’s radius. We define the magnetic and electron density anomalies statistically in terms of specific thresholds with respect to the same statistical quantity along the whole residual satellite track (|geomagnetic latitude| ≤ 50°, quiet geomagnetic conditions). Once normalized by the analysed satellite tracks, the anomalies associated to all earthquakes resemble a linear dependence with earthquake magnitude, so supporting the statistical correlation with earthquakes and excluding a relationship by chance.

Description: An original source apportionment study was conducted on atmospheric particles (PM10) collected in Metz, one of the largest cities of Eastern France. A Positive matrix factorization (PMF) analysis was applied to a sampling filter-based chemical dataset obtained for the April 2015 to January 2017 period. Nine factors were clearly identified, showing mainly contributions from anthropogenic sources of primary PM (19.2% and 16.1% for traffic and biomass burning, respectively) as well as secondary aerosols (12.3%, 14.5%, 21.8% for sulfate-, nitrate-, and oxalate-rich factors, respectively). Wood-burning aerosols exhibited strong temporal variations and contributed up to 30% of the PM mass fraction during winter, while primary traffic concentrations remained relatively constant throughout the year. These two sources are also the main contributors during observed PM10 pollution episodes. Furthermore, the dominance of the oxalate-rich factor among other secondary aerosol factors underlines the role of atmospheric processing to secondary organic aerosol loadings which are still poorly characterized in this region. Finally, Concentration-Weighted Trajectory (CWT) analysis were performed to investigate the geographical origins of the apportioned sources, notably illustrating a significant transport of both nitrate-rich and sulfate-rich factors from Northeastern Europe but also from the Balkan region.

Description: We have developed a simple measuring system prototype that uses an unmanned aerial vehicle (UAV) and a non-dispersive infrared (NDIR) analyzer to detect regional carbon dioxide (CO2) concentrations and obtain vertical CO2 distributions. Here, we report CO2 measurement results for the lower troposphere above Ogata Village, Akita Prefecture, Japan (about 40° N, 140° E, approximately −1 m amsl), obtained with this UAV system. The actual flight observations were conducted at 500, 400, 300, 200, 100, and 10 m above the ground, at least once a month during the daytime from February 2018 to February 2019. The raw CO2 values from the NDIR were calibrated by two different CO2 standard gases and high-purity nitrogen (N2) gas (as a CO2 zero gas; 0 ppm). During the observation period, the maximum CO2 concentration was measured in February 2019 and the minimum in August 2018. In all seasons, CO2 concentrations became higher as the flight altitude was increased. The monthly pattern of observed CO2 changes is similar to that generally observed in the Northern Hemisphere as well as to surface CO2 changes simulated by an atmospheric transport model of the Japan Meteorological Agency. It is highly probable that these changes reflect the vegetation distribution around the study area.

Description: Coastal rural areas can be a source of elemental mercury, but the potential influence of their topographic and climatic particularities on gaseous elemental mercury (GEM) fluxes have not been investigated extensively. In this study, gaseous elemental mercury was measured over Mediterranean coastal grassland located in Northern Greece from 2014 to 2015 and GEM fluxes were evaluated utilizing Monin–Obukhov similarity theory. The GEM fluxes ranged from –50.30 to 109.69 ng m−2 h−1 with a mean value equal to 10.50 ± 19.14 ng m−2 h−1. Concerning the peak events, with high positive and low negative GEM fluxes, those were recorded from the morning until the evening. Rain events were a strong contributing factor for enhanced GEM fluxes. The enhanced turbulent mixing under daytime unstable conditions led to greater evasion and positive GEM fluxes, while, during nighttime periods, the GEM evasion is lower, indicating the effect of atmospheric stability on GEM fluxes. The coastal grassland with its specific characteristics influences the GEM fluxes and this area could be characterized as a source of elemental mercury. This study is one of the rare efforts in the research community to estimate GEM fluxes in a coastal natural site based on aerodynamic gradient method.

Description: Volatile organic compounds (VOCs) are important precursors to the formation of ozone and fine particulate matter, the two pollutants of most concern in Sydney, Australia. Despite this importance, there are very few published measurements of ambient VOC concentrations in Australia. In this paper, we present mole fractions of several important VOCs measured during the campaign known as MUMBA (Measurements of Urban, Marine and Biogenic Air) in the Australian city of Wollongong (34°S). We particularly focus on measurements made during periods when clean marine air impacted the measurement site and on VOCs of biogenic origin. Typical unpolluted marine air mole fractions during austral summer 2012-2013 at latitude 34°S were established for CO2 (391.0 ± 0.6 ppm), CH4 (1760.1 ± 0.4 ppb), N2O (325.04 ± 0.08 ppb), CO (52.4 ± 1.7 ppb), O3 (20.5 ± 1.1 ppb), acetaldehyde (190 ± 40 ppt), acetone (260 ± 30 ppt), dimethyl sulphide (50 ± 10 ppt), benzene (20 ± 10 ppt), toluene (30 ± 20 ppt), C8H10 aromatics (23 ± 6 ppt) and C9H12 aromatics (36 ± 7 ppt). The MUMBA site was frequently influenced by VOCs of biogenic origin from a nearby strip of forested parkland to the east due to the dominant north-easterly afternoon sea breeze. VOCs from the more distant densely forested escarpment to the west also impacted the site, especially during two days of extreme heat and strong westerly winds. The relative amounts of different biogenic VOCs observed for these two biomes differed, with much larger increases of isoprene than of monoterpenes or methanol during the hot westerly winds from the escarpment than with cooler winds from the east. However, whether this was due to different vegetation types or was solely the result of the extreme temperatures is not entirely clear. We conclude that the clean marine air and biogenic signatures measured during the MUMBA campaign provide useful information about the typical abundance of several key VOCs and can be used to constrain chemical transport model simulations of the atmosphere in this poorly sampled region of the world.

Description: In this study, AERONET (Aerosol Robotic Network) and EARLINET (European Aerosol Research Lidar Network) data from 17 collocated lidar and sun photometer stations were used to characterize the optical properties of aerosol and their types for the 2008–2018 period in various regions of Europe. The analysis was done on six cluster domains defined using circulation types around each station and their common circulation features. As concluded from the lidar photometer measurements, the typical aerosol particles observed during 2008–2018 over Europe were medium-sized, medium absorbing particles with low spectral dependence. The highest mean values for the lidar ratio at 532 nm were recorded over Northeastern Europe and were associated with Smoke particles, while the lowest mean values for the Angstrom exponent were identified over the Southwest cluster and were associated with Dust and Marine particles. Smoke (37%) and Continental (25%) aerosol types were the predominant aerosol types in Europe, followed by Continental Polluted (17%), Dust (10%), and Marine/Cloud (10%) types. The seasonal variability was insignificant at the continental scale, showing a small increase in the percentage of Smoke during spring and a small increase of Dust during autumn. The aerosol optical depth (AOD) slightly decreased with time, while the Angstrom exponent oscillated between “hot and smoky” years (2011–2015) on the one hand and “dusty” years (2008–2010) and “wet” years (2017–2018) on the other hand. The high variability from year to year showed that aerosol transport in the troposphere became more and more important in the overall balance of the columnar aerosol load.

Description: A long-lived supercell developed in Northwest Bulgaria on 15 May 2018 and inflicted widespread damage along its track. The first part of this article presents a detailed overview of the observed storm evolution. Doppler radar observations reveal that the storm acquired typical supercellular signatures and maintained reflectivity values in excess of 63 dBZ for more than 4 h. The thunderstorm was also analyzed through lightning observations that highlighted important characteristics of the overall supercell dynamics. In its second part, the study investigates the predictability of the severe weather outbreak. In the medium forecast ranges, the global European Centre for Medium-Range Weather Forecasts (ECMWF) ensemble indicated the presence of favorable conditions for the development of deep moist convection 4 days prior to the event. A set of three convection-allowing ensemble simulations also demonstrated that the practical predictability of the supercell was approximately 12 h, which is considerably higher than some previously reported estimates. Nevertheless, the skill of the convective forecasts appears to be limited by the presence of typical model errors, such as the timing of convection initiation and the development of spurious convective activity. The relevance of these errors to the optimal ensemble size and to the design of future convection-allowing numerical weather prediction (NWP) systems is further discussed.

Description: The increasing frequency, intensity and duration of heat waves seem to follow the observed global warming in recent decades. Vulnerability to heat waves is expected to increase in urban environments mainly due to population density and the effect of the urban heat island that make cities hotter than surrounding non-urban areas. The present study focuses on a vulnerable area of the eastern Mediterranean, already characterized as a ‘hot spot’ with respect to heat-related risk and investigates the change in heat stress levels during heat wave compared to non-heat wave conditions as well as the way that heat stress levels respond to heat waves in urban, compared to non-urban, environments. The adoption of a metric accounting for both the intensity and duration of the hot event yielded a total of 46 heat wave episodes over a nearly 60-year period, but with very rare occurrence until the late 1990s and a profound increased frequency thereafter. The results reveal a difference of at least one thermal stress category between heat wave and non-heat wave periods, which is apparent across the entire range of the thermal stress distribution. The analysis demonstrates a robust intensification of nighttime heat stress conditions in urban, compared to non-urban, sites during severe heat waves. Nevertheless, severe heat waves almost equalize heat stress conditions between urban and non-urban sites during midday.

Description: The modular structure of the boundary layer and convection parameterizations in atmospheric models have long been affecting the numerical representation of subgrid-scale motions and their mutual interactions. A promising alternative, the eddy-diffusivity/mass-flux approach (EDMF), has the potential for unifying the existing formulations into a consistent scheme and improving some of the long-standing issues. This study documents a step towards developing such a unified approach by implementing a stochastic multi-plume EDMF scheme into the Community Atmosphere Model (Version 5.0). Its performance in single-column mode is evaluated against the control parameterization and large-eddy simulation (LES) for two benchmark cases: marine and continental shallow convection. Overall, the results for the two parameterizations agree well with each other and with LES in terms of mean profiles of moist conserved variables and their vertical fluxes, as well as the updraft properties. However, systematic differences between the two schemes, especially for transient continental convection, are also documented. Using EDMF helps improve some of the parameterized features of shallow convection. In particular, for the highest tested vertical resolution, the EDMF cloud base and top heights and the vertical fluxes of energy and water are remarkably close to LES.

Description: In this study, using Aerosol Robotic Network aerosol optical depth (AOD) products at three stations in the North China Plain (NCP)—a heavily polluted region in China—the AOD products from six satellite-borne radiometers: the Moderate Resolution Imagining Spectroradiometer (MODIS), the Multiangle Imaging Spectroradiometer (MISR), Ozone Mapping Imaging (OMI), the Visible Infrared Imaging Radiometer (VIIRS), the Sea-Viewing Wide Field-of-View Sensor (SeaWiFS), and Polarization and Directionality of the Earth’s Reflectances (POLDER), were thoroughly validated, shedding new light on their advantages and disadvantages. The MODIS Deep Blue (DB) products provide more accurate retrievals than the MODIS Dark Target (DT) and other satellite products at the Beijing site (BJ, a megacity), with higher correlations with AERONET (R 〉 0.93), lower mean absolute bias (MB 〈 0.012), and higher percentages (〉68%) falling within the expected error (EE). All MODIS DT and DB products perform better than the other satellite products at the Xianghe site (XH, a suburb). The MODIS/Aqua DT products at both 3-km and 10-km resolutions performed better than the other space-borne AOD products at the Xinglong site (XL, a rural area at the top of a mountain). MISR, VIIRS, and SeaWiFS tend to underestimate high AOD values and overestimate AOD values under very low AOD conditions in the NCP. Both OMI and POLDER significantly underestimate the AOD. In terms of data volume, MISR with the limited swath width of 380 km has less data volume than the other satellite sensors. MODIS products have the highest sampling rate, especially the MODIS DT and DB merged products, and can be used for various climate study and air-quality monitoring.

Description: The Asia Pacific Mercury Monitoring Network (APMMN) cooperatively measures mercury in precipitation in a network of sites operating in Asia and the Western Pacific region. The network addresses significant data gaps in a region where mercury emission estimates are the highest globally, and available measurement data are limited. The reduction of mercury emissions under the Minamata Convention on Mercury also justifies the need for continent-wide and consistent observations that can help determine the magnitude of the problem and assess the efficacy of reductions over time. The APMMN’s primary objectives are to monitor wet deposition and atmospheric concentrations of mercury and assist partners in developing their own monitoring capabilities. Network planning began in 2012 with wet deposition sampling starting in 2014. Currently, eight network sites measure mercury in precipitation following standardized procedures adapted from the National Atmospheric Deposition Program. The network also has a common regional analytical laboratory (Taiwan), and quality assurance and data flagging procedures, which ensure the network makes scientifically valid and consistent measurements. Results from our ongoing analytical and field quality assurance measurements show minimal contamination in the network and accurate analytical analyses. We are continuing to monitor a potential concentration and precipitation volume bias under certain conditions. The average mercury concentration in precipitation was 11.3 (+9.6) ng L−1 for 139 network samples in 2018. Concentrations for individual sites vary widely. Low averages compare to the low concentrations observed on the U.S. West Coast; while other sites have average concentrations similar to the high values reported from many urban areas in China. Future APMMN goals are to (1) foster new network partnerships, (2) continue to collect, quality assure, and distribute results on the APMMN website, (3) provide training and share best monitoring practices, and (4) establish a gaseous concentration network for estimating dry deposition.

Description: Studies in actual urban settings that integrate chemical reaction modeling, radiation, and particular emissions are mandatory to evaluate the effects of traffic-related air pollution on street canyons. In this paper, airflow patterns and reactive pollutant behavior for over 24 h, in a realistic urban canyon in Osaka City, Japan, was conducted using a computational fluid dynamics (CFD) model coupled with a chemical reaction model (CBM-IV). The boundary conditions for the CFD model were obtained from mesoscale meteorological and air quality models. Inherent street canyon processes, such as ground and wall radiation, were evaluated using a surface energy budget model of the ground and a building envelope model, respectively. The CFD-coupled chemical reaction model surpassed the mesoscale models in describing the NO, NO2, and O3 transport process, representing pollutants concentrations more accurately within the street canyon since the latter cannot capture the local phenomena because of coarse grid resolution. This work showed that the concentration of pollutants in the urban canyon is heavily reliant on roadside emissions and airflow patterns, which, in turn, is strongly affected by the heterogeneity of the urban layout. The CFD-coupled chemical reaction model characterized better the complex three-dimensional site and hour-dependent dispersion of contaminants within an urban canyon.

Description: Ocean salinity variability provides a new way to study the evolution of the the El Niño-Southern Oscillation (ENSO). Comparisons between the salinity variation and related processes responsible for sea surface temperature anomaly (SSTA) were extensively examined for the two strong El Niño (EN) events in 1997/1998 and 2015/2016, and a special EN event in 2014/2015. The results show that the development of EN is significantly correlated with a sea surface salinity anomaly (SSSA) in the tropical western-central Pacific. In the spring of 1997 and 2015 with strong EN events, the western-central equatorial Pacific exhibited significant negative SSSA that propagated eastward to the west of the dateline. The negative SSSA induced increased barrier layer thickness (BLT) which enhanced sea surface temperature (SST) warming in the tropical central Pacific. In contrast, although a negative SSSA occurred during April of the 2014/2015 weak EN event in the western-central equatorial Pacific, this SSSA was mainly confined to between 160° E and 180° E without significant eastward movement, resulting in a weakened BLT thickening process and a weak modulation effect on SST. We also confirm that the surface forcing associated with fresh water flux (FWF: evaporation (E) minus precipitation (P)) plays a prominent role in the surface salinity tendency in the tropical Pacific during EN events. Moreover, the negative FWF anomaly leads a strong negative SSSA by two months. Compared with the two strong ENs, the early negative FWF anomaly in the weak 2014/2015 EN did not present distinct development and eastward propagation and weakened rapidly in the summer of 2015. We demonstrate that change in salinity can modulate the ENSO, and the variation of SSSA and associated physical processes in the tropical western-central Pacific and could be used as an indicator for predicting the development of ENSO.

Description: Contribution of extra-tropical synoptic cyclones to the formation of mean summer atmospheric circulation patterns in the Arctic domain (≥60° N) was investigated by clustering dominant Arctic circulation patterns based on daily mean sea-level pressure using self-organizing maps (SOMs). Three SOM patterns were identified; one pattern had prevalent low-pressure anomalies in the Arctic Circle (SOM1), while two exhibited opposite dipoles with primary high-pressure anomalies covering the Arctic Ocean (SOM2 and SOM3). The time series of their occurrence frequencies demonstrated the largest inter-annual variation in SOM1, a slight decreasing trend in SOM2, and the abrupt upswing after 2007 in SOM3. Analyses of synoptic cyclone activity using the cyclone track data confirmed the vital contribution of synoptic cyclones to the formation of large-scale patterns. Arctic cyclone activity was enhanced in the SOM1, which was consistent with the meridional temperature gradient increases over the land–Arctic ocean boundaries co-located with major cyclone pathways. The composite daily synoptic evolution of each SOM revealed that all three SOMs persisted for less than five days on average. These evolutionary short-term weather patterns have substantial variability at inter-annual and longer timescales. Therefore, the synoptic-scale activity is central to forming the seasonal-mean climate of the Arctic.

Description: The divergence problem, which manifests as an unstable response relationship between tree-ring growth and climatic factors under the background of global warming, poses a challenge to both the traditional theory of dendroclimatology and the reliability of climatic reconstructions based on tree-ring data. Although Schrenk spruce, as the dominant tree species in the Tianshan Mountains, is frequently applied in the dendrochronological studies, the understanding of the divergence problem of this tree species is still limited. This study conducted correlation analysis between climatic factors and tree-ring width chronologies from 51 living and healthy specimens of Schrenk spruce at sites of high and low elevation in the Alatau Mountains to determine the stability of the response. The results revealed that the tree-ring width of the spruce specimens was correlated positively with precipitation and correlated negatively with temperature. Although the variations of the two tree-ring chronologies were similar, the radial growth of the spruce at the low elevation was found more sensitive to climatic factors. Furthermore, the sensitivity of tree growth to climate demonstrated an obvious increase after an abrupt change of climate under the background of the recent warming and wetting trend. Increased drought stress, calculated based on climatic data, was regarded as the main reason for this phenomenon. The results supply the gap of the stability of climatic response of tree growth in Central Asia to some extent.

Description: The measured concentrations of inorganic pollutants, such as ozone (2015–2018), sulfur, and nitrogen oxides (2012–2018) at air monitoring sites in the south of Eastern Siberia were sampled, following the passive sampling method, and analyzed. The spatial inhomogeneity of atmospheric gas concentrations is presented. The ozone concentration is lower in urban areas than those in rural areas and the background level. However, the nitrogen and sulfur oxide concentrations are higher in the atmosphere over the city site. The seasonal dependence of the ozone concentration was determined using its maximum (March–April) and minimum (September–October) levels. The dynamics of the nitrogen and sulfur oxide concentrations indicate that they are at their highest in December–June and their lowest in July–August. To verify the validity of the pollutant concentration measurements sampled by passive sampling, we compared our results with those obtained following the automatic and filter pack methods. A linear regression analysis and a pairwise modification of Student’s t test evaluated the concentrations of the air pollutant, sampled and measured using different methods, and they correlate well (r = 0.7–0.9). Full validation of the passive sampling method is not possible for some sites; therefore it is necessary to remove the remaining systematic errors in future work.

Description: This study aims to address hydrological processes and impacts of an atmospheric river (AR) event that occurred during 15–18 February 2004 in the Russian River basin in California. The National Water Model (NWM), a fully distributed hydrologic model, was used to evaluate the hydrological processes including soil moisture flux, overland flow, and streamflow. Observed streamflow and volumetric soil water content data were used to evaluate the performance of the NWM using various error metrics. The simulation results showed that this AR event (15–18 February 2004) with a long duration of precipitation could cause not only deep soil saturation, but also high direct runoff depth. Taken together, the analysis revealed the complex interaction between precipitation and land surface response to the AR event. The results emphasize the significance of a change of water contents in various soil layers and suggest that soil water content monitoring could aid in improving flood forecasting accuracy caused by the extreme events such as the AR.

Description: In this study, based on daily precipitation records during 1958–2017 from 28 meteorological stations in the Beijing-Tianjin-Hebei (BTH) region, the spatio-temporal variations in precipitation extremes defined by twelve indices are analyzed by the methods of linear regression, Mann-Kendall test and continuous wavelet transform. The results showed that the spatial patterns of all the indices except for consecutive dry days (CDD) and consecutive wet days (CWD) were similar to that of annual total precipitation with the high values in the east and the low value in the west. Regionally averaged precipitation extremes were characterized by decreasing trends, of which five indices (i.e., very heavy precipitation days (R50), very wet precipitation (R95p), extreme wet precipitation (R99p), max one-day precipitation (R × 1day), and max five-day precipitation (R × 5day)) exhibited significantly decreasing trends at 5% level. From monthly and seasonal scale, almost all of the highest values in R × 1day and R × 5day occurred in summer, especially in July and August due to the impacts of East Asian monsoon climate on inter-annual uneven distribution of precipitation. The significant decreasing trends in annual R×1day and R×5day were mainly caused by the significant descend in summer. Besides, the possible associations between precipitation extremes and large-scale climate anomalies (e.g., ENSO (El Niño Southern Oscillation), NAO (North Atlantic Oscillation), IOD (Indian Ocean Dipole), and PDO (Pacific Decadal Oscillation)) were also investigated using the correlation analysis. The results showed that the precipitation extremes were significantly influenced by ENSO with one-year ahead, and the converse correlations between the precipitation extremes and climate indices with one-year ahead and 0-year ahead were observed. Moreover, all the indices show significant two- to four-year periodic oscillation during the entire period of 1958–2017, and most of indices show significant four- to eight-year periodic oscillation during certain periods. The influences of climate anomalies on precipitation extremes were composed by different periodic components, with most of higher correlations occurring in low-frequency components.

Description: The recent development of high-resolution climate models offers a promising approach in improving the simulation of precipitation, clouds and temperature. However, higher grid spacing is also a promising feature to improve the simulation of snow cover. In particular, it provides a refined representation of topography and allows for an explicit simulation of convective precipitation processes. In this study we analyze the snow cover in a set of decade-long high-resolution climate simulation with horizontal grid spacing of 2.2 km over the greater Alpine region. Results are compared against observations and lower resolution models (12 and 50 km), which use parameterized convection. The simulations are integrated using the COSMO (Consortium for Small-Scale Modeling) model. The evaluation of snow water equivalent (SWE) in the simulation of present-day climate, driven by the ERA-Interim reanalysis, against an observational dataset, reveals that the high-resolution simulation clearly outperforms simulations with grid spacing of 12 and 50 km. The latter simulations underestimate the cumulative amount of SWE over Switzerland over the whole annual cycle by 33% (12 km simulation) and 56% (50 km simulation) while the high-resolution simulation shows a spatially and temporally averaged difference of less than 1%. Scenario simulations driven by GCM MPI-ESM-LR (2081–2090 RCP8.5 vs. 1991–2000) reveal a strong decrease of SWE over the Alps, consistent with previous studies. Previous studies had found that the relative decrease becomes gradually smaller with elevation, but this finding was limited to low and intermediate altitudes (as a 12 km simulation resolves the topography up to 2500 m). In the current study we find that the height gradient reverses sign, and relative reductions in snow cover increases above 3000 m asl, where important parts of the cryosphere are present. In addition, the simulations project a transition from permanent to seasonal snow cover at high altitudes, with potentially important impacts to Alpine permafrost. This transition and the more pronounced decline of SWE emphasize the value of the higher grid spacing. Overall, we show that high-resolution climate models offer a promising approach in improving the simulation of snow cover in Alpine terrain.

Description: China is experiencing severe PM2.5 (fine particles with a diameter of 2.5 μg or smaller) pollution problem. Little is known, however, about how the increasing concentration trend is spatially distributed, nor whether there are some areas that experience a stable or decreasing concentration trend. Managers and policymakers require such information to make strategic decisions and monitor progress towards management objectives. Here, we present a pixel-based linear trend analysis of annual PM2.5 concentration variation in China during the period 1999–2016, and our results provide guidance about where to prioritize management efforts and affirm the importance of controlling coal energy consumption. We show that 87.9% of the whole China area had an increasing trend. The drastic increasing trends of PM2.5 concentration during the last 18 years in the Beijing–Tianjin–Hebei region, Shandong province, and the Three Northeastern Provinces are discussed. Furthermore, by exploring regional PM2.5 pollution, we find that Tarim Basin endures a high PM2.5 concentration, and this should have some relationship with oil exploration. The relationship between PM2.5 pollution and energy consumption is also discussed. Not only energy structure reconstruction should be repeatedly emphasized, the amount of coal burned should be strictly controlled.

Description: A 10-year study of elevated severe thunderstorms was performed using The National Centers for Environmental Information Storm Events Database. A total of 80 elevated thunderstorm cases were identified, verified, and divided into “Prolific” and “Marginal” classes. These severe cases occurred at least 80 km away from, and on the cold side of, a surface boundary. The downdraft convective available potential energy (DCAPE), downdraft convective inhibition (DCIN), and their ratio are tools to help estimate the potential for a downdraft to penetrate through the depth of a stable surface layer. The hypothesis is that as the DCIN/DCAPE ratio decreases, there exists enhanced possibility of severe surface winds. Using the initial fields from the Rapid Refresh numerical weather prediction model, datasets of DCIN, DCAPE, and their ratio were created. Mann-Whitney U tests on the Prolific versus Marginal case sets were undertaken to determine if the DCAPE and DCIN values come from different populations for the two different case sets. Results show that the Prolific cases have values of DCIN closer to zero, suggesting the downdraft is able to penetrate to the surface causing severe winds. Thus, comparing DCIN and DCAPE is a viable tool in determining if downdrafts will reach the surface from elevated thunderstorms.

Description: This research aims to find the most ecological itineraries for urban mobility in a small city (eco-routes), where distances are rather short, but car dependence is really high. A real life citywide survey was carried out in the city of Caceres (Spain) with almost 100,000 inhabitants. Research was done on alternating routes, traffic, times of day, and weather conditions. The output of the study was to assess fuel consumption, CO2, and regulated pollutant emissions for different type of vehicles, routes, and drivers. The results show that in the case studied, urban roads had fewer emissions (CO2 and pollutants) but there was an increase in the population affected by pollutants. On the contrary, bypasses reduced travel time and congestion but increased fuel consumption and emissions. Traffic conditions had a greater influence on fuel consumption in petrol vehicles than diesel ones. Therefore, there must be a balanced distribution of traffic in order to minimize congestion, and at the same time to reduce emissions and the number of people affected by harmful pollution levels. There should be a combination of regulatory measures in traffic policies in order to achieve that balance by controlling access to city centres, limiting parking spaces, pedestrianization, and lowering traffic speeds in sensitive areas.

Description: The impact of absorbing aerosols on climate is complex, with their potential positive or negative forcing, depending on many factors, including their height distribution and reflective properties of the underlying background. Measurement data is very limited, due to insufficient remote sensing methods dedicated to the retrieval of their vertical distribution. Columnar values of absorbing aerosol optical depth (AAOD) and single scattering albedo (SSA) are retrieved by the Aerosol Robotic Network (AERONET). However, the number of available results is low due to sky condition and aerosol optical depth (AOD) limitation. Presented research describes results of field campaigns in Strzyżów (South-East Poland, Eastern Europe) dedicated to the comparison of the absorption coefficient and SSA measurements performed with on-ground in-situ devices (aethalomter, nephelometer), small unmanned aerial system (UAS) carrying micro-aethalometer, as well as with lidar/ceilometer. An important aspect is the comparison of measurement results with those delivered by AERONET. Correlation of absorption to scattering coefficients measured on ground (0.79) and correlation of extinction on ground to AOD measured by AERONET (0.77) was visibly higher than correlation between AOD and AAOD retrieved by AERONET (0.56). Columnar SSA was weakly correlated with ground SSA (higher values of columnar SSA), which were mainly explained by hygroscopic effects, increasing scattering coefficient in ambient (wet conditions), and partly high uncertainty of SSA retrieval. AAOD derived with the use of profiles from UAS up to PBL height, was estimated to contribute in average to 37% of the total AAOD. A method of AAOD estimation, in the whole troposphere, with use of measured vertical profiles of absorption coefficient and extinction coefficient profiles from lidars was proposed. AAOD measured with this method has poor correlation with AERONET data, however for some measurements, within PBL, AAOD was higher than reported by AERONET, suggesting potential underestimation in photometric measurement under particular conditions. Correlation of absorption coefficient in profile to on ground measurements decrease with altitude. Measurements of SSA from drones agree well with ground measurements and are lower than results from AERONET, which suggests a larger contribution of absorbing aerosols. As an alternative for AAOD estimation in case of lack of AERONET AAOD data simple models are proposed, which base on AOD scaling with SSA measured with different methods. Proposed solution increase potential of absorption coefficient measurements in vertical profiles and columns of the atmosphere. Presented solutions make measurements of absorption coefficients in vertical profiles more affordable and allow rough estimation of columnar values for the whole atmosphere.

Description: Nature has been adversely affected by increasing industrialization, especially during the latter part of the last century, as a result of accelerating technological development, unplanned urbanization, incorrect agricultural policies and deforestation, which have contributed to the elevated concentration of the greenhouse gases (GHGs) in the environment. GHG accumulation has an adverse impact on meteorological and hydro-meteorological parameters, particularly temperature. Temperature plays a prominent and well-known role in evaporation, transpiration and changes in water demand, and thus significantly affects both water availability and food security. Therefore, a systematic understanding of temperature is important for fighting food insecurity and household poverty. Variations in temperature are often assessed and characterized through trend analysis. Hence, the objective of this paper is to determine long-term trends in mean monthly maximum and minimum air temperatures for the Maloti-Drakensberg region. The Mann–Kendall test, a non-parametric test, was applied on mean air temperature for the 1960–2016 period. A significant rising trend (p 〈 0.001) was detected with a yearly change in the long term annual mean maximum and mean minimum temperature by 0.03 °C/annum and 0.01 °C/annum, respectively. This knowledge has important implications for both the state of the environment and livelihoods in the region, since its use can be useful in planning and policymaking in water resource management, biodiversity conservation, agriculture, tourism and other sectors of the economy within the region.

Description: Subway air pollution mainly refers to inhalable particulate matter (PM) pollution, organic pollution, and microbial pollution. Based on the investigation and calculation of the existing researches, this paper summarizes the sources of air pollutants, chemical compositions, and driving factors of PM variations in subway. It evaluates the toxicity and health risks of pollutants. In this paper, the problems and challenges during the deployment of air pollution governance are discussed. Results show that the global PM compliance rate of subway is about 30%. Subway air pollution is endogenous, which means that pollutants mainly come from mechanical wear and building materials erosions. Particles are mainly metal particles, black carbon, and floating dust. The health risks of some chemical elements in the subway have reached critical levels. The variations of PM concentrations show spatial-temporal characteristics, which are mainly controlled by train age, brakes types, and environmental control systems. The authors then analyze the dynamics of interactions among government, companies and public during the air pollution governance by adding the following questions: (a) who pays the bill; (b) how to evaluate the cost-effectiveness of policies; (c) how the public moves from risk perception to actions; (d) how to develop clean air technology better so as to ultimately incentivize stakeholders and to facilitate the implementation of subway clean air programme in a resilient mode.

Description: Prescribed pasture burning plays a critical role in ecosystem maintenance in tallgrass prairie ecosystems and may contribute to agricultural productivity but can also have negative impacts on air quality. Volatile organic compound (VOC) concentrations were measured immediately downwind of prescribed tallgrass prairie fires in the Flint Hills region of Kansas, United States. The VOC mixture is dominated by alkenes and oxygenated VOCs, which are highly reactive and can drive photochemical production of ozone downwind of the fires. The computed emission factors are comparable to those previous measured from pasture maintenance fires in Brazil. In addition to the emission of large amounts of particulate matter, hazardous air pollutants such as benzene and acrolein are emitted in significant amounts and could contribute to adverse health effects in exposed populations.

Description: The same prescribed anthropogenic aerosol forcing was implemented into three climate models. The atmosphere components of these participating climate models were the GAMIL, ECHAM, and CAM models. Ensemble simulations were carried out to obtain a reliable estimate of anthropogenic aerosol effective radiative forcing (ERF). The ensemble mean ERFs from these three participating models with this aerosol forcing were −0.27, −0.63, and −0.54 W∙m−2. The model diversity in ERF is clearly reduced as compared with those based on the models’ own default approaches (−1.98, −0.21, and −2.22 W∙m−2). This is consistent with the design of this aerosol forcing. The modeled ERF can be decomposed into two basic components, i.e., the instantaneous radiative forcing (RF) from aerosol–radiation interactions (RFari) and the aerosol-induced changes in cloud forcing (△Fcloud*). For the three participating models, the model diversity in RFari (−0.21, −0.33, and −0.29 W∙m−2) could be constrained by reducing the differences in natural aerosol radiative forcings. However, it was difficult to figure out the reason for the model diversity in △Fcloud* (−0.05, −0.28, and −0.24 W∙m−2), which was the dominant source of the model diversity in ERF. The variability of modeled ERF was also studied. Ensemble simulations showed that the modeled RFs were very stable. The rapid adjustments (ERF − RF) had an important role to play in the quantification of the perturbation of ERF. Fortunately, the contribution from the rapid adjustments to the mean ERF was very small. This study also showed that we should pay attention to the difference between the aerosol climate effects we want and the aerosol climate effects we calculate.

Description: Residential coal combustion is one of the most significant sources of carbonaceous aerosols in the Highveld region of South Africa, significantly affecting the local and regional climate. This study investigated single coal-burning particles emitted when using different fire-ignition techniques (top-lit up-draft versus bottom-lit up-draft) and air ventilation rates (defined by the number of air holes above and below the fire grate) in selected informal braziers. Aerosol samples were collected on nucleopore filters at the Sustainable Energy Technology and Research Centre Laboratory, University of Johannesburg. The individual particles (~700) were investigated using a scanning electron microscope equipped with energy-dispersive X-ray spectroscopy (EDX). Two distinct forms of spherical organic particles (SOPs) were identified, one less oxidized than the other. The particles were further classified into electronically dark and bright. The EDX analysis showed that 70% of the dark spherical organic particles had higher (~60%) relative oxygen content than in the bright SOPs. The morphology of spherical organic particles were quantified and classified into four categories: ~50% were bare single particles; ~35% particles were aggregated and formed diffusion accretion chains; 10% had inclusions, and 5% were deformed due to impaction on filter material during sampling. This study concludes that there are two distinct kinds of coal burning spherical organic particles and that dark SOPs are less volatile than bright SOPs. The authors also show that these spherical organic particles are similar in nature and characteristics to tar balls observed in biomass combustion and that they have the potential to absorb sunlight thereby affecting the earth’s radiative budget and climate. This study provides insights into the mixing states, morphology, and possible formation mechanisms of these organic particles from residential coal combustion in informal stoves.

Description: We present the framework of the modeling package IonoSeis. This software models Global Navigation Satellite System (GNSS) derived slant total electron content (sTEC) perturbations in the ionosphere due to the interaction of the neutral atmosphere and charged particles in the ionosphere. We use a simplified model to couple the neutral particle momentum into the ionosphere and reconstruct time series of sTEC perturbations that match observed data in both arrival time and perturbation shape. We propagate neutral atmosphere disturbances to ionospheric heights using a three-dimensional ray-tracing code in spherical coordinates called Windy Atmospheric Sonic Propagation (WASP3D), which works for a stationary or non-stationary atmospheric models. The source of the atmosphere perturbation can be an earthquake or volcanic eruption; both couple significant amounts of energy into the atmosphere in the frequency range of a few Millihertz. We demonstrate the output of the code by comparing modeled sTEC perturbation data to the observed perturbation recorded at GNSS station BTNG (Bitung, Indonesia) immediately following the 28 September 2018, Sulawesi-Palu earthquake. With this framework, we provide a software to couple the lithosphere, atmosphere, and ionosphere that can be used to study post-seismic ionospherically-derived signals.

Description: The present study explores the cloud microphysics (MPs) impact on the simulation of two convective rainfall events (CREs) over the complex topography of Andes mountains, using the Weather Research and Forecasting- Advanced Research (WRF-ARW) model. The events occurred on December 29 2015 (CRE1) and January 7 2016 (CRE2). Six microphysical parameterizations (MPPs) (Thompson, WSM6, Morrison, Goddard, Milbrandt and Lin) were tested, which had been previously applied in complex orography areas. The one-way nesting technique was applied to four domains, with horizontal resolutions of 18, 6, and 3 km for the outer ones, in which cumulus and MP parameterizations were applied, while for the innermost domain, with a resolution of 0.75 km, only MP parameterization was used. It was integrated for 36 h with National Centers for Environmental Prediction (NCEP Final Operational Global Analysis (NFL) initial conditions at 00:00 UTC (Coordinated Universal Time). The simulations were verified using Geostationary Operational Environmental Satellites (GOES) brightness temperature, Ka band cloud radar, and surface meteorology variables observed at the Huancayo Observatory. All the MPPs detected the surface temperature signature of the CREs, but for CRE2, it was underestimated during its lifetime in its vicinity, matching well after the simulated event. For CRE1, all the schemes gave good estimations of 24 h precipitation, but for CRE2, Goddard and Milbrandt underestimated the 24 h precipitation in the inner domain. The Morrison and Lin configurations reproduced the general dynamics of the development of cloud systems for the two case studies. The vertical profiles of the hydrometeors simulated by different schemes showed significant differences. The best performance of the Morrison scheme for both case studies may be related to its ability to simulate the role of graupel in precipitation formation. The analysis of the maximum reflectivity field, cloud top distribution, and vertical structure of the simulated cloud field also shows that the Morrison parameterization reproduced the convective systems consistently with observations.

Description: In recent years, new and strict air quality regulations have been implemented in China. Therefore, it is of great significance to evaluate the current air pollution situation and effectiveness of actions. In this study, Ozone Monitoring Instrument (OMI) satellite data were used to detect the spatiotemporal characteristics of tropospheric NO2 columns over China from 2005 to 2018, including spatial distribution, seasonal cycles and long-term trends. The averaged NO2 pollution is higher in southeastern China and lower in the northwest, which are well delineated by the Heihe–Tengchong line. Furthermore, the NO2 loadings are highest in the North China Plain, with vertical column density (VCD) exceeding 13 × 1015 molec cm−2. Regarding the seasonal cycle, the NO2 loadings in eastern China is highest in winter and lowest in summer, while the western region shows the opposite feature. The amplitude of annual range increase gradually from the south to the north. If the entire period of 2005–2018 is taken into account, China has experienced little change in NO2. In fact, however, there appears to be significant trends of an increase followed by a downward tendency, with the turning point in the year 2012. In the former episode of 2005–2012, increasing trends overwhelm nearly the whole nation, especially in the Jing–Jin–Tang region, Shandong Province, and Northern Henan and Southern Hebei combined regions, where the rising rates were as high as 1.0–1.8 × 1015 molec cm−2 year−1. In contrast, the latter episode of 2013–2018 features remarkable declines in NO2 columns over China. Particularly, the regions where the decreased degree was remarkable in 2013–2018 were consistent with the regions where the upward trend was obvious in 2005–2012. Overall, this upward–downward pattern is true for most parts of China. However, some of the largest metropolises, such as Beijing, Shanghai and Guangzhou, witnessed a continuous decrease in the NO2 amounts, indicating earlier and more stringent measures adopted in these areas. Finally, it can be concluded that China’s recent efforts to cut NO2 pollution are successful, especially in mega cities.

Description: The use of low-cost sensors for air quality measurements is expanding rapidly, with an associated rise in the number of citizens measuring air quality themselves. This has major implications for traditional air quality monitoring as performed by Environmental Protection Agencies. Here we reflect on the experiences of the Dutch Institute for Public Health and the Environment (RIVM) with the use of low-cost sensors, particularly NO2 and PM10/PM2.5-sensors, and related citizen science, over the last few years. Specifically, we discuss the Dutch Innovation Program for Environmental Monitoring, which comprises the development of a knowledge portal and sensor data portal, new calibration approaches for sensors, and modelling and assimilation techniques for incorporating these uncertain sensor data into air pollution models. Finally, we highlight some of the challenges that come with the use of low-cost sensors for air quality monitoring, and give some specific use-case examples. Our results show that low-cost sensors can be a valuable addition to traditional air quality monitoring, but so far, their use in official monitoring has been limited. More research is needed to establish robust calibration methods while ongoing work is also aimed at a better understanding of the public’s needs for air quality information to optimize the use of low-cost sensors.

Description: Prokaryotic microbes can become aerosolized and deposited into new environments located thousands of kilometers away from their place of origin. The Mediterranean Sea is an oligotrophic to ultra-oligotrophic marginal sea, which neighbors northern Africa (a major source of natural aerosols) and Europe (a source of mostly anthropogenic aerosols). Previous studies demonstrated that airborne bacteria deposited during dust events over the Mediterranean Sea may significantly alter the ecology and function of the surface seawater layer, yet little is known about their abundance and diversity during ‘background’ non-storm conditions. Here, we describe the abundance and genetic diversity of airborne bacteria in 16 air samples collected over an East-West transect of the entire Mediterranean Sea during non-storm conditions in April 2011. The results show that airborne bacteria represent diverse groups with the most abundant bacteria from the Firmicutes (Bacilli and Clostridia) and Proteobacteria (Alphaproteobacteria, Betaproteobacteria, and Gammaproteobacteria) phyla. Most of the bacteria in our samples have previously been observed in the air at other open ocean locations, in the air over the Mediterranean Sea during dust storms, and in the Mediterranean seawater. Airborne bacterial abundance ranged from 0.7 × 104 to 2.5 × 104 cells m−3 air, similar to abundances at other oceanic regimes. Our results demonstrate that airborne bacterial diversity is positively correlated with the mineral dust content in the aerosols and was spatially separated between major basins of the Mediterranean Sea. To our knowledge, this is the first comprehensive biogeographical dataset to assess the diversity and abundance of airborne microbes over the Mediterranean Sea. Our results shed light on the spatiotemporal distribution of airborne microbes and may have implications for dispersal and distribution of microbes (biogeography) in the ocean.

Description: This work focuses on the interaction of a Coronal Mass Ejection (CME) with the ambient solar corona, by studying the spatial and temporal evolution of the density fluctuations observed by the SOHO/UV Coronagraph Spectrometer (UVCS) during the CME. The investigation is performed by applying a wavelet analysis to the HI Ly α 1216 Å line intensity fluctuations observed with UVCS during the CME occurred on 24 December 2006. Strong and coherent fluctuations, with a significant spatial periodicity of about 84 Mm ≃ 0.12 R ⊙ , are shown to develop in about an hour along the front of the CME. The results seem to indicate the Rayleigh-Taylor (RT) instability, susceptible to the deceleration of the heavier fluid of the CME front into the lighter surrounding coronal plasma, as the likely mechanism underlying the generation of the observed plasma fluctuations. This could be the first inference of the RT instability in the outer solar corona in UV, due to the transit of a CME front in the quiet coronal plasma; this interpretation is also supported by a linear magnetohydrodynamic analysis of the RT instability.

Description: Exploring the mechanism that influences the choice of urban public travel mode is an important policy research topic that can promote urban residents’ pro-environment travel (PET) behaviour and relieve the pressure on urban traffic and environmental problems. By expanding the theory of planned behaviour by considering the effects of the quality of public transport service and individual behaviour, this paper establishes a mixed PET behaviour model. Grounded theory is used to analyse data obtained from in-depth interviews, with the aim of determining the relationships among different attributes of the quality of public transport service and PET. An empirical examination in the form of a questionnaire was conducted in Changsha, China, to obtain the intensity and mechanism of various factors influencing pro-environment behaviour decision-making. The results reveal three new pieces of information. First, the influence of many psychological variables (except subjective norms) is consistent with the prediction results of the theory of planned behaviour (TPB), and the predictions of the model are accurate. More specifically, intention (0.535) and habit (0.354) are key factors in PET behaviour, while attitude (0.527) has the most significant effect on intention towards PET behaviour. Second, the perceived service quality of public transport has a direct and significant impact on intention towards PET behaviour. Satisfaction with public transport service quality exerts a mediating effect on perceived service quality and PET behaviour. More specifically, operation and management (0.808) and vehicle environment (0.809) have the most important influence on intention towards PET behaviour. Last but not least, the extent of the influence of PET behaviour varies based on travellers’ demographic characteristics. The driving age, income and ownership of private cars show the greatest impact. The perceived service quality of public transport and travellers’ social and economic characteristics all play roles in the psychology of travel decisions, and are associated with PET behaviour on several distinct levels. From the perspectives of passenger psychology, public transport service quality and personal attributes of passengers, this paper provides a scientific basis for decision-making in transportation systems and the formulation of traffic intervention strategies to promote voluntary public reductions in carbon-intensive travel behaviour.

Description: Four cross-coupled models were used to investigate the relative contributions of atmospheric and oceanic components to the asymmetry of the El Niño–Southern Oscillation (ENSO). Strong El Niño and La Niña events related to the negative heat flux feedbacks were found to be determined mainly by the atmospheric component, and the stronger sea surface temperature (SST) anomalies in the warm phase did not lead to an increased SST asymmetry. The skewness of the four models could be affected by both atmospheric and oceanic components; the atmospheric component determines the strength of positive and negative SST anomalies, and the oceanic component affects the strength of the negative SST anomalies in the cold phase under the same atmospheric component group. The Bjerknes stability index (BJ index) of warm and cold phases contributed to the El Niño–La Niña SSTA asymmetries in observation, but the BJ index did not necessarily explain the El Niño–La Niña SSTA asymmetries in climate model simulations. The SST asymmetries in these four models were closely associated with convective precipitation and wind stress asymmetries, which are also determined by both the atmospheric and oceanic components.

Description: In ultrasonic equipment (anemometers and thermometers), for the measurement of parameters of atmospheric turbulence, a standard algorithm that calculates parameters from temporary structural functions constructed on the registered data is usually used. The algorithm is based on the Kolmogorov–Obukhov law. The experience of using ultrasonic meters shows that such an approach can lead to significant errors. Therefore, an improved algorithm for calculating the parameters is developed, which allows more accurate estimation of the structural characteristics of turbulent fluctuations, with an error that is not more than 10%. The algorithm was used in the development of a new ultrasonic hardware-software complex, autonomous meteorological complex AMK-03-4, which differs from similar measuring instruments of turbulent atmosphere parameters by the presence of four identical ultrasonic anemometers. The design of the complex allows not only registration of the characteristics of turbulence, but also measurement of the statistical characteristics of the spatial derivatives of turbulent temperature fluctuations and orthogonal components of wind speed along each of the axes of the Cartesian coordinate system. This makes it possible to investigate the space–time structure of turbulent meteorological fields of the surface layer of the atmosphere for subsequent applications in the Monin–Obukhov similarity theory and to study turbulent coherent structures. The new measurement data of the spatial derivatives of temperature at stable stratification (at positive Monin–Obukhov parameters) were obtained, at which the behavior of the derivatives was been investigated earlier. In the most part of the interval of positive Monin–Obukhov parameters, the vertical derivative of the temperature is close to a constant value. This fact can be considered as a new significant result in similarity theory.

Description: In this study, the atmospheric moisture transport involved in the East Asian summer monsoon (EASM) water cycle is examined. Observational estimates are contrasted with the Met Office Unified Model (MetUM) climate simulations to evaluate the model’s ability to capture this transport. We explore the role of large circulation in determining the regional water cycle by analyzing key systematic errors in the model. MetUM exhibits robust errors in its representation of the summer Asian-Pacific monsoon system, including dry biases in the Indian peninsula and wet biases in the tropical Indian Ocean and tropical West Pacific. Such errors are consistent with errors in the atmospheric moisture convergence in the area. Diabatic heating biases in the Maritime Continent domain are shown, via nudging sensitivity experiments, to play a crucial role in remotely forcing the model circulation and moisture transport errors in the East Asian area. We also examine changes in the regional water cycle in response to interannual variability of the West Pacific subtropical high (WPSH). It is shown by water budget analysis that, although the model in general is not able to faithfully reproduce the response on a month to month basis, it gives comparable seasonal trends in regional moisture convergence and precipitation associated with shifts of the WPSH.

Description: The present paper aims to elucidate impacts of climate change on the availability and security of water and energy in the Middle East and North Africa region (MENA region; including the Eastern Mediterranean) in the context of the water–energy nexus. It largely builds on existing knowledge and understanding and aims to present a review of existing information on this topic. The region is particularly challenged by a number of factors, including the large variability of bio-geographical characteristics, extreme population growth over the last few decades, and substantial societal and economical transitions, as well as armed conflicts in some of the countries in the region. Anticipated changes in climate conditions will exacerbate the challenges regarding water and energy security in the region. Major impacts of climate change include a significant increase in summer temperatures, which will lead to a growing number of heat waves, primarily in urban structures. A general decrease in precipitation in many of the MENA countries is foreseen, resulting in enhanced droughts and a growing number of dry spells. In addressing energy and water scarcities and their mutual interrelationships, an integrated water–energy nexus concept offers promising prospects to improve environmental, climate, human, and political security. However, only very few countries in the MENA region have presently implemented such a concept. Mitigation and adaptation strategies addressing water and energy scarcity include enhanced efficiency of resource use, integrated technology assessments regarding electricity generation, and a stronger reliance on renewable/solar technologies. While looking at the MENA region as a whole, some emphasis will be given to Cyprus and the Eastern Mediterranean.

Description: Transparent exopolymer particles (TEP) act as a major transport mechanism for organic matter (OM) to the sea surface microlayer (SML) via bubble scavenging, and into the atmosphere via bubble bursting. However; little is known about the effects of bubble scavenging on TEP enrichment in the SML. This study examined the effects of several bubbling conditions and algae species on the enrichment of TEP in the SML. TEP enrichment in the SML was enhanced by bubbling, with a larger impact from bubbling rate than bubble size and increasing enrichment over time. Depth profiles showed that any TEP aggregates formed in the underlying water (ULW) were rapidly (〈2 min) transported to the SML, and that TEP was entrained in the SML by bubbling. Species experiments determined that the presence of different phytoplankton species and their subsequent release of precursor material further enhance the effectiveness of TEP enrichment via bubble scavenging.

Description: A better knowledge of the local and regional sources of the atmospheric particulate matter provides policy makers with the proper awareness when acting to improve air quality, in order to protect public health. A source apportionment study of the carbonaceous aerosol in Naples (Italy) is presented here, in order to improve this understanding in a vulnerable urban area. The aim of this study is quantifying directly fossil and non-fossil contributions to carbonaceous aerosol, by means of radiocarbon measurements. This is the first time that such an approach is implemented in this area. Fine particles with diameter ≤ 2.5 µm (PM2.5) were collected daily on top of a building in the city center, from November 2016 until January 2017. The carbonaceous aerosol was separated into organic carbon (OC) and elemental carbon (EC), by a two-step thermal desorption method. Subsequent radiocarbon analysis enabled the partitioning of the major sources of carbonaceous aerosol into fossil and non-fossil ones by applying radiocarbon isotopic mass balance. The PM2.5 concentration was on average 29 ± 3 µg⁄m3 (mean ± standard error; n = 18), with a maximum of 68.6 ± 0.7 µg⁄m3 on a day when air masses back-trajectories suggest a local origin and stagnant airflow conditions in the region. The carbonaceous component accounts for roughly half of the PM2.5 mass. Fossil fuel emissions are a minor source of OC (23%), but the dominant source of EC (66%), which is directly emitted during combustion processes. However, overall only 30% of the total carbon is of fossil origin, accounting for 14% of PM2.5 mass. Surprisingly, a comparable contribution is due to primary biomass burning carbon, which accounts in total for 15% of PM2.5 mass. Traffic pollution, the main cause of fossil fuel emissions in urban areas, is a significant, but not the predominant source of carbonaceous particle concentration. These findings support the conclusion of a predominant contribution from non-fossil sources to the carbon in airborne particulate matter, which policy makers should take into account when planning mitigation strategies to improve urban air quality.

Description: We present results of the studies of internal gravity waves based on altitude-time dependences of the temperature and the density of the neutral component and the velocity of the vertical plasma motion at altitudes of the lower ionosphere (60–130 km). The vertical plasma velocity, which in the specified altitude range is equal to the velocity of the neutral component, the temperature, and the density of the neutral atmosphere are determined by the method of the resonant scattering of radio waves by artificial periodic irregularities (APIs) of the ionosphere plasma. We have developed an API technique and now we are evolving it for studying the ionosphere and the neutral atmosphere using the Sura heating facility (56.1 N; 46.1 E), Nizhny Novgorod, Russia. An advantage of the API technique is the opportunity to determine the parameters of the undisturbed natural environment under a disturbance of the ionosphere by a field of powerful high frequency radio waves. Analysis of altitude-time variations of the neutral temperature, the density, and the vertical plasma velocity allows one to estimate periods of atmospheric waves propagation. Wavelike variations with a period from 5 min to 3 h and more are clearly determined.

Description: Desert dust is now recognized as a major health hazard. However, there still exists a lack of measurements of desert dust atmospheric loads in regions located in the vicinity of the major desert areas, where a growing part of the world population is living. Dust emission is a sporadic and intense phenomenon so that the classical air quality standards should not be appropriate to reflect the real population exposure to desert dust. In order to give some insight to answer this question, PM10 concentrations were continuously measured at a five-minute time step in southern Tunisia from February 2014 to February 2019. The daily and annual PM10 concentrations were first discussed according to the Tunisian air quality standards: In this relatively remote area, close to dust source, these standards were respected at the annual, but no longer at the daily scale. Measurements performed at a high temporal resolution (five minutes) allowed to discriminate the different situations that led to exceed Tunisian daily standards in southern Tunisia. In particular, intense (five-minute PM10 concentrations up to more than 1500 µg m−3) and short-lived (a few hours) dust outbreaks were regularly observed. This result raises the question of the necessity of defining specific air quality standards at the sub-daily scale in countries affected by intense desert dust outbreaks.

Description: Based on the ensemble empirical mode decomposition method, this study explores the differences and similarities in multiple time-scale characteristics of summer air temperature (T) and equivalent temperature (Te) over China during 1961–2017, using daily meteorological observations collected at 412 stations in China. Their relationships to global sea surface temperature variations is also discussed. Results show that both T and Te can be decomposed into five components, which includes multiple timescales, from interannual to long-term trends. The spatial patterns of each timescale’s leading mode show that the variations of Te are generally larger than that of T. Meanwhile, both T and Te are dominated by their inter-annual, multi-decadal variations and the non-linear trend. High correlations of T and Te can also be found in these major scales. The related sea surface temperature variations in these major scales also show consistent patterns, which correspond to El Niño–Southern Oscillation, Atlantic Multidecadal Oscillation and the global warming trend in the sea, respectively. In other scales, both spatial patterns of T and Te and the corresponding correlation patterns with sea surface temperature are distinguishable. The current results explore the compound changes of surface temperature-humidity during the past five decades from a new perspective, which provides some insights for a better understanding of the possible causes of climate change over China.

Description: Polar mesocyclones (MCs) are small marine atmospheric vortices. The class of intense MCs, called polar lows, are accompanied by extremely strong surface winds and heat fluxes and thus largely influencing deep ocean water formation in the polar regions. Accurate detection of polar mesocyclones in high-resolution satellite data, while challenging, is a time-consuming task, when performed manually. Existing algorithms for the automatic detection of polar mesocyclones are based on the conventional analysis of patterns of cloudiness and they involve different empirically defined thresholds of geophysical variables. As a result, various detection methods typically reveal very different results when applied to a single dataset. We develop a conceptually novel approach for the detection of MCs based on the use of deep convolutional neural networks (DCNNs). As a first step, we demonstrate that DCNN model is capable of performing binary classification of 500 × 500 km patches of satellite images regarding MC patterns presence in it. The training dataset is based on the reference database of MCs manually tracked in the Southern Hemisphere from satellite mosaics. We use a subset of this database with MC diameters falling in the range of 200–400 km. This dataset is further used for testing several different DCNN setups, specifically, DCNN built “from scratch”, DCNN based on VGG16 pre-trained weights also engaging the Transfer Learning technique, and DCNN based on VGG16 with Fine Tuning technique. Each of these networks is further applied to both infrared (IR) and a combination of infrared and water vapor (IR + WV) satellite imagery. The best skills (97% in terms of the binary classification accuracy score) is achieved with the model that averages the estimates of the ensemble of different DCNNs. The algorithm can be further extended to the automatic identification and tracking numerical scheme and applied to other atmospheric phenomena that are characterized by a distinct signature in satellite imagery.

Description: Rainfall during typhoons is one of the most important water resources in Taiwan, but heavy typhoon rainfall often leads to serious disasters and consequently results in loss of lives and property. Hence, accurate forecasts of typhoon rainfall are always required as important information for water resources management and rainfall-induced disaster warning system. In this study, a methodology is proposed for providing quantitative forecasts of 24 h cumulative rainfall during typhoons. Firstly, ensemble forecasts of typhoon rainfall are obtained from an ensemble numerical weather prediction (NWP) system. Then, an evolutionary algorithm, i.e., genetic algorithm (GA), is adopted to real-time decide the weights for optimally combining these ensemble forecasts. That is, the novelty of this proposed methodology is the effective integration of the NWP-based ensemble forecasts through an evolutionary algorithm-based strategy. An actual application is conducted to verify the forecasts resulting from the proposed methodology, namely NWP-based ensemble forecasts with a GA-based integration strategy. The results confirm that the forecasts from the proposed methodology are in good agreement with observations. Besides, the results from the GA-based strategy are more accurate as compared to those by simply averaging all ensemble forecasts. On average, the root mean square error decreases about 7%. In conclusion, more accurate typhoon rainfall forecasts are obtained by the proposed methodology, and they are expected to be useful for disaster warning system and water resources management during typhoons.

Description: We analyse and quantify the recurrences of European temperature extremes using 32 historical simulations (1900–1999) of the fifth Coupled Model Intercomparison Project (CMIP5) and 8 historical simulations (1971–2005) from the EUROCORDEX experiment. We compare the former simulations to the 20th Century Reanalysis (20CRv2c) dataset to compute recurrence spectra of temperature in Europe. We find that, (1) the spectra obtained by the model ensemble mean are generally consistent with those of 20CR; (2) spectra biases have a strong regional dependence; (3) the resolution does not change the order of magnitude of spectral biases between models and reanalysis, (4) the spread in recurrence biases is larger for cold extremes. Our analysis of biases provides a new way of selecting a subset of the CMIP5 ensemble to obtain an optimal estimate of temperature recurrences for a range of time-scales.

Description: This study evaluates the predictability of the formation phase of a radiation fog event observed during the night of 31 October 2015 to 01 November 2015 in the north-east of France at three sites managed by OPE (Observatoire Pérenne de l’Environnement). The fog layer shows significantly different behaviors at the three areas, which are located only a few kilometers apart. Three fog life cycles were observed: the formation of a dense adiabatic fog, the formation of a thin patchy fog, or no fog formation despite favorable conditions. This event was studied with the Meso-NH numerical mesoscale model at two horizontal resolutions, 500 m and 50 m. Simulations at 50 m allow estimation of the spread of the predicted parameters over the heterogeneous terrain studied. These numerical simulations strongly suggest that this event involved numerous interactions and complex circulations. The wind above the nocturnal boundary layer greatly affects the transition of shallow patchy fog into thick adiabatic fog. These numerical simulations also show that the occurrence and type of fog could be very different over a small but heterogeneous area. It is also interesting to note that the spread of the simulated parameters was very high during the transition from shallow fog to a deep fog layer. The spread was concentrated during the regime transition between the fog formation and its maturity. This appeared to be the result of the complex interplay of processes at numerous ranges of scale. A new concept called “pseudo-process diagram” is presented. These pseudo-process diagrams are very good tools to analyze fog, and allow a good illustration of the spread of fog during this chaotic phase. This kind of concept seems a promising tool to analyze fog predictability in depth.

Description: In this research, for the first time, we present the evaluation of a semi-continuous pit manure recharge system on the mitigation of ammonia (NH3) and hydrogen sulfide (H2S) emissions from a swine finisher barn. The pit recharge system is practiced on many swine farms in the Republic of Korea, primarily for improving air quality in the barn. It consists of an integrated waste management system where the fraction of stored manure is pumped out (10× of the daily production of manure, 3× a day); solids are separated and composted, while the aerobically treated liquid fraction is then returned to the pit. We compared emissions from two 240-pig rooms, one equipped with a pit recharge system, and the other operating a conventional slurry pit under the slatted floor. Mean reduction of NH3 and H2S emissions were 49 ± 6% and 82 ± 7%, respectively, over 14 days of measurements. The removal efficiency of H2S was higher than NH3, likely because the pH of aerobically treated liquid manure remained slightly above 8. More work is warranted to assess the N balance in this system and the emissions of odor and greenhouse gasses (GHGs). It is also expected that it will be possible to control the NH3 and H2S removal rates by controlling the nitrification level of the liquid manure in the aerobic treatment system.

Description: We present observations of X-rays from laboratory sparks created in the air at atmospheric pressure by applying an impulse voltage with long (250 µs) rise-time. X-ray production in 35 and 46 cm gaps for three different electrode configurations was studied. The results demonstrate, for the first time, the production of X-rays in gaps subjected to switching impulses. The low rate of rise of the voltage in switching impulses does not significantly reduce the production of X-rays. Additionally, the timing of the X-ray occurrence suggests the possibility that the mechanism of X-ray production by sparks is related to the collision of streamers of opposite polarity.

Description: Aerosol mixing state significantly affects concentrations of cloud condensation nuclei (CCN), wet removal rates, thermodynamic properties, heterogeneous chemistry, and aerosol optical properties, with implications for human health and climate. Over the last two decades, significant research effort has gone into finding computationally-efficient methods for representing the most important aspects of aerosol mixing state in air pollution, weather prediction, and climate models. In this review, we summarize the interactions between mixing-state and aerosol hygroscopicity, optical properties, equilibrium thermodynamics and heterogeneous chemistry. We focus on the effects of simplified assumptions of aerosol mixing state on CCN concentrations, wet deposition, and aerosol absorption. We also summarize previous approaches for representing aerosol mixing state in atmospheric models, and we make recommendations regarding the representation of aerosol mixing state in future modelling studies.

Description: This paper reports estimated maintenance-cleaning costs, cost savings and cleaning interval increases for structural surfaces and windows in Europe obtainable by reducing the air pollution. Methodology and data from the ICP-materials project were used. The average present (2018) cleaning costs for sheltered white painted steel surfaces and modern glass due to air pollution over background, was estimated to be ~2.5 Euro/m2∙year. Hypothetical 50% reduction in the air pollution was found to give savings in these cleaning costs of ~1.5 Euro/m2∙year. Observed reduction in the air pollution, from 2002–2005 until 2011–2014, have probably increased the cleaning interval for white painted steel with ~100% (from 12 to 24 years), representing reductions in the single intervention cleaning costs from 7 to 4%/year (= % of one cleaning investment, per year during the cleaning interval) and for the modern glass with ~65% (from 0.85 to 1.3 years), representing reductions in the cleaning cost from 124 to 95%/year. The cleaning cost reductions, obtainable by 50% reduction in air pollution, would have been ~3 %/year for white painted steel and ~60%/year for the modern glass, representing ~100 and 50% additional cleaning interval increases. These potential cleaning cost savings are significantly higher than previously reported for the weathering of Portland limestone ornament and zinc monuments.

Description: Crop residue burning negatively impacts both the environment and human health, whether in the aspect of air pollution, regional and global climate change, or transboundary air pollution. Accordingly, this study aims to assess the level of air pollutant emissions caused by the rice residue open burning activities in 2018, by analyzing the remote sensing information and country specific data. This research also aims to analyze the trend of particulate matter 10 microns or less in diameter (PM10) concentration air quality sites in provinces with large paddy rice planting areas from 2010–2017. According to the results, 61.87 megaton (Mt) of rice residue were generated, comprising 21.35 Mt generated from the irrigated fields and 40.53 Mt generated from the rain-fed field. Only 23.0% of the total rice residue generated were subject to open burning—of which nearly 32% were actually burned in the fields. The emissions from such rice residue burning consisted of: 5.34 ± 2.33 megaton (Mt) of CO2, 44 ± 14 kiloton (kt) of CH4, 422 ± 179 kt of CO, 2 ± 2 kt of NOX, 2 ± 2 kt of SO2, 38 ± 22 kt of PM2.5, 43 ± 29 kt of PM10, 2 ± 1 kt of black carbon (BC), and 14 ± 5 kt of organic carbon (OC). According to the air quality trends, the results shows the higher level of PM10 concentration was due to the agricultural burning activities, as reflected in the higher monthly averages of the months with the agricultural burning, by around 1.9–2.1 times. The result also shows the effect of government’s policy for farmers on the crop burning activities and air quality trends.

Description: European beech (Fagus sylvatica L.) is the most important deciduous tree species in Europe. According to different climate scenarios, there is a relatively high probability of a massive decline in and loss of beech forests in southern Europe and in the southern part of central Europe. Thus, the authors of this study explored the dynamics of tree diameter increments and the influence of extremely dry years on the width of tree rings. This study used dendroecological methods to analyze the growth and diameter increments of European beech trees at locations in Serbia and the Republic of Srpska. The sampling was conducted along the vertical distribution of beech forests, at five sites at the lower limit of the distribution, at five optimal sites of the distribution, and at five sites at the upper limit of the distribution. Long-term analyses indicate that dry conditions during a growing season can reduce tree-ring width, but a reduction in tree growth can be expected as a result of more than one season of unfavorable conditions. Low temperatures in autumn and winter and prolonged winters can strongly affect upcoming vegetation and reduce tree development even under normal thermal conditions during a growing season.

Description: Several simulations of the surface climate and energy balance of Vatnajökull ice cap, Iceland, are used to estimate the glacier runoff for the period 1980–2015 and the sensitivity of runoff to the spring conditions (e.g., snow thickness). The simulations are calculated using the snow pack scheme from the regional climate model HIRHAM5, forced with incoming mass and energy fluxes from the numerical weather prediction model HARMONIE-AROME. The modeled runoff is compared to available observations from two outlet glaciers to assess the quality of the simulations. To test the sensitivity of the runoff to spring conditions, simulations are repeated for the spring conditions of each of the years 1980–2015, followed by the weather of all summers in the same period. We find that for the whole ice cap, the variability in runoff as a function of varying spring conditions was on average 31% of the variability due to changing summer weather. However, some outlet glaciers are very sensitive to the amount of snow in the spring, as e.g., the variation in runoff from Brúarjökull due to changing spring conditions was on average 50% of the variability due to varying summer weather.

Description: More than half of the people in the world now live in cities and this proportion is expected to continue to increase. [...]

Description: A prequel study showed that dynamic downscaling using a regional climate model (RCM) over Africa improved the Goddard Institute for Space Studies Atmosphere-Ocean Global Climate Model (GISS AOGCM: ModelE) simulation of June–September rainfall patterns over Africa. The current study applies bias corrections to the lateral and lower boundary data from the AOGCM driving the RCM, based on the comparison of a 30-year simulation to the actual climate. The analysis examines the horizontal pattern of June–September total accumulated precipitation, the time versus latitude evolution of zonal mean West Africa (WA) precipitation (showing monsoon onset timing), and the latitude versus altitude cross-section of zonal winds over WA (showing the African Easterly Jet and the Tropical Easterly Jet). The study shows that correcting for excessively warm AOGCM Atlantic sea-surface temperatures (SSTs) improves the simulation of key features, whereas applying 30-year mean bias corrections to atmospheric variables driving the RCM at the lateral boundaries does not improve the RCM simulations. We suggest that AOGCM climate projections for Africa should benefit from downscaling by nesting an RCM that has demonstrated skill in simulating African climate, driven with bias-corrected SST.

Description: The magnitude and timing of seasonal rainfall is vitally important to the health and vitality of key agro-ecological and social-economic systems of the Niger River Basin. Given this unique context, knowledge concerning how climate change is likely to impact future rainfall characteristics and patterns is critically needed for adaptation and mitigation planning. Using nine ensemble bias-corrected climate model projection results under RCP4.5 and RCP8.5 (RCP—Representative Concentration Pathway) emissions scenarios at the mid-future time period, 2021/2025-2050 from the Coordinated Regional Climate Downscaling Experiments (CORDEX) dataset; this study provides a comprehensive analysis of the projected changes in rainfall characteristics in three agro-ecological zones of the Niger River Basin. The results show an increase in the average rainfall of about 5%, 10–20% and 10–15% for the Southern Guinea, Northern Guinea and Sahelian zones, respectively, relative to the baseline, 1981/1985–2005. On the other hand, the change in future rainfall intensities are largely significant and the frequency of rainfall at the low, heavy and extreme rainfall events in the future decrease at most locations in the Niger River Basin. The results also showed an increase in the frequency of moderate rainfall events at all locations in the basin. However, in the Northern Guinea and Sahel locations, there is an increase in the frequency of projected heavy and extreme rainfall events. The results reveal a shift in the future onset/cessation and a shortening of the duration of the rainy season in the basin. Specifically, the mean date of rainfall onset will be delayed by between 10 and 32 days. The mean onset of cessation will also be delayed by between 10 and 21 days. It is posited that the projected rainfall changes pose serious risks for food security of the region and may require changes in the cropping patterns and management.

Description: Micro-scale cold-air flow along a gentle slope was analyzed using thermal infrared imaging (TIR), focusing exclusively on the lowermost 2 m above ground. Cold-air pulses were analyzed with regard to their vertical temperature stratification as well as flow characteristics, such as flow speed. Analyses on the transition zone between the near-surface very stable inversion layer and the less stable, warmer air above highlight turbulent situations and detrainment effects at the cold-air inversion top. Using thermal imaging in a high spatiotemporal resolution with up to 90 vertical data points and TIR pixels for 1.5 m cold-air depth, a high-precision cold-air flow analysis was realized.

Description: The rain statistics of 0–45° N area including equatorial, Sahelian, and mid-latitude regions, are studied using the probability distributions of the duration of rainy and dry events. Long time daily data set from ground measurements and satellite observations of rain fields are used. This technique highlights a sharp latitudinal transition of the statistics between equatorial and all other regions (Sahel, mid-latitude). The probability distribution of the 8° S to 8° N latitude band shows a large-scale organization with a slow decreasing (power law decrease) distributions for the time and space size of rain events. This observation is in agreement with a scaling, or macro turbulent, behavior of the equatorial regions rain fields. For the Sahelian and mid-latitude regions, our observations are clearly not in agreement with this behavior. They show that the largest rain systems have a limited time and space size (well described with a decreasing exponential distribution). For these non-equatorial regions it is possible to define a local characteristic duration and a characteristic horizontal size of the large rain events. These characteristics time and space scales of observed mesoscale convective systems could be a sensible indicator for the detection of the possible trend of rain distribution properties due to anthropogenic influence.

Description: Managed grassland is occasionally renovated to maintain plant productivity by killing old vegetation, ploughing, and reseeding. This study aimed to investigate the combined effect of grassland renovation and long-term manure application on the temporal dynamics of nitrous oxide (N2O) emission and nitrate nitrogen (NO3−–N) leaching. The study was conducted from September 2013 to September 2016 in a managed grassland renovated in September 2013. In this grassland, two treatments were managed—chemical fertilizer application (F treatment) and the combined application of chemical fertilizer and beef cattle manure (MF treatment)—for eight years before the renovation. The control treatment without fertilization (CT) was newly established in the F treatment. The soil N2O flux was measured using a closed chamber method. A leachate sample was collected using a tension-free lysimeter that was installed at the bottom of the Ap horizon (25 cm deep), and total NO3−–N leaching was calculated from leachate NO3−–N concentration and drainage volume was estimated by the water balance method. In the first year after renovation, the absence of plant nitrogen uptake triggered NO3−–N leaching following rainfall during renovation and increased drainage water after thawing. NO3−–N movement from topsoil to deeper soil enhanced N2O production and emission from the soil. N2O emission in MF treatment was 1.6–2.0 times larger than those of CT and F treatments, and NO3−–N leaching in MF treatment was 2.3–2.6 times larger than those of CT and F treatments in the first year. Mineral nitrogen release derived from long-term manure application increased NO3−–N leaching and N2O emission. In the second year, N2O emission and NO3−–N leaching significantly decreased from the first year because of increased plant N uptake and decreased mineral nitrogen surplus, and no significant differences in N2O emission and NO3−–N leaching were observed among the treatments. In the second and third years, NO3−–N leaching was regulated by plant nitrogen uptake. There were no significant differences in NO3−–N leaching among the treatments, but N2O emission in MF treatment was significantly smaller than in the F treatment. Long-term manure application could be a possible option to mitigate N2O emission in permanent grassland; however, the risk of increased NO3−–N leaching and N2O emission in the renovation year induced by manure nitrogen release should be noted.

Description: This study explored the influence of choosing a nonhydrostatic dynamical core or a hydrostatic dynamical core in the weather research and forecasting (WRF) model on the intensity and structure of simulated tropical cyclones (TCs). A comparison of cloud-resolving simulations using each core revealed significant differences in the TC simulations. In comparison with the nonhydrostatic simulation, the hydrostatic simulation produced a stronger and larger TC, associated with stronger convective activity. A budget analysis of the vertical momentum equation was conducted to investigate the underlying mechanisms. Although the hydrostatic dynamical core was used, the vertical motion was not in strict hydrostatic balance because of the existence of the vertical perturbation pressure gradient force, local buoyancy force, water loading, and sum of the Coriolis and diffusion effects. The contribution of the enhanced vertical perturbation pressure gradient force was found to be more important for stronger upward acceleration in the eyewall in the hydrostatic simulation than in the nonhydrostatic simulation. This is because it leads to intensified convection in the eyewall that releases more latent heat, which induces a larger low-level radial pressure gradient and inflow motion, and eventually leads to a stronger storm.

Description: This paper describes experiments on the ageing of a monolayer model for the air–water interface of marine aerosols composed of a typical glycolipid, galactocerebroside (GCB). Lipopolysaccharides have been observed in marine aerosols, and GCB is used as a proxy for these more complex lipopolysaccharides. GCB monolayers are investigated as pure films, as mixed films with palmitic acid, which is abundant in marine aerosols and forms a stable attractively mixed film with GCB, particularly with divalent salts present in the subphase, and as mixed films with palmitoleic acid, an unsaturated analogue of palmitic acid. Such mixed films are more realistic models of atmospheric aerosols than simpler single-component systems. Neutron reflectometry (NR) has been combined in situ with Fourier transform infra-red reflection absorption spectroscopy (IRRAS) in a pioneering analysis and reaction setup designed by us specifically to study mixed organic monolayers at the air–water interface. The two techniques in combination allow for more sophisticated observation of multi-component monolayers than has previously been possible. The structure at the air–water interface was also investigated by complementary Brewster angle microscopy (BAM). This study looks specifically at the oxidation of the organic films by nitrate radicals (NO3•), the key atmospheric oxidant present at night. We conclude that NO3• oxidation cannot fully remove a cerebroside monolayer from the surface on atmospherically relevant timescales, leaving its saturated tail at the interface. This is true for pure and salt water subphases, as well as for single- and two-component films. The behaviour of the unsaturated tail section of the molecule is more variable and is affected by interactions with co-deposited species. Most surprisingly, we found that the presence of CaCl2 in the subphase extends the lifetime of the unsaturated tail substantially—a new explanation for longer residence times of materials in the atmosphere compared to lifetimes based on laboratory studies of simplified model systems. It is thus likely that aerosols produced from the sea-surface microlayer at night will remain covered in surfactant molecules on atmospherically relevant timescales with impact on the droplet’s surface tension and on the transport of chemical species across the air–water interface.

Description: In this study, we investigated the effects of grid and spectral nudging in regional hydroclimatic simulations over the Eastern Mediterranean climate change hot-spot. We performed year-long simulations for the hydrological year October 2001–September 2002 using the Weather Research and Forecasting (WRF) model at 12-km resolution, driven by the ERA-Interim reanalyses. Six grid and three spectral nudging options were tested using a number of model configurations. Due to the large uncertainty of regional observations, we compared the model with various satellite- and station-based meteorological datasets. The effect of nudging was tested for mean weather conditions and precipitation characteristics and extremes. For certain parts of the study domain, WRF was found to reproduce both aspects of rainfall over the Eastern Mediterranean reasonably well. Our findings highlighted that, for the WRF modeling system, nudging is critical for the simulation of rainfall; however, the application of interior constraint methods was found to have different impacts on various locations and climatic regimes. For the hyperarid parts of the domain, nudging did not improve the simulation of precipitation amounts (about 20% additional drying was introduced), while it added much value for the wetter rainfall regimes of the Eastern Mediterranean (corrections of about 30%). Improvements in the simulated precipitation were mostly introduced by spectral nudging; however, this option required significant computational resources. For these ERA-Interim-driven simulations, grid nudging that involves specific humidity within the planetary boundary layer is not recommended for the simulation of precipitation since it introduces dry biases up to 75–80%.

Description: Extreme warm stratospheric events during polar winters from ERA-Interim reanalysis and CMIP5-ESM-LR runs were separated by duration and strength of the polar-night jet oscillation (PJO) using a high statistical confidence level of three standard deviations (strong-PJO events). With a composite analysis, we demonstrate that strong-PJO events show a significantly stronger downward propagating signal in both, northern annular mode (NAM) and zonal mean zonal wind anomaly in the stratosphere in comparison with non-PJO events. The lower stratospheric EP-flux-divergence difference in ERA-Interim was stronger in comparison to long-term CMIP5-ESM-LR runs (by a factor of four). This suggests that stratosphere–troposphere coupling is stronger in ERA-Interim than in CMIP5-ESM-LR. During the 60 days following the central date (CD), the Arctic oscillation signal was more intense during strong-PJO events than during non-PJO events in ERA-Interim data in comparison to CMIP5-ESM-LR runs. During the 15-day phase after CD, strong PJO events had a significant increase in stratospheric ozone, upper tropospheric zonally asymmetric impact, and a regional surface impact in ERA-Interim. Finally, we conclude that the applied high statistical threshold gives a clearer separation of extreme warm stratospheric events into strong-PJO events and non-PJO events including their different downward propagating NAM signal and tropospheric impacts.

Description: A meridional Northern Hemisphere (NH) circulation epoch, which began in 1957, is marked by changes in the temperature and precipitation regimes over southwest Russia and central USA depending on the occurrence of NH atmospheric circulation regimes. A classification scheme proposed in 1968, and studied later put forth 13 NH circulation types, fitting more broadly into four groups, two of which are more zonal type flows and two of which are more meridional flows. Using the results of a previous study that showed four distinct sub-periods during the 1957–2017 epoch, the temperature and precipitation regimes of both regions were studied across all seasons in order to characterize modern day climate variability and their suitability for vegetation growth. Then the Hydrologic Coefficient, which combined the temperature and precipitation variables, was briefly studied. The most optimal conditions for vegetation growth, positive temperature and precipitation anomalies, were noted during the period 1970–1980 for southwest Russia, which was dominated by an increasingly more zonal flow regime in the Belgorod region and NH in general. For the central USA, the HTC showed more ideal conditions for agriculture in recent years due to favorable precipitation occurrence. In southwest Russia, variable precipitation regimes were noted during the meridional flow periods, and with the increase in temperature (since 1998), these can adversely affect the hydrothermal characteristics of the growing season. Finally, a comparison of the 13 NH circulation types with several teleconnection indexes demonstrated the robustness of the NH flow regime classification scheme used here.

Description: This work documents the main features of six subtropical cyclones occurred between the years 2010 and 2016 over the southwestern South Atlantic Ocean, near the Brazilian coast, which received names (with the exception of one) from the Brazilian Navy Hydrographic Center. The fine-resolution ERA5 reanalysis and rainfall estimates from the Tropical Rainfall Measuring Mission (TRMM) were used to describe the synoptic environment and the adverse weather conditions during the six events. The support of a small-amplitude trough at mid-levels or a cut-off low, weak vertical wind shear, and moisture flux convergence are the main features contributing to the subtropical cyclogenesis at the surface. On the other hand, sea surface temperature (SST) presents a secondary contribution since the cyclones develop over the ocean with a wide range of SST values (from 22.5 °C to 28.6 °C in the initial phase of cyclones). The six subtropical cyclones are less deep in the atmosphere column than the tropical ones and, unlike the extratropical cyclones, they have little or no westward tilt with an increase in height. The studied subtropical cyclones produced adverse weather conditions such as (a) strong winds (reaching 17 m·s−1 at 10 m high) for a long period occurring east/southeastward of the cyclone center, and (b) high amounts of rainfall along the southeastern coast of Brazil, where the accumulated rainfall varied between 170 to 350 mm, being in most cases higher than the monthly climatology. Over the continent, the Brazilian states of Rio de Janeiro and Espírito Santo were the most affected by the intense rainfall associated with the cyclones.

Description: Sea spray aerosol (SSA) is highly enriched in marine-derived organic compounds during seasons of high biological productivity, and saturated fatty acids comprise one of the most abundant classes of molecules. Fatty acids and other organic compounds form a film on SSA surfaces, and SSA particle surface-area-to-volume ratios are altered during aging in the marine boundary layer (MBL). To understand SSA surface organization and its role during dynamic atmospheric conditions, an SSA proxy fatty acid film and its individual components stearic acid (SA), palmitic acid (PA), and myristic acid (MA) are studied separately using surface pressure–area ( Π − A ) isotherms and Brewster angle microscopy (BAM). The films were spread on an aqueous NaCl subphase at pH 8.2, 5.6, and 2.0 to mimic nascent to aged SSA aqueous core composition in the MBL, respectively. We show that the individual fatty acid behavior differs from that of the SSA proxy film, and at nascent SSA pH the mixture yields a monolayer with intermediate rigidity that folds upon film compression to the collapse state. Acidification causes the SSA proxy film to become more rigid and form 3D nuclei. Our results reveal film morphology alterations, which are related to SSA reflectivity, throughout various stages of SSA aging and provide a better understanding of SSA impacts on climate.

Description: The current wealth of spaceborne passive and active measurements from ultraviolet to the infrared wavelengths provides an unprecedented opportunity to construct ice cloud bulk optical property models that lead to consistent ice cloud property retrievals across multiple sensors and platforms. To infer the microphysical and radiative properties of ice clouds from these satellite measurements, the general approach is to assume an ice cloud optical property model that implicitly assumes the habit (shape) and size distributions of the ice particles in these clouds. The assumption is that this ice optical property model will be adequate for global retrievals. In this review paper, we first summarize the key optical properties of individual particles and then the bulk radiative properties of their ensemble, followed by a review of the ice cloud models developed for application to satellite remote sensing. We illustrate that the random orientation condition assumed for ice particles is arguably justified for passive remote sensing applications based on radiometric measurements. The focus of the present discussion is on the ice models used by the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Clouds and Earth’s Radiant Energy System (CERES) science teams. In addition, we briefly review the ice cloud models adopted by the Polarization and Directionality of the Earth’s Reflectance (POLDER) and the Himawari-8 Advanced Himawari Imager (AHI) for ice cloud retrievals. We find that both the MODIS Collection 6 ice model and the CERES two-habit model result in spectrally consistent retrievals.

Description: Observational data from the Global Precipitation Measurement (GPM) Core Observatory during four summers (2014–2017) has been used to investigate deep convection systems (DCSs) over the Tibetan Plateau (TP) and its south slope (SS). The frequency, geographical distribution diurnal variation, and vertical structure of DCSs over the TP and SS are compared among these two regions. The frequency of DCSs over the SS (0.98%) was far higher than over the TP (0.15%), suggesting that stronger DCSs occur to the east and south of the TP. The maximum number of DCS occurred in July and August. A clear diurnal variation in DCS was found over the whole region, DCSs over the TP and SS both have a greatest amplitude in the afternoon. The probability of DCSs from 1200 to 1800 local time (LT) was 76.3% and 44.1% over TP and SS respectively, whereas the probability of DCSs being generated from 2200 (LT) to 0600 on the next day LT was 0.03% and 33.1% over the TP and SS respectively. There was a very low frequency of DCSs over the TP during the night. Five special echo top heights were used to investigate the vertical structure of DCSs. DCSs over the TP were both weaker and smaller than those over the SS.

Description: In this work, we have characterized the iron local structure in samples of two different types of atmospheric dust using X-ray absorption spectroscopy and selective leaching experiments. Specifically, we have investigated samples of long-range transported Saharan dust and freshly emitted steel plant fumes with the aim of individuating possible fingerprints of iron in the two cases. Findings include (1) prevalence of octahedral coordinated Fe 3 + for all samples; (2) presence of 6-fold coordinated Fe 3 + , aluminosilicates and iron oxy(hydr)oxides in Saharan dust and (3) of Fe-bearing spinel-like structures in the industrial fumes; (4) general predominance of the residual insoluble fraction with a notable difference: 69% for Saharan dust and 93% for steel production emissions, associated with aluminosilicates and non-reducible iron oxy(hydr)oxides, and Fe spinels, respectively. The remarkable differences between the X-ray absorption spectroscopy (XAS) spectra and leaching test results for the two sample types suggest the possibility to exploit the present approach in more complex cases. To this aim, two additional case studies of mixed aerosol samples are presented and discussed.

Description: Large-Eddy Simulations (LES) corresponding to four convective intensive observation periods of Sagebrush Phase 1 tracer experiment were conducted with realistic boundary conditions using Weather Research and Forecast model (WRF). Multiple nested domains were used to dynamically downscale the conditions from domain with grid size of 24 km to local scales with grid size of 150 m. Sensitivity analysis of mesoscale model was conducted using three boundary layer, three surface layer and two micro-physics schemes. Model performance was evaluated by comparing the surface meteorological variables and boundary layer height from the mesoscale runs and observed values during tracer experiment. Output from mesoscale simulations was used to drive the LES domains. Effect of vertical resolution and sub-grid scale parameterizations were studied by comparing the wind speed and direction profiles along with turbulent kinetic energy at two different heights. Atmospheric stability estimated using the Richardson number and shear exponent evaluated between 8- and 60-m levels was found to vary between weakly unstable to unstable. Comparing the wind direction standard deviations coupled with the wind speeds showed that the WRF-LES underestimated the wind direction fluctuations for wind speeds smaller than 3-ms − 1 . Based on the strengths of convection and shear, WRF-LES was able to simulate horizontal convection roll and convective cell type features.

Description: The variability and mechanisms of multi-decadal megadroughts over eastern China during the last millennium were investigated using a control, full-forcing, and four sensitivity experiments from the Community Earth System Model (CESM) Last Millennium Ensemble (LME) archive. The model simulated megadroughts have comparable magnitudes and durations with those derived from reconstructed proxy data, although the megadroughts are not temporally synchronous. In all experiments, the megadroughts exhibit similar spatial structures, corresponding to a weakening of the East Asia summer monsoon (EASM) and a strengthening of the East Asia winter monsoon (EAWM). The results show that internal climate variability within the coupled climate system plays an essential role in triggering megadroughts, while different external forcings may contribute to persistence and modify the anomaly patterns of megadroughts. A pattern of meridional tripolar (warm-cold-warm) sea surface temperature (SST) anomalies in the western Pacific stretching from the equator to high latitude is responsible for the EASM weakening and EAWM strengthening. The weakening of the EASM and strengthening of the EAWM are essentially caused by negative SST anomalies over the northwestern Pacific and positive SST anomalies over the equatorial western Pacific, which are associated with a La Niña-like SST gradient across the tropical Pacific. The external forcings prolong the megadroughts through maintenance of the meridional tripolar SST anomalies and enlarge the megadrought spatial extent by magnifying the meridional tripolar SST anomalies.

Description: The establishment of sown pasture is an important agricultural practice in many landscapes. Although both native grassland and sown pasture play a key role in the global carbon cycle, due to lack of data and field experiments, our understanding of grassland CH4 fluxes and CO2 emissions remains limited, especially when it comes to sown pasture. We measured ecosystem respiration and CH4 fluxes in response to a variety of potential drivers (soil temperature, soil moisture, ammonium nitrogen, nitrate nitrogen and dissolved organic carbon) in CG (continuous grazing), RG (rotational grazing) and UG (ungrazed) plots in sown grassland for one year in Inner Mongolia. Fluxes of CH4 and ecosystem respiration were measured using static opaque chambers and gas chromatography. Grazing significantly reduced ecosystem respiration (p 〈 0.01), and grazing pattern significantly influenced respiration in CG and RG plots (p 〈 0.01). We find that the sown grassland is a net sink for atmospheric CH4. No influence of grazing pattern was observed on CH4 flux in CG, RG and UG (p 〉 0.05). Soil temperature is the most important factor influencing ecosystem respiration and CH4 flux in the sown grassland, with soil moisture playing a secondary role to soil temperature. Variation in levels of ammonium nitrogen, nitrate nitrogen and dissolved organic carbon had little influence on ecosystem respiration or CH4 flux (except in UG plots). The values obtained for ecosystem respiration of grasslands have a large uncertainty range, which may be due to spatial variability as well as differences in research methods. Mean CH4 fluxes measured only during the growing season were much higher than the annual mean CH4 fluxes.

Description: The range and time of the environmental effects of Asian dust are closely dependent on the pathways and the speed of dust plume movement. In this study, the occurrence and movement of two dust storms in China in May 2017 were examined by using open space- and ground-based measurement data and the backward trajectories of dust plumes. Results from the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis data showed that the dust storms were caused by the rapid coupling development of Mongolian cyclones and Asian highs. After the dust plumes arrived at the Southeastern China in the first dust event, the stable weather conditions and the Asian high slowed down the movement of the plumes, leading to the gradual diffusion of dust particles. Moreover, the Asian high in the first event and the Huabei low (a low-pressure system in North China Plain) in the second altered the movement direction of the dust plumes from southward to northward, which we denote as the “dust reverse transport (DRT)”. The DRT occurred only within the lower troposphere even though dust plumes could extended to 5–10 km in vertical direction. Statistical results of 28 spring dust events occurred in 2015–2018 showed that all these dust storms were triggered by Mongolian cyclones and/or Asian highs, and approximately 39% moved as the DRT, indicating about one third of severe spring dust storms could influence larger areas or longer time than the remained ones.

Description: Recently, NCAR (the National Center for Atmospheric Research) released the Community Earth System Model’s low-warming simulations, which provided long-term climate data for stabilization pathways at 1.5 °C and 2.0 °C above pre-industrial levels. Based on these data, six extreme low temperature indices—TXn (coldest day), TNn (coldest night), TX10p (cool days), TN10p (cool nights), CSDI (cold spell duration indicator), and DTR (diurnal temperature range)—were calculated to assess the changes in extreme low temperature over Northern China under 1.5 °C and 2.0 °C warmer future. The results indicate that compared to the preindustrial level, the whole of China will experience 0.32–0.46 °C higher minimum surface air temperature (SAT) warming than the global average, and the winter temperature increase in Northern China will be the most pronounced over the country. In almost all the regions of Northern China, especially Northeast and Northwest China, extreme low temperature events will occur with lower intensity, frequency, and duration. Compared with the present day, the intensity of low temperature events will decrease most in Northeast China, with TXn increasing by 1.9 °C/2.0 °C and TNn increasing by 2.0 °C/2.5 °C under 1.5 °C/2.0 °C global warming, respectively. The frequency of low temperature events will decrease relatively more in North China, with TX10p decreasing by 8 days/11 days and TN10p decreasing by 7 days/9 days under 1.5 °C/2.0 °C warming. CSDI will decrease most in Northwest China, with decreases of 7 days/10 days with 1.5 °C/2.0 °C warming. DTR will decrease in the Northwest and Northeast but increase in North China, with −0.9 °C/−2.0 °C in the Northwest, −0.4 °C/−1.5 °C in the Northeast, and 1.7 °C/2.0 °C in North China in the 1.5 °C/2.0 °C warming scenarios. For temperatures lower than the 5th percentile, the PRs (probability ratios) will be 0.68 and 0.55 of that of the present day under 1.5 °C and 2.0 °C warmer futures, respectively. Global warming of 2.0 °C instead of 1.5 °C will lead to extreme low temperature events decreasing by 6–56% in regard to intensity, frequency, and duration over Northern China, and the maximal values of decrease (24–56%) will be seen in Northeast China.

Description: The effects of conversion from staple rice to forage rice on carbon and greenhouse gas (GHG) balances in a paddy field were evaluated. A staple rice plot without the application of livestock manure compost (LMC, S − M plot) and forage rice plots with and without the application of LMC, derived mainly from cattle (2 kg−FW m−2, F + M and F − M plots, respectively), were established. CH4 and N2O fluxes and CO2 flux from a bare soil plot for organic matter decomposition (OMD) were measured. The carbon budget was calculated by subtracting the OMD, CH4 emission, and harvested grain and straw (forage rice only) from the net primary production and LMC. The net GHG balance was calculated by integrating them as CO2 equivalents. There were no significant differences in GHG flux among the plots. Compared to the carbon loss in the S − M plot, the loss increased by harvesting straw and was mitigated by LMC application. The net GHG emission in the F + M plot was significantly lower than that in other plots (1.78 and 2.63−2.77 kg CO2-eq m−2 year−1, respectively). There is a possibility that GHG emissions could be suppressed by forage rice cultivation with the application of LMC.

Description: Precise source identification for ambient pollution incidents in industrial parks were often difficult due to limited measurements. Source area analysis method was one of the applicable source identification methods, which could provide potential source areas under these circumstances. However, a source area usually covered several sources and the method was unable to identify the real one. This article introduces a case study on the statistical source identification of methyl mercaptan based on the long-term measurements, in 2014, in an industrial park. A procedure for statistical source area analysis was established, which contains independent pollution episode extraction, source area calculation scenario definition, meteorological data selection, and source area statistical analysis. A total of 414 violation records were detected by five monitors inside the park. Three kinds of calculation scenarios were found and, finally, three key source areas were revealed. The typical scenarios of source area calculations were described in detail. The characteristics of the statistical source areas for all pollution episodes were examined. Finally, the applicability of the method, as well as the source of uncertainties, was discussed. This study shows that more concentrated source areas can be identified through the statistical source area method if several excessive emission sources exist in an industrial park.

Description: We present findings from the Measurements of Urban, Marine and Biogenic Air (MUMBA) campaign, which took place in the coastal city of Wollongong in New South Wales, Australia. We focus on a few key air quality indicators, along with a comparison to regional scale chemical transport model predictions at a spatial resolution of 1 km by 1 km. We find that the CSIRO chemical transport model provides accurate simulations of ozone concentrations at most times, but underestimates the ozone enhancements that occur during extreme temperature events. The model also meets previously published performance standards for fine particulate matter less than 2.5 microns in diameter (PM2.5), and the larger aerosol fraction (PM10). We explore the observed composition of the atmosphere within this urban air-shed during the MUMBA campaign and discuss the different influences on air quality in the city. Our findings suggest that further improvements to our ability to simulate air quality in this coastal city can be made through more accurate anthropogenic and biogenic emissions inventories and better understanding of the impact of extreme temperatures on air quality. The challenges in modelling air quality within the urban air-shed of Wollongong, including difficulties in accurate simulation of the local meteorology, are likely to be replicated in many other coastal cities in the Southern Hemisphere.

Description: Besides local emissions, long-range transportation of polluted air masses also has a huge impact on haze pollution. In this study, the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model was used to determine the transport paths and potential sources of haze pollution in the Yangtze River Delta Urban Agglomeration. Haze days were determined by setting the threshold of meteorological elements. Shanghai, Hangzhou, Nanjing and Hefei were selected as four representative cities to calculate the −72 h backward transport trajectory of haze air mass; thus, the main transport path was obtained after clustering. A potential source contribution function and concentration weighted field were used to identify potential pollution sources of the study. The results showed that the number of haze days in the northern Yangtze River Delta Urban Agglomeration is much higher than that in the south. Haze days and Fine particulate matter (PM2.5) concentration showed a downward trend. The transport paths could be summarized as long-range transports from the northwest and coastal direction during the dry season and short-distance transports from all directions. −72 h air flow trajectories come from the higher altitudes in dry season than these in wet season. The main sources of potential pollution are Hebei, Shandong, Anhui and northern Jiangsu.

Description: The author would like to correct a published article by Teixeira [1], in which there is a factor of 2 missing from his Equation (24).[...]

Description: Sulfate aerosol (SO42−) is a major component of particulate matter in Japan. The Japanese model intercomparison study, J-STREAM, found that although SO42− is well captured by models, it is underestimated during winter. In the first phase of J-STREAM, we refined the Fe- and Mn-catalyzed oxidation and partly improved the underestimation. The winter haze in December 2016 was a target period in the second phase. The results from the Community Multiscale Air Quality (CMAQ) and Comprehensive Air quality Model with eXtentions (CAMx) regional chemical transport models were compared with observations from the network over Japan and intensive observations at Nagoya and Tokyo. Statistical analysis showed both models satisfied the suggested model performance criteria. CMAQ sensitivity simulations explained the improvements in model performance. CMAQ modeled lower SO42− concentrations than CAMx, despite increased aqueous oxidation via the metal catalysis pathway and NO2 reaction in CMAQ. Deposition explained this difference. A scatter plot demonstrated that the lower SO42− concentration in CMAQ than in CAMx arose from the lower SO2 concentration and higher SO42− wet deposition in CMAQ. The dry deposition velocity caused the difference in SO2 concentration. These results suggest the importance of deposition in improving our understanding of ambient concentration behavior.

Description: The aim of this research is to present a weather-based forecasting system for apple fire blight (Erwinia amylovora) and downy mildew of grapevine (Plasmopara viticola) under Serbian agroecological conditions and test its efficacy. The weather-based forecasting system contains Numerical Weather Prediction (NWP) model outputs and a disease occurrence model. The weather forecast used is a product of the high-resolution forecast (HRES) atmospheric model by the European Centre for Medium-Range Weather Forecasts (ECMWF). For disease modelling, we selected a biometeorological system for messages on the occurrence of diseases in fruits and vines (BAHUS) because it contains both diseases with well-known and tested algorithms. Several comparisons were made: (1) forecasted variables for the fifth day are compared against measurements from the agrometeorological network at seven locations for three months (March, April, and May) in the period 2012–2018 to determine forecast efficacy; (2) BAHUS runs driven with observed and forecast meteorology were compared to test the impact of forecasted meteorological data; and (3) BAHUS runs were compared with field disease observations to estimate system efficacy in plant disease forecasts. The BAHUS runs with forecasted and observed meteorology were in good agreement. The results obtained encourage further development, with the goal of fully utilizing this weather-based forecasting system.

Description: Comparisons between bin and bulk cloud microphysics schemes are conducted by simulating a heavy precipitation case using a bin microphysics scheme and four double-moment bulk microphysics schemes in the Weather Research and Forecasting (WRF) model. For this, we implemented an updated bin microphysics scheme in the WRF model. All of the microphysics schemes underestimate observed strong precipitation, but the bin microphysics scheme yields the result that is closest to observations. The differences among the schemes are more pronounced in terms of hydrometeor number concentration than in terms of hydrometeor mixing ratio. In this case, the bin scheme exhibits remarkably more latent heat release by deposition and riming than the bulk schemes. This causes stronger updrafts and more upward transport of water vapor, which leads to more deposition, and again, increases the latent heat release. An additional simulation using the bin scheme but excluding the riming of cloud droplets on ice crystals, which is not or poorly treated in the examined bulk schemes, shows that surface precipitation is slightly weakened and moved farther downwind compared to that of the control simulation. This implies that the more appropriate representation of microphysical processes in the bin microphysics scheme contributes to the more accurate prediction of precipitation in this case.

Description: The performance of recent reanalysis products (i.e., ERA-Interim, NCEP2, MERRA, CFSR, and JRA-55) was evaluated based on in situ observations from nine automatic weather stations and one stake network to investigate the monthly and seasonal variability of the surface mass balance in Antarctica. Synoptic precipitation simulations were also evaluated by an investigation of high precipitation events. The seasonal variations showed large fluctuations and were inconsistent at each station, probably owing to the large interannual variability of snow accumulation based on the short temporal coverage of the data. The ERA-Interim and JRA-55 datasets revealed better simulated precision, with the other three models presenting similar simulations at monthly and seasonal timescales. The JRA-55 dataset captured a greater number of synoptic high precipitation events at four of the nine stations. Such events at the other five stations were mainly captured by ERA and CFSR. The NCEP2 dataset was more weakly correlated with each station on all timescales. These results indicate that significant monthly or seasonal correlations between in situ observations and the models had little effect on the capability of the reanalyses to capture high precipitation events. The precision of the five reanalysis datasets widely fluctuated in specific regions or at specific stations at different timescales. Great caution is needed when using a single reanalysis dataset to assess the surface mass balance over all of Antarctica.

Description: The main goal of the study is to determine the pressure conditions that cause waves of hot nights in Central Europe. The goal was implemented on the basis of data from 1966 to 2015, made available by the Institute of Meteorology and Water Management—National Research Institute, Deutscher Wetterdienst and the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR). A hot night was defined as a day with the minimum air temperature above 18 °C. In the analysed years, an increase in the number of hot nights was found, which was predominantly statistically significant within the studied area. The study shows that the occurrence of waves of hot nights in Central Europe was associated on average with the ridge of high pressure, within which a local high-pressure area developed. During the waves of hot nights, there were positive anomalies of heights of isobaric surfaces over the study area with a maximum in the upper troposphere.

Description: A limited-area kilometre scale numerical weather prediction system is applied to evaluate the effect of refined surface data assimilation on short-range heavy precipitation forecasts. The refinements include a spatially dependent background error representation, use of a flow-dependent data assimilation technique, and use of data from a satellite-based scatterometer instrument. The effect of the enhancements on short-term prediction of intense precipitation events is confirmed through a number of case studies. Verification scores and subjective evaluation of one particular case points at a clear impact of the enhanced surface data assimilation on short-range heavy precipitation forecasts and suggest that it also tends to slightly improve them. Although this is not strictly statistically demonstrated, it is consistent with the expectation that a better surface state should improve rainfall forecasts.

Description: Strong El Niño events alter tropical climates and may lead to a negative carbon balance in tropical forests and consequently a disruption to the global carbon cycle. The complexity of tropical forests and the lack of data from these regions hamper the assessment of the spatial distribution of El Niño impacts on these ecosystems. Typically, maps of climate anomaly are used to detect areas of greater risk, ignoring baseline climate conditions and forest cover. Here, we integrated climate anomalies from the 1982–1983, 1997–1998, and 2015–2016 El Niño events with baseline climate and forest edge extent, using a risk assessment approach to hypothetically assess the spatial and temporal distributions of El Niño risk over tropical forests under several risk scenarios. The drivers of risk varied temporally and spatially. Overall, the relative risk of El Niño has been increasing driven mainly by intensified forest fragmentation that has led to a greater chance of fire ignition and increased mean annual air temperatures. We identified areas of repeated high risk, where conservation efforts and fire control measures should be focused to avoid future forest degradation and negative impacts on the carbon cycle.

Description: In this study, aridity data and tree ring data were collected in Northern Serbia, in Southeast (SE) Banat, a subregion within Vojvodina, and Vojvodina at large. They were each investigated independently. The De Martonne Aridity Index and the Forestry Aridity Index are derived from examining the relationship between precipitation and surface air temperature data sets sourced from seven meteorological stations in SE Banat, and from 10 meteorological stations located in Vojvodina as a whole. Vojvodina is a large territory and used as the control area, for the period 1949–2017. The Palmer Drought Severity Index was derived for the period 1927–2016, for both SE Banat and the totality of Vojvodina. The results of the Tree Ring Width Index were obtained from samples collected in or around the villages of Vlajkovac and Šušara, both located in SE Banat, for the period 1927–2017. These tree ring records were compared with three previous aridity and drought indices, and the meteorological data on the surface air temperature and the precipitation, with the objective being to evaluate the response of tree growth to climate dynamics in the SE Banat subregion. It was noted that the significant positive temperature trends recorded in both areas were too insufficient to trigger any trends in aridity or the Tree Ring Width Index, as neither displayed any change. Instead, it appears that these climatic parameters only changed in response to the precipitation trend, which remained unchanged during the investigated period, rather than in response to the temperature trend. It appears that the forest vegetation in the investigated areas was not affected significantly by climate change in response to the dominant temperature increase.

Description: Rainfall patterns in the San Francisco Bay Area (SFBA) are highly influenced by local topography. It has been a forecasting challenge for the main US forecast models. This study investigates the ability of the Weather Research and Forecasting (WRF) model to improve upon forecasts, with particular emphasis on the rain shadow common to the southern end of the SFBA. Three rain events were evaluated: a mid-season atmospheric river (AR) event with copious rains; a typical non-AR frontal passage rain event; and an area-wide rain event in which zero rain was recorded in the southern SFBA. The results show that, with suitable choices of parameterizations, the WRF model with a resolution around 1 km can forecast the observed rainfall patterns with good accuracy, and would be suitable for operational use, especially to water and emergency managers. Additionally, the three synoptic situations were investigated for further insight into the common ingredients for either flooding rains or strong rain shadow events.