Search Results

Now showing 1 - 10 of 15
  • 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
    Three-dimensional evolution of Saharan dust transport towards Europe based on a 9-year EARLINET-optimized CALIPSO dataset
    (Katlenburg-Lindau : EGU, 2017) Marinou, Eleni; Amiridis, Vassilis; Binietoglou, Ioannis; Tsikerdekis, Athanasios; Solomos, Stavros; Proestakis, Emannouil; Konsta, Dimitra; Papagiannopoulos, Nikolaos; Tsekeri, Alexandra; Vlastou, Georgia; Zanis, Prodromos; Balis, Dimitrios; Wandinger, Ulla; Ansmann, Albert
    In this study we use a new dust product developed using CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) observations and EARLINET (European Aerosol Research Lidar Network) measurements and methods to provide a 3-D multiyear analysis on the evolution of Saharan dust over North Africa and Europe. The product uses a CALIPSO L2 backscatter product corrected with a depolarization-based method to separate pure dust in external aerosol mixtures and a Saharan dust lidar ratio (LR) based on long-term EARLINET measurements to calculate the dust extinction profiles. The methodology is applied on a 9-year CALIPSO dataset (2007-2015) and the results are analyzed here to reveal for the first time the 3-D dust evolution and the seasonal patterns of dust over its transportation paths from the Sahara towards the Mediterranean and Continental Europe. During spring, the spatial distribution of dust shows a uniform pattern over the Sahara desert. The dust transport over the Mediterranean Sea results in mean dust optical depth (DOD) values up to 0.1. During summer, the dust activity is mostly shifted to the western part of the desert where mean DOD near the source is up to 0.6. Elevated dust plumes with mean extinction values between 10 and 75 Mm-1 are observed throughout the year at various heights between 2 and 6 km, extending up to latitudes of 40° N. Dust advection is identified even at latitudes of about 60° N, but this is due to rare events of episodic nature. Dust plumes of high DOD are also observed above the Balkans during the winter period and above northwest Europe during autumn at heights between 2 and 4 km, reaching mean extinction values up to 50 Mm-1. The dataset is considered unique with respect to its potential applications, including the evaluation of dust transport models and the estimation of cloud condensation nuclei (CCN) and ice nuclei (IN) concentration profiles. Finally, the product can be used to study dust dynamics during transportation, since it is capable of revealing even fine dynamical features such as the particle uplifting and deposition on European mountainous ridges such as the Alps and Carpathian Mountains.
  • Item
    Compiling geophysical and geological information into a 3-D model of the glacially-affected island of Föhr
    (Munich : EGU, 2012) Burschil, T.; Scheer, W.; Kirsch, R.; Wiederhold, H.
    Within the scope of climatic change and associated sea level rise, coastal aquifers are endangered and are becoming more a focus of research to ensure the future water supply in coastal areas. For groundwater modelling a good understanding of the geological/hydrogeological situation and the aquifer behavior is necessary. In preparation of groundwater modelling and assessment of climate change impacts on coastal water resources, we setup a geological/hydrogeological model for the North Sea Island of Föhr. Data from different geophysical methods applied from the air, the surface and in boreholes contribute to the 3-D model, e.g. airborne electromagnetics (SkyTEM) for spatial mapping the resistivity of the entire island, seismic reflections for detailed cross-sections in the groundwater catchment area, and geophysical borehole logging for calibration of these measurements. An iterative and integrated evaluation of the results from the different geophysical methods contributes to reliable data as input for the 3-D model covering the whole island and not just the well fields. The complex subsurface structure of the island is revealed. The local waterworks use a freshwater body embedded in saline groundwater. Several glaciations reordered the youngest Tertiary and Quaternary sediments by glaciotectonic thrust faulting, as well as incision and refill of glacial valleys. Both subsurface structures have a strong impact on the distribution of freshwater-bearing aquifers. A digital geological 3-D model reproduces the hydrogeological structure of the island as a base for a groundwater model. In the course of the data interpretation, we deliver a basis for rock identification. We demonstrate that geophysical investigation provide petrophysical parameters and improve the understanding of the subsurface and the groundwater system. The main benefit of our work is that the successful combination of electromagnetic, seismic and borehole data reveals the complex geology of a glacially-affected island. A sound understanding of the subsurface structure and the compilation of a 3-D model is imperative and the basis for a groundwater flow model to predict climate change effects on future water resources.
  • Item
    Three-dimensional monitoring of soil water content in a maize field using Electrical Resistivity Tomography
    (Munich : EGU, 2013) Beff, L.; Günther, T.; Vandoorne, B.; Couvreur, V.; Javaux, M.
    A good understanding of the soil water content (SWC) distribution at the field scale is essential to improve the management of water, soil and crops. Recent studies proved that Electrical Resistivity Tomography (ERT) opens interesting perspectives in the determination of the SWC distribution in 3 dimensions (3-D). This study was conducted (i) to check and validate how ERT is able to monitor SWC distribution in a maize field during the late growing season; and (ii) to investigate how maize plants and rainfall affect the dynamics of SWC distribution. Time Domain Reflectometry (TDR) measurements were used to validate ERT-inverted SWC values. Evolution of water mass balance was also calculated to check whether ERT was capable of giving a reliable estimate of soil water stock evolution. It is observed that ERT was able to give the same average SWC as TDR (R2 = 0.98). In addition, ERT gives better estimates of the water stock than TDR thanks to its higher spatial resolution. The high resolution of ERT measurements also allows for the discrimination of SWC heterogeneities. The SWC distribution showed that alternation of maize rows and inter-rows was the main influencing factor of the SWC distribution. The drying patterns were linked to the root profiles, with drier zones under the maize rows. During short periods, with negligible rainfall, the SWC decrease took place mainly in the two upper soil horizons and in the inter-row area. In contrast, rainfall increased the SWC mostly under the maize rows and in the upper soil layer. Nevertheless, the total amount of rainfall during the growing season was not sufficient to modify the SWC patterns induced by the maize rows. During the experimental time, there was hardly any SWC redistribution from maize rows to inter-rows. Yet, lateral redistribution from inter-rows to maize rows induced by potential gradient generates SWC decrease in the inter-row area and in the deeper soil horizons.
  • 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
    Can VHF radars at polar latitudes measure mean vertical winds in the presence of PMSE?
    (Göttingen : Copernicus GmbH, 2019) Gudadze, N.; Stober, G.; Chau, J.L.
    Mean vertical velocity measurements obtained from radars at polar latitudes using polar mesosphere summer echoes (PMSEs) as an inert tracer have been considered to be non-representative of the mean vertical winds over the last couple of decades. We used PMSEs observed with the Middle Atmosphere Alomar Radar System (MAARSY) over Andøya, Norway (69.30°N, 16.04°E), during summers of 2016 and 2017 to derive mean vertical winds in the upper mesosphere. The 3-D vector wind components (zonal, meridional and vertical) are based on a Doppler beam swinging experiment using five beam directions (one vertical and four oblique). The 3-D wind components are computed using a recently developed wind retrieval technique. The method includes full non-linear error propagation, spatial and temporal regularisation, and beam pointing corrections and angular pointing uncertainties. Measurement uncertainties are used as weights to obtain seasonal weighted averages and characterise seasonal mean vertical velocities. Weighted average values of vertical velocities reveal a weak upward behaviour at altitudes ∼ 84-87 km after eliminating the influence of the speed of falling ice. At the same time, a sharp decrease (increase) in the mean vertical velocities at the lower (upper) edges of the summer mean altitude profile, which are attributed to the sampling issues of the PMSE due to disappearance of the target corresponding to the certain regions of motions and temperatures, prevails. Thus the mean vertical velocities can be biased downwards at the lower edge, and the mean vertical velocities can be biased upwards at the upper edge, while at the main central region the obtained mean vertical velocities are consistent with expected upward values of mean vertical winds after considering ice particle sedimentation. © 2019 Author(s). This work is distributed under the Creative Commons Attribution 4.0 License.