2011, Hartmann, S., Niedermeier, D., Voigtländer, J., Clauss, T., Shaw, R.A., Wex, H., Kiselev, A., Stratmann, F.

At the Leipzig Aerosol Cloud Interaction Simulator (LACIS) experiments investigating homogeneous and heterogeneous nucleation of ice (particularly immersion freezing in the latter case) have been carried out. Here both the physical LACIS setup and the numerical model developed to design experiments at LACIS and interpret their results are presented in detail. Combining results from the numerical model with experimental data, it was found that for the experimental parameter space considered, classical homogeneous ice nucleation theory is able to predict the freezing behavior of highly diluted ammonium sulfate solution droplets, while classical heterogeneous ice nucleation theory, together with the assumption of a constant contact angle, fails to predict the immersion freezing behavior of surrogate mineral dust particles (Arizona Test Dust, ATD). The main reason for this failure is the compared to experimental data apparently overly strong temperature dependence of the nucleation rate coefficient. Assuming, in the numerical model, Classical Nucleation Theory (CNT) for homogeneous ice nucleation and a CNT-based parameterization for the nucleation rate coefficient in the immersion freezing mode, recently published by our group, it was found that even for a relatively effective ice nucleating agent such as pure ATD, there is a temperature range where homogeneous ice nucleation is dominant. The main explanation is the apparently different temperature dependencies of the two freezing mechanisms. Finally, reviewing the assumptions made during the derivation of the CNT-based parameterization for immersion freezing, it was found that the assumption of constant temperature during ice nucleation and the chosen ice nucleation time were justified, underlining the applicability of the method to determine the fitting coefficients in the parameterization equation.

2011, Meinshausen, M., Raper, S.C.B., Wigley, T.M.L.

Current scientific knowledge on the future response of the climate system to human-induced perturbations is comprehensively captured by various model intercomparison efforts. In the preparation of the Fourth Assessment Report (AR4) of the Intergovernmental Panel on Climate Change (IPCC), intercomparisons were organized for atmosphere-ocean general circulation models (AOGCMs) and carbon cycle models, named "CMIP3" and "C4MIP", respectively. Despite their tremendous value for the scientific community and policy makers alike, there are some difficulties in interpreting the results. For example, radiative forcings were not standardized across the various AOGCM integrations and carbon cycle runs, and, in some models, key forcings were omitted. Furthermore, the AOGCM analysis of plausible emissions pathways was restricted to only three SRES scenarios. This study attempts to address these issues. We present an updated version of MAGICC, the simple carbon cycle-climate model used in past IPCC Assessment Reports with enhanced representation of time-varying climate sensitivities, carbon cycle feedbacks, aerosol forcings and ocean heat uptake characteristics. This new version, MAGICC6, is successfully calibrated against the higher complexity AOGCMs and carbon cycle models. Parameterizations of MAGICC6 are provided. The mean of the emulations presented here using MAGICC6 deviates from the mean AOGCM responses by only 2.2% on average for the SRES scenarios. This enhanced emulation skill in comparison to previous calibrations is primarily due to: making a "like-with-like comparison" using AOGCM-specific subsets of forcings; employing a new calibration procedure; as well as the fact that the updated simple climate model can now successfully emulate some of the climate-state dependent effective climate sensitivities of AOGCMs. The diagnosed effective climate sensitivity at the time of CO2 doubling for the AOGCMs is on average 2.88 °C, about 0.33 °C cooler than the mean of the reported slab ocean climate sensitivities. In the companion paper (Part 2) of this study, we examine the combined climate system and carbon cycle emulations for the complete range of IPCC SRES emissions scenarios and the new RCP pathways.

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.

2011, Winkelmann, R., Martin, M.A., Haseloff, M., Albrecht, T., Bueler, E., Khroulev, C., Levermann, A.

We present the Potsdam Parallel Ice Sheet Model (PISM-PIK), developed at the Potsdam Institute for Climate Impact Research to be used for simulations of large-scale ice sheet-shelf systems. It is derived from the Parallel Ice Sheet Model (Bueler and Brown, 2009). Velocities are calculated by superposition of two shallow stress balance approximations within the entire ice covered region: the shallow ice approximation (SIA) is dominant in grounded regions and accounts for shear deformation parallel to the geoid. The plug-flow type shallow shelf approximation (SSA) dominates the velocity field in ice shelf regions and serves as a basal sliding velocity in grounded regions. Ice streams can be identified diagnostically as regions with a significant contribution of membrane stresses to the local momentum balance. All lateral boundaries in PISM-PIK are free to evolve, including the grounding line and ice fronts. Ice shelf margins in particular are modeled using Neumann boundary conditions for the SSA equations, reflecting a hydrostatic stress imbalance along the vertical calving face. The ice front position is modeled using a subgrid-scale representation of calving front motion (Albrecht et al., 2011) and a physically-motivated calving law based on horizontal spreading rates. The model is tested in experiments from the Marine Ice Sheet Model Intercomparison Project (MISMIP). A dynamic equilibrium simulation of Antarctica under present-day conditions is presented in Martin et al. (2011).

2011, Meinshausen, M., Raper, S.C.B., Wigley, T.M.L.

Intercomparisons of coupled atmosphere-ocean general circulation models (AOGCMs) and carbon cycle models are important for galvanizing our current scientific knowledge to project future climate. Interpreting such intercomparisons faces major challenges, not least because different models have been forced with different sets of forcing agents. Here, we show how an emulation approach with MAGICC6 can address such problems. In a companion paper (Meinshausen et al., 2011a), we show how the lower complexity carbon cycle-climate model MAGICC6 can be calibrated to emulate, with considerable accuracy, globally aggregated characteristics of these more complex models. Building on that, we examine here the Coupled Model Intercomparison Project's Phase 3 results (CMIP3). If forcing agents missed by individual AOGCMs in CMIP3 are considered, this reduces ensemble average temperature change from pre-industrial times to 2100 under SRES A1B by 0.4 °C. Differences in the results from the 1980 to 1999 base period (as reported in IPCC AR4) to 2100 are negligible, however, although there are some differences in the trajectories over the 21st century. In a second part of this study, we consider the new RCP scenarios that are to be investigated under the forthcoming CMIP5 intercomparison for the IPCC Fifth Assessment Report. For the highest scenario, RCP8.5, relative to pre-industrial levels, we project a median warming of around 4.6 °C by 2100 and more than 7 °C by 2300. For the lowest RCP scenario, RCP3-PD, the corresponding warming is around 1.5 °C by 2100, decreasing to around 1.1 °C by 2300 based on our AOGCM and carbon cycle model emulations. Implied cumulative CO2 emissions over the 21st century for RCP8.5 and RCP3-PD are 1881 GtC (1697 to 2034 GtC, 80% uncertainty range) and 381 GtC (334 to 488 GtC), when prescribing CO2 concentrations and accounting for uncertainty in the carbon cycle. Lastly, we assess the reasons why a previous MAGICC version (4.2) used in IPCC AR4 gave roughly 10% larger warmings over the 21st century compared to the CMIP3 average. We find that forcing differences and the use of slightly too high climate sensitivities inferred from idealized high-forcing runs were the major reasons for this difference.

2010, Spracklen, D.V., Carslaw, K.S., Merikanto, J., Mann, G.W., Reddington, C.L., Pickering, S., Ogren, J.A., Andrews, E., Baltensperger, U., Weingartner, E., Boy, M., Kulmala, M., Laakso, L., Lihavainen, H., Kivekäs, N., Komppula, M., Mihalopoulos, N., Kouvarakis, G., Jennings, S.G., O'Dowd, C., Birmili, W., Wiedensohler, A., Weller, R., Gras, J., Laj, P., Sellegri, K., Bonn, B., Krejci, R., Laaksonen, A., Hamed, A., Minikin, A., Harrison, R.M., Talbot, R., Sun, J.

We synthesised observations of total particle number (CN) concentration from 36 sites around the world. We found that annual mean CN concentrations are typically 300–2000 cm−3 in the marine boundary layer and free troposphere (FT) and 1000–10 000 cm−3 in the continental boundary layer (BL). Many sites exhibit pronounced seasonality with summer time concentrations a factor of 2–10 greater than wintertime concentrations. We used these CN observations to evaluate primary and secondary sources of particle number in a global aerosol microphysics model. We found that emissions of primary particles can reasonably reproduce the spatial pattern of observed CN concentration (R2=0.46) but fail to explain the observed seasonal cycle (R2=0.1). The modeled CN concentration in the FT was biased low (normalised mean bias, NMB=−88%) unless a secondary source of particles was included, for example from binary homogeneous nucleation of sulfuric acid and water (NMB=−25%). Simulated CN concentrations in the continental BL were also biased low (NMB=−74%) unless the number emission of anthropogenic primary particles was increased or a mechanism that results in particle formation in the BL was included. We ran a number of simulations where we included an empirical BL nucleation mechanism either using the activation-type mechanism (nucleation rate, J, proportional to gas-phase sulfuric acid concentration to the power one) or kinetic-type mechanism (J proportional to sulfuric acid to the power two) with a range of nucleation coefficients. We found that the seasonal CN cycle observed at continental BL sites was better simulated by BL particle formation (R2=0.3) than by increasing the number emission from primary anthropogenic sources (R2=0.18). The nucleation constants that resulted in best overall match between model and observed CN concentrations were consistent with values derived in previous studies from detailed case studies at individual sites. In our model, kinetic and activation-type nucleation parameterizations gave similar agreement with observed monthly mean CN concentrations.

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.

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.

2011, Fountoukis, C., Racherla, P.N., Denier van der Gon, H.A.C., Polymeneas, P., Charalampidis, P.E., Pilinis, C., Wiedensohler, A., Dall'Osto, M., O'Dowd, C., Pandis, S.N.

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.

2011, Albrecht, T., Martin, M., Haseloff, M., Winkelmann, R., Levermann, A.

In order to explore the response of the Greenland ice sheet (GIS) to climate change on long (centennial to multi-millennial) time scales, a regional energy-moisture balance model has been developed. This model simulates seasonal variations of temperature and precipitation over Greenland and explicitly accounts for elevation and albedo feedbacks. From these fields, the annual mean surface temperature and surface mass balance can be determined and used to force an ice sheet model. The melt component of the surface mass balance is computed here using both a positive degree day approach and a more physically-based alternative that includes insolation and albedo explicitly. As a validation of the climate model, we first simulated temperature and precipitation over Greenland for the prescribed, present-day topography. Our simulated climatology compares well to observations and does not differ significantly from that of a simple parameterization used in many previous simulations. Furthermore, the calculated surface mass balance using both melt schemes falls within the range of recent regional climate model results. For a prescribed, ice-free state, the differences in simulated climatology between the regional energy-moisture balance model and the simple parameterization become significant, with our model showing much stronger summer warming. When coupled to a three-dimensional ice sheet model and initialized with present-day conditions, the two melt schemes both allow realistic simulations of the present-day GIS.