0000000001192713
AUTHOR
Peter Spichtinger
Thermodynamic correction of particle concentrations measured by underwing probes on fast-flying aircraft
Particle concentration measurements with underwing probes on aircraft are impacted by air compression upstream of the instrument body as a function of flight velocity. In particular, for fast-flying aircraft the necessity arises to account for compression of the air sample volume. Hence, a correction procedure is needed to invert measured particle number concentrations to ambient conditions that is commonly applicable to different instruments to gain comparable results. In the compression region where the detection of particles occurs (i.e. under factual measurement conditions), pressure and temperature of the air sample are increased compared to ambient (undisturbed) conditions in certain …
Stochastic Galerkin method for cloud simulation
AbstractWe develop a stochastic Galerkin method for a coupled Navier-Stokes-cloud system that models dynamics of warm clouds. Our goal is to explicitly describe the evolution of uncertainties that arise due to unknown input data, such as model parameters and initial or boundary conditions. The developed stochastic Galerkin method combines the space-time approximation obtained by a suitable finite volume method with a spectral-type approximation based on the generalized polynomial chaos expansion in the stochastic space. The resulting numerical scheme yields a second-order accurate approximation in both space and time and exponential convergence in the stochastic space. Our numerical results…
IMEX Finite Volume Methods for Cloud Simulation
We present new implicit-explicit (IMEX) finite volume schemes for numerical simulation of cloud dynamics. We use weakly compressible equations to describe fluid dynamics and a system of advection-diffusion-reaction equations to model cloud dynamics. In order to efficiently resolve slow dynamics we split the whole nonlinear system in a stiff linear part governing the acoustic and gravitational waves as well as diffusive effects and a non-stiff nonlinear part that models nonlinear advection effects. We use a stiffly accurate second order IMEX scheme for time discretization to approximate the stiff linear operator implicitly and the non-stiff nonlinear operator explicitly. Fast microscale clou…
Tropical tropopause ice clouds: a dynamic approach to the mystery of low crystal numbers
Abstract. The occurrence of high, persistent ice supersaturation inside and outside cold cirrus in the tropical tropopause layer (TTL) remains an enigma that is intensely debated as the "ice supersaturation puzzle". However, it was recently confirmed that observed supersaturations are consistent with very low ice crystal concentrations, which is incompatible with the idea that homogeneous freezing is the major method of ice formation in the TTL. Thus, the tropical tropopause "ice supersaturation puzzle" has become an "ice nucleation puzzle". To explain the low ice crystal concentrations, a number of mainly heterogeneous freezing methods have been proposed. Here, we reproduce in situ measure…
Reappraising the appropriate calculation of a common meteorological quantity: Potential Temperature
Abstract. The potential temperature is a widely used quantity in atmospheric science since it is conserved for air's adiabatic changes of state. Its definition involves the specific heat capacity of dry air, which is traditionally assumed as constant. However, the literature provides different values of this allegedly constant parameter, which are reviewed and discussed in this study. Furthermore, we derive the potential temperature for a temperature-dependent parameterization of the specific heat capacity of dry air, thus providing a new reference potential temperature with a more rigorous basis. This new reference shows different values and vertical gradients in the upper troposphere and …
In situ detection of stratosphere-troposphere exchange of cirrus particles in the midlatitudes
Airborne trace gas, microphysical, and radiation measurements were performed during the AIRcraft TOwed Sensor Shuttle - Inhomogeneous Cirrus Experiment over northern Germany in 2013. Based on high-precision nitrous oxide (N2O) and carbon monoxide (CO) in situ data, stratospheric air could be identified, which contained cirrus cloud particles. Consistent with the stratospheric N2O data, backward trajectories indicate that the sampled air masses crossed the dynamical tropopause in the last 3 h before the measurement. These air masses contained cirrus particles, which were formed during slow ascent in the troposphere and subsequently mixed with stratospheric air. From the CO-N2O correlation th…
Pattern formation in clouds via Turing instabilities
Pattern formation in clouds is a well-known feature, which can be observed almost every day. However, the guiding processes for structure formation are mostly unknown, and also theoretical investigations of cloud patterns are quite rare. From many scientific disciplines the occurrence of patterns in non-equilibrium systems due to Turing instabilities is known, i.e. unstable modes grow and form spatial structures. In this study we investigate a generic cloud model for the possibility of Turing instabilities. For this purpose, the model is extended by diffusion terms. We can show that for some cloud models, i.e special cases of the generic model, no Turing instabilities are possible. However,…
New particle formation inside ice clouds: In-situ observations in the tropical tropopause layer of the 2017 Asian Monsoon Anticyclone
Abstract. From 27 July to 10 August 2017 the airborne StratoClim mission took place in Kathmandu, Nepal where eight mission flights were conducted with the M-55 Geophysica up to altitudes of 20 km. New Particle Formation (NPF) was identified by the abundant presence of ultrafine aerosols, with particle diameters dp smaller than 15 nm, which were in-situ detected by means of condensation nuclei counting techniques. NPF fields in clear-skies as well as in the presence of cloud ice particles (dp > 3 µm) were encountered at upper troposphere/lowermost stratosphere (UT/LS) levels and within the Asian Monsoon Anticyclone (AMA). NPF-generated ultrafine particles in elevated concentrations (Nuf)…
Automated detection and classification of synoptic-scale fronts from atmospheric data grids
Automatic determination of fronts from atmospheric data is an important task for weather prediction as well as for research of synoptic-scale phenomena. In this paper we introduce a deep neural network to detect and classify fronts from multi-level ERA5 reanalysis data. Model training and prediction is evaluated using two different regions covering Europe and North America with data from two weather services. We apply label deformation within our loss function, which removes the need for skeleton operations or other complicated post-processing steps as used in other work, to create the final output. We obtain good prediction scores with a critical success index higher than 66.9 % and an obj…
ML-CIRRUS: The Airborne Experiment on Natural Cirrus and Contrail Cirrus with the High-Altitude Long-Range Research Aircraft HALO
Abstract The Midlatitude Cirrus experiment (ML-CIRRUS) deployed the High Altitude and Long Range Research Aircraft (HALO) to obtain new insights into nucleation, life cycle, and climate impact of natural cirrus and aircraft-induced contrail cirrus. Direct observations of cirrus properties and their variability are still incomplete, currently limiting our understanding of the clouds’ impact on climate. Also, dynamical effects on clouds and feedbacks are not adequately represented in today’s weather prediction models. Here, we present the rationale, objectives, and selected scientific highlights of ML-CIRRUS using the G-550 aircraft of the German atmospheric science community. The first combi…
Formulation and test of an ice aggregation scheme for two-moment bulk microphysics schemes
A simple formulation of aggregation for 2-moment bulk microphysical models is de-rived. The solution involves the evaluation of a double integral of the collection kernelweighted with the crystal size (or mass) distribution. This quantity is to be inserted intothe differential equation for the crystal number concentration which has classical form. The double integrals are evaluated numerically for log-normal size distributions overa large range of geometric mean masses. A polynomial fit of the results is given thatyields good accuracy. Various tests of the new parameterization are described: aggre-gation as stand-alone process, in a box-model, and in 2-D simulations of a cirrostratuscloud. …
Tropopause level Rossby wave breaking in the Northern Hemisphere: a feature-based validation of the ECHAM5-HAM climate model
Breaking synoptic-scale Rossby waves (RWB) at the tropopause level are central to the daily weather evolution in the extratropics and the subtropics. RWB leads to pronounced meridional transport of heat, moisture, momentum, and chemical constituents. RWB events are manifest as elongated and narrow structures in the tropopause-level potential vorticity (PV) field. A feature-based validation approach is used to assess the representation of Northern Hemisphere RWB in present-day climate simulations carried out with the ECHAM5-HAM climate model at three different resolutions (T42L19, T63L31, and T106L31) against the ERA-40 reanalysis data set. An objective identification algorithm extracts RWB …
Spectral optical layer properties of cirrus from collocated airborne measurements – a feasibility study
Abstract. Spectral optical layer properties of cirrus are derived from simultaneous and vertically collocated measurements of spectral upward and downward solar irradiance above and below the cloud layer and concurrent in situ microphysical sampling. From the irradiance data spectral transmissivity, absorptivity, reflectivity, and cloud top albedo of the observed cirrus layer are obtained. At the same time microphysical properties of the cirrus were sampled. The close collocation of the radiative and microphysical measurements, above, beneath and inside the cirrus, is obtained by using a research aircraft (Learjet 35A) in tandem with a towed platform called AIRTOSS (AIRcraft TOwed Sensor Sh…
Balloon-borne match measurements of midlatitude cirrus clouds
Observations of high supersaturations with respect to ice inside cirrus clouds with high ice water content (> 0.01 g kg−1) and high crystal number densities (> 1 cm−3) are challenging our understanding of cloud microphysics and of climate feedback processes in the upper troposphere. However, single measurements of a cloudy air mass provide only a snapshot from which the persistence of ice supersaturation cannot be judged. We introduce here the "cirrus match technique" to obtain information about the evolution of clouds and their saturation ratio. The aim of these coordinated balloon soundings is to analyze the same air mass twice. To this end the standard radiosonde equipment is complemente…
A statistical subgrid-scale algorithm for precipitation formation in stratiform clouds in the ECHAM5 single column model
Abstract. Cloud properties are usually assumed to be homogeneous within the cloudy part of the grid-box, i.e. subgrid-scale inhomogeneities in cloud cover and/or microphysical properties are often neglected. However, precipitation formation is initiated by large particles. Thus mean values are not representative and could lead to a delayed onset of precipitation. For a more physical description of the subgrid-scale structure of clouds we introduce a new statistical sub-column algorithm to study the impact of cloud inhomogeneities on stratiform precipitation. Each model column is divided into N independent sub-columns with sub-boxes in each layer, which are completely clear or cloudy. The cl…
Upper tropospheric water vapour and its interaction with cirrus clouds as seen from IAGOS long-term routine in-situ observations
IAGOS (In-service Aircraft for a Global Observing System) performs long-term routinein situobservations of atmospheric chemical composition (O3, CO, NOx, NOy, CO2, CH4), water vapour, aerosols, clouds, and temperature on a global scale by operating compact instruments on board of passenger aircraft. The unique characteristics of the IAGOS data set originate from the global scale sampling on air traffic routes with similar instrumentation such that the observations are truly comparable and well suited for atmospheric research on a statistical basis. Here, we present the analysis of 15 months of simultaneous observations of relative humidity with respect to ice (RHice) and ice crystal number …
A case study on the formation and evolution of ice supersaturation in the vicinity of a warm conveyor belt\'s outflow region
A case study is presented on the formation and evolution of an ice-supersaturated region (ISSR) that was detected by a radiosonde in NE Germany at 06:00 UTC 29 November 2000. The ISSR was situated in the vicinity of the outflow region of a warm conveyor belt associated with an intense event of cyclogenesis in the eastern North Atlantic. Using ECMWF analyses and trajectory calculations it is determined when the air parcels became supersaturated and later subsaturated again. In the case considered, the state of air parcel supersaturation can last for longer than 24h. The ISSR was unusually thick: while the mean vertical extension of ISSRs in NE Germany is about 500m, the one investigated here…
Explicit cloud-top entrainment parameterization in the global climate model ECHAM5-HAM
New developments in the turbulence parameterization in the general circulation model ECHAM5-HAM are presented. They consist mainly of an explicit entrainment closure at the top of stratocumulus-capped boundary layers and the addition of an explicit contribution of the radiative divergence in the buoyancy production term. The impact of the new implementations on a single column model study and on the global scale is presented here. The parameterization has a "smoothing" effect: the abnormally high values of turbulence kinetic energy are reduced, both in the single column and in the Californian stratocumulus region. A sensitivity study with prescribed droplet concentration shows a reduction i…
Ice supersaturated regions: properties and validation of ERA-Interim reanalysis with IAGOS in situ water vapour measurements
Cirrus clouds and their potential formation regions, so-called ice supersaturated regions (ISSRs), with values of relative humidity with respect to ice exceeding 100 %, occur frequently in the tropopause region. It is assumed that ISSRs and cirrus clouds can change the tropopause structure by diabatic processes, driven by latent heating due to phase transition and interaction with radiation. For many research questions, a three-dimensional picture including a sufficient temporal resolution of the water vapour fields in the tropopause region is required. This requirement is fulfilled nowadays by reanalysis products such as the European Centre for Medium-Range Weather Forecasts (ECMWF)…
Algorithmic differentiation for cloud schemes (IFS Cy43r3) using CoDiPack (v1.8.1)
Abstract. Numerical models in atmospheric sciences not only need to approximate the flow equations on a suitable computational grid, they also need to include subgrid effects of many non-resolved physical processes. Among others, the formation and evolution of cloud particles is an example of such subgrid processes. Moreover, to date there is no universal mathematical description of a cloud, hence many cloud schemes have been proposed and these schemes typically contain several uncertain parameters. In this study, we propose the use of algorithmic differentiation (AD) as a method to identify parameters within the cloud scheme, to which the output of the cloud scheme is most sensitive. We il…
Influence of heterogeneous freezing on the microphysical and radiative properties of orographic cirrus clouds
The influence of heterogeneous freezing on the microphysical and optical properties of orographic cirrus clouds has been simulated with the large eddy simulation model EULAG. Idealised simulations with different concentrations of ice nuclei (IN) in a dynamically dominated regime with high vertical velocities have been performed. Furthermore the temperature at cloud formation as well as the critical supersaturation for initiation of heterogenous freezing have been varied. The shortwave, longwave and net cloud forcing has been calculated under the assumption that the clouds form between 06:00 and 12:00 local time (LT) or between 12:00 and 18:00 LT. In general it can be seen that the onset of …
A tandem approach for collocated in-situ measurements of microphysical and radiative cirrus properties
Abstract. Microphysical and radiation measurements were collected with the unique AIRcraft TOwed Sensor Shuttle (AIRTOSS) – Learjet tandem platform. It is a combination of a Learjet 35A research aircraft and an instrumented aerodynamic body, which can be detached from and retracted back to the aircraft during flight. Both platforms are equipped with radiative, cloud microphysical, trace gas (CO, N2O, O3 and H2O) and meteorological instruments to study the inhomogeneity of cirrus as well as other layer clouds. Sophisticated numerical flow simulations were conducted in advance in order to optimally integrate a Cloud Combination Probe (CCP) inside the AIRTOSS. The tandem platform was used for …
Automatic shape detection of ice crystals
Abstract Clouds have a crucial impact on the energy balance of the Earth-Atmosphere system. They can cool the system by partly reflecting or scattering of the incoming solar radiation (albedo effect); moreover, thermal radiation as emitted from the Earth's surface can be absorbed and partly re-emitted by clouds leading to a warming of the atmosphere (greenhouse effect). The effectiveness of both effects crucially depends on the size and the shape of a cloud's particulate constituents, i.e. liquid water droplets or solid ice crystals. For studying cloud microphysics, in situ measurements on board of aircraft are commonly used. An important class of measurement techniques comprises optical ar…
A tandem approach for collocated measurements of microphysical and radiative cirrus properties
Microphysical and radiation measurements were collected with the novel AIRcraft TOwed Sensor Shuttle (AIRTOSS) – Learjet tandem platform. The platform is a combination of an instrumented Learjet 35A research aircraft and an aerodynamic bird, which is detached from and retracted back to the aircraft during flight via a steel wire with a length of 4000 m. Both platforms are equipped with radiative, cloud microphysical, trace gas, and meteorological instruments. The purpose of the development of this tandem set-up is to study the inhomogeneity of cirrus as well as other stratiform clouds. Sophisticated numerical flow simulations were conducted in order to optimally integrate an axially asymmet…
Impact of formulations of the homogeneous nucleation rate on ice nucleation events in cirrus
Abstract. Homogeneous freezing of solution droplets is an important pathway of ice formation in the tropopause region. The nucleation rate can be parameterized as a function of water activity, based on empirical fits and some assumptions on the underlying properties of super-cooled water, although a general theory is missing. It is not clear how nucleation events are influenced by the exact formulation of the nucleation rate or even their inherent uncertainty. In this study we investigate the formulation of the nucleation rate of homogeneous freezing of solution droplets (1) to link the formulation to the nucleation rate of pure water droplets, (2) to derive a robust and simple formulation …
A trajectory-based classification of ERA-Interim ice clouds in the region of the North Atlantic storm track
A two-type classification of ice clouds (cirrus) is introduced, based on the liquid and ice water content, LWC and IWC, along air parcel backward trajectories from the clouds. In situ cirrus has no LWC along the trajectory segment containing IWC; it forms via nucleation from the gas phase. In contrast, liquid-origin cirrus has both LWC and IWC along their backward trajectories; it forms via lifting from the lower troposphere and freezing of mixed-phase clouds. This classification is applied to 12 years of ERA-Interim ice clouds in the North Atlantic region. Between 400 and 500 hPa more than 50% are liquid-origin cirrus, whereas this frequency decreases strongly with altitude (<10% at 200 hP…
Ice supersaturations and cirrus cloud crystal numbers
Upper tropospheric observations outside and inside of cirrus clouds indicate water vapour mixing ratios sometimes exceeding water saturation. Relative humidities over ice (RHice) of up to and more than 200% have been reported from aircraft and balloon measurements in recent years. From these observations a lively discussion continues on whether there is a lack of understanding of ice cloud microphysics or whether the water measurements are tainted with large uncertainties or flaws. Here, RHice in clear air and in ice clouds is investigated. Strict quality-checked aircraft in situ observations of RHice were performed during 28 flights in tropical, mid-latitude and Arctic field experiments in…
Algorithmic Differentiation for Cloud Schemes
&lt;p&gt;Numerical models in atmospheric sciences do not only need to approximate the flow equations on a suitable computational grid, they also need to include subgrid effects of many non-resolved physical processes. Among others, the formation and evolution of cloud particles is an example of such subgrid processes. Moreover, to date there is no universal mathematical description of a cloud, hence many cloud schemes were proposed and these schemes typically contain several uncertain parameters. In this study, we propose the use of algorithmic differentiation (AD) as a method to identify parameters within the cloud scheme, to which the output of the cloud scheme is most sensitive.…
Automated detection and classification of synoptic scale fronts from atmospheric data grids
&lt;p&gt;Automatic determination of fronts from atmospheric data is an important task for weather prediction as well as for research of synoptic scale phenomena. We developed a deep neural network to detect and classify fronts from multi-level ERA5 reanalysis data. Model training and prediction is evaluated using two different regions covering Europe and North America with data from two weather services. Due to a label deformation step performed during training we are able to directly generate frontal lines with no further thinning during post processing. Our network compares well against the weather service labels with a Critical Success Index higher than 66.9% and a Object Detecti…
Microphysical and radiative changes in cirrus clouds by geoengineering the stratosphere
[1] In the absence of tangible progress in reducing greenhouse gas emissions, the implementation of solar radiation management has been suggested as measure to stop global warming. Here we investigate the impacts on northern midlatitude cirrus from continuous SO2emissions of 2–10 Mt/a in the tropical stratosphere. Transport of geoengineering aerosols into the troposphere was calculated along trajectories based on ERA Interim reanalyses using ozone concentrations to quantify the degree of mixing of stratospheric and tropospheric air termed “troposphericity”. Modeled size distributions of the geoengineered H2SO4-H2O droplets have been fed into a cirrus box model with spectral microphysics. Th…