The applicability of physical optics in the millimetre and sub-millimetre spectral region. : Part I: The ray tracing with diffraction on facets method

This document is the Accepted Manuscript version of the following article: A. J. Baran, Evelyn Hesse, and Odran Sourdeval, ‘ The applicability of physical optics in the millimetre and sub-millimetre spectral region. Part I: The ray tracing with diffraction on facets method’, Journal of Quantitative Spectroscopy & Radiative Transfer, Vol. 190, (2017), pp. 13-25, first published online on 3 January 2017. The version of record is available online at doi: http://dx.doi.org/10.1016/j.jqsrt.2016.12.030 Crown Copyright © 2016 Published by Elsevier Ltd. All rights reserved.Future satellite missions, from 2022 onwards, will obtain near-global measurements of cirrus at microwave and sub-millimetre frequencies. To realise the potential of these observations, fast and accurate light-scattering methods are required to calculate scattered millimetre and sub-millimetre intensities from complex ice crystals. Here, the applicability of the ray tracing with diffraction on facets method (RTDF) in predicting the bulk scalar optical properties and phase functions of randomly oriented hexagonal ice columns and hexagonal ice aggregates at millimetre frequencies is investigated. The applicability of RTDF is shown to be acceptable down to size parameters of about 18, between the frequencies of 243 and 874 GHz. It is demonstrated that RTDF is generally well within about 10% of T-matrix solutions obtained for the scalar optical properties assuming hexagonal ice columns. Moreover, on replacing electromagnetic scalar optical property solutions obtained for the hexagonal ice aggregate with the RTDF counterparts at size parameter values of about 18 or greater, the bulk scalar optical properties can be calculated to generally well within ±5% of an electromagnetic-based database. The RTDF-derived bulk scalar optical properties result in brightness temperature errors to generally within about ±4 K at 874 GHz. Differing microphysics assumptions can easily exceed such errors. Similar findings are found for the bulk scattering phase functions. This finding is owing to the scattering solutions being dominated by the processes of diffraction and reflection, both being well described by RTDF. The impact of centimetre-sized complex ice crystals on interpreting cirrus polarisation measurements at sub-millimetre frequencies is discussed.Peer reviewe

Constraining the aerosol influence on cloud liquid water path

The impact of aerosols on cloud properties is one of the largest uncertainties in the anthropogenic radiative forcing of the climate. In recent years, significant progress has been made in constraining this forcing using observations, but uncertainty still remains, particularly in the adjustments of cloud properties to aerosol perturbations. Cloud liquid water path (LWP) is the leading control on liquid-cloud albedo, making it important to observationally constrain the aerosol impact LWP. Previous modelling and observational studies have shown that multiple processes play a role in determining the LWP response to aerosol perturbations, but that the aerosol effect can be difficult to isolate. Following previous studies using mediating variables, this work investigates use of the relationship between cloud droplet number concentration (Nd) and LWP for constraining the role of aerosols. Using joint probability histograms to account for the non-linear relationship, this work finds a relationship that is broadly consistent with previous studies. There is significant geographical variation in the relationship, partly due to role of meteorological factors (particularly relative humidity) in the relationship. However, the Nd-LWP relationship is negative in the majority of regions, suggesting that aerosol induced LWP reductions could offset a significant fraction of the radiative forcing from aerosol-cloud interactions (RFaci). However, variations in the Nd-LWP relationship in response to volcanic and shipping aerosol perturbations indicate that the Nd-LWP relationship overestimates the Nd impact on LWP. As such, the estimate of LWP changes due to aerosol in this work provides an upper bound to the radiative forcing from aerosol-induced changes in the LWP

Ice crystal number concentration estimates from lidar-radar satellite retrievals. Part 2: Controls on the ice crystal number concentration

The ice crystal number concentration (Ni) is a keyproperty of ice clouds, both radiatively and microphysically.Due to sparse in situ measurements of ice cloud properties,the controls on theNihave remained difficult to determine.As more advanced treatments of ice clouds are included inglobal models, it is becoming increasingly necessary to de-velop strong observational constraints on the processes in-volved.This work uses the DARDAR-NiceNiretrieval describedin Part 1 to investigate the controls on theNiat a globalscale. The retrieved clouds are separated by type. The ef-fects of temperature, proxies for in-cloud updraft and aerosolconcentrations are investigated. Variations in the cloud topNi(Ni(top)) consistent with both homogeneous and hetero-geneous nucleation are observed along with differing rela-tionships between aerosol andNi(top)depending on the pre-vailing meteorological situation and aerosol type. Away fromthe cloud top, theNidisplays a different sensitivity to thesecontrolling factors, providing a possible explanation for thelowNisensitivity to temperature and ice nucleating particles(INP) observed in previous in situ studies.This satellite dataset provides a new way of investigat-ing the response of cloud properties to meteorological andaerosol controls. The results presented in this work increaseour confidence in the retrievedNiand will form the basis for further study into the processes influencing ice and mixedphase clouds

Cloud base height retrieval from multi-angle satellite data

Clouds are a key modulator of the Earth energy budget at the top of the atmosphere and at the surface. While the cloud top height is operationally retrieved with global coverage, only few methods have been proposed to determine cloud base height (zbase) from satellite measurements. This study presents a new approach to retrieve cloud base heights using the Multi-angle Imaging SpectroRadiometer (MISR) on the Terra satellite. It can be applied if some cloud gaps occur within the chosen distance of typically 10&thinsp;km. The MISR cloud base height (MIBase) algorithm then determines zbase from the ensemble of all MISR cloud top heights retrieved at a 1.1&thinsp;km horizontal resolution in this area. MIBase is first calibrated using 1 year of ceilometer data from more than 1500 sites within the continental United States of America. The 15th percentile of the cloud top height distribution within a circular area of 10&thinsp;km radius provides the best agreement with the ground-based data. The thorough evaluation of the MIBase product zbase with further ceilometer data yields a correlation coefficient of about 0.66, demonstrating the feasibility of this approach to retrieve zbase. The impacts of the cloud scene structure and macrophysical cloud properties are discussed. For a 3-year period, the median zbase is generated globally on a 0.25∘&thinsp;×&thinsp;0.25∘ grid. Even though overcast cloud scenes and high clouds are excluded from the statistics, the median zbase retrievals yield plausible results, in particular over ocean as well as for seasonal differences. The potential of the full 16 years of MISR data is demonstrated for the southeast Pacific, revealing interannual variability in zbase in accordance with reanalysis data. The global cloud base data for the 3-year period (2007–2009) are available at https://doi.org/10.5880/CRC1211DB.19.</p

Analysis and quantification of ENSO-linked changes in the tropical Atlantic cloud vertical distribution using 14 years of MODIS observations

A total of 14 years (September 2002 to September 2016) of Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) monthly mean cloud data are used to quantify possible changes in the cloud vertical distribution over the tropical Atlantic. For the analysis multiple linear regression techniques are used. For the investigated time period significant linear changes were found in the domain-averaged cloud-top height (CTH) (−178 m per decade), the high-cloud fraction (HCF) (−0.0006 per decade), and the low-cloud amount (0.001 per decade). The interannual variability of the time series (especially CTH and HCF) is highly influenced by the El Niño–Southern Oscillation (ENSO). Separating the time series into two phases, we quantified the linear change associated with the transition from more La Niña-like conditions to a phase with El Niño conditions (Phase 2) and vice versa (Phase 1). The transition from negative to positive ENSO conditions was related to a decrease in total cloud fraction (TCF) (−0.018 per decade; not significant) due to a reduction in the high-cloud amount (−0.024 per decade; significant). Observed anomalies in the mean CTH were found to be mainly caused by changes in HCF rather than by anomalies in the height of cloud tops themselves. Using the large-scale vertical motion ω at 500 hPa (from ERA-Interim ECMWF reanalysis data), the observed anomalies were linked to ENSO-induced changes in the atmospheric large-scale dynamics. The most significant and largest changes were found in regions with strong large-scale upward movements near the Equator. Despite the fact that with passive imagers such as MODIS it is not possible to vertically resolve clouds, this study shows the great potential for large-scale analysis of possible changes in the cloud vertical distribution due to the changing climate by using vertically resolved cloud cover and linking those changes to large-scale dynamics using other observations or model data

Ice crystal number concentration estimates from lidar–radar satellite remote sensing – Part 2: Controls on the ice crystal number concentration

The ice crystal number concentration (Ni) is a key property of ice clouds, both radiatively and microphysically. Due to sparse in situ measurements of ice cloud properties, the controls on the Ni have remained difficult to determine. As more advanced treatments of ice clouds are included in global models, it is becoming increasingly necessary to develop strong observational constraints on the processes involved. This work uses the DARDAR-Nice Ni retrieval described in Part 1 to investigate the controls on the Ni at a global scale. The retrieved clouds are separated by type. The effects of temperature, proxies for in-cloud updraft and aerosol concentrations are investigated. Variations in the cloud top Ni (Ni(top)) consistent with both homogeneous and heterogeneous nucleation are observed along with differing relationships between aerosol and Ni(top) depending on the prevailing meteorological situation and aerosol type. Away from the cloud top, the Ni displays a different sensitivity to these controlling factors, providing a possible explanation for the low Ni sensitivity to temperature and ice nucleating particles (INP) observed in previous in situ studies. This satellite dataset provides a new way of investigating the response of cloud properties to meteorological and aerosol controls. The results presented in this work increase our confidence in the retrieved Ni and will form the basis for further study into the processes influencing ice and mixed phase clouds

Ice crystal number concentration estimates from lidar–radar satellite remote sensing – Part 2: Controls on the ice crystal number concentration

The ice crystal number concentration (Ni) is a key property of ice clouds, both radiatively and microphysically. Due to sparse in situ measurements of ice cloud properties, the controls on the Ni have remained difficult to determine. As more advanced treatments of ice clouds are included in global models, it is becoming increasingly necessary to develop strong observational constraints on the processes involved.This work uses the DARDAR-Nice Ni retrieval described in Part 1 to investigate the controls on the Ni at a global scale. The retrieved clouds are separated by type. The effects of temperature, proxies for in-cloud updraft and aerosol concentrations are investigated. Variations in the cloud top Ni (Ni(top)) consistent with both homogeneous and heterogeneous nucleation are observed along with differing relationships between aerosol and Ni(top) depending on the prevailing meteorological situation and aerosol type. Away from the cloud top, the Ni displays a different sensitivity to these controlling factors, providing a possible explanation for the low Ni sensitivity to temperature and ice nucleating particles (INP) observed in previous in situ studies.This satellite dataset provides a new way of investigating the response of cloud properties to meteorological and aerosol controls. The results presented in this work increase our confidence in the retrieved Ni and will form the basis for further study into the processes influencing ice and mixed phase clouds.</p

A variational approach for retrieving ice cloud properties from infrared measurements: application in the context of two IIR validation campaigns

Cirrus are cloud types that are recognized to have a strong impact on the Earth-atmosphere radiation balance. This impact is however still poorly understood, due to the difficulties in describing the large variability of their properties in global climate models. Consequently, numerous airborne and space borne missions have been dedicated to their study in the last decades. The satellite constellation A-Train has proven to be particularly helpful to study cirrus on global scale due to such instruments as the Infrared Imaging Radiometer (IIR), which shows great sensitivity to the radiative and microphysical properties of these clouds. This study presents an algorithm that uses thermal infrared measurements to retrieve the optical thickness of cirrus and the effective size of their ice crystals. This algorithm is based on an optimal estimation scheme, which possesses the advantage of attributing precise uncertainties to the retrieved parameters. Two IIR airborne validation campaigns have been chosen as case studies. It is observed that optical thicknesses could be accurately retrieved but that large uncertainties may occur on the effective diameters. Strong agreements have been found between the products of our algorithm when separately applied to the measurements of IIR and of the airborne radiometer CLIMAT-AV, which comforts the results of previous validations of IIR level-1 measurements. Comparisons with in situ observations and with operational products of IIR also show confidence in our results. However, we have found that the quality of our retrievals can be strongly impacted by uncertainties related to the choice of a pristine crystal model and by poor constraints on the properties of possible liquid cloud layers underneath cirrus. Simultaneous retrievals of liquid clouds radiative and microphysical properties or the use of different ice crystal models should therefore be considered to improve the quality of the results

Characterisation of low-base and mid-base clouds and their thermodynamic phase over the Southern Ocean and Arctic marine regions

The thermodynamic phase of clouds in low and middle levels over the Southern Ocean and the Arctic marine regions is poorly known, leading to uncertainties in the radiation budget in weather and climate models. To improve the knowledge of the cloud phase, we analyse 2 years of the raDAR-liDAR (DARDAR) dataset based on active satellite instruments. We classify clouds according to their base and top height and focus on low-, mid-, and mid- to low-level clouds as they are the most frequent in the mixed-phase temperature regime. Low-level single-layer clouds occur in 8 %–15 % of all profiles, but single-layer clouds spanning the mid-level also amount to approx. 15 %. Liquid clouds show mainly a smaller vertical extent but a horizontally larger extent compared to ice clouds. The results show the highest liquid fractions for low-level and mid-level clouds. Two local minima in the liquid fraction are observed around cloud top temperatures of −15 and −5 °C. Mid-level and mid- to low-level clouds over the Southern Ocean and low-level clouds in both polar regions show higher liquid fractions if they occur over sea ice compared to the open ocean. Low-level clouds and mid- to low-level clouds with high sea salt concentrations, used as a proxy for sea spray, show reduced liquid fractions. In mid-level clouds, dust shows the largest correlations with liquid fraction, with a lower liquid fraction for a higher dust aerosol concentration. Low-level clouds clearly show the largest contribution to the shortwave cloud radiative effect in both polar regions, followed by mid- to low-level clouds.</p

Cold cloud microphysical process rates in a global chemistry–climate model

Microphysical processes in cold clouds which act as sources or sinks of hydrometeors below 0 degrees C control the ice crystal number concentrations (ICNCs) and in turn the cloud radiative effects. Estimating the relative importance of the cold cloud microphysical process rates is of fundamental importance to underpin the development of cloud parameterizations for weather, atmospheric chemistry, and climate models and to compare the output with observations at different temporal resolutions. This study quantifies and investigates the ICNC rates of cold cloud microphysical processes by means of the chemistry-climate model EMAC (ECHAM/MESSy Atmospheric Chemistry) and defines the hierarchy of sources and sinks of ice crystals. Both microphysical process rates, such as ice nucleation, aggregation, and secondary ice production, and unphysical correction terms are presented. Model ICNCs are also compared against a satellite climatology. We found that model ICNCs are in overall agreement with satellite observations in terms of spatial distribution, although the values are overestimated, especially around high mountains. The analysis of ice crystal rates is carried out both at global and at regional scales. We found that globally the freezing of cloud droplets and convective detrainment over tropical land masses are the dominant sources of ice crystals, while aggregation and accretion act as the largest sinks. In general, all processes are characterized by highly skewed distributions. Moreover, the influence of (a) different ice nucleation parameterizations and (b) a future global warming scenario on the rates has been analysed in two sensitivity studies. In the first, we found that the application of different parameterizations for ice nucleation changes the hierarchy of ice crystal sources only slightly. In the second, all microphysical processes follow an upward shift in altitude and an increase by up to 10 % in the upper troposphere towards the end of the 21st century