The Lagrangian particle dispersion model FLEXPART was originally designed for calculating long-range and mesoscale dispersion of air pollutants from point sources, such that occurring after an accident in a nuclear power plant. In the meantime, FLEXPART has evolved into a comprehensive tool for atmospheric transport modeling and analysis at different scales. A need for further multiscale modeling and analysis has encouraged new developments in FLEXPART. In this paper, we present a FLEXPART version that works with the Weather Research and Forecasting (WRF) mesoscale meteorological model. We explain how to run this new model and present special options and features that differ from those of the preceding versions. For instance, a novel turbulence scheme for the convective boundary layer has been included that considers both the skewness of turbulence in the vertical velocity as well as the vertical gradient in the air density. To our knowledge, FLEXPART is the first model for which such a scheme has been developed. On a more technical level, FLEXPART-WRF now offers effective parallelization, and details on computational performance are presented here. FLEXPART-WRF output can either be in binary or Network Common Data Form (NetCDF) format, both of which have efficient data compression. In addition, test case data and the source code are provided to the reader as a Supplement. This material and future developments will be accessible at http://www.flexpart.eu.Peer ReviewedPostprint (published version

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 ...

No abstract

The structure of turbulent flows along a transition between tall-forested canopies and forest clearings continues to be an active research topic in canopy turbulence. The difficulties in describing the turbulent flow along these transitions stem from the fact that the vertical structure of the canopy and its leaf area distribution cannot be ignored or represented by an effective roughness length. Large-eddy simulation (LES) runs were performed to explore the effect of a homogeneous variation in the forest leaf area index (LAI) on the turbulent flow across forest edges. A nested grid numerical method was used to ensure the development of a deep boundary layer above the forest while maintaining a sufficiently high resolution in the region close to the ground. It was demonstrated that the LES here predicted first-order and second-order mean velocity statistics within the canopy that agree with reported Reynolds-Averaged Navier-Stokes (RANS) model results, field and laboratory experiments. In the simulations reported here, the LAI was varied between 2 and 8 spanning a broad range of observed LAI in terrestrial ecosystems. By increasing the forest LAI, the mean flow properties both within the forest and in the clearing near the forest edge were altered in two fundamental ways: near the forest edge and into the clearing, the flow statistical properties resembled the so-called back-facing step (BFS) flow with a mean recirculation zone near the edge. Another recirculation zone sets up downstream of the clearing as the flow enters the tall forest canopy. The genesis of this within-forest recirculation zone can be primarily described using the interplay between the mean pressure gradients (forcing the flow) and the drag force (opposing the flow). Using the LES results, a simplified analytical model was also proposed to explain the location of the recirculation zone inside the canopy and its dependence on the forest LAI. Furthermore, a simplified scaling argument that decomposes the mean velocity at the outflow edge into a superposition of 'exit flow' and BFS-like flow with their relative importance determined by LAI was explored.

Abstract. In this study we have analyzed whether tourist cruise ships have an influence on measured sulfur dioxide (SO2), ozone (O3), Aitken mode particle and equivalent black carbon (EBC) concentrations at Ny Ålesund and Zeppelin Mountain on Svalbard in the Norwegian Arctic during summer. We separated the measurement data set into periods when ships were present and periods when ships were not present in the Kongsfjord area, according to a long-term record of the number of passengers visiting Ny Ålesund. We show that when ships with more than 50 passengers cruise in the Kongsfjord, measured daytime mean concentrations of 60 nm particles and EBC in summer show enhancements of 72 and 45%, respectively, relative to values when ships are not present. Even larger enhancements of 81 and 72% were found for stagnant conditions. In contrast, O3 concentrations were 5% lower on average and 7% lower under stagnant conditions, due to titration of O3 with the emitted nitric oxide (NO). The differences between the two data subsets are largest for the highest measured percentiles, while relatively small differences were found for the median concentrations, indicating that ship plumes are sampled relatively infrequently even when ships are present although they carry high pollutant concentrations. We estimate that the ships increased the total summer mean concentrations of SO2, 60 nm particles and EBC by 15, 18 and 11%, respectively. Our findings have two important implications. Firstly, even at such a remote Arctic observatory as Zeppelin, the measurements can be influenced by tourist ship emissions. Careful data screening is recommended before summertime Zeppelin data is used for data analysis or for comparison with global chemistry transport models. However, Zeppelin remains as one of the most valuable Arctic observatories, as most other Arctic observatories face even larger local pollution problems. Secondly, given landing statistics of tourist ships on Svalbard, it is suspected that large parts of the Svalbard archipelago are affected by cruise ship emissions. Thus, our results may be taken as a warning signal of future pan-Arctic conditions if Arctic shipping becomes more frequent and emission regulations are not strict enough.