Research Facilities
We use the latest atmospheric measurement techniques to study the properties of the aerosols particles in the most arduous areas of the globe, e.g. the Alps, the Himalaya, the Bolivian Andes. Measurements data are combined with Cluster Analysis and Chemical Transport Model simulations to obtain an holistic view of the underlying influence of air mass origin and chemical processes.

Particle size magnifier (PSM)

The Airmodus A10 PSM determines the size distribution of particles between 1 and 3 nm diameter. Particles grow inside the PSM with diethylene glycol as a working fluid to about 90 nm in diameter, after which counting is performed with an Airmodus A20 butanol condensation particle counter (CPC). The cut-off size of the PSM can be varied by altering the mixing ratio of the sample and diethylene glycol saturated flows.

Neutral cluster and Air Ion Spectrometer (NAIS)

The NAIS is an ion mobility spectrometer that measures the number size distribution of ions and total (charged and neutral) particles with electrical mobilities that range from 3.2 to 0.0013 cm2 V−1 s−1, or from 0.8 to 42 nm mobility diameters for ions and from 2 to 42 nm for neutral particles. Two differential mobility analysers classified the positive and negative ions.

Atmospheric Pressure interface Time of Flight Mass Spectrometer (APi-TOF), CI-APi-TOF and cluster chemical composition

The APi-TOF mass spectrometer measures the chemical composition of the naturally charged ions (TOFWERK AG and Aerodyne Research). In the APi section, the sample air passed through a critical orifice with a flow rate of 0.8 l min−1, after which the ions were guided by two quadrupoles and an ion lens to the time-of-flight mass spectrometer.


Atmospheric measurements are largely influenced by synoptic meteorology and local topography of the measurements site. Thus, accurate information on air mass origin and transport pathways to high-altitude sites is required. The ABC group have developed a new method, based on the source–receptor relationship (SRR) obtained from backwards WRF-FLEXPART simulations and a k-means clustering approach, to identify source regions of air masses arriving at measurement sites.

A key aspect of our method is that it is probabilistic, and for each observation time, more than one air mass (cluster) can influence the station, and the percentage influence of each air mass can be quantified. This is in contrast to binary methods, which label each observation time as influenced by either boundary layer or free-troposphere air masses.

The method was successfully applied in areas characterized by complex topography (i.e. the Chacaltaya Observatory) and for stations influenced by both local and long-range sources.

High-resolution meteorological modelling (step 1)

To generate a high-resolution, gridded dataset of meteorological variables that can be used to drive a Lagrangian dispersion model, we use the Advanced Research WRF model version 4.0.3 (Skamarock et al., 2019). One-way nesting is used: the outer domain provides boundary data for the inner nest, but the inner nest does not provide any feedback to the outer domain. To ensure that the long simulation remains close to reality throughout the 6-month period, the outer domain is nudged (i.e. analysis nudging) to the boundary conditions every 6 h.

Backward dispersion simulations (step 2)

The FLEXible PARTicle dispersion model (FLEXPART) is a Lagrangian transport and dispersion model which can be used for both forward and backward simulations. We use version FLEXPART-WRF_v3.3.2 to perform backward simulations and thus to determine the source regions of air masses arriving at a specific location. The FLEXPART simulations are driven using the meteorological output from the WRF simulation. Output from all the WRF output domains are used at a temporal resolution of 15 min. This high temporal resolution is a clear advantage over using reanalysis data, which at best are only available once per hour.


Emissions of both anthropogenic and biogenic chemical compounds are highly dependent on their geographical location, meteorological condition and human habits. Once released from the local source into the various layers of the atmosphere, those compounds undergoes hundred of thousands of chemical reactions - eventually leading to the formation of new particles - and are subject to removal process such as dry and wet deposition. These - non-linear - processes are highly depended on the physical and chemical characteristics of the aerosols population.

Additionally, aerosol particles can affect climate by interacting with the incoming solar radiation (direct effect) and by acting as cloud condensation nuclei (CCN), therefore changing the microphysics of cloud.    

We use the WRF-CHIMERE Chemical Transport model ( to simulate physical and chemical processes taking place into the atmosphere, from the injection of emissions in the planetary boundary layer (PBL), to chemical reaction of hundreds of chemical compounds to dry and wet deposition processes. The model routinely participates in intercomparison exercises and it is a member of the Copernicus Atmosphere Monitoring Service ( operational ensemble. It can be operated in a so call “online” configuration with the Weather Research and Forecast (WRFv3.71) model to include the exchange of the aerosol size distribution between CHIMERE and the meteorological model, i.e. WRF, which allows to investigate the direct and indirect effects of aerosol particles. It can be applied at various horizontal resolutions in a nested configuration, and it is therefore suitable for both global (hundreds of kilometers) and urban (1 km) scale applications. 

Meteorological data

We simulate meteorological input data using the Weather Research and Forecast (WRF) regional model. We nudge the simulations on the coarser grid towards reanalysis data from the National Centers for Environmental Prediction (NCEP) Climate Forecast System Version 2 ( with a temporal resolution of 6 h (i.e. surface grid nudging). For most of the application we use the Rapid Radiative Transfer Model (RRTMG) radiation scheme, the Thompson aerosol-aware MP scheme to treat the microphysics the Monin–Obukhov surface layer scheme, and the NOAA Land Surface Model scheme for land surface physics.

Initial and Boundary conditions

Initial and boundary conditions of aerosols and gas-phase constituents were retrieved from the climatological simulations of LMDz-INCA3 ( and the Goddard Chemistry Aerosol Radiation and Transport ( model. For aerosol constituent, in particular, the model includes inorganic species such as fine and coarse nitrate, ammonium, sulfate, dust, as well as organic carbon and black carbon.

Anthropogenic Emissions

Annual anthropogenic emissions of black carbon (BC), organic carbon (OC), methane (CH4), carbon dioxide (CO2), carbon monoxide (CO), ammonia (NH3), non-methane volatile organic compounds (NMVOCs), nitrogen oxides (NOx) and sulfur dioxide (SO2) are retrieved from CAMS ( for the year of interest at 0.1 x 0.1 degree resolution (corresponding to a resolution of about 10 km at the equator) and hourly distributed over the period of interest. NMVOCs are split into 221 chemical compounds and lumped for the chemical mechanism of interest (i.e. MELCHIOR or SAPRC). The total emission fluxes are additionally downscaled using a top-down approach (down to 1 km) with source-dependent proxy such as road maps, population density and landuse characteristic data. A vertical distribution of the emission is also used depending on the specific source category.

Biogenic Emissions

Biogenic emissions of NO, isoprene, limonene, α-pinene, β-pinene, ocimene, and humulene are prepared using the MEGAN ( model. Emission rates of 15 plant functional types (PFTs), at an original horizontal resolution of 0.008º × 0.008º, are re-gridded to match the resolution of both the coarse and high-resolution nested domains (down to 1 km). Standard emissions rate are adjusted based on several environmental factors, based on local radiation and temperature values (among others variables). 

Emission of seasalt, dust (including detailed mineralogy), NOx from lightning and resuspension material are included.


The chemical mechanism used for the gas-phase chemistry is the MELCHIOR2 scheme (alike the SAPRC scheme), including up to about 120 reactions with updated reaction rates. The ISORROPIA thermodynamic model is used to calculate the partitioning of the inorganic aerosol constituents and a logarithmic sectional distribution approach is deployed to treat the size distribution of aerosol particles using tens of adaptable bins (from 1nm up to 40µm). The model additionally account for coagulation process as well as binary nucleation of sulfuric acid (H2SO4) and water with the addition of HOMs formation and nucleation from biogenic compounds (i.e. α-pinene). The treatment of the organic aerosol chemistry in the model is treated within in the framework of the Volatility Basis Set (VBS) scheme.