Search Results

Now showing 1 - 8 of 8
  • Item
    Organic aerosol concentration and composition over Europe: Insights from comparison of regional model predictions with aerosol mass spectrometer factor analysis
    (München : European Geopyhsical Union, 2014) Fountoukis, C.; Megaritis, A.G.; Skyllakou, K.; Charalampidis, P.E.; Pilinis, C.; van der Gon, H.A.C. Denier; Crippa, M.; Canonaco, F.; Mohr, C.; Prévôt, A.S.H.; Allan, J.D.; Poulain, L.; Petäjä, T.; Tiitta, P.; Carbone, S.; Kiendler-Scharr, A.; Nemitz, E.; O'Dowd, C.; Swietlicki, E.; Pandis, S.N.
    A detailed three-dimensional regional chemical transport model (Particulate Matter Comprehensive Air Quality Model with Extensions, PMCAMx) was applied over Europe, focusing on the formation and chemical transformation of organic matter. Three periods representative of different seasons were simulated, corresponding to intensive field campaigns. An extensive set of AMS measurements was used to evaluate the model and, using factor-analysis results, gain more insight into the sources and transformations of organic aerosol (OA). Overall, the agreement between predictions and measurements for OA concentration is encouraging, with the model reproducing two-thirds of the data (daily average mass concentrations) within a factor of 2. Oxygenated OA (OOA) is predicted to contribute 93% to total OA during May, 87% during winter and 96% during autumn, with the rest consisting of fresh primary OA (POA). Predicted OOA concentrations compare well with the observed OOA values for all periods, with an average fractional error of 0.53 and a bias equal to −0.07 (mean error = 0.9 μg m−3, mean bias = −0.2 μg m−3). The model systematically underpredicts fresh POA at most sites during late spring and autumn (mean bias up to −0.8 μg m−3). Based on results from a source apportionment algorithm running in parallel with PMCAMx, most of the POA originates from biomass burning (fires and residential wood combustion), and therefore biomass burning OA is most likely underestimated in the emission inventory. The sensitivity of POA predictions to the corresponding emissions' volatility distribution is discussed. The model performs well at all sites when the Positive Matrix Factorization (PMF)-estimated low-volatility OOA is compared against the OA with saturation concentrations of the OA surrogate species C* ≤ 0.1 μg m−3 and semivolatile OOA against the OA with C* > 0.1 μg m−3.
  • Item
    Simulating ultrafine particle formation in Europe using a regional CTM: Contribution of primary emissions versus secondary formation to aerosol number concentrations
    (München : European Geopyhsical Union, 2012) Fountoukis, C.; Riipinen, I.; Denier van der Gon, H.A.C.; Charalampidis, P.E.; Pilinis, C.; Wiedensohler, A.; O'Dowd, C.; Putaud, J.P.; Moerman, M.; Pandis, S.N.
    A three-dimensional regional chemical transport model (CTM) with detailed aerosol microphysics, PMCAMx-UF, was applied to the European domain to simulate the contribution of direct emissions and secondary formation to total particle number concentrations during May 2008. PMCAMx-UF uses the Dynamic Model for Aerosol Nucleation and the Two-Moment Aerosol Sectional (TOMAS) algorithm to track both aerosol number and mass concentration using a sectional approach. The model predicts nucleation events that occur over scales of hundreds up to thousands of kilometers especially over the Balkans and Southeast Europe. The model predictions were compared against measurements from 7 sites across Europe. The model reproduces more than 70% of the hourly concentrations of particles larger than 10 nm (N10) within a factor of 2. About half of these particles are predicted to originate from nucleation in the lower troposphere. Regional nucleation is predicted to increase the total particle number concentration by approximately a factor of 3. For particles larger than 100 nm the effect varies from an increase of 20% in the eastern Mediterranean to a decrease of 20% in southern Spain and Portugal resulting in a small average increase of around 1% over the whole domain. Nucleation has a significant effect in the predicted N50 levels (up to a factor of 2 increase) mainly in areas where there are condensable vapors to grow the particles to larger sizes. A semi-empirical ternary sulfuric acid-ammonia-water parameterization performs better than the activation or the kinetic parameterizations in reproducing the observations. Reducing emissions of ammonia and sulfur dioxide affects certain parts of the number size distribution.
  • Item
    Simulating the formation of carbonaceous aerosol in a European Megacity (Paris) during the MEGAPOLI summer and winter campaigns
    (München : European Geopyhsical Union, 2016) Fountoukis, Christos; Megaritis, Athanasios G.; Skyllakou, Ksakousti; Charalampidis, Panagiotis E.; van der Gon, Hugo A.C.Denier; Crippa, Monica; Prévôt, André S.H.; Fachinger, Friederike; Wiedensohler, Alfred; Pilinis, Christodoulos; Pandis, Spyros N.
    We use a three-dimensional regional chemical transport model (PMCAMx) with high grid resolution and high-resolution emissions (4 × 4 km2) over the Paris greater area to simulate the formation of carbonaceous aerosol during a summer (July 2009) and a winter (January/February 2010) period as part of the MEGAPOLI (megacities: emissions, urban, regional, and global atmospheric pollution and climate effects, and Integrated tools for assessment and mitigation) campaigns. Model predictions of carbonaceous aerosol are compared against Aerodyne aerosol mass spectrometer and black carbon (BC) high time resolution measurements from three ground sites. PMCAMx predicts BC concentrations reasonably well reproducing the majority (70 %) of the hourly data within a factor of two during both periods. The agreement for the summertime secondary organic aerosol (OA) concentrations is also encouraging (mean bias = 0.1 µg m−3) during a photochemically intense period. The model tends to underpredict the summertime primary OA concentrations in the Paris greater area (by approximately 0.8 µg m−3) mainly due to missing primary OA emissions from cooking activities. The total cooking emissions are estimated to be approximately 80 mg d−1 per capita and have a distinct diurnal profile in which 50 % of the daily cooking OA is emitted during lunch time (12:00–14:00 LT) and 20 % during dinner time (20:00–22:00 LT). Results also show a large underestimation of secondary OA in the Paris greater area during wintertime (mean bias =  −2.3 µg m−3) pointing towards a secondary OA formation process during low photochemical activity periods that is not simulated in the model.
  • Item
    Albedo and heat transport in 3-D model simulations of the early Archean climate
    (München : European Geopyhsical Union, 2013) Kienert, H.; Feulner, G.; Petoukhov, V.
    At the beginning of the Archean eon (ca. 3.8 billion years ago), the Earth's climate state was significantly different from today due to the lower solar luminosity, smaller continental fraction, higher rotation rate and, presumably, significantly larger greenhouse gas concentrations. All these aspects play a role in solutions to the "faint young Sun paradox" which must explain why the ocean surface was not fully frozen at that time. Here, we present 3-D model simulations of climate states that are consistent with early Archean boundary conditions and have different CO2 concentrations, aiming at an understanding of the fundamental characteristics of the early Archean climate system. In order to do so, we have appropriately modified an intermediate complexity climate model that couples a statistical-dynamical atmosphere model (involving parameterizations of the dynamics) to an ocean general circulation model and a thermodynamic-dynamic sea-ice model. We focus on three states: one of them is ice-free, one has the same mean surface air temperature of 288 K as today's Earth and the third one is the coldest stable state in which there is still an area with liquid surface water (i.e. the critical state at the transition to a "snowball Earth"). We find a reduction in meridional heat transport compared to today, which leads to a steeper latitudinal temperature profile and has atmospheric as well as oceanic contributions. Ocean surface velocities are largely zonal, and the strength of the atmospheric meridional circulation is significantly reduced in all three states. These aspects contribute to the observed relation between global mean temperature and albedo, which we suggest as a parameterization of the ice-albedo feedback for 1-D model simulations of the early Archean and thus the faint young Sun problem.
  • Item
    Heinrich event 1: An example of dynamical ice-sheet reaction to oceanic changes
    (München : European Geopyhsical Union, 2011) Álvarez-Solas, J.; Montoya, M.; Ritz, C.; Ramstein, G.; Charbit, S.; Dumas, C.; Nisancioglu, K.; Dokken, T.; Ganopolski, A.
    Heinrich events, identified as enhanced ice-rafted detritus (IRD) in North Atlantic deep sea sediments (Heinrich, 1988; Hemming, 2004) have classically been attributed to Laurentide ice-sheet (LIS) instabilities (MacAyeal, 1993; Calov et al., 2002; Hulbe et al., 2004) and assumed to lead to important disruptions of the Atlantic meridional overturning circulation (AMOC) and North Atlantic deep water (NADW) formation. However, recent paleoclimate data have revealed that most of these events probably occurred after the AMOC had already slowed down or/and NADW largely collapsed, within about a thousand years (Hall et al., 2006; Hemming, 2004; Jonkers et al., 2010; Roche et al., 2004), implying that the initial AMOC reduction could not have been caused by the Heinrich events themselves. Here we propose an alternative driving mechanism, specifically for Heinrich event 1 (H1; 18 to 15 ka BP), by which North Atlantic ocean circulation changes are found to have strong impacts on LIS dynamics. By combining simulations with a coupled climate model and a three-dimensional ice sheet model, our study illustrates how reduced NADW and AMOC weakening lead to a subsurface warming in the Nordic and Labrador Seas resulting in rapid melting of the Hudson Strait and Labrador ice shelves. Lack of buttressing by the ice shelves implies a substantial ice-stream acceleration, enhanced ice-discharge and sea level rise, with peak values 500–1500 yr after the initial AMOC reduction. Our scenario modifies the previous paradigm of H1 by solving the paradox of its occurrence during a cold surface period, and highlights the importance of taking into account the effects of oceanic circulation on ice-sheets dynamics in order to elucidate the triggering mechanism of Heinrich events.
  • Item
    Simulation of the last glacial cycle with a coupled climate ice-sheet model of intermediate complexity
    (München : European Geopyhsical Union, 2010) Ganopolski, A.; Calov, R.; Claussen, M.
    A new version of the Earth system model of intermediate complexity, CLIMBER-2, which includes the three-dimensional polythermal ice-sheet model SICOPOLIS, is used to simulate the last glacial cycle forced by variations of the Earth's orbital parameters and atmospheric concentration of major greenhouse gases. The climate and ice-sheet components of the model are coupled bi-directionally through a physically-based surface energy and mass balance interface. The model accounts for the time-dependent effect of aeolian dust on planetary and snow albedo. The model successfully simulates the temporal and spatial dynamics of the major Northern Hemisphere (NH) ice sheets, including rapid glacial inception and strong asymmetry between the ice-sheet growth phase and glacial termination. Spatial extent and elevation of the ice sheets during the last glacial maximum agree reasonably well with palaeoclimate reconstructions. A suite of sensitivity experiments demonstrates that simulated ice-sheet evolution during the last glacial cycle is very sensitive to some parameters of the surface energy and mass-balance interface and dust module. The possibility of a considerable acceleration of the climate ice-sheet model is discussed.
  • Item
    ASAMgpu V1.0 - A moist fully compressible atmospheric model using graphics processing units (GPUs)
    (München : European Geopyhsical Union, 2012) Horn, S.
    In this work the three dimensional compressible moist atmospheric model ASAMgpu is presented. The calculations are done using graphics processing units (GPUs). To ensure platform independence OpenGL and GLSL are used, with that the model runs on any hardware supporting fragment shaders. The MPICH2 library enables interprocess communication allowing the usage of more than one GPU through domain decomposition. Time integration is done with an explicit three step Runge-Kutta scheme with a time-splitting algorithm for the acoustic waves. The results for four test cases are shown in this paper. A rising dry heat bubble, a cold bubble induced density flow, a rising moist heat bubble in a saturated environment, and a DYCOMS-II case.
  • Item
    The Potsdam Parallel Ice Sheet Model (PISM-PIK) - Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet
    (München : European Geopyhsical Union, 2011) Martin, M.A.; Winkelmann, R.; Haseloff, M.; Albrecht, T.; Bueler, E.; Khroulev, C.; Levermann, A.
    We present a dynamic equilibrium simulation of the ice sheet-shelf system on Antarctica with the Potsdam Parallel Ice Sheet Model (PISM-PIK). The simulation is initialized with present-day conditions for bed topography and ice thickness and then run to steady state with constant present-day surface mass balance. Surface temperature and sub-shelf basal melt distribution are parameterized. Grounding lines and calving fronts are free to evolve, and their modeled equilibrium state is compared to observational data. A physically-motivated calving law based on horizontal spreading rates allows for realistic calving fronts for various types of shelves. Steady-state dynamics including surface velocity and ice flux are analyzed for whole Antarctica and the Ronne-Filchner and Ross ice shelf areas in particular. The results show that the different flow regimes in sheet and shelves, and the transition zone between them, are captured reasonably well, supporting the approach of superposition of SIA and SSA for the representation of fast motion of grounded ice. This approach also leads to a natural emergence of sliding-dominated flow in stream-like features in this new 3-D marine ice sheet model.