Southern Africa produces almost a third of the Earth's biomass burning (BB) aerosol particles, yet the fate of these particles and their influence on regional and global climate is poorly understood. ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) is a 5-year NASA EVS-2 (Earth Venture Suborbital-2) investigation with three intensive observation periods designed to study key atmospheric processes that determine the climate impacts of these aerosols. During the Southern Hemisphere winter and spring (June–October), aerosol particles reaching 3–5 km in altitude are transported westward over the southeast Atlantic, where they interact with one of the largest subtropical stratocumulus (Sc) cloud decks in the world. The representation of these interactions in climate models remains highly uncertain in part due to a scarcity of observational constraints on aerosol and cloud properties, as well as due to the parameterized treatment of physical processes. Three ORACLES deployments by the NASA P-3 aircraft in September 2016, August 2017, and October 2018 (totaling ∼350 science flight hours), augmented by the deployment of the NASA ER-2 aircraft for remote sensing in September 2016 (totaling ∼100 science flight hours), were intended to help fill this observational gap. ORACLES focuses on three fundamental science themes centered on the climate effects of African BB aerosols: (a) direct aerosol radiative effects, (b) effects of aerosol absorption on atmospheric circulation and clouds, and (c) aerosol–cloud microphysical interactions. This paper summarizes the ORACLES science objectives, describes the project implementation, provides an overview of the flights and measurements in each deployment, and highlights the integrative modeling efforts from cloud to global scales to address science objectives. Significant new findings on the vertical structure of BB aerosol physical and chemical properties, chemical aging, cloud condensation nuclei, rain and precipitation statistics, and aerosol indirect effects are emphasized, but their detailed descriptions are the subject of separate publications. The main purpose of this paper is to familiarize the broader scientific community with the ORACLES project and the dataset it produced.

The total effect of aerosols, both directly and on cloud properties, remains the biggest source of uncertainty in anthropogenic radiative forcing on the climate. Correct characterization of intensive aerosol optical properties, particularly in conditions where absorbing aerosol is present, is a crucial factor in quantifying these effects. The southeast Atlantic Ocean (SEA), with seasonal biomass burning smoke plumes overlying and mixing with a persistent stratocumulus cloud deck, offers an excellent natural laboratory to make the observations necessary to understand the complexities of aerosol–cloud–radiation interactions. The first field deployment of the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) campaign was conducted in September of 2016 out of Walvis Bay, Namibia. Data collected during ORACLES-2016 are used to derive aerosol properties from an unprecedented number of simultaneous measurement techniques over this region. Here, we present results from six of the eight independent instruments or instrument combinations, all applied to measure or retrieve aerosol absorption and single-scattering albedo. Most but not all of the biomass burning aerosol was located in the free troposphere, in relative humidities typically ranging up to 60 %. We present the single-scattering albedo (SSA), absorbing and total aerosol optical depth (AAOD and AOD), and absorption, scattering, and extinction Ångström exponents (AAE, SAE, and EAE, respectively) for specific case studies looking at near-coincident and near-colocated measurements from multiple instruments, and SSAs for the broader campaign average over the month-long deployment. For the case studies, we find that SSA agrees within the measurement uncertainties between multiple instruments, though, over all cases, there is no strong correlation between values reported by one instrument and another. We also find that agreement between the instruments is more robust at higher aerosol loading (AOD400>0.4). The campaign-wide average and range shows differences in the values measured by each instrument. We find the ORACLES-2016 campaign-average SSA at 500 nm (SSA500) to be between 0.85 and 0.88, depending on the instrument considered (4STAR, AirMSPI, or in situ measurements), with the interquartile ranges for all instruments between 0.83 and 0.89. This is consistent with previous September values reported over the region (between 0.84 and 0.90 for SSA at 550nm). The results suggest that the differences observed in the campaign-average values may be dominated by instrument-specific spatial sampling differences and the natural physical variability in aerosol conditions over the SEA, rather than fundamental methodological differences.

In the southeast Atlantic, well-defined smoke plumes from Africa advect over marine boundary layer cloud decks; both are most extensive around September, when most of the smoke resides in the free troposphere. A framework is put forth for evaluating the performance of a range of global and regional atmospheric composition models against observations made during the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) airborne mission in September 2016. A strength of the comparison is a focus on the spatial distribution of a wider range of aerosol composition and optical properties than has been done previously. The sparse airborne observations are aggregated into approximately 2^∘ grid boxes and into three vertical layers: 3–6 km, the layer from cloud top to 3 km, and the cloud-topped marine boundary layer. Simulated aerosol extensive properties suggest that the flight-day observations are reasonably representative of the regional monthly average, with systematic deviations of 30 % or less. Evaluation against observations indicates that all models have strengths and weaknesses, and there is no single model that is superior to all the others in all metrics evaluated. Whereas all six models typically place the top of the smoke layer within 0–500 m of the airborne lidar observations, the models tend to place the smoke layer bottom 300–1400 m lower than the observations. A spatial pattern emerges, in which most models underestimate the mean of most smoke quantities (black carbon, extinction, carbon monoxide) on the diagonal corridor between 16^∘ S, 6^∘ E, and 10^∘ S, 0^∘ E, in the 3–6 km layer, and overestimate them further south, closer to the coast, where less aerosol is present. Model representations of the above-cloud aerosol optical depth differ more widely. Most models overestimate the organic aerosol mass concentrations relative to those of black carbon, and with less skill, indicating model uncertainties in secondary organic aerosol processes. Regional-mean free-tropospheric model ambient single scattering albedos vary widely, between 0.83 and 0.93 compared with in situ dry measurements centered at 0.86, despite minimal impact of humidification on particulate scattering. The modeled ratios of the particulate extinction to the sum of the black carbon and organic aerosol mass concentrations (a mass extinction efficiency proxy) are typically too low and vary too little spatially, with significant inter-model differences. Most models overestimate the carbonaceous mass within the offshore boundary layer. Overall, the diversity in the model biases suggests that different model processes are responsible. The wide range of model optical properties requires further scrutiny because of their importance for radiative effect estimates.

The southeast Atlantic (SEA) region is host to a climatologically significant biomass burning aerosol layer overlying marine stratocumulus. We present the first results of the directly measured above-cloud aerosol optical depth (ACAOD) from the recent ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) airborne field campaign during August and September 2016. In our analysis, we use data from the Spectrometers for Sky-Scanning Sun-Tracking Atmospheric Research (4STAR) instrument and found an average ACAOD of 0.32 at 501 nm (range of 0.02 to 1.04), with an average Ångström exponent (AE) above clouds of 1.71. The AE is much lower at 1.25 for the full column (including below-cloud-level aerosol, with an average of 0.36 at 501 nm and a range of 0.02 to 0.74), indicating the presence of large aerosol particles, likely marine aerosol, in the lower atmospheric column. The ACAOD is observed from 4STAR to be highest near the coast at about 12^∘ S, whereas its variability is largest at the southern edge of the average aerosol plume, as indicated by 12 years of MODIS observations. In comparison to MODIS-derived ACAOD and long-term fine-mode plume-average AOD along a diagonal routine track extending out from the coast of Namibia, the directly measured ACAOD from 4STAR is slightly lower than the ACAOD product from MODIS. The peak ACAOD expected from MODIS AOD retrievals averaged over a long term along the routine diagonal flight track (peak of 0.5) was measured to be closer to coast in 2016 at about 1.5–4^∘ E, with 4STAR ACAOD averages showing a peak of 0.42. When considering the full observation set over the SEA, by spatially binning each sampled AOD, we obtain a geographically representative mean ACAOD of 0.37. Vertical profiles of AOD showcase the variability in the altitude of the aerosol plume and its separation from the cloud top. We measured larger AOD at a high altitude near the coast than farther from the coast, while generally observing a larger vertical gap farther from the coast. Changes in AOD with altitude are correlated with carbon monoxide, a gas tracer of the biomass burning aerosol plume. Vertical extent of gaps between aerosol and cloud show a wide distribution, with a near-zero gap being most frequent. The gap distribution with longitude is observed to be largest at about 7^∘ E, farther from coast than expected from previous studies.

In southern Africa, widespread agricultural fires produce substantial biomass burning (BB) emissions over the region. The seasonal smoke plumes associated with these emissions are then advected westward over the persistent stratocumulus cloud deck in the southeast Atlantic (SEA) Ocean, resulting in aerosol effects which vary with time and location. Much work has focused on the effects of these aerosol plumes, but previous studies have also described an elevated free tropospheric water vapor signal over the SEA. Water vapor influences climate in its own right, and it is especially important to consider atmospheric water vapor when quantifying aerosol–cloud interactions and aerosol radiative effects. Here we present airborne observations made during the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) campaign over the SEA Ocean. In observations collected from multiple independent instruments on the NASA P-3 aircraft (from near-surface to 6–7 km), we observe a strongly linear correlation between pollution indicators (carbon monoxide (CO) and aerosol loading) and atmospheric water vapor content, seen at all altitudes above the boundary layer. The focus of the current study is on the especially strong correlation observed during the ORACLES-2016 deployment (out of Walvis Bay, Namibia), but a similar relationship is also observed in the August 2017 and October 2018 ORACLES deployments. Using reanalyses from the European Centre for Medium-Range Weather Forecasts (ECMWF) and Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), and specialized WRF-Chem simulations, we trace the plume–vapor relationship to an initial humid, smoky continental source region, where it mixes with clean, dry upper tropospheric air and then is subjected to conditions of strong westward advection, namely the southern African easterly jet (AEJ-S). Our analysis indicates that air masses likely left the continent with the same relationship between water vapor and carbon monoxide as was observed by aircraft. This linear relationship developed over the continent due to daytime convection within a deep continental boundary layer (up to ∼5–6 km) and mixing with higher-altitude air, which resulted in fairly consistent vertical gradients in CO and water vapor, decreasing with altitude and varying in time, but this water vapor does not originate as a product of the BB combustion itself. Due to a combination of conditions and mixing between the smoky, moist continental boundary layer and the dry and fairly clean upper-troposphere air above (∼6 km), the smoky, humid air is transported by strong zonal winds and then advected over the SEA (to the ORACLES flight region) following largely isentropic trajectories. Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) back trajectories support this interpretation. This work thus gives insights into the conditions and processes which cause water vapor to covary with plume strength. Better understanding of this relationship, including how it varies spatially and temporally, is important to accurately quantify direct, semi-direct, and indirect aerosol effects over this region.

All-sky direct aerosol radiative effects (DARE) play a significant yet still uncertain role in climate. This is partly due to poorly quantified radiative properties of aerosol above clouds (AAC). We compute global estimates of shortwave top-of-atmosphere DARE over opaque water clouds (OWCs), DARE_OWC, using observation-based aerosol and cloud radiative properties from a combination of A-Train satellite sensors and a radiative transfer model. There are three major differences between our DARE_OWC calculations and previous studies: (1) we use the depolarization ratio method (DR) on CALIOP (Cloud–Aerosol Lidar with Orthogonal Polarization) Level 1 measurements to compute the AAC frequencies of occurrence and the AAC aerosol optical depths (AODs), thus introducing fewer uncertainties compared to using the CALIOP standard product; (2) we apply our calculations globally, instead of focusing exclusively on regional AAC hotspots such as the southeast Atlantic; and (3) instead of the traditional look-up table approach, we use a combination of satellite-based sensors to obtain AAC intensive radiative properties. Our results agree with previous findings on the dominant locations of AAC (south and northeast Pacific, tropical and southeast Atlantic, northern Indian Ocean and northwest Pacific), the season of maximum occurrence and aerosol optical depths (a majority in the 0.01–0.02 range and that can exceed 0.2 at 532 nm) across the globe. We find positive averages of global seasonal DARE_OWC between 0.13 and 0.26 W m⁻² (i.e., a warming effect on climate). Regional seasonal DARE_OWC values range from −0.06 W m⁻² in the Indian Ocean offshore from western Australia (in March–April–May) to 2.87 W m⁻² in the southeast Atlantic (in September–October–November). High positive values are usually paired with high aerosol optical depths ( > 0.1) and low single scattering albedos ( < 0.94), representative of, for example, biomass burning aerosols. Because we use different spatial domains, temporal periods, satellite sensors, detection methods and/or associated uncertainties, the DARE_OWC estimates in this study are not directly comparable to previous peer-reviewed results. Despite these differences, we emphasize that the DARE_OWC estimates derived in this study are generally higher than previously reported. The primary reasons for our higher estimates are (i) the possible underestimate of the number of dust-dominated AAC cases in our study; (ii) our use of Level 1 CALIOP products (instead of CALIOP Level 2 products in previous studies) for the detection and quantification of AAC aerosol optical depths, which leads to larger estimates of AOD above OWC; and (iii) our use of gridded 4° × 5° seasonal means of aerosol and cloud properties in our DARE_OWC calculations instead of simultaneously derived aerosol and cloud properties from a combination of A-Train satellite sensors. Each of these areas is explored in depth with detailed discussions that explain both the rationale for our specific approach and the subsequent ramifications for our DARE calculations.