In recent years, marine, freshwater and terrestrial pollution with microplastics has been discussed extensively, whereas atmospheric microplastic transport has been largely overlooked. Here, we present global simulations of atmospheric transport of microplastic particles produced by road traffic (TWPstire wear particles and BWPsbrake wear particles), a major source that can be quantified relatively well. We find a high transport efficiencies of these particles to remote regions. About 34% of the emitted coarse TWPs and 30% of the emitted coarse BWPs (100 kt yr −1 and 40 kt yr −1 respectively) were deposited in the World Ocean. These amounts are of similar magnitude as the total estimated direct and riverine transport of TWPs and fibres to the ocean (64 kt yr −1). We suggest that the Arctic may be a particularly sensitive receptor region, where the light-absorbing properties of TWPs and BWPs may also cause accelerated warming and melting of the cryosphere.

The Lagrangian particle dispersion model FLEX-PART in its original version in the mid-1990s was designed for calculating the long-range and mesoscale dispersion of hazardous substances from point sources, such as those released after an accident in a nuclear power plant. Over the past decades, the model has evolved into a comprehensive tool for multi-scale atmospheric transport modeling and analysis and has attracted a global user community. Its application fields have been extended to a large range of atmospheric gases and aerosols, e.g., greenhouse gases, short-lived climate forcers like black carbon and volcanic ash, and it has also been used to study the atmospheric branch of the water cycle. Given suitable meteorological input data, it can be used for scales from dozens of meters to global. In particular, inverse modeling based on source-receptor relationships from FLEXPART has become widely used. In this paper, we present FLEXPART version 10.4, which works with meteorological input data from the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS) and data from the United States National Cen-ters of Environmental Prediction (NCEP) Global Forecast System (GFS). Since the last publication of a detailed FLEX-PART description (version 6.2), the model has been improved in different aspects such as performance, physicochemical parameterizations, input/output formats, and available preprocessing and post-processing software. The model code has also been parallelized using the Message Passing Interface (MPI). We demonstrate that the model scales well up to using 256 processors, with a parallel efficiency greater than 75 % for up to 64 processes on multiple nodes in runs with very large numbers of particles. The deviation from 100 % efficiency is almost entirely due to the remaining nonparallelized parts of the code, suggesting large potential for further speedup. A new turbulence scheme for the convective boundary layer has been developed that considers the skewness in the vertical velocity distribution (updrafts and downdrafts) and vertical gradients in air density. FLEXPART is the only model available considering both effects, making it highly accurate for small-scale applications, e.g., to quantify dispersion in the vicinity of a point source. The wet deposi-Published by Copernicus Publications on behalf of the European Geosciences Union. 4956 I. Pisso et al.: FLEXPART version 10.4tion scheme for aerosols has been completely rewritten and a new, more detailed gravitational settling parameterization for aerosols has also been implemented. FLEXPART has had the option of running backward in time from atmospheric concentrations at receptor locations for many years, but this has now been extended to also work for deposition values and may become useful, for instance, for the interpretation of ice core measurements. To our knowledge, to date FLEX-PART is the only model with that capability. Furthermore, the temporal variation and temperature dependence of chemical reactions with ...

Abstract. Sea-spray aerosols (SSA) are an important part of the climate system because of their effects on the global radiative budget -both directly as scatterers and absorbers of solar and terrestrial radiation, and indirectly as cloud condensation nuclei (CCN) influencing cloud formation, lifetime, and precipitation. In terms of their global mass, SSA have the largest uncertainty of all aerosols. In this study we review 21 SSA source functions from the literature, several of which are used in current climate models. In addition, we propose a new function. Even excluding outliers, the global annual SSA mass produced spans roughly 3-70 Pg yr −1 for the different source functions, for particles with dry diameter D p < 10 µm, with relatively little interannual variability for a given function. The FLEXPART Lagrangian particle dispersion model was run in backward mode for a large global set of observed SSA concentrations, comprised of several station networks and ship cruise measurement campaigns. FLEX-PART backward calculations produce gridded emission sensitivity fields, which can subsequently be multiplied with gridded SSA production fluxes in order to obtain modeled SSA concentrations. This allowed us to efficiently and simultaneously evaluate all 21 source functions against the measurements. Another advantage of this method is that sourceregion information on wind speed and sea surface temperatures (SSTs) could be stored and used for improving the SSA source function parameterizations. The best source functions reproduced as much as 70 % of the observed SSA concentration variability at several stations, which is comparable with "state of the art" aerosol models. The main driver of SSA production is wind, and we found that the best fit to the observation data could be obtained when the SSA production is proportional to U 3.5 10 , where U 10 is the source region averaged 10 m wind speed. A strong influence of SST on SSA production, with higher temperatures leading to higher production, could be detected as well, although the underlying physical mechanisms of the SST influence remains unclear. Our new source function with wind speed and temperature dependence gives a global SSA production for particles smaller than D p < 10 µm of 9 Pg yr −1 , and is the best fit to the observed concentrations.