Search Results

Now showing 1 - 10 of 10
Loading...
Thumbnail Image
Item

3+2 + X : what is the most useful depolarization input for retrieving microphysical properties of non-spherical particles from lidar measurements using the spheroid model of Dubovik et al. (2006)?

2019, Tesche, Matthias, Kolgotin, Alexei, Haarig, Moritz, Burton, Sharon P., Ferrare, Richard A., Hostetler, Chris A., Müller, Detlef

The typical multiwavelength aerosol lidar data set for inversion of optical to microphysical parameters is composed of three backscatter coefficients (β) at 355, 532, and 1064 nm and two extinction coefficients (α) at 355 and 532 nm. This data combination is referred to as a 3β C 2α or 3 + 2 data set. This set of data is sufficient for retrieving some important microphysical particle parameters if the particles have spherical shape. Here, we investigate the effect of including the particle linear depolarization ratio (δ) as a third input parameter for the inversion of lidar data. The inversion algorithm is generally not used if measurements show values of d that exceed 0.10 at 532 nm, i.e. in the presence of nonspherical particles such as desert dust, volcanic ash, and, under special circumstances, biomass-burning smoke. We use experimental data collected with instruments that are capable of measuring d at all three lidar wavelengths with an inversion routine that applies the spheroidal light-scattering model of Dubovik et al. (2006) with a fixed axis-ratio distribution to replicate scattering properties of non-spherical particles. The inversion gives the fraction of spheroids required to replicate the optical data as an additional output parameter. This is the first systematic test of the effect of using all theoretically possible combinations of d taken at 355, 532, and 1064 nm as input in the lidar data inversion. We find that depolarization information of at least one wavelength already provides useful information for the inversion of optical data that have been collected in the presence of non-spherical mineral dust particles. However, any choice of d will give lower values of the single-scattering albedo than the traditional 3 + 2 data set. We find that input data sets that include d355 give a spheroid fraction that closely resembles the dust ratio we obtain from using β532 and d532 in a methodology applied in aerosol-type separation. The use of d355 in data sets of two or three d? reduces the spheroid fraction that is retrieved when using d532 and d1064. Use of the latter two parameters without accounting for d355 generally leads to high spheroid fractions that we consider not trustworthy. The use of three d instead of two δ, including the constraint that one of these is measured at 355 nm does not provide any advantage over using 3 + 2 + d355 for the observations with varying contributions of mineral dust considered here. However, additional measurements at wavelengths different from 355 nm would be desirable for application to a wider range of aerosol scenarios that may include non-spherical smoke particles, which can have values of d355 that are indistinguishable from those found for mineral dust. We therefore conclude that - depending on measurement capability - the future standard input for inversion of lidar data taken in the presence of mineral dust particles and using the spheroid model of Dubovik et al. (2006) might be 3+2Cδ355 or 3 + 2 + δ355 + δ532. © 2019 The Author(s).

Loading...
Thumbnail Image
Item

Biomass burning and urban emission impacts in the Andes Cordillera region based on in situ measurements from the Chacaltaya observatory, Bolivia (5240a.s.l.)

2019, Chauvigné, Aurélien, Aliaga, Diego, Sellegri, Karine, Montoux, Nadège, Krejci, Radovan, Močnik, Griša, Moreno, Isabel, Müller, Thomas, Pandolfi, Marco, Velarde, Fernando, Weinhold, Kay, Ginot, Patrick, Wiedensohler, Alfred, Andrade, Marcos, Laj, Paolo

This study documents and analyses a 4-year continuous record of aerosol optical properties measured at the Global Atmosphere Watch (GAW) station of Chacaltaya (CHC; 5240a.s.l.), in Bolivia. Records of particle light scattering and particle light absorption coefficients are used to investigate how the high Andean Cordillera is affected by both long-range transport and by the fast-growing agglomeration of La Paz-El Alto, located approximately 20km away and 1.5km below the sampling site. The extended multi-year record allows us to study the properties of aerosol particles for different air mass types, during wet and dry seasons, also covering periods when the site was affected by biomass burning in the Bolivian lowlands and the Amazon Basin. The absorption, scattering, and extinction coefficients (median annual values of 0.74, 12.14, and 12.96Mm-1 respectively) show a clear seasonal variation with low values during the wet season (0.57, 7.94, and 8.68Mm-1 respectively) and higher values during the dry season (0.80, 11.23, and 14.51Mm-1 respectively). The record is driven by variability at both seasonal and diurnal scales. At a diurnal scale, all records of intensive and extensive aerosol properties show a pronounced variation (daytime maximum, night-time minimum), as a result of the dynamic and convective effects. The particle light absorption, scattering, and extinction coefficients are on average 1.94, 1.49, and 1.55 times higher respectively in the turbulent thermally driven conditions than the more stable conditions, due to more efficient transport from the boundary layer. Retrieved intensive optical properties are significantly different from one season to the other, reflecting the changing aerosol emission sources of aerosol at a larger scale. Using the wavelength dependence of aerosol particle optical properties, we discriminated between contributions from natural (mainly mineral dust) and anthropogenic (mainly biomass burning and urban transport or industries) emissions according to seasons and local circulation. The main sources influencing measurements at CHC are from the urban area of La Paz-El Alto in the Altiplano and from regional biomass burning in the Amazon Basin. Results show a 28% to 80% increase in the extinction coefficients during the biomass burning season with respect to the dry season, which is observed in both tropospheric dynamic conditions. From this analysis, long-term observations at CHC provide the first direct evidence of the impact of biomass burning emissions of the Amazon Basin and urban emissions from the La Paz area on atmospheric optical properties at a remote site all the way to the free troposphere. © Author(s) 2019.

Loading...
Thumbnail Image
Item

The effect of rapid relative humidity changes on fast filter-based aerosol-particle light-absorption measurements: Uncertainties and correction schemes

2019, Düsing, Sebastian, Wehner, Birgit, Müller, Thomas, Stöcker, Almond, Wiedensohler, Alfred

Measuring vertical profiles of the particle light-absorption coefficient by using absorption photometers may face the challenge of fast changes in relative humidity (RH). These absorption photometers determine the particle light-absorption coefficient due to a change in light attenuation through a particle-loaded filter. The filter material, however, takes up or releases water with changing relative humidity (RH in %), thus influencing the light attenuation. A sophisticated set of laboratory experiments was therefore conducted to investigate the effect of fast RH changes (dRH/dt) on the particle light-absorption coefficient (σabs in Mm-1) derived with two absorption photometers. The RH dependence was examined based on different filter types and filter loadings with respect to loading material and areal loading density. The Single Channel Tricolor Absorption Photometer (STAP) relies on quartz-fiber filter, and the microAeth® MA200 is based on a polytetrafluoroethylene (PTFE) filter band. Furthermore, three cases were investigated: clean filters, filters loaded with black carbon (BC), and filters loaded with ammonium sulfate. The filter areal loading densities (ρ∗) ranged from 3.1 to 99.6 mg m-2 in the case of the STAP and ammonium sulfate and 1.2 to 37.6 mg m-2 in the case the MA200. Investigating BC-loaded cases, M8 scroll mrow miBCm 15pt was in the range of 2.9 to 43.0 and 1.1 to 16.3 mg m-2 for the STAP and MA200, respectively.

Both instruments revealed opposing responses to relative humidity changes ("RH) with different magnitudes. The STAP shows a linear dependence on relative humidity changes. The MA200 is characterized by a distinct exponential recovery after its filter was exposed to relative humidity changes. At a wavelength of 624 nm and for the default 60 s running average output, the STAP reveals an absolute change in σabs per absolute change of RH ("σabsĝ•"RH) of 0.14 Mm-1 %-1 in the clean case, 0.29 Mm-1 %-1 in the case of BC-loaded filters, and 0.21 Mm-1 %-1 in the case filters loaded with ammonium sulfate. The 60 s running average of the particle light-absorption coefficient at 625 nm measured with the MA200 revealed a response of around -0.4 Mm-1 %-1 for all three cases. Whereas the response of the STAP varies over the different loading materials, in contrast, the MA200 was quite stable. The response was, for the STAP, in the range of 0.17 to 0.24 Mm-1 %-1 and, in the case of ammonium sulfate loading and in the BC-loaded case, 0.17 to 0.62 Mm-1 %-1. In the ammonium sulfate case, the minimum response shown by the MA200 was -0.42 with a maximum of -0.36 Mm-1 %-1 and a minimum of -0.42 and maximum -0.37 Mm-1 %-1 in the case of BC.

A linear correction function for the STAP was developed here. It is provided by correlating 1 Hz resolved recalculated particle light-absorption coefficients and RH change rates. The linear response is estimated at 10.08 Mm-1 s-1 %-1. A correction approach for the MA200 is also provided; however, the behavior of the MA200 is more complex. Further research and multi-instrument measurements have to be conducted to fully understand the underlying processes, since the correction approach resulted in different correction parameters across various experiments. However, the exponential recovery after the filter of the MA200 experienced a RH change could be reproduced. However, the given correction approach has to be estimated with other RH sensors as well, since each sensor has a different response time. And, for the given correction approaches, the uncertainties could not be estimated, which was mainly due to the response time of the RH sensor. Therefore, we do not recommend using the given approaches. But they point in the right direction, and despite the imperfections, they are useful for at least estimating the measurement uncertainties due to relative humidity changes.

Due to our findings, we recommend using an aerosol dryer upstream of absorption photometers to reduce the RH effect significantly. Furthermore, when absorption photometers are used in vertical measurements, the ascending or descending speed through layers of large relative humidity gradients has to be low to minimize the observed RH effect. But this is simply not possible in some scenarios, especially in unmixed layers or clouds. Additionally, recording the RH of the sample stream allows correcting for the bias during post-processing of the data. This data correction leads to reasonable results, according to the given example in this study. © Author(s) 2019.

Loading...
Thumbnail Image
Item

Profiles of cloud condensation nuclei, dust mass concentration, and ice-nucleating-particle-relevant aerosol properties in the Saharan Air Layer over Barbados from polarization lidar and airborne in situ measurements

2019, Haarig, Moritz, Walser, Adrian, Ansmann, Albert, Dollner, Maximilian, Althausen, Dietrich, Sauer, Daniel, Farrell, David, Weinzierl, Bernadett

The present study aims to evaluate lidar retrievals of cloud-relevant aerosol properties by using polarization lidar and coincident airborne in situ measurements in the Saharan Air Layer (SAL) over the Barbados region. Vertical profiles of the number concentration of cloud condensation nuclei (CCN), large particles (diameter d > 500 nm), surface area, mass, and ice-nucleating particle (INP) concentration are derived from the lidar measurements and compared with CCN concentrations and the INP-relevant aerosol properties measured in situ with aircraft. The measurements were performed in the framework of the Saharan Aerosol Longrange Transport and Aerosol-Cloud-Interaction Experiment (SALTRACE) in summer 2013. The CCN number concentrations derived from lidar observations were up to a factor of 2 higher than the ones measured in situ aboard the research aircraft Falcon. Possible reasons for the difference are discussed. The number concentration of particles with a dry radius of more than 250 nm and the surface-area concentration obtained from the lidar observations and used as input for the INP parameterizations agreed well (< 30 %-50 % deviation) with the aircraft measurements. In a pronounced lofted dust layer during summer (10 July 2013), the lidar retrieval yielded 100-300 CCN per cubic centimeter at 0.2 % water supersaturation and 10-200 INPs per liter at-25?C. Excellent agreement was also obtained in the comparison of mass concentration profiles. During the SALTRACE winter campaign (March 2014), the dust layer from Africa was mixed with smoke particles which dominated the CCN number concentration. This example highlights the unique lidar potential to separate smoke and dust contributions to the CCN reservoir and thus to identify the sensitive role of smoke in trade wind cumuli developments over the tropical Atlantic during the winter season. © 2017 Georg Thieme Verlag. All rights reserved.

Loading...
Thumbnail Image
Item

Retrieval of ice-nucleating particle concentrations from lidar observations and comparison with UAV in situ measurements

2019, Marinou, Eleni, Tesche, Matthias, Nenes, Athanasios, Ansmann, Albert, Schrod, Jann, Mamali, Dimitra, Tsekeri, Alexandra, Pikridas, Michael, Baars, Holger, Engelmann, Ronny, Voudouri, Kalliopi-Artemis, Solomos, Stavros, Sciare, Jean, Groß, Silke, Ewald, Florian, Amiridis, Vassilis

Aerosols that are efficient ice-nucleating particles (INPs) are crucial for the formation of cloud ice via heterogeneous nucleation in the atmosphere. The distribution of INPs on a large spatial scale and as a function of height determines their impact on clouds and climate. However, in situ measurements of INPs provide sparse coverage over space and time. A promising approach to address this gap is to retrieve INP concentration profiles by combining particle concentration profiles derived by lidar measurements with INP efficiency parameterizations for different freezing mechanisms (immersion freezing, deposition nucleation). Here, we assess the feasibility of this new method for both ground-based and spaceborne lidar measurements, using in situ observations collected with unmanned aerial vehicles (UAVs) and subsequently analyzed with the FRIDGE (FRankfurt Ice nucleation Deposition freezinG Experiment) INP counter from an experimental campaign at Cyprus in April 2016. Analyzing five case studies we calculated the cloud-relevant particle number concentrations using lidar measurements (n250,dry with an uncertainty of 20 % to 40 % and Sdry with an uncertainty of 30 % to 50 %), and we assessed the suitability of the different INP parameterizations with respect to the temperature range and the type of particles considered. Specifically, our analysis suggests that our calculations using the parameterization of Ullrich et al. (2017) (applicable for the temperature range −50 to −33 ∘C) agree within 1 order of magnitude with the in situ observations of nINP; thus, the parameterization of Ullrich et al. (2017) can efficiently address the deposition nucleation pathway in dust-dominated environments. Additionally, our calculations using the combination of the parameterizations of DeMott et al. (2015, 2010) (applicable for the temperature range −35 to −9 ∘C) agree within 2 orders of magnitude with the in situ observations of INP concentrations (nINP) and can thus efficiently address the immersion/condensation pathway of dust and nondust particles. The same conclusion is derived from the compilation of the parameterizations of DeMott et al. (2015) for dust and Ullrich et al. (2017) for soot.

Loading...
Thumbnail Image
Item

Sun photometer retrievals of Saharan dust properties over Barbados during SALTRACE

2019, Toledano, Carlos, Torres, Benjamín, Velasco-Merino, Cristian, Althausen, Dietrich, Groß, Silke, Wiegner, Matthias, Weinzierl, Bernadett, Gasteiger, Josef, Ansmann, Albert, González, Ramiro, Mateos, David, Farrel, David, Müller, Thomas, Haarig, Moritz, Cachorro, Victoria E.

The Saharan Aerosol Long-Range Transport and Aerosol-Cloud-Interaction Experiment (SALTRACE) was devoted to the investigation of Saharan dust properties over the Caribbean. The campaign took place in June-July 2013. A wide set of ground-based and airborne aerosol instrumentation was deployed at the island of Barbados for a comprehensive experiment. Several sun photometers performed measurements during this campaign: two AERONET (Aerosol Robotic Network) Cimel sun photometers and the Sun and Sky Automatic Radiometer (SSARA). The sun photometers were co-located with the ground-based multi-wavelength lidars BERTHA (Backscatter Extinction lidar Ratio Temperature Humidity profiling Apparatus) and POLIS (Portable Lidar System). Aerosol properties derived from direct sun and sky radiance observations are analyzed, and a comparison with the co-located lidar and in situ data is provided. The time series of aerosol optical depth (AOD) allows identifying successive dust events with short periods in between in which the marine background conditions were observed. The moderate aerosol optical depth in the range of 0.3 to 0.6 was found during the dust periods. The sun photometer infrared channel at the 1640nm wavelength was used in the retrieval to investigate possible improvements to aerosol size retrievals, and it was expected to have a larger sensitivity to coarse particles. The comparison between column (aerosol optical depth) and surface (dust concentration) data demonstrates the connection between the Saharan Air Layer and the boundary layer in the Caribbean region, as is shown by the synchronized detection of the successive dust events in both datasets. However the differences of size distributions derived from sun photometer data and in situ observations reveal the difficulties in carrying out a column closure study. © 2019 All rights reserved.

Loading...
Thumbnail Image
Item

Vertical aerosol distribution in the southern hemispheric midlatitudes as observed with lidar in Punta Arenas, Chile (53.2° and 70.9° W), during ALPACA

2019, Foth, Andreas, Kanitz, Thomas, Engelmann, Ronny, Baars, Holger, Radenz, Martin, Seifert, Patric, Barja, Boris, Fromm, Michael, Kalesse, Heike, Ansmann, Albert

Within this publication, lidar observations of the vertical aerosol distribution above Punta Arenas, Chile (53.2 S and 70.9 W), which have been performed with the Raman lidar PollyXT from December 2009 to April 2010, are presented. Pristine marine aerosol conditions related to the prevailing westerly circulation dominated the measurements. Lofted aerosol layers could only be observed eight times during the whole measurement period. Two case studies are presented showing long-range transport of smoke from biomass burning in Australia and regionally transported dust from the Patagonian Desert, respectively. The aerosol sources are identified by trajectory analyses with the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) and FLEXible PARTicle dispersion model (FLEXPART). However, seven of the eight analysed cases with lofted layers show an aerosol optical thickness of less than 0.05. From the lidar observations, a mean planetary boundary layer (PBL) top height of 1150 350m was determined. An analysis of particle backscatter coefficients confirms that the majority of the aerosol is attributed to the PBL, while the free troposphere is characterized by a very low background aerosol concentration. The ground-based lidar observations at 532 and 1064 nm are supplemented by the Aerosol Robotic Network (AERONET) Sun photometers and the space-borne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) aboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO). The averaged aerosol optical thickness (AOT) determined by CALIOP was 0:02 0:01 in Punta Arenas from 2009 to 2010. © Author(s) 2019. This work is distributed under the Creative Commons Attribution 4.0 License.

Loading...
Thumbnail Image
Item

EARLINET evaluation of the CATS Level 2 aerosol backscatter coefficient product

2019, Proestakis, Emmanouil, Amiridis, Vassilis, Marinou, Eleni, Binietoglou, Ioannis, Ansmann, Albert, Wandinger, Ulla, Hofer, Julian, Yorks, John, Nowottnick, Edward, Makhmudov, Abduvosit, Papayannis, Alexandros, Pietruczuk, Aleksander, Gialitaki, Anna, Apituley, Arnoud, Szkop, Artur, Muñoz Porcar, Constantino, Bortoli, Daniele, Dionisi, Davide, Althausen, Dietrich, Mamali, Dimitra, Balis, Dimitris, Nicolae, Doina, Tetoni, Eleni, Liberti, Gian Luigi, Baars, Holger, Mattis, Ina, Stachlewska, Iwona Sylwia, Voudouri, Kalliopi Artemis, Mona, Lucia, Mylonaki, Maria, Perrone, Maria Rita, Costa, Maria João, Sicard, Michael, Papagiannopoulos, Nikolaos, Siomos, Nikolaos, Burlizzi, Pasquale, Pauly, Rebecca, Engelmann, Ronny, Abdullaev, Sabur, Pappalardo, Gelsomina

We present the evaluation activity of the European Aerosol Research Lidar Network (EARLINET) for the quantitative assessment of the Level 2 aerosol backscatter coefficient product derived by the Cloud-Aerosol Transport System (CATS) aboard the International Space Station (ISS; Rodier et al., 2015). The study employs correlative CATS and EARLINET backscatter measurements within a 50km distance between the ground station and the ISS overpass and as close in time as possible, typically with the starting time or stopping time of the EARLINET performed measurement time window within 90min of the ISS overpass, for the period from February 2015 to September 2016. The results demonstrate the good agreement of the CATS Level 2 backscatter coefficient and EARLINET. Three ISS overpasses close to the EARLINET stations of Leipzig, Germany; Évora, Portugal; and Dushanbe, Tajikistan, are analyzed here to demonstrate the performance of the CATS lidar system under different conditions. The results show that under cloud-free, relative homogeneous aerosol conditions, CATS is in good agreement with EARLINET, independent of daytime and nighttime conditions. CATS low negative biases are observed, partially attributed to the deficiency of lidar systems to detect tenuous aerosol layers of backscatter signal below the minimum detection thresholds; these are biases which may lead to systematic deviations and slight underestimations of the total aerosol optical depth (AOD) in climate studies. In addition, CATS misclassification of aerosol layers as clouds, and vice versa, in cases of coexistent and/or adjacent aerosol and cloud features, occasionally leads to non-representative, unrealistic, and cloud-contaminated aerosol profiles. Regarding solar illumination conditions, low negative biases in CATS backscatter coefficient profiles, of the order of 6.1%, indicate the good nighttime performance of CATS. During daytime, a reduced signal-to-noise ratio by solar background illumination prevents retrievals of weakly scattering atmospheric layers that would otherwise be detectable during nighttime, leading to higher negative biases, of the order of 22.3%. © Author(s) 2019.

Loading...
Thumbnail Image
Item

EURODELTA III exercise: An evaluation of air quality models’ capacity to reproduce the carbonaceous aerosol

2019, Mircea, Mihaela, Bessagnet, Bertrand, D'Isidoro, Massimo, Pirovano, Guido, Aksoyoglu, Sebnem, Ciarelli, Giancarlo, Tsyro, Svetlana, Manders, Astrid, Bieser, Johannes, Stern, Rainer, Vivanco, Marta García, Cuvelier, Cornelius, Aas, Wenche, Prévôt, André S.H., Aulinger, Armin, Briganti, Gino, Calori, Giuseppe, Cappelletti, Andrea, Colette, Augustin, Couvidat, Florian, Fagerli, Hilde, Finardi, Sandro, Kranenburg, Richard, Rouïl, Laurence, Silibello, Camillo, Spindler, Gerald, Poulain, Laurent, Herrmann, Hartmut, Jimenez, Jose L., Day, Douglas A., Tiitta, Petri, Carbone, Samara

The carbonaceous aerosol accounts for an important part of total aerosol mass, affects human health and climate through its effects on physical and chemical properties of the aerosol, yet the understanding of its atmospheric sources and sinks is still incomplete. This study shows the state-of-the-art in modelling carbonaceous aerosol over Europe by comparing simulations performed with seven chemical transport models (CTMs) currently in air quality assessments in Europe: CAMx, CHIMERE, CMAQ, EMEP/MSC-W, LOTOS-EUROS, MINNI and RCGC. The simulations were carried out in the framework of the EURODELTA III modelling exercise and were evaluated against field measurements from intensive campaigns of European Monitoring and Evaluation Programme (EMEP) and the European Integrated Project on Aerosol Cloud Climate and Air Quality Interactions (EUCAARI). Model simulations were performed over the same domain, using as much as possible the same input data and covering four seasons: summer (1–30 June 2006), winter (8 January – 4 February 2007), autumn (17 September- 15 October 2008) and spring (25 February - 26 March 2009). The analyses of models’ performances in prediction of elemental carbon (EC) for the four seasons and organic aerosol components (OA) for the last two seasons show that all models generally underestimate the measured concentrations. The maximum underestimation of EC is about 60% and up to about 80% for total organic matter (TOM). The underestimation of TOM outside of highly polluted area is a consequence of an underestimation of secondary organic aerosol (SOA), in particular of its main contributor: biogenic secondary aerosol (BSOA). This result is independent on the SOA modelling approach used and season. The concentrations and daily cycles of total primary organic matter (TPOM) are generally better reproduced by the models since they used the same anthropogenic emissions. However, the combination of emissions and model formulation leads to overestimate TPOM concentrations in 2009 for most of the models. All models capture relatively well the SOA daily cycles at rural stations mainly due to the spatial resolution used in the simulations. For the investigated carbonaceous aerosol compounds, the differences between the concentrations simulated by different models are lower than the differences between the concentrations simulated with a model for different seasons. © 2019 The Authors

Loading...
Thumbnail Image
Item

The EMEP Intensive Measurement Period campaign, 2008-2009: Characterizing carbonaceous aerosol at nine rural sites in Europe

2019, Yttri, Karl Espen, Simpson, David, Bergström, Robert, Kiss, Gyula, Szidat, Sönke, Ceburnis, Darius, Eckhardt, Sabine, Hueglin, Christoph, Nøjgaard, Jacob Klenø, Perrino, Cinzia, Pisso, Ignazio, Prevot, Andre Stephan Henry, Putaud, Jean-Philippe, Spindler, Gerald, Vana, Milan, Zhang, Yan-Lin, Aas, Wenche

Carbonaceous aerosol (total carbon, TCp) was source apportioned at nine European rural background sites, as part of the European Measurement and Evaluation Programme (EMEP) Intensive Measurement Periods in fall 2008 and winter/spring 2009. Five predefined fractions were apportioned based on ambient measurements: elemental and organic carbon, from combustion of biomass (ECbb and OCbb) and from fossil-fuel (ECff and OCff) sources, and remaining non-fossil organic carbon (OCrnf), dominated by natural sources.OCrnf made a larger contribution to TCp than anthropogenic sources (ECbb, OCbb, ECff, and OCff) at four out of nine sites in fall, reflecting the vegetative season, whereas anthropogenic sources dominated at all but one site in winter/spring. Biomass burning (OCbb + ECbb) was the major anthropogenic source at the central European sites in fall, whereas fossil-fuel (OCff + ECff) sources dominated at the southernmost and the two northernmost sites. Residential wood burning emissions explained 30 %-50 % of TCp at most sites in the first week of sampling in fall, showing that this source can be the dominant one, even outside the heating season. In winter/spring, biomass burning was the major anthropogenic source at all but two sites, reflecting increased residential wood burning emissions in the heating season. Fossil-fuel sources dominated EC at all sites in fall, whereas there was a shift towards biomass burning for the southernmost sites in winter/spring.Model calculations based on base-case emissions (mainly officially reported national emissions) strongly underpredicted observational derived levels of OCbb and ECbb outside Scandinavia. Emissions based on a consistent bottom-up inventory for residential wood burning (and including intermediate volatility compounds, IVOCs) improved model results compared to the base-case emissions, but modeled levels were still substantially underestimated compared to observational derived OCbb and ECbb levels at the southernmost sites.Our study shows that natural sources are a major contributor to carbonaceous aerosol in Europe, even in fall and in winter/spring, and that residential wood burning emissions are equally as large as or larger than that of fossil-fuel sources, depending on season and region. The poorly constrained residential wood burning emissions for large parts of Europe show the obvious need to improve emission inventories, with harmonization of emission factors between countries likely being the most important step to improve model calculations for biomass burning emissions, and European PM2.5 concentrations in general. © Author(s) 2019.