A novel bulk microphysics scheme that predicts the evolution of ice properties, including aspect ratio (shape), mass, number, size, and density is described, tested, and demonstrated. The scheme is named the Ice-Spheroids Habit Model with Aspect-Ratio Evolution (ISHMAEL). Ice is modeled as spheroids and is nucleated as one of two species depending on nucleation temperature. Microphysical process rates determine how shape and other ice properties evolve. A third aggregate species is also employed, diversifying ice properties in the model. Tests of ice shape evolution during vapor growth and riming are verified against wind tunnel data, revealing that the model captures habit-dependent riming and its effect on fall speed. Lagrangian parcel studies demonstrate that the bulk model captures ice property evolution during riming and melting compared with a bin model. Finally, the capabilities of ISHMAEL are shown in a 2D kinematic framework with a simple updraft. A direct result of predicting ice shape evolution is that various states of ice from unrimed to lightly rimed to densely rimed can be modeled without converting ice mass between predefined ice categories (e.g., snow and graupel). This leads to a different spatial precipitation distribution compared with the traditional method of separating snow and graupel and converting between the two categories, because ice in ISHMAEL sorts in physical space based on the amount of rime, which controls the thickness and therefore fall speed. Predicting these various states of rimed ice leads to a reduction in vapor growth rate and an increase in riming rate in a simple updraft compared with the traditional approach.

This paper describes and tests a single-particle ice growth model that evolves both ice crystal mass and shape as a result of vapor growth and riming. Columnar collision efficiencies in the model are calculated using a new theoretical method derived from spherical collision efficiencies. The model is able to evolve mass, shape, and fall speed of growing ice across a range of temperatures, and it compares well with wind tunnel data. The onset time of riming and the effects of riming on mass and fall speed between −3° and −16°C are modeled, as compared with wind tunnel data for a liquid water content of 0.4 g m−3. Under these conditions, riming is constrained to the more isometric habits near −10° and −4°C. It is shown that the mass and fall speed of riming dendrites depend on the liquid drop distribution properties, leading to a range of mass–size and fall speed–size relationships. Riming at low liquid water contents is shown to be sensitive to ice crystal habit and liquid drop size. Moreover, very light riming can affect the shape of ice crystals enough to reduce vapor growth and suppress overall mass growth, as compared with those same ice crystals if they were unrimed.

[1] Ozone profiles from balloon-borne ozonesondes are used for development of satellite algorithms and in chemistry-climate model initialization, assimilation and evaluation. An important issue in the application of these profiles is how best to treat variations where varying photochemical and dynamical influences can cause the ozone mixing ratio in the tropospheric segments of the profile to change by of a factor of 2-3 within a day. Clustering techniques are an ideal way to approach the statistical classification of profile data and we apply self-organizing maps to tropical tropospheric SHADOZ data, hypothesizing that the data will sort according to various influences on ozone, namely anthropogenic sources like biomass burning, meteorological conditions, and stratospheric or extra-tropical intrusions. Self-organizing maps, that use a learning algorithm to reveal the most prominent features of a data set according to a specified number of clusters, have been determined for the 1998-2009 SHADOZ profiles over Ascension Island (512 profiles, 7.98°S, 14.42°W) and Natal, Brazil (425 profiles, 5.42°S, 35.38°W). The 2 Â 2 self-organizing map, which creates 4 clusters, reveals that deviations from the average ozone in the free troposphere include both increased ozone resulting from seasonal biomass burning in Africa and locally reduced ozone brought about by convective lifting of unpolluted boundary-layer air. Expanding to a 4 Â 4 self-organizing map shows how biomass burning influences the yearly cycle of tropospheric ozone at Ascension Island and captures the seasonality of ozone at both Ascension Island and Natal. Comparing Ascension Island and Natal using a 4 Â 4 self-organizing map at each site reveals similarities in mid-tropospheric ozone, but shows differences in lower-tropospheric ozone due to Ascension Island being closer to African biomass burning and more affected by descent from the mean Walker circulation, with less convective activity, than Natal.

Concentrations of airborne chemical and biological agents from a hazardous release are not spread uniformly. Instead, there are regions of higher concentration, in part due to local atmospheric flow conditions which can attract agents. We equipped a ground station and two rotary-wing unmanned aircraft systems (UASs) with ultrasonic anemometers. Flights reported here were conducted 10 to 15 m above ground level (AGL) at the Leach Airfield in the San Luis Valley, Colorado as part of the Lower Atmospheric Process Studies at Elevation—a Remotely-Piloted Aircraft Team Experiment (LAPSE-RATE) campaign in 2018. The ultrasonic anemometers were used to collect simultaneous measurements of wind speed, wind direction, and temperature in a fixed triangle pattern; each sensor was located at one apex of a triangle with ∼100 to 200 m on each side, depending on the experiment. A WRF-LES model was used to determine the wind field across the sampling domain. Data from the ground-based sensors and the two UASs were used to detect attracting regions (also known as Lagrangian Coherent Structures, or LCSs), which have the potential to transport high concentrations of agents. This unique framework for detection of high concentration regions is based on estimates of the horizontal wind gradient tensor. To our knowledge, our work represents the first direct measurement of an LCS indicator in the atmosphere using a team of sensors. Our ultimate goal is to use environmental data from swarms of sensors to drive transport models of hazardous agents that can lead to real-time proper decisions regarding rapid emergency responses. The integration of real-time data from unmanned assets, advanced mathematical techniques for transport analysis, and predictive models can help assist in emergency response decisions in the future.

Because unmanned aircraft systems (UAS) offer new perspectives on the atmosphere, their use in atmospheric science is expanding rapidly. In support of this growth, the International Society for Atmospheric Research Using Remotely-Piloted Aircraft (ISARRA) has been developed and has convened annual meetings and “flight weeks.” The 2018 flight week, dubbed the Lower Atmospheric Profiling Studies at Elevation–A Remotely-Piloted Aircraft Team Experiment (LAPSE-RATE), involved a 1-week deployment to Colorado’s San Luis Valley. Between 14 and 20 July 2018 over 100 students, scientists, engineers, pilots, and outreach coordinators conducted an intensive field operation using unmanned aircraft and ground-based assets to develop datasets, community, and capabilities. In addition to a coordinated “Community Day” which offered a chance for groups to share their aircraft and science with the San Luis Valley community, LAPSE-RATE participants conducted nearly 1,300 research flights totaling over 250 flight hours. The measurements collected have been used to advance capabilities (instrumentation, platforms, sampling techniques, and modeling tools), conduct a detailed system intercomparison study, develop new collaborations, and foster community support for the use of UAS in atmospheric science.

A quasi-idealized 3D squall-line case is simulated using a novel bulk microphysics scheme called the Ice-Spheroids Habit Model with Aspect-ratio Evolution (ISHMAEL). In ISHMAEL, the evolution of ice particle properties (e.g., mass, shape, maximum diameter, density, and fall speed) are predicted during vapor growth, sublimation, riming, and melting, allowing ice properties to evolve from various microphysical processes without needing separate unrimed and rimed ice categories. ISHMAEL produces both a transition zone and an enhanced stratiform precipitation region, and ice particle properties are analyzed to determine the characteristics of ice that lead to the development of these squall-line features. Rimed particles advected rearward from the convective region produce the enhanced stratiform precipitation region. The transition zone results from hydrometeor sorting; the evolution of ice particle properties in the convective region leads to fall speeds that favor ice advecting rearward of the transition zone before reaching the melting level, causing a local minimum in precipitation rate and reflectivity there. Sensitivity studies show that the fall speed of ice particles largely determines the location of the enhanced stratiform precipitation region and whether or not a transition zone forms. The representation of microphysical processes, such as rime splintering and aggregation, and ice size distribution shape can impact the mean ice particle fall speeds enough to significantly impact the location of the enhanced stratiform precipitation region and the existence of the transition zone.

The self-organizing map (SOM) statistical technique is applied to vertical profiles of thermodynamic and kinematic parameters from a Rapid Update Cycle-2 (RUC-2) proximity sounding dataset with the goal of better distinguishing and predicting supercell and tornadic environments. An SOM is a topologically ordered mapping of input data onto a two-dimensional array of nodes that can be used to classify large datasets into meaningful clusters. The relative ability of SOMs derived from each parameter to separate soundings in a way that is useful in discriminating between storm type, location, and time of year is discussed. Sensitivity to SOM configuration is also explored. Simple skill scores are computed for each SOM to evaluate the relative potential of each variable for future development as a method of probabilistic forecasting. It is found that variance in SOM nodes is reduced compared to the overall dataset, indicating that this is a viable classification method. SOMs of profiles of wind-derived variables are more effective in discriminating between storm type than thermodynamic variables. The SOM method also identifies meteorological, geographic, and temporal regimes within the dataset. In general, conditional probabilities of storm-type occurrence generated using SOMs have higher skill when wind-derived variables are considered and when forecasting nonsupercell events. Stormrelative wind variables tend to have better skill than ground-relative wind variables when forecasting nonsupercells, whereas ground-relative variables become more important when forecasting tornadoes.

Abstract. Unmanned aircraft systems (UASs) offer innovative capabilities for providing new perspectives on the atmosphere, and therefore atmospheric scientists are rapidly expanding their use, particularly for studying the planetary boundary layer. In support of this expansion, from 14 to 20 July 2018 the International Society for Atmospheric Research using Remotely piloted Aircraft (ISARRA) hosted a community flight week, dubbed the Lower Atmospheric Profiling Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE; de Boer et al., 2020a). This field campaign spanned a 1-week deployment to Colorado's San Luis Valley, involving over 100 students, scientists, engineers, pilots, and outreach coordinators. These groups conducted intensive field operations using unmanned aircraft and ground-based assets to develop comprehensive datasets spanning a variety of scientific objectives, including a total of nearly 1300 research flights totaling over 250 flight hours. This article introduces this campaign and lays the groundwork for a special issue on the LAPSE-RATE project. The remainder of the special issue provides detailed overviews of the datasets collected and the platforms used to collect them. All of the datasets covered by this special issue have been uploaded to a LAPSE-RATE community set up at the Zenodo data archive (https://zenodo.org/communities/lapse-rate/, last access: 3 December 2020).