Aircraft Observations Research Articles

Consistency evaluation of tropospheric ozone from ozonesonde and IAGOS (In-service Aircraft for a Global Observing System) observations: vertical distribution, ozonesonde types, and station–airport distance

Abstract. The vertical distribution of tropospheric O3 from ozonesondes is compared with that from In-service Aircraft for a Global Observing System (IAGOS) measurements collected at 23 pairs of sites between about 30° S and 55° N from 1995 to 2021. Profiles of tropospheric O3 from IAGOS are generally in good agreement with ozonesonde observations from electrochemical concentration cells (ECCs), Brewer–Mast sondes, and carbon–iodine sensors, with average biases of 2.58, −0.28, and 0.67 ppb and correlation coefficients (R) of 0.72, 0.82, and 0.66, respectively. Agreement between aircraft and Indian-sonde observations is poor, with an average bias of 15.32 ppb and an R value of 0.44. The O3 concentration observed by ECC sondes is, on average, 5 %–10 % higher than that observed by IAGOS, and the relative bias increases modestly with altitude. For other sonde types, there are some seasonal and altitudinal variations in the relative bias with respect to the IAGOS measurements, but these appear to be caused by local differences. The distance between the station and airport, when within 4° (latitude and longitude), has little effect on the comparison results. For the ECC ozonesondes, the overall bias with respect to the IAGOS measurements varies from 5.7 to 9.8 ppb when the station pairs are grouped by station–airport distances of <1° (latitude and longitude), 1–2°, and 2–4°. Correlations for these groups correspond to R=0.8, 0.9, and 0.7. These comparison results provide important information for merging ozonesonde and IAGOS measurement datasets. They can also be used to evaluate the relative biases of different sonde types in the troposphere, using the aircraft as a transfer standard.

Troposphere–stratosphere-integrated bromine monoxide (BrO) profile retrieval over the central Pacific Ocean

Abstract. Bromine monoxide (BrO) is relevant to atmospheric oxidative capacity, affecting the lifetime of greenhouse gases (i.e., methane, dimethylsulfide) and mercury oxidation. However, measurements of BrO radical vertical profiles are rare, and BrO is highly variable. As a result, the few available aircraft observations in different regions of the atmosphere are not easily reconciled. Autonomous multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments placed at remote mountaintop observatories (MT-DOAS) present a cost-effective alternative to aircraft, with the potential to probe the climate-relevant yet understudied free troposphere more routinely. Here, we describe an innovative full-atmosphere BrO and formaldehyde (HCHO) profile retrieval algorithm using MT-DOAS measurements at Mauna Loa Observatory (MLO – 19.536° N, 155.577° W; 3401 m a.s.l.). The retrieval is based on time-dependent optimal estimation and simultaneously inverts 190+ individual BrO (and formaldehyde, HCHO) SCDs (slant column densities; SCD = dSCD + SCDRef) from solar stray light spectra measured in the zenith and off-axis geometries at high and low solar zenith angles (92° > SZA > 30°) to derive BrO concentration profiles from 1.9 to 35 km with 7.5 degrees of freedom (DoFs). Two case study days are characterized by the absence (26 April 2017, base case) and presence of a Rossby-wave-breaking double tropopause (29 April 2017, RW-DT case). Stratospheric-BrO vertical columns are nearly identical on both days (VCD = (1.5 ± 0.2) × 1013 molec. cm−2), and the stratospheric-BrO profile peaks at a lower altitude during the RW-DT (1.6–2.0 DoFs). Tropospheric-BrO VCDs increase from (0.70 ± 0.14) × 1013 molec. cm−2 (base case) to (1.00 ± 0.14) × 1013 molec. cm−2 (RW-DT) owing to a 3-fold increase in BrO in the upper troposphere (1.7–1.9 DoFs). BrO at MLO increases from (0.23 ± 0.03) pptv (base case) to (0.46 ± 0.03) pptv (RW-DT) and is characterized by an added time resolution (∼ 3.8 DoFs). Up to (0.9 ± 0.1) pptv BrO is observed above MLO in the lower free troposphere in the absence of the double tropopause. We validate the retrieval using aircraft BrO profiles and in situ HCHO measurements aboard the NSF/NCAR GV aircraft above MLO (11 January 2014) that establish BrO peaks around 2.4 pptv above 13 km in the upper troposphere–lower stratosphere (UTLS) during a similar RW-DT event (0.83 × 1013 molec. cm2 tropospheric-BrO VCD above 2 km). The tropospheric-BrO profile measured using MT-DOAS (RW-DT case) and using the aircraft agree well (after averaging-kernel smoothing). Furthermore, these tropospheric-BrO profiles over the central Pacific Ocean are found to closely resemble those over the eastern Pacific Ocean (2–14 km) and are in contrast to those over the western Pacific Ocean, where a C-shaped tropospheric-BrO profile shape has been observed.

Analysis of noctilucent clouds’ fields according to ground-based network and airborne photography data

The article analyzes the fields of noctilucent clouds over the territory of the Russian Federation, recorded by a ground-based network of cameras using also aircraft photography, over the two nights in June 2021. It is demonstrated that aircraft photography can significantly improve the coverage of the territory of noctilucent clouds’ probable appearance. The detected noctilucent cloud fields are compared with model regions of water vapor condensation derived from satellite measurements of temperature and water vapor mixing ratio. Practical steps are proposed for the development of aircraft observations of noctilucent clouds.

Aircraft Observations Reveal the Relationship Between Cumulus Entrainment Rate and Aerosol Loading

AbstractThe influence of entrainment, a key process characterized by the entrainment rate in cumulus parameterization, on aerosol‐cloud interactions has been widely recognized. However, despite qualitative links established between entrainment and aerosol loading, a quantitative relationship based on observational evidence remains elusive. This study utilizes aircraft observations of cumulus clouds during two field campaigns to determine the quantitative relationship between entrainment rate and aerosol loading. In both campaigns, the entrainment rate is negatively correlated with aerosol loading. It is speculated that increased aerosol loading enhances cloud edge droplet evaporation, which leads to increased buoyancy and vertical velocity within the cloud, thereby reducing the entrainment rate. Further analysis shows that the response of entrainment rate to aerosol perturbations is more significant in smaller cumulus clouds with weak buoyancy and less pronounced under opposite conditions. These findings shed new light on improving the description of aerosol‐cloud interactions in cumulus parameterizations.

The Dual Nature of Entrainment-Mixing Signatures Revealed Through Large-Eddy Simulations of a Convection-Cloud Chamber

Abstract Entrainment of subsaturated air into a cloud can influence its optical and microphysical properties in various ways, depending on the droplet evaporation and turbulent mixing time scales. Previous experiments in the Pi Convection-Cloud Chamber have revealed that, given a fixed entrained air property, the mixing of entrained subsaturated air results in complete evaporation of some cloud droplets, with the rest remaining unchanged. This is a signature of inhomogeneous mixing. While comparing the results of entrainment with varying air properties, the mixing signature appears as if the subsaturated air is well mixed with the cloud to evenly reduce the droplets’ size. In other words, taken together, the experiments appear to have the signature of homogeneous mixing. To explore these results in a greater depth, we conduct large-eddy simulations combined with a bin microphysics scheme. Our results reproduce the similar signatures of inhomogeneous and homogeneous mixing, implying that LES can resolve the inhomogeneous mixing when the grid spacing is smaller than the entrained air parcel. Additionally, we observe that increasing aerosol injection rate enhances the signature of inhomogeneous mixing, while coarser grid spacing diminishes it. Finally, the change in wall fluxes in response to various entrained air properties confirms that the homogeneous signature seen in the analysis of an ensemble of simulations is the result of various equilibrium states. This further strengthens the suggestion that the homogeneous mixing signature found in aircraft observations near the cloud top may result from combining entrainment events of different intensities, possibly caused by various-sized eddies.

Evidence of a New Population of Weak Terrestrial Gamma‐Ray Flashes Observed From Aircraft Altitude

AbstractTerrestrial Gamma‐ray Flashes (TGFs) are ten‐to‐hundreds of microsecond bursts of gamma‐rays produced when electrons in strong electric fields in thunderclouds are accelerated to relativistic energies. Space instruments have observed TGFs with source photon brightness down to ∼1017–1016. Based on space and aircraft observations, TGFs have been considered rare phenomena produced in association with very few lightning discharges. Space observations associated with lightning ground observations in the radio band have indicated that there exists a population of dimmer TGFs. Here we show observations of TGFs from aircraft altitude that were not detected by a space instrument viewing the same area. The TGFs were found through Monte Carlo modeling to be associated with 1015–1012 photons at source, which is several orders of magnitude below what can be seen from space. Our results suggest that there exists a significant population of TGFs that are too weak to be observed from space.

Evaluating CHASER V4.0 global formaldehyde (HCHO) simulations using satellite, aircraft, and ground-based remote-sensing observations

Abstract. Formaldehyde (HCHO), a precursor to tropospheric ozone, is an important tracer of volatile organic compounds (VOCs) in the atmosphere. Two years (2019–2020) of HCHO simulations obtained from the global chemistry transport model CHASER at a horizontal resolution of 2.8° × 2.8° have been evaluated using the Tropospheric Monitoring Instrument (TROPOMI) and multi-axis differential optical absorption spectroscopy (MAX-DOAS) observations. In situ measurements from the Atmospheric Tomography Mission (ATom) in 2018 were used to evaluate the HCHO simulations for 2018. CHASER reproduced the TROPOMI-observed global HCHO spatial distribution with a spatial correlation (r) of 0.93 and a negative bias of 7 %. The model showed a good capability to reproduce the observed magnitude of the HCHO seasonality in different regions, including the background conditions. The discrepancies between the model and satellite in the Asian regions were related mainly to the underestimated and missing anthropogenic emission inventories. The maximum difference between two HCHO simulations based on two different nitrogen oxide (NOx) emission inventories was 20 %. TROPOMI's finer spatial resolution than that of the Ozone Monitoring Instrument (OMI) sensor reduced the global model–satellite root-mean-square error (RMSE) by 20 %. The OMI- and TROPOMI-observed seasonal variations in HCHO abundances were consistent. The simulated seasonality showed better agreement with TROPOMI in most regions. The simulated HCHO and isoprene profiles correlated strongly (R=0.81) with the ATom observations. However, CHASER overestimated HCHO mixing ratios over dense vegetation areas in South America and the remote Pacific region (background condition), mainly within the planetary boundary layer (< 2 km). The simulated seasonal variations in the HCHO columns showed good agreement (R>0.70) with the MAX-DOAS observations and agreed within the 1σ standard deviation of the observed values. However, the temporal correlation (R∼0.40) was moderate on a daily scale. CHASER underestimated the HCHO levels at all sites, and the peak occurrences in the observed and simulated HCHO seasonality differed. The coarseness of the model's resolution could potentially lead to such discrepancies. Sensitivity studies showed that anthropogenic emissions were the highest contributor (up to ∼ 35 %) to the wintertime regional HCHO levels.

The algorithm of microphysical-parameter profiles of aerosol and small cloud droplets based on the dual-wavelength lidar data

Abstract. This study proposed an inversion method for atmospheric-aerosol or cloud microphysical parameters based on dual-wavelength lidar data. The matching characteristics between aerosol and cloud particle size distributions and gamma distributions were studied using aircraft observation data. The feasibility of the retrieval of the particle effective radius from lidar ratios and backscatter ratios was simulated and studied. A method for inverting the effective radius and number concentration of atmospheric aerosols or small cloud droplets using the backscatter ratio was proposed, and the error sources and applicability of the algorithm were analyzed. This algorithm was suitable for the inversion of uniformly mixed and single-property aerosol layers or small cloud droplets. Compared with the previous study, this algorithm could quickly obtain the microphysical parameters of atmospheric particles and has good robustness. For aerosol particles, the inversion range that this algorithm can achieve is 0.3–1.7 µm. For cloud droplets, it is 1.0–10 µm. An atmospheric-observation experiment was conducted using the multi-wavelength lidar developed by Xi'an University of Technology, and a thin cloud layer was captured. The microphysical parameters of aerosol and clouds during this process were retrieved. The results clearly demonstrate the growth of the effective radius and number concentration.

Pathways of ice multiplication in nimbostratus clouds during the Indian summer monsoon

The present study illustrates the microphysical parameters of nimbostratus clouds during the Indian summer monsoon using airborne and radar observations conducted as part of the Cloud Aerosol Interaction and Precipitation Enhancement Experiment (CAIPEEX), and simulations with the Weather Research and Forecasting (WRF) model. The study also analyzes the impact of secondary ice production (emphasizing the microphysical pathways of Hallett-Mossop (HM) as represented in the models) on cloud processes using model sensitivity simulations. The effect of possible pathways of riming during the HM process is the focus of the sensitivity simulations. Aircraft observations showed higher ice particle concentrations in the Hallett–Mossop zone (−3 °C to - 8 °C) along with the existence of smaller and larger cloud droplets, rimed needles, columns, and snow particles. Observations strongly suggested the active presence of the HM process in this cloud. The observed mean values of microphysical parameters including liquid water content, ice number concentrations, and maximum vertical velocity, agreed well with model simulations. The number concentration and water content of ice crystals and snow decreased in the HM zone during the HM inactive case. A reduction in rainfall is also observed in the absence of HM. HM increases convection as it causes a sudden increase in the number concentration of small ice particles, which causes rainwater to freeze quickly and water vapour to be consumed by diffusional growth. Latent heating is produced by both effects, which energizes the convection.

Evaluation of a New Approach for Entrainment and Detrainment Rate Estimation

AbstractEntrainment and detrainment rates (ε and δ) constitute the most critical free parameters in mass flux schemes commonly employed for cumulus parameterizations. Recently, Zhu et al. (2021) introduced a new approach that utilizes aircraft observations to simultaneously estimate ε and δ for cumulus clouds, overcoming the limitation of other observation‐based approaches that solely yield ε without offering insights into δ. This study aims to comprehensively evaluate the reliability of this new approach. First, evaluation using an Explicit Mixing Parcel Model demonstrates the capability of the new approach to back‐calculate predetermined ε and δ based on the physical properties before and after the entrainment mixing. Second, evaluation using large‐eddy simulations illustrates that the new approach yields consistent ε and δ profiles compared to the traditional approach. Sensitivity tests indicate a weak sensitivity of the estimated δ with the new approach to the entrained air source. A decrease in the proportion of cloudy air in the assumed detrained air leads to a reduction in the estimated δ, while ε remains unaffected. Finally, the most appropriate assumptions for entrained and detrained air are discussed. Estimating ε for cumulus parameterizations involves acquiring ambient air more than 500 m away from the cloud edge as entrained air. Due to implicit mean field approximations in the traditional approach, determining the optimal assumption for detrained air properties proves challenging. This study confirms the reliability of the new approach in estimating ε and δ, providing confidence in its application to extensive observational data and advancement in parameterization.

A global surface CO2 flux dataset (2015–2022) inferred from OCO-2 retrievals using the GONGGA inversion system

Abstract. Accurate assessment of the size and distribution of carbon dioxide (CO2) sources and sinks is important for efforts to understand the carbon cycle and support policy decisions regarding climate mitigation actions. Satellite retrievals of the column-averaged dry-air mole fractions of CO2 (XCO2) have been widely used to infer spatial and temporal variations in carbon fluxes through atmospheric inversion techniques. In this study, we present a global spatially resolved terrestrial and ocean carbon flux dataset for 2015–2022. The dataset was generated by the Global ObservatioN-based system for monitoring Greenhouse GAses (GONGGA) atmospheric inversion system through the assimilation of Orbiting Carbon Observatory-2 (OCO-2) XCO2 retrievals. We describe the carbon budget, interannual variability, and seasonal cycle for the global scale and a set of TransCom regions. The 8-year mean net biosphere exchange and ocean carbon fluxes were −2.22 ± 0.75 and −2.32 ± 0.18 Pg C yr−1, absorbing approximately 23 % and 24 % of contemporary fossil fuel CO2 emissions, respectively. The annual mean global atmospheric CO2 growth rate was 5.17 ± 0.68 Pg C yr−1, which is consistent with the National Oceanic and Atmospheric Administration (NOAA) measurement (5.24 ± 0.59 Pg C yr−1). Europe has the largest terrestrial sink among the 11 TransCom land regions, followed by Boreal Asia and Temperate Asia. The dataset was evaluated by comparing posterior CO2 simulations with Total Carbon Column Observing Network (TCCON) retrievals as well as Observation Package (ObsPack) surface flask observations and aircraft observations. Compared with CO2 simulations using the unoptimized fluxes, the bias and root mean square error (RMSE) in posterior CO2 simulations were largely reduced across the full range of locations, confirming that the GONGGA system improves the estimates of spatial and temporal variations in carbon fluxes by assimilating OCO-2 XCO2 data. This dataset will improve the broader understanding of global carbon cycle dynamics and their response to climate change. The dataset can be accessed at https://doi.org/10.5281/zenodo.8368846 (Jin et al., 2023a).

Third‐Order Structure Functions of Zonal Winds in the Thermosphere Using CHAMP and GOCE Observations

AbstractWe use multi‐year observations of cross‐track winds (u) from the CHAllenging Minisatellite Payload (CHAMP) and the Gravity Field and Steady State Ocean Circulation Explorer (GOCE) to calculate third‐order structure functions in the thermosphere as a function of horizontal separation (s). They are computed using the mean (〈δu3〉) and the median and implemented over non‐polar satellite paths in both hemispheres. On height averages, 〈δu3〉 is shown to scale with s2 for s ≃ 80–1,000 km, in agreement with equivalent estimates in the lower atmosphere from aircraft observations. Conversely, follows an s3 power law for almost the whole s range, consistent with the two‐dimensional turbulence scaling law for a direct enstrophy cascade. These scaling laws appear independent of winds in distinct atmospheric regions. Furthermore, the functions are predominantly positive, indicating a preferential cyclonic motion for the wind.

Lightning NOx in the 29–30 May 2012 Deep Convective Clouds and Chemistry (DC3) Severe Storm and Its Downwind Chemical Consequences

AbstractA cloud‐resolved storm and chemistry simulation of a severe convective system in Oklahoma constrained by anvil aircraft observations of NOx was used to estimate the mean production of NOx per flash in this storm. An upward ice flux scheme was used to parameterize flash rates in the model. Model lightning was also constrained by observed lightning flash types and the altitude distribution of flash channel segments. The best estimate of mean NOx production by lightning in this storm was 80–110 mol per flash, which is smaller than many other literature estimates. This result is likely due to the storm having been a high flash rate event in which flash extents were relatively small. Over the evolution of this storm a moderate negative correlation was found between the total flash rate and flash extent and energy per flash. A longer‐term simulation at 36‐km horizontal resolution with parameterized convection was used to simulate the downwind transport and chemistry of the anvil outflow from the same storm. Convective transport of low‐ozone air from the boundary layer decreased ozone in the anvil outflow by up to 20–40 ppbv compared with the initial conditions, which contained stratospheric influence. Photochemical ozone production in the lightning‐NOx enhanced convective plume proceeded at a rate of 10–11 ppbv per day in the 9–11 km outflow layer over the 24‐hr period of downwind transport to the Southern Appalachians. Photochemical production plays a large role in the restoration of upper tropospheric ozone following deep convection.

Combined assimilation of NOAA surface and MIPAS satellite observations to constrain the global budget of carbonyl sulfide

Abstract. Carbonyl sulfide (COS), a trace gas in our atmosphere that leads to the formation of aerosols in the stratosphere, is largely taken up by terrestrial ecosystems. Quantifying the biosphere uptake of COS could provide a useful quantity to estimate gross primary productivity (GPP). Some COS sources and sinks still contain large uncertainties, and several top-down estimates of the COS budget point to an underestimation of sources, especially in the tropics. We extended the inverse model TM5-4DVAR to assimilate Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) satellite data, in addition to National Oceanic and Atmospheric Administration (NOAA) surface data as used in a previous study. To resolve possible discrepancies among the two observational data sets, a bias correction scheme is necessary and implemented. A set of inversions is presented that explores the influence of the different measurement streams and the settings of the prior fluxes. To evaluate the performance of the inverse system, the HIAPER Pole-to-Pole Observations (HIPPO) aircraft observations and NOAA airborne profiles are used. All inversions reduce the COS biosphere uptake from a prior value of 1053 GgS a−1 to much smaller values, depending on the inversion settings. These large adjustments of the biosphere uptake often turn parts of Amazonia into a COS source. Only inversions that exclusively use MIPAS observations, or strongly reduce the prior errors on the biosphere flux, maintain the Amazon as a COS sink. Inclusion of MIPAS data in the inversion leads to a better separation of land and ocean fluxes. Over the Amazon, these inversions reduce the biosphere uptake from roughly 300 to 100 GgS a−1, indicating a strongly overestimated prior uptake in this region. Although a recent study also reported reduced COS uptake over the Amazon, we emphasise that a careful construction of prior fluxes and their associated errors remains important. For instance, an inversion that gives large freedom to adjust the anthropogenic and ocean fluxes of CS2, an important COS precursor, also closes the budget satisfactorily with much smaller adjustments to the biosphere. We achieved better characterisation of biosphere prior and uncertainty, better characterisation of combined ocean and land fluxes, and better constraint of both by combining surface and satellite observations. We recommend more COS observations to characterise biosphere and ocean fluxes, especially over the data-poor tropics.

Spatial downscaling of SMAP soil moisture estimation using multiscale geographically weighted regression during SMAPVEX16

Passive microwave remote sensing missions such as Soil Moisture Active Passive (SMAP) are widely used for global soil moisture (SM) estimation. But the low spatial resolution of passive microwave SM products limits their applications at regional and local scales (normally 1–10 km). Using high-resolution auxiliary data to downscale passive microwave SM products is the primary way to obtain high-resolution SM data. However, the existing SM downscaling methods do not entirely account for the scale differences in the impact of various variables on the distribution of SM. To solve this problem, this paper introduced MGWR (Multiscale Geographically Weighted Regression) as a novel way to analyze the scale differences of normalized difference vegetation index (NDVI) and land surface temperature (LST) on the spatial pattern of SM. Accordingly, a new spatial downscaling algorithm for SMAP SM products is proposed and compared with the geographically weighted regression (GWR) algorithm and moving window algorithm (MW). The algorithm was applied to data taken during the SMAP Validation Experiment 2016 (SMAPVEX16), and the downscaled SM evaluated with in situ SM and aircraft observations. These results show that: 1) the SM downscaling conversion function, based on the MGWR, effectively reveals how different surface parameters relate to SM, with LST influencing SM globally and NDVI affecting it locally, 2) compared with GWR and MW, the 1 km SM obtained by the MGWR-based downscaling method showed better spatial agreement with airborne SM in the semiarid region with the ubRMSE decreasing by 0.023 m3/m3 and 0.026 m3/m3 for GWR and MW respectively, 3) when evaluating the 1 km SM with in situ SM in the semiarid region, the SM obtained through MGWR exhibited lower ubRMSE values as compared to both GWR and MW. Specifically, the median ubRMSE values for MGWR, GWR, and MW were 0.043 m3/m3, 0.053 m3/m3, and 0.094 m3/m3, respectively. In conclusion, this study demonstrates that MGWR can effectively improve the accuracy of downscaled SM in semiarid regions.

Integration of surface-based and space-based atmospheric CO2 measurements for improving carbon flux estimates using a new developed 3-GAS inversion model

The ongoing necessity and developmental trajectory in atmospheric CO2 assimilation system revolve around enhancing the optimization efficacy of assimilation inversions through the synergistic integration of distinct observational datasets. In this study, we developed an new inversion system-Global Greenhouse-Gases Assimilation System (3-GAS), by coupling the global atmospheric transport model GEOS-Chem with the four-dimensional local ensemble transform Kalman filter (4D-LETKF) algorithm. We compiled an integrated dataset utilizing the widely recognized carbon observation data from ObsPack (Observation Package), GOSAT (The Greenhouse Gases Observing Satellite), and OCO-2 (Orbiting Carbon Observatory-2) to estimate global surface fluxes. The integration of multiple observational data significantly amplifies spatial coverage, enhances the variability in flux inversion, and fosters improvements in localized assimilation efficacy. The spatial distribution of the carbon flux from the combined assimilation emerges as a comprehensive outcome, born from the assimilation with individual observational data. This distribution is primarily influenced by the coverage attributes of the original observations and the meticulous selection of observations during the assimilation process. The validation with surface and aircraft observations reveals that the results of combined assimilation possess composite advantages, which are manifested in alleviating certain intrinsic errors inherent to individual assimilation. The combined assimilation improves the accuracy based on root-mean-square error reduced by 9%, 22% and 33% in comparison to surface-only, GOSAT-only and OCO-2-only assimilation. Eventually, predicated upon preliminary inversion results and regional characteristics, we endeavor to propose a more refined filtration scheme and error determination method for the integrated dataset, as the initiative for the efficient assimilation of multi-source data.

Deep Convective Microphysics Experiment (DCMEX) coordinated aircraft and ground observations: microphysics, aerosol, and dynamics during cumulonimbus development

Abstract. Cloud feedbacks associated with deep convective anvils remain highly uncertain. In part, this uncertainty arises from a lack of understanding of how microphysical processes influence the cloud radiative effect. In particular, climate models have a poor representation of microphysics processes, thereby encouraging the collection and study of observation data to enable better representation of these processes in models. As such, the Deep Convective Microphysics Experiment (DCMEX) undertook an in situ aircraft and ground-based measurement campaign of New Mexico deep convective clouds during July–August 2022. The campaign coordinated a broad range of instrumentation measuring aerosol, cloud physics, radar, thermodynamics, dynamics, electric fields, and weather. This paper introduces the potential data user to DCMEX observational campaign characteristics, relevant instrument details, and references to more detailed instrument descriptions. Also included is information on the structure and important files in the dataset in order to aid the accessibility of the dataset to new users. Our overview of the campaign cases illustrates the complementary operational observations available and demonstrates the breadth of the campaign cases observed. During the campaign, a wide selection of environmental conditions occurred, ranging from dry, northerly air masses with low wind shear to moist, southerly air masses with high wind shear. This provided a wide range of different convective growth situations. Of 19 flight days, only 2 d lacked the formation of convective cloud. The dataset presented (https://doi.org/10.5285/B1211AD185E24B488D41DD98F957506C; Facility for Airborne Atmospheric Measurements et al., 2024) will help establish a new understanding of processes on the smallest cloud- and aerosol-particle scales and, once combined with operational satellite observations and modelling, can support efforts to reduce the uncertainty of anvil cloud radiative impacts on climate scales.

High Elevation Radiation Array (HERA) detectors for airborne thunderstorm investigations

A high-energy atmospheric physics phenomenon, referred to as a terrestrial gamma ray flash (TGF), is associated with lightning and produces large bursts of energetic photon radiation. TGFs will be investigated using a suite of gamma-ray instruments designed and constructed to fly on ten United States Air Force (USAF) WC-130J Hurricane Hunter aircraft as part of an aircrew ionization study led by the Air Force Institute of Technology (AFIT) and the United States Air Force School of Aerospace Medicine (USAFSAM), in cooperation with the 53rd Weather Reconnaissance Squadron (WRS). Each instrument consists of one NaI and one plastic detector, a GPS timing device, and an instrument computer that performs data acquisition. High Elevation Radiation Array (HERA) detectors will be employed to maximize the chances of observing TGFs near their source and to gain a better understanding of their origin, mechanism, ubiquity, and to assess potential hazards posed to military and commercial aircrew and passengers. The HERA program, deployed on 10 separate Air Force aircraft over a multi-year campaign, will result in thousands of observational flight hours and be the largest concerted effort to date to observe TGFs in situ through aircraft observations. In this paper, we give an overview of the scientific goals of this campaign and how the HERA instruments have been designed to meet those goals. We include a detailed description of the HERA instrument, along with mass model and signal processing simulations.

Aircraft observations of aerosol and BC during the East Asian dust storm event: Vertical profiles, size distribution and mixing state

The vertical structure of meteorological elements, black carbon (BC), and aerosols were in situ characterized by the King Air 350 aircraft measurements over the northeast of China during an East Asian dust storm event from April 10 to 13, 2023. The main causes of this dust storm event were the strong winds generated by the Mongolian cyclone, the cold front on the ground provided dynamic conditions, and the unstable atmospheric state provided thermal conditions and less precipitation in the dust source areas led to exposed the ground surface. The number concentrations of BC and aerosol had significant vertical distribution. The number concentration of aerosol in the dust storm (DS) layer increased rapidly with height, while the number and mass concentrations of BC were higher and varied slightly with height. The BC core mass median diameters (MMD) first increased with height, reaching a maximum value of 0.24 μm at 3000 m; and it rapidly decreased to 0.18 μm at 3200 m and then varied slightly with height above 3200 m. BC's bulk relative coating thickness (coated diameter divided by uncoated core diameter Dp/Dc) first decreased with height, reaching a minimum value of 1.16 at 3000 m, then increased with height, reaching a maximum of 1.38 (4250 m). A wind shear layer (7.9 m s−1 at 1850 m) with a sudden drop in temperature (the cooling rate can reach 7.4 °C·(100 m)−1) at 1900 m hinders the transportation of dust particles to the ground, allowing the dust particles to be transported at long-term distances. There were two stratifications in the dust storm layer: DS1 (1700–3000 m, the vertical wind speed was downward, the horizontal wind speed was small) and DS2 (3000–4500 m, the vertical wind speed was upward, the horizontal wind speed was large). The distinct atmospheric flow field conditions between DS1 and DS2 led to different BC and aerosol pollution characteristics. In DS2, the number concentrations of BC and aerosol were higher, the aging ratio of BC particles was higher, and the size of the core particle of BC was smaller. The concentrations of BC, aerosol, Dp/Dc, and MMD in DS2 were 1.1, 59.0, 1.02, and 0.84 times higher than those in DS1, respectively. Dust particles transported over long distances were primarily small. The number concentrations of aerosols at 0.1–0.12 μm and 0.1–0.24 μm were higher in DS1 and DS2, respectively. Still, the number concentration of aerosol at 0.12–0.28 μm in DS1 was lower than in the planetary boundary layer (PBL) and residual layer (RL).

Quantifying riming from airborne data during the HALO-(AC)3 campaign

Abstract. Riming is a key precipitation formation process in mixed-phase clouds which efficiently converts cloud liquid to ice water. Here, we present two methods to quantify riming of ice particles from airborne observations with the normalized rime mass, which is the ratio of rime mass to the mass of a size-equivalent spherical graupel particle. We use data obtained during the HALO-(AC)3 aircraft campaign, where two aircraft collected radar and in situ measurements that were closely spatially and temporally collocated over the Fram Strait west of Svalbard in spring 2022. The first method is based on an inverse optimal estimation algorithm for the retrieval of the normalized rime mass from a closure between cloud radar and in situ measurements during these collocated flight segments (combined method). The second method relies on in situ observations only, relating the normalized rime mass to optical particle shape measurements (in situ method). We find good agreement between both methods during collocated flight segments with median normalized rime masses of 0.024 and 0.021 (mean values of 0.035 and 0.033) for the combined and in situ method, respectively. Assuming that particles with a normalized rime mass smaller than 0.01 are unrimed, we obtain average rimed fractions of 88 % and 87 % over all collocated flight segments. Although in situ measurement volumes are in the range of a few cubic centimeters and are therefore much smaller than the radar volume (about 45 m footprint diameter at an altitude of 500 m above ground, with a vertical resolution of 5 m), we assume they are representative of the radar volume. When this assumption is not met due to less homogeneous conditions, discrepancies between the two methods result. We show the performance of the methods in a case study of a collocated segment of cold-air outbreak conditions and compare normalized rime mass results with meteorological and cloud parameters. We find that higher normalized rime masses correlate with streaks of higher radar reflectivity. The methods presented improve our ability to quantify riming from aircraft observations.