Aviation impact on the cirrus ice crystal number using satellite observation

<p align="justify">Aviation outflow is the only anthropogenic source of pollution that is directly emitted into the upper troposphere. This emission has the potential to modify the cloudiness directly by forming linear contrails and indirectly by injecting aerosols, which can act as cloud condensation nuclei (CCN) and ice nucleating particles (INP). Contrail cirrus can persist either in cloud-free supersaturated air, increasing high-cloud cover or inside natural cirrus cloud, and therefore modifying the microphysical properties of already existing cirrus clouds. Even though the situation that an aircraft flies through a natural cirrus is one of the highly probable situations in the upper troposphere, its subsequent impact is unclear with the present state of knowledge. Quantifying such impact is necessary if we are to properly account for the influence of aviation on climate. One main limitation preventing us to better identify these impacts is the lack of height resolved measurements inside the cirrus clouds.</p> <p align="justify">In this study, we used new retrievals from combined satellite cloud radar and lidar (Cloud- Sat/CALIPSO; DARDAR-Nice algorithm), which provide height resolved information of ice crystal number concentration, at intercepts between the CALIPSO ground track and the position of civil aircraft operating between the west coast of the continental United States (Seattle, San Francisco and Los Angeles) and Hawaii during 2010 and 2011 from an earlier study.</p> <p align="justify">Comparing cloudy air behind the aircraft inside the flight track to the adjacent regions and to ahead of the aircraft revealed a notable difference in ice number concentration at 300 m to 540 m beneath the flight height. These differences are derived from the reduction of ice number concentrations as we proceed toward the cloud base in regions unaffected by aviation and the increase of ice crystals as we distance a few hundreds of meters beneath the flight level in the regions affected by aviation.</p>

Exploring the uncertainties in the aviation soot–cirrus effect

A global aerosol–climate model, including a two-moment cloud microphysical scheme and a parametrization for aerosol-induced ice formation in cirrus clouds, is applied in order to quantify the impact of aviation soot on natural cirrus clouds. Several sensitivity experiments are performed to assess the uncertainties in this effect related to (i) the assumptions on the ice nucleation abilities of aviation soot, (ii) the representation of vertical updrafts in the model, and (iii) the use of reanalysis data to relax the model dynamics (the so-called nudging technique). Based on the results of the model simulations, a radiative forcing from the aviation soot–cirrus effect in the range of −35 to 13 mW m−2 is quantified, depending on the assumed critical saturation ratio for ice nucleation and active fraction of aviation soot but with a confidence level below 95 % in several cases. Simple idealized experiments with prescribed vertical velocities further show that the uncertainties on this aspect of the model dynamics are critical for the investigated effect and could potentially add a factor of about 2 of further uncertainty to the model estimates of the resulting radiative forcing. The use of the nudging technique to relax model dynamics is proved essential in order to identify a statistically significant signal from the model internal variability, while simulations performed in free-running mode and with prescribed sea-surface temperatures and sea-ice concentrations are shown to be unable to provide robust estimates of the investigated effect. A comparison with analogous model studies on the aviation soot–cirrus effect show a very large model diversity, with a conspicuous lack of consensus across the various estimates, which points to the need for more in-depth analyses on the roots of such discrepancies.

Enhanced soot particle ice nucleation ability induced by aggregate compaction and densification

Soot particles, acting as ice nucleating particles (INPs), can contribute to cirrus cloud formation which has an important influence on climate. Aviation activities emitting soot particles in the upper troposphere can potentially impact ice nucleation (IN) in cirrus clouds. Pore condensation and freezing (PCF) is an important ice formation pathway for soot particles in the cirrus regime, which requires the soot INP to have specific morphological properties, i.e. mesopore structures. In this study, the morphology and pore size distribution of two kinds of soot samples were modified by a physical agitation method without any chemical modification, by which more compacted soot sample aggregates could be produced compared to the unmodified sample. The IN activities of both fresh and compacted soot particles with different sizes, 60, 100, 200 and 400 nm, were systematically tested by the Horizontal Ice Nucleation Chamber (HINC) under mixed-phase and cirrus clouds relevant temperatures (T). Our results show that soot particles are unable to form ice crystals at T > 235 K (homogeneous nucleation temperature, HNT) but IN was observed for compacted and larger size soot aggregates (> 200 nm) well below homogeneous freezing relative humidity (RHhom) at T < HNT, demonstrating PCF as the dominating mechanism for soot IN. We also observed that mechanically compacted soot particles can reach a higher particle activation fraction (AF) value for the same T and RH condition, compared to the same aggregate size fresh soot particles. The results also reveal a clear size dependence for the IN activity of soot particles with the same agitation degree, showing that compacted soot particles with large sizes (200 and 400 nm) are more active INPs and can convey the single importance of soot aggregate morphology for the IN ability. In order to understand the role of soot aggregate morphology for its IN activity, both fresh and compacted soot samples were characterized systematically using particle mass and size measurements, comparisons from TEM (transmission electron microscopy) images, soot porosity characteristics from argon (Ar) and nitrogen (N2) physisorption measurements, as well as soot-water interaction results from DVS (dynamic vapor sorption) measurements. Considering the soot particle physical properties along with its IN activities, the enhanced IN abilities of compacted soot particles are attributed to decreasing mesopore width and increasing mesopore occurrence probability due to the compaction process.

Empirical model of multiple scattering effect on single-wavelength lidar data of aerosols and clouds

We performed extensive Monte Carlo (MC) simulations of single-wavelength lidar signals from a plane-parallel homogeneous layer of atmospheric particles and developed an empirical model to account for the multiple scattering in the lidar signals. The simulations have taken into consideration four types of lidar configurations (the ground based, the airborne, the CALIOP, and the ATLID) and four types of particles (coarse aerosol, water cloud, jet-stream cirrus and cirrus). Most of simulations were performed with the spatial resolution of 20 m and the particles extinction coefficient εp between 0.06 km−1 and 1.0 km−1. The resolution was of 5 m for high values of εp (up to 10.0 km−1). The majority of simulations for ground-based and airborne lidars were performed at two values of the receiver field-of-view (RFOV): 0.25 mrad and 1.0 mrad. The effect of the width of the RFOV was studied for the values up to 50 mrad. The proposed empirical model is a function that has only three free parameters and approximates the multiple-scattering relative contribution to lidar signals. It is demonstrated that the empirical model has very good quality of MC data fitting for all considered cases. Special attention was given to the usual operational conditions, i.e., low distances to a particles layer, small optical depths and quite narrow receiver field-of-views. It is demonstrated that multiple scattering effects cannot be neglected when the distance to a particles layer is about 8 km or higher and the full RFOV is of 1.0 mrad. As for the full RFOV of 0.25 mrad, the single scattering approximation is acceptable for aerosols (εp ≲ 1.0 km−1), water clouds (εp ≲ 0.5 km−1), and cirrus clouds (εp ≤ 0.1 km−1). When the distance to a particles layer is of 1 km, the single scattering approximation is acceptable for aerosols and water clouds (εp ≲ 1.0 km−1, both RFOV = 0.25 and RFOV = 1 mrad). As for cirrus clouds, the effect of multiple scattering cannot be neglected even at such low distance when εp ≳ 0.5 km−1.

Properties of Cirrus Clouds over the European Arctic (Ny-Ålesund, Svalbard)

Cirrus is the only cloud type capable of inducing daytime cooling or heating at the top of the atmosphere (TOA) and the sign of its radiative effect highly depends on its optical depth. However, the investigation of its geometrical and optical properties over the Arctic is limited. In this work the long-term properties of cirrus clouds are explored for the first time over an Arctic site (Ny-Ålesund, Svalbard) using lidar and radiosonde measurements from 2011 to 2020. The optical properties were quality assured, taking into account the effects of specular reflections and multiple-scattering. Cirrus clouds were generally associated with colder and calmer wind conditions compared to the 2011–2020 climatology. However, the dependence of cirrus properties on temperature and wind speed was not strong. Even though the seasonal cycle was not pronounced, the winter-time cirrus appeared under lower temperatures and stronger wind conditions. Moreover, in winter, geometrically- and optically-thicker cirrus were found and their ice particles tended to be more spherical. The majority of cirrus was associated with westerly flow and westerly cirrus tended to be geometrically-thicker. Overall, optically-thinner layers tended to comprise smaller and less spherical ice crystals, most likely due to reduced water vapor deposition on the particle surface. Compared to lower latitudes, the cirrus layers over Ny-Ålesund were more absorbing in the visible spectral region and they consisted of more spherical ice particles.

Monte Carlo simulation of halos in the crystal clouds

In this paper, we try to answer the question: how the multiple scattering, the sun elevation, shape and orientation of ice crystals in the cirrus clouds affect a halo pattern. To study the radiation transfer in optically anisotropic clouds, we have developed the software based on Monte Carlo method and ray tracing. In addition to halos, this software enables one to simulate “anti-halos”, which above the cloud layer can be seen by observers. We present the visualization of halos and anti-halos generated by the cirrus clouds for different shapes and orientations of ice crystals.

Vertical Structure of Ice Clouds and Vertical Air Motion from Vertically Pointing Cloud Radar Measurements

The vertical structure of ice clouds and vertical air motion (Vair) were investigated using vertically pointing Ka-band cloud radar. The distributions of reflectivity (Z), Doppler velocity (VD), and spectrum width (SW) were analyzed for three ice cloud types, namely, cirrus, anvil, and stratiform clouds. The radar parameters of the cirrus clouds showed narrower distributions than those of the stratiform and anvil clouds. In the vertical structures, the rapid growth of Z and VD occurred in the layer between 8 and 12 km (roughly a layer of −40 °C to −20 °C) for all ice clouds. The prominent feature in the stratiform clouds was an elongated “S” shape in the VD near 7–7.5 km (at approximately −16 °C to −13 °C) due to a significant decrease in an absolute value of VD. The mean terminal fall velocity (Vt) and Vair in the ice clouds were estimated using pre-determined Vt–Z relationships (Vt = aZb) and the observed VD. Although the cirrus clouds demonstrated wide distributions in coefficients a and exponents b depending on cloud heights, they showed a smaller change in Z and Vt values compared to that of the other cloud types. The anvil clouds had a larger exponent than that of the stratiform clouds, indicating that the ice particle density of anvil clouds increases at a faster rate compared with the density of stratiform clouds for the same Z increment. The significant positive Vair appeared at the top of all ice clouds in range up to 0.5 m s−1, and the anvil clouds showed the deepest layer of upward motion. The stratiform and anvil clouds showed a dramatic increase in vertical air motion in the layer of 6–8 km as shown by the rapid decrease of VD. This likely caused increase of supersaturation above. A periodic positive Vair linked with a significant reduction in VD appeared at the height of 7–8 km (approximately −15 °C) dominantly in the stratiform clouds. This layer exhibited a bi-modal power spectrum produced by pre-existing larger ice particles and newly formed numerous smaller ice particles. This result raised a question on the origins of smaller ice particles such as new nucleation due to increased supersaturation by upward motion below or the seeder-feeder effect. In addition, the retrieved Vair with high-resolution data well represented a Kelvin-Helmholtz wave development.