The modern era of polarimetric radar begins with radiowave propagation research starting in the early 1970s with applications to measurement and modeling of wave attenuation in rain and depolarization due to ice particles along satellite-earth links. While there is a rich history of radar in meteorology after World War II, the impetus provided by radiowave propagation requirements lead to high quality antennas and feeds. Our journey starts by describing the key institutions and personnel responsible for development of weather radar polarimetry. The early period was dominated by circularly polarized radars for propagation research and at S-band for hail detection. By the mid-to late 70s, a paradigm shift occurred which led to the dominance of linear polarizations with applications to slant path attenuation prediction as well as estimation of rain rates and inferences of precipitation physics. The period from early 1980s to 1995 can be considered as the “golden” period of rapid research that brought in meteorologists, cloud physicists and hydrologists. This article describes the evolution of this technology from the vantage point of the authors. Their personal reflections and “behind the scenes” descriptions offer a glimpse into the inner workings at several key institutions which cannot be found elsewhere.

This paper characterizes the influence of synoptic-scale air mass conditions on spatial and temporal patterns of precipitation in North Carolina over a 16-year period (2003-2018). National Center for Environmental Prediction Stage IV multi-sensor precipitation estimates were used to describe seasonal variations in precipitation in the context of prevailing air mass conditions classified using the spatial synoptic classification system. Spatial analyses identified significant clustering of high daily precipitation amounts distributed along the east side of the Appalachian Mountains and along the coastal plains. Significant and heterogeneous clustering was prevalent in summer months and tended to coincide with land cover boundaries and complex terrain. The summer months were dominated by maritime tropical air mass conditions whereas dry moderate air mass conditions prevailed in the winter, spring, and fall. Between the three geographic regions of North Carolina, highest precipitation amounts were received in western North Carolina during the winter and spring, and in eastern North Carolina in the summer and fall. Central North Carolina received the least amount of precipitation; however, there was substantial variability between regions due to prevailing air mass conditions. There was an observed shift toward warmer and more humid air mass conditions in the winter, spring, and fall months throughout the study period (2003-2018), indicating a shift toward air mass conditions conducive to higher daily average rain rates in North Carolina.

Chemical processing of organic material in aqueous atmospheric aerosols and cloudwater is known to form secondary organic aerosols (SOA), although the extent to which each of these processes contributes to total aerosol mass is unclear. In this study, we use GAMMA 5.0, a photochemical box model with coupled gas and aqueous-phase chemistry, to consider the impact of aqueous organic reactions in both aqueous aerosols and clouds on isoprene epoxydiol (IEPOX) SOA over a range of pH for both aqueous phases, including cycling between cloud and aerosol within a single simulation. Low-pH aqueous aerosol, in the absence of organic coatings or other morphology which may limit uptake of IEPOX, is found to be an efficient source of IEPOX SOA, consistent with previous work. Cloudwater at pH 4 or lower is also found to be a potentially significant source of IEPOX SOA. This phenomenon is primarily attributed to the relatively high uptake of IEPOX to clouds as a result of higher water content in clouds as compared to aerosol. For more acidic cloudwater, the aqueous organic material is comprised primarily of IEPOX SOA and lower-volatility organic acids. For both cloudwater and aqueous aerosol, pH is the most significant factor considered in this study in determining the mass of aqueous phase organic acids and IEPOX SOA. Other factors, such as the time of day or sequence of aerosol-to-cloud or cloud-to-aerosol transitions, contribute to less than 15% difference in the final aqSOA fractional composition. The potential significance of cloud processing as a contributor to IEPOX SOA production could account for discrepancies between predicted IEPOX SOA mass from atmospheric models and measured ambient IEPOX SOA mass, or observations of IEPOX SOA in locations where mass transfer limitations are expected in aerosol particles.

An open source computational fluid dynamics (CFD) solver has been incorporated into the WindNinja modeling framework widely used by wildland fire managers as well as researchers and practitioners in other fields, such as wind energy, wind erosion, and search and rescue. Here we describe incorporation of the CFD solver and evaluate its performance compared to the conservation of mass (COM) solver in WindNinja and previously published large-eddy simulations (LES) for three field campaigns conducted over isolated terrain obstacles of varying terrain complexity: Askervein Hill, Bolund Hill, and Big Southern Butte. We also compare the effects of two important model settings in the CFD solver and provide guidance on model sensitivity to these settings. Additionally, we investigate the computational mesh and difficulties regarding terrain representation. Two important findings from this work are: (1) the choice of discretization scheme for advection has a significantly larger effect on the simulated winds than the choice of turbulence model and (2) CFD solver predictions are significantly better than the COM solver predictions at windward and lee side observation locations, but no difference was found in predicted speed-up at ridgetop locations between the two solvers.

The 24 September 2001 College Park, Maryland, tornado was a long-track and strong tornado that passed within a close range of two Doppler radars. It was the third in a series of three tornadoes associated with a supercell storm that developed in Stafford County, Virginia, and initiated 3 - 4 km southwest of College Park and dissipated near Columbia, Howard County. The supercell tracked approximately 120 km and lasted for about 126 minutes. This study presents a synoptic and mesoscale overview of favorable conditions and forcing mechanisms that resulted in the severe convective outbreak associated with the College Park tornado. Results show many critical elements of the tornadic event, including a negative-tilted upper-level trough over the Ohio Valley, a jet stream with moderate vertical shear, a low-level warm, moist tongue of the air associated with strong southerly flow over south-central Maryland and Virginia, and significantly increased convective available potential energy (CAPE) during the late afternoon hours. A possible role of the urban heat island effects from Washington, DC in increasing CAPE for the development of the supercell is discussed. Satellite imagery reveals banded convective morphology with high cloud tops associated with the supercell that produced the College Park tornado. Operational WSR-88D data exhibits a high reflectivity “debris ball” or tornadic debris signature (TDS) within the hook echo, the evolution of the parent storm from a supercell structure to a bow echo, and a tornado cyclone signature (TCS). Many of the mesoscale features could be captured by contemporary numerical model analyses. This study concludes with a discussion of the effectiveness of the coordinated use of satellite and radar observations in the operational environment of nowcasting severe convection.

Land surface temperature (LST) is a key variable in surface-atmosphere energy and water exchanges. The main goals of this study are to (i) evaluate the LST of the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim and ERA5 reanalyses over Iberian Peninsula using the Satellite Application Facility on Land Surface Analysis (LSA-SAF) product and to (ii) understand the main drivers of the LST errors in the reanalysis. Simulations with the ECMWF land-surface model in offline mode (uncoupled) were carried out over the Iberian Peninsula and compared with the reanalysis data. Several sensitivity simulations were performed in a confined domain centered in Southern Portugal to investigate potential sources of the LST errors. The Copernicus Global Land Service (CGLS) fraction of green vegetation cover (FCover) and the European Space Agency’s Climate Change Initiative (ESA-CCI) Land Cover dataset were explored. We found a general underestimation of daytime LST and slightly overestimation at night-time. The results indicate that there is still room for improvement in the simulation of LST in ECMWF products. Still, ERA5 presents an overall higher quality product in relation to ERA-Interim. Our analysis suggested a relation between the large daytime cold bias and vegetation cover differences between (ERA5 and CGLS FCocver) with a correlation of -0.45. The replacement of the low and high vegetation cover by those of ESA-CCI provided an overall reduction of the large Tmax biases during summer. The increased vertical resolution of the soil at the surface, has a positive impact, but much smaller when compared with the vegetation changes. The sensitivity of the vegetation density parameter, that currently depends on the vegetation type, provided further proof for a needed revision of the vegetation in the model, as there is a reasonable correlation between this parameter and the Tmax mean errors when using the ESA-CCI vegetation cover (while the same correlation cannot be reproduced with the original model vegetation). Our results support the hypothesis that vegetation cover is one of the main drivers of the LST summertime cold bias in ERA5 over Iberian Peninsula.

The present study uses the 3 years of data from a vertically pointed profile Radar (VPPR) at the mountain site (Huancayo at 3.2 km msl; HYO) in Peru to investigate the precipitation characteristics/patterns including diurnal variation, bright band (BB) characteristics and vertical structure of rain (VSR). VPPR provides continuous 1­min profiles of radar reflectivity (Ze), rain rate (RR), liquid water content (LWC) and Doppler velocity. At the same time, field campaign data are used to investigate the role of the surface and atmospheric variables in generating the rainfall and organization of the cloud systems over HYO. The precipitation shows the bimodal distribution; frequent and higher RR during afternoon and nighttime. The BB height also shows a diurnal pattern and BB height is higher during the afternoon and nighttime compared; and follows the diurnal heating of the Andes. More than 70% BB height lies between ~4.3­4.7 km and only 20% BB has altitude higher than 5 km. The austral summer monsoon (December to March months) have higher intense vertical profiles (higher Z e ) compared to pre monsoon seasons (September to November) and shows the nega­tive gradient for most of the altitude. The RR and LWC show the opposite characteristics, and both have a positive gradient below the 6 km altitude and then negative gradient above 6 km altitude. The raindrop size distribution (DSD) parameters show most of the variation below the freezing level, and a higher concentra­tion of large sized of hydrometeors are observed for higher RR, however the dominant modes of Dm are less than 1 mm. The changes in the VSR in DSD parameters are significant for the RR>20 mm/h, whereas for RR<2 mm/h the vertical structure in DSD parameters do not have much differences during monsoon and premonsoon seasons. Satellite and reanalysis data reveal the short periods of convective activity with higher ac­ cumulated rainfall over HYO compared to stratiform precipitation, which is more common in the nighttime and sustain for many hours. Wet spells (higher rainy days) have low pressure circulation, which favours the transport of moisture from the Amazon to the central Sierra of Peru, while the anticyclonic circulation at high levels favours the divergence at higher pressure levels and, enhances the convective in the region. During the dry spells, low level weaker circulation at the west of Brazil, restricts the transport of moisture to the central Sierra, while the circulation at high levels does not favor rain processes. The improved understanding of the tropical Andes precipitation would be very important for assessing climate variability and changes as well as to forecast precipitation using the numerical models.

The gap in knowledge regarding the radiative effects of aerosols in the UV region of the solar spectrum is large, mainly due to the lack of systematic measurements of the aerosol single scattering albedo (SSA) and absorption optical depth (AAOD). In the present study, spectral UV measurements performed in Thessaloniki, Greece by a double monochromator Brewer spectrophotometer in the period 1998 - 2017 are used for the calculation of the aforementioned optical properties. The main uncertainty factors have been described and there is an effort to quantify the overall uncertainties in SSA and AAOD. Analysis of the results suggests that the absorption by aerosols is much stronger in the UV relative to the visible. SSA follows a clear annual pattern ranging from ~0.7 in winter to ~0.85 in summer at wavelengths 320 – 360 nm, while AAOD peaks in summer and winter. The average AAOD for 2009 – 2011 is ~50% above the 2003 – 2006 average, possibly due to increased emissions of absorbing aerosols related to the economical crisis and the metro-railway construction works in the city center. A detailed analysis of the uncertainties in the retrieval of the SSA and the AAOD from the Brewer spectrophotometer has been also performed.

Tornadoes remain an active subject of observational and numerical research due to the damage and fatalities they cause worldwide as well as poor understanding of their behavior, such as what processes result in their genesis and what determines their longevity and morphology. A numerical model executed on a supercomputer run at high resolution can serve as a powerful tool for exploring the rapidly evolving tornado-scale features within a simulated storm, but saving large amounts data for meaningful analysis can result in unacceptably slow model performance, an unwieldy amount of saved data, and saved data spread across millions of files. In this paper, a system for efficiently saving and managing hundreds of terabytes of compressed model output is described in order to support a supercomputer simulation of a tornadic supercell thunderstorm. The challenges of managing a simulation spanning over a quarter trillion grid volumes across the Blue Waters supercomputer are also described. The simulated supercell produces a long-track EF5 tornado, and the near-tornado environment is described during tornadogenesis, where single upward-growing vortex undergoes several vortex mergers before transitioning into a multiple vortex tornado.

Ocean surface heat fluxes play a significant role in the genesis and evolution of various marine-based atmospheric phenomena, from the synoptic scale down to the microscale. While in-situ measurements from buoys and flux towers will continue to be the standard in regards to surface heat flux estimates, they commonly have significant gaps in temporal and spatial coverage. Previous and current satellite missions have filled these gaps; though they may not observe the fluxes directly, they can measure the variables needed (wind speed, temperature, and humidity) to estimate latent and sensible heat fluxes. However, current remote sensing instruments have their own limitations, such as infrequent coverage, signals attenuated by precipitation, or both. The Cyclone Global Navigation Satellite System (CYGNSS) mission overcomes these limitations over the tropical and subtropical oceans by providing improved coverage in nearly all weather conditions. While CYGNSS (Level 2) primarily estimates surface winds, when coupled with observations or estimates of temperature and humidity from reanalysis data, it can provide estimates of latent and sensible heat fluxes along its orbit. This paper describes the development of the Surface Heat Flux Product for the CYGNSS mission, its current results, and expected improvements and changes in future releases.

This paper characterizes the influence of synoptic-scale airmass conditions on spatial and temporal patterns of precipitation in North Carolina over a 16-year period (2003-2018). National Center for Environmental Prediction Stage IV multi-sensor precipitation estimates were used to describe seasonal variations in precipitation in the context of prevailing airmass conditions classified using the spatial synoptic classification system. Spatial analyses identified significant clustering of high daily average precipitation amounts distributed along the lee side of the Appalachian Mountains and along the coastal plains. Significant and heterogenous clustering was prevalent in summer months and tended to coincide with land cover boundaries and complex terrain. Between the three geographic regions of North Carolina, highest precipitations amounts were received in western North Carolina during the winter and spring, but this signal shifted to eastern North Carolina in the summer and fall. Central North Carolina received the least amount of precipitation; however, there was substantial variability between regions due to prevailing airmass conditions. The summer months were dominated by maritime tropical airmass conditions with no clear shift in summertime airmass trends over the study period. Most days with recorded precipitation in the winter, spring, and fall occurred under dry moderate airmass conditions; however, the highest daily average precipitation and total precipitation occurred under the influence of maritime moderate airmasses. Importantly, there was an observed shift toward warmer and more humid airmass conditions in the winter, spring, and fall months throughout the study period (2003-2018), indicating a shift toward airmass conditions conducive to higher daily average rain rates in North Carolina.

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 at 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.695 ng m-2 h-1 with a mean value equal to 10.501 ng m-2 h-1 ± 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 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.

Improvements that can be attained in seasonal climate predictions in various parts of Africa using the multimodel supersensemble scheme are presented in this study. The synthetic superensemble (SSE) used follows the approach originally developed at Florida State University (FSU). The technique takes more advantage of the skill in the climate forecast data sets from atmosphere-ocean general circulation models running at many centres worldwide including the WMO global producing centers (GPCs). The module used in this work drew data sets from the Four versions of FSU coupled model system, seven models from the DEMETER project which is the forerun to the current European Ensembles Forecast System, the NCAR Model, and the Predictive Ocean Atmosphere Model for Australia (POAMA), all making a set of 13 individual models. An archive consisting of monthly simulations of precipitation was available over all the 5 regions of Africa, namely Eastern, Central, Northern, Southern, and Western Africa. The results showed that the SSE forecast for precipitation carries a higher skill compared to each of the member models and the ensemble mean. Relative to the ensemble mean (EM), the SSE provides an improvement of 18% in simulation of season cycle of precipitation climatology. In Eastern Africa, during December-February season, a north-south gradient of precipitation prevails between Tropical East Africa and the sector of the region towards Southern Africa. This regional scale climate pattern is a direct influence of the Intertropical Convergence Zone (ITZC) across the African continent during this time of the year. The SSE emerges with superior skill scores such as lowest root mean square error above the EM and the member models, for example in the prediction of spatial location and precipitation magnitudes that characterize the see-saw precipitation pattern in Eastern Africa. In all parts of Africa, and especially Eastern Africa where seasonal precipitation variability is a frequent cause huge human suffering in due to droughts and famine, the multimodel superensemble and its subsequent improvements will always provide a forecast that out weighs the best Atmosphere-Ocean Climate Model.

Atmospheric Motion Vectors (AMVs) calculated by six different institutions (Brazil Center for Weather Prediction and Climate Studies/CPTEC/INPE, European Organization for the Exploitation of Meteorological Satellites/EUMETSAT, Japan Meteorological Agency/JMA, Korea Meteorological Administration/KMA, Unites States National Oceanic and Atmospheric Administration/NOAA and the Satellite Application Facility on Support to Nowcasting/NWCSAF) with JMA’s Himawari-8 satellite data and other common input data are here compared. The comparison is based on two different AMV input datasets, calculated with two different image triplets for 21 July 2016, and the use of a prescribed and a specific configuration. The main results of the study are summarized as follows: (1) the differences in the AMV datasets depend very much on the “AMV height assignment” used and much less on the use of a prescribed or specific configuration; (2) the use of the “Common Quality Indicator (CQI)” has a quantified skill in filtering collocated AMVs for an improved statistical agreement between centers; (3) JMA AMV algorithm has the best overall performance considering all validation metrics, most likely due to its height assignment: “optimal estimation using observed radiance and NWP wind vertical profile”.

This study proposes a new algorithm termed rain cell identification and tracking (RCIT) to identify and track rain cells from high resolution weather radar data. Previous algorithms have limitations when tracking non-consequent rain cells owing to their use of maximum correlation coefficient methods and their lack of an alternative way to handle the variation stages of rain cells during their life cycles. To address these deficiencies, various methods are implemented in the new algorithm. These include the particle image velocimetry (PIV) method for motion estimation and the rain cell matching rule to obtain the stage changes of rain cells. High resolution (5-min and 1-km) radar reflectivity data from three rainy days over the German federal state North Rhine Westphalia (NRW) are used to evaluate the proposed algorithm. The performance of the new algorithm is compared with a radar reflectivity map and verified by two object-oriented methods: structure–amplitude–location (SAL) and geometric index. The verification results suggest that the performance of the new algorithm is good. Application of the RCIT algorithm to the selected cases shows that the inner structure of rainfall events in the experimental region present extreme value distributions, with most rainfall events having a short duration with less intensity. The new algorithm can effectively capture the stage changes of rain cells during their life cycles. The proposed algorithm can serve as the basis for further hydro-meteorological applications such as spatial and temporal analysis of rainfall events and short-term flood forecasting.

In the context of the French-German space lidar mission MERLIN dedicated to the determination of the atmospheric methane content, an end-to-end mission simulator is being developed. In order to check whether the instrument design meets the performance requirements, simulations have to count all the sources of noise on the measurements like the optical energy variability induced by speckle. Speckle is due to interference as the lidar beam are quasi monochromatic. Speckle contribution to the error budget has to be estimated but also simulated. In this paper, the speckle theory is revisited and applied to MERLIN double pulsed IPDA lidar and also to the DLR demonstrator CHARM-F. Results show: on the signal path, speckle noise depends mainly on the size of the illuminating area on ground; on the solar flux, speckle is fully negligible both because the pixel size and the optical filter spectral width; on energy monitoring path a decorrelation mechanism is needed to reduce speckle noise on averaged data. Speckle noises for MERLIN and CHARM-F can be simulated by Gaussian noises with only one random draw by shot separately for energy monitoring and signal paths.

Visualization is an essential tool for analysis of data and communication of findings in the sciences, and the Earth System Science (ESS) are no exception. However, within ESS specialized visualization requirements and data models --- particularly for those data arising from numerical models --- often make general-purpose visualization packages difficult, if not impossible, to effectively use. This paper presents VAPOR: a domain-specific visualization package that targets the specialized needs of ESS modelers, particularly those working in research settings where highly interactive exploratory visualization is beneficial. We specifically describe VAPOR’s ability to handle ESS simulation data from a wide variety of numerical models, as well as a multi-resolution representation that enables interactive visualization on very large data while using only commodity computing resources. We also describe VAPOR’s visualization capabilities, paying particular attention to features for geo-referenced data and advanced rendering algorithms suitable for time-varying, 3D data. Finally, we illustrate VAPOR's utility in the study of a numerically simulated tornado. Our results demonstrate both ease-of-use and the rich capabilities of VAPOR in such a use case.

The “NWCSAF High Resolution Winds (NWC/GEO-HRW)” software is developed by the EUMETSAT’s “Satellite Application Facility on support to Nowcasting and very short range forecasting (NWCSAF)”, inside its stand-alone software package for calculation of meteorological products with geostationary satellite data (NWC/GEO). The whole NWC/GEO software package can be obtained after registration at the NWCSAF Helpdesk, www.nwcsaf.org. It is easy to get, install and use. The code is easy to read and fully documented. And in the NWCSAF Helpdesk, users find support and help for its use. “NWCSAF High Resolution Winds” provides a detailed calculation of Atmospheric Motion Vectors (AMVs) and Trajectories, locally and in near real time, using as input NWP model data and geostationary satellite image data. The latest version of the software, v2018, is able to process MSG, Himawari-8/9, GOES-N and GOES-R satellite series images, so that AMVs and Trajectories can be calculated all throughout the planet Earth with the same algorithm and quality. In the “2014 and 2018 AMV Intercomparison Studies”, “NWCSAF High Resolution Winds” has shown to be one of the two best AMV algorithms for both MSG and Himawari-8/9 satellites. And the “Coordination Group for Meteorological Satellites (CGMS)” has recognized in its “2012 Meeting Report”: 1. “NWCSAF High Resolution Winds” fulfills the requirements to be a portable stand-alone AMV calculation software due to its easy installation and usability. 2. It has been successfully adapted by some CGMS members and serves as an important tool for development. It is modular, well documented, and well suited as stand-alone AMV software. 3. Although alternatives exist as portable stand-alone AMV calculation software, they are not as advanced in terms of documentation and do not have an existing Helpdesk. Considering this, a full description and validation of the “NWCSAF/High Resolution Winds” is shown here for the first time in a peer-reviewed paper. The procedure to obtain the software for operational meteorology and research is also explained.

The contribution of extra-tropical synoptic cyclones to the formation of summer-mean atmospheric circulation patterns in the Arctic is investigated by clustering the dominant Arctic circulation patterns by the self-organizing maps (SOMs) using the daily mean sea level pressure (MSLP) in the Arctic domain (≥ 60°N). Three SOM patterns are identified: one with prevalent low pressure anomalies in the Arctic Circle (SOM1) and two opposite dipoles with primary high pressure anomalies covering the Arctic Ocean (SOM2 and SOM3). The time series of summertime occurrence frequencies demonstrate the largest inter-annual variation in the SOM1, the slight decreasing trend in the SOM2, and the abrupt upswing after 2007 in the SOM3. The relevant analyses with produced cyclone track data confirm that the vital contribution. The Arctic cyclone activity is enhanced in the SOM1 because the meridional temperature gradient increases over the land–Arctic Ocean boundaries co-located with major extra-tropical cyclone pathways. The composite daily synoptic evolutions for each SOM reveal that the persistence of all the three SOMs is less than 5 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.

The 24 September 2001 College Park, Maryland, tornado was a long track and strong tornado that passed within a close range of two Doppler radars. It was the third in a series of three tornadoes associated with a supercell storm that developed in Stafford County, Virginia, and initiated 3 - 4 km southwest of College Park and dissipated near Columbia, Howard County. The supercell tracked approximately 120 km and lasted for about 126 minutes. This study presents a synoptic and mesoscale overview of favorable conditions and forcing mechanisms that resulted in the severe convective outbreak associated with the College Park tornado. Results show many critical elements of the tornadic event, including a negative-tilted upper-level trough over the Ohio Valley, a jet stream with moderate vertical shear, a warm, moist tongue of the air associated with strong southerly flow over south-central Maryland and Virginia, and significantly increased convective available potential energy during the late afternoon hours. Satellite imagery reveals banded convective morphology with high cloud tops associated with the supercell that produced the College Park tornado. Operational WSR-88D data exhibits a high reflectivity “debris ball” or tornadic debris signature (TDS) within the hook echo, the evolution of the parent storm from a supercell structure to a bow echo, and a tornado cyclone signature (TCS). Many of the mesoscale environmental features could be captured by contemporary numerical model analyses. This study concludes with a discussion of the effectiveness of the coordinated use of satellite and radar observations in the operational environment of nowcasting severe convection.

This study analyses uncertainties associated with the annual cycle of West African rainfall characteristics in 15 simulations of the Coupled Model Intercomparison Project phase 5 (CMIP5) over the Sahel and Guinean regions. Indices based on daily rainfall such as the frequency and the ntensity of wet days, the consecutive dry days (CDD) and wet days (CWD), the 95th percentile of daily rainfall (R95) and its contribution to the umulative monsoon rainfall (R95PTOT) have been assessed. Over both regions, TRMM, GPCP and CHIRPS observational datasets provide very consistent results on the annual cycle of precipitation but less so on the frequency of wet days. Conversely, higher uncertainties are noted on the intensity of wet days over both study areas, particularly over the Guinean region. Overall, CMIP5 simulations present much higher uncertainties in the representation of the mean precipitation climatology, often provide too early (late) onset dates over the Sahel (the Guinean region) and overestimate rainfall during the early and late monsoon phases. These errors do not compensate at the annual scale nor when considering West Africa as a hole. Results also reveal that over the Guinean region, the difficulty of models to represent the annual structure of the mean precipitation strongly involves biases in the representation of the annual cycle of the frequency of wet days. We found strong uncertainties in the simulation of the CWD and he CDD over both areas. Conversely for R95 and R95PTOT, the ncertainties in CMIP5 models appear somewhat weaker, but the magnitude of R95 is largely underestimated in most models.

This study analyses uncertainties associated with the annual cycle of West African rainfall characteristics in 15 simulations of the Coupled Model Intercomparison Project phase 5 (CMIP5) over the Sahel and Guinean regions. Indices based on daily rainfall such as the frequency and the ntensity of wet days, the consecutive dry days (CDD) and wet days (CWD), the 95th percentile of daily rainfall (R95) and its contribution to the umulative monsoon rainfall (R95PTOT) have been assessed. Over both regions, TRMM, GPCP and CHIRPS observational datasets provide very consistent results on the annual cycle of precipitation but less so on the frequency of wet days. Conversely, higher uncertainties are noted on the intensity of wet days over both study areas, particularly over the Guinean region. Overall, CMIP5 simulations present much higher uncertainties in the representation of the mean precipitation climatology, often provide too early (late) onset dates over the Sahel (the Guinean region) and overestimate rainfall during the early and late monsoon phases. These errors do not compensate at the annual scale nor when considering West Africa as a hole. Results also reveal that over the Guinean region, the difficulty of models to represent the annual structure of the mean precipitation strongly involves biases in the representation of the annual cycle of the frequency of wet days. We found strong uncertainties in the simulation of the CWD and he CDD over both areas. Conversely for R95p and R95PTOT, the ncertainties in CMIP5 models appear somewhat weaker, but the magnitude f R95 is largely underestimated in most models.

AIUS (Atmospheric Infrared Ultraspectral Sounder) is an infrared occultation spectrometer onboard the Chinese GaoFen-5 satellite, which covers a spectral range of 2.4--13.3 μm (750--4100 cm-1) with a spectral resolution of about 0.02 cm-1. AIUS was designed to measure and to study chemical processes of ozone (O3) and other trace gases in the upper troposphere and stratosphere over Antarctic. In this study, the corresponding retrieval methodology is described. The comparison between AIUS measurements and simulated spectra illustrates that AIUS measurements agree well to the simulated spectra. To first evaluate the reliability of our retrieval algorithm, three retrieval O3 experiments are performed based on ACE-FTS observation spectra. A comparison between our retrieved profiles and the ACE-FTS official products shows that the relative difference of these three retrieval experiments is mostly within 10% between 20 km and 70 km. These retrieval experiments demonstrate that the retrieval algorithm described in this study work fine and reliable. Furthermore, O3, H2O and HCl profiles are retrieved from eight orbits of AIUS measurements and compared with the official AURA/MLS level-2 v4.2 profiles. Comparison experiments show that the relative difference is mostly within 10% (about 0.02 - 0.4 ppm) between 18 and 58 km for O3 retrieval, within 10% (0-0.5 ppm) between 15 and 80 km for H2O retrieval, and within 10% (about 0.1 ppb) between 30 and 60 km for HCl retrieval. There is a good agreement in the retrieved trace gas profiles obtained from AIUS and from coincident profiles from MLS.

Radiation is the major driver of snowmelt, and hence its estimation is critically important. Net radiation reaching the forest floor is influenced by vegetation density. Previous studies in mid-latitude conifer forests have confirmed that net radiation decreases and then subsequently increases with increasing vegetation density, for clear sky conditions. This leads to existence of a net radiation minimum at an intermediate vegetation density. With increasing cloud cover, the minimum radiation shifts toward lower densities, sometimes resulting in a monotonically increasing radiation with vegetation density. The net radiation trend, however, is expected to change across sites, affecting the magnitude and timing of individual radiation components. This research explores the variability of net radiation on snow-covered forest floor for different vegetation densities along a latitudinal gradient. We especially investigate how the magnitude of minimum/maximum radiation and the corresponding vegetation density change with the site geographical location. To evaluate these, the net radiation is evaluated using the Forest Radiation Model at six different locations in predominantly white spruce (Picea glauca) canopy cover across North America, ranging from 45 to 66°N latitudes. Results show that the variation of net radiation with vegetation density considerably varies with latitude. In higher latitude forests, the magnitude of net radiation is generally smaller, and the minimum radiation is exhibited at relatively sparser vegetation densities, under clear sky conditions. For interspersed cloudy sky conditions, net radiation non-monotonically varies with latitude across the sites, depending on the seasonal sky cloudiness and air temperature. Latitudinal sensitivity of net radiation is lower on north-facing hillslopes than on south-facing sites.

The Sentinel Application Platform (SNAP) architecture facilitates Earth Observation data processing (http://step.esa.int/main/toolboxes/snap/). In this work we present results from a new Snow Processor for SNAP. We also describe physical principles behind the developed snow property retrieval technique based on the analysis of Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A/B measurements over clean and polluted snow fields. Using OLCI spectral reflectance measurements in the range 400-1020nm, we derive important snow properties such as spectral and broadband albedo, snow specific surface area, snow extent and grain size on the spatial grid of 300m. The algorithm also incorporates cloud screening and atmospheric correction procedures over snow surfaces. We present validation results using ground measurements from Antarctica, the Greenland ice sheet and the French Alps. We find the spectral albedo retrieved with accuracy of better than 3% on average, making our retrievals sufficient for a variety of applications. Broadband albedo is retrieved with the average accuracy of about 5% over snow. Therefore, the uncertainties of satellite retrievals are close to experimental errors of ground measurements. The retrieved surface grain size shows good agreement with ground observations. Snow specific surface area observations are also consistent with our OLCI retrievals. We present snow albedo and grain size mapping over the inland ice sheet of Greenland for areas including dry snow, melted/melting snow and impurity rich bare ice. The algorithm can be applied to OLCI Sentinel-3 measurements providing an opportunity for creation of long – term snow property records essential for climate monitoring and data assimilation studies - especially in the Arctic region, where we face rapid environmental changes including reduction of snow/ice extent and, therefore, planetary albedo.

This manuscript examines from the diurnal convection cycle (CDC) to the interdecadal variability in the region of the Peruvian Altiplano (RAP). Currently, estimating precipitation using satellites is an alternative which can be used to study the spatio-temporal evolution of precipitation systems. Herein CPC data Morphing technique - CMORPH (Joyce et al, 2004) was used between 2002 and 2014 to analyze the CDC in RAP. The CMOPRH data were compared with rainfall data series measured by rain gauges of meteorological stations (EMS) in the RAP. The results indicate that the CDC shows high variability in the Titicaca Basin and is associated with patterns of lake breeze (day), land breeze (night) and mountain - valley circulation. The CDC starts at 1800 HL (local time) in the northern region of Lake Titicaca, lasting between 2 h and 6 h, and most of 2000 HL. The CDC over the dry surface (ST) of Titicaca Basin starts early at around 1200 HL, lasting 4 h to 7 h, and maximum at 1800 HL.

The increasing spatial and spectral resolution of hyperspectral imagers yields detailed spectroscopy measurements from both space-based and airborne platforms. Machine learning algorithms have achieved state-of-the-art material classification performance on benchmark hyperspectral data sets; however, these techniques often do not consider varying atmospheric conditions experienced in a real-world detection scenario. To reduce the impact of atmospheric effects in the at-sensor signal, atmospheric compensation must be performed. Radiative Transfer (RT) modeling can generate high-fidelity atmospheric estimates at detailed spectral resolutions, but is often too time-consuming for real-time detection scenarios. This research utilizes machine learning methods to perform dimension reduction on the transmittance, upwelling radiance, and downwelling radiance (TUD) data to create high accuracy atmospheric estimates with lower computational cost than RT modeling. The utility of this approach is investigated using the instrument line shape for the Mako long-wave infrared hyperspectral sensor. This study employs physics-based metrics and loss functions to identify promising dimension reduction techniques. As a result, TUD vectors can be produced in real-time allowing for atmospheric compensation across diverse remote sensing scenarios.

Among other applications, radar-rainfall (RR) and QPE (Quantitative Precipitation Estimation) based on radar reflectivity, dual polarization variables, and multi-sensor information, provide important information for land surface hydrology, such as flood forecasting. Therefore, we developed a flood alert system using rainfall-runoff model forced with RR and QPE, and tipping-bucket observations to forecast river water levels (using rating-curves). In this study, we used an hourly dataset from an S-Band dual-polarimetric radar with two tropical R(Z) relations based distrometer data, a polarimetric R(Z,ZDR) algorithm from the literature and a multi-sensor approach using radar, satellite and rain gauge. Two hydrological models were used and calibrated using observed discharge time-series. Although our previous studies indicated accurate RR-based simulations, in some cases floods were not detected when using catchment-lumped rainfall derived from multi-sensor QPE. In this study, we advance further in this subject using improved R(Z,ZDR) relations and QPE for the period of 2016-2017 and flood event-based rainfall-runoff calibration. Thus, we focused on the development (and timing) of floods in the Marrecas River can be complex and strongly related to storms spatiotemporal distribution. To explore this aspect, we also perform a first analysis in using RR in rainfall-runoff model with a nested catchment discretization.

The 24 September 2001 College Park, Maryland, tornado was a long-tracked and violent tornado that passed within a close range of two Doppler radars. It was the third in a series of three tornadoes associated with a supercell storm that developed in Stafford County, Virginia, and initiated 3 - 4 km southwest of College Park and dissipated near Columbia, Howard County. The supercell tracked approximately 120 km and lasted for about 126 minutes. This study presents a synoptic and mesoscale overview of favorable conditions and forcing mechanisms that resulted in the severe convective outbreak associated with the College Park tornado. Results show many critical elements of the tornadic event, including a negative-tilted upper-level trough over the Ohio Valley, a jet stream with moderate vertical shear, a warm, moist tongue of the air associated with strong southerly flow over south-central Maryland and Virginia, and significantly increased convective available potential energy during the late afternoon hours. Satellite imagery reveals banded convective morphology with high cloud tops associated with the supercell that produced the College Park tornado. Operational WSR-88D data exhibits a high reflectivity “debris ball” or tornadic debris signature (TDS) within the hook echo, the evolution of the parent storm from a supercell structure to a bow echo, and a tornado cyclone signature (TCS). Many of the mesoscale environmental features could be captured by contemporary numerical model analyses. This study concludes with a discussion of the effectiveness of the coordinated use of satellite and radar observations in the operational environment of nowcasting severe convection.

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. In this study, we investigated single coal-burning particles emitted when using different fire-ignition techniques (top-lit up-draft vs 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. Individual particles (~700) were investigated using a scanning electron microscope equipped with energy-dispersive X-ray spectroscopy (EDS). 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. EDS analysis showed that 70% of the dark spherical organic particles had higher (~60%) relative oxygen content than in the bright SOPs. We quantified the morphology of spherical organic particles and classified them into four categories: ~50% are bare single particles; ~35% particles are aggregated and form diffusion accretion chains; 10% have inclusions, and 5% are deformed due to impaction on filter material during sampling. We conclude that there are two distinct kinds of coal burning spherical organic particles and that dark SOPs are less volatile than bright SOPs. We 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.

The assimilation of Soil Moisture and Ocean Salinity (SMOS) data into the ECMWF (European Centre for Medium Range Weather Forecasts) H-TESSEL (Hydrology revised - Tiled ECMWF Scheme for Surface Exchanges over Land) model is presented. SMOS soil moisture (SM) estimates have been produced specifically by training a neural network with SMOS brightness temperatures as input and H-TESSEL model SM simulations as reference. This can help the assimilation of SMOS information in several ways: (1) the neural network soil moisture (NNSM) data have a similar climatology to the model, (2) no global bias is present with respect to the model even if regional differences can exist. Experiments performing joint data assimilation (DA) of NNSM, 2 metre air temperature and relative humidity or NNSM-only DA are discussed. The resulting SM was evaluated against a large number of in situ measurements of SM obtaining similar results to those of the model with no assimilation, even if significant differences were found from site to site. In addition, atmospheric forecasts initialized with H-TESSEL runs (without DA) or with the analysed SM were compared to measure of the impact of the satellite information. Although, NNSM DA has an overall neutral impact in the forecast in the Tropics, a significant positive impact was found in other areas and periods, especially in regions with limited in situ information. The joint NNSM, T2m and RH2m DA improves the forecast for all the seasons in the Southern Hemisphere. The impact is mostly due to T2m and RH2m, but SMOS NN DA alone also improves the forecast in July- September. In the Northern Hemisphere, the joint NNSM, T2m and RH2m DA improves the forecast in April-September, while NNSM alone has a significant positive effect in July-September. Furthermore, forecasting skill maps show that SMOS NNSM improves the forecast in North America and in Northern Asia for up to 72 hours lead time.

A convenient practical model for accurately estimating the total entropy (ΣSi) of atmospheric gases based on physical action is proposed. This realistic approach is fully consistent with statistical mechanics, but reinterprets its partition functions as measures of translational, rotational and vibrational action or quantum states, to estimate the entropy. With all kinds of molecular action expressed as logarithmic functions, the total heat required for warming a chemical system from 0o K (ΣSiT) to a given temperature and pressure can be computed, yielding results identical with published experimental third law values of entropy. All thermodynamic properties of gases including entropy, enthalpy, Gibbs energy and Helmholtz energy are directly estimated using simple algorithms based on simple molecular and physical properties, without resource to tables of standard values; both free energies are measures of quantum field states and of minimal statistical degeneracy, decreasing with temperature and declining density. We propose that this more realistic approach has heuristic value for thermodynamic computation of atmospheric profiles, based on steady state heat flows equilibrating with gravity. Potentially, this application of an action principle can provide better understanding of emergent properties of many natural or evolving complex systems, including modelling of predictions for global warming.

Stretching along the border of North Dakota and Minnesota, The Red River Valley (RRV) of the North has the highest frequency of reported blizzards within the contiguous United States. Despite the numerous impacts these events have, few systematic studies exist discussing the meteorological properties of blizzards. As a result, forecasting these events and lesser blowing snow events is an ongoing forecast challenge. This study presents a climatology of atmospheric patterns associated with RRV blizzards for the winter seasons of 1979-1980 to 2017-2018. Patterns were identified using subjective and objective techniques using meteorological fields from the North American Regional Reanalysis (NARR). The RRV experiences on average, 2.6 events per year. Blizzard frequency is bimodal with peaks occurring in December and March. The events can largely be typed into four meteorological categories dependent on the forcing that drives the blizzard: Alberta Clippers, Arctic Fronts, Colorado Lows, and Hybrids. Objective classification of these blizzards using a competitive neural network known as the Self-Organizing Map (SOM) demonstrates that gross segregation of the events can be achieved with a small (8-class) map. This implies that objective analysis techniques can be used to identify these events in weather and climate model output that may aid future forecasting and risk assessment projects.

Relationships between drought and fire danger indices are examined to 1) incorporate fire risk information into the National Integrated Drought Information System California-Nevada Drought Early Warning System and 2) provide a baseline analysis for application of drought indices into a fire risk management framework. We analyzed four drought indices that incorporate precipitation and evaporative demand (E₀) and three fire indices that reflect fuel moisture and potential fire intensity. Seasonally averaged fire danger indices were most strongly correlated to multi-scalar drought indices that use E₀ (the Evaporative Demand Drought Index [EDDI] and Standardized Precipitation Evapotranspiration Index [SPEI]) at approximately annual time scales that reflect buildup of antecedent drought conditions. Results indicate that EDDI and SPEI can inform seasonal fire potential outlooks at the beginning of summer. An E₀ decomposition case study of conditions prior to the Tubbs Fire in Northern California indicate high E₀ (97th percentile) driven predominantly by low humidity signaled increased fire potential several days before the start of the fire. Initial use of EDDI by fire management groups during summer and fall 2018 highlights several value-added applications, including seasonal fire potential outlooks, funding fire severity level requests, and assessing set-up conditions prior to large, explosive fire cases.

We evaluate the complexity of Colombian climate from extreme behavior of gauge temperature and precipitation, using the the novel Tsallis' non-extensive entropy principle based on physical information through the q-index. We find the spatial structure of non additive universal categories (q-index) and compare with some complex systems with the potential to have some degree of dynamical affinity. Our results evidence the great dynamical variability of regional climate expressed in the large range of values of $q$-index, and the high degree of non-extensitivity for both temperature and precipitation.

Relationships between drought and fire danger indices are examined to 1) incorporate fire risk information into the National Integrated Drought Information System California-Nevada Drought Early Warning System and 2) provide a baseline analysis for application of drought indices into a fire risk management framework. We analyzed four drought indices that incorporate precipitation and evaporative demand (E₀) and three fire indices that reflect fuel moisture and potential fire intensity. Seasonally averaged fire danger indices were most strongly correlated to multi-scalar drought indices that use E₀ (the Evaporative Demand Drought Index [EDDI] and Standardized Precipitation Evapotranspiration Index [SPEI]) at approximately annual time scales that reflect buildup of antecedent drought conditions. Results indicate that EDDI and SPEI can inform seasonal fire potential outlooks at the beginning of summer. An E₀ decomposition case study of conditions prior to the Tubbs Fire in Northern California indicate high E₀ (97th percentile) driven predominantly by low humidity signaled increased fire potential several days before the start of the fire. Initial use of EDDI by fire management groups during summer and fall 2018 highlights several value-added applications, including seasonal fire potential outlooks, funding fire severity level requests, and assessing set-up conditions prior to large, explosive fire cases.

Upper Blue Nile basin (UBNB) is the water tower of Ethiopia and downstream countries. It contributes significant moistures to the surrounding atmosphere. However, the contribution of the moisture from the basin to the precipitation in the area is not well documented. Therefore, this paper is aimed at seasonal variation of upper Blue Nile basin moisture budget and the global moistures in the role of temporal and spatial precipitation variability. To this end, we used European Centre for Medium-range Weather Forecast (ECMWF) data from 1979-2017. The UBNB moisture contributed precipitation in the central parts of the study area during the summer season, while in spring; it contributed in southern part of the study area. Northwest part of the study area got precipitation from the basin moistures during autumn season. The recycling ratios for four seasons (summer, autumn, spring and winter) were 9.70%, 16.33%, 19.01%, and 35.30% respectively. The maximum amount of precipitation is extracted from the local moistures during winter season. The annual average value of recycling ratio was found 20.11%. Hence, we concluded that UBNB moisture budget had lesser contribution of precipitation over the study area. It rather contributed a significant precipitation to the neighboring countries such as Egypt and Sudan. Further studies on moisture budget are required to explain this phenomenon in the context of Ethiopia.

Upper Blue Nile basin (UBNB) is the water tower of Ethiopia and downstream countries. It contributes significant moistures to the surrounding atmosphere. However, the contribution of the moisture from the basin to the precipitation in the area is not well documented. Therefore, this paper is aimed at seasonal variation of upper Blue Nile basin moisture budget and the global moistures in the role of temporal and spatial precipitation variability. To this end, we used European Centre for Medium-range Weather Forecast (ECMWF) data from 1979-2017. The UBNB moisture contributed precipitation in the central parts of the study area during the summer season, while in spring; it contributed in southern part of the study area. Northwest part of the study area got precipitation from the basin moistures during autumn season. The recycling ratios for four seasons (summer, autumn, spring and winter) were 9.70%, 16.33%, 19.01%, and 35.30% respectively. The maximum amount of precipitation is extracted from the local moistures during winter season. The annual average value of recycling ratio was found 20.11%. Hence, we concluded that UBNB moisture budget had lesser contribution of precipitation over the study area. It rather contributed a significant precipitation to the neighboring countries such as Egypt and Sudan. Further studies on moisture budget are required to explain this phenomenon in the context of Ethiopia.

The Convective Rainfall Rate from Cloud Physical Properties (CRPh) for Meteosat Second Generation Satellites is a day-only precipitation algorithm developed at the Spanish Meteorological Agency (AEMET) for EUMETSAT’ Satellite Application Facility in support to Nowcasting and Very Short Range Forecasting (NWC SAF). It is therefore mainly intended to provide input for monitoring and near-real-time forecasts for the next few hours. This paper critically discusses the theoretical basis of the algorithm with special emphasis in the empirical values and assumptions in the microphysics of precipitation and compares the performances of the CRPh with its antecessor, the Convective Rainfall Rate algorithm (CRR), using an object-based method. The analyses show that AEMET’s CRPh is physically consistent and that outperforms the CRR. The applicability of the algorithm for nowcasting and the challenges to evolve the product to an all-day algorithm are also presented.

The main objective of this study is to search better prediction result of rainy seasonal rainfall (15 June-15 August). A correlation between rainfall of Bengali rainy seasons at Rangpur, Dhaka, Barisal and Sylhet and global sea surface temperature (SST) of different areas of the world was studied by using the both data of 1975- 2008 years with the help of the Climate Predictability Tool (CPT) to find more positive correlated SST with observed rainfall and use as predictor for giving the prediction of the year 2009. Using SST of one month before rainy season as predictor, the positive deviation of predicted rainfall from observed rainfall was 1.34 mm/day at Sylhet and 0.9 mm/day at Dhaka. The negative deviation of mean rainfall was 1.16 mm/day at Rangpur and 1.10 mm/day at Barisal. Again, using of starting one month SST of rainy season as predictor, positive deviation of predicted rainfall from observed rainfall was 4.03 mm/day at Sylhet. The positive deviation of daily mean rainfall was found 6.58 mm/day at Dhaka and 6.23 mm/day over southern Bangladesh. The study reveals that sea surface temperature (SST) of one month before rainy season was better predictor than SST of starting month of rainy season.

Electric and/or magnetic fields are generated by stationary charges, uniformly moving charges and accelerating charges. These field components are described in the literature as static fields, velocity fields (or generalized Coulomb field) and radiation fields (or acceleration fields), respectively. In the literature, the electromagnetic fields generated by lightning return strokes are presented using the field components associated with short dipoles and in this description the one to one association of the electromagnetic field terms with the physical process that gives rise to them is lost. In this paper, we will derive expressions for the electromagnetic fields using field equations associated with accelerating (and moving) charges and separate the resulting fields into static, velocity and radiation fields. The results illustrates how the radiation fields emanating from the lightning channel give rise to field terms varying as inverse of distance and distance squired, the velocity fields generating field terms varying as inverse of distance squired and the static fields generating fields components varying as inverse of distance squired and distance cube. These field components depend explicitly on the speed of propagation of the current pulse. However, the total field does not depend explicitly on the speed of propagation of the current pulse. It is shown that these field components can be combined to generate the field components pertinent to the dipole technique. However, in this conversion process the connection of the field components to the physical process taking place at the source that generate these fields (i.e. static charges, uniformly moving charges and accelerating charges) is lost.    

We propose a new technique to prepare statistically-robust benchmarking data for evaluating chemical transport model meteorology and air quality parameters within the urban boundary layer. The approach employs atmospheric class-typing, using nocturnal radon measurements to assign atmospheric mixing classes, and can be applied temporally (across the diurnal cycle), or spatially (to create angular distributions of pollutants as a top-down constraint on emissions inventories). In this study only a short (<1-month) campaign is used, but grouping of the relative mixing classes based on nocturnal mean radon concentrations can be adjusted according to dataset length (i.e., number of days per category), or desired range of within-class variability. Calculating hourly distributions of observed and simulated values across diurnal composites of each class-type helps to: (i) bridge the gap between scales of simulation and observation, (ii) represent the variability associated with spatial and temporal heterogeneity of sources and meteorology without being confused by it, and (iii) provide an objective way to group results over whole diurnal cycles that separates ‘natural complicating factors’ (synoptic non-stationarity, rainfall, mesoscale motions, extreme stability, etc.) from problems related to parameterizations, or between-model differences. We demonstrate the utility of this technique using output from a suite of seven contemporary regional forecast and chemical transport models. Meteorological model skill varied across the diurnal cycle for all models, with an additional dependence on the atmospheric mixing class that varied between models. From an air quality perspective, model skill regarding the duration and magnitude of morning and evening “rush hour” pollution events varied strongly as a function of mixing class. Model skill was typically the lowest when public exposure would have been the highest, which has important implications for assessing potential health risks in new and rapidly evolving urban regions, and also for prioritizing the areas of model improvement for future applications.

Earth atmosphere is almost opaque in the infrared: about 374 W/m² is absorbed by the atmosphere out of 396 W/m² surface upward longwave radiation, and only about 22 W/m² leaves the system unabsorbed in the atmospheric window. This makes rise to the idea to approximate the annual global mean energy flow system from a simple idealized greenhouse model, where the surface is surrounded by a single-layer shortwave (SW) transparent, longwave (LW) opaque, non-turbulent atmosphere. The energy flows in this geometry can be described by elementary arithmetic relationships. Starting from this model, the realistic Earth’s atmosphere can be achieved by introducing partial atmospheric SW opacity, partial atmospheric LW transparency and turbulent fluxes during the course of the deduction. The resulted global mean energy flow system is then compared to several data sets such as satellite observations from the CERES mission; estimates using direct surface observations and climate models; global energy and water cycle assessments; and independent detailed clear-sky radiative transfer computations. We find that the deduction from this idealized model approximates the real values in Earth energy budget with reasonable accuracy: the deduced fluxes and the observed ones are consistent within the acknowledged error of observations; while fundamental features of the initial geometry like special ratios and definite relationships between the fluxes are preserved.

This paper contributes to the understating of tornadoes in South Africa using case study analysis. In South Africa tornadoes are the recurring phenomenon (the climatology) but so far they have received less attention. Damages from storms itself (tornadoes inclusive) are significant in South Africa relative to other weather-related disasters for example floods, heat waves, and droughts. For their understanding, a case study approach was adopted in the current study. Data were in courtesy of the following, National Oceanic and Atmospheric Administration (NOAA), National Centers for Environmental Predictions (NCEP), Eumetsat Germany, and South African Weather Service (SAWS). The aim of the study was to provide an overview of the occurrence of tornadoes in South Africa using a Klerksdorp tornado, which occurred on March 4, 2007, Northwest Province in South Africa. From the case study analysis, the tornado was associated with the cold front and cut-off low (both are extratropical circulation) which were the dominant weather systems of the day. Therefore we conclude that, a case study approach may be the best way to study events of these nature for a more informed decision, for example, issuing an early warning system.  In future, case studies, for example, involving interaction between extratropical and tropical circulation will also be an interesting study.

This paper contributes to the understating of tornadoes in South Africa using case study analysis. In South Africa tornadoes are the recurring phenomenon (the climatology) but they have received less attention. Damages from storms itself (tornadoes inclusive) are significant in South Africa relative to other weather-related disasters for example floods, heat waves, and droughts. Case study approach was adopted in the current study. Data were in the courtesy of the following, National Oceanic and Administration (NOOA), National Centers for Environmental Protections (NCEP), Eumetsat, and South African weather Service. The aim of the study was to provide an overview of the occurrence of tornadoes in South Africa using a case study. From the case study analysis, the tornadoes at Klerksdorp on March 4, 2007, was associated with the cold frontal systems and the cut-ff low (extratropical circulation) which were the dominant weather systems of the day. Case study approach may be the best way to study events of these nature for a more informed decision, for example, issuing an early warning system.  Case studies, for example, involving interaction between extratropical and tropical circulation may be also more informative.

Understanding the distribution in drought and floods plays an important role in disaster risk management. The present study aims to explore the trends in the standardized precipitation index and extreme precipitation days in China, as well as to estimate the economic losses they cause. We found that in the Northeast China, northern of North China and northeast of Northwest China were severely affected by drought disasters (average damaged areas were 6.44 million hectares) and the most severe drought trend was located in West China. However, in the north of East China and Central China, the northeastern of the Southwest China was severely affected by flood disasters (average damaged areas were 3.97 million hectares) and the extreme precipitation trend is increasing in the northeastern of the Southwest China. In the Yangtze River basin, there were increasing trends in terms of drought and extreme precipitation, especially in the northeastern of the Southwest China, where accompanied by severe disaster losses. By combining the trends in drought and extreme precipitation days with the distribution of damaged areas, we found that the increasing trend in droughts shifted gradually from north to south, especially in the Southwest China, and the increasing trend in extreme precipitation gradually shifted from south to north.

Thunderstorm electrification has been studied for hundreds of years. Several mechanisms have been proposed to elucidate the electrification, including convective charging, inductive precipitation charging, and ice crystal-graupel collision charging. Field observations and model studies have demonstrated the vital roles that graupel and ice crystals play in the electrification, but the mechanism of the collision charging is still unclear. The fundamental essence of relative growth rate theory used for explaining the tripole charge structure in thunderclouds also needs a further exploration. We analyze the processes of ice crystal-graupel collision charging from charge migration inside hydrometeors to charge separation between two hydrometeors. The driving effects of temperature gradient and chemical potential gradient in charge migration are clarified, as well as the applicability of the relative growth rate theory, thermoelectric effect and surface tension gradient in different humidities. Based on the understanding from these electrification mechanisms, we propose that the essence of charge separation is driven by non-thermal equilibrium, and future studies on thunderstorm electrification should focus on the dynamical non-thermal equilibrium of cloud particles.

Analysis of the hourly Canadian Prairie data for the past 60 years has transformed our quantitative understanding of land-atmosphere-cloud coupling. The key reason is that trained observers made hourly estimates of opaque cloud fraction that obscures the sun, moon or stars, following the same protocol for 60 years at all stations. These 24 daily estimates of opaque cloud data are of sufficient quality that they can be calibrated against Baseline Surface Radiation Network data to give the climatology of the daily short-wave, longwave and total cloud forcing (SWCF, LWCF and CF). This key radiative forcing has not been available previously for climate datasets. Net cloud radiative forcing reverses sign from negative in the warm season to positive in the cold season, when reflective snow reduces the negative SWCF below the positive LWCF. This in turn leads to a large climate discontinuity with snow cover, with a systematic cooling of 10°C or more with snow cover. In addition, snow cover transforms the coupling between cloud cover and the diurnal range of temperature. In the warm season, maximum temperature increases with decreasing cloud, while minimum temperature barely changes; while in the cold season with snow cover, maximum temperature decreases with decreasing cloud and minimum temperature decreases even more. In the warm season, the diurnal ranges of temperature, relative humidity, equivalent potential temperature and the pressure height of the lifting condensation level are all tightly coupled to opaque cloud cover. Given over 600 station-years of hourly data, we are able to extract, perhaps for the first time, the coupling between cloud forcing and the warm season imbalance of the diurnal cycle; which changes monotonically from a warming and drying under clear skies to a cooling and moistening under cloudy skies with precipitation. Because we have the daily cloud radiative forci, which is large, we are able to show that the memory of water storage anomalies, from precipitation and the snowpack, goes back many months. The spring climatology shows the memory of snowfall back through the entire winter, and the memory in summer goes back to the months of snowmelt. Lagged precipitation anomalies modify the thermodynamic coupling of the diurnal cycle to the cloud forcing, and shift the diurnal cycle of mixing ratio which has a double peak. The seasonal extraction of the surface total water storage is a large damping of the interannual variability of precipitation anomalies in the growing season. The large land-use change from summer fallow to intensive cropping, which peaked in the early 1990s, has led to a coupled climate response that has cooled and moistened the growing season, lowering cloud-base, increasing equivalent potential temperature, and increasing precipitation. We show a simplified energy balance of the Prairies during the growing season and its dependence on reflective cloud.

Mixing layer height (MLH) is a crucial parameter for air quality modelling that is still not routinely measured. Common methods for MLH determination use atmospheric profiles recorded by radiosonde but they suffer from course temporal resolution since balloon launching is only twice a day. Recently cheap ceilometers are gaining popularity in the retrieval of MLH diurnal evolution based on aerosol profiles. This study presents a comparison of a proprietary (Jenoptik) and a free available (STRAT) algorithms to retrieve MLH diurnal cycle. The comparison is accomplished in summer season over urban area and radiosonde data is used to estimate MLHs according to parcel, lapse rate, and Richardson methods (the last algorithm is used as a reference in the study) in addition. It was found that STRAT and Jenoptik give lower MLH than radiosonde with an underestimation of about 150m and 650m respectively. Additionally, STRAT showed reasonable performance in tracking of MLH diurnal evolution. Daily MLH maximum of about 2000m was found in the late afternoon (18-19 LT). In contrast, Jenoptik algorithm showed more weaknesses, mainly attributed to its real-time operation and independent processing of a single profile. At night and during morning transition period, both lidar-based methods showed difficulties as MLH was often in the ceilometer’s incomplete overlapping zone so residual or advected aerosol layers aloft were misleadingly reported as mixing layer (ML).

For water-stressed regions like Iran improving the effectiveness and productivity of agricultural water-use is of utmost importance due to climate change and unsustainable demands. Therefore, a hydro-economic model has been developed here for the Zarrine River Basin with the central concept of that demands are value-sensitive functions, where quantities of water-uses at different locations and times have a changeable economic benefits. To do this, the potential crop yields and the surface and groundwater resources, especially Boukan Dam inflow are simulated using the hydrologic model, SWAT, based on predicted climatic scenarios i.e. quantile mapping-downscaled projections. Then, to allocate the agricultural water based on the agro- economic crop water productivity (AEWP) of crops, a basin-wide water management tool, MODSIM, is customized. Next, a simulation- optimization model has been developed using a coupled CSPSO-MODSIM, to optimize the total AEWP, considering climatic impact and crop pattern scenarios, for 2020-2038, 2050-2068 and 2080-2098 periods. Finally, the optimum crop pattern and crop water irrigation depths are presented for different RCPs and periods. The results indicated that this approach will improve considerably the AEWPs and decrease the agricultural water-use up to 40%. Thus, this integrated model is able to support water authorities and other stakeholder in a water-scarce basin, as is the study area.

This study presents a novel approach for the early detection of developing thunderstorms. To date, methods for the detection of developing thunderstorms usually rely on accurate atmospheric motion vectors (AMVs) for the estimation of cooling rate of convective clouds, which corresponds to the updraft strength of the cloud objects. In this study we present a method for the estimation of the updraft strength that does not rely on AMVs. The updraft strength is derived directly from the satellite observations in the SEVIRI water vapor channels. For this purpose the absolute value of the vector product of two vectors derived from the observed radiances emitted by the water vapor channels in 2 subsequent satellite images is calculated. The x and y components of the vectors consists of the temporal change of the spatial brightness temperature gradients in x- and y-direction, the z component consists of the temporal variation of the brightness temperature at pixel level. The absolute value of the vector product is referred to as normalized updraft strength (NUS). The performance of the method has been investigated for 2 summer periods in 2016 and 2017 by validation with lightning data. Values of the Critical Success Index (CSI) of about 66 % for the 2016 period and 62 % for the 2017 period demonstrate the good performance of the method. The POD values for the experiments with the highest CSI values are 90.3 for 2016 period and 87.6 for the 2017 period, respectively The corresponding FAR values are 28.2% (2016) and 32.5 % (2017), respectively. In summary, the method has the potential to reduce the forecast lead time significantly, and can be quite useful in regions without a well maintained radar network.

Data on the atmospheric gas concentrations can be received with high efficiency and on a large spatial scale only from remote laser sounding instruments. The remote laser techniques with the use of lidars are widely used in the study of the atmosphere and control of its state. The aims of this work are the design and test in numerical and field experiments of a DIAL OPO lidar system based on KTA and KTP crystals for gas analysis of the atmosphere. Lidar measurements of atmospheric gases in the near/mid infrared region have been numerically simulated. The differential absorption lidar system based on optical parametric oscillators with nonlinear KTA and KTP crystals which allow laser radiation tuning both in the near and in the middle IR spectral region is described; it allows tuning laser radiation in the near/mid-IR wavelength regions. Lidar echo signals have been experimentally recorded in the 1.8–2.5 and 3–4 m wavelength ranges. The results of H2O and CO2 profile measurements along the surface sounding path are presented.

A forecasting scheme of the thermal updraft velocity based on a theoretical model and data collected from flights records at gliding competitions, is presented. The forecasting scheme was based on the hypothesis that there is linear relationship between the overheat function at ground surface and the temperature difference between soil and air. The proportionality factor of this relationship was determined experimentally using observations recorded during gliding flights. The results showed that based on this simple scheme forecasting thermal convection is possible at any geographical location.

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 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 500x500km 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 characterized by a distinct signature in satellite imagery.

Turbulent flows over and within forest canopies have recently been modeled with success using Large Eddy Simulations (LES). Validation exercises against experimental data suggest that models can be applied with a high degree of confidence for many applications, mechanical and physiological plant/atmosphere interaction analysis, seed or pollen dispersal, wildfire spread and firebrand transport, or investigation of causes of eddy-covariance technique bias. Long distances required for shear-induced turbulence to equilibrate, result in the widespread use of cyclic boundary conditions in LES atmospheric boundary layer studies. Vegetation drag dissipates air momentum in the atmosphere, but equilibrium is often achieved through compensatory momentum source, supplied by macro-scale pressure gradient forcing. Unfortunately, both classical Ekman balance or simple spatially-constant pressure gradient techniques for implementing this forcing have major drawbacks in the context of cyclic boundary conditions for the applications listed above. Among them, it is difficult to specify aspects of the mean velocity profile such as a specific desired wind velocity and direction at a reference height. In the present paper, we propose a new technique for capturing the effects of  a large-scale pressure gradient force (LSPGF) that can be used at stand scale and enables  simulation of realistic and specifiable wind fields. Several variants of this LSPGF are developed and analyzed here and validated against experimental data. Although this LSPGF technique is developed in the context of HIGRAD/FIRETEC wildfire simulations, LSPGF can be used for any LES wind modeling application aimed at generating detailed stand-scale wind fields with resolved turbulence and shear profiles consistent with vegetation structure in the boundary layer.

Improving air quality and reducing human exposure to unhealthy levels of airborne chemicals are important global missions, particularly in China. Satellite remote sensing offers a powerful tool to examine regional trends in NO2, thus providing a direct measure of key parameters that strongly affect surface air quality. To accurately resolve spatial gradients in NO2 concentration using satellite observations and thus understand local and regional aspects of air quality, a priori input data at sufficiently high spatial and temporal resolution to account for pixel-to-pixel variability in the characteristics of the land and atmosphere are required. In this paper, we adapt the Berkeley High Resolution product (BEHR v3.0A, v3.0B and v3.0C) and meteorological outputs from the Weather Research and Forecasting (WRF) model to describe column NO2 in southern China. The BEHR approach is particularly useful for places with large spatial variabilities and terrain height differences such as China. We retrieved tropospheric NO2 vertical column density (TVCD) within part of southern China, for four seasons of 2015, based upon satellite datasets from Ozone Monitoring Instrument (OMI). Retrieval results are validated by comparing with MAX-DOAS tropospheric column measurements conducted in Guangzhou. BEHR retrieval algorithms are more consistent with MAX-DOAS measurements than OMI-NASA retrieval, opening new windows into research questions that require high spatial resolution, for example retrieving NO2 vertical column and ground pollutant concentration in China and other countries.

Towards the investigation and further understanding of the development and propagation of Medicanes, this study explores the forecasting capability of WRF model in case of cyclone “Cleopatra” which affected with extreme rainfall and strong winds Sardinia and Calabria, Italy, in November 2013.  This cyclone was unusual in that it developed a warm core but did not fulfill its transformation into a tropical-like cyclone because its core did not expand high enough in the tropospnere. The ERA5 reanalysis dataset was dynamically downscaled from 31 km spatial horizontal resolution to 9 km using WRF model. The methodology consists of; firstly, an extensive physical parameterization schemes sensitivity test and consequently, a short-range ensemble forecasting implementation based on the highest statistical scored physics configuration. All simulation results were validated against surface observations and remote sensing products. Subsequently, the modeled cyclone trajectories are compared to satellite imagery derived from EUMETSAT-SEVIRI gridded data. The findings of the conducted analysis illustrate that ensemble average displays significant difference in performance compared to any of the deterministic runs individually, suggesting that ensemble forecasts will be beneficial in studies assessing cyclonic events in the Mediterranean region.

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 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 present a conceptually novel approach for the detection of MCs based on the use of deep convolutional neural networks (DCNNs). We demonstrate that DCNN model is capable of performing binary classification of 500x500km 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. 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 characterized by a distinct signature in satellite imagery.

Here we present an overview of speleothem δ18O records from South America, which mostly are available in the Speleothem Isotopes Synthesis and Analysis (SISAL_v1) database. South American tropical and subtropical δ18O time series are primarily interpreted as being driven by the amount effect and, consequently show the past history of the convection intensity of convergence zones such as the Intertropical Convergence Zone and the South America Monsoon System. We investigate past hydroclimate scenarios in South America related to the South American Monsoon System in three different time scales: Late Pleistocene, Holocene and the last two millennia. The precession driven insolation is the main driver of convective variability over the continent during the last 250 kyrs, including the Holocene period. However a dipole is observed between the west and east portions of the continent. Records located in the central region of Brazil appear to be weakly affected by insolation driven variability and more susceptible to the South Atlantic Convergence Zone. Cold episodic events in Northern Hemisphere increase the activity of the South American Monsoon System on all time scales, in turn increasing rainfall amounts in South America, as was documented during Heinrich events in the late Pleistocene and Bond events in the Holocene, as well as during the Little Ice Age.

Global wind observations are fundamental for studying weather and climate dynamics. Most wind measurements come from atmospheric motion vectors (AMVs) by tracking the displacement of cloud or water vapor features. These AMVs generally rely on thermal infrared (IR) techniques for their height assignments, which are subject to large uncertainties in the presence of weak or reversed vertical temperature gradients around the planetary boundary layer (PBL) and with tropopause folding. Stereo imaging can overcome the height assignment problem using geometric parallax for feature height determination. In this study we develop a stereo 3D-Wind algorithm to simultaneously retrieve AMV and height from geostationary (GEO) and low Earth orbit (LEO) satellite imagery and apply it to collocated Geostationary Operational Environmental Satellite (GOES) and Multi-angle Imaging SpectroRadiometer (MISR) imagery. The new algorithm improves AMV and height relative to products from GOES or MISR alone, with an estimated accuracy of <0.5 m/s in AMV and <200 m in height with 2.2 km sampling. The algorithm can be generalized to other LEO-GEO or GEO-GEO combinations for greater spatiotemporal coverage. The technique demonstrated with MISR and GOES has important implications for future high-quality AMV observations, for which a low-cost constellation of CubeSats can play a vital role.

This study focusses on identifying a set of representative future climate projections for the Upper Indus Basin (UIB). Although a large number of GCM’s predictor sets are nowadays available in the CMIP5 archive, the issue of their reliability for specific regions must still be confronted. This situation makes it imperative to sort out the most appropriate, single or small-ensemble set of GCMs for the assessment of climate change impacts in a region. Here a set of different approaches is adopted and applied for a step-wise shortlist and selection of appropriate climate models for the UIB under two RCPs: RCP 4.5 and RCP 8.5, based on, a) range of projected mean changes, b) range of projected extreme changes, and c) skill in reproducing the past climate. Furthermore, because of higher uncertainties in climate projection for high mountainous regions like the UIB, a wider range of future GCM climate projections is considered by using all possible future extreme scenarios (wet-warm, wet-cold, dry-warm, dry-cold). Based on this two-fold procedure, a limited number of climate models is pre-selected, out of which the final selection is done by assigning ranks to the weighted score for each of the mentioned selection criteria. The dynamically downscaled climate projections from the Coordinated Regional Downscaling Experiment (CORDEX) available for the top-ranked GCMs are further statistically downscaled (bias-corrected) over the UIB. The downscaled projections up to year 2100 indicate temperature increases ranging between 2.3 °C and 9.0 °C and precipitation changes that range, from a slight annual increase of 2.2% under the drier scenarios,  to as high as 15.9% for the wet scenarios.  Moreover, for all scenarios, the future precipitation will be more extreme, as the probability of wet days will decrease, while, at the same time, the precipitation intensities will increase. The spatial distribution of the downscaled predictors across the UIB also shows similar patterns for all scenarios, with a distinct precipitation decrease over the south-eastern parts of the basin, but an increase in the northeastern parts. These two features are particularly intense for the “Dry-Warm” and the “Median” scenarios over the late 21st century.

Nowadays a recognized need for accurate observations of atmospheric aerosols (AEs) and reactive gases (RGs) exists in the framework of regional, national and global near-surface networks based on permanent or mobile measurement stations. In this context, a paramount and not-trivial issue is related to the correct execution of continuous sampling of ambient air and its subsequent distribution to measurement analyzers hosted inside the stations. Sampling artifacts must be minimized for obtaining reliable pictures of ambient air composition. To respond to this need, a suite of novel “smart” and relatively low-cost systems for the continuous sampling of ambient air was developed in the framework of the Project I-AMICA (2012 – 2015, www.i-amica.eu). These systems were designed to execute AE and RG measurements according with WMO/GAW and ACTRIS recommendations and standard operation procedures. A particular attention was dedicated to the stabilization and control of the sampling flow rates and temperatures. The analysis of one full year of operations at the WMO/GAW regional station of Capo Granitola (GAW ID: CGR, Italy), allowed to conclude that these systems are effective in meeting the technical requirements for correct execution of AE and RG measurements

The impact of snow darkening by deposition of light absorbing aerosols (LAAs) on snow cover over the Himalaya-Tibetan-Plateau (HTP) and influence on the Asian monsoon are investigated using the NASA Goddard Earth Observing System Model Version 5 (GEOS-5). We find that during April-May-June, deposition of LAAs on snow leads to a reduction in surface albedo, initiating a sequence of feedback processes, starting with increased surface solar radiation, rapid snowmelt in HTP and warming of the surface and upper troposphere, followed by enhanced low-level southwesterlies and increased dust loading over the Himalayas-Indo-Gangetic Plain. The warming is amplified by increased dust aerosol heating, and subsequently amplified by latent heating from enhanced precipitation over the Himalaya foothills and northern India, via the Elevated Heat Pump (EHP) effect during June-July-August. The reduced snow cover in the HTP anchors the enhanced heating over the Tibetan Plateau and its southern slopes, in conjunction with an enhancement of the Tibetan Anticyclone, and the development of an anomalous Rossby wavetrain over East Asia, leading to weakening of the subtropical westerly jet, and northward displacement and intensification of the Mei-Yu rainbelt. Our results suggest that atmosphere-land heating by LAAs, particularly desert dust play a fundamental role in physical processes underpinning the snow-monsoon relationship proposed by Blandford more than a century ago.

Abstract: Characterization the 3-D structure of clouds is needed for a more complete understanding of the Earth's radiative and latent heat fluxes. Here we develop and explore a “ray casting” algorithm applied to the Multi-angle Imaging SpectroRadiometer (MISR) on board the Terra satellite, to reconstruct 3-D cloud volumes for observed clouds. The ray casting algorithm is first applied to geometrically simple synthetic clouds to show that, under the assumption of perfect, clear-conservative cloud masks, the reconstruction method yields overestimation whose magnitude depends on the cloud geometry and the resolution of the reconstruction grid relative to the image pixel resolution. The method is then applied to two select MISR scenes, fully accounting for MISR’s viewing geometry for reconstructions over the Earth’s ellipsoidal surface. The MISR Radiometric Camera-by-camera Cloud Masks at 1.1 km resolution and custom cloud masks at 275 m resolution independently derived from MISR RGB channels are used as input cloud masks. A wind correction method, termed “cloud spreading”, is devised and applied to the cloud masks to offset potential cloud movements over short time intervals (around 7 minutes at maximum) between the cameras. The MISR cloud top height product is used as a constraint to reduce the overestimation at the cloud top. The reconstruction results show that their uncertainty is significant when the wind correction is applied, and that they have more refined structures when the input cloud mask has a higher resolution. Recommendations for improving the presented cloud volume reconstructions as well as for future passive remote sensing satellite missions are discussed.

 This study shows the application of self-organizing maps (SOMs) to probabilistic forecasts of wind power generation and ramps in Japan. SOMs are applied to atmospheric variables obtained from atmospheric reanalysis over the region, thus deriving classified weather patterns (WPs). Probabilistic relationships are established between the synoptic-scale atmospheric variables over East Japan and the generation of regionally integrated wind power in East Japan. Medium-range probabilistic wind power predictions are derived by SOM, as analog ensembles based on the WPs of the multi-center ensemble forecasts. As this analog approach handles stochastic uncertainties effectively, probabilistic wind power forecasts are rapidly generated from a very large number of forecast ensembles. The use of a multi-model ensemble provides better results than a one-forecast model. The hybrid ensemble forecasts further improve the probabilistic predictability skill of wind power generation, as compared with non-hybrid methods. It is expected that long-term wind forecasts will provide better guidance to transmission grid operators. The advantage of this method is that it can include an interpretative analysis of meteorological factors for variations in renewable energy.

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 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 present a conceptually novel approach for the detection of MCs based on the use of deep convolutional neural networks (DCNNs). We demonstrate that DCNN model is capable of performing binary classification of 500x500km 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. 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 characterized by a distinct signature in satellite imagery.

The current study applied a new approach for the interpolation and regionalization of observed precipitation series to a smaller spatial scale (0.125° by 0.125° grid) across the Upper Indus Basin (UIB), with appropriate adjustments for the orographic effect and changes in glacier storage. The approach is evaluated and validated through reverse hydrology, guided by observed flows and available knowledge base. More specifically, the generated corrected precipitation data is validated by means of SWAT-modelled responses of the observed flows to the different input precipitation series (original and corrected ones). The results show that the SWAT- simulated flows using the corrected, regionalized precipitation series as input are much more in line with the observed flows than those using the uncorrected observed precipitation input for which significant underestimations are obtained.

It has been an outstanding challenge for global climate models to simulate and predict East Asia (EA) summer monsoon (EASM) rainfall. This study evaluates the dynamical hindcast skills with the newly developed Nanjing University of Information Science and Technology Earth System Model version 3.0 (NESM3.0). To improve the poor prediction of an earlier version of NESM3.0, we have modified convective parameterization schemes to suppress excessive deep convection and enhance insufficient shallow and stratiform clouds. The new version of NESM3.0 with modified parameterizations (MOD hereafter) yields significantly improved rainfall prediction in the northern and southern China but not over the Yangtze River Valley. The improved prediction is primarily attributed to the improvements in the predicted climatological summer mean rainfall and circulations, seasonal march of the subtropical rain belt, Nino 3.4 SST anomaly, and the rainfall anomalies associated with the development and decay of El Nino events. However, the MOD still has notable biases in the predicted leading mode of interannual variability of precipitation. The leading mode captures the dry (wet) anomalies over the South China Sea (northern EA) but misplaced precipitation anomalies over the Yangtze River Valley. The model can capture the interannual variation of the circulation indices very well, but the bias in the circulation-rainfall connection caused predicted rainfall errors. The results here suggest that over EA land regions, the skillful rainfall prediction relies on not only model’s capability in predicting better summer mean and seasonal march of rainfall and ENSO teleconnection with EASM, but also accurate prediction of the leading modes of interannual variability.

The present study aims to evaluate the capability of the TRMM-3B42-(V7) precipitation product to estimate appropriate precipitation rates in the Upper Indus basin (UIB) and the analysis of the dependency of the estimates’ accuracies on the time scale. To that avail statistical analyses and comparison of the TMPA- products with gauge measurements in the UIB are carried out. The dependency of the TMPA estimates’ quality on the time scale is analysed by comparisons of daily, monthly, seasonal and annual sums for the UIB. The results show considerable biases in the TMPA- (TRMM) precipitation estimates for the UIB, as well as high false alarms and miss ratios. The correlation of the TMPA- estimates with ground-based gauge data increases considerably and almost in a linear fashion with increasing temporal aggregation, i.e. time scale. The BIAS is mostly positive for the summer season, while for the winter season it is predominantly negative, thereby showing a slight over-estimation of the precipitation in summer and under-estimation in winter. The results of the study suggest that, in spite of these discrepancies between TMPA- estimates and gauge data, the use of the former in hydrological watershed modelling, endeavoured presently by the authors, may be a valuable alternative in data- scarce regions, like the UIB, but still must be taken with a grain of salt.

This study presents a pre-processing approach adopted for the radar reflectivity data assimilation and results of simulations with the Harmonie numerical weather prediction model. The method shows an improvement of precipitation prediction within the radar location area in both the rain rates and spatial pattern presentation. With the assimilation of radar data, the model simulates larger water content in the middle troposphere within the layer from 1 to 6 km, with major variations at 2.5–3 km; it also reproduces better the mesoscale belt and cell patterns of precipitation fields.

Passive microwave measurements have been available on satellites dating back to the 1970s on research satellites flown by the National Aeronautics and Space Administration (NASA). Since then, several other sensors have been flown to retrieve hydrological products for both operational weather applications (e.g., the Special Sensor Microwave/Imager–SSM/I; the Advanced Microwave Sounding Unit–AMSU) and climate applications (e.g., the Advanced Microwave Scanning Radiometer–AMSR; the Tropical Rainfall Measurement Mission Microwave Imager–TMI; the Global Precipitation Mission Microwave Imager–GMI). Here the focus is on measurements from the AMSU-A, AMSU-B and Microwave Humidity Sounder (MHS). These sensors have been in operation since 1998 with the launch of NOAA-15, and are also on board NOAA-16, -17, -18, -19 and the MetOp-A and -B satellites. A data set called the “Hydrological Bundle” is a Climate Data Record (CDR) that utilizes brightness temperatures from Fundamental CDRs to generate Thematic CDRs (TCDR). The TCDR’s include: Total Precipitable Water (TPW), Cloud Liquid Water (CLW), Sea-Ice concentration (SIC), Land surface temperature (LST), Land surface emissivity (LSE) for 23, 31, 50 GHz, rain rate (RR), snow cover (SC), ice water path (IWP), and snow water equivalent (SWE). The TCDR’s are shown to be in general good agreement with similar products from other sources such as the Global Precipitation Climatology Project (GPCP) and the Modern-Era Retrospective Analysis for Research and Applications (MERRA-2). Because of the careful intercalibration of the FCDR’s, little bias is found among the different TCDR’s produced from individual NOAA and MetOp satellites, except for normal diurnal cycle differences.

The data assimilation method to improve sea fog forecast over the Yellow Sea is usually three-dimensional variational assimilation (3DVAR), whereas ensemble Kalman filter (EnKF) has not yet been applied on this weather phenomenon. In this paper, two sea fog cases over the Yellow sea, one spread widely and the other spread narrowly along the coastal area, are studied in detail by a series of numerical experiments with 3DVAR and EnKF based on the Grid-point Statistical Interpolation (GSI) system and the Weather Research and Forecasting (WRF) model. The results show that the assimilation effect of EnKF outperforms that of 3DVAR: for the widespread-fog case, the probability of detection and equitable threat scores of the forecasted sea fog area get improved respectively by ~57.9% and ~55.5%; the sea fog of the other case completely mis-forecasted by 3DVAR is produced successfully by EnKF. These improvements of EnKF relative to 3DVAR are benefited from its flow-dependent background error, resulting in more realistic depiction of sea surface wind for the widespread-fog case and better moisture distribution for the other case in the initial conditions. More importantly, the correlation between temperature and humidity in the background error of EnKF plays a vital role in the response of moisture to the assimilation of temperature, which leads to a great improvement on the initial moisture conditions for sea fog forecast.

We review the state of dynamical downscaling with scale-constrained regional and global models. The methodology, in particular spectral nudging, has become a routine and well-researched tool for hindcasting climatologies of sub-synoptic atmospheric disturbances in coastal regions. At present, the spectrum of applications is expanding to other phenomena, but also to ocean dynamics and to extended forecasting. Also new diagnostic challenges are appearing such as spatial characteristics of small-scale phenomena such as Low Level Jets.

The predictability of wind information in a given location is essential for the evaluation of a wind power project. Predicting wind speed accurately improves the planning of wind power generation, reducing costs and improving the use of resources. This paper seeks to predict the mean hourly wind speed in anemometric towers (at a height of 50 meters) at two locations: a coastal region and one with complex terrain characteristics. To this end, the Holt-Winters (HW), Artificial Neural Networks (ANN) and Hybrid time-series models were used. Observational data evaluated by the Modern-Era Retrospective analysis for Research and Applications-Version 2 (MERRA-2) reanalysis at the same height of the towers. The results show that the hybrid model had a better performance in relation to the others, including when compared to the evaluation with MERRA-2. For example, in terms of statistical residuals, RMSE and MAE were 0.91 and 0.62 m/s, respectively. As such, the hybrid models are a good method to forecast wind speed data for wind generation.

Large water bodies such as the Laurentian Great Lakes have significant influences on local and regional climate through their unique physical features. Due to the coarse spatial resolution of general circulation models (GCMs), the Great Lakes are geometrically ignored in most GCMs. Thus, the dynamical downscaling technique serves as a necessary and feasible solution to bridge the gap. The Weather Research and Forecasting model (WRF) with an updated lake scheme is employed to downscale from a GCM, GFDL-CM3. The WRF-Lake’s performance is evaluated against observations, the GCM, as well as 23 other GCMs. Results show that the coupled air-lake model, with a fine spatial resolution and realistic lake bathymetries, reproduced a more reasonable spatiotemporal climatology than GCMs, as well as the lake-induced characteristics that were missed in GCMs. With lakes present, the seasonal variability of air temperature was reduced in WRF-Lake relative to GFDL-CM3, especially in summer. A reduced air temperature trend, about 4.5 °C/100 year in the 21st century, was projected in WRF-Lake. The seasonal evolutions of lake surface temperature and lake ice coverage were well captured by the lake model. The lake surface temperature was projected to be warming by 3.5-4 °C and the lake ice diminishing by 58.9% - 86%. Those results brought by the WRF-Lake model suggest that a fine resolution of the topography and the incorporation of the lake-atmosphere interactions are crucial to improve the understanding of the climate and climate change in the Great Lakes region.

A diagnostic method is presented for analysing the large-scale behaviour of the Met Office Unified Model, which is a comprehensive atmospheric model used for weather and climate prediction. Outside the boundary layer, on scales larger than the radius of deformation, semigeostrophic theory will give an accurate approximation to the model evolution. In particular, the ageostrophic circulation required to maintain geostrophic and hydrostatic balance against prescribed forcing and a rate of change of the geostrophic pressure can be calculated. In the tropics the balance condition degenerates to the weak temperature gradient approximation. Within the boundary layer the semigeostriptic approximation has to be used because friction and rotation are equally important. Assuming the calculated pressure tendency and ageotriptic circulation match the observed model behaviour, the influence of the large-scale state and the nature of the forcing on the model response can be deduced in a straightforward way. This process is illustrated by comparing predictions of the ageotriptic circulation from the theory and the model. It is then used to show that the effects of latent heat release can be included by modifying the static stability, and to show the effect of an idealised tropical heat source on the subtropical jet. Finally the response of the ageotriptic flow to boundary layer heating in the tropics is demonstrated. These illustrations show that the model behaviour on large scales conforms with theoretical expectations, so that the results of the diagnostic can be used to aid the development of further improvements to the model.

El Niño Southern Oscillation (ENSO) is an important of mode of climatic variability that through altering climate patterns exerts a discernible impact on ecosystems and society. For this reason, ENSO has attracted much interest in the climate and health science community with many analysts investigating ENSO health links through considering the degree of dependency of the incidence of a range of climate diseases on the occurrence of El Niño events. Because of the mounting interest in the relationship between ENSO as a major mode of climatic variability and health this paper presents an overview of the basic characteristics of the ENSO phenomenon and its climate impacts, discusses the use of ENSO indices in climate and health research and outlines the present understanding of ENSO health associations. Also touched upon are ENSO-based seasonal health forecasting and the possible impacts of climate change on ENSO and the implications this holds for future assessments of ENSO health associations. The review concludes that there is still some way to go before a thorough understanding of the association between ENSO and health is achieved with a need to move beyond analyses undertaken through a purely statistical lens with due acknowledgement that ENSO as a complex non-canonical phenomenon and that simple ENSO health associations should not be expected.

The Absorbing/scattering nature of aerosols affects the total radiative forcing and this absorption to total extinction ratio is quantified by single scattering albedo (SSA). Effect of SSA in the Ultraviolet (UV) irradiance is less studied and limited measurements are available. SSA retrieved at Athens, Greece during 2009-2014 from Ultraviolet Multifilter Radiometer (UVMFR) at 332 and 368 nm, were used to calculate incoming UV irradiance, alongside with ones from AERONET at visible wavelengths, from OMI satellite at 342.5 nm and from AEROCOM climatological database at 300 nm. UVA and UVB irradiances were estimated using a Radiative Transfer Model and we found that relative differences could be as high as 20%, while average relative differences varied from 2% to 8.7 % for the whole experimental period. Both UVA and UVB drop by a rate of ~12% for 0.05 aerosol absorption optical depth compared to ones estimated using SSA at visible range. Brewer irradiance measurements at 324nm were used to validate simulated irradiances and a better agreement was found when UVMFR SSAs were used with an average difference of 0.86%, while when using visible or climatological input, relative differences were estimated +4.91 and +4.15% accordingly.

Land Surface Temperature (LST) is a basic parameter in energy exchange between the land and atmosphere and is frequently used in many sciences such as climatology, hydrology, agriculture, ecology, etc. LST time series data have usually deficient, missing and unacceptable data caused by the presence of clouds in images, presence of dust in atmosphere and sensor failure. In this study, Singular Spectrum Analysis (SSA) algorithm was used to resolve the problem of missing and outlier data caused by cloud cover. The region studied in the present research included an image frame of MODIS with horizontal number 22 and vertical number 05 (h22v05). This image involved a large part of Iran and Turkmenistan and Caspian Sea. In this study, MODIS LST sensor (MOD11A1) was used during 2015 with 1×1 Km spatial resolution and day/night LST data (daily temporal resolution). The results of the data quality showed that cloud cover caused 36.37% of missing data in the studied time series with 730 day/night LST images. Further, the results of SSA algorithm in reconstruction of LST images indicated the Root Mean Square Error (RMSE) of 2.95 K between the original and reconstructed data in LST time series in the study region. In general, the findings showed that SSA algorithm using spatio-temporal interpolation in LST time series can be effectively used to resolve the problem of missing data caused by cloud cover.

With the continual global increase of PV power systems and the inherent weather-induced volatility of their power output, understanding the underlying variability of solar radiation in both space and time is important for the planning and reliable operation of future power grids. This paper concisely reviews recent advances in the characterization of irradiance variability with an emphasis on small spatial and temporal scales (respectively less than about 10\,km and 1\,min), for which comprehensive data sets are time consuming and expensive to collect, and although needed have not been available until recently. Special attention is given to studies dealing with the quantification of variability using such unique data, the analysis and modeling of spatial smoothing, and the evaluation of temporal averaging.

Maritime Continent (MC) positions in between Asian and Australian summer monsoons zone. Its complex topography and shallow seas around it is a major challenge for the climate researchers to model and understand it. Monsoon in this area is affected by inter-scale ocean-atmospheric interactions like El-Niño Southern Oscillation (ENSO), Indian Ocean Dipole (IOD) and Madden-Julian Oscillation. Monsoon rainfall in MC (especially in Indonesia and Malaysia) profoundly exhibits its variability dependency on ocean-atmospheric phenomena in this region. This monsoon shift often introduces to dreadful events like biomass burning (BB) in Southeast Asia (SEA) which sometimes leads to severe trans-boundary haze pollution. In this study, the episode of BB in 2015 of SEA is highlighted and discussed. Observational satellite datasets are tested by performing simulations with numerical weather prediction (NWP) model using WRF-ARW (Advanced research WRF). Observed and model datasets are compared to study the sea surface temperature (SST) and precipitation (rainfall) anomalies influenced by ENSO, IOD and MJO. Correlations have been recognised which explains the delayed rainfall of regular monsoon in MC due to the influence of ENSO, IOD and MJO during 2015 BB episode, eventually leading to intensification of fire and severe haze.

Several long-term monitoring of aerosol datasets from the Moderate Resolution Imaging Spectroradiometer (MODIS) on board Terra/Aqua, Multi-angle Imaging SpectroRadiometer (MISR), Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) are used to derive the dust aerosol optical depth (DOD) in Central Asia based on the Angstrom exponent parameter and/or the particle shape. All sensors agree very well on the interannual variability of DOD. The seasonal analysis of DOD and dust occurrences identified the largest dust loading and the most frequent dust occurrence in the spring and summer, respectively. No significant trend is found during the research period in terms of both DOD and the dust occurrence. Further analysis of Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) aerosol products on a case-by-case basis in most dust months of 2007 suggested that the vertical structure is varying in terms of the extension and the dust loading from one event to another, although dust particles of most episodes have similar physical characteristics (the particle shape and size). Our analysis on the vertical structure of dust plumes, the layer-integrated color ratio and depolarization ratio indicates a varied climate effect (e.g., the direct radiative impact) by mineral dust, dependent on the event being observed in Central Asia.

Traffic source emissions inventories for the rapidly growing West African urban cities are necessary for better local characterization of vehicle emissions released into these cities atmosphere. This study based on local field campaign in a representative site of anthropogenic activities over West African cities such as Yopougon (Abidjan, Côte d'Ivoire) during 2016, provided useful information on vehicle type and age, traveling time, fuel type and amount for fuel consumption estimation. Also, high traffic flow of personal car were recorded on highway, boulevard and backstreet whereas high flows of intra-communal sedan taxi are recorded on main and secondary roads. In addition, the highest daily fuel consumption value of 56 L.day^-1 was recorded in heavy vehicle while the lowest value of 15 L.day^-1 is recorded for personal car using gasoline. This study will be useful for the improvement of uncertainties related to the different databases used to estimate inventories emissions either national or international reports.   This work provides useful information for future studies on urban air quality, climate and health impacts assessment in African cities. It may also be useful for policy makers to support implementation of emission reduction policy in West African cities.

The characteristics of biogenic aerosols in urban area were explored by determining the composition, temporal distribution of saccharides in PM_2.5 in Shanghai. The total saccharides showed a wide range of 15.2 ng/m³ to 1752.8 ng/m³, with the averaged concentrations were 169.8 ng/m³，300.5 ng/m³，288.4 ng/m³，688 ng/m³ in spring, summer, autumn, and winter, respectively. The concerned saccharides include anhydrosaccharides (levoglucosan and mannosan), which were higher in cold seasons due to the increased biomass burning, saccharide alcohols (mannitol, arabitol, sorbitol) and monosaccharides (fructose, glucose), which showed more abundant in warm seasons attributed to the biological emissions. By PMF analysis, four emission sources of saccharides were demonstrated, including biomass burning, fungal spores, soil suspension and plant pollens. Resolution of backward trajectory and fire points showed a process of high concentrations of levoglucosan. We found that concentrations of anhydrosaccharides showed relatively stable under different pollution levels while saccharide alcohols exhibited an obvious decrease, indicated that biomass burning was not the core reason of the heavy haze pollution, however, and high level PM_2.5 pollution might inhibit effects of biological activities.

In the lightning rods categorized as ESE an intermittent voltage impulse is applied to the lightning rod to modulate the electric field at its tip in an attempt to speed up the initiation of a connecting leader from the lightning rod when it is under the influence of a down coming stepped leader. In this paper it is shown that due to the stepping nature of the down coming stepped leader, there is a natural modulation of the electric field at the tip of any lightning rod exposed to the lightning stepped leaders and this modulation is much more intense than any artificial modulation that is possible under practical conditions. Based on the results it is claimed that artificial modulation of the electric field at the tip of lightning rods by applying voltage pulses to it is an unnecessary endeavor because the nature itself has endowed the tip of the lightning rod with a modulating electric field. Thus, as far as the effectiveness of artificial modulation of the tip electric field is concerned, there could be no difference in the lightning attachment efficiency between ESE and Franklin lightning rods.

The Urban Heat Island (UHI) is one of the more serious consequences of urbanization resulting in impacts on thermal comfort levels, heat stress, and even mortality. For Municipal Beirut, implementation of “cool” surface materials and green spaces have been recommended to counterbalance the UHI. This paper builds on previous findings on the topic of non-constructible parcels within the district of Bachoura in Municipal Beirut and examines the possibility of implementing “cool” surface or paving materials and urban vegetation which can improve thermal conditions especially during the summer period and with the viewto project the positive findings of this case study to the entire Municipal Beirut area. A numerical analysis using ENVI-met 4.0 investigates the thermal performance of these non-constructibles further to implementation of high reflective surfaces and urban vegetation within a broad neighborhood scale in Bachoura. Results show reductions in ambient temperatures up to 1K on a summer day.. Within the framework of an integrated approach to planning, this form of urban acupuncture aims for substantial UHI reduction. Energy performance of buildings further to implementation of these mitigation measures is also recommended for future studies and to validate the findings in this paper.

Diablo winds are dry and gusty north-northeasterly downslope windstorms that affect the Northern California. On the evening of October 8, 2017, Diablo winds contributed to the ignitions and rapid spread of the “Northern California Firestorm”, including the Tubbs Fire, which burned 2,800 homes in Santa Rosa, resulted in 22 fatalities and $1.2B USD in damages. We analyzed 18 years of data from a network of surface meteorological stations and show that Diablo winds tend to occur overnight through early morning in Fall, Winter and Spring. Conditions similar to Diablo winds occur commonly in the Sierra Nevada Mountains in addition to the area North of the San Francisco Bay Area including the Northern Coastal Ranges. Diablo winds are characterized by high wind speeds and low relative humidity, but they tend to neither be warmer than climatology nor do they have a higher gust factor, or ratio of wind gusts to mean wind speeds, than climatology.

On 5-7 April 2018 a landfalling atmospheric river resulted in widespread heavy precipitation in the Sierra Nevada of California and Nevada. Observed snow levels during this event were among the highest snow levels recorded since observations began in 2002 and exceeded 2.75 km for 31 hours in the northern Sierra Nevada and 3.75 km for 12 hours in the southern Sierra Nevada. The anomalously high snow levels and over 80 mm of precipitation caused flooding, debris flows, and wet snow avalanches in the upper elevations of the Sierra Nevada. The origin of this atmospheric river was super typhoon Jelawat, whose moisture remnants were entrained and maintained by an extratropical cyclone in the northeast Pacific. This event was notable due to its April occurrence, as six other typhoon remnants that caused heavy precipitation with high snow levels (mean = 2.92 km) in the northern Sierra Nevada all occurred during October.

A record-breaking dust episode took place in Crete on 22 March 2018. The event was characterized by surface concentrations exceeding 1 mg m-3 for a period of 4-7 hours, reaching record values higher than 6 mg m-3 at the background station of Finokalia. We present here a detailed analysis of the atmospheric dynamical processes during this period, to identify the main reasons that resulted in such extreme dust advection over Crete. At synoptic scale, the weakening of the polar vortex and the meridional transport of polar air masses at upper tropospheric layers resulted in a strong jet streak over N. Africa and Central Mediterranean and corresponding tropospheric folding that brought cold stratospheric air in mid and upper troposphere. Cyclogenesis occurred at the Gulf of Sirte in Libya, resulting on strong winds over the NE parts of Libya, enhancing particle emissions. The dust plume traveled at low heights (0.5-3 km) along the warm conveyor belt preceding the depression cold front. This type of dusty southerly wind is commonly known as “Khamsin”. As the flow approached Crete, Foehn winds at the lee side of the island favored the downward mixing of dust towards the surface, resulting in local maxima of PM10 in Heraklion and Finokalia. The analysis is based on the synergy of high resolution WRF-Chem simulations reaching up to 1×1 km grid space over Crete, ground-based and satellite remote sensing of the dust plumes (PollyXT lidar, MSG-SEVIRI, MODIS) and detailed surface aerosol in-situ measurements at urban (Heraklion, Chania) and background (Finokalia) stations in Crete.

In this paper we illustrate a new, simple and complementary ground-based methodology to retrieve the vertically resolved atmospheric precipitation intensity through a synergy between measurements from the National Aeronautics and Space Administration (NASA) Micropulse Lidar network (MPLNET), an analytical model solution and ground-based disdrometer measurements. The presented results are obtained at two mid-latitude MPLNET permanent observational sites, located respectively at NASA Goddard Space Flight Center, USA, and at the Universitat Politècnica de Catalunya, Barcelona, Spain. The methodology is suitable to be applied to existing and/or future lidar/ceilometer networks with the main objective of either providing near-real time (3h latency) rainfall intensity measurements and/or to validate satellite missions, especially for critical light precipitation (<3 mm hr⁻¹).

The air quality monitoring network in Alaska is currently limited to ground-based observations in urban areas and national parks leaving a large proportion of the state unmonitored. The use of MODIS aerosol optical depth (AOD) to estimate ground-level particulate pollution concentrations has been successfully demonstrated around the world, and could potentially be used in Alaska. In this work, MODIS AOD measurements at 550 nm were validated against AOD derived from AERONET ground-based sunphotometers in Barrow and Bonanza Creek to determine if MODIS AOD from the Terra and Aqua satellites could be used to estimate ground-level particulate pollution concentrations. The MODIS AOD was obtained from MODIS collection 6 using the dark target Land and Ocean algorithms from 2000 to 2014. MODIS data could only be obtained between the months of April and October; therefore, it could only be validated for those months. Individual and combined Terra and Aqua MODIS data were considered. The results showed that MODIS collection 6 products at 10 km resolution for Terra and Aqua combined are not valid over land but are valid over the ocean. On the other hand, the individual Terra and Aqua MODIS collection 6 AOD products at 10 km resolution are valid over land individually but not when combined. Results also suggest the MODIS collection 6 AOD products at 3 km resolution are valid over land and ocean and perform better over land than the 10-km product. These findings indicate that MODIS collection 6 AOD products can be used quantitatively in air quality applications in Alaska during the summer months.

Aerosol optical depth (AOT) was measured at Xavier University of Louisiana (29.96ᵒ N, 90.11ᵒ W and about 3m above sea level) using a hand-held sun photometer. AOT was measured at two different wavelengths (green at 505nm and red at 625nm) during the period from Sept-2017 to Feb-2018. In this study, we investigate the relationship between AOT and temperature, precipitation, barometric pressure and relative humidity. The maximum monthly average value for green AOT was 0.176 (measured in September) and the minimum monthly value is 0.040 (measured in January). For the red AOT the maximum monthly average value was 0.123 (measured in September) and the minimum is 0.034 (measured in December). The AOT–temperature relationship was predominantly positive, meaning that high AOT values correspond to high temperatures and low AOT values correspond to low temperatures. Relationship between AOT and rainfall is negative, meaning high rainfall averages are associated with low AOT values. AOT and atmospheric pressure have a predominantly negative relationship. The relationship between relative humidity and AOT is a complicated one and is hard to qualify as relative humidity varies very little during the study period.

The increasing use of renewable energies as a source of electricity has led to a fundamental transition of the power supply system. The integration of fluctuating weather-dependent energy sources into the grid already has a major impact on the load flows of the grid. As a result, the interest in forecasting wind and solar radiation with a sufficient accuracy over short time horizons grew. In this study the short-term forecast of the effective cloud albedo based on optical flow estimation methods are investigated. The optical flow method utilized here is TV-L¹ from the open source library OpenCV. This method uses a multi-scale-approach to capture cloud motions on various spatial scales. After the clouds are displaced the solar surface radiation will be calculated with SPECMAGIC NOW which computes the global irradiation spectrally resolved from satellite imagery. Due to a high temporal and spatial resolution of satellite measurements the effective cloud albedo and thus solar radiation can be forecasted from 5 minutes up to 4 hours with a resolution of 0.05°. In the following there will be a brief description of the method for the short-term forecast of the effective cloud albedo. Subsequently evaluation results will be presented and discussed. Finally an outlook of further developments will be given.

Typhoon Durian forming over the Western North Pacific Ocean and entering into the South China Sea (SCS), caused extreme and widespread damages in 2006. In this research, sensitivity analyses on different physical parameterization schemes of the Weather Research and Forecasting Atmospheric Model (WRF-ATM) have been utilized to study typhoon Durian. Model accuracy and performance testing were investigated with different initial conditions during the tropical cyclone simulation in the SCS. The initial and boundary conditions (IBCs) for all experiments were derived from the European Centre for Medium Range Weather Forecasts (ECMWF), Re-Analysis Interim (ERAI), and the National Centers for Environmental Prediction (NCEP) with Final (FNL) analysis data compiled through the WRF-ATM model. The sensitivity analysis results indicated a major improvement for the cumulus scheme by using the Grell-Devenyi scheme along with the PBL scheme of Yonsei University, mixed-phase microphysics scheme of the WRF Single Moment 5-class and IBCs for ECMWF-ERAI of TC simulation under the context of Wind-Pressure Relationships. This predicted better track and intensity comparing with these of the Joint Typhoon Warning Center. The results revealed that the TC track and intensity were well simulated by the WSM5-GD combination for the WRF-ATM model with an intensity error of 1.69 hPa for minimum surface level pressure, maximum wind speed of 1.83 knots and average track error of 25 km in 72 hours. The simulations showed that the potential track and intensity error decreased with the delayed IBCs, suggesting that the model simulation is more dependable when the coast is approached by the TC.

Site environments and instrumental characteristics of carbon dioxide (CO₂) measurements operated by local governments in the Kanto Plain, the center of which is Tokyo were summarized for this study. The observation sites were classified into environments of three types: urban, suburban, and woodland. Based on a few decades of accumulated hourly data, the diurnal and seasonal variations of CO₂ concentrations were analyzed as a composite of anomalies from annual means recorded for each site. In urban areas, the highest concentrations appear before midnight in winter. The second peak corresponds to the morning rush hour and strengthening of the inversion layer. Suburban areas can be characterized as having the highest concentration before dawn and the lowest concentration during the daytime in summer in association with the activation of respiration and photosynthesis of vegetation. In these areas, concentration peaks also appear during the morning rush hour. Woodland areas show background features, with the highest concentration in early spring: higher than the global background by about 5 ppmv.

AIUS (Atmospheric Infrared Ultraspectral Sounder) is an infrared occultation spectrometer onboard the Chinese GaoFen-5 satellite, which covers a spectral range of 2.4–13.3 μm (750–4100 cm⁻¹) with a spectral resolution of about 0.02 cm⁻¹. AIUS is designed to measure and study chemical processes of ozone (O₃) and other trace gases in the upper troposphere and stratosphere around Antarctic. In this study, the corresponding retrieval methodology is described. The retrieval simulations based on the simulated spectra of AIUS have been carried out, with a focus on O₃. The relative difference between the retrieved and the true O₃ profiles is within 5% from the 15 km to 70 km and about 10% below 15 km. The corresponding averaging kernels illustrate that the overall retrieval information mainly come from the spectra, not the a priori. The retrieval experiments also demonstrate that the shape of the retrieved profiles resembles the shape of the true profile even if the shape of the a priori profile is different from that of the true profile. Further, we perform the O₃ retrieval from the real ACE-FTS (Atmospheric Chemistry Experiment-Fourier Transform Spectrometer) measurements and compare the results with the official ACE-FTS Level-2 products. Overall, both profiles agree well in the stratosphere where the retrieval sensitivity is high. The relative difference between both profiles is about 15% below 70 km, which may due to the measurement errors and different forward model parameters.