In situ particle sampling relationships to surface and turbulent fluxes using large eddy simulations with Lagrangian particles

Abstract

Abstract. Source functions for mechanically driven coarse-mode sea spray and dust aerosol particles span orders of magnitude owing to a combination of physical sensitivity in the system and large measurement uncertainty. Outside special idealized settings (such as wind tunnels), aerosol particle fluxes are largely inferred from a host of methods, including local eddy correlation, gradient methods, and dry deposition methods. In all of these methods, it is difficult to relate point measurements from towers, ships, or aircraft to a general representative flux of aerosol particles. This difficulty is from the particles' inhomogeneous distribution due to multiple spatiotemporal scales of an evolving marine environment. We hypothesize that the current representation of a point in situ measurement of sea spray or dust particles is a likely contributor to the unrealistic range of flux and concentration outcomes in the literature. This paper aims to help the interpretation of field data: we conduct a series of high-resolution, cloud-free large eddy simulations (LESs) with Lagrangian particles to better understand the temporal evolution and volumetric variability of coarse- to giant-mode marine aerosol particles and their relationship to turbulent transport. The study begins by describing the Lagrangian LES model framework and simulates flux measurements that were made using numerical analogs to field practices such as the eddy covariance method. Using these methods, turbulent flux sampling is quantified based on key features such as coherent structures within the marine atmospheric boundary layer (MABL) and aerosol particle size. We show that for an unstable atmospheric stability, the MABL exhibits large coherent eddy structures, and as a consequence, the flux measurement outcome becomes strongly tied to spatial length scales and relative sampling of crosswise and streamwise sampling. For example, through the use of ogive curves, a given sampling duration of a fixed numerical sampling instrument is found to capture 80 % of the aerosol flux given a sampling rate of zf/w∗∼ 0.2, whereas a spanwise moving instrument results in a 95 % capture. These coherent structures and other canonical features contribute to the lack of convergence to the true aerosol vertical flux at any height. As expected, sampling all of the flow features results in a statistically robust flux signal. Analysis of a neutral boundary layer configuration results in a lower predictive range due to weak or no vertical roll structures compared to the unstable boundary layer setting. Finally, we take the results of each approach and compare their surface flux variability: a baseline metric used in regional and global aerosol models.