Search Results

Now showing 1 - 10 of 63
  • Item
    Implementation of multirate time integration methods for air pollution modelling
    (München : European Geopyhsical Union, 2012) Schlegel, M.; Knoth, O.; Arnold, M.; Wolke, R.
    Explicit time integration methods are characterised by a small numerical effort per time step. In the application to multiscale problems in atmospheric modelling, this benefit is often more than compensated by stability problems and step size restrictions resulting from stiff chemical reaction terms and from a locally varying Courant-Friedrichs-Lewy (CFL) condition for the advection terms. Splitting methods may be applied to efficiently combine implicit and explicit methods (IMEX splitting). Complementarily multirate time integration schemes allow for a local adaptation of the time step size to the grid size. In combination, these approaches lead to schemes which are efficient in terms of evaluations of the right-hand side. Special challenges arise when these methods are to be implemented. For an efficient implementation, it is crucial to locate and exploit redundancies. Furthermore, the more complex programme flow may lead to computational overhead which, in the worst case, more than compensates the theoretical gain in efficiency. We present a general splitting approach which allows both for IMEX splittings and for local time step adaptation. The main focus is on an efficient implementation of this approach for parallel computation on computer clusters.
  • Item
    Meteorological and trace gas factors affecting the number concentration of atmospheric Aitken (DP Combining double low line 50 nm) particles in the continental boundary layer: Parameterization using a multivariate mixed effects model
    (München : European Geopyhsical Union, 2011) Mikkonen, S.; Korhonen, H.; Romakkaniemi, S.; Smith, J.N.; Joutsensaari, J.; Lehtinen, K.E.J.; Hamed, A.; Breider, T.J.; Birmili, W.; Spindler, G.; Plass-Duelmer, C.; Facchini, M.C.; Laaksonen, A.
    Measurements of aerosol size distribution and different gas and meteorological parameters, made in three polluted sites in Central and Southern Europe: Po Valley, Italy, Melpitz and Hohenpeissenberg in Germany, were analysed for this study to examine which of the meteorological and trace gas variables affect the number concentration of Aitken (Dp= 50 nm) particles. The aim of our study was to predict the number concentration of 50 nm particles by a combination of in-situ meteorological and gas phase parameters. The statistical model needs to describe, amongst others, the factors affecting the growth of newly formed aerosol particles (below 10 nm) to 50 nm size, but also sources of direct particle emissions in that size range. As the analysis method we used multivariate nonlinear mixed effects model. Hourly averages of gas and meteorological parameters measured at the stations were used as predictor variables; the best predictive model was attained with a combination of relative humidity, new particle formation event probability, temperature, condensation sink and concentrations of SO2, NO2 and ozone. The seasonal variation was also taken into account in the mixed model structure. Model simulations with the Global Model of Aerosol Processes (GLOMAP) indicate that the parameterization can be used as a part of a larger atmospheric model to predict the concentration of climatically active particles. As an additional benefit, the introduced model framework is, in theory, applicable for any kind of measured aerosol parameter.
  • Item
    The regional aerosol-climate model REMO-HAM
    (München : European Geopyhsical Union, 2012) Pietikäinen, J.-P.; O'Donnell, D.; Teichmann, C.; Karstens, U.; Pfeifer, S.; Kazil, J.; Podzun, R.; Fiedler, S.; Kokkola, H.; Birmili, W.; O'Dowd, C.; Baltensperger, U.; Weingartner, E.; Gehrig, R.; Spindler, G.; Kulmala, M.; Feichter, J.; Jacob, D.; Laaksonen, A.
    REMO-HAM is a new regional aerosol-climate model. It is based on the REMO regional climate model and includes most of the major aerosol processes. The structure for aerosol is similar to the global aerosol-climate model ECHAM5-HAM, for example the aerosol module HAM is coupled with a two-moment stratiform cloud scheme. On the other hand, REMO-HAM does not include an online coupled aerosol-radiation nor a secondary organic aerosol module. In this work, we evaluate the model and compare the results against ECHAM5-HAM and measurements. Four different measurement sites were chosen for the comparison of total number concentrations, size distributions and gas phase sulfur dioxide concentrations: Hyytiälä in Finland, Melpitz in Germany, Mace Head in Ireland and Jungfraujoch in Switzerland. REMO-HAM is run with two different resolutions: 50 × 50 km2 and 10 × 10 km2. Based on our simulations, REMO-HAM is in reasonable agreement with the measured values. The differences in the total number concentrations between REMO-HAM and ECHAM5-HAM can be mainly explained by the difference in the nucleation mode. Since we did not use activation nor kinetic nucleation for the boundary layer, the total number concentrations are somewhat underestimated. From the meteorological point of view, REMO-HAM represents the precipitation fields and 2 m temperature profile very well compared to measurement. Overall, we show that REMO-HAM is a functional aerosol-climate model, which will be used in further studies.
  • Item
    Number size distributions and seasonality of submicron particles in Europe 2008–2009
    (München : European Geopyhsical Union, 2011) Asmi, A.; Wiedensohler, A.; Laj, P.; Fjaeraa, A.-M.; Sellegri, K.; Birmili, W.; Weingartner, E.; Baltensperger, U.; Zdimal, V.; Zikova, N.; Putaud, J.-P.; Marinoni, A.; Tunved, P.; Hansson, H.-C.; Fiebig, M.; Kivekäs, N.; Lihavainen, H.; Asmi, E.; Ulevicius, V.; Aalto, P.P.; Swietlicki, E.; Kristensson, A.; Mihalopoulos, N.; Kalivitis, N.; Kalapov, I.; Kiss, G.; de Leeuw, G.; Henzing, B.; Harrison, R.M.; Beddows, D.; O'Dowd, C.; Jennings, S.G.; Flentje, H.; Weinhold, K.; Meinhardt, F.; Ries, L.; Kulmala, M.
    Two years of harmonized aerosol number size distribution data from 24 European field monitoring sites have been analysed. The results give a comprehensive overview of the European near surface aerosol particle number concentrations and number size distributions between 30 and 500 nm of dry particle diameter. Spatial and temporal distribution of aerosols in the particle sizes most important for climate applications are presented. We also analyse the annual, weekly and diurnal cycles of the aerosol number concentrations, provide log-normal fitting parameters for median number size distributions, and give guidance notes for data users. Emphasis is placed on the usability of results within the aerosol modelling community. We also show that the aerosol number concentrations of Aitken and accumulation mode particles (with 100 nm dry diameter as a cut-off between modes) are related, although there is significant variation in the ratios of the modal number concentrations. Different aerosol and station types are distinguished from this data and this methodology has potential for further categorization of stations aerosol number size distribution types. The European submicron aerosol was divided into characteristic types: Central European aerosol, characterized by single mode median size distributions, unimodal number concentration histograms and low variability in CCN-sized aerosol number concentrations; Nordic aerosol with low number concentrations, although showing pronounced seasonal variation of especially Aitken mode particles; Mountain sites (altitude over 1000 m a.s.l.) with a strong seasonal cycle in aerosol number concentrations, high variability, and very low median number concentrations. Southern and Western European regions had fewer stations, which decreases the regional coverage of these results. Aerosol number concentrations over the Britain and Ireland had very high variance and there are indications of mixed air masses from several source regions; the Mediterranean aerosol exhibit high seasonality, and a strong accumulation mode in the summer. The greatest concentrations were observed at the Ispra station in Northern Italy with high accumulation mode number concentrations in the winter. The aerosol number concentrations at the Arctic station Zeppelin in Ny-\AA lesund in Svalbard have also a strong seasonal cycle, with greater concentrations of accumulation mode particles in winter, and dominating summer Aitken mode indicating more recently formed particles. Observed particles did not show any statistically significant regional work-week or weekday related variation in number concentrations studied. Analysis products are made for open-access to the research community, available in a freely accessible internet site. The results give to the modelling community a reliable, easy-to-use and freely available comparison dataset of aerosol size distributions.
  • Item
    Evaluation on the role of sulfuric acid in the mechanisms of new particle formation for Beijing case
    (München : European Geopyhsical Union, 2011) Wang, Z.B.; Hu, M.; Yue, D.L.; Zheng, J.; Zhang, R.Y.; Wiedensohler, A.; Wu, Z.J.; Nieminen, T.; Boy, M.
    New particle formation (NPF) is considered as an important mechanism for gas-to-particle transformation, and gaseous sulfuric acid is believed as a crucial precursor. Up to now few field-based studies on nucleation mechanisms and the role of sulfuric acid were conducted in China. In this study, simultaneously measurements of particle number size distributions and gaseous sulfuric acid concentrations were performed from July to September in 2008. Totally, 22 new particle formation events were observed during the entire 85 campaign days. The results show that in the case of both higher source and sink values, the result of the competition between source and sink is more likely the key limiting factor to determine the observation of NPF events in Beijing. The concentrations of gaseous sulfuric acid show good correlations with freshly nucleated particles (N3-6 and formation rates (J3 and J1.5. The power-law relationship between H2SO4 concentration and N3-6 or J is adopted to explore the nucleation mechanism. The exponents are showed a great range (from 1 to 7). More than half of the NPF events exhibit an exponent larger than 2.5. For these cases, the thermodynamic process works better than the activation or kinetic nucleation theories to explain the nucleation events in urban atmosphere of Beijing.
  • Item
    First long-term study of particle number size distributions and new particle formation events of regional aerosol in the North China Plain
    (München : European Geopyhsical Union, 2011) Shen, X.J.; Sun, J.Y.; Zhang, Y.M.; Wehner, B.; Nowak, A.; Tuch, T.; Zhang, X.C.; Wang, T.T.; Zhou, H.G.; Zhang, X.L.; Dong, F.; Birmili, W.; Wiedensohler, A.
    Atmospheric particle number size distributions (size range 0.003–10 μm) were measured between March 2008 and August 2009 at Shangdianzi (SDZ), a rural research station in the North China Plain. These measurements were made in an attempt to better characterize the tropospheric background aerosol in Northern China. The mean particle number concentrations of the total particle, as well as the nucleation, Aitken, accumulation and coarse mode were determined to be 1.2 ± 0.9 × 104, 3.6 ± 7.9 × 103, 4.4 ± 3.4 × 103, 3.5 ± 2.8 × 103 and 2 ± 3 cm−3, respectively. A general finding was that the particle number concentration was higher during spring compared to the other seasons. The air mass origin had an important effect on the particle number concentration and new particle formation events. Air masses from northwest (i.e. inner Asia) favored the new particle formation events, while air masses from southeast showed the highest particle mass concentration. Significant diurnal variations in particle number were observed, which could be linked to new particle formation events, i.e. gas-to-particle conversion. During particle formation events, the number concentration of the nucleation mode rose up to maximum value of 104 cm−3. New particle formation events were observed on 36% of the effective measurement days. The formation rate ranged from 0.7 to 72.7 cm−3 s−1, with a mean value of 8.0 cm−3 s−1. The value of the nucleation mode growth rate was in the range of 0.3–14.5 nm h−1, with a mean value of 4.3 nm h−1. It was an essential observation that on many occasions the nucleation mode was able to grow into the size of cloud condensation nuclei (CCN) within a matter of several hours. Furthermore, the new particle formation was regularly followed by a measurable increase in particle mass concentration and extinction coefficient, indicative of a high abundance of condensable vapors in the atmosphere under study.
  • Item
    Primary versus secondary contributions to particle number concentrations in the European boundary layer
    (München : European Geopyhsical Union, 2011) Reddington, C.L.; Carslaw, K.S.; Spracklen, D.V.; Frontoso, M.G.; Collins, L.; Merikanto, J.; Minikin, A.; Hamburger, T.; Coe, H.; Kulmala, M.; Aalto, P.; Flentje, H.; Plass-Dülmer, C.; Birmili, W.; Wiedensohler, A.; Wehner, B.; Tuch, T.; Sonntag, A.; O'Dowd, C.D.; Jennings, S.G.; Dupuy, R.; Baltensperger, U.; Weingartner, E.; Hansson, H.-C.; Tunved, P.; Laj, P.; Sellegri, K.; Boulon, J.; Putaud, J.-P.; Gruening, C.; Swietlicki, E.; Roldin, P.; Henzing, J.S.; Moerman, M.; Mihalopoulos, N.; Kouvarakis, G.; Ždímal, V.; Zíková, N.; Marinoni, A.; Bonasoni, P.; Duchi, R.
    It is important to understand the relative contribution of primary and secondary particles to regional and global aerosol so that models can attribute aerosol radiative forcing to different sources. In large-scale models, there is considerable uncertainty associated with treatments of particle formation (nucleation) in the boundary layer (BL) and in the size distribution of emitted primary particles, leading to uncertainties in predicted cloud condensation nuclei (CCN) concentrations. Here we quantify how primary particle emissions and secondary particle formation influence size-resolved particle number concentrations in the BL using a global aerosol microphysics model and aircraft and ground site observations made during the May 2008 campaign of the European Integrated Project on Aerosol Cloud Climate Air Quality Interactions (EUCAARI). We tested four different parameterisations for BL nucleation and two assumptions for the emission size distribution of anthropogenic and wildfire carbonaceous particles. When we emit carbonaceous particles at small sizes (as recommended by the Aerosol Intercomparison project, AEROCOM), the spatial distributions of campaign-mean number concentrations of particles with diameter >50 nm (N50) and >100 nm (N100) were well captured by the model (R2≥0.8) and the normalised mean bias (NMB) was also small (−18% for N50 and −1% for N100). Emission of carbonaceous particles at larger sizes, which we consider to be more realistic for low spatial resolution global models, results in equally good correlation but larger bias (R2≥0.8, NMB = −52% and −29%), which could be partly but not entirely compensated by BL nucleation. Within the uncertainty of the observations and accounting for the uncertainty in the size of emitted primary particles, BL nucleation makes a statistically significant contribution to CCN-sized particles at less than a quarter of the ground sites. Our results show that a major source of uncertainty in CCN-sized particles in polluted European air is the emitted size of primary carbonaceous particles. New information is required not just from direct observations, but also to determine the "effective emission size" and composition of primary particles appropriate for different resolution models.
  • Item
    Modelling of sea salt concentrations over Europe: Key uncertainties and comparison with observations
    (München : European Geopyhsical Union, 2011) Tsyro, S.; Aas, W.; Soares, J.; Sofiev, M.; Berge, H.; Spindler, G.
    Sea salt aerosol can significantly affect the air quality. Sea salt can cause enhanced concentrations of particulate matter and change particle chemical composition, in particular in coastal areas, and therefore should be accounted for in air quality modelling. We have used an EMEP Unified model to calculate sea salt concentrations and depositions over Europe, focusing on studying the effects of uncertainties in sea salt production and lifetime on calculation results. Model calculations of sea salt have been compared with EMEP observations of sodium concentrations in air and precipitation for a four year period, from 2004 to 2007, including size (fine/coarse) resolved EMEP intensive measurements in 2006 and 2007. In the presented calculations, sodium air concentrations are between 8% and 46% overestimated, whereas concentrations in precipitation are systematically underestimated by 65–70% for years 2004–2007. A series of model tests have been performed to investigate the reasons for this underestimation, but further studies are needed. The model is found to reproduce the spatial distribution of Na+ in air and precipitation over Europe fairly well, and to capture most of sea salt episodes. The paper presents the main findings from a series of tests in which we compare several different sea spray source functions and also look at the effects of meteorological input and the efficiency of removal processes on calculated sea salt concentrations. Finally, sea salt calculations with the EMEP model have been compared with results from the SILAM model and observations for 2007. While the models produce quite close results for Na+ at the majority of 26 measurement sites, discrepancies in terms of bias and temporal correlation are also found. Those differences are believed to occur due to differences in the representation of source function and size distribution of sea salt aerosol, different meteorology used for model runs and the different models' resolution. This study contributes to getting a better insight on uncertainties associated with sea salt calculations and thus facilitates further improvement of aerosol modelling on both regional and global scales.
  • Item
    Changes in the production rate of secondary aerosol particles in Central Europe in view of decreasing SO2 emissions between 1996 and 2006
    (München : European Geopyhsical Union, 2010) Hamed, A.; Birmili, W.; Joutsensaari, J.; Mikkonen, S.; Asmi, A.; Wehner, B.; Spindler, G.; Jaatinen, A.; Wiedensohler, A.; Korhonen, H.; Lehtinen, K.E.J.; Laaksonen, A.
    In anthropogenically influenced atmospheres, sulphur dioxide (SO2) is the main precursor of gaseous sulphuric acid (H2SO4), which in turn is a main precursor for atmospheric particle nucleation. As a result of socio-economic changes, East Germany has seen a dramatic decrease in anthropogenic SO2 emissions between 1989 and present, as documented by routine air quality measurements in many locations. We have attempted to evaluate the influence of changing SO2 concentrations on the frequency and intensity of new particle formation (NPF) using two different data sets (1996–1997; 2003–2006) of experimental particle number size distributions (diameter range 3–750 nm) from the atmospheric research station Melpitz near Leipzig, Germany. Between the two periods SO2 concentrations decreased by 65% on average, while the frequency of NPF events dropped by 45%. Meanwhile, the average formation rate of 3 nm particles decreased by 68% on average. The trends were statistically significant and therefore suggest a connection between the availability of anthropogenic SO2 and freshly formed new particles. In contrast to the decrease in new particle formation, we found an increase in the mean growth rate of freshly nucleated particles (+22%), suggesting that particle nucleation and subsequent growth into larger sizes are delineated with respect to their precursor species. Using three basic parameters, the condensation sink for H2SO4, the SO2 concentration, and the global radiation intensity, we were able to define the characteristic range of atmospheric conditions under which particle formation events take place at the Melpitz site. While the decrease in the concentrations and formation rates of the new particles was rather evident, no similar decrease was found with respect to the generation of cloud condensation nuclei (CCN; particle diameter >100 nm) as a result of atmospheric nucleation events. On the contrary, the production of CCN following nucleation events appears to have increased by tens of percents. Our aerosol dynamics model simulations suggest that such an increase can be caused by the increased particle growth rate.
  • Item
    On the sub-micron aerosol size distribution in a coastal-rural site at El Arenosillo Station (SW – Spain)
    (München : European Geopyhsical Union, 2011) Sorribas, M.; de la Morena, B.A.; Wehner, B.; López, J.F.; Prats, N.; Mogo, S.; Wiedensohler, A.; Cachorro, V.E.
    This study focuses on the analysis of the sub-micron aerosol characteristics at El Arenosillo Station, a rural and coastal environment in South-western Spain between 1 August 2004 and 31 July 2006 (594 days). The mean total concentration (NT) was 8660 cm−3 and the mean concentrations in the nucleation (NNUC), Aitken (NAIT) and accumulation (NACC) particle size ranges were 2830 cm−3, 4110 cm−3 and 1720 cm−3, respectively. Median size distribution was characterised by a single-modal fit, with a geometric diameter, median number concentration and geometric standard deviation of 60 nm, 5390 cm−3 and 2.31, respectively. Characterisation of primary emissions, secondary particle formation, changes to meteorology and long-term transport has been necessary to understand the seasonal and annual variability of the total and modal particle concentration. Number concentrations exhibited a diurnal pattern with maximum concentrations around noon. This was governed by the concentrations of the nucleation and Aitken modes during the warm seasons and only by the nucleation mode during the cold seasons. Similar monthly mean total concentrations were observed throughout the year due to a clear inverse variation between the monthly mean NNUC and NACC. It was related to the impact of desert dust and continental air masses on the monthly mean particle levels. These air masses were associated with high values of NACC which suppressed the new particle formation (decreasing NNUC). Each day was classified according to a land breeze flow or a synoptic pattern influence. The median size distribution for desert dust and continental aerosol was dominated by the Aitken and accumulation modes, and marine air masses were dominated by the nucleation and Aitken modes. Particles moved offshore due to the land breeze and had an impact on the particle burden at noon, especially when the wind was blowing from the NW sector in the morning during summer time. This increased NNUC and NAIT by factors of 3.1 and 2.4, respectively. Nucleation events with the typical "banana" shape were characterised by a mean particle nucleation rate of 0.74 cm−3 s−1, a mean growth rate of 1.96 nm h−1 and a mean total duration of 9.25 h (starting at 10:55 GMT and ending at 20:10 GMT). They were observed for 48 days. Other nucleation events were identified as those produced by the emissions from the industrial areas located at a distance of 35 km. They were observed for 42 days. Both nucleation events were strongly linked to the marine air mass origin.