ALBERT

Description: Influence of surface morphology on the immersion mode ice nucleation efficiency of hematite particles Atmospheric Chemistry and Physics Discussions, 13, 23757-23780, 2013 Author(s): N. Hiranuma, N. Hoffmann, A. Kiselev, A. Dreyer, K. Zhang, G. Kulkarni, T. Koop, and O. Möhler In this paper, the effect of the morphological modification of aerosol particles with respect to heterogeneous ice nucleation is comprehensively investigated for laboratory-generated hematite particles as a model substrate for atmospheric dust particles. The surface area-scaled ice nucleation efficiencies of monodisperse cubic hematite particles and milled hematite particles were measured with a series of expansion cooling experiments using the Aerosol Interaction and Dynamics in the Atmosphere (AIDA) cloud simulation chamber. Complementary off-line characterization of physico-chemical properties of both hematite subsets were also carried out with scanning electron microscopy (SEM), energy dispersive X-ray (EDX) spectroscopy, dynamic light scattering (DLS), and an electro-kinetic particle charge detector to further constrain droplet-freezing measurements of hematite particles. Additionally, an empirical parameterization derived from our laboratory measurements was implemented in the single-column version of the Community Atmospheric Model version 5 (CAM5) to investigate the model sensitivity in simulated ice crystal number concentration on different ice nucleation efficiencies. From an experimental perspective, our results show that the immersion mode ice nucleation efficiency of milled hematite particles is almost an order of magnitude higher at −35.2 °C 〈 T 〈 −33.5 °C than that of the cubic hematite particles, indicating a substantial effect of morphological irregularities on immersion mode freezing. Our modeling results similarly show that the increased droplet-freezing rates of milled hematite particles lead to about one order magnitude higher ice crystal number in the upper troposphere than cubic hematite particles. Overall, our results suggest that the surface irregularities and associated active sites lead to greater ice activation through droplet-freezing.

Description: Fracture-induced softening for large-scale ice dynamics The Cryosphere Discussions, 7, 4501-4544, 2013 Author(s): T. Albrecht and A. Levermann Floating ice shelves can exert a retentive and hence stabilizing force onto the inland ice sheet of Antarctica. However, this effect has been observed to diminish by fracture-coupled dynamic processes within the protective ice shelves leading to accelerated ice flow and hence to a sea-level contribution. In order to better understand the role of fractures in ice dynamics we apply a large-scale continuum representation of fractures and related fracture growth into the prognostic Parallel Ice Sheet Model (PISM). To this end we introduce a higher-order accuracy advection scheme for the transport of the two-dimensional fracture density across the regular computational grid. Dynamic coupling of fractures and ice flow is attained by a reduction of effective ice viscosity proportional to the inferred fracture density. This formulation implies the possibility of a non-linear threshold behavior due to self-amplified fracturing in shear regions triggered by small variations in damage threshold. As a result of prognostic flow simulations, flow patterns with realistically large across-flow velocity gradients in fracture-weakened regions as seen in observations are reproduced. This model framework is expandable to grounded ice streams and accounts for climate-induced effects on fracturing and hence on the ice-flow dynamics. It further allows for an enhanced fracture-based calving parameterization.

Description: Magnetic field and dynamic pressure ULF fluctuations in coronal-mass-ejection-driven sheath regions Annales Geophysicae, 31, 1559-1567, 2013 Author(s): E. K. J. Kilpua, H. Hietala, H. E. J. Koskinen, D. Fontaine, and L. Turc Compressed sheath regions form ahead of interplanetary coronal mass ejections (ICMEs) that are sufficiently faster than the preceding solar wind. The turbulent sheath regions are important drivers of magnetospheric activity, but due to their complex internal structure, relatively little is known on the distribution of the magnetic field and plasma variations in them. In this paper we investigate ultra low frequency (ULF) fluctuations in the interplanetary magnetic field (IMF) and in dynamic pressure ( P dyn ) using a superposed epoch analysis of 41 sheath regions observed during solar cycle 23. We find strongest fluctuation power near the shock and in the vicinity of the ICME leading edge. The IMF and P dyn ULF power have different profiles within the sheath; the former is enhanced in the leading part of the sheath, while the latter is increased in the trailing part of the sheath. We also find that the ICME properties affect the level and distribution of the ULF power in sheath regions. For example, sheath regions associated with strong or fast ICMEs, or those that are crossed at intermediate distances from the center, have strongest ULF power and large variation in the power throughout the sheath region. The weaker or slower ICMEs, or those that are crossed centrally, have in general considerably weaker ULF power with relatively smooth profiles. The strong and abrupt decrease of the IMF ULF power at the ICME leading edge could be used to distinguish the ICME from the preceding sheath plasma.

Description: Key chemical NO x sink uncertainties and how they influence top-down emissions of nitrogen oxides Atmospheric Chemistry and Physics, 13, 9057-9082, 2013 Author(s): T. Stavrakou, J.-F. Müller, K. F. Boersma, R. J. van der A, J. Kurokawa, T. Ohara, and Q. Zhang Triggered by recent developments from laboratory and field studies regarding major NO x sink pathways in the troposphere, this study evaluates the influence of chemical uncertainties in NO x sinks for global NO x distributions calculated by the IMAGESv2 chemistry-transport model, and quantifies their significance for top-down NO x emission estimates. Our study focuses on five key chemical parameters believed to be of primary importance, more specifically, the rate of the reaction of NO 2 with OH radicals, the newly identified HNO 3 -forming channel in the reaction of NO with HO 2 , the reactive uptake of N 2 O 5 and HO 2 by aerosols, and the regeneration of OH in the oxidation of isoprene. Sensitivity simulations are performed to estimate the impact of each source of uncertainty. The model calculations show that, although the NO 2 +OH reaction is the largest NO x sink globally accounting for ca. 60% of the total sink, the reactions contributing the most to the overall uncertainty are the formation of HNO 3 in NO+HO 2 , leading to NO x column changes exceeding a factor of two over tropical regions, and the uptake of HO 2 by aqueous aerosols, in particular over East and South Asia. Emission inversion experiments are carried out using model settings which either minimise (MINLOSS) or maximise (MAXLOSS) the total NO x sink, both constrained by one year of OMI NO 2 column data from the DOMINO v2 KNMI algorithm. The choice of the model setup is found to have a major impact on the top-down flux estimates, with 75% higher emissions for MAXLOSS compared to the MINLOSS inversion globally. Even larger departures are found for soil NO (factor of 2) and lightning (1.8). The global anthropogenic source is better constrained (factor of 1.57) than the natural sources, except over South Asia where the combined uncertainty primarily associated to the NO+HO 2 reaction in summer and HO 2 uptake by aerosol in winter lead to top-down emission differences exceeding a factor of 2. Evaluation of the emission optimisation is performed against independent satellite observations from the SCIAMACHY sensor, with airborne NO 2 measurements of the INTEX-A and INTEX-B campaigns, as well as with two new bottom-up inventories of anthropogenic emissions in Asia (REASv2) and China (MEIC). Neither the MINLOSS nor the MAXLOSS setup succeeds in providing the best possible match with all independent datasets. Whereas the minimum sink assumption leads to better agreement with aircraft NO 2 profile measurements, consistent with the results of a previous analysis (Henderson et al., 2012), the same assumption leads to unrealistic features in the inferred distribution of emissions over China. Clearly, although our study addresses an important issue which was largely overlooked in previous inversion exercises, and demonstrates the strong influence of NO x loss uncertainties on top-down emission fluxes, additional processes need to be considered which could also influence the inferred source.

Description: Validation of stratospheric and mesospheric ozone observed by SMILES from International Space Station Atmospheric Measurement Techniques, 6, 2311-2338, 2013 Author(s): Y. Kasai, H. Sagawa, D. Kreyling, E. Dupuy, P. Baron, J. Mendrok, K. Suzuki, T. O. Sato, T. Nishibori, S. Mizobuchi, K. Kikuchi, T. Manabe, H. Ozeki, T. Sugita, M. Fujiwara, Y. Irimajiri, K. A. Walker, P. F. Bernath, C. Boone, G. Stiller, T. von Clarmann, J. Orphal, J. Urban, D. Murtagh, E. J. Llewellyn, D. Degenstein, A. E. Bourassa, N. D. Lloyd, L. Froidevaux, M. Birk, G. Wagner, F. Schreier, J. Xu, P. Vogt, T. Trautmann, and M. Yasui We observed ozone (O 3 ) in the vertical region between 250 and 0.0005 hPa (~ 12–96 km) using the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) on the Japanese Experiment Module (JEM) of the International Space Station (ISS) between 12 October 2009 and 21 April 2010. The new 4 K superconducting heterodyne receiver technology of SMILES allowed us to obtain a one order of magnitude better signal-to-noise ratio for the O 3 line observation compared to past spaceborne microwave instruments. The non-sun-synchronous orbit of the ISS allowed us to observe O 3 at various local times. We assessed the quality of the vertical profiles of O 3 in the 100–0.001 hPa (~ 16–90 km) region for the SMILES NICT Level 2 product version 2.1.5. The evaluation is based on four components: error analysis; internal comparisons of observations targeting three different instrumental setups for the same O 3 625.371 GHz transition; internal comparisons of two different retrieval algorithms; and external comparisons for various local times with ozonesonde, satellite and balloon observations (ENVISAT/MIPAS, SCISAT/ACE-FTS, Odin/OSIRIS, Odin/SMR, Aura/MLS, TELIS). SMILES O 3 data have an estimated absolute accuracy of better than 0.3 ppmv (3%) with a vertical resolution of 3–4 km over the 60 to 8 hPa range. The random error for a single measurement is better than the estimated systematic error, being less than 1, 2, and 7%, in the 40–1, 80–0.1, and 100–0.004 hPa pressure regions, respectively. SMILES O 3 abundance was 10–20% lower than all other satellite measurements at 8–0.1 hPa due to an error arising from uncertainties of the tangent point information and the gain calibration for the intensity of the spectrum. SMILES O 3 from observation frequency Band-B had better accuracy than that from Band-A. A two month period is required to accumulate measurements covering 24 h in local time of O 3 profile. However such a dataset can also contain variation due to dynamical, seasonal, and latitudinal effects.

Description: Stable isotopes dissect food webs from top to the bottom Biogeosciences Discussions, 10, 14923-14952, 2013 Author(s): J. J. Middelburg Stable isotopes have been used extensively to study food web functioning, i.e. the flow of energy and matter among organisms. Traditional food-web studies are based on the natural variability of carbon and nitrogen isotopes and are limited to larger organisms that can be physically separated from their environment. Recent developments allow isotope ratio measurements of microbes and this in turn allows then measurement of entire food webs, i.e. from small producers at the bottom to large consumers at the top. Here, I provide a concise review on the use and potential of stable isotope to reconstruct end-to-end food webs. I will first discuss food web reconstruction based on natural abundances isotope data and will then show that the use of stable isotopes as deliberately added tracers provides complementary information. Finally, challenges and opportunities for end-to-end food web reconstructions in a changing world are discussed.

Description: Study of the unknown HONO daytime source at an European suburban site during the MEGAPOLI summer and winter field campaigns Atmospheric Chemistry and Physics Discussions, 13, 23639-23690, 2013 Author(s): V. Michoud, A. Colomb, A. Borbon, K. Miet, M. Beekmann, M. Camredon, B. Aumont, S. Perrier, P. Zapf, G. Siour, W. Ait-Helal, C. Afif, A. Kukui, M. Furger, J. C. Dupont, M. Haeffelin, and J. F. Doussin Nitrous acid measurements were carried out during the MEGAPOLI summer and winter field campaigns at SIRTA observatory in Paris surroundings. Highly variable HONO levels were observed during the campaigns, ranging from 10 ppt to 500 ppt in summer and from 10 ppt to 1.7 ppb in winter. Significant HONO mixing ratios have also been measured during daytime hours, comprised between some tenth of ppt and 200 ppt for the summer campaign and between few ppt and 1 ppb for the winter campaign. Ancillary measurements, such as NO x , O 3 , photolysis frequencies, meteorological parameters (pressure, temperature, relative humidity, wind speed and wind direction), black carbon concentration, total aerosol surface area, boundary layer height and soil moisture, were conducted during both campaigns. In addition, for the summer period, OH radical measurements were made with a CIMS (Chemical Ionisation Mass Spectrometer). This large dataset has been used to investigate the HONO budget in a suburban environment. To do so, calculations of HONO concentrations using PhotoStationary State (PSS) approach have been performed, for daytime hours. The comparison of these calculations with measured HONO concentrations revealed an underestimation of the calculations making evident a missing source term for both campaigns. This unknown HONO source exhibits a bell shaped like average diurnal profile with a maximum around noon of approximately 0.7 ppb h −1 and 0.25 ppb h −1 , during summer and winter respectively. This source is the main HONO source during daytime hours for both campaigns. In both cases, this source shows a slight positive correlation with J (NO 2 ) and the product between J (NO 2 ) and soil moisture. This original approach had, thus, indicated that this missing source is photolytic and might be heterogeneous occurring on ground surface and involving water content available at the ground.

Description: Is blue intensity ready to replace maximum latewood density as a strong temperature proxy? A tree-ring case study on Scots pine from northern Sweden Climate of the Past Discussions, 9, 5227-5261, 2013 Author(s): J. A. Björklund, B. E. Gunnarson, K. Seftigen, J. Esper, and H. W. Linderholm At high latitudes, where low temperatures mainly limit tree-growth, measurements of wood density (e.g. Maximum Latewood Density, MXD) using the X-Ray methodology provide a temperature proxy that is superior to that of TRW. Density measurements are however costly and time consuming and have lead to experimentation with optical flatbed scanners to produce Maximum Blue Intensity (BI max ). BI max is an excellent proxy for density on annual scale but very limited in skill on centennial scale. Discolouration between samples is limiting BI max where specific brightnesses can have different densities. To overcome this, the new un-exploited parameter Δ blue intensity (ΔBI) was constructed by using the brightness in the earlywood (BI EW ) as background, (BI max − BI EW = ΔBI). This parameter was tested on X-Ray material (MXD − earlywood density = ΔMXD) and showed great potential both as a quality control and as a booster of climate signals. Unfortunately since the relationship between grey scale and density is not linear, and between-sample brightness can differ tremendously for similar densities, ΔBI cannot fully match ΔMXD in skill as climate proxy on centennial scale. For ΔBI to stand alone, the range of brightness/density offset must be reduced. Further studies are needed to evaluate this possibility, and solutions might include heavier sample treatment (reflux with chemicals) or image-data treatment (digitally manipulating base-line levels of brightness).

Description: An efficient method to generate a perturbed parameter ensemble of a fully coupled AOGCM without flux-adjustment Geoscientific Model Development, 6, 1447-1462, 2013 Author(s): P. J. Irvine, L. J. Gregoire, D. J. Lunt, and P. J. Valdes We present a simple method to generate a perturbed parameter ensemble (PPE) of a fully-coupled atmosphere-ocean general circulation model (AOGCM), HadCM3, without requiring flux-adjustment. The aim was to produce an ensemble that samples parametric uncertainty in some key variables and gives a plausible representation of the climate. Six atmospheric parameters, a sea-ice parameter and an ocean parameter were jointly perturbed within a reasonable range to generate an initial group of 200 members. To screen out implausible ensemble members, 20 yr pre-industrial control simulations were run and members whose temperature responses to the parameter perturbations were projected to be outside the range of 13.6 ± 2 °C, i.e. near to the observed pre-industrial global mean, were discarded. Twenty-one members, including the standard unperturbed model, were accepted, covering almost the entire span of the eight parameters, challenging the argument that without flux-adjustment parameter ranges would be unduly restricted. This ensemble was used in 2 experiments; an 800 yr pre-industrial and a 150 yr quadrupled CO 2 simulation. The behaviour of the PPE for the pre-industrial control compared well to ERA-40 reanalysis data and the CMIP3 ensemble for a number of surface and atmospheric column variables with the exception of a few members in the Tropics. However, we find that members of the PPE with low values of the entrainment rate coefficient show very large increases in upper tropospheric and stratospheric water vapour concentrations in response to elevated CO 2 and one member showed an implausible nonlinear climate response, and as such will be excluded from future experiments with this ensemble. The outcome of this study is a PPE of a fully-coupled AOGCM which samples parametric uncertainty and a simple methodology which would be applicable to other GCMs.

Description: Relationship between Amazon biomass burning aerosols and rainfall over La Plata Basin Atmospheric Chemistry and Physics Discussions, 13, 23995-24021, 2013 Author(s): G. Camponogara, M. A. F. Silva Dias, and G. G. Carrió High aerosol loads are discharged into the atmosphere by biomass burning in Amazon and Central Brazil during the dry season. These particles can interact with clouds as cloud condensation nuclei (CCN) changing cloud microphysics and radiative properties and, thereby, affecting the radiative budget of the region. Furthermore, the biomass burning aerosols can be transported by the low level jet (LLJ) to La Plata Basin where many mesoscale convective systems (MCS) are observed during spring and summer. This work proposes to investigate whether the aerosols from biomass burning may affect the MCS in terms of rainfall over La Plata Basin during spring. Since the aerosol effect is very difficult to isolate because convective clouds are very sensitive to small environment disturbances, detailed analyses using different techniques are used. The binplot, 2D histograms and combined empirical orthogonal function (EOF) methods are used to separate certain environment conditions with the possible effects of aerosol loading. Reanalysis 2, TRMM-3B42 and AERONET data are used from 1999 up to 2012 during September-December. The results show that there are two patterns associated to rainfall-aerosol interaction in La Plata Basin: one in which the dynamic conditions are more important than aerosols to generate rain; and a second one where the aerosol particles have a role in rain formation, acting mainly to suppress rainfall over La Plata Basin.

Description: Size-resolved aerosol composition and link to hygroscopicity at a forested site in Colorado Atmospheric Chemistry and Physics Discussions, 13, 23817-23843, 2013 Author(s): E. J. T. Levin, A. J. Prenni, B. Palm, D. Day, P. Campuzano-Jost, P. M. Winkler, S. M. Kreidenweis, P. J. DeMott, J. Jimenez, and J. N. Smith Aerosol hygroscopicity describes the ability of a particle to take up water and form a cloud droplet. Modeling studies have shown sensitivity of precipitation-producing cloud systems to the availability of aerosol particles capable of serving as cloud condensation nuclei (CCN), and hygroscopicity is a key parameter controlling the number of available CCN. Continental aerosol is typically assumed to have a representative hygroscopicity parameter, κ, of 0.3; however, in remote locations this value can be lower due to relatively large mass fractions of organic components. To further our understanding of aerosol properties in remote areas, we measured size-resolved aerosol chemical composition and hygroscopicity in a forested, mountainous site in Colorado during the six-week BEACHON-RoMBAS campaign. This campaign followed a year-long measurement period at this site, and results from the intensive campaign shed light on the previously reported seasonal cycle in aerosol hygroscopicity. New particle formation events were observed routinely at this site and nucleation mode composition measurements indicated that the newly formed particles were predominantly organic. These events likely contribute to the dominance of organic species at smaller sizes, where aerosol organic mass fractions of non-refractory components were between 70–90%. Corresponding aerosol hygroscopicity was observed to range from κ = 0.15–0.22, with hygroscopicity increasing with particle size. Aerosol chemical composition measured by an Aerosol Mass Spectrometer and calculated from hygroscopicity measurements agreed very well during the intensive study with an assumed value of κ org = 0.13 resulting in the best agreement.

Description: The impact of initial spread calibration on the RELO ensemble and its application to Lagrangian dynamics Nonlinear Processes in Geophysics, 20, 621-641, 2013 Author(s): M. Wei, G. Jacobs, C. Rowley, C. N. Barron, P. Hogan, P. Spence, O. M. Smedstad, P. Martin, P. Muscarella, and E. Coelho A number of real-time ocean model forecasts were carried out successfully at Naval Research Laboratory (NRL) to provide modeling support and numerical guidance to the CARTHE GLAD at-sea experiment during summer 2012. Two RELO ensembles and three single models using NCOM and HYCOM with different resolutions were carried out. A calibrated ensemble system with enhanced spread and reliability was developed to better support this experiment. The calibrated ensemble is found to outperform the un-calibrated ensemble in forecasting accuracy, skill, and reliability for all the variables and observation spaces evaluated. The metrics used in this paper include RMS error, anomaly correlation, PECA, Brier score, spread reliability, and Talagrand rank histogram. It is also found that even the un-calibrated ensemble outperforms the single forecast from the model with the same resolution. The advantages of the ensembles are further extended to the Lagrangian framework. In contrast to a single model forecast, the RELO ensemble provides not only the most likely Lagrangian trajectory for a particle in the ocean, but also an uncertainty estimate that directly reflects the complicated ocean dynamics, which is valuable for decision makers. The examples show that the calibrated ensemble with more reliability can capture trajectories in different, even opposite, directions, which would be missed by the un-calibrated ensemble. The ensembles are applied to compute the repelling and attracting Lagrangian coherent structures (LCSs), and the uncertainties of the LCSs, which are hard to obtain from a single model forecast, are estimated. It is found that the spatial scales of the LCSs depend on the model resolution. The model with the highest resolution produces the finest, small-scale, LCS structures, while the model with lowest resolution generates only large-scale LCSs. The repelling and attracting LCSs are found to intersect at many locations and create complex mesoscale eddies. The fluid particles and drifters in the middle of these tangles are subject to attraction and repulsion simultaneously from these two kinds of LCSs. As a result, the movements of particles near the Deepwater Horizon (DWH) location are severely limited. This is also confirmed by the Lagrangian trajectories predicted by the ensembles.

Description: Multifractal properties of embedded convective structures in orographic precipitation: toward subgrid-scale predictability Nonlinear Processes in Geophysics, 20, 605-620, 2013 Author(s): M. Nogueira, A. P. Barros, and P. M. A. Miranda Rain and cloud fields produced by fully nonlinear idealized cloud resolving numerical simulations of orographic convective precipitation display statistical multiscaling behavior, implying that multifractal diagnostics should provide a physically robust basis for the downscaling and sub-grid scale parameterizations of moist processes. Our results show that the horizontal scaling exponent function (and respective multiscaling parameters) of the simulated rainfall and cloud fields varies with atmospheric and terrain properties, particularly small-scale terrain spectra, atmospheric stability, and advective timescale. This implies that multifractal diagnostics of moist processes for these simulations are fundamentally transient, exhibiting complex nonlinear behavior depending on atmospheric conditions and terrain forcing at each location. A particularly robust behavior found here is the transition of the multifractal parameters between stable and unstable cases, which has a clear physical correspondence to the transition from stratiform to organized (banded and cellular) convective regime. This result is reinforced by a similar behavior in the horizontal spectral exponent. Finally, our results indicate that although nonlinearly coupled fields (such as rain and clouds) have different scaling exponent functions, there are robust relationships with physical underpinnings between the scaling parameters that can be explored for hybrid dynamical-statistical downscaling.

Description: Cluster observation of few-hour-scale evolution of structured plasma in the inner magnetosphere Annales Geophysicae, 31, 1569-1578, 2013 Author(s): M. Yamauchi, I. Dandouras, H. Rème, R. Lundin, and L. M. Kistler Using Cluster Ion Spectrometry (CIS) data from the spacecraft-4 perigee traversals during the 2001–2006 period (nearly 500 traversals after removing those that are highly contaminated by radiation belt particles), we statistically examined the local time distribution of structured trapped ions at sub- to few-keV range as well as inbound–outbound differences of these ion signatures in intensities and energy–latitude dispersion directions. Since the Cluster orbit during this period was almost constant and approximately north–south symmetric at nearly constant local time near the perigee, inbound–outbound differences are attributed to temporal developments in a 1–2 h timescale. Three types of structured ions at sub- to few keV range that are commonly found in the inner magnetosphere are examined: – Energy–latitude dispersed structured ions at less than a few keV, – Short-lived dispersionless ion stripes at wide energy range extending 0.1–10 keV, – Short-lived low-energy ion bursts at less than a few hundred eV. The statistics revealed that the wedge-like dispersed ions are most often observed in the dawn sector (60% of traversals), and a large portion of them show significant enhancement during the traversals at all local times. The short-lived ion stripes are predominantly found near midnight, where most stripes are significantly enhanced during the traversals and are associated with substorm activities with geomagnetic AL 〈 −300 nT. The low-energy bursts are observed at all local times and under all geomagnetic conditions, with moderate peak of the occurrence rate in the afternoon sector. A large portion of them again show significant enhancement or decay during the traversals.

Description: Atmospheric parameters in a subtropical cloud regime transition derived by AIRS+MODIS – observed statistical variability compared to ERA-Interim Atmospheric Chemistry and Physics Discussions, 13, 24051-24085, 2013 Author(s): M. M. Schreier, B. H. Kahn, K. Sušelj, J. Karlsson, S. C. Ou, Q. Yue, and S. L. Nasiri Cloud occurrence, microphysical and optical properties and atmospheric profiles within a subtropical cloud regime transition in the northeastern Pacific Ocean are obtained from a synergistic combination of the Atmospheric Infrared Sounder (AIRS) and the MODerate resolution Imaging Spectroradiometer (MODIS). The observed cloud parameters and atmospheric thermodynamic profile retrievals are binned by cloud type and analyzed based on their probability density functions (PDFs). Comparison of the PDFs to data from the European Center for Medium Range Weather Forecasting Re-analysis (ERA-Interim) shows a strong difference in the occurrence of the different cloud types compared to clear sky. An increasing non-Gaussian behavior is observed in cloud optical thickness (τ c ), effective radius ( r e ) and cloud top temperature ( T c ) distributions from Stratocumulus to Trade Cumulus, while decreasing values of lower tropospheric stability are seen. However, variations in the mean, width and shape of the distributions are found. The AIRS potential temperature (θ) and water vapor ( q ) profiles in the presence of varying marine boundary layer (MBL) cloud types show overall similarities to the ERA-Interim in the mean profiles, but differences arise in the higher moments at some altitudes. The differences between the PDFs from AIRS+MODIS and ERA-Interim make it possible to pinpoint systematic errors in both systems and helps to understand joint PDFs of cloud properties and coincident thermodynamic profiles from satellite observations.

Description: The Atmospheric radiation measurement (ARM) program network of microwave radiometers: instrumentation, data, and retrievals Atmospheric Measurement Techniques, 6, 2359-2372, 2013 Author(s): M. P. Cadeddu, J. C. Liljegren, and D. D. Turner The Climate Research Facility of the US Department of Energy's Atmospheric Radiation Measurement (ARM) Program operates a network of ground-based microwave radiometers. Data and retrievals from these instruments have been available to the scientific community for almost 20 yr. In the past five years the network has expanded to include a total of 22 microwave radiometers deployed in various locations around the world. The new instruments cover a frequency range between 22 and 197 GHz and are consistently and automatically calibrated. The latest addition to the network is a new generation of three-channel radiometers, currently in the early stage of deployment at all ARM sites. The network has been specifically designed to achieve increased accuracy in the retrieval of precipitable water vapor (PWV) and cloud liquid water path (LWP) with the long-term goal of providing the scientific community with reliable, calibrated radiometric data and retrievals of important geophysical quantities with well-characterized uncertainties. The radiometers provide high-quality, continuous datasets that can be utilized in a wealth of applications and scientific studies. This paper presents an overview of the microwave instrumentation, calibration procedures, data, and retrievals that are available for download from the ARM data archive.

Description: GAPPARD: a computationally efficient method of approximating gap-scale disturbance in vegetation models Geoscientific Model Development, 6, 1517-1542, 2013 Author(s): M. Scherstjanoi, J. O. Kaplan, E. Thürig, and H. Lischke Models of vegetation dynamics that are designed for application at spatial scales larger than individual forest gaps suffer from several limitations. Typically, either a population average approximation is used that results in unrealistic tree allometry and forest stand structure, or models have a high computational demand because they need to simulate both a series of age-based cohorts and a number of replicate patches to account for stochastic gap-scale disturbances. The detail required by the latter method increases the number of calculations by two to three orders of magnitude compared to the less realistic population average approach. In an effort to increase the efficiency of dynamic vegetation models without sacrificing realism, we developed a new method for simulating stand-replacing disturbances that is both accurate and faster than approaches that use replicate patches. The GAPPARD (approximating GAP model results with a Probabilistic Approach to account for stand Replacing Disturbances) method works by postprocessing the output of deterministic, undisturbed simulations of a cohort-based vegetation model by deriving the distribution of patch ages at any point in time on the basis of a disturbance probability. With this distribution, the expected value of any output variable can be calculated from the output values of the deterministic undisturbed run at the time corresponding to the patch age. To account for temporal changes in model forcing (e.g., as a result of climate change), GAPPARD performs a series of deterministic simulations and interpolates between the results in the postprocessing step. We integrated the GAPPARD method in the vegetation model LPJ-GUESS, and evaluated it in a series of simulations along an altitudinal transect of an inner-Alpine valley. We obtained results very similar to the output of the original LPJ-GUESS model that uses 100 replicate patches, but simulation time was reduced by approximately the factor 10. Our new method is therefore highly suited for rapidly approximating LPJ-GUESS results, and provides the opportunity for future studies over large spatial domains, allows easier parameterization of tree species, faster identification of areas of interesting simulation results, and comparisons with large-scale datasets and results of other forest models.

Description: Climate information based streamflow and rainfall forecasts for Huai River Basin using Hierarchical Bayesian Modeling Hydrology and Earth System Sciences Discussions, 10, 11559-11584, 2013 Author(s): X. Chen, Z. Hao, N. Devineni, and U. Lall A Hierarchal Bayesian model for forecasting regional summer rainfall and streamflow season-ahead using exogenous climate variables for East Central China is presented. The model provides estimates of the posterior forecasted probability distribution for 12 rainfall and 2 streamflow stations considering parameter uncertainty, and cross-site correlation. The model has a multilevel structure with regression coefficients modeled from a common multivariate normal distribution results in partial-pooling of information across multiple stations and better representation of parameter and posterior distribution uncertainty. Covariance structure of the residuals across stations is explicitly modeled. Model performance is tested under leave-10-out cross-validation. Frequentist and Bayesian performance metrics used include Receiver Operating Characteristic, Reduction of Error, Coefficient of Efficiency, Rank Probability Skill Scores, and coverage by posterior credible intervals. The ability of the model to reliably forecast regional summer rainfall and streamflow season-ahead offers potential for developing adaptive water risk management strategies.

Description: Halocline water modification and along slope advection at the Laptev Sea continental margin Ocean Science Discussions, 10, 1581-1617, 2013 Author(s): D. Bauch, S. Torres-Valdes, I. Polyakov, A. Novikhin, I. Dmitrenko, J. McKay, and A. Mix A general pattern in water mass distribution and potential shelf-basin exchanges is revealed at the Laptev Sea continental slope based on hydrochemical and stable oxygen isotope data from summers 2005–2009. Despite considerable interannual variations, a frontal system can be inferred between shelf, continental slope and central Eurasian Basin waters in the upper 100 m of the water column along the continental slope. Net sea-ice melt is consistently found at the continental slope; however the sea-ice meltwater signal is independent from the local retreat of the ice cover and appears to be advected from upwind locations. In addition to the along-slope frontal system at the continental shelf break a strong gradient is identified on the Laptev Sea shelf between 122 and 126° E with an eastward increase of riverine and sea-ice related brine water contents. These waters cross the shelf break at ~ 140° E and feed the Low Salinity Halocline Water (LSHW, salinity S 〈 33) in the upper 50 m of the water column. Extremely high silicate concentrations in Laptev Sea bottom waters may lead to speculation on a link to the local silicate maximum found within the salinity range of ~ 33 to 34.5, typical for the Lower Halocline Water (LHW) at the continental slope. But brine signatures and nutrient ratios from the central Laptev Sea differ from those at the continental slope. Thus a significant contribution of Laptev Sea bottom waters to the LHW at the continental slope can be excluded. The silicate maximum within the LHW at the continental slope may be formed locally or at the outer Laptev Sea shelf. Similar to the advection of the sea-ice melt signal along the Laptev Sea continental slope the nutrient signal at 50–70 m water depth within the LHW might also be fed by advection parallel to the slope. Thus, our analyses suggest that advective processes from upwind locations play a significant role in the halocline formation in the northern Laptev Sea.

Description: The effect of phase partitioning of semivolatile compounds on the measured CCN activity of aerosol particles Atmospheric Measurement Techniques Discussions, 6, 8413-8433, 2013 Author(s): S. Romakkaniemi, A. Jaatinen, A. Laaksonen, A. Nenes, and T. Raatikainen The effect of inorganic semivolatile aerosol compounds on the CCN activity of aerosol particles was studied by using a computational model for a DMT-CCN counter, a cloud parcel model for condensation kinetics and experiments to quantify the modelled results. Concentrations of water vapour and semivolatiles as well as aerosol trajectories in the CCN column were calculated by a computational fluid dynamics model. These trajectories and vapour concentrations were then used as an input for the cloud parcel model to simulate mass transfer kinetics of water and semivolatiles between aerosol particles and the gas phase. Two different questions were studied: (1) how big fraction of semivolatiles is evaporated from particles before activation in the CCN counter? (2) How much the CCN activity can be increased due to condensation of semivolatiles prior to the maximum water supersaturation in the case of high semivolatile concentration in the gas phase? The results show that, to increase the CCN activity of aerosol particles, a very high gas phase concentration (as compared to typical ambient conditions) is needed. We used nitric acid as a test compound. A concentration of several ppb or higher is needed for measurable effect. In the case of particle evaporation, we used ammonium nitrate as a test compound and found that it partially evaporates before maximum supersaturation is reached in the CCN counter, thus causing an underestimation of CCN activity. The effect of evaporation is clearly visible in all supersaturations, leading to an underestimation of the critical dry diameter by 10 to 15 nanometres in the case of ammonium nitrate particles in different supersaturations. This result was also confirmed by measurements in supersaturations between 0.1 and 0.7%.

Description: Responses of carbon dioxide flux and plant biomass to drought in a treed peatland in northern Alberta: a climate change perspective Biogeosciences Discussions, 10, 14999-15031, 2013 Author(s): T. M. Munir, B. Xu, M. Perkins, and M. Strack Northern peatland ecosystems represent large carbon stocks that are susceptible to changes such as accelerated mineralization due to water table lowering expected under a climate change scenario. During the growing seasons of 2011 and 2012 we monitored CO 2 fluxes and plant biomass along a microtopographic gradient (hummocks-hollows) in an undisturbed dry continental boreal treed bog (control) and a nearby site that was drained (drained) in 2001. Ten years of drainage in the bog significantly increased coverage of shrubs at hummocks and lichens at hollows. Considering measured hummock coverage and including tree incremental growth, we estimate that the control site was a larger sink in 2011 of −40 than that of −13 g C m −2 in 2012 while the drained site was a source of 144 and 140 g C m −2 over the same years. We infer that, drainage induced changes in vegetation growth led to increased biomass to counteract a portion of soil carbon losses. These results suggest that spatial variability (microtopography) and changes in vegetation community in boreal peatlands will affect how these ecosystems respond to lowered water table potentially induced by climate change.

Description: The coccolithophores Emiliania huxleyi and Coccolithus pelagicus : extant populations from the Norwegian-Iceland Sea and Fram Strait Biogeosciences Discussions, 10, 15077-15106, 2013 Author(s): C. V. Dylmer, J. Giraudeau, V. Hanquiez, and K. Husum Extant coccolithophores and their relation to the governing oceanographic features in the northern North Atlantic were investigated along two zonal transects of surface water sampling, both conducted during summer 2011 and fall 2007. The northern transects crossed Fram Strait and its two opposing boundary currents (West Spitsbergen Current and East Greenland Current), while the southern transects sampled the Norwegian and Iceland Seas (passing the island Jan Mayen) from the Lofoten Islands to the continental margin off Eastern Greenland. The distribution of the dominant coccolithophore species Emiliania huxleyi and Coccolithus pelagicus is discussed in view of both the surface hydrology at the time of sampling and the structure of the surface mixed layer. Remote-sensing images as well as CTD and ARGO profiles are used to constrain the physico-chemical state of the surface water at the time of sampling. Both transects were characterized by strong seasonal differences in bulk coccolithophore standing stocks with maximum values of 53 × 10 3 cells L −1 for the northern transect and 72 × 10 3 cells L −1 for the southern transect in fall and summer, respectively. The highest recorded bulk cell densities are essentially explained by E. huxleyi . This species shows a zonal shift in peak abundance in the Norwegian-Iceland Seas from a summer maximum in the Lofoten gyre to peak cell densities around the island Jan Mayen in fall. Vertical mixing of Atlantic waters west of Lofoten Island, a phenomenom related to pervasive summer large scale atmospheric changes in the eastern Nordic Seas, on one hand, and strengthened influence of melt-water and related surface water stratification around the island Jan Mayen during fall, on the other hand, explains the observed seasonal migration of the E. huxleyi peak production area, as well as the seasonal change in dominating species within the Iceland Sea. In addition our datasets are indicative of a well-defined maximum boundary temperature of 6 °C for the production of C. pelagicus in the northern North Atlantic. The Fram Strait transects provides, to our knowledge, a first view of the zonal distribution of extant coccolithophores in this remote setting during summer and fall. Our datasets are indicative of a seasonal change in the species community from an E. huxleyi -dominated assemblage during summer to a C. pelagicus -rich population during fall. Here, higher irradiance and increased Atlantic water influence during summer favored the production of the opportunistic species E. huxleyi close to the Arctic Front, whereas the peak production area during fall, with high concentrations of C. pelagicus , lays in true Arctic/Polar waters.

Description: AnaWEGE: a weather generator based on analogues of atmospheric circulation Geoscientific Model Development Discussions, 6, 4745-4774, 2013 Author(s): P. Yiou This paper presents a stochastic weather generator based on analogues of circulation (AnaWEGE). Analogues of circulation have been a promising paradigm to analyse climate variability and its extremes. The weather generator uses precomputed analogues of sea-level pressure over the North Atlantic. The stochastic rules of the generator constrain the continuity in time of the simulations. The generator then simulates spatially coherent time series of a climate variable, drawn from meteorological observations. The weather generator is tested for European temperatures, and for winter and summer seasons. The biases in temperature quantiles and autocorrelation are rather small compared to observed variability. The ability of simulating extremely hot summers and cold winters is also assessed.

Description: Methodological aspects of a pattern-scaling approach to produce global fields of monthly means of daily maximum and minimum temperature Geoscientific Model Development Discussions, 6, 4833-4882, 2013 Author(s): S. Kremser, G. E. Bodeker, and J. Lewis A Climate Pattern-Scaling Model (CPSM) that simulates global patterns of climate change, for a prescribed emissions scenario, is described. A CPSM works by quantitatively establishing the statistical relationship between a climate variable at a specific location (e.g. daily maximum surface temperature, T max ) and one or more predictor time series (e.g. global mean surface temperature, T global ) – referred to as the "training" of the CPSM. This training uses a regression model to derive fit-coefficients that describe the statistical relationship between the predictor time series and the target climate variable time series. Once that relationship has been determined, and given the predictor time series for any greenhouse gas (GHG) emissions scenario, the change in the climate variable of interest can be reconstructed – referred to as the "application" of the CPSM. The advantage of using a CPSM rather than a typical atmosphere-ocean global climate model (AOGCM) is that the predictor time series required by the CPSM can usually be generated quickly using a simple climate model (SCM) for any prescribed GHG emissions scenario and then applied to generate global fields of the climate variable of interest. The training can be performed either on historical measurements or on output from an AOGCM. Using model output from 21st century simulations has the advantage that the climate change signal is more pronounced than in historical data and therefore a more robust statistical relationship is obtained. The disadvantage of using AOGCM output is that the CPSM training might be compromised by any AOGCM inadequacies. For the purposes of exploring the various methodological aspects of the CPSM approach, AOGCM output was used in this study to train the CPSM. These investigations of the CPSM methodology focus on monthly mean fields of daily temperature extremes ( T max and T min ). Key conclusions are: (1) overall, the CPSM trained on simulations based on the Representative Concentration Pathway (RCP) 8.5 emissions scenario is able to reproduce AOGCM simulations of T max and T min based on predictor time series from an RCP 4.5 emissions scenario; (2) access to hemisphere average land and ocean temperatures as predictors improves the variance that can be explained, particularly over the oceans; (3) regression model fit-coefficients derived from individual simulations based on the RCP 2.6, 4.5 and 8.5 emissions scenarios agree well over most regions of the globe (the Arctic is the exception); (4) training the CPSM on concatenated time series from an ensemble of simulations does not result in fit-coefficients that explain significantly more of the variance than an approach that weights results based on single simulation fits; and (5) the inclusion of a linear time dependence in the regression model fit-coefficients improves the variance explained, primarily over the oceans.

Description: 3-D reflection seismic imaging of the Hontomín structure in the Basque-Cantabrian Basin (Spain) Solid Earth Discussions, 5, 1575-1614, 2013 Author(s): J. Alcalde, D. Martí, C. Juhlin, A. Malehmir, D. Sopher, E. Saura, I. Marzán, P. Ayarza, A. Calahorrano, A. Pérez-Estaún, and R. Carbonell The Basque-Cantabrian Basin of the Northern Iberia peninsula constitutes a unique example of a major deformation system, featuring a dome structure developed by extensional tectonics followed by compressional reactivation. The occurrence of natural resources in the area and the possibility of establishing a geological storage site for carbon dioxide motivated the acquisition of a 3-D seismic reflection survey in 2010, centered on the Jurassic Hontomín dome. The objectives of this survey were to obtain a geological model of the overall structure and to establish a baseline model for a possible geological CO 2 storage site. The 36 km 2 survey included approximately 5000 mixed (Vibroseis and explosives) source points recorded with a 25 m inline source and receiver spacing. The target reservoir is a saline aquifer, at approximately 1450 m depth, encased and sealed by carbonate formations. Acquisition and processing parameters were influenced by the rough topography and relatively complex geology. A strong near surface velocity inversion is evident in the data, affecting the quality of the data. The resulting 3-D image provides constraints on the key features of the geologic model. The Hontomín structure is interpreted to consist of an approximately 10 7 m 2 large elongated dome with two major W–E and NW–SE striking faults bounding it.

Description: Sea ice detection with space-based LIDAR The Cryosphere Discussions, 7, 4681-4701, 2013 Author(s): S. Rodier, Y. Hu, and M. Vaughan Monitoring long-term climate change in the Polar Regions relies on accurate, detailed and repeatable measurements of geophysical processes and states. These regions are among the Earth's most vulnerable ecosystems, and measurements there have shown rapid changes in the seasonality and the extent of snow and sea ice coverage. The authors have recently developed a promising new technique that uses lidar surface measurements from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission to infer ocean surface ice-water phase. CALIPSO's 532 nm depolarization ratio measurements of the ocean surface are uniquely capable of providing information about the ever-changing sea surface state within the Polar Regions. With the finer resolution of the CALIPSO footprint (90 m diameter, spaced 335 m apart) and its ability to acquire measurements during both daytime and nighttime orbit segments and in the presence of clouds, the CALIPSO sea ice product provides fine-scale information on mixed phase scenes and can be used to assess/validate the estimates of sea-ice concentration currently provided by passive sensors. This paper describes the fundamentals of the CALIPSO sea-ice detection and classification technique. We present retrieval results from a six-year study, which are compared to existing data sets obtained by satellite-based passive remote sensors.

Description: Airborne hydrogen cyanide measurements using a chemical ionisation mass spectrometer for the plume identification of biomass burning forest fires Atmospheric Chemistry and Physics, 13, 9217-9232, 2013 Author(s): M. Le Breton, A. Bacak, J. B. A. Muller, S. J. O'Shea, P. Xiao, M. N. R. Ashfold, M. C. Cooke, R. Batt, D. E. Shallcross, D. E. Oram, G. Forster, S. J.-B. Bauguitte, and C. J. Percival A chemical ionisation mass spectrometer (CIMS) was developed for measuring hydrogen cyanide (HCN) from biomass burning events in Canada using I − reagent ions on board the FAAM BAe-146 research aircraft during the BORTAS campaign in 2011. The ionisation scheme enabled highly sensitive measurements at 1 Hz frequency through biomass burning plumes in the troposphere. A strong correlation between the HCN, carbon monoxide (CO) and acetonitrile (CH 3 CN) was observed, indicating the potential of HCN as a biomass burning (BB) marker. A plume was defined as being 6 standard deviations above background for the flights. This method was compared with a number of alternative plume-defining techniques employing CO and CH 3 CN measurements. The 6-sigma technique produced the highest R 2 values for correlations with CO. A normalised excess mixing ratio (NEMR) of 3.68 ± 0.149 pptv ppbv −1 was calculated, which is within the range quoted in previous research (Hornbrook et al., 2011). The global tropospheric model STOCHEM-CRI incorporated both the observed ratio and extreme ratios derived from other studies to generate global emission totals of HCN via biomass burning. Using the ratio derived from this work, the emission total for HCN from BB was 0.92 Tg (N) yr −1 .

Description: Modeling of very low frequency (VLF) radio wave signal profile due to solar flares using the GEANT4 Monte Carlo simulation coupled with ionospheric chemistry Atmospheric Chemistry and Physics, 13, 9159-9168, 2013 Author(s): S. Palit, T. Basak, S. K. Mondal, S. Pal, and S. K. Chakrabarti X-ray photons emitted during solar flares cause ionization in the lower ionosphere (~60 to 100 km) in excess of what is expected to occur due to a quiet sun. Very low frequency (VLF) radio wave signals reflected from the D-region of the ionosphere are affected by this excess ionization. In this paper, we reproduce the deviation in VLF signal strength during solar flares by numerical modeling. We use GEANT4 Monte Carlo simulation code to compute the rate of ionization due to a M-class flare and a X-class flare. The output of the simulation is then used in a simplified ionospheric chemistry model to calculate the time variation of electron density at different altitudes in the D-region of the ionosphere. The resulting electron density variation profile is then self-consistently used in the LWPC code to obtain the time variation of the change in VLF signal. We did the modeling of the VLF signal along the NWC (Australia) to IERC/ICSP (India) propagation path and compared the results with observations. The agreement is found to be very satisfactory.

Description: Establishing the contribution of lawn mowing to atmospheric aerosol levels in American suburbs Atmospheric Chemistry and Physics Discussions, 13, 24435-24480, 2013 Author(s): R. M. Harvey, J. Zahardis, and G. A. Petrucci Green leaf volatiles (GLVs) are a class of wound-induced volatile organic compounds emitted by several plant species. Turfgrasses emit a complex profile of GLVs upon mowing, as evidenced by the "freshly cut grass" smell, some of which are readily oxidized in the atmosphere to contribute to secondary organic aerosol (SOA). The contribution of lawn mowing-induced SOA production may be especially impactful at the urban/suburban interface, where urban hubs provide a source of anthropogenic oxidants and SOA while suburban neighborhoods have the potential to emit large quantities of reactive, mow-induced GLVs. This interface provides a unique opportunity to study aerosol formation in a multi-component system and at a regionally relevant scale. Freshly cut grass was collected from a study site in Essex Junction, Vermont and was placed inside a 775 L Teflon experimental chamber. Thermal desorption gas chromatography mass spectrometry (TD-GC/MS) was used to characterize the emitted GLV profile. Ozone was introduced to the experimental chamber and TD-GC/MS was used to monitor the consumption of these GLVs and the subsequent evolution of gas phase products while a scanning mobility particle sizer was used to continuously measure aerosol size distributions and mass loadings as a result of grass clipping ozonolysis. Freshly cut grass found to emit a complex mixture of GLVs, dominated by cis -3-hexenyl acetate and cis -3-hexenol, which were released at an initial rate of 1.8 (±0.5) μg and 0.07 (±0.03) μg per square meter of lawn mowed with each mowing. Chamber studies using pure standards of cis -3-hexenyl acetate (CHA) and cis -3-hexenol (HXL) were found to have aerosol yields of 1.2 (±1.1)% and 3.3 (±3.1)%, respectively. Using these aerosol yields and the emission rate of these CHA and HXL by grass, SOA evolution by ozonolysis of grass clippings was predicted. However, the measured SOA mass produced from the ozonolysis of grass clippings exceeded the predicted amount, by upwards of ~ 150%. The ozonolysis of a mixture of CHA and HXL representative of environmental mixing ratios also failed to accurately model the SOA mass produced by grass clippings. Aerial photographs and geospatial analysis were used to determine the turfgrass coverage in a suburban neighborhood, which was then used along with measured SOA production as a function of grass mowed to determine that lawn mowing has the potential to contribute 47 μg m −2 SOA to the atmosphere per mowing event by ozonolysis, which cannot be modeled solely by the ozonolysis of CHA, HXL or a representative mixture of the two.

Description: Global distributions and trends of atmospheric ammonia (NH 3 ) from IASI satellite observations Atmospheric Chemistry and Physics Discussions, 13, 24301-24342, 2013 Author(s): M. Van Damme, L. Clarisse, C. L. Heald, D. Hurtmans, Y. Ngadi, C. Clerbaux, A. J. Dolman, J. W. Erisman, and P. F. Coheur Ammonia (NH 3 ) emissions in the atmosphere have strongly increased in the past decades, largely because of the intensive livestock production and use of fertilizers. As a short-lived species, NH 3 is highly variable in the atmosphere and its concentration is generally small, except in and close to local source areas. While ground-based measurements are possible, they are challenging and sparse. Advanced infrared sounders in orbit have recently demonstrated their capability to measure NH 3 , offering a new tool to refine global and regional budgets. In this paper we describe an improved retrieval scheme of NH 3 total columns from the measurements of the Infrared Atmospheric Sounding Interferometer (IASI). It exploits the hyperspectral character of this instrument by using an extended spectral range (800–1200 cm −1 ) where NH 3 is optically active. This scheme consists of the calculation of a dimensionless spectral index from the IASI level1C radiances, which is subsequently converted to a total NH 3 column using look-up-tables built from forward radiative transfer model simulations. We show how to retrieve the NH 3 total columns from IASI quasi-globally and twice daily, above both land and sea, without large computational resources and with an improved detection limit. The retrieval also provides error characterization on the retrieved columns. Five years of IASI measurements (1 November 2007 to 31 October 2012) have been processed to acquire the first global and multiple-year dataset of NH 3 total columns, which are evaluated and compared to similar products from other retrieval methods. Spatial distributions from the five years dataset are provided and analyzed at global and regional scales. We show in particular the ability of this method to identify smaller emission sources than those reported previously, as well as transport patterns above sea. The five year time series is further examined in terms of seasonality and inter-annual variability (in particular as a function of fire activity) separately for the Northern and Southern Hemispheres.

Description: Absorption properties of Mediterranean aerosols obtained from multi-year ground-based remote sensing observations Atmospheric Chemistry and Physics, 13, 9195-9210, 2013 Author(s): M. Mallet, O. Dubovik, P. Nabat, F. Dulac, R. Kahn, J. Sciare, D. Paronis, and J. F. Léon Aerosol absorption properties are of high importance to assess aerosol impact on regional climate. This study presents an analysis of aerosol absorption products obtained over the Mediterranean basin or land stations in the region from multi-year ground-based AERONET observations with a focus on the Absorbing Aerosol Optical Depth (AAOD), Single Scattering Albedo (SSA) and their spectral dependence. The AAOD and Absorption Angström Exponent (AAE) dataset is composed of daily averaged AERONET level 2 data from a total of 22 Mediterranean stations having long time series, mainly under the influence of urban-industrial aerosols and/or soil dust. This dataset covers the 17-yr period 1996–2012 with most data being from 2003–2011 (~89% of level-2 AAOD data). Since AERONET level-2 absorption products require a high aerosol load (AOD at 440 nm 〉 0.4), which is most often related to the presence of desert dust, we also consider level-1.5 SSA data, despite their higher uncertainty, and filter out data with an Angström exponent 〈 1.0 in order to study absorption by carbonaceous aerosols. The SSA dataset includes AERONET level-2 products. Sun-photometer observations show that values of AAOD at 440 nm vary between 0.024 ± 0.01 (resp. 0.040 ± 0.01) and 0.050 ± 0.01 (0.055 ± 0.01) for urban (dusty) sites. Analysis shows that the Mediterranean urban-industrial aerosols appear "moderately" absorbing with values of SSA close to ~0.94–0.95 ± 0.04 (at 440 nm) in most cases except over the large cities of Rome and Athens, where aerosol appears more absorbing (SSA ~0.89–0.90 ± 0.04). The aerosol Absorption Angström Exponent (AAE, estimated using 440 and 870 nm) is found to be larger than 1 for most sites over the Mediterranean, a manifestation of mineral dust (iron) and/or brown carbon producing the observed absorption. AERONET level-2 sun-photometer data indicate a possible East-West gradient, with higher values over the eastern basin (AAE East = 1.39/AAE West = 1.33). The North-South AAE gradient is more pronounced, especially over the western basin. Our additional analysis of AERONET level-1.5 data also shows that organic absorbing aerosols significantly affect some Mediterranean sites. These results indicate that current climate models treating organics as nonabsorbing over the Mediterranean certainly underestimate the warming effect due to carbonaceous aerosols.

Description: Airborne sun photometer PLASMA: concept, measurements, comparison of aerosol extinction vertical profile with lidar Atmospheric Measurement Techniques, 6, 2383-2389, 2013 Author(s): Y. Karol, D. Tanré, P. Goloub, C. Vervaerde, J. Y. Balois, L. Blarel, T. Podvin, A. Mortier, and A. Chaikovsky A 15-channel airborne sun-tracking photometer has been developed. The instrument provides aerosol optical depths over a wide spectral range (0.34–2.25 μm) with an accuracy (ΔAOD) of approximately 0.01. Taking measurements at different altitudes allow us to derive the aerosol extinction vertical profile. Thanks to the wide spectral range of the instrument, information on the aerosol size distribution along the vertical is also available.

Description: Modeling different freeze/thaw processes in heterogeneous landscapes of the Arctic polygonal tundra using an ecosystem model Geoscientific Model Development Discussions, 6, 4883-4932, 2013 Author(s): given_name prefix surname suffix, S. Yi, K. Wischnewski, M. Langer, S. Muster, and J. Boike Freeze/thaw (F/T) processes can be quite different under the various land surface types found in the heterogeneous polygonal tundra of the Arctic. Proper simulation of these different processes is essential for accurate prediction of the release of greenhouse gases under a warming climate scenario. In this study we have modified the dynamic organic soil version of the Terrestrial Ecosystem Model (DOS-TEM) to simulate F/T processes beneath the polygon rims, polygon centers (with and without water), and lakes that are common features in Arctic lowland regions. We first verified the F/T algorithm in the DOS-TEM against analytical solutions, and then compared the results with in situ measurements from Samoylov Island, Siberia. In the final stage, we examined the different responses of the F/T processes for different water levels at the various land surface types. The simulations revealed that (1) the DOS-TEM was very efficient and its results compared very well with analytical solutions for idealized cases, (2) the simulations compared reasonably well with in situ measurements although there were a number of model limitations and uncertainties, (3) the DOS-TEM was able to successfully simulate the differences in F/T dynamics under different land surface types, and (4) permafrost beneath water bodies was found to respond highly sensitive to changes in water depths between 1 and 2 m. Our results indicate that water is very important in the thermal processes simulated by the DOS-TEM; the heterogeneous nature of the landscape and different water depths therefore need to be taken into account when simulating methane emission responses to a warming climate.

Description: Anthropogenic and natural methane fluxes in Switzerland synthesized within a spatially-explicit inventory Biogeosciences Discussions, 10, 15181-15224, 2013 Author(s): R. V. Hiller, D. Bretscher, T. DelSontro, T. Diem, W. Eugster, R. Henneberger, S. Hobi, E. Hodson, D. Imer, M. Kreuzer, T. Künzle, L. Merbold, P. A. Niklaus, B. Rihm, A. Schellenberger, M. H. Schroth, C. J. Schubert, H. Siegrist, J. Stieger, N. Buchmann, and D. Brunner We present the first high-resolution (500 m × 500 m) gridded methane (CH 4 ) emission inventory for Switzerland, which integrates the national emission totals reported to the United Nations Framework Convention on Climate Change (UNFCCC) and recent CH 4 flux studies conducted by research groups across Switzerland. In addition to anthropogenic emissions, we also include natural and semi-natural CH 4 fluxes, i.e., emissions from lakes and reservoirs, wetlands, wild animals as well as uptake by forest soils. National CH 4 emissions were disaggregated using detailed geostatistical information on source locations and their spatial extent and process- or area-specific emission factors. In Switzerland, the highest CH 4 emissions in 2011 originated from the agricultural sector (150 Gg CH 4 yr −1 ), mainly produced by ruminants and manure management, followed by emissions from waste management (15 Gg CH 4 yr −1 ) mainly from landfills and the energy sector (12 Gg CH 4 yr −1 ), which was dominated by emissions from natural gas distribution. Compared to the anthropogenic sources, emissions from natural and semi-natural sources were relatively small (6 Gg CH 4 yr −1 ), making up only 3 % of the total emissions in Switzerland. CH 4 fluxes from agricultural soils were estimated to be not significantly different from zero (between −1.5 and 0 Gg CH 4 yr −1 ), while forest soils are a CH 4 sink (approx. −2.8 Gg CH 4 yr −1 ), partially offsetting other natural emissions. Estimates of uncertainties are provided for the different sources, including an estimate of spatial disaggregation errors deduced from a comparison with a global (EDGAR v4.2) and a European CH 4 inventory (TNO/MACC). This new spatially-explicit emission inventory for Switzerland will provide valuable input for regional scale atmospheric modeling and inverse source estimation.

Description: Seasonal trends of dry and bulk concentration of nitrogen compounds over a rain forest in Ghana Biogeosciences Discussions, 10, 15225-15255, 2013 Author(s): F. Fattore, T. Bertolini, S. Materia, S. Gualdi, A. Thongo M'Bou, G. Nicolini, R. Valentini, A. De Grandcourt, D. Tedesco, and S. Castaldi African tropical forests of the equatorial belt might receive significant input of extra nitrogen derived from biomass burning occurring in the north savanna belt and transported equator wards by NE winds. In order to test this hypothesis an experiment was set up in a tropical rain forest in the National park of Ankasa (Ghana) aiming at: quantifying magnitude and seasonal variability of concentrations of N compounds, present as gas and aerosol (dry nitrogen) or in the rainfall (bulk nitrogen), over the studied forest; relating their seasonal variability to trends of local and regional winds and rainfall and to variations of fire events in the region. Three Delta systems, implemented for monthly measurements of NO 2 , were mounted over a tower at 45 m height, 20 m above forest canopy to sample gas (NH 3 , NO 2 , HNO 3 , HCl, SO 2 ) and aerosol (NH 4 + , NO 3 − , and several ions), together with three tanks for bulk rainfall collection (to analyze NH 4 + , NO 3 − and ion concentration). The tower was provided with a sonic anemometer to estimate local wind data. The experiment started in October 2011 and data up to October 2012 are presented. To interpret the observed seasonal trends of measured compounds, local and regional meteo data and regional satellite fire data were analyzed. The concentration of N compounds significantly increased from December to April, during the drier period, peaking in December-February when North Eastern winds (Harmattan) were moving dry air masses over the West central African region and the inter tropical convergence zone (ITCZ) was at its minimum latitude over the equator. This period also coincided with peaks of fire in the whole region. On the contrary, N concentration in gas, aerosol and rain decreased from May to October when prevalent winds arrived from the sea (South-East), during the Monsoon period. Both ionic compositions of rain and analysis of local wind direction showed a significant and continuous presence of see-breeze at site. The ionic composition of rain water resulted much closer to see water and poorer in N compounds from May to October.

Description: Aircraft measurements of polar organic tracer compounds in tropospheric particles (PM 10 ) over Central China Atmospheric Chemistry and Physics Discussions, 13, 24481-24516, 2013 Author(s): P. Q. Fu, K. Kawamura, Y. F. Cheng, S. Hatakeyama, A. Takami, H. Li, and W. Wang Atmospheric aerosol samples were collected by aircraft at low to middle altitudes (0.8–3.5 km a.g.l.) over Central East to West China during summer 2003 and spring 2004. The samples were analyzed for polar organic compounds using a technique of solvent extraction/BSTFA derivatization/gas chromatography-mass spectrometry. Biogenic secondary organic aerosol (SOA) tracers from the oxidation of isoprene were found to be more abundant in summer (3.3–138 ng m −3 , mean 39 ng m −3 ) than in spring (3.2–42 ng m −3 , 15 ng m −3 ), while α/β-pinene and β-caryophyllene SOA tracers showed similar abundance between these two seasons. A strong positive correlation ( R 2 =0.83) between levoglucosan and β-caryophyllinic acid was found in the spring samples versus a weak correlation ( R 2 =0.17) in the summer samples, implying substantial contributions from biomass burning to the β-caryophyllinic acid production in spring. Two organic nitrogen species (oxamic acid and carbamide) were detected in the aircraft aerosol samples and their concentrations were comparable to those of biogenic SOA tracers. Most of the POA and SOA tracers were less abundant at higher altitudes, suggesting they are of ground surface origin, either being directly emitted from anthropogenic/natural sources on the ground surface, or rapidly formed through photooxidation of their precursors emitted from the ground surface and then diluted during uplifting into the troposphere. This study demonstrates that primary biological aerosols, biogenic SOA, and organic nitrogen species are important components of organic aerosols in the troposphere over Central China.

Description: McClear: a new model estimating downwelling solar radiation at ground level in clear-sky conditions Atmospheric Measurement Techniques, 6, 2403-2418, 2013 Author(s): M. Lefèvre, A. Oumbe, P. Blanc, B. Espinar, B. Gschwind, Z. Qu, L. Wald, M. Schroedter-Homscheidt, C. Hoyer-Klick, A. Arola, A. Benedetti, J. W. Kaiser, and J.-J. Morcrette A new fast clear-sky model called McClear was developed to estimate the downwelling shortwave direct and global irradiances received at ground level under clear skies. It is a fully physical model replacing empirical relations or simpler models used before. It exploits the recent results on aerosol properties, and total column content in water vapour and ozone produced by the MACC project (Monitoring Atmosphere Composition and Climate). It accurately reproduces the irradiance computed by the libRadtran reference radiative transfer model with a computational speed approximately 10 5 times greater by adopting the abaci, or look-up table, approach combined with interpolation functions. It is therefore suited for geostationary satellite retrievals or numerical weather prediction schemes with many pixels or grid points, respectively. McClear irradiances were compared to 1 min measurements made in clear-sky conditions at several stations within the Baseline Surface Radiation Network in various climates. The bias for global irradiance comprises between −6 and 25 W m −2 . The RMSE ranges from 20 W m −2 (3% of the mean observed irradiance) to 36 W m −2 (5%) and the correlation coefficient ranges between 0.95 and 0.99. The bias for the direct irradiance comprises between −48 and +33 W m −2 . The root mean square error (RMSE) ranges from 33 W m −2 (5%) to 64 W m −2 (10%). The correlation coefficient ranges between 0.84 and 0.98. This work demonstrates the quality of the McClear model combined with MACC products, and indirectly the quality of the aerosol properties modelled by the MACC reanalysis.

Description: Optical, microphysical, mass and geometrical properties of aged volcanic particles observed over Athens, Greece, during the Eyjafjallajökull eruption in April 2010 through synergy of Raman lidar and sunphotometer measurements Atmospheric Chemistry and Physics, 13, 9303-9320, 2013 Author(s): P. Kokkalis, A. Papayannis, V. Amiridis, R. E. Mamouri, I. Veselovskii, A. Kolgotin, G. Tsaknakis, N. I. Kristiansen, A. Stohl, and L. Mona Vertical profiles of the optical (extinction and backscatter coefficients, lidar ratio and Ångström exponent), microphysical (mean effective radius, mean refractive index, mean number concentration) and geometrical properties as well as the mass concentration of volcanic particles from the Eyjafjallajökull eruption were retrieved at selected heights over Athens, Greece, using multi-wavelength Raman lidar measurements performed during the period 21–24 April 2010. Aerosol Robotic Network (AERONET) particulate columnar measurements along with inversion schemes were initialized together with lidar observations to deliver the aforementioned products. The well-known FLEXPART (FLEXible PARTicle dispersion model) model used for volcanic dispersion simulations is initiated as well in order to estimate the horizontal and vertical distribution of volcanic particles. Compared with the lidar measurements within the planetary boundary layer over Athens, FLEXPART proved to be a useful tool for determining the state of mixing of ash with other, locally emitted aerosol types. The major findings presented in our work concern the identification of volcanic particles layers in the form of filaments after 7-day transport from the volcanic source (approximately 4000 km away from our site) from the surface and up to 10 km according to the lidar measurements. Mean hourly averaged lidar signals indicated that the layer thickness of volcanic particles ranged between 1.5 and 2.2 km. The corresponding aerosol optical depth was found to vary from 0.01 to 0.18 at 355 nm and from 0.02 up to 0.17 at 532 nm. Furthermore, the corresponding lidar ratios ( S ) ranged between 60 and 80 sr at 355 nm and 44 and 88 sr at 532 nm. The mean effective radius of the volcanic particles estimated by applying inversion scheme to the lidar data found to vary within the range 0.13–0.38 μm and the refractive index ranged from 1.39+0.009 i to 1.48+0.006 i . This high variability is most probably attributed to the mixing of aged volcanic particles with other aerosol types of local origin. Finally, the LIRIC (LIdar/Radiometer Inversion Code) lidar/sunphotometric combined inversion algorithm has been applied in order to retrieve particle concentrations. These have been compared with FLEXPART simulations of the vertical distribution of ash showing good agreement concerning not only the geometrical properties of the volcanic particles layers but also the particles mass concentration.

Description: Review of some significant claimed irregularities in Scandinavian postglacial uplift in time scales from tens to thousands of years: earthquakes? Solid Earth Discussions, 5, 1615-1640, 2013 Author(s): S. Gregersen The postglacial uplift/subsidence in Scandinavia is regular. And the phenomenon is similar in time scales of tens, hundreds and thousands of years studied via geodesy, seismology and geology. Searches for irregularities in the form of earthquakes claimed in the scientific literature have disclosed many earthquakes right after the Ice Age and some later cases for further evaluation. In a previous report the present author has mentioned doubts about the validity of some of the most significant claimed irregularities. In the present paper a review is made of these significant claimed irregularities in the southwestern flank of the Scandinavian postglacial uplift/subsidence via literature studies of geodetic and geological claims of earthquakes as well as discussions in the field. Geodetic observations exist for all of Scandinavia. Those describe the phenomenon in 10s–100s of years scale. Earthquake observations in seismology are of relevance in the same time scales. Geological studies of dated shore lines describe the postglacial vertical earth-surface motion in a quite different time scale of 100s–1000s of years. There is a need for integration of these observations geographically. This is happening in the various time scales in the DynaQlim project. The review finds the claims improbable about the following: (1) geodynamical motion in the Copenhagen area, (2) a paleo-earthquake in Læsø and (3) the recently proposed water level discrepancy in the southern part of Denmark. The assessment is less certain, but falls to improbable concerning (4) proposed paleo-earthquakes by Hallandsåsen in southwestern Sweden.

Description: Using the Firefly optimization method to weight an ensemble of rainfall forecasts from the Brazilian developments on the Regional Atmospheric Modeling System (BRAMS) Advances in Geosciences, 35, 123-136, 2013 Author(s): A. F. dos Santos, S. R. Freitas, J. G. Z. de Mattos, H. F. de Campos Velho, M. A. Gan, E. F. P. da Luz, and G. A. Grell In this paper we consider an optimization problem applying the metaheuristic Firefly algorithm (FY) to weight an ensemble of rainfall forecasts from daily precipitation simulations with the Brazilian developments on the Regional Atmospheric Modeling System (BRAMS) over South America during January 2006. The method is addressed as a parameter estimation problem to weight the ensemble of precipitation forecasts carried out using different options of the convective parameterization scheme. Ensemble simulations were performed using different choices of closures, representing different formulations of dynamic control (the modulation of convection by the environment) in a deep convection scheme. The optimization problem is solved as an inverse problem of parameter estimation. The application and validation of the methodology is carried out using daily precipitation fields, defined over South America and obtained by merging remote sensing estimations with rain gauge observations. The quadratic difference between the model and observed data was used as the objective function to determine the best combination of the ensemble members to reproduce the observations. To reduce the model rainfall biases, the set of weights determined by the algorithm is used to weight members of an ensemble of model simulations in order to compute a new precipitation field that represents the observed precipitation as closely as possible. The validation of the methodology is carried out using classical statistical scores. The algorithm has produced the best combination of the weights, resulting in a new precipitation field closest to the observations.

Description: Instrument self-noise and sensor misalignment Advances in Geosciences, 36, 17-20, 2013 Author(s): A. Gerner and G. Bokelmann In this study we investigate self-noise of RefTek TM 151-60A "Observer" broadband seismometers (flat to velocity between 50 Hz down to T 0 =60 s, f 0 ≈17 mHz) using the coherency analysis method introduced by Sleeman et al. (2006). We present a self-noise model for this type of sensor and compare it to the self-noise models of the standard observatory sensor STS-2 (Streckeisen) and RefTek's 151-120 broadband seismometer, which both have natural periods T 0 of 120 s. We further report on the sensitivity of this technique to sensor misalignment and our success of eliminating leakage of the omnipresent microseism noise into self-noise estimates by numerically rotating seismic traces in order to find real self-noise.

Description: Review Article: Potential geomorphic consequences of a future great ( M w = 8.0+) Alpine Fault earthquake, South Island, New Zealand Natural Hazards and Earth System Science, 13, 2279-2299, 2013 Author(s): T. R. Robinson and T. R. H. Davies The Alpine Fault in New Zealand's South Island has not sustained a large magnitude earthquake since ca. AD 1717. The time since this rupture is close to the average inferred recurrence interval of the fault (~300 yr). The Alpine Fault is therefore expected to generate a large magnitude earthquake in the near future. Previous ruptures of this fault are inferred to have generated M w = 8.0 or greater earthquakes and to have resulted in, amongst other geomorphic hazards, large-scale landslides and landslide dams throughout the Southern Alps. There is currently 85% probability that the Alpine Fault will cause a M w = 8.0+ earthquake within the next 100 yr. While the seismic hazard is fairly well understood, that of the consequential geomorphic activity is less well studied, and these consequences are explored herein. They are expected to include landsliding, landslide damming, dam-break flooding, debris flows, river aggradation, liquefaction, and landslide-generated lake/fiord tsunami. Using evidence from previous events within New Zealand as well as analogous international examples, we develop first-order estimates of the likely magnitude and possible locations of the geomorphic effects associated with earthquakes. Landsliding is expected to affect an area 〉 30 000 km 2 and involve 〉 1billion m 3 of material. Some tens of landslide dams are expected to occur in narrow, steep-sided gorges in the affected region. Debris flows will be generated in the first long-duration rainfall after the earthquake and will continue to occur for several years as rainfall (re)mobilises landslide material. In total more than 1000 debris flows are likely to be generated at some time after the earthquake. Aggradation of up to 3 m will cover an area 〉 125 km 2 and is likely to occur on many West Coast alluvial fans and floodplains. The impact of these effects will be felt across the entire South Island and is likely to continue for several decades.

Description: Corrigendum to "The magnitude and causes of uncertainty in global model simulations of cloud condensation nuclei" published in Atmos. Chem. Phys., 13, 8879–8914, 2013 Atmospheric Chemistry and Physics, 13, 9375-9377, 2013 Author(s): L. A. Lee, K. J. Pringle, C. L. Reddington, G. W. Mann, P. Stier, D. V. Spracklen, J. R. Pierce, and K. S. Carslaw No abstract available.

Description: Black carbon emissions from in-use ships: a California regional assessment Atmospheric Chemistry and Physics Discussions, 13, 24675-24712, 2013 Author(s): G. M. Buffaloe, D. A. Lack, E. J. Williams, D. Coffman, K. L. Hayden, B. M. Lerner, S-M. Li, I. Nuaaman, P. Massoli, T. B. Onasch, P. K. Quinn, and C. D. Cappa Black carbon (BC) mass emission factors (EF BC ; g-BC (kg-fuel) −1 ) from a variety of ocean going vessels have been determined from measurements of BC and carbon dioxide (CO 2 ) concentrations in ship plumes intercepted by the R/V Atlantis during the 2010 California Nexus (CalNex) campaign. The ships encountered were all operating within 24 nautical miles of the California coast and were utilizing relatively low sulphur fuels. Black carbon concentrations within the plumes, from which EF BC values are determined, were measured using four independent instruments: a photoacoustic spectrometer and a particle soot absorption photometer, which measure light absorption, and a single particle soot photometer and soot particle aerosol mass spectrometer, which measure the mass concentration of refractory BC directly. The measured EF BC have been divided into vessel type categories and engine type categories, from which averages have been determined. The geometric average EF BC , determined from over 71 vessels and 135 plumes encountered, was 0.31 g-BC (kg-fuel) −1 . The most frequent engine type encountered was the slow speed diesel (SSD), and the most frequent SSD vessel type was the cargo ship sub-category. Average and median EF BC values from the SSD category are compared with previous observations from the Texas Air Quality Study (TexAQS) in 2006, during which the ships encountered were predominately operating on high sulphur fuels. There is a statistically significant difference between the EF BC values from CalNex and TexAQS for SSD vessels and for the cargo and tanker ship types within this engine category. The CalNex EF BC values are lower than those from TexAQS, suggesting that operation on lower sulphur fuels is associated with smaller EF BC values.

Description: A case study into the measurement of ship emissions from plume intercepts of the NOAA Ship Miller Freeman Atmospheric Chemistry and Physics Discussions, 13, 24635-24674, 2013 Author(s): C. D. Cappa, E. J. Williams, D. A. Lack, G. M. Buffaloe, D. Coffman, K. L. Hayden, S. C. Herndon, B. M. Lerner, S-M. Li, P. Massoli, R. McLaren, I. Nuaaman, T. B. Onasch, and P. K. Quinn Emissions factors (EFs) for gas and sub-micron particle-phase species were measured in intercepted plumes as a function of vessel speed from an underway research vessel, the NOAA Ship Miller Freeman , operating a medium-speed diesel engine on low-sulfur marine gas oil. For many of the particle-phase species, EFs were determined using multiple measurement methodologies, allowing for an assessment of how well EFs from different techniques agree. The total sub-micron PM (PM 1 ) was dominated by particulate black carbon (BC) and particulate organic matter (POM), with an average POM / BC ratio of 1.3. Consideration of the POM / BC ratios observed here with literature studies suggests that laboratory and in-stack measurement methods may over-estimate primary POM EFs relative to those observed in emitted plumes. Comparison of four different methods for black carbon measurement indicates that careful attention must be paid to instrument limitations and biases when assessing EF BC . Particulate sulfate (SO 4 2− ) EFs were extremely small and the particles emitted by Miller Freeman were inefficient as cloud condensation nuclei (CCN), even at high super saturations, consistent with the use of very low sulfur fuel and the overall small emitted particle sizes. All measurement methodologies consistently demonstrate that the measured EFs (fuel mass basis) for PM 1 mass, BC and POM decreased as the ship slowed. Particle number EFs were approximately constant across the speed change, with a shift towards smaller particles being emitted at slower speeds. Emissions factors for gas-phase CO and formaldehyde (HCHO) both increased as the vessel slowed, while EFs for NO x decreased and SO 2 EFs were approximately constant.

Description: Ten years of MIPAS measurements with ESA Level 2 processor V6 – Part 1: Retrieval algorithm and diagnostics of the products Atmospheric Measurement Techniques, 6, 2419-2439, 2013 Author(s): P. Raspollini, B. Carli, M. Carlotti, S. Ceccherini, A. Dehn, B. M. Dinelli, A. Dudhia, J.-M. Flaud, M. López-Puertas, F. Niro, J. J. Remedios, M. Ridolfi, H. Sembhi, L. Sgheri, and T. von Clarmann The MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) instrument on the Envisat (Environmental satellite) satellite has provided vertical profiles of the atmospheric composition on a global scale for almost ten years. The MIPAS mission is divided in two phases: the full resolution phase, from 2002 to 2004, and the optimized resolution phase, from 2005 to 2012, which is characterized by a finer vertical and horizontal sampling attained through a reduction of the spectral resolution. While the description and characterization of the products of the ESA processor for the full resolution phase has been already described in previous papers, in this paper we focus on the performances of the latest version of the ESA (European Space Agency) processor, named ML2PP V6 (MIPAS Level 2 Prototype Processor), which has been used for reprocessing the entire mission. The ESA processor had to perform the operational near real time analysis of the observations and its products needed to be available for data assimilation. Therefore, it has been designed for fast, continuous and automated analysis of observations made in quite different atmospheric conditions and for a minimum use of external constraints in order to avoid biases in the products. The dense vertical sampling of the measurements adopted in the second phase of the MIPAS mission resulted in sampling intervals finer than the instantaneous field of view of the instrument. Together with the choice of a retrieval grid aligned with the vertical sampling of the measurements, this made ill-conditioned the retrieval problem of the MIPAS operational processor. This problem has been handled with minimal changes to the original retrieval approach but with significant improvements nonetheless. The Levenberg–Marquardt method, already present in the retrieval scheme for its capability to provide fast convergence for nonlinear problems, is now also exploited for the reduction of the ill-conditioning of the inversion. An expression specifically designed for the regularizing Levenberg–Marquardt method has been implemented for the computation of the covariance matrices and averaging kernels of the retrieved products. The regularization of the Levenberg–Marquardt method is controlled by the convergence criteria and is deliberately kept weak. The resulting oscillations of the retrieved profile are a posteriori damped by an innovative self-adapting Tikhonov regularization. The convergence criteria and the weakness of the self-adapting regularization ensure that minimum constraints are used and the best vertical resolution obtainable from the measurements is achieved in all atmospheric conditions. Random and systematic errors, as well as vertical and horizontal resolution are compared in the two phases of the mission for all products, namely: temperature, H 2O , O 3 , HNO 3 , CH 4 , N 2O , NO 2 , CFC-11, CFC-12, N 2 O 5 and ClONO 2 . The use in the two phases of the mission of different optimized sets of spectral intervals ensures that, despite the different spectral resolutions, comparable performances are obtained in the whole MIPAS mission in terms of random and systematic errors, while the vertical resolution and the horizontal resolution are significantly better in the case of the optimized resolution measurements.

Description: Ground-based stratospheric O 3 and HNO 3 measurements at Thule, Greenland: an intercomparison with Aura MLS observations Atmospheric Measurement Techniques, 6, 2441-2453, 2013 Author(s): I. Fiorucci, G. Muscari, L. Froidevaux, and M. L. Santee In response to the need for improving our understanding of the evolution and the interannual variability of the winter Arctic stratosphere, in January 2009 a Ground-Based Millimeter-wave Spectrometer (GBMS) was installed at the Network for the Detection of Atmospheric Composition Change (NDACC) site in Thule (76.5° N, 68.8° W), Greenland. In this work, stratospheric GBMS O 3 and HNO 3 vertical profiles obtained from Thule during the winters 2010 (HNO 3 only), 2011 and 2012 are characterized and intercompared with co-located measurements of the Aura Microwave Limb Sounder (MLS) experiment. Using a recently developed algorithm based on Optimal Estimation, we find that the GBMS O 3 retrievals show good sensitivity (〉 80%) to atmospheric variations between ~ 17 and ~ 50 km, where their 1σ uncertainty is estimated to be the larger of ~ 11% or 0.2 ppmv. Similarly, HNO 3 profiles can be considered for scientific use between ~ 17 and ~ 45 km altitude, with a 1σ uncertainty that amounts to the larger of 15% or 0.2 ppbv. Comparisons with Aura MLS version 3.3 observations show that, on average, GBMS O 3 mixing ratios are biased negatively with respect to MLS throughout the stratosphere, with differences ranging between ~ 0.3 ppmv (8%) and 0.9 ppmv (18%) in the 17–50 km vertical range. GBMS HNO 3 values display instead a positive bias with respect to MLS up to 26 km, reaching a maximum of ~ 1 ppbv (10%) near the mixing ratio profile peak. O 3 and HNO 3 values from the two datasets prove to be well correlated at all altitudes, although their correlations worsen at the lower end of the altitude ranges considered. Column contents of GBMS and MLS O 3 (from 20 km upwards) and HNO 3 (from 17 km upwards) correlate very well and indicate that GBMS measurements can provide valuable estimates of column interannual and seasonal variations for these compounds.

Description: Inferences from CO 2 and CH 4 concentration profiles at the Zotino Tall Tower Observatory (ZOTTO) on local summer-time ecosystem fluxes Biogeosciences Discussions, 10, 15337-15372, 2013 Author(s): J. Winderlich, C. Gerbig, O. Kolle, and M. Heimann The Siberian region is still sparsely covered by ecosystem observatories, which motivates to exploit existing datasets to gain spatially and temporally better-resolved carbon fluxes. The Zotino Tall Tower Observatory (ZOTTO, 60°48' N, 89°21' E) observations of CO 2 and CH 4 mole fractions as well as meteorological parameters from six different heights up to 301 m allow for an additional estimate of surface-atmosphere fluxes of CO 2 and CH 4 for the Middle-Siberian region since 2009. The total carbon flux is calculated from the storage and the turbulent flux component. The gradients between the different tower levels determine the storage flux component, which dominates the local fluxes, especially during night. As a correction term, the turbulent flux component was estimated by the modified Bowen ratio method based on the sensible heat flux measurements at the top of the tower. The gained average night time fluxes (23:00 to 04:00 local time) are 2.7 ± 1.1 μmol (m 2 s) −1 for CO 2 and 5.6 ± 4.5 nmol (m 2 s) −1 for CH 4 during the summer months June-September in 2009 and 2011. During day, the method is limited due to numeric instabilities from vanishing vertical gradients; however, the derived CO 2 fluxes exhibit reasonable diurnal shape and magnitude compared to the eddy covariance technique, which become available at the site in 2012. Therefore, the tall tower data facilitates the extension of the new eddy covariance flux dataset back in time. The diurnal signal of the CH 4 flux is predominantly characterized by a strong morning transition, which is explained by local topographic effects.

Description: Live foraminiferal faunas (Rose Bengal stained) from the northern Arabian Sea: links with bottom-water oxygenation Biogeosciences Discussions, 10, 15257-15304, 2013 Author(s): C. Caulle, K. A. Koho, M. Mojtahid, G. J. Reichart, and F. J. Jorissen Live (Rose Bengal stained) benthic foraminifera from the Murray Ridge, within and below the northern Arabian Sea Oxygen Minimum Zone (OMZ), were studied in order to determine the relationship between faunal composition, bottom-water oxygenation (BWO), pore-water chemistry and organic matter (organic carbon and phytopigment) distribution. A series of multicores were recovered from a ten-station oxygen (BWO: 2–78 μM) and bathymetric (885–3010 m depth) transect during the winter monsoon in January 2009. Foraminifera were investigated from three different size fractions (63–125 μm, 125–150 μm and 〉 150 μm). The larger foraminifera (〉 125 μm) were strongly dominated by agglutinated species (e.g. Reophax spp.). In contrast, in the 63–125 μm fraction, calcareous taxa were more abundant, especially in the core of the OMZ, suggesting an opportunistic behaviour. On the basis of a Principal Component Analysis, three foraminiferal groups were identified, reflecting the environmental parameters along the study transect. The faunas from the shallowest stations, in the core of the OMZ (BWO: 2 μM), were composed of "low oxygen" species, typical of the Arabian Sea OMZ (e.g., Rotaliatinopsis semiinvoluta , Praeglobobulimina spp. , Bulimina exilis, Uvigerina peregrina type parva ). These taxa are adapted to the very low BWO conditions and to high phytodetritus supplies. The transitional group, typical for the lower part of the OMZ (BWO: 5–16 μM), is composed of more cosmopolitan taxa tolerant to low-oxygen concentrations ( Globocassidulina subglobosa , Ehrenbergina trigona ). Below the OMZ (BWO: 26–78 μM), where food availability is more limited and becomes increasingly restricted to surficial sediments, more cosmopolitan calcareous taxa were present, such as Bulimina aculeata, Melonis barleeanus, Uvigerina peregrina and Epistominella exigua . Miliolids were uniquely observed in this last group, reflecting the higher BWO. At these deeper sites, the faunas exhibit a clear depth succession of superficial, intermediate and deep-infaunal microhabitats, because of the deeper oxygen and nitrate penetration into the sediment.

Description: Evidence of meltwater retention within the Greenland ice sheet The Cryosphere, 7, 1433-1445, 2013 Author(s): A. K. Rennermalm, L. C. Smith, V. W. Chu, J. E. Box, R. R. Forster, M. R. Van den Broeke, D. Van As, and S. E. Moustafa Greenland ice sheet mass losses have increased in recent decades with more than half of these attributed to surface meltwater runoff. However, the magnitudes of englacial storage, firn retention, internal refreezing and other hydrologic processes that delay or reduce true water export to the global ocean remain less understood, partly due to a scarcity of in situ measurements. Here, ice sheet surface meltwater runoff and proglacial river discharge between 2008 and 2010 near Kangerlussuaq, southwestern Greenland were used to establish sub- and englacial meltwater storage for a small ice sheet watershed (36–64 km 2 ). This watershed lacks significant potential meltwater storage in firn, surface lakes on the ice sheet and in the proglacial area, and receives limited proglacial precipitation. Thus, ice sheet surface runoff not accounted for by river discharge can reasonably be attributed to retention in sub- and englacial storage. Evidence for meltwater storage within the ice sheet includes (1) characteristic dampened daily river discharge amplitudes relative to ice sheet runoff; (2) three cold-season river discharge anomalies at times with limited ice sheet surface melt, demonstrating that meltwater may be retained up to 1–6 months; (3) annual ice sheet watershed runoff is not balanced by river discharge, and while near water budget closure is possible as much as 54% of melting season ice sheet runoff may not escape to downstream rivers; (4) even the large meltwater retention estimate (54%) is equivalent to less than 1% of the ice sheet volume, which suggests that storage in en- and subglacial cavities and till is plausible. While this study is the first to provide evidence for meltwater retention and delayed release within the Greenland ice sheet, more information is needed to establish how widespread this is along the Greenland ice sheet perimeter.

Description: Heterogeneous reaction of N 2 O 5 with illite and Arizona Test Dust particles Atmospheric Chemistry and Physics Discussions, 13, 24855-24884, 2013 Author(s): M. J. Tang, G. Schuster, and J. N. Crowley The heterogeneous reaction of N 2 O 5 with airborne illite and Arizona Test Dust particles was investigated at room temperature and at different relative humidities using an atmospheric pressure aerosol flow tube. N 2 O 5 at concentrations in the range 8 to 24×10 12 molecule cm −3 was monitored using thermal-dissociation cavity ring-down spectroscopy at 662 nm. At zero relative humidity a large uptake coefficient of N 2 O 5 to illite was obtained, γ(N 2 O 5 ) = 0.09, which decreased to 0.04 as relative humidity was increased to 67%. In contrast, the uptake coefficient derived for ATD is much lower (~ 0.006) and, within experimental uncertainty, independent of relative humidity (0–67%). Potential explanations are given for the significant differences between the uptake behaviour for ATD and illite and the results are compared with uptake coefficients for N 2 O 5 on other mineral surfaces.

Description: A case study of sea breeze blocking regulated by sea surface temperature along the English south coast Atmospheric Chemistry and Physics Discussions, 13, 24785-24807, 2013 Author(s): J. K. Sweeney, J. M. Chagnon, and S. L. Gray The sensitivity of sea breeze structure to sea surface temperature (SST) and coastal orography is investigated in convection-permitting Met Office Unified Model simulations of a case study along the south coast of England. Changes in SST of 1 K are shown to significantly modify the structure of the sea breeze. On the day of the case study the sea breeze was partially blocked by coastal orography, particularly within Lyme Bay. The extent to which the flow is blocked depends strongly on the static stability of the marine boundary layer. In experiments with colder SST, the marine boundary layer is more stable, and the degree of blocking is more pronounced. The implications of prescribing fixed SST from climatology in numerical weather prediction model forecasts of the sea breeze are discussed.

Description: Using a WRF simulation to examine regions where convection impacts the Asian summer monsoon anticyclone Atmospheric Chemistry and Physics Discussions, 13, 24809-24853, 2013 Author(s): N. K. Heath and H. E. Fuelberg The Asian summer monsoon is a prominent feature of the global circulation that is associated with an upper-level anticyclone (ULAC) that stands out vividly in satellite observations of trace gases. The ULAC also is an important region of troposphere-to-stratosphere transport. We ran the Weather Research and Forecasting (WRF) model at convective-permitting scales (4 km grid spacing) between 10–20 August 2012 to understand the role of convection in transporting boundary layer air into the upper-level anticyclone. Such high-resolution modeling of the Asian ULAC previously has not been documented in the literature. Comparison of our WRF simulation with reanalysis and satellite observations showed that WRF simulated the atmosphere sufficiently well to be used to study convective transport into the ULAC. A back-trajectory analysis based on hourly WRF output showed that 〉 90% of convectively influenced parcels reaching the ULAC came from the Tibetan Plateau (TP) and the southern slope (SS) of the Himalayas. A distinct diurnal cycle is seen in the convective trajectories, with their greatest impact occurring between 1600–2300 local solar time. This finding highlights the role of "everyday" diurnal convection in transporting boundary layer air into the ULAC. WRF output at 15 min intervals was produced for 16 August to examine the convection in greater detail. This high-temporal output revealed that the weakest convection in the study area occurred over the TP. However, because the TP is at 3000–5000 m a.m.s.l., its convection does not have to be as strong to reach the ULAC as in lower altitude regions. In addition, because the TP's elevated heat source is a major cause of the ULAC, we propose that convection over the TP and the neighboring SS is ideally situated geographically to impact the ULAC. The vertical mass flux of water vapor into the ULAC also was calculated. Results show that the TP and SS regions dominate other Asian regions in transporting moisture vertically into the ULAC. Because convection reaching the ULAC is more widespread over the TP than nearby, we propose that the abundant convection partially explains the TP's dominant water vapor fluxes. In addition, greater outgoing longwave radiation reaches the upper levels of the TP due to its elevated terrain. This creates a warmer ambient upper level environment, allowing parcels with greater saturation mixing ratios to enter the ULAC. Lakes in the Tibetan Plateau are shown to provide favorable conditions for deep convection during the night.

Description: A standard test case suite for two-dimensional linear transport on the sphere: results from a collection of state-of-the-art schemes Geoscientific Model Development Discussions, 6, 4983-5076, 2013 Author(s): P. H. Lauritzen, P. A. Ullrich, C. Jablonowski, P. A. Bosler, D. Calhoun, A. J. Conley, T. Enomoto, L. Dong, S. Dubey, O. Guba, A. B. Hansen, E. Kaas, J. Kent, J.-F. Lamarque, M. J. Prather, D. Reinert, V. V. Shashkin, W. C. Skamarock, B. Sørensen, M. A. Taylor, and M. A. Tolstykh Recently, a standard test case suite for 2-D linear transport on the sphere was proposed to assess important aspects of accuracy in geophysical fluid dynamics with a "minimal" set of idealized model configurations/runs/diagnostics. Here we present results from 19 state-of-the-art transport scheme formulations based on finite-difference/finite-volume methods as well as emerging (in the context of atmospheric/oceanographic sciences) Galerkin methods. Discretization grids range from traditional regular latitude-longitude grids to more isotropic domain discretizations such as icosahedral and cubed-sphere tessellations of the sphere. The schemes are evaluated using a wide range of diagnostics in idealized flow environments. Accuracy is assessed in single- and two-tracer configurations using conventional error norms as well as novel diagnostics designed for climate and climate-chemistry applications. In addition, algorithmic considerations that may be important for computational efficiency are reported on. The latter is inevitably computing platform dependent, The ensemble of results from a wide variety of schemes presented here helps shed light on the ability of the test case suite diagnostics and flow settings to discriminate between algorithms and provide insights into accuracy in the context of global atmospheric/ocean modeling. A library of benchmark results is provided to facilitate scheme intercomparison and model development. Simple software and data-sets are made available to facilitate the process of model evaluation and scheme intercomparison.

Description: Spatial and seasonal variations in evapotranspiration over Canada's landmass Hydrology and Earth System Sciences, 17, 3561-3575, 2013 Author(s): S. Wang, Y. Yang, Y. Luo, and A. Rivera A 30 yr (1979–2008) dataset of actual evapotranspiration (ET) at 1 km resolution was generated over Canada's landmass by integrating remote sensing land surface data and gridded climate data using the EALCO model run at a 30 min time step. This long-term high-resolution dataset was used to characterize the spatiotemporal variations in ET across Canada. The results show that annual ET varied from 600 mm yr −1 over several regions in the south to less than 100 mm yr −1 in the northern Arctic. Nationally, ET in summer (i.e., June to August) comprised 65% of the annual total amount. ET in the cold season remained mostly below 10 mm month −1 over the country. Negative monthly ET was obtained over the Arctic region in winter, indicating EALCO simulated a larger amount of condensation than ET. Overall, the mean ET over the entire Canadian landmass for the 30 yr was 239 mm yr −1 , or 44% of its corresponding precipitation. Comparisons of available ET studies in Canada revealed large uncertainties in ET estimates associated with using different approaches. The scarcity of ET measurements for the diverse ecosystems in Canada remains a significant challenge for reducing the uncertainties; this gap needs to be addressed in future studies to improve capabilities in climate/weather modeling and water resource management.

Description: A spatial bootstrap technique for parameter estimation of rainfall annual maxima distribution Hydrology and Earth System Sciences Discussions, 10, 11755-11794, 2013 Author(s): F. Uboldi, A. N. Sulis, C. Lussana, M. Cislaghi, and M. Russo Estimation of extreme event distributions and depth-duration-frequency (DDF) curves is achieved at any target site by repeated sampling among all available raingauge data in the surrounding area. The estimate is computed over a gridded domain in Northern Italy, using precipitation time series from 1929 to 2011, including data from historical analog stations and from the present-day automatic observational network. The presented local regionalisation naturally overcomes traditional station-point methods, with their demand of long historical series and their sensitivity to very rare events occurring at very few stations, possibly causing unrealistic spatial gradients in DDF relations. At the same time, the presented approach allows for spatial dependence, necessary in a geographical domain such as Lombardy, complex for both its topography and its climatology. The bootstrap technique enables evaluating uncertainty maps for all estimated parameters and for rainfall depths at assigned return periods.

Description: Overview of the first HyMeX Special Observation Period over Italy: observations and model results Hydrology and Earth System Sciences Discussions, 10, 11643-11710, 2013 Author(s): R. Ferretti, E. Pichelli, S. Gentile, I. Maiello, D. Cimini, S. Davolio, M. M. Miglietta, G. Panegrossi, L. Baldini, F. Pasi, F. S. Marzano, A. Zinzi, S. Mariani, M. Casaioli, G. Bartolini, N. Loglisci, A. Montani, C. Marsigli, A. Manzato, A. Pucillo, M. E. Ferrario, V. Colaiuda, and R. Rotunno During the first Hymex campaign (5 September–6 November 2012) referred to as Special Observation Period (SOP-1), dedicated to heavy precipitation events and flash floods in Western Mediterranean, three Italian hydro-meteorological monitoring sites were activated: Liguria-Tuscany, North-Eastern Italy and Central Italy. The extraordinary deployment of advanced instrumentation, including instrumented aircrafts, and the use of several different operational weather forecast models has allowed an unprecedented monitoring and analysis of high impact weather events around the Italian hydro-meteorological sites. This activity has seen the strict collaboration between the Italian scientific and operational communities. In this paper, an overview of the Italian organization during the SOP-1 is provided, and selected Intensive Observation Periods (IOPs) are described. A significant event for each Italian target area is chosen for this analysis: IOP2 (12–13 September 2012) in North-Eastern Italy, IOP13 (15–16 October 2012) in Central Italy and IOP19 (3–5 November 2012) in Liguria and Tuscany. For each IOP the meteorological characteristics, together with special observations and weather forecasts, are analyzed with the aim of highlighting strengths and weaknesses of the forecast modeling systems. Moreover, using one of the three events, the usefulness of different operational chains is highlighted.

Description: Correcting basin-scale snowfall in a mountainous basin using a distributed snowmelt model and remote sensing data Hydrology and Earth System Sciences Discussions, 10, 11711-11753, 2013 Author(s): M. Shrestha, L. Wang, T. Koike, H. Tsutsui, Y. Xue, and Y. Hirabayashi Adequate estimation of the spatial distribution of snowfall is critical in hydrologic modeling. However, this is a well-known problem in estimating basin-scale snowfall, especially in mountainous basins with data scarcity. This study focuses on correction and estimation of this spatial distribution, which considers topographic effects within the basin. A method is proposed that optimizes an altitude-based snowfall correction factor ( C fsnow ). This is done through multi-objective calibration of a spatially distributed, multilayer energy and water balance-based snowmelt model (WEB-DHM-S) with observed discharge and remotely sensed snow cover data from the Moderate Resolution Imaging Spectroradiometer (MODIS). The Shuffled Complex Evolution – University of Arizona automatic search algorithm is used to obtain the optimal value of C fsnow for minimum cumulative error in discharge and snow cover simulations. Discharge error is quantified by Nash–Sutcliffe efficiency and relative volume deviation, and snow cover error was estimated by pixel-by-pixel analysis. The study region is the heavily snow-fed Yagisawa Basin of the Upper Tone River in northeast Japan. First, the system was applied to one snow season (2002–2003), obtaining an optimized C fsnow of 0.0007 m −1 . For validation purposes, the optimized C fsnow was implemented to correct snowfall in 2004, 2002 and 2001. Overall, the system was effective, implying improvements in correlation of simulated vs. observed discharge and snow cover. The 4 yr mean of basin-average snowfall for the corrected spatial snowfall distribution was 1160 mm (780 mm before correction). Execution of sensitivity runs against other model input and parameters indicated that C fsnow could be affected by uncertainty in shortwave radiation and setting of the threshold air temperature parameter. Our approach is suitable to correct snowfall and estimate its distribution in poorly-gauged basins, where elevation dependence of snowfall amount is strong.

Description: High-resolution provenance of desert dust deposited on Mt. Elbrus, Caucasus in 2009–2012 using snow pit and firn core records The Cryosphere, 7, 1481-1498, 2013 Author(s): S. Kutuzov, M. Shahgedanova, V. Mikhalenko, P. Ginot, I. Lavrentiev, and S. Kemp The first record of dust deposition events on Mt. Elbrus, Caucasus Mountains derived from a snow pit and a shallow firn core is presented for the 2009–2012 period. A combination of isotopic analysis, SEVIRI red-green-blue composite imagery, MODIS atmospheric optical depth fields derived using the Deep Blue algorithm, air mass trajectories derived using the HYSPLIT model and analyses of meteorological data enabled identification of dust source regions with high temporal (hours) and spatial (ca. 20–100 km) resolution. Seventeen dust deposition events were detected; fourteen occurred in March–June, one in February and two in October. Four events originated in the Sahara, predominantly in northeastern Libya and eastern Algeria. Thirteen events originated in the Middle East, in the Syrian Desert and northern Mesopotamia, from a mixture of natural and anthropogenic sources. Dust transportation from Sahara was associated with vigorous Saharan depressions, strong surface winds in the source region and mid-tropospheric southwesterly flow with daily winds speeds of 20–30 m s −1 at 700 hPa level. Although these events were less frequent than those originating in the Middle East, they resulted in higher dust concentrations in snow. Dust transportation from the Middle East was associated with weaker depressions forming over the source region, high pressure centred over or extending towards the Caspian Sea and a weaker southerly or southeasterly flow towards the Caucasus Mountains with daily wind speeds of 12–18 m s −1 at 700 hPa level. Higher concentrations of nitrates and ammonium characterised dust from the Middle East deposited on Mt. Elbrus in 2009 indicating contribution of anthropogenic sources. The modal values of particle size distributions ranged between 1.98 μm and 4.16 μm. Most samples were characterised by modal values of 2.0–2.8 μm with an average of 2.6 μm and there was no significant difference between dust from the Sahara and the Middle East.

Description: Simulating the role of gravel on the dynamics of permafrost on the Qinghai-Tibetan Plateau The Cryosphere Discussions, 7, 4703-4740, 2013 Author(s): S. Yi, J. Chen, Q. Wu, and Y. Ding Gravel (particle size ≥ 2 mm) is common in soil profiles of the Qinghai-Tibetan Plateau (QTP). It has different thermal and hydrological properties than other fine mineral soils (particle size 〈 2 mm), which may have significant impacts on the thermal and hydrological processes of soil. However, few models have considered gravel. In this study, we implemented the thermal and hydraulic properties of gravel into the Dynamic Organic Soil-Terrestrial Ecosystem Model to develop new schemes to simulate the dynamics of permafrost on the QTP. Results showed that: (1) the widely used Farouki thermal scheme always simulated higher thermal conductivity of frozen soils than unfrozen soils with the same soil water content; therefore it tends to overestimate permafrost thickness strongly; (2) there exists a soil moisture threshold, below which the new set of schemes with gravel simulated smaller thermal conductivity of frozen soils than unfrozen soils; (3) soil with gravel has higher hydraulic conductivity and poorer water retention capability; and simulations with gravel were usually drier than those without gravel; and (4) the new schemes simulated faster upward degradation than downward degradation; and the simulated permafrost thicknesses were sensitive to the fraction of gravel, the gravel size, the thickness of soil with gravel, and the subsurface drainage. To reduce the uncertainties in the projection of permafrost degradation on the QTP, more effort should be made to: (1) developing robust relationships between soil thermal and hydraulic properties and gravel characteristics based on laboratory work; and (2) compiling spatial datasets of the vertical distribution of gravel content based on measurements during drilling or the digging of soil pits.

Description: Changes in glacier equilibrium-line altitude in the western Alps from 1984 to 2010: evaluation by remote sensing and modeling of the morpho-topographic and climate controls The Cryosphere, 7, 1455-1471, 2013 Author(s): A. Rabatel, A. Letréguilly, J.-P. Dedieu, and N. Eckert We present time series of equilibrium-line altitude (ELA) measured from the end-of-summer snow line altitude computed using satellite images, for 43 glaciers in the western Alps over the 1984–2010 period. More than 120 satellite images acquired by Landsat, SPOT and ASTER were used. In parallel, changes in climate variables, summer cumulative positive degree days (CPDD) and winter precipitation, were analyzed over the same time period using 22 weather stations located inside and around the study area. Assuming a continuous linear trend over the study period: (1) the average ELA of the 43 glaciers increased by about 170 m; (2) summer CPDD increased by about 150 PDD at 3000 m a.s.l.; and (3) winter precipitation remained rather stationary. Summer CPDD showed homogeneous spatial and temporal variability; winter precipitation showed homogeneous temporal variability, but some stations showed a slightly different spatial pattern. Regarding ELAs, temporal variability between the 43 glaciers was also homogeneous, but spatially, glaciers in the southern part of the study area differed from glaciers in the northern part, mainly due to a different precipitation pattern. A sensitivity analysis of the ELAs to climate and morpho-topographic variables (elevation, aspect, latitude) highlighted the following: (1) the average ELA over the study period of each glacier is strongly controlled by morpho-topographic variables; and (2) the interannual variability of the ELA is strongly controlled by climate variables, with the observed increasing trend mainly driven by increasing temperatures, even if significant nonlinear, low-frequency fluctuations appear to be driven by winter precipitation anomalies. Finally, we used an expansion of Lliboutry's approach to reconstruct fluctuations in the ELA of any glacier of the study area with respect to morpho-topographic and climate variables, by quantifying their respective weight and the related uncertainties in a consistent manner within a hierarchical Bayesian framework. This method was tested and validated using the ELA measured on the satellite images.

Description: Identification of response and timing issues at permanent European broadband stations from automated data analysis Advances in Geosciences, 36, 21-25, 2013 Author(s): C. Weidle, R. A. Soomro, L. Cristiano, and T. Meier We apply automated two-station broadband phase velocity dispersion measurements to all available broadband data from permanent seismic stations in Europe, as available through the European Integrated Data Archive (EIDA, http://www.orfeus-eu.org/eida/ ) infrastructure. As part of our quality control we detect several typical patterns in our measurements that can be related to technical problems, incorrect metadata information or uncover inconsistencies in data processing routines. These effects include timing and various response issues, most prominently erroneous response information. Our procedure is thus able to identify potentially problematic (meta)data from a large set of seismic data and offers an applicable way to increase data quality at data centers.

Description: Introduction to special issue: the TransBrom Sonne expedition in the tropical West Pacific Atmospheric Chemistry and Physics, 13, 9439-9446, 2013 Author(s): K. Krüger and B. Quack This special section of Atmospheric Chemistry and Physics gives an overview of scientific results, collected during a West Pacific ship expedition in October 2009 with the Research Vessel (R/V) Sonne . The cruise focussed on chemical interactions between the ocean surface and the atmosphere above the tropical West Pacific and was planned within the national research project TransBrom ( www.geomar.de/~transbrom ). TransBrom aimed to particularly investigate very short lived bromine compounds in the ocean and their transport to and relevance for the stratosphere. For this purpose, chemical and biological parameters were analysed in the ocean and in the atmosphere, accompanied by a high frequency of meteorological measurements, to derive new insights into the multidisciplinary research field. This introduction paper presents the scientific goals and the meteorological and oceanographic background. The main research findings of the TransBrom Sonne expedition are highlighted.

Description: On the analytic approximation of bulk collision rates of non-spherical hydrometeors Geoscientific Model Development Discussions, 6, 5077-5116, 2013 Author(s): A. Seifert, U. Blahak, and R. Buhr Analytic approximations of the binary collision rates of hydrometeors are derived for use in bulk microphysical parameterizations. Special attention is given to non-spherical hydrometeors like raindrops and snowflakes. The terminal fall velocity of these particles cannot be sufficiently well approximated by power law relations which are used in most microphysical parameterizations and therefore an improved formulation is needed. The analytic approximations of the bulk collision rates given in this paper are an alternative to look-up tables and can replace the Wisner approximation which is used in many atmospheric models.

Description: Attribution of hydrologic forecast uncertainty within scalable forecast windows Hydrology and Earth System Sciences Discussions, 10, 11795-11828, 2013 Author(s): L. Yang, F. Tian, Y. Sun, X. Yuan, and H. Hu Hindcasts based on the Extended Streamflow Prediction (ESP) approach are carried out in a typical rainfall-dominated basin in China, aiming to examine the roles of initial condition (IC), future atmospheric forcing (FC) and hydrologic model uncertainty (MU) in the streamflow forecast skill. The combined effects of IC and FC are explored within the framework of a forecast window. By implementing virtual numerical simulations without the consideration of MU, it is found that the dominance of IC could last up to 90 days in dry season, while its impact gives way to FC for lead times exceeding 30 days in the wet season. The combined effects of IC and FC on the forecast skill are further investigated by proposing a dimensionless parameter ( β ) that represents the ratio of the total amount of initial water storage and the incoming rainfall. The forecast skill increases exponentially with β , and varies greatly in different forecast windows. Moreover, the influence of MU on forecast skill is examined by focusing on the uncertainty of model parameters. Two different hydrologic model calibration strategies are carried out. The results indicate that the uncertainty of model parameters exhibits a more significant influence on the forecast skill in the dry season than in the wet season. The ESP approach is more skillful in monthly streamflow forecast during the transition period from wet to dry than otherwise. For the transition period from dry to wet, the low skill of the forecasts could be attributed to the combined effects of IC and FC, but less to the biases in the hydrologic model parameters. For the forecasting in dry season, the usefulness of the ESP approach is heavily dependent on the strategy of the model calibration.

Description: Integrating spatial, temporal, and size probabilities for the annual landslide hazard maps in the Shihmen watershed, Taiwan Natural Hazards and Earth System Science, 13, 2353-2367, 2013 Author(s): C. Y. Wu and S. C. Chen Landslide spatial, temporal, and size probabilities were used to perform a landslide hazard assessment in this study. Eleven intrinsic geomorphological, and two extrinsic rainfall factors were evaluated as landslide susceptibility related factors as they related to the success rate curves, landslide ratio plots, frequency distributions of landslide and non-landslide groups, as well as probability–probability plots. Data on landslides caused by Typhoon Aere in the Shihmen watershed were selected to train the susceptibility model. The landslide area probability, based on the power law relationship between the landslide area and a noncumulative number, was analyzed using the Pearson type 5 probability density function. The exceedance probabilities of rainfall with various recurrence intervals, including 2, 5, 10, 20, 50, 100 and 200 yr, were used to determine the temporal probabilities of the events. The study was conducted in the Shihmen watershed, which has an area of 760 km 2 and is one of the main water sources for northern Taiwan. The validation result of Typhoon Krosa demonstrated that this landslide hazard model could be used to predict the landslide probabilities. The results suggested that integration of spatial, area, and exceedance probabilities to estimate the annual probability of each slope unit is feasible. The advantage of this annual landslide probability model lies in its ability to estimate the annual landslide risk, instead of a scenario-based risk.

Description: A top-down model to generate ensembles of runoff from a large number of hillslopes Nonlinear Processes in Geophysics, 20, 683-704, 2013 Author(s): P. R. Furey, V. K. Gupta, and B. M. Troutman We hypothesize that total hillslope water loss for a rainfall–runoff event is inversely related to a function of a lognormal random variable, based on basin- and point-scale observations taken from the 21 km 2 Goodwin Creek Experimental Watershed (GCEW) in Mississippi, USA. A top-down approach is used to develop a new runoff generation model both to test our physical-statistical hypothesis and to provide a method of generating ensembles of runoff from a large number of hillslopes in a basin. The model is based on the assumption that the probability distributions of a runoff/loss ratio have a space–time rescaling property. We test this assumption using streamflow and rainfall data from GCEW. For over 100 rainfall–runoff events, we find that the spatial probability distributions of a runoff/loss ratio can be rescaled to a new distribution that is common to all events. We interpret random within-event differences in runoff/loss ratios in the model to arise from soil moisture spatial variability. Observations of water loss during events in GCEW support this interpretation. Our model preserves water balance in a mean statistical sense and supports our hypothesis. As an example, we use the model to generate ensembles of runoff at a large number of hillslopes for a rainfall–runoff event in GCEW.

Description: A mechanism for catastrophic filter divergence in data assimilation for sparse observation networks Nonlinear Processes in Geophysics, 20, 705-712, 2013 Author(s): G. A. Gottwald and A. J. Majda We study catastrophic filter divergence in data assimilation procedures whereby the forecast model develops severe numerical instabilities leading to a blow-up of the solution. Catastrophic filter divergence can occur in sparse observational grids with small observational noise for intermediate observation intervals and finite ensemble sizes. Using a minimal five-dimensional model, we establish that catastrophic filter divergence is a numerical instability of the underlying forecast model caused by the filtering procedure producing analyses which are not consistent with the true dynamics, and stiffness caused by the fast attraction of the inconsistent analyses towards the attractor during the forecast step.

Description: Quantitative determination of carbonaceous particle mixing state in Paris using single-particle mass spectrometer and aerosol mass spectrometer measurements Atmospheric Chemistry and Physics, 13, 9479-9496, 2013 Author(s): R. M. Healy, J. Sciare, L. Poulain, M. Crippa, A. Wiedensohler, A. S. H. Prévôt, U. Baltensperger, R. Sarda-Estève, M. L. McGuire, C.-H. Jeong, E. McGillicuddy, I. P. O'Connor, J. R. Sodeau, G. J. Evans, and J. C. Wenger Single-particle mixing state information can be a powerful tool for assessing the relative impact of local and regional sources of ambient particulate matter in urban environments. However, quantitative mixing state data are challenging to obtain using single-particle mass spectrometers. In this study, the quantitative chemical composition of carbonaceous single particles has been determined using an aerosol time-of-flight mass spectrometer (ATOFMS) as part of the MEGAPOLI 2010 winter campaign in Paris, France. Relative peak areas of marker ions for elemental carbon (EC), organic aerosol (OA), ammonium, nitrate, sulfate and potassium were compared with concurrent measurements from an Aerodyne high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS), a thermal–optical OCEC analyser and a particle into liquid sampler coupled with ion chromatography (PILS-IC). ATOFMS-derived estimated mass concentrations reproduced the variability of these species well ( R 2 = 0.67–0.78), and 10 discrete mixing states for carbonaceous particles were identified and quantified. The chemical mixing state of HR-ToF-AMS organic aerosol factors, resolved using positive matrix factorisation, was also investigated through comparison with the ATOFMS dataset. The results indicate that hydrocarbon-like OA (HOA) detected in Paris is associated with two EC-rich mixing states which differ in their relative sulfate content, while fresh biomass burning OA (BBOA) is associated with two mixing states which differ significantly in their OA / EC ratios. Aged biomass burning OA (OOA 2 -BBOA) was found to be significantly internally mixed with nitrate, while secondary, oxidised OA (OOA) was associated with five particle mixing states, each exhibiting different relative secondary inorganic ion content. Externally mixed secondary organic aerosol was not observed. These findings demonstrate the range of primary and secondary organic aerosol mixing states in Paris. Examination of the temporal behaviour and chemical composition of the ATOFMS classes also enabled estimation of the relative contribution of transported emissions of each chemical species and total particle mass in the size range investigated. Only 22% of the total ATOFMS-derived particle mass was apportioned to fresh, local emissions, with 78% apportioned to regional/continental-scale emissions.

Description: Inter-calibration of polar imager solar channels using SEVIRI Atmospheric Measurement Techniques, 6, 2495-2508, 2013 Author(s): J. F. Meirink, R. A. Roebeling, and P. Stammes Accurate calibration of satellite imagers is a prerequisite for using their measurements in climate applications. Here we present a method for the inter-calibration of geostationary and polar-orbiting imager solar channels based on regressions of collocated near-nadir reflectances. Specific attention is paid to correcting for differences in spectral response between instruments. The method is used to calibrate the solar channels of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on the geostationary Meteosat satellite with corresponding channels of the Moderate Resolution Imaging Spectroradiometer (MODIS) on the polar-orbiting Aqua satellite. The SEVIRI operational calibration is found to be stable during the years 2004 to 2009, but offset by −8, −6, and +3.5 % for channels 1 (0.6 μm), 2 (0.8 μm), and 3 (1.6 μm), respectively. These results are robust for a range of choices that can be made regarding data collocation and selection, as long as the viewing and illumination geometries of the two instruments are matched. Uncertainties in the inter-calibration method are estimated to be 1 % for channel 1 and 1.5 % for channels 2 and 3. A specific application of our method is the inter-calibration of polar imagers using SEVIRI as a transfer instrument. This offers an alternative to direct inter-calibration, which in general has to rely on high-latitude collocations. Using this method we have tied MODIS-Terra and Advanced Very High Resolution Radiometer (AVHRR) instruments on National Oceanic and Atmospheric Administration (NOAA) satellites 17 and 18 to MODIS-Aqua for the years 2007 to 2009. While reflectances of the two MODIS instruments differ less than 2 % for all channels considered, deviations of an existing AVHRR calibration from MODIS-Aqua reach −3.5 and +2.5 % for the 0.8 and 1.6 μm channels, respectively.

Description: An assessment of the performance of the Monitor for AeRosols and GAses in ambient air (MARGA): a semi-continuous method for soluble compounds Atmospheric Chemistry and Physics Discussions, 13, 25067-25124, 2013 Author(s): I. C. Rumsey, K. A. Cowen, J. T. Walker, T. J. Kelly, E. A. Hanft, K. Mishoe, C. Rogers, R. Proost, G. M. Beachley, G. Lear, T. Frelink, and R. P. Otjes Ambient air monitoring as part of the US Environmental Protection Agency's (US EPA's) Clean Air Status and Trends Network (CASTNet) currently uses filter packs to measure weekly integrated concentrations. The US EPA is interested in supplementing CASTNet with semi-continuous monitoring systems at select sites to characterize atmospheric chemistry and deposition of nitrogen and sulfur compounds at higher time resolution than the filter pack. The Monitor for AeRosols and GAses in ambient air (MARGA) measures water-soluble gases and aerosols at hourly temporal resolution. The performance of the MARGA was assessed under the US EPA Environmental Technology Verification (ETV) program. The assessment was conducted in Research Triangle Park, North Carolina from 8 September–8 October 2010 and focused on gaseous SO 2 , HNO 3 and NH 3 and aerosol SO 4 − , NO 3 − and NH 4 + . Precision of the MARGA was evaluated by calculating the median absolute relative percent difference (MARPD) between paired hourly results from duplicate MARGA units (MUs), with a performance goal of

Description: Using ocean-glint scattered sunlight as a diagnostic tool for satellite remote sensing of greenhouse gases Atmospheric Measurement Techniques, 6, 2509-2520, 2013 Author(s): A. Butz, S. Guerlet, O. P. Hasekamp, A. Kuze, and H. Suto Spectroscopic measurements of sunlight backscattered by the Earth's surface is a technique widely used for remote sensing of atmospheric constituent concentrations from space. Thereby, remote sensing of greenhouse gases poses particularly challenging accuracy requirements for instrumentation and retrieval algorithms which, in general, suffer from various error sources. Here, we investigate a method that helps disentangle sources of error for observations of sunlight backscattered from the glint spot on the ocean surface. The method exploits the backscattering characteristics of the ocean surface, which is bright for glint geometry but dark for off-glint angles. This property allows for identifying a set of clean scenes where light scattering due to particles in the atmosphere is negligible such that uncertain knowledge of the lightpath can be excluded as a source of error. We apply the method to more than 3 yr of ocean-glint measurements by the Thermal And Near infrared Sensor for carbon Observation (TANSO) Fourier Transform Spectrometer (FTS) onboard the Greenhouse Gases Observing Satellite (GOSAT), which aims at measuring carbon dioxide (CO 2 ) and methane (CH 4 ) concentrations. The proposed method is able to clearly monitor recent improvements in the instrument calibration of the oxygen (O 2 ) A-band channel and suggests some residual uncertainty in our knowledge about the instrument. We further assess the consistency of CO 2 retrievals from several absorption bands between 6400 cm −1 (1565 nm) and 4800 cm −1 (2100 nm) and find that the absorption bands commonly used for monitoring of CO 2 dry air mole fractions from GOSAT allow for consistency better than 1.5 ppm. Usage of other bands reveals significant inconsistency among retrieved CO 2 concentrations pointing at inconsistency of spectroscopic parameters.

Description: Landscape-scale changes in forest structure and functional traits along an Andes-to-Amazon elevation gradient Biogeosciences Discussions, 10, 15415-15454, 2013 Author(s): G. P. Asner, C. Anderson, R. E. Martin, D. E. Knapp, R. Tupayachi, T. Kennedy-Bowdoin, F. Sinca, and Y. Malhi Elevation gradients provide opportunities to explore environmental controls on forest structure and functioning, but plot-based studies have proven highly variable due to limited geographic scope. We used airborne imaging spectroscopy and LiDAR (light detection and ranging) to quantify changes in three-dimensional forest structure and canopy functional traits in a series of 25 ha landscapes distributed along a 3300 m elevation gradient from lowland Amazonia to treeline in the Peruvian Andes. Canopy greenness, photosynthetic fractional cover and exposed non-photosynthetic vegetation varied as much across lowland forests (100–200 m) as they did from the lowlands to the Andean treeline (3400 m). Elevation was positively correlated with canopy gap density and understory vegetation cover, and negatively related to canopy height and vertical profile. Increases in gap density were tightly linked to increases in understory plant cover, and larger gaps (20–200 m 2 produced 25–30 times the response in understory cover than did smaller gaps ( 〈 5 m 2 . Scaling of gap size to gap frequency was, however, relatively constant along the elevation gradient, which when combined with other canopy structural information, indicates equilibrium turnover patterns from the lowlands to treeline. Our results provide a first landscape-scale quantification of forest structure and canopy functional traits with changing elevation, thereby improving our understanding of disturbance, demography and ecosystem processes in the Andes-to-Amazon corridor.

Description: Post-processing rainfall forecasts from numerical weather prediction models for short-term streamflow forecasting Hydrology and Earth System Sciences, 17, 3587-3603, 2013 Author(s): D. E. Robertson, D. L. Shrestha, and Q. J. Wang Sub-daily ensemble rainfall forecasts that are bias free and reliably quantify forecast uncertainty are critical for flood and short-term ensemble streamflow forecasting. Post-processing of rainfall predictions from numerical weather prediction models is typically required to provide rainfall forecasts with these properties. In this paper, a new approach to generate ensemble rainfall forecasts by post-processing raw numerical weather prediction (NWP) rainfall predictions is introduced. The approach uses a simplified version of the Bayesian joint probability modelling approach to produce forecast probability distributions for individual locations and forecast lead times. Ensemble forecasts with appropriate spatial and temporal correlations are then generated by linking samples from the forecast probability distributions using the Schaake shuffle. The new approach is evaluated by applying it to post-process predictions from the ACCESS-R numerical weather prediction model at rain gauge locations in the Ovens catchment in southern Australia. The joint distribution of NWP predicted and observed rainfall is shown to be well described by the assumed log-sinh transformed bivariate normal distribution. Ensemble forecasts produced using the approach are shown to be more skilful than the raw NWP predictions both for individual forecast lead times and for cumulative totals throughout all forecast lead times. Skill increases result from the correction of not only the mean bias, but also biases conditional on the magnitude of the NWP rainfall prediction. The post-processed forecast ensembles are demonstrated to successfully discriminate between events and non-events for both small and large rainfall occurrences, and reliably quantify the forecast uncertainty. Future work will assess the efficacy of the post-processing method for a wider range of climatic conditions and also investigate the benefits of using post-processed rainfall forecasts for flood and short-term streamflow forecasting.

Description: On the puzzling feature of the silence of precursory electromagnetic emissions Natural Hazards and Earth System Science, 13, 2381-2397, 2013 Author(s): K. Eftaxias, S. M. Potirakis, and T. Chelidze It has been suggested that fracture-induced MHz–kHz electromagnetic emissions (EME), which emerge from a few days up to a few hours before the main seismic shock occurrence permit a real-time monitoring of the damage process during the last stages of earthquake preparation, as it happens at the laboratory scale. Despite fairly abundant evidence, electromagnetic (EM) precursors have not been adequately accepted as credible physical phenomena. These negative views are enhanced by the fact that certain "puzzling features" are repetitively observed in candidate fracture-induced pre-seismic EME. More precisely, EM silence in all frequency bands appears before the main seismic shock occurrence, as well as during the aftershock period. Actually, the view that "acceptance of "precursive" EM signals without convincing co-seismic signals should not be expected" seems to be reasonable. In this work we focus on this point. We examine whether the aforementioned features of EM silence are really puzzling ones or, instead, reflect well-documented characteristic features of the fracture process, in terms of universal structural patterns of the fracture process, recent laboratory experiments, numerical and theoretical studies of fracture dynamics, critical phenomena, percolation theory, and micromechanics of granular materials. Our analysis shows that these features should not be considered puzzling.

Description: A note on the water budget of temperate glaciers The Cryosphere, 7, 1557-1564, 2013 Author(s): J. Oerlemans In this note, the total dissipative melting in temperate glaciers is studied. The analysis is based on the notion that the dissipation is determined by the loss of potential energy due to the downward motion of mass (ice, snow, meltwater and rain). A mathematical formulation of the dissipation is developed and applied to a simple glacier geometry. In the next step, meltwater production resulting from enhanced ice motion during a glacier surge is calculated. The amount of melt energy available follows directly from the lowering of the centre of gravity of the glacier. To illustrate the concept, schematic calculations are presented for a number of glaciers with different geometric characteristics. Typical dissipative melt rates, expressed as water-layer depth averaged over the glacier, range from a few centimetres per year for smaller glaciers to half a metre per year for Franz Josef Glacier, one of the most active glaciers in the world (in terms of mass turnover). The total generation of meltwater during a surge is typically half a metre. For Variegated Glacier a value of 70 cm is found, for Kongsvegen 20 cm. These values refer to water layer depth averaged over the entire glacier. The melt \textit{rate} depends on the duration of the surge. It is generally an order of magnitude greater than water production by `normal' dissipation. On the other hand, the additional basal melt rate during a surge is comparable in magnitude with the water input from meltwater and precipitation. This suggests that enhanced melting during a surge does not grossly change the total water budget of a glacier. Basal water generated by enhanced sliding is an important ingredient in many theories of glacier surges. It provides a positive feedback mechanism that actually makes the surge happen. The results found here suggest that this can only work if water generated by enhanced sliding accumulates in a part of the glacier base where surface meltwater and rain have no or very limited access. This finding seems compatible with the fact that, on many glaciers, surges are initiated in the lower accumulation zone.

Description: Migrating subtropical front and Agulhas Return Current affect the southwestern Indian Ocean during the late Quaternary Climate of the Past Discussions, 9, 5521-5551, 2013 Author(s): D. K. Naik, R. Saraswat, N. Khare, A. C. Pandey, and R. Nigam The position of sub-tropical front (STF), Agulhas Current (AC) and Agulhas Return Current (ARC) controls the hydrography of southwestern Indian Ocean. Although, equator-ward migration of STF and reduction in Agulhas leakage has been reported during the last glacial period, the fate of ARC during the last glacial–interglacial cycle is not clear. Therefore, in order to understand changes in the position and strength of ARC during the last glacial–interglacial cycle, here we reconstruct hydrographic changes in the southwestern Indian Ocean from temporal variation in planktic foraminiferal abundance, stable isotopic ratio (δ 18 O) and trace metal ratio (Mg / Ca) of planktic foraminifera Globigerina bulloides in a core collected from the Agulhas Retroflection Region (ARR) in the southwestern Indian Ocean. Increased abundance of G. bulloides suggests that the productivity in the southwestern Indian Ocean increased during glacial period which confirms previous reports of high glacial productivity in the Southern Ocean. The increased productivity was likely driven by a combination of equator-ward migration of subtropical front and westerlies. Increase in relative abundance of Neogloboquadrina pachyderma Dextral suggests warming of ARR leading to strong thermocline in the southwestern Indian Ocean during the last glacial period. We suggest that the warming of Agulhas Retroflection Region was driven by strengthened ARC which shifted to the east of its present location, thus bringing warmer and saltier water to the southwestern Indian Ocean. Therefore, it is inferred that over the last glacial–interglacial cycle, the hydrography of southwestern Indian Ocean was driven by an eastward shift of retroflection region as well as migrating subtropical front.

Description: The response of the Peruvian Upwelling Ecosystem to centennial-scale global change during the last two millennia Climate of the Past Discussions, 9, 5479-5519, 2013 Author(s): R. Salvatteci, D. Gutiérrez, D. Field, A. Sifeddine, L. Ortlieb, I. Bouloubassi, M. Boussafir, H. Boucher, and F. Cetin The Tropical Pacific ocean-atmosphere system influences global climate on interannual, decadal, as well as at longer timescales. Given the uncertainties in the response of the Tropical Pacific to the ongoing greenhouse effect, it is important to assess the natural range of the Tropical Pacific climate variability in response to global natural changes, and to understand the underlying mechanisms. The Peruvian Upwelling Ecosystem (PUE) represents an ideal area to reconstruct past changes in ocean-atmosphere systems because productivity and subsurface oxygenation are strongly linked to changes in the strength of the Walker circulation. Throughout the last 2000 yr, warmer (the Roman Warm Period [RWP], the Medieval Climate Anomaly [MCA] and the Current Warm Period [CWP]), and colder (the Dark Ages Cold Period [DACP] and Little Ice Age [LIA]) intervals occurred with considerable changes around the globe. In order to reconstruct the PUE response to these climatic periods and reveal the underlying mechanisms, we use a multi-proxy approach including organic and inorganic proxies in finely laminated sediments retrieved off Pisco (~ 14° S), Peru. Our results indicate that the PUE exhibited a La Niña-like mean state during the warm periods, characterized by an intense OMZ and high marine productivity. During cold periods the PUE exhibited an El Niño-like mean state, characterized by a weak OMZ and low marine productivity. Comparing our results with other relevant paleoclimatic reconstructions revealed that changes in the strength of the Walker circulation and the expansion/contraction of the South Pacific Sub-tropical High controlled productivity and subsurface oxygenation in the PUE during the last two millennia. This indicate that large scale circulation changes are the driving forces in maintaining productivity and subsurface oxygenation off Peru at centennial time scales during the past two millennia.

Description: A conceptual model of check dam hydraulics for gully control Hydrology and Earth System Sciences Discussions, 10, 11901-11941, 2013 Author(s): C. Castillo, R. Pérez, and J. A. Gómez There is little information in scientific literature regarding the modifications induced by check dam systems in flow regimes in restored gully reaches, despite it being a crucial issue for the design of conservation measures. Here, we develop a conceptual model to classify flow regimes in straight rectangular channels for initial and dam-filling conditions as well as a method of estimating efficiency in order to provide guidelines for optimal design. The model integrates several previous mathematical approaches for assessing the main processes involved (hydraulic jump HJ, impact flow, gradually varied flows). Its performance was compared with the simulations obtained from IBER, a bi-dimensional hydrodynamic model. The impact of check dam spacing (defined by the geometric factor of influence c ) on efficiency was explored. Eleven main classifications of flow regimes were identified depending on the element and level of influence. The model produced similar results when compared with IBER, but led to higher estimations of HJ and impact lengths. Total influence guaranteed maximum efficiency and HJ control defining the location of the optimal c . Geometric total influence ( c = 1) was a valid criterion for the different stages of the structures in a wide range of situations provided that hydraulic roughness conditions remained high within the gully, e.g. through revegetation. Our total influence criterion involved shorter spacing than that habitually recommended in technical manuals for restoration, but was in line with those values found in spontaneous and stable step-pools systems, which might serve as a reference for man-made interventions.

Description: Comparison of climate change signals in CMIP3 and CMIP5 multi-model ensembles and implications for Central Asian glaciers Hydrology and Earth System Sciences, 17, 3661-3677, 2013 Author(s): A. F. Lutz, W. W. Immerzeel, A. Gobiet, F. Pellicciotti, and M. F. P. Bierkens Central Asian water resources largely depend on melt water generated in the Pamir and Tien Shan mountain ranges. To estimate future water availability in this region, it is necessary to use climate projections to estimate the future glacier extent and volume. In this study, we evaluate the impact of uncertainty in climate change projections on the future glacier extent in the Amu and Syr Darya river basins. To this end we use the latest climate change projections generated for the upcoming IPCC report (CMIP5) and, for comparison, projections used in the fourth IPCC assessment (CMIP3). With these projections we force a regionalized glacier mass balance model, and estimate changes in the basins' glacier extent as a function of the glacier size distribution in the basins and projected temperature and precipitation. This glacier mass balance model is specifically developed for implementation in large scale hydrological models, where the spatial resolution does not allow for simulating individual glaciers and data scarcity is an issue. Although the CMIP5 ensemble results in greater regional warming than the CMIP3 ensemble and the range in projections for temperature as well as precipitation is wider for the CMIP5 than for the CMIP3, the spread in projections of future glacier extent in Central Asia is similar for both ensembles. This is because differences in temperature rise are small during periods of maximum melt (July–September) while differences in precipitation change are small during the period of maximum accumulation (October–February). However, the model uncertainty due to parameter uncertainty is high, and has roughly the same importance as uncertainty in the climate projections. Uncertainty about the size of the decline in glacier extent remains large, making estimates of future Central Asian glacier evolution and downstream water availability uncertain.

Description: Prediction of dissolved reactive phosphorus losses from small agricultural catchments: calibration and validation of a parsimonious model Hydrology and Earth System Sciences, 17, 3679-3693, 2013 Author(s): C. Hahn, V. Prasuhn, C. Stamm, P. Lazzarotto, M. W. H. Evangelou, and R. Schulin Eutrophication of surface waters due to diffuse phosphorus (P) losses continues to be a severe water quality problem worldwide, causing the loss of ecosystem functions of the respective water bodies. Phosphorus in runoff often originates from a small fraction of a catchment only. Targeting mitigation measures to these critical source areas (CSAs) is expected to be most efficient and cost-effective, but requires suitable tools. Here we investigated the capability of the parsimonious Rainfall-Runoff-Phosphorus (RRP) model to identify CSAs in grassland-dominated catchments based on readily available soil and topographic data. After simultaneous calibration on runoff data from four small hilly catchments on the Swiss Plateau, the model was validated on a different catchment in the same region without further calibration. The RRP model adequately simulated the discharge and dissolved reactive P (DRP) export from the validation catchment. Sensitivity analysis showed that the model predictions were robust with respect to the classification of soils into "poorly drained" and "well drained", based on the available soil map. Comparing spatial hydrological model predictions with field data from the validation catchment provided further evidence that the assumptions underlying the model are valid and that the model adequately accounts for the dominant P export processes in the target region. Thus, the parsimonious RRP model is a valuable tool that can be used to determine CSAs. Despite the considerable predictive uncertainty regarding the spatial extent of CSAs, the RRP can provide guidance for the implementation of mitigation measures. The model helps to identify those parts of a catchment where high DRP losses are expected or can be excluded with high confidence. Legacy P was predicted to be the dominant source for DRP losses and thus, in combination with hydrologic active areas, a high risk for water quality.

Description: Senstitivity of water balance components to environmental changes in a mountainous watershed: uncertainty assessment based on models comparison Hydrology and Earth System Sciences Discussions, 10, 11983-12026, 2013 Author(s): E. Morán-Tejeda, J. Zabalza, K. Rahman, A. Gago-Silva, J. I. López-Moreno, S. Vicente-Serrano, A. Lehmann, C. L. Tague, and M. Beniston This paper evaluates the response of stream flow and other components of the water balance to changes in climate and land-use in a Pyrenean watershed. It further provides a measure of uncertainty in water resources forecasts by comparing the performance of two hydrological models: Soil and Water Assessment Tool (SWAT) and Regional Hydro-Ecological Simulation System (RHESSys). Regional Climate Model outputs for the 2021–2050 time-frame, and hypothetical (but plausible) land-use scenarios considering re-vegetation and wildfire processes were used as inputs to the models. Results indicate an overall decrease in river flows when the scenarios are considered, except for the post-fire vegetation scenario, in which stream flows are simulated to increase. However the magnitude of these projections varies between the two models used, as SWAT tends to produce larger hydrological changes under climate change scenarios, and RHESSys shows more sensitivity to changes in land-cover. The final prediction will therefore depend largely on the combination of the land-use and climate scenarios, and on the model utilized.

Description: Ice volume estimates for the Himalaya–Karakoram region: evaluating different methods The Cryosphere Discussions, 7, 4813-4854, 2013 Author(s): H. Frey, H. Machguth, M. Huss, C. Huggel, S. Bajracharya, T. Bolch, A. Kulkarni, A. Linsbauer, N. Salzmann, and M. Stoffel Ice volume estimates are crucial for assessing water reserves stored in glaciers. A variety of different methodologies exist but there is a lack of systematic comparative analysis thereof. Due to its large glacier coverage, such estimates are of particular interest for the Himalayan-Karakoram (HK) region. Here, three volume–area (V–A) relations, a slope-dependent estimation method, and two ice-thickness distribution models are applied to a complete glacier inventory of the HK region. An uncertainty and sensitivity assessment is performed to investigate the influence of the input glacier areas, and model approaches and parameters on the resulting total ice volumes. Results of the two ice-thickness distribution models are validated with local ice-thickness measurements at six glaciers. The resulting ice volumes for the entire HK region range from 2955 km 3 to 6455 km 3 , depending on the approach. Results from the ice thickness distribution models and the slope-dependent thickness estimations agree well with measured local ice thicknesses while V–A relations show stronger deviations. The study provides evidence on the significant effect of the selected method on results and underlines the importance of a careful and critical evaluation. More ice-thickness measurements are needed to improve models and results in the future.

Description: Sensitivity of simulated climate to latitudinal distribution of solar insolation reduction in SRM geoengineering methods Atmospheric Chemistry and Physics Discussions, 13, 25387-25415, 2013 Author(s): A. Modak and G. Bala Solar radiation management (SRM) geoengineering has been proposed as a potential option to counteract climate change. We perform a set of idealized geoengineering simulations to understand the global hydrological implications of varying the latitudinal distribution of solar insolation reduction in SRM methods. We find that for a fixed total mass of sulfate aerosols (12.6 Mt of SO 4 ), relative to a uniform distribution which mitigates changes in global mean temperature, global mean radiative forcing is larger when aerosol concentration is maximum at the poles leading to a warmer global mean climate and consequently an intensified hydrological cycle. Opposite changes are simulated when aerosol concentration is maximized in the tropics. We obtain a range of 1 K in global mean temperature and 3% in precipitation changes by varying the distribution pattern: this range is about 50% of the climate change from a doubling of CO 2 . Hence, our study demonstrates that a range of global mean climate states, determined by the global mean radiative forcing, are possible for a fixed total amount of aerosols but with differing latitudinal distribution, highlighting the need for a careful evaluation of SRM proposals.

Description: Source apportionment of PM 10 in a North-Western Europe regional urban background site (Lens, France) using Positive Matrix Factorization and including primary biogenic emissions Atmospheric Chemistry and Physics Discussions, 13, 25325-25385, 2013 Author(s): A. Waked, O. Favez, L. Y. Alleman, C. Piot, J.-E. Petit, T. Delaunay, E. Verlinden, B. Golly, J.-L. Besombes, J.-L. Jaffrezo, and E. Leoz-Garziandia In this work, the source of ambient particulate matter (PM 10 ) collected over a one year period at an urban background site in Lens (France) were determined and investigated using a~Positive Matrix Factorization receptor model (US EPA PMF v3.0). In addition, a Potential Source Contribution Function (PSCF) was performed by means of the Hysplit v4.9 model to assess prevailing geographical origins of the identified sources. A selective iteration process was followed for the qualification of the more robust and meaningful PMF solution. Components measured and used in the PMF include inorganic and organic species: soluble ionic species, trace elements, elemental carbon (EC), sugars alcohols, sugar anhydride, and organic carbon (OC). The mean PM 10 concentration measured from March 2011 to March 2012 was about 21 μg m −3 with typically OM, nitrate and sulfate contributing to most of the mass and accounting respectively for 5.8, 4.5 and 2.3 μg m −3 on a yearly basis. Accordingly, PMF outputs showed that the main emission sources were (in a decreasing order of contribution): secondary inorganic aerosols (28% of the total PM 10 mass), aged marine emissions (19%), with probably predominant contribution of shipping activities, biomass burning (13%), mineral dust (13%), primary biogenic emissions (9%), fresh sea salts (8%), primary traffic emissions (6%) and heavy oil combustion (4%). Significant temporal variations were observed for most of the identified sources. In particular, biomass burning emissions were negligible in summer but responsible for about 25% of total PM 10 and 50% of total OC at wintertime. Conversely, primary biogenic emissions were found to be negligible in winter but to represent about 20% of total PM 10 and 40% of total OC in summer. The latter result calls for more investigations of primary biogenic aerosols using source apportionment studies, which quite usually disregards this type of sources. This study furthermore underlines the major influence of secondary processes during daily threshold exceedances. Finally, apparent discrepancies that could be generally observed between filter-based studies (such as the present one) and Aerosol Mass Spectrometer-based PMF analyses (organic fractions) are also discussed here.

Description: Spectral Aerosol Extinction Monitoring System (SÆMS): setup, observational products, and comparisons Atmospheric Measurement Techniques Discussions, 6, 8647-8677, 2013 Author(s): A. Skupin, A. Ansmann, R. Engelmann, and H. Baars A Spectral Aerosol Extinction Monitoring System (SÆMS) is presented that allows us to continuously measure the spectral extinction coefficient of atmospheric aerosol particles along an about 2.7 km long optical path at 30–50 m height above ground at Leipzig (51.3° N, 12.4° E), Germany. The fully automated instrument measures the ambient aerosol extinction coefficients from 300–1000 nm. The main goal of SÆMS observations are long-term studies of the relationship between particle extinction and relative humidity from below 40 % to almost 100 %. The setup is presented and observations (a case study and statistical results for 2009) are discussed in terms of time series of 550 nm particle optical depth, Ångström exponent, and particle size distribution retrieved from the spectrally resolved extinction. The SÆMS measurements are compared with simultaneously performed EARLINET lidar, AERONET photometer, and in situ aerosol observations of particle size distribution and related extinction coefficients at the roof of our institute. Consistency between the different measurements is found which corroborates the quality of the SÆMS observations.

Description: Wind extraction potential from 4D-Var assimilation of O 3 , N 2 O, and H 2 O using a global shallow water model Atmospheric Chemistry and Physics Discussions, 13, 25291-25323, 2013 Author(s): D. R. Allen, K. W. Hoppel, and D. D. Kuhl The wind extraction due to assimilation of trace gas (tracer) data is examined using a 4D-Var data assimilation system based on the shallow water equations coupled to the tracer continuity equation. The procedure is outlined as follows. First, a Nature Run is created, simulating middle stratospheric winter conditions. Second, ozone (O 3 ), nitrous oxide (N 2 O), and water vapor (H 2 O) (treated in this study as passive tracers) are initialized using Microwave Limb Sounder (MLS) mixing ratios at 850 K potential temperature and advected by the Nature Run winds. Third, the initial dynamical conditions are perturbed by using a 6 h offset. Fourth, observations based on the simulated tracer data are then assimilated with a 4D-Var system in which the tracer and winds are coupled via the adjoint of the tracer continuity equation. Finally, the wind extraction potential (WEP) is calculated as the reduction of the Root Mean Square (RMS) vector wind error due to tracer assimilation relative to the total possible reduction from the initial perturbed conditions. For a single 6 h assimilation cycle of "perfect" tracer (unbiased and no imposed random errors), WEP values are 70% for O 3 , 49% for N 2 O and 16% for H 2 O. O 3 and N 2 O provide more wind information than H 2 O due to stronger background gradients relative to the tracer precisions. 10 day multi-cycle simulations with "perfect" tracer result in WEP of 98% for O 3 , 97% for N 2 O, and 90% for H 2 O. There is therefore sufficient information in these fields to nearly completely specify the dynamics, even without assimilation of dynamical information. For assimilation of tracer observations with realistic random noise (based on MLS precision at 10 hPa), the WEP after 10 days decreases to 57% for O 3 , 42% for N 2 O, and 28% for H 2 O. The root-mean-square (RMS) wind errors level out at ~ 1–2 m s −1 for these runs, suggesting a limit to which realistic tracers could constrain the winds, given complete global coverage. With higher observation noise levels, the WEP values decrease further, with negative WEP occurring in cases of very large errors for H 2 O, indicating that assimilation of very noisy observations may worsen the wind fields.

Description: Springtime phytoplankton dynamics in the Arctic Krossfjorden and Kongsfjorden (Spitsbergen) as a function of glacier proximity Biogeosciences Discussions, 10, 15519-15557, 2013 Author(s): A. M.-T. Piquet, W. H. van de Poll, R. J. W. Visser, C. Wiencke, H. Bolhuis, and A. G. J. Buma The hydrographic properties of the Kongsfjorden – Krossfjorden system (79° N, Spitsbergen) are affected by Atlantic water incursions as well as glacier meltwater runoff. This results in strong physical gradients (temperature, salinity and irradiance) within the fjords. Here, we tested the hypothesis that glaciers affect phytoplankton dynamics as early as the productive spring bloom period. During two campaigns in 2007 (late spring) and 2008 (early spring) we studied hydrographic characteristics and phytoplankton variability along 2 transects in both fjords, using HPLC-CHEMTAX pigment fingerprinting, molecular fingerprinting (DGGE) and sequencing of 18S rRNA genes. The sheltered inner fjord locations remained colder during spring as opposed to the outer locations. Vertical light attenuation coefficients increased from early spring onwards, at all locations, but in particular at the inner locations. During the end of spring, meltwater input had stratified surface waters throughout the fjords. The inner fjord locations were characterized by overall lower phytoplankton biomass. Furthermore HPLC-CHEMTAX data revealed that diatoms and Phaeocystis sp. were replaced by small nano- and picophytoplankton during late spring, coinciding with low nutrient availability. The innermost stations showed higher relative abundances of nano- and picophytoplankton throughout, notably of cyanophytes and cryptophytes. Molecular fingerprinting revealed a high similarity between inner fjord samples from early spring and late spring samples from all locations, while outer samples from early spring clustered separately. We conclude that glacier influence, mediated by early meltwater input, modifies phytoplankton biomass and composition already during the spring bloom period, in favor of low biomass and small cell size communities. This may affect higher trophic levels especially when regional warming further increases the period and volume of meltwater.

Description: Validation of SCIAMACHY O 2 A band cloud heights using Cloudnet radar/lidar measurements Atmospheric Measurement Techniques Discussions, 6, 8603-8645, 2013 Author(s): P. Wang and P. Stammes For the first time two SCIAMACHY O 2 A band cloud height products are validated using ground-based radar/lidar measurements between January 2003 and December 2011. The products are the ESA Level 2 (L2) version 5.02 cloud top height and the FRESCO (Fast Retrieval Scheme for Clouds from the Oxygen A band) version 6 cloud height. The radar/lidar profiles are obtained at the Cloudnet sites of Cabauw and Lindenberg, and are averaged for one hour centered at the SCIAMACHY overpass time to achieve an optimal temporal and spatial match. In total we have about 220 cases of single layer clouds and 200 cases of multi-layer clouds. The FRESCO cloud height and ESA L2 cloud top height are compared with the Cloudnet cloud top height and Cloudnet cloud middle height. We find that the ESA L2 cloud top height has a better agreement with the Cloudnet cloud top height than the Cloudnet cloud middle height. The ESA L2 cloud top height is on average 0.44 km higher than the Cloudnet cloud top height, with a standard deviation of 3.07 km. The FRESCO cloud height is closer to the Cloudnet cloud middle height than the Cloudnet cloud top height. The mean difference between the FRESCO cloud height and the Cloudnet cloud middle height is −0.14 km with a standard deviation of 1.88 km. The SCIAMACHY cloud height products are further compared to the Cloudnet cloud top height and the Cloudnet cloud middle height in 1 km bins. For single layer clouds, the difference between the ESA L2 cloud top height and the Cloudnet cloud top height is less than 1 km for each cloud bin at 3–7 km, which is 24 % percent of the data. The difference between the FRESCO cloud height and the Cloudnet cloud middle height is less than 1 km for each cloud bin at 0–6 km, which is 85 % percent of the data. The results are similar for multi-layer clouds, but the percentage of cases having a bias within 1 km is smaller than for single layer clouds. Since globally about 60 % of all clouds are low clouds and 42 % are single-layer low clouds, we expect that globally for a large percentage of cases the FRESCO cloud height would be close to the cloud middle height.

Description: Simulating microbial degradation of organic matter in a simple porous system using the 3-D diffusion based model MOSAIC Biogeosciences Discussions, 10, 15613-15640, 2013 Author(s): O. Monga, P. Garnier, V. Pot, E. Coucheney, N. Nunan, W. Otten, and C. Chenu This paper deals with the simulation of microbial degradation in soil within pore space at microscopic scale. Pore space was described using sphere network coming from a geometrical modeling algorithm. The biological model was improved regarding previous work in order to include transformation of dissolved organic compounds and diffusion processes. Our model was tested using experimental results of a simple substrate decomposition (Fructose) within a simple media (the sand). Diverse microbial communities were inoculated. Separated incubations in microcosms were carried out using 5 different bacterial communities at 2 different water potentials of −10 cm and −100 cm of water. We calibrated the biological parameters by means of experimental data obtained at high water content and we tested the model without any parameters change at low water content. Same as for experimental data, our simulation results showed the decrease in water content involved the decrease of mineralisation. The model was able to simulate the decrease of connectivity between substrate and microorganism due the decrease of water content.

Description: Influence of temperature and CO 2 on the strontium and magnesium composition of coccolithophore calcite Biogeosciences Discussions, 10, 15559-15586, 2013 Author(s): M. N. Müller, M. Lebrato, U. Riebesell, J. Barcelos e Ramos, K. G. Schulz, S. Blanco-Ameijeiras, S. Sett, A. Eisenhauer, and H. M. Stoll Marine calcareous sediments provide a fundamental basis for paleoceanographic studies aiming to reconstruct past oceanic conditions and understand key biogeochemical element cycles. Calcifying unicellular phytoplankton (coccolithophores) are a major contributor to both carbon and calcium cycling by photosynthesis and the production of calcite (coccoliths) in the euphotic zone and the subsequent long-term deposition and burial into marine sediments. Here we present data from controlled laboratory experiments on four coccolithophore species and elucidate the relation between the divalent cation (Sr, Mg and Ca) partitioning in coccoliths and cellular physiology (growth, calcification and photosynthesis). Coccolithophores were cultured under different seawater temperature and carbonate chemistry conditions. The partition coefficient of strontium ( D Sr ) was positively correlated with both carbon dioxide ( p CO 2 ) and temperature but displayed no coherent relation to particulate organic and inorganic carbon production rates. Furthermore, D Sr correlated positively with cellular growth rates when driven by temperature but no correlation was present when changes in growth rates were p CO 2 -induced. The results demonstrate the complex interaction between environmental forcing and physiological control on the strontium partitioning in coccolithophore calcite. The partition coefficient of magnesium ( D Mg ) displayed species-specific differences and elevated values under nutrient limitation. No conclusive correlation between coccolith D Mg and temperature was observed but p CO 2 induced a rising trend in coccolith D Mg . Interestingly, the best correlation was found between coccolith D Mg and chlorophyll a production suggesting that chlorophyll a and calcite associated Mg originate from the same intracellular pool. These results give an extended insight into the driving factors that lead to variations in the coccolith Mg / Ca ratio and can be used for Sr / Ca and Mg / Ca paleoproxy calibration.

Description: Forest NEP is significantly driven by previous year's weather Biogeosciences Discussions, 10, 15587-15611, 2013 Author(s): S. Zielis, S. Etzold, R. Zweifel, W. Eugster, M. Haeni, and N. Buchmann Understanding the response of forest net ecosystem productivity (NEP) to environmental drivers under climate change is highly relevant for predictions of annual forest carbon (C) flux budgets. Modeling annual forest NEP with soil–vegetation–atmosphere transfer models (SVATs), however, remains challenging due to unknown responses of forests to weather of the previous year. In this study, we addressed the influence of previous year's weather on the inter-annual variability of NEP for a subalpine spruce forest in Switzerland. Analysis of long-term (1997–2011) eddy covariance measurements showed that the Norway spruce forest Davos Seehornwald was a consistent sink for atmospheric CO 2 , sequestering 210 ± 88 g C m −2 per year on average. Previous year's weather strongly affected inter-annual variability of NEP, increasing the explained variance in linear models to 53% compared to 20% without previous year's weather. Thus, our results highlight the need to consider previous year's weather in modeling annual C budgets of forests. Furthermore, soil temperature in the current year's spring played a major role controlling annual NEP, mainly by influencing gross primary productivity early in the year, with spring NEP accounting for 56% of annual NEP. Consequently, we expect an increase in net CO 2 uptake with future climate warming, as long as no other resources become limiting.

Description: Black carbon concentrations from a Tibetan Plateau ice core spanning 1843–1982: recent increases due to emissions and glacier melt The Cryosphere Discussions, 7, 4855-4880, 2013 Author(s): M. Jenkins, S. Kaspari, S. Kang, B. Grigholm, and P. A. Mayewski Black carbon (BC) deposited on snow and glacier surfaces can reduce albedo and lead to accelerated melt. An ice core recovered from Guoqu glacier on Mt. Geladaindong and analyzed using a Single Particle Soot Photometer provides the first long-term (1843–1982) record of BC concentrations from the Central Tibetan Plateau. The highest concentrations are observed from 1975–1982, which corresponds to a 2.0-fold and 2.4-fold increase in average and median values, respectively, relative to 1843–1940. BC concentrations post-1940 are also elevated relative to the earlier portion of the record. Causes for the higher BC concentrations include increased regional BC emissions and subsequent deposition, and melt induced enrichment of BC, with the melt potentially accelerated due to the presence of BC at the glacier surface. A qualitative comparison of the BC and Fe (used as a dust proxy) records suggests that if changes in the concentrations of absorbing impurities at the glacier surface have influenced recent glacial melt, the melt may be due to the presence of BC rather than dust. Guoqu glacier has received no net ice accumulation since the 1980s, and is a potential example of a glacier where an increase in the equilibrium line altitude is exposing buried high impurity layers. That BC concentrations in the uppermost layers of the Geladaindong ice core are not substantially higher relative to deeper in the ice core suggests that some of the BC that must have been deposited on Guoqu glacier via wet or dry deposition between 1983 and 2005 has been removed from the surface of the glacier, potentially via supraglacial or englacial meltwater.

Description: Cumulated insolation: a simple explanation of Milankovitch's forcing on climate changes Climate of the Past Discussions, 9, 5553-5568, 2013 Author(s): F. Marra The occurrence of the sudden melting of the ice sheets during the glacial terminations is explained in this paper as the consequence of the combined role of the minima and the maxima of mean summer insolation on the Northern Hemisphere, providing a new contribution to understand the mechanisms ruling glacial forcing. Indeed, no satisfactory answer has been provided so far to the question why one specific maximum, after a series of consecutive maxima of insolation, has the potentiality to trigger a deglaciation. The explanation proposed in this paper accounts for a pre-conditioning factor, represented by "mild" (warmer) minimum, followed by a sufficiently warm maximum as the conditions that cause the end of a glacial cycle. These conditions are realized whenever the sum of the values of each consecutive minima and maxima ("cumulated insolation") on the curve of mean summer insolation at 65° N exceeds 742 Watt m −2 . The comparison of the succession of these cumulated insolation values with the astronomically tuned Oxygen isotopes record provides a satisfactory match with the occurrence of all the glacial terminations in the last 800 ka.

Description: Global and regional sea surface temperature trends during Marine Isotope Stage 11 Climate of the Past, 9, 2231-2252, 2013 Author(s): Y. Milker, R. Rachmayani, M. F. G. Weinkauf, M. Prange, M. Raitzsch, M. Schulz, and M. Kučera The Marine Isotope Stage (MIS) 11 (424–374 ka) was characterized by a protracted deglaciation and an unusually long climatic optimum. It remains unclear to what degree the climate development during this interglacial reflects the unusually weak orbital forcing or greenhouse gas trends. Previously, arguments about the duration and timing of the MIS11 climatic optimum and about the pace of the deglacial warming were based on a small number of key records, which appear to show regional differences. In order to obtain a global signal of climate evolution during MIS11, we compiled a database of 78 sea surface temperature (SST) records from 57 sites spanning MIS11, aligned these individually on the basis of benthic ( N = 28) or planktonic ( N = 31) stable oxygen isotope curves to a common time frame and subjected 48 of them to an empirical orthogonal function (EOF) analysis. The analysis revealed a high commonality among all records, with the principal SST trend explaining almost 49% of the variability. This trend indicates that on the global scale, the surface ocean underwent rapid deglacial warming during Termination V, in pace with carbon dioxide rise, followed by a broad SST optimum centered at ~410 kyr. The second EOF, which explained ~18% of the variability, revealed the existence of a different SST trend, characterized by a delayed onset of the temperature optimum during MIS11 at ~398 kyr, followed by a prolonged warm period lasting beyond 380 kyr. This trend is most consistently manifested in the mid-latitude North Atlantic and Mediterranean Sea and is here attributed to the strength of the Atlantic meridional overturning circulation. A sensitivity analysis indicates that these results are robust to record selection and to age-model uncertainties of up to 3–6 kyr, but more sensitive to SST seasonal attribution and SST uncertainties 〉1 °C. In order to validate the CCSM3 (Community Climate System Model, version 3) predictive potential, the annual and seasonal SST anomalies recorded in a total of 74 proxy records were compared with runs for three time slices representing orbital configuration extremes during the peak interglacial of MIS11. The modeled SST anomalies are characterized by a significantly lower variance compared to the reconstructions. Nevertheless, significant correlations between proxy and model data are found in comparisons on the seasonal basis, indicating that the model captures part of the long-term variability induced by astronomical forcing, which appears to have left a detectable signature in SST trends.

Description: Influence of the terrestrial magnetic field geometry on the cutoff rigidity of cosmic ray particles Annales Geophysicae, 31, 1637-1643, 2013 Author(s): K. Herbst, A. Kopp, and B. Heber Studies of the propagation of charged energetic particles in the Earth's magnetic field go back to Carl Størmer. In the end, his investigations finally lead to the definition of the so-called cutoff rigidity R C ; that is, the minimum momentum per charge a particle must have in order to reach a certain geographical location. Employing Monte Carlo simulations with the PLANETOCOSMICS code we investigate the correlation between the geomagnetic field structure and the cutoff rigidity. We show that the geometry of the magnetic field has a considerable influence on the resulting cutoff rigidity distribution. Furthermore, we will present a simple geometry-based parameter, δ B , which is able to reflect the location-dependent cutoff rigidity. We show that this correlation is also visible in the temporal evolution of the Earth's magnetic field, at least over the last 100 yr. Using latitude scans with neutron monitors, changes of the relative counting rates at different positions are calculated, showing small variations for, e.g., Kiel and Moscow, while large ones occur at Mexico City as well as on the British Virgin Islands.

Description: Electron acceleration at Jupiter: input from cyclotron-resonant interaction with whistler-mode chorus waves Annales Geophysicae, 31, 1619-1630, 2013 Author(s): E. E. Woodfield, R. B. Horne, S. A. Glauert, J. D. Menietti, and Y. Y. Shprits Jupiter has the most intense radiation belts of all the outer planets. It is not yet known how electrons can be accelerated to energies of 10 MeV or more. It has been suggested that cyclotron-resonant wave-particle interactions by chorus waves could accelerate electrons to a few MeV near the orbit of Io. Here we use the chorus wave intensities observed by the Galileo spacecraft to calculate the changes in electron flux as a result of pitch angle and energy diffusion. We show that, when the bandwidth of the waves and its variation with L are taken into account, pitch angle and energy diffusion due to chorus waves is a factor of 8 larger at L -shells greater than 10 than previously shown. We have used the latitudinal wave intensity profile from Galileo data to model the time evolution of the electron flux using the British Antarctic Survey Radiation Belt (BAS) model. This profile confines intense chorus waves near the magnetic equator with a peak intensity at ∼5° latitude. Electron fluxes in the BAS model increase by an order of magnitude for energies around 3 MeV. Extending our results to L = 14 shows that cyclotron-resonant interactions with chorus waves are equally important for electron acceleration beyond L = 10. These results suggest that there is significant electron acceleration by cyclotron-resonant interactions at Jupiter contributing to the creation of Jupiter's radiation belts and also increasing the range of L -shells over which this mechanism should be considered.

Description: The coupling between the solar wind and proton fluxes at GEO Annales Geophysicae, 31, 1631-1636, 2013 Author(s): R. J. Boynton, S. A. Billings, O. A. Amariutei, and I. Moiseenko The relationship between the solar wind and the proton flux at geosynchronous Earth orbit (GEO) is investigated using the error reduction ratio (ERR) analysis. The ERR analysis is able to search for the most appropriate inputs that control the evolution of the system. This approach is a black box method and is able to derive a mathematical model of a system from input-output data. This method is used to analyse eight energy ranges of the proton flux at GEO from 80 keV to 14.5 MeV. The inputs to the algorithm were solar wind velocity, density and pressure; the Dst index; the solar energetic proton (SEP) flux; and a function of the interplanetary magnetic field (IMF) tangential magnitude and clock angle. The results show that for lowest five energy channels (80 to 800 keV) the GEO proton fluxes are controlled by the solar wind velocity with a lag of two to three days. However, above 350 keV, the SEP fluxes, accounts for a significant portion of the GEO proton flux variance. For the highest three energy channels (0.74 to 14.5 MeV), the SEPs account for the majority of the ERR. The results also show an anisotropy of protons with gyrocenters inside GEO and outside GEO, where the protons inside GEO are controlled partly by the Dst index and also an IMF-clock angle function.

Description: Ensemble initialization of the oceanic component of a coupled model through bred vectors at seasonal-to-interannual time scales Geoscientific Model Development Discussions, 6, 5189-5214, 2013 Author(s): J. Baehr and R. Piontek We evaluate the ensemble spread at seasonal-to-interannual time scales for two perturbation techniques implemented into the ocean component of a coupled model: (1) lagged initial conditions as commonly used for decadal predictions, (2) bred vectors as commonly used for weather and seasonal forecasting. We show that relative to an uninitialized reference simulation the implementation for bred vectors can improve the ensemble spread compared to lagged initialization at time scales from one months up to three years. As bred vectors have so far mostly been used at short time scales, we initially focus on the implementation of the bred vectors into the ocean component. We introduce a depth-dependent vertical rescaling norm, accounting for the vertical dependence of the variability, and extending the commonly used upper-ocean rescaling norm to the full water column. We further show that it is sufficient for the (sub-surface) ocean to breed temperature and salinity (i.e., scalar quantities), and rely on the governing physics to carry the temperature and salinity perturbations to the flow field. Using these bred vectors with a rescaling interval of 12 months, we initialize hindcast simulations and compare them to hindcast simulations initialized with lagged initial conditions. We quantify the ensemble spread by analyzing Talagrand diagrams and spread-error ratios. For both temperature and salinity, the lagged initialized ensemble is particularly under-dispersive for the first few months of predictable lead time. The ensemble initialized with bred vectors improves the spread for temperature and salinity for the 0–700 m and 1000–3500 m means, compared to the lagged ensemble at lead times of several months to one year. As the lead time increases to years, the differences between the two ensemble initialization techniques becomes more difficult to discern. While the results need to be confirmed in an initialized framework, the present analysis represents a first step towards an improved ensemble generation at the transition from seasonal-to-interannual time scales, in particular at lead times up to one year.

Description: Understanding the performance of the FLake model over the African Great Lakes Geoscientific Model Development Discussions, 6, 5141-5187, 2013 Author(s): W. Thiery, A. Martynov, F. Darchambeau, J.-P. Descy, P.-D. Plisnier, L. Sushama, and N. P. M. van Lipzig The ability of the one-dimensional lake model FLake to represent the mixolimnion temperatures for tropical conditions was tested for three locations in East Africa: Lake Kivu, Lake Tanganyika's northern and southern basins. Meteorological observations from surrounding Automatic Weather Stations were corrected and used to drive FLake, whereas a comprehensive set of water temperature profiles served to evaluate the model at each site. Careful forcing data correction and model configuration allowed to reproduce the observed mixed layer seasonality at Lake Kivu and Lake Tanganyika (northern and southern basins), with correct representation of both the mixed layer depth and temperature structure. At Lake Kivu, mixolimnion temperatures predicted by FLake were found sensitive both to minimal variations in the external parameters (lake depth and water transparency) as to small changes in the meteorological driving data, in particular wind velocity. In each case, small modifications may already lead to a regime switch from the correctly represented seasonal mixed layer deepening to either completely mixed (down to the model lake bottom) or permanently stratified (from ~10 m downwards) conditions. In contrast, model temperatures are found robust close to the surface, with acceptable predictions of near-surface water temperatures even when the seasonal mixing regime is not reproduced. FLake can thus be a suitable tool to parameterize tropical lake water surface temperatures within atmospheric prediction models, but may be less appropriate, in its current form, to study complex limnological processes within tropical lakes. Furthermore, a study of different initial conditions showed that for tropical lakes lacking reliable initial data, a fully mixed, artificially warm initialisation is to be preferred, but only if the model is allowed to spin up until convergence is reached. Finally, FLake was used to attribute the seasonal mixing cycle at Lake Kivu to variations in the near-surface meteorological conditions. It was found that the annual mixing down to 60 m during the main dry season is primarily due to enhanced lake evaporation and secondarily due to the decreased incoming long wave radiation, both causing a significant heat loss from the lake surface and associated mixolimnion cooling.