ALBERT

Description: Source attribution and process analysis for atmospheric mercury in eastern China simulated by CMAQ-Hg Atmospheric Chemistry and Physics, 15, 8767-8779, 2015 Author(s): J. Zhu, T. Wang, J. Bieser, and V. Matthias The contribution from different emission sources and atmospheric processes to gaseous elemental mercury (GEM), gaseous oxidized mercury (GOM), particulate bound mercury (PBM) and mercury deposition in eastern China were quantified using the Community Multi-scale Air Quality (CMAQ-Hg) modeling system run with a nested domain. Natural sources (NAT) and six categories of anthropogenic mercury sources (ANTH) including cement production (CEM), domestic life (DOM), industrial boilers (IND), metal production (MET), coal-fired power plants (PP) and traffic (TRA) were considered for source apportionment. NAT were responsible for 36.6 % of annual averaged GEM concentration, which was regarded as the most important source for GEM in spite of obvious seasonal variation. Among ANTH, the influence of MET and PP on GEM were most evident especially in winter. ANTH dominated the variations of GOM and PBM concentrations with contributions of 86.7 and 79.1 %, respectively. Among ANTH, IND were the largest contributor for GOM (57.5 %) and PBM (34.4 %) so that most mercury deposition came from IND. The effect of mercury emitted from out of China was indicated by a 〉 30 % contribution to GEM concentration and wet deposition. The contributions from nine processes – consisting of emissions (EMIS), gas-phase chemical production/loss (CHEM), horizontal advection (HADV), vertical advection (ZADV), horizontal advection (HDIF), vertical diffusion (VDIF), dry deposition (DDEP), cloud processes (CLDS) and aerosol processes (AERO) – were calculated for process analysis with their comparison in urban and non-urban regions of the Yangtze River delta (YRD). EMIS and VDIF affected surface GEM and PBM concentrations most and tended to compensate each other all the time in both urban and non-urban areas. However, DDEP was the most important removal process for GOM with 7.3 and 2.9 ng m −3 reduced in the surface of urban and non-urban areas, respectively, in 1 day. The diurnal profile variation of processes revealed the transportation of GOM from urban area to non-urban areas and the importance of CHEM/AERO in higher altitudes which partly caused diffusion of GOM downwards to non-urban areas. Most of the anthropogenic mercury was transported and diffused away from urban areas by HADV and VDIF and increased mercury concentrations in non-urban areas by HADV. Natural emissions only influenced CHEM and AERO more significantly than anthropogenic. Local emissions in the YRD contributed 8.5 % more to GEM and ~ 30 % more to GOM and PBM in urban areas compared to non-urban areas.

Description: Seasonal variation of secondary organic aerosol tracers in Central Tibetan Plateau Atmospheric Chemistry and Physics, 15, 8781-8793, 2015 Author(s): R.-Q. Shen, X. Ding, Q.-F. He, Z.-Y. Cong, Q.-Q. Yu, and X.-M. Wang Secondary organic aerosol (SOA) affects the earth's radiation balance and global climate. High-elevation areas are sensitive to global climate change. However, at present, SOA origins and seasonal variations are understudied in remote high-elevation areas. In this study, particulate samples were collected from July 2012 to July 2013 at the remote Nam Co (NC) site, Central Tibetan Plateau and analyzed for SOA tracers from biogenic (isoprene, monoterpenes and β-caryophyllene) and anthropogenic (aromatics) precursors. Among these compounds, isoprene SOA (SOA I ) tracers represented the majority (26.6 ± 44.2 ng m −3 ), followed by monoterpene SOA (SOA M ) tracers (0.97 ± 0.57 ng m −3 ), aromatic SOA (SOA A ) tracer (2,3-dihydroxy-4-oxopentanoic acid, DHOPA, 0.25 ± 0.18 ng m −3 ) and β-caryophyllene SOA tracer (β-caryophyllenic acid, 0.09 ± 0.10 ng m −3 ). SOA I tracers exhibited high concentrations in the summer and low levels in the winter. The similar temperature dependence of SOA I tracers and isoprene emission suggested that the seasonal variation of SOA I tracers at the NC site was mainly influenced by the isoprene emission. The ratio of high-NO x to low-NO x products of SOA I (2-methylglyceric acid to 2-methyltetrols) was highest in the winter and lowest in the summer, due to the influence of temperature and relative humidity. The seasonal variation of SOA M tracers was impacted by monoterpenes emission and gas-particle partitioning. During the summer to the fall, temperature effect on partitioning was the dominant process influencing SOA M tracers' variation; while the temperature effect on emission was the dominant process influencing SOA M tracers' variation during the winter to the spring. SOA M tracer levels did not elevate with increased temperature in the summer, probably resulting from the counteraction of temperature effects on emission and partitioning. The concentrations of DHOPA were 1–2 orders of magnitude lower than those reported in the urban regions of the world. Due to the transport of air pollutants from the adjacent Bangladesh and northeastern India, DHOPA presented relatively higher levels in the summer. In the winter when air masses mainly came from northwestern India, mass fractions of DHOPA in total tracers increased, although its concentrations declined. The SOA-tracer method was applied to estimate secondary organic carbon (SOC) from these four precursors. The annual average of SOC was 0.22 ± 0.29 μgC m −3 , with the biogenic SOC (sum of isoprene, monoterpenes and β-caryophyllene) accounting for 75 %. In the summer, isoprene was the major precursor with its SOC contributions of 81 %. In the winter when the emission of biogenic precursors largely dropped, the contributions of aromatic SOC increased. Our study implies that anthropogenic pollutants emitted in the Indian subcontinent could be transported to the TP and have an impact on SOC over the remote NC.

Description: Influence of biomass burning on CCN number and hygroscopicity during summertime in the eastern Mediterranean Atmospheric Chemistry and Physics Discussions, 15, 21539-21582, 2015 Author(s): A. Bougiatioti, S. Bezantakos, I. Stavroulas, N. Kalivitis, P. Kokkalis, G. Biskos, N. Mihalopoulos, A. Papayannis, and A. Nenes This study investigates the CCN activity and hygroscopic properties of particles influenced by biomass burning in the eastern Mediterranean. Air masses sampled were subject to a range of atmospheric processing (several hours up to 3 days). Values of the hygroscopicity parameter, κ, were derived from cloud condensation nuclei (CCN) measurements and a Hygroscopic Tandem Differential Mobility Analyzer (HTDMA). An Aerosol Chemical Speciation Monitor (ACSM) was also used to determine the chemical composition and mass concentration of non-refractory components of the submicron aerosol fraction. During fire events, the increased organic content (and lower inorganic fraction) of the aerosol decreases the hygroscopicity parameter, κ, for all particle sizes. The reason, however, for this decrease was not the same for all size modes; smaller particle sizes appeared to be richer in less hygroscopic, less CCN-active components due to coagulation processes while larger particles become less hygroscopic during the biomass burning events due to condensation of less hygroscopic gaseous compounds. In addition, smaller particles exhibited considerable chemical dispersion (where hygroscopicity varied up to 100 % for particles of same size); larger particles, however, exhibited considerably less dispersion owing to the effects of aging and retained high levels of CCN activity. These conclusions are further supported by the observed mixing state determined by the HTDMA measurements. ACSM measurements indicate that the bulk composition reflects the hygroscopicity and chemical nature of the largest particles and a large fraction of the CCN concentrations sampled. Based on Positive Matrix Factorization (PMF) analysis of the organic ACSM spectra, CCN concentrations follow a similar trend with the BBOA component, with enhancements of CCN in biomass burning plumes ranging between 65 and 150 %, for supersaturations ranging between 0.2 and 0.7 %. Using multilinear regression, we determine the hygroscopicity of the prime organic aerosol components (BBOA, OOA-BB and OOA); it is found that the total organic hygroscopicity is very close to the inferred hygroscopicity of the oxygenated organic aerosol components. Finally, the transformation of freshly-emitted biomass burning (BBOA) to more oxidized organic aerosol (OOA-BB) can result in a two-fold increase of the inferred organic hygroscopicity. Almost 10 % of the total aerosol hygroscopicity is related to the two biomass burning components (BBOA and OOA-BB), which in turn contribute almost 35 % to the fine-particle organic water of the aerosol. This is important as organic water can contribute to the atmospheric chemistry and the direct radiative forcing.

Description: Measurement-based direct radiative effect by brown carbon over Indo-Gangetic Plain Atmospheric Chemistry and Physics Discussions, 15, 21583-21606, 2015 Author(s): A. Arola, G. L. Schuster, M. R. A. Pitkänen, O. Dubovik, H. Kokkola, A. V. Lindfors, T. Mielonen, T. Raatikainen, S. Romakkaniemi, S. N. Tripathi, and H. Lihavainen The importance of light absorbing organic aerosols, often called brown carbon (BrC), has become evident in recent years. However, there are relatively few measurement-based estimates for the direct radiative effect of BrC so far. In those earlier studies, the AErosol RObotic NETwork (AERONET) measured Aerosol Absorption Optical Depth (AAOD) and Absorption Angstrom Exponent (AAE) have been exploited. However, these two pieces of information are clearly not sufficient to separate properly carbonaceous aerosols from dust, while imaginary indices of refraction would contain more and better justified information for this purpose. This is first time that the direct radiative effect (DRE) of BrC is estimated by exploiting the AERONET-retrieved imaginary indices. We estimated it for four sites in Indo-Gangetic Plain (IGP), Karachi, Lahore, Kanpur and Gandhi College. We found a distinct seasonality, which was generally similar among all the sites, but with slightly different strengths. The monthly warming effect up to 0.5 W m -2 takes place during spring season. On the other hand, BrC results in overall cooling effect in the winter season, which can reach levels close to −1W m -2 . We then estimated similarly also DRE of black carbon and total aerosol, in order to assess the relative significance of BrC radiative effect in the radiative effects of other components. Even though BrC impact seems minor in this context, we demonstrated that it is not insignificant and moreover that it is crucial to perform spectrally resolved radiative transfer calculations to obtain good estimates for DRE of BrC.

Description: Contrail life cycle and properties from 1 year of MSG/SEVIRI rapid-scan images Atmospheric Chemistry and Physics, 15, 8739-8749, 2015 Author(s): M. Vázquez-Navarro, H. Mannstein, and S. Kox The automatic contrail tracking algorithm (ACTA) – developed to automatically follow contrails as they age, drift and spread – enables the study of a large number of contrails and the evolution of contrail properties with time. In this paper we present a year's worth of tracked contrails, from August 2008 to July 2009 in order to derive statistically significant mean values. The tracking is performed using the 5 min rapid-scan mode of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on board the Meteosat Second Generation (MSG) satellites. The detection is based on the high spatial resolution of the images provided by the Moderate Resolution Imaging Spectroradiometer on board the Terra satellite (Terra/MODIS), where a contrail detection algorithm (CDA) is applied. The results show the satellite-derived average lifetimes of contrails and contrail-cirrus along with the probability density function (PDF) of other geometric characteristics such as mean coverage, distribution and width. In combination with specifically developed algorithms (RRUMS; Rapid Retrieval of Upwelling irradiance from MSG/SEVIRI and COCS (Cirrus Optical properties derived from CALIOP and SEVIRI), explained below) it is possible to derive the radiative forcing (RF), energy forcing (EF), optical thickness (τ) and altitude of the tracked contrails. Mean values here retrieved are duration, 1 h; length, 130 km; width, 8 km; altitude, 11.7 km; optical thickness, 0.34. Radiative forcing and energy forcing are shown for land/water backgrounds in day/night situations.

Description: Accounting for the effects of sastrugi in the CERES clear-sky Antarctic shortwave angular distribution models Atmospheric Measurement Techniques, 8, 3163-3175, 2015 Author(s): J. Corbett and W. Su The Cloud and the Earth's Radiant Energy System (CERES) instruments on NASA's Terra, Aqua and Soumi NPP satellites are used to provide a long-term measurement of Earth's energy budget. To accomplish this, the radiances measured by the instruments must be inverted to fluxes by the use of a scene-type-dependent angular distribution model (ADM). For permanent snow scenes over Antarctica, shortwave (SW) ADMs are created by compositing radiance measurements over the full viewing zenith and azimuth range. However, the presence of small-scale wind blown roughness features called sastrugi cause the BRDF (bidirectional reflectance distribution function) of the snow to vary significantly based upon the solar azimuth angle and location. This can result in monthly regional biases between −12 and 7.5 Wm −2 in the inverted TOA (top-of-atmosphere) SW flux. The bias is assessed by comparing the CERES shortwave fluxes derived from nadir observations with those from all viewing zenith angles, as the sastrugi affect fluxes inverted from the oblique viewing angles more than for the nadir viewing angles. In this paper we further describe the clear-sky Antarctic ADMs from Su et al. (2015). These ADMs account for the sastrugi effect by using measurements from the Multi-Angle Imaging Spectro-Radiometer (MISR) instrument to derive statistical relationships between radiance from different viewing angles. We show here that these ADMs reduce the bias and artifacts in the CERES SW flux caused by sastrugi, both locally and Antarctic-wide. The regional monthly biases from sastrugi are reduced to between −5 and 7 Wm −2 , and the monthly-mean biases over Antarctica are reduced by up to 0.64 Wm −2 , a decrease of 74 %. These improved ADMs are used as part of the Edition 4 CERES SSF (Single Scanner Footprint) data.

Description: Space-borne observation of methane from atmospheric infrared sounder version 6: validation and implications for data analysis Atmospheric Measurement Techniques Discussions, 8, 8563-8597, 2015 Author(s): X. Xiong, F. Weng, Q. Liu, and E. Olsen Atmospheric Methane (CH 4 ) is generated as a standard product in recent version of the hyperspectral Atmospheric Infrared Sounder (AIRS-V6) aboard NASA's Aqua satellite at the NASA Goddard Earth Sciences Data and Information Services Center (NASA/GES/DISC). Significant improvements in AIRS-V6 was expected but without a thorough validation. This paper first introduced the improvements of CH 4 retrieval in AIRS-V6 and some characterizations, then presented the results of validation using ~ 1000 aircraft profiles from several campaigns spread over a couple of years and in different regions. It was found the mean biases of AIRS CH 4 at layers 343–441 and 441–575 hPa are −0.76 and −0.05 % and the RMS errors are 1.56 and 1.16 %, respectively. Further analysis demonstrates that the errors in the spring and in the high northern latitudes are larger than in other seasons or regions. The error is correlated with Degree of Freedoms (DOFs), particularly in the tropics or in the summer, and cloud amount, suggesting that the "observed" spatiotemporal variation of CH 4 by AIRS is imbedded with some artificial impact from the retrieval sensitivity in addition to its variation in reality, so the variation of information content in the retrievals needs to be taken into account in data analysis of the retrieval products. Some additional filtering (i.e. rejection of profiles with obvious oscillation as well as those deviating greatly from the norm) for quality control is recommended for the users to better utilize AIRS-V6 CH 4 , and their implementation in the future versions of the AIRS retrieval algorithm is under consideration.

Description: Impact of aerosols on the OMI tropospheric NO 2 retrievals over industrialized regions: how accurate is the aerosol correction of cloud-free scenes via a simple cloud model? Atmospheric Measurement Techniques Discussions, 8, 8385-8437, 2015 Author(s): J. Chimot, T. Vlemmix, J. P. Veefkind, J. F. de Haan, and P. F. Levelt The Ozone Monitoring Instrument (OMI) instrument has provided daily global measurements of tropospheric NO 2 for more than a decade. Numerous studies have drawn attention to the complexities related to measurements of tropospheric NO 2 in the presence of aerosols. Fine particles affect the OMI spectral measurements and the length of the average light path followed by the photons. However, they are not explicitly taken into account in the current OMI tropospheric NO 2 retrieval chain. Instead, the operational OMI O 2 -O 2 cloud retrieval algorithm is applied both to cloudy scenes and to cloud free scenes with aerosols present. This paper describes in detail the complex interplay between the spectral effects of aerosols, the OMI O 2 -O 2 cloud retrieval algorithm and the impact on the accuracy of the tropospheric NO 2 retrievals through the computed Air Mass Factor (AMF) over cloud-free scenes. Collocated OMI NO 2 and MODIS Aqua aerosol products are analysed over East China, in industrialized area. In addition, aerosol effects on the tropospheric NO 2 AMF and the retrieval of OMI cloud parameters are simulated. Both the observation-based and the simulation-based approach demonstrate that the retrieved cloud fraction linearly increases with increasing Aerosol Optical Thickness (AOT), but the magnitude of this increase depends on the aerosol properties and surface albedo. This increase is induced by the additional scattering effects of aerosols which enhance the scene brightness. The decreasing effective cloud pressure with increasing AOT represents primarily the absorbing effects of aerosols. The study cases show that the actual aerosol correction based on the implemented OMI cloud model results in biases between −20 and −40 % for the DOMINO tropospheric NO 2 product in cases of high aerosol pollution (AOT ≥ 0.6) and elevated particles. On the contrary, when aerosols are relatively close to the surface or mixed with NO 2 , aerosol correction based on the cloud model results in overestimation of the DOMINO tropospheric NO 2 product, between 10 and 20 %. These numbers are in line with comparison studies between ground-based and OMI tropospheric NO 2 measurements under conditions with high aerosol pollution and elevated particles. This highlights the need to implement an improved aerosol correction in the computation of tropospheric NO 2 AMFs.

Description: Sensitivity of thermal infrared sounders to the chemical and micro-physical properties of UTLS secondary sulphate aerosols Atmospheric Measurement Techniques Discussions, 8, 8439-8481, 2015 Author(s): P. Sellitto and B. Legras Monitoring upper tropospheric-lower stratospheric (UTLS) secondary sulphate aerosols and their chemical and micro-physical properties from satellite nadir observations is crucial to better understand their formation and evolution processes and then to estimate their impact to the UTLS chemistry, and on regional and global radiative balance. Here we present a study aimed at the evaluation of the sensitivity of thermal infrared (TIR) satellite nadir observations to the chemical composition and the size distribution of idealized UTLS sulphate aerosol layers. The extinction properties of sulphuric acid/water droplets, for different sulphuric acid mixing ratios and temperatures, are systematically analysed. The extinction coefficients are derived by means of a Mie code, using refractive indexes taken from the GEISA (Gestion et Étude des Informations Spectroscopiques Atmosphériques: Management and Study of Spectroscopic Information) spectroscopic database and log-normal size distributions with different effective radii and number concentrations. IASI (Infrared Atmospheric Sounding Interferometer) pseudo-observations are generated using forward radiative transfer calculations performed with the 4A (Automatized Atmospheric Absorption Atlas) radiative transfer model, to estimate the impact of the extinction of idealized aerosol layers, at typical UTLS conditions, on the brightness temperature spectra observed by this satellite instrument. We found a marked and typical spectral signature of these aerosol layers between 700 and 1200 cm −1 , due to the absorption bands of the sulphate and bi-sulphate ions and the undissociated sulphuric acid, with the main absorption peaks at 1170 and 905 cm −1 . The dependence of the aerosol spectral signature to the sulphuric acid mixing ratio, and effective number concentration and radius, as well as the role of interferring parameters like the ozone, sulphur dioxide, carbon dioxide and ash absorption, and temperature and water vapour profile uncertainties, are analyzed and critically discussed.

Description: Global cloud top height retrieval using SCIAMACHY limb spectra: model studies and first results Atmospheric Measurement Techniques Discussions, 8, 8295-8352, 2015 Author(s): K.-U. Eichmann, L. Lelli, C. von Savigny, H. Sembhi, and J. P. Burrows Cloud top heights (CTH) were retrieved for the period 1 January 2003 to 7 April 2012 using height-resolved limb spectra measured with the Scanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) on board ENVISAT (ENVIronmental SATellite). In this study, we tested the sensitivity of the colour index method used in the retrieval code SCODA (SCIAMACHY Cloud Detection Algorithm) and the accuracy of the retrieved CTHs in comparison to other methods. Sensitivity studies using the radiative transfer model SCIATRAN showed that the method is capable of generally detecting cloud tops down to about 5 km and very thin cirrus clouds even up to the tropopause. Volcanic particles can also be detected that occasionally reach the lower stratosphere. Low clouds at 2–3 km can only be retrieved under very clean atmospheric conditions, as light scattering of aerosols interferes with the cloud retrieval. Upper tropospheric ice clouds are detectable for cloud optical depths down to about τ N = 0.005, which is in the subvisual range. The detection sensitivity decreases towards the surface. An optical thickness of roughly 0.1 was the lower detection limit for water cloud top heights at 5 km. This value is much lower than thresholds reported for the passive cloud detection in nadir viewing direction. Comparisons with SCIAMACHY nadir cloud top heights, calculated with the Semi-Analytical CloUd Retrieval Algorithm (SACURA), showed a good agreement in the global cloud field distribution. But only opaque clouds (τ N 〉 5) are detectable with the nadir passive retrieval technique in the UV-visible and infrared wavelength range. So due to the frequent occurrence of thin and sub-visual cirrus clouds in the tropics, large cloud top height deviations were detected between both viewing geometries. Also the land/sea contrast seen in nadir retrievals was not detected in limb mode. Co-located cloud top height measurements of the limb viewing Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on ENVISAT for the period from January 2008 to March 2012 were compared, showing good agreement to within 1 km, which is smaller than the vertical field of view of both instruments. Lower stratospheric aerosols from volcanic eruptions occasionally interfered with the cloud retrieval and inhibited detection of tropospheric clouds. Examples of the impact of these events are shown for the volcanoes Kasatochi in August 2008, Sarychev Peak in June 2009, and Nabro in June 2010. Long-lasting aerosol layers were detected after these events in the Northern Hemisphere down to the tropics. Particle top heights up to about 22 km were retrieved in 2009, when the enhanced lower stratospheric aerosol layer persisted for about 7 months. Up to about 82 % of the Northern hemispheric lower stratosphere between 30° and 70° was covered by scattering particles in August 2009 and nearly half in October 2008.

Description: Operation of the Airmodus A11 nano Condensation Nucleus Counter at various inlet pressures, various operation temperatures and design of a new inlet system Atmospheric Measurement Techniques Discussions, 8, 8483-8508, 2015 Author(s): J. Kangasluoma, A. Franchin, J. Duplissy, L. Ahonen, F. Korhonen, M. Attoui, J. Mikkilä, K. Lehtipalo, J. Vanhanen, M. Kulmala, and T. Petäjä Measuring sub-3 nm particles outside of controlled laboratory conditions is a challenging task, as many of the instruments are operated at their limits and are subjected to changing ambient conditions. In this study, we advance the current understanding on the operation of Airmodus A11 nano Condensation Nucleus Counter (nCNC), which consists of a A10 Particle Size Magnifier (PSM) and A20 condensation particle counter (CPC). We explore the effect of the inlet line pressure on the measured particle concentration. We identify two different regions inside the instrument where supersaturation of working fluid can take place. We show the possibility of varying the cut-off of the instrument from 1 to 6 nm, a wider size range than the one usually covered by the PSM. We also present a new inlet system, which allows automated measurements of the background, minimizes the diffusion losses in the sampling line and is equipped with an electrostatic filter to remove ions. Finally, our view of the guidelines for optimal use of the Airmodus nCNC are provided.

Description: Looking beyond stratification: a model-based analysis of the biological drivers of oxygen depletion in the North Sea Biogeosciences Discussions, 12, 12543-12610, 2015 Author(s): F. Große, N. Greenwood, M. Kreus, H. J. Lenhart, D. Machoczek, J. Pätsch, L. A. Salt, and H. Thomas The problem of low oxygen conditions, often referred to as hypoxia, occurs regularly in the North Sea, a temperate European shelf sea. Stratification represents a major process regulating the seasonal dynamics of bottom oxygen. However, lowest oxygen conditions in the North Sea do not occur in the regions of strongest stratification. This suggests that stratification is an important prerequisite for hypoxia, but that the complex interaction between hydrodynamics and the biological processes drives its development. In this study we use the ecosystem model HAMSOM-ECOHAM5 to provide a general characteristic of the different North Sea oxygen regimes, and to quantify the impact of the different physical and biological factors driving the oxygen dynamics below the thermocline and in the bottom layer. We show that the North Sea can be subdivided into three different regimes in terms of oxygen dynamics: (1) a highly productive, non-stratified coastal regime, (2) a productive, seasonally stratified regime with a small sub-thermocline volume, and (3) a productive, seasonally stratified regime with a large sub-thermocline volume, with regime 2 being highly susceptible to hypoxic conditions. Our analysis of the different processes driving the oxygen development reveals that inter-annual variations in the oxygen conditions are caused by variations in primary production, while spatial differences can be attributed to differences in stratification and water depth. In addition, we show that benthic bacteria represent the main oxygen consumers in the bottom layer, consistently accounting for more than 50 % of the overall consumption. By providing these valuable insights, we show that ecosystem models can be a useful tool for the interpretation of observations and the estimation of the impact of anthropogenic drivers on the North Sea oxygen conditions.

Description: Impact of the oceanic geothermal heat flux on a glacial ocean state Climate of the Past Discussions, 11, 3597-3624, 2015 Author(s): M. Ballarotta, F. Roquet, S. Falahat, Q. Zhang, and G. Madec The oceanic geothermal heating (OGH) has a significant impact on the present-day ocean state, but its role during glacial periods, when the ocean circulation and stratification were different from those of today, remains poorly known. In the present study, we analyzed the response of the glacial ocean to OGH, by comparing ocean simulations of the Last Glacial Maximum (LGM, ∼ 21 ka ago) including or not geothermal heating. We found that applying the OGH warmed the Antarctic Bottom Waters (AABW) by ∼ 0.4 °C and increased the abyssal circulation by 15 to 30 % north of 30° S in the deep Pacific and Atlantic basins. The geothermally heated deep waters were then advected toward the Southern Ocean where they upwelled to the surface due to the Ekman transport. The extra heat transport towards Antarctica acted to reduce the amount of sea ice contributing to the freshening of the whole AABW overturning cell. The global amount of salt being conserved, this bottom freshening induced a salinification of the North Atlantic and North Pacific surface and intermediate waters, contributing to the deepening of the North Atlantic Deep Water. This indirect mechanism is responsible for the largest observed warming, found in the North Atlantic deep western boundary current between 2000 and 3000 m (up to 2 °C). The characteristic time scale of the ocean response to the OGH corresponds to an advective time scale (associated with the overturning of the AABW cell) rather than a diffusive time scale. The OGH might facilitate the transition from a glacial to an inter-glacial state but its effect on the deep stratification seems insufficient to drive alone an abrupt climate change.

Description: Woodland Survey of Great Britain 1971–2001 Earth System Science Data, 7, 203-214, 2015 Author(s): C. M. Wood, S. M. Smart, and R. G. H. Bunce The Woodland Survey of Great Britain is a unique data set, consisting of a detailed range of ecological measurements at a national scale, covering a time span of 30 years. A set of 103 woods spread across Britain were first surveyed in 1971, which were again surveyed in 2000–2003 (for convenience referred to subsequently as the "2001 survey"). Standardised methods of describing the trees, shrubs, ground flora, soils and general habitats present were used for both sets of surveys. The sample of 1648 plots spread through 103 woodland sites located across Britain makes it probably the most extensive quantitative ecological woodland survey undertaken in Britain; it is also notable for the range of sites that have been revisited after such a long interval. The data set provides a unique opportunity to explore the effects of a range of potential drivers of woodland change that operated between 1971 and 2001. The data set is available in four discrete parts, which have been assigned the following DOIs: doi:10.5285/4d93f9ac-68e3-49cf-8a41-4d02a7ead81a (Kirby et al., 2013b), doi:10.5285/d6409d40-58fe-4fa7-b7c8-71a105b965b4 (Kirby et al., 2013d), doi:10.5285/fb1e474d-456b-42a9-9a10-a02c35af10d2 (Kirby et al., 2013c), doi:10.5285/2d023ce9-6dbe-4b4f-a0cd-34768e1455ae (Kirby et al., 2013a).

Description: Estimating greenhouse gas emissions from travel – a GIS-based study Geographica Helvetica, 70, 185-192, 2015 Author(s): S. Kuonen Conferences, meetings and congresses are an important part of today's economic and scientific world. But the environmental impact, especially from greenhouse gas emissions associated with travel, can be extensive. Anthropogenic greenhouse gas (GHG) emissions account for the warming of the atmosphere and oceans. This study draws on the need to quantify and reduce greenhouse gas emissions associated with travel activities and aims to give suggestions for organizers and participants on possible ways to reduce greenhouse gas emissions, demonstrated on the example of the European Geography Association (EGEA) Annual Congress 2013 in Wasilkow, Poland. The lack of a comprehensive methodology for the estimation of greenhouse gas emissions from travel led to an outline of a methodology that uses geographic information systems (GIS) to calculate travel distances. The calculation of travel distances in GIS is adapted from actual transportation infrastructure, derived from the open-source platform OpenStreetMap. The methodology also aims to assess the possibilities to reduce GHG emissions by choosing different means of transportation and a more central conference location. The results of the participants of the EGEA congress, who shared their travel data for this study, show that the total travel distance adds up to 238 000 km, with average travel distance of 2429 km per participant. The travel activities of the participants in the study result in total GHG emissions of 39 300 kg CO 2 -eq including both outward and return trip. On average a participant caused GHG emissions of 401 kg CO 2 -eq. In addition, the analysis of the travel data showed differences in travel behaviour depending on the distance between conference site and point of origin. The findings on travel behaviour have then been used to give an estimation of total greenhouse gas emissions from travel for all participants of the conference, which result in a total amount of 79 711 kg CO 2 -eq. The potential for reducing greenhouse gas emissions by substituting short flights with train rides and car rides with bus and train rides is limited. Only 6 % of greenhouse gas emissions could be saved by applying these measures. Further considerable savings could only be made by substituting longer flights (32.6 %) or choosing a more central conference location (26.3 %).

Description: Singular vector based targeted observations of chemical constituents: description and first application of the EURAD-IM-SVA Geoscientific Model Development Discussions, 8, 6267-6307, 2015 Author(s): N. Goris and H. Elbern Measurements of the large dimensional chemical state of the atmosphere provide only sparse snapshots of the state of the system due to their typically insufficient temporal and spatial density. In order to optimize the measurement configurations despite those limitations, the present work describes the identification of sensitive states of the chemical system as optimal target areas for adaptive observations. For this purpose, the technique of singular vector analysis (SVA), which has been proved effective for targeted observations in numerical weather predication, is implemented into the chemical transport model EURAD-IM (EURopean Air pollution and Dispersion – Inverse Model) yielding the EURAD-IM-SVA. Besides initial values, emissions are investigated as critical simulation controlling targeting variables. For both variants, singular vectors are applied to determine the optimal placement for observations and moreover to quantify which chemical compounds have to be observed with preference. Based on measurements of the airship based ZEPTER-2 campaign, the EURAD-IM-SVA has been evaluated by conducting a comprehensive set of model runs involving different initial states and simulation lengths. Since the considered cases are restricted in terms of considered chemical compounds and selected areas, they allow for a retracing of the results and a confirmation of their correctness. Our analysis shows that the optimal placement for observations of chemical species is not entirely determined by mere transport and mixing processes. Rather, a combination of initial chemical concentrations, chemical conversions, and meteorological processes determine the influence of chemical compounds and regions. We furthermore demonstrate that the optimal placement of observations of emission strengths is highly dependent on the location of emission sources and that the benefit of including emissions as target variables outperforms the value of initial value optimisation with growing simulation length. The obtained results confirm the benefit of considering both initial values and emission strengths as target variables and of applying the EURAD-IM-SVA for measurement decision guidance with respect to chemical compounds.

Description: A global scale mechanistic model of the photosynthetic capacity Geoscientific Model Development Discussions, 8, 6217-6266, 2015 Author(s): A. A. Ali, C. Xu, A. Rogers, R. A. Fisher, S. D. Wullschleger, N. G. McDowell, E. C. Massoud, J. A. Vrugt, J. D. Muss, J. B. Fisher, P. B. Reich, and C. J. Wilson Although plant photosynthetic capacity as determined by the maximum carboxylation rate (i.e., V c, max25 ) and the maximum electron transport rate (i.e., J max25 ) at a reference temperature (generally 25 °C) is known to vary substantially in space and time in response to environmental conditions, it is typically parameterized in Earth system models (ESMs) with tabulated values associated to plant functional types. In this study, we developed a mechanistic model of leaf utilization of nitrogen for assimilation (LUNA V1.0) to predict the photosynthetic capacity at the global scale under different environmental conditions, based on the optimization of nitrogen allocated among light capture, electron transport, carboxylation, and respiration. The LUNA model was able to reasonably well capture the observed patterns of photosynthetic capacity in view that it explained approximately 55 % of the variation in observed V c, max25 and 65 % of the variation in observed J max25 across the globe. Our model simulations under current and future climate conditions indicated that V c, max25 could be most affected in high-latitude regions under a warming climate and that ESMs using a fixed V c, max25 or J max25 by plant functional types were likely to substantially overestimate future global photosynthesis.

Description: Uncertainties in calculating precipitation climatology in East Asia Hydrology and Earth System Sciences Discussions, 12, 7765-7783, 2015 Author(s): J. Kim and S. K. Park This study examines the uncertainty in calculating the fundamental climatological characteristics of precipitation in the East Asia region from multiple fine-resolution gridded analysis datasets based on in-situ rain gauge observations. Five observation-based gridded precipitation datasets are used to derive the long-term means, standard deviations in lieu of interannual variability and linear trends over the 28-year period from 1980 to 2007. Both the annual and summer (June–July–August) mean precipitation is examined. The agreement amongst these precipitation datasets are examined using multiple metrics including the signal-to-noise ratio (SNR) defined as the ratio between long-term means and the corresponding standard deviations, and Taylor diagrams which allows examinations of the pattern correlation, the standard deviation, and the centered root mean square error. It is found that the five gauge-based precipitation analysis datasets agree well in the long-term mean and interannual variability in most of the East Asia region including eastern China, Manchuria, South Korea, and Japan, which are densely populated and have fairly high density observation networks. The regions of large inter-dataset variations include Tibetan Plateau, Mongolia, northern Indo-China, and North Korea. The regions of large uncertainties are typically lightly populated and are characterized by severe terrain and/or extreme high elevations. Unlike the long-term mean and interannual variability, agreements between datasets in the linear trend is weak, both for the annual and summer mean values. In most of the East Asia region, the SNR for the linear trend is below 0.5, i.e., the inter-dataset variability exceeds the multi-data ensemble mean. The uncertainty in the spatial distribution of long-term means among these datasets occurs both in the spatial pattern and variability, but the uncertainty for the interannual variability and time trend is much larger in the variability than in the pattern correlation. Thus, care must be taken in using long-term trends calculated from gridded precipitation analysis data for climate studies over the East Asia region.

Description: Qualitative soil moisture assessment in semi-arid Africa – the role of experience and training on inter-rater reliability Hydrology and Earth System Sciences, 19, 3505-3516, 2015 Author(s): M. Rinderer, H. C. Komakech, D. Müller, G. L. B. Wiesenberg, and J. Seibert Soil and water management is particularly relevant in semi-arid regions to enhance agricultural productivity. During periods of water scarcity, soil moisture differences are important indicators of the soil water deficit and are traditionally used for allocating water resources among farmers of a village community. Here we present a simple, inexpensive soil wetness classification scheme based on qualitative indicators which one can see or touch on the soil surface. It incorporates the local farmers' knowledge on the best soil moisture conditions for seeding and brick making in the semi-arid environment of the study site near Arusha, Tanzania. The scheme was tested twice in 2014 with farmers, students and experts (April: 40 persons, June: 25 persons) for inter-rater reliability, bias of individuals and functional relation between qualitative and quantitative soil moisture values. During the test in April farmers assigned the same wetness class in 46 % of all cases, while students and experts agreed on about 60 % of all cases. Students who had been trained in how to apply the method gained higher inter-rater reliability than their colleagues with only a basic introduction. When repeating the test in June, participants were given improved instructions, organized in small subgroups, which resulted in a higher inter-rater reliability among farmers. In 66 % of all classifications, farmers assigned the same wetness class and the spread of class assignments was smaller. This study demonstrates that a wetness classification scheme based on qualitative indicators is a robust tool and can be applied successfully regardless of experience in crop growing and education level when an in-depth introduction and training is provided. The use of a simple and clear layout of the assessment form is important for reliable wetness class assignments.

Description: Nonlinear effects of locally heterogeneous hydraulic conductivity fields on regional stream–aquifer exchanges Hydrology and Earth System Sciences Discussions, 12, 7727-7764, 2015 Author(s): J. Zhu, C. L. Winter, and Z. Wang Computational experiments are performed to evaluate the effects of locally heterogeneous conductivity fields on regional exchanges of water between stream and aquifer systems in the Middle Heihe River Basin (MHRB) of northwestern China. The effects are found to be nonlinear in the sense that simulated discharges from aquifers to streams are systematically lower than discharges produced by a base model parameterized with relatively coarse effective conductivity. A similar, but weaker, effect is observed for stream leakage. The study is organized around three hypotheses: (H1) small-scale spatial variations of conductivity significantly affect regional exchanges of water between streams and aquifers in river basins, (H2) aggregating small-scale heterogeneities into regional effective parameters systematically biases estimates of stream–aquifer exchanges, and (H3) the biases result from slow-paths in groundwater flow that emerge due to small-scale heterogeneities. The hypotheses are evaluated by comparing stream–aquifer fluxes produced by the base model to fluxes simulated using realizations of the MHRB characterized by local (grid-scale) heterogeneity. Levels of local heterogeneity are manipulated as control variables by adjusting coefficients of variation. All models are implemented using the MODFLOW simulation environment, and the PEST tool is used to calibrate effective conductivities defined over 16 zones within the MHRB. The effective parameters are also used as expected values to develop log-normally distributed conductivity ( K ) fields on local grid scales. Stream-aquifer exchanges are simulated with K fields at both scales and then compared. Results show that the effects of small-scale heterogeneities significantly influence exchanges with simulations based on local-scale heterogeneities always producing discharges that are less than those produced by the base model. Although aquifer heterogeneities are uncorrelated at local scales, they appear to induce coherent slow-paths in groundwater fluxes that in turn reduce aquifer–stream exchanges. Since surface water–groundwater exchanges are critical hydrologic processes in basin-scale water budgets, these results also have implications for water resources management.

Description: Near–surface air temperature and snow skin temperature comparison from CREST-SAFE station data with MODIS land surface temperature data Hydrology and Earth System Sciences Discussions, 12, 7665-7687, 2015 Author(s): C. L. Pérez Díaz, T. Lakhankar, P. Romanov, J. Muñoz, R. Khanbilvardi, and Y. Yu Land Surface Temperature (LST) is a key variable (commonly studied to understand the hydrological cycle) that helps drive the energy balance and water exchange between the Earth's surface and its atmosphere. One observable constituent of much importance in the land surface water balance model is snow. Snow cover plays a critical role in the regional to global scale hydrological cycle because rain-on-snow with warm air temperatures accelerates rapid snow-melt, which is responsible for the majority of the spring floods. Accurate information on near-surface air temperature ( T -air) and snow skin temperature ( T -skin) helps us comprehend the energy and water balances in the Earth's hydrological cycle. T -skin is critical in estimating latent and sensible heat fluxes over snow covered areas because incoming and outgoing radiation fluxes from the snow mass and the air temperature above make it different from the average snowpack temperature. This study investigates the correlation between MODerate resolution Imaging Spectroradiometer (MODIS) LST data and observed T -air and T -skin data from NOAA-CREST-Snow Analysis and Field Experiment (CREST-SAFE) for the winters of 2013 and 2014. LST satellite validation is imperative because high-latitude regions are significantly affected by climate warming and there is a need to aid existing meteorological station networks with the spatially continuous measurements provided by satellites. Results indicate that near-surface air temperature correlates better than snow skin temperature with MODIS LST data. Additional findings show that there is a negative trend demonstrating that the air minus snow skin temperature difference is inversely proportional to cloud cover. To a lesser extent, it will be examined whether the surface properties at the site are representative for the LST properties within the instrument field of view.

Description: Does drought alter hydrological functions in forest soils? An infiltration experiment Hydrology and Earth System Sciences Discussions, 12, 7689-7725, 2015 Author(s): K. F. Gimbel, H. Puhlmann, and M. Weiler The water cycle is expected to change in future and severely affect precipitation patterns across central Europe and in other parts of the world, leading to more frequent and severe droughts. Usually, it is assumed that system properties, like soil properties, remain stable and will not be affected by drought events. To study if this assumption is appropriate, we address the effects of drought on the infiltration behavior of forest soils using dye tracer experiments on six sites in three regions across Germany, which were forced into drought conditions. The sites cover clayey, loamy and sandy textured soils. In each region, we compared a deciduous and a coniferous forest stand to address differences between the main tree species. The results of the dye tracer experiments show clear evidence for changes in infiltration behavior at the sites. The infiltration changed at the clayey plots from regular and homogeneous flow to fast preferential flow. Similar behavior was observed at the loamy plots, where large areas in the upper layers remained dry, displaying signs of strong water repellency. This was confirmed by WDPT tests, which revealed, in all except one plot, moderate to severe water repellency. Water repellency was also accountable for the change of regular infiltration to fingered flow in the sandy soils. The results of this study suggest that the "drought-history" or generally the climatic conditions in the past of a soil are more important than the actual antecedent soil moisture status regarding hydrophobicity and infiltration behavior; and also, that drought effects on infiltration need to be considered in hydrological models to obtain realistic predictions concerning water quality and quantity in runoff and groundwater recharge.

Description: State determination of catalytic converters based on an ultra-wideband communication system Journal of Sensors and Sensor Systems, 4, 255-262, 2015 Author(s): I. Motroniuk, R. Stöber, and G. Fischerauer A novel microwave-based approach for monitoring the state of aftertreatment systems such as diesel particulate filters (DPFs), three-way catalytic converters (TWCs), and selective catalytic reduction (SCR) catalysts is proposed. The volume inside the metallic housing of the DPF, TWC, or SCR is considered as a wireless communication channel between two terminals of a communication system. It is shown that, depending on the transmission channel characteristics, the properties of the catalyst, such as the catalyst state, can be inferred. This is done by means of an ultra-wideband (UWB) measurement and the subsequent evaluation and processing of the waveform in the time and frequency domains.

Description: A comparison of chemical mechanisms using tagged ozone production potential (TOPP) analysis Atmospheric Chemistry and Physics, 15, 8795-8808, 2015 Author(s): J. Coates and T. M. Butler Ground-level ozone is a secondary pollutant produced photochemically from reactions of NO x with peroxy radicals produced during volatile organic compound (VOC) degradation. Chemical transport models use simplified representations of this complex gas-phase chemistry to predict O 3 levels and inform emission control strategies. Accurate representation of O 3 production chemistry is vital for effective prediction. In this study, VOC degradation chemistry in simplified mechanisms is compared to that in the near-explicit Master Chemical Mechanism (MCM) using a box model and by "tagging" all organic degradation products over multi-day runs, thus calculating the tagged ozone production potential (TOPP) for a selection of VOCs representative of urban air masses. Simplified mechanisms that aggregate VOC degradation products instead of aggregating emitted VOCs produce comparable amounts of O 3 from VOC degradation to the MCM. First-day TOPP values are similar across mechanisms for most VOCs, with larger discrepancies arising over the course of the model run. Aromatic and unsaturated aliphatic VOCs have the largest inter-mechanism differences on the first day, while alkanes show largest differences on the second day. Simplified mechanisms break VOCs down into smaller-sized degradation products on the first day faster than the MCM, impacting the total amount of O 3 produced on subsequent days due to secondary chemistry.

Description: Size distribution and optical properties of mineral dust aerosols transported in the western Mediterranean Atmospheric Chemistry and Physics Discussions, 15, 21607-21669, 2015 Author(s): C. Denjean, F. Cassola, A. Mazzino, S. Triquet, S. Chevaillier, N. Grand, T. Bourrianne, G. Momboisse, K. Sellegri, A. Schwarzenbock, E. Freney, M. Mallet, and P. Formenti This study presents in situ aircraft measurements of Saharan mineral dust transported over the western Mediterranean basin in June–July 2013 during the ChArMEx/ADRIMED (the Chemistry-Aerosol Mediterranean Experiment/Aerosol Direct Radiative Impact on the regional climate in the MEDiterranean region) airborne campaign. Dust events differing in terms of source region (Algeria, Tunisia and Morocco), time of tranport (1–5 days) and height of transport were sampled. Mineral dust were transported above the marine boundary layer, which conversely was dominated by pollution and marine aerosols. The dust vertical structure was extremely variable and characterized by either a single layer or a more complex and stratified structure with layers originating from different source regions. Mixing of mineral dust with pollution particles was observed depending on the height of transport of the dust layers. Dust layers carried higher concentration of pollution particles at intermediate altitude (1–3 km) than at elevated altitude (〉 3 km), resulting in scattering Angstrom exponent up to 2.2 within the intermediate altitude. However, the optical properties of the dust plumes remained practically unchanged with respect to values previously measured over source regions, regardless of the altitude. Moderate light absorption of the dust plumes was observed with values of aerosol single scattering albedo at 530 nm ranging from 0.90 to 1.00 ± 0.04. Concurrent calculations from the aerosol chemical composition revealed a negligible contribution of pollution particles to the absorption properties of the dust plumes that was due to a low contribution of refractory black carbon in regards to the fraction of dust and sulfate particles. This suggests that, even in the presence of moderate pollution, likely a persistent feature in the Mediterranean, the optical properties of the dust plumes could be assimilated to those of native dust in radiative transfer simulations, modeling studies and satellite retrievals over the Mediterranean. Measurements also showed that the coarse mode of mineral dust was conserved even after 5 days of transport in the Mediterranean, which contrasts with the gravitational depletion of large particles observed during the transport of dust plumes over the Atlantic. Simulations with the WRF mesoscale meteorological model highlighted a strong vertical turbulence within the dust layers that could prevent deposition of large particles during their atmospheric transport. This has important implications for the dust radiative effects due to surface dimming, atmospheric heating and cloud formation. The results presented here add to the observational dataset necessary for evaluating the role of mineral dust on the regional climate and rainfall patterns in the western Mediterranean basin.

Description: Identification of particulate organosulfates in three megacities at the middle and lower reaches of the Yangtze River Atmospheric Chemistry and Physics Discussions, 15, 21415-21448, 2015 Author(s): X. K. Wang, S. Rossignol, Y. Ma, L. Yao, M. Y. Wang, J. M. Chen, C. George, and L. Wang PM 2.5 filter samples have been collected in three megacities i.e., Wuhan (WH), Nanjing (NJ), and Shanghai (SH) at the middle and lower reaches of the Yangtze River, respectively. Analysis of those samples using an ultra-high performance liquid chromatography (UHPLC) coupled to an orbitrap mass spectrometer (MS) allowed detection of about two hundred particulate organosulfates (OSs), including dozens of nitrooxy-organosulfates, at each location. While aliphatic OSs represented more than 78 % of the detected OSs at the three locations, aromatic OSs were much less abundant. OSs with two to four isomers accounted for about 50 % of the total OSs on average in these megacities, and the percentage of OSs with six and more isomers in WH was more significant than those in SH and NJ. The average molecular weight, and the degrees of oxidation and saturation of OSs in the WH summer samples were greater than those in WH winter samples. In SH, the average molecular weight and the degree of oxidation of OSs in summer samples were greater than those in winter samples, but the degree of saturation was similar between the two seasons. In summer, the average molecular weight, and the degrees of oxidation and unsaturation of OSs were smallest in WH among the three cities. Between NJ and SH, the average molecular weight and the degree of saturation of OSs were close and the degree of oxidation of OSs in NJ was smaller. Kendrick mass defect diagrams and Van Krevelen diagrams indicated that the characteristics of identified OSs between in NJ and in SH shared better similarity. In addition, the identity and abundance of OSs in SH showed clear seasonal and diurnal variations. OSs in summer were more abundant than they were in winter due to stronger photochemical reactions in summer. The relative abundance of OSs at night was greater than that in the daytime and more nitrooxy-OSs existed at night, probably because of active NO 3 radical chemistry at night. In SH summer samples, OSs with 5 and 10 carbons (C 5 and C 10 ) were the most abundant, indicating the importance of isoprene and monoterpenes as precursors of OSs, whereas the relative abundances of OSs with 8, 9, and more than 14 carbons (C 8 , C 9 , and C 14+ ) were also high in SH winter samples, urging the need to further understand the precursors of OSs.

Description: Trace metal characterization of aerosol particles and cloud water during HCCT 2010 Atmospheric Chemistry and Physics, 15, 8751-8765, 2015 Author(s): K. W. Fomba, D. van Pinxteren, K. Müller, Y. Iinuma, T. Lee, J. L. Collett Jr., and H. Herrmann Trace metal characterization of bulk and size-resolved aerosol and cloud water samples were performed during the Hill Cap Cloud Thuringia (HCCT) campaign. Cloud water was collected at the top of Mt. Schmücke while aerosol samples were collected at two stations upwind and downwind of Mt. Schmücke. Fourteen trace metals including Ti, V, Fe, Mn, Co, Zn, Ni, Cu, As, Sr, Rb, Pb, Cr, and Se were investigated during four full cloud events (FCEs) that fulfilled the conditions of a continuous air mass flow through the three stations. Aerosol particle trace metal concentrations were found to be lower than those observed in the same region during previous field experiments but were within a similar range to those observed in other rural regions in Europe. Fe and Zn were the most abundant elements with concentration ranges of 0.2–111.6 and 1.1–32.1 ng m −3 , respectively. Fe, Mn, and Ti were mainly found in coarse mode aerosols while Zn, Pb, and As were mostly found in the fine mode. Correlation and enrichment factor analysis of trace metals revealed that trace metals such as Ti and Rb were mostly of crustal origin while trace metals such as Zn, Pb, As, Cr, Ni, V, and Cu were of anthropogenic origin. Trace metals such as Fe and Mn were of mixed origins including crustal and combustion sources. Trace metal cloud water concentration decreased from Ti, Mn, Cr, to Co with average concentrations of 9.18, 5.59, 5.54, and 0.46 μg L −1 , respectively. A non-uniform distribution of soluble Fe, Cu, and Mn was observed across the cloud drop sizes. Soluble Fe and Cu were found mainly in cloud droplets with diameters between 16 and 22 μm, while Mn was found mostly in larger drops greater than 22 μm. Fe(III) was the main form of soluble Fe especially in the small and larger drops with concentrations ranging from 2.2 to 37.1 μg L −1 . In contrast to other studies, Fe(II) was observed mainly in the evening hours, implying its presence was not directly related to photochemical processes. Aerosol–cloud interaction did not lead to a marked increase in soluble trace metal concentrations; rather it led to differences in the chemical composition of the aerosol due to preferential loss of aerosol particles through physical processes including cloud drop deposition to vegetative surfaces.

Description: Known and unknown unknowns: the application of ensemble techniques to uncertainty estimation in satellite remote sensing data Atmospheric Measurement Techniques Discussions, 8, 8509-8562, 2015 Author(s): A. C. Povey and R. G. Grainger This paper discusses a best-practice representation of uncertainty in satellite remote sensing data. An estimate of uncertainty is necessary to make appropriate use of the information conveyed by a measurement. Traditional error propagation quantifies the uncertainty in a measurement due to well-understood perturbations in a measurement and auxiliary data – known, quantified "unknowns". The underconstrained nature of most satellite remote sensing observations requires the use of various approximations and assumptions that produce non-linear systematic errors that are not readily assessed – known, unquantifiable "unknowns". Additional errors result from the inability to resolve all scales of variation in the measured quantity – unknown "unknowns". The latter two categories of error are dominant in underconstrained remote sensing retrievals and the difficulty of their quantification limits the utility of existing uncertainty estimates, degrading confidence in such data. This paper proposes the use of ensemble techniques to present multiple self-consistent realisations of a data set as a means of depicting unquantified uncertainties. These are generated using various systems (different algorithms or forward models) believed to be appropriate to the conditions observed. Benefiting from the experience of the climate modelling community, an ensemble provides a user with a more complete representation of the uncertainty as understood by the data producer and greater freedom to consider different realisations of the data.

Description: Improved stratospheric aerosol extinction profiles from SCIAMACHY: validation and sample results Atmospheric Measurement Techniques Discussions, 8, 8353-8383, 2015 Author(s): C. von Savigny, F. Ernst, A. Rozanov, R. Hommel, K.-U. Eichmann, V. Rozanov, J. P. Burrows, and L. W. Thomason Stratospheric aerosol extinction profiles have been retrieved from SCIAMACHY/Envisat measurements of limb-scattered solar radiation. The retrieval is an improved version of an algorithm presented earlier. The retrieved aerosol extinction profiles are compared to co-located aerosol profile measurements with the SAGE II solar occultation instrument at a wavelength of 525 nm. Comparisons were carried out with two versions of the SAGE II data set (version 6.2 and the new version 7.0). In a global average sense the SCIAMACHY and the SAGE II version 7.0 extinction profiles agree to within about 10 % for altitudes above 15 km. Larger relative differences (up to 40 %) are observed at specific latitudes and altitudes. We also find differences between the two SAGE II data versions of up to 40 % for specific latitudes and altitudes. Sample results on the latitudinal and temporal variability of stratospheric aerosol extinction and optical depth during the SCIAMACHY mission period are presented. The results indicate that a series of volcanic eruptions is responsible for the increase in stratospheric aerosol optical depth from 2002 to 2012. Above about 28 km altitude volcanic eruptions are found to have negligible impact in the period 2002 to 2012.

Description: The tropopause inversion layer in baroclinic life cycles experiments: the role of diabatic and mixing processes Atmospheric Chemistry and Physics Discussions, 15, 21495-21537, 2015 Author(s): D. Kunkel, P. Hoor, and V. Wirth Recent studies on the formation of a quasi-permanent layer of enhanced static stability above the thermal tropopause revealed the contributions of dynamical and radiative processes. Dry dynamics lead to the evolution of a tropopause inversion layer (TIL) which is, however, too weak compared to observations and thus diabatic contributions are required. In this study we aim to assess the importance of diabatic as well as mixing processes in the understanding of TIL formation at midlatitudes. The non-hydrostatic model COSMO is applied in an idealized mid-latitude channel configuration to simulate baroclinic life cycles. The effect of individual diabatic, i.e. related to humidity and radiation, and turbulent processes is studied first to estimate the additional contribution of these processes to dry dynamics. In a second step these processes are stepwise included in the model to increase the complexity and finally estimate the relative importance of each process. The results suggest that including turbulence leads to a weaker TIL than in a dry reference simulation. In contrast, the TIL evolves stronger when radiation is included but the temporal occurrence is still comparable to the reference. Using various cloud schemes in the model shows that latent heat release and consecutive increased vertical motions foster an earlier and stronger appearance of the TIL than in all other life cycles. Furthermore, updrafts moisten the upper troposphere and as such increase the radiative effect from water vapor. Particularly, this process becomes more relevant for maintaining the TIL during later stages of the life cycles. Increased convergence of the vertical wind induced by updrafts and by propagating and potentially dissipating inertia-gravity waves further contributes to the enhanced stability of the lower stratosphere. Furthermore, radiative feedback of ice clouds reaching up to the tropopause is identified to potentially further affect the strength of the TIL in the region of the cloud.

Description: Impact of model developments on present and future simulations of permafrost in a global land-surface model The Cryosphere, 9, 1505-1521, 2015 Author(s): S. E. Chadburn, E. J. Burke, R. L. H. Essery, J. Boike, M. Langer, M. Heikenfeld, P. M. Cox, and P. Friedlingstein There is a large amount of organic carbon stored in permafrost in the northern high latitudes, which may become vulnerable to microbial decomposition under future climate warming. In order to estimate this potential carbon–climate feedback it is necessary to correctly simulate the physical dynamics of permafrost within global Earth system models (ESMs) and to determine the rate at which it will thaw. Additional new processes within JULES, the land-surface scheme of the UK ESM (UKESM), include a representation of organic soils, moss and bedrock and a modification to the snow scheme; the sensitivity of permafrost to these new developments is investigated in this study. The impact of a higher vertical soil resolution and deeper soil column is also considered. Evaluation against a large group of sites shows the annual cycle of soil temperatures is approximately 25 % too large in the standard JULES version, but this error is corrected by the model improvements, in particular by deeper soil, organic soils, moss and the modified snow scheme. A comparison with active layer monitoring sites shows that the active layer is on average just over 1 m too deep in the standard model version, and this bias is reduced by 70 cm in the improved version. Increasing the soil vertical resolution allows the full range of active layer depths to be simulated; by contrast, with a poorly resolved soil at least 50 % of the permafrost area has a maximum thaw depth at the centre of the bottom soil layer. Thus all the model modifications are seen to improve the permafrost simulations. Historical permafrost area corresponds fairly well to observations in all simulations, covering an area between 14 and 19 million km 2 . Simulations under two future climate scenarios show a reduced sensitivity of permafrost degradation to temperature, with the near-surface permafrost loss per degree of warming reduced from 1.5 million km 2 °C −1 in the standard version of JULES to between 1.1 and 1.2 million km 2 °C −1 in the new model version. However, the near-surface permafrost area is still projected to approximately half by the end of the 21st century under the RCP8.5 scenario.

Description: Numerical simulations of the Cordilleran ice sheet through the last glacial cycle The Cryosphere Discussions, 9, 4147-4203, 2015 Author(s): J. Seguinot, I. Rogozhina, A. P. Stroeven, M. Margold, and J. Kleman Despite more than a century of geological observations, the Cordilleran ice sheet of North America remains poorly understood in terms of its former extent, volume and dynamics. Although geomorphological evidence is abundant, its complexity is such that whole ice-sheet reconstructions of advance and retreat patterns are lacking. Here we use a numerical ice sheet model calibrated against field-based evidence to attempt a quantitative reconstruction of the Cordilleran ice sheet history through the last glacial cycle. A series of simulations is driven by time-dependent temperature offsets from six proxy records located around the globe. Although this approach reveals large variations in model response to evolving climate forcing, all simulations produce two major glaciations during marine oxygen isotope stages 4 (61.9–56.5 ka) and 2 (23.2–16.8 ka). The timing of glaciation is better reproduced using temperature reconstructions from Greenland and Antarctic ice cores than from regional oceanic sediment cores. During most of the last glacial cycle, the modelled ice cover is discontinuous and restricted to high mountain areas. However, widespread precipitation over the Skeena Mountains favours the persistence of a central ice dome throughout the glacial cycle. It acts as a nucleation centre before the Last Glacial Maximum and hosts the last remains of Cordilleran ice until the middle Holocene (6.6–6.2 ka).

Description: Stimulated Brillouin scattering during electron gyro-harmonic heating at EISCAT Annales Geophysicae, 33, 983-990, 2015 Author(s): H. Y. Fu, W. A. Scales, P. A. Bernhardt, S. J. Briczinski, M. J. Kosch, A. Senior, M. T. Rietveld, T. K. Yeoman, and J. M. Ruohoniemi Observations of secondary radiation, stimulated electromagnetic emission (SEE), produced during ionospheric modification experiments using ground-based, high-power, high-frequency (HF) radio waves are considered. The High Frequency Active Auroral Research Program (HAARP) facility is capable of generating narrowband SEE in the form of stimulated Brillouin scatter (SBS) and stimulated ion Bernstein scatter (SIBS) in the SEE spectrum. Such narrowband SEE spectral lines have not been reported using the European Incoherent Scatter (EISCAT) heater facility before. This work reports the first EISCAT results of narrowband SEE spectra and compares them to SEE previously observed at HAARP during electron gyro-harmonic heating. An analysis of experimental SEE data shows observations of emission lines within 100 Hz of the pump frequency, interpreted as SBS, during the 2012 July EISCAT campaign. Experimental results indicate that SBS strengthens as the pump frequency approaches the third electron gyro-harmonic. Also, for different heater antenna beam angles, the CUTLASS radar backscatter induced by HF radio pumping is suppressed near electron gyro-harmonics, whereas electron temperature enhancement weakens as measured by EISCAT/UHF radar. The main features of these new narrowband EISCAT observations are generally consistent with previous SBS measurements at HAARP.

Description: Optimizing the impact of temperature on bio-hydrogen production from food waste and its derivatives under no pH control using statistical modelling Biogeosciences Discussions, 12, 12823-12850, 2015 Author(s): A. Sattar, C. Arslan, C. Ji, S. Sattar, K. Yousaf, and S. Hashim The effect of temperature on bio-hydrogen production by co-digestion of sewerage sludge with food waste and its two derivatives, i.e. noodle waste and rice waste, was investigated by statistical modelling. Experimental results showed that increasing temperature from mesophilic (37 °C) to thermophilic (55 °C) was an effective mean for increasing bio-hydrogen production from food waste and noodle waste, but it caused a negative impact on bio-hydrogen production from rice waste. The maximum cumulative bio-hydrogen production of 650 mL was obtained from noodle waste under mesophilic temperature condition. Most of the production was observed during 48 h of incubation that continued till 72 h of incubation, and a decline in pH during this interval was 4.3 and 4.4 from a starting value of 7 under mesophilic and thermophilic conditions, respectively. Most of glucose consumption was also observed during 72 h of incubation and the maximum consumption was observed during the first 24 h, which was the same duration where the maximum pH drop occurred. The maximum hydrogen yields of 82.47 mL VS −1 , 131.38 mL COD −1 , and 44.90 mL glucose −1 were obtained from mesophilic food waste, thermophilic noodle waste and mesophilic rice waste respectively. The production of volatile fatty acids increased with an increase in time and temperature from food waste and noodle waste reactors whereas it decreased with temperature in rice waste reactors. The statistical modelling returned good results with high values of coefficient of determination ( R 2 ) for each waste type when it was opted for the study of cumulative hydrogen production, glucose consumption and volatile fatty acid production. The 3-D response surface plots developed by the statistical models helped a lot in developing better understanding of the impact of temperature and incubation time.

Description: Characterization and source apportionment of organic aerosol using offline aerosol mass spectrometry Atmospheric Measurement Techniques Discussions, 8, 8599-8644, 2015 Author(s): K. R. Daellenbach, C. Bozzetti, A. Křepelová, F. Canonaco, R. Wolf, P. Zotter, P. Fermo, M. Crippa, J. G. Slowik, Y. Sosedova, Y. Zhang, R.-J. Huang, L. Poulain, S. Szidat, U. Baltensperger, A. S. H. Prévôt, and I. El Haddad Field deployments of the Aerodyne Aerosol Mass Spectrometer (AMS) have significantly advanced real-time measurements and source apportionment of non-refractory particulate matter. However, the cost and complex maintenance requirements of the AMS make impractical its deployment at sufficient sites to determine regional characteristics. Furthermore, the negligible transmission efficiency of the AMS inlet for supermicron particles significantly limits the characterization of their chemical nature and contributing sources. In this study, we utilize the AMS to characterize the water-soluble organic fingerprint of ambient particles collected onto conventional quartz filters, which are routinely sampled at many air quality sites. The method was applied to 256 particulate matter (PM) filter samples (PM 1 , PM 2.5 , PM 10 ) collected at 16 urban and rural sites during summer and winter. We show that the results obtained by the present technique compare well with those from co-located online measurements, e.g. AMS or Aerosol Chemical Speciation Monitor (ACSM). The bulk recoveries of organic aerosol (60–91 %) achieved using this technique, together with low detection limits (0.8 μg of organic aerosol on the analyzed filter fraction) allow its application to environmental samples. We will discuss the recovery variability of individual hydrocarbon, oxygen containing and other ions. The performance of such data in source apportionment is assessed in comparison to ACSM data. Recoveries of organic components related to different sources as traffic, wood burning and secondary organic aerosol are presented. This technique, while subjected to the limitations inherent to filter-based measurements (e.g. filter artifacts and limited time resolution) may be used to enhance the AMS capabilities in measuring size-fractionated, spatially-resolved long-term datasets.

Description: Statistical detection and modeling of the over-dispersion of winter storm occurrence Natural Hazards and Earth System Science, 15, 1757-1761, 2015 Author(s): M. Raschke In this communication, I improve the detection and modeling of the over-dispersion of winter storm occurrence. For this purpose, the generalized Poisson distribution and the Bayesian information criterion are introduced; the latter is used for statistical model selection. Moreover, I replace the frequently used dispersion statistics by an over-dispersion parameter which does not depend on the considered return period of storm events. These models and methods are applied in order to properly detect the over-dispersion in winter storm data for Germany, carrying out a joint estimation of the distribution models for different samples.

Description: Impacts of beaver dams on hydrologic and temperature regimes in a mountain stream Hydrology and Earth System Sciences, 19, 3541-3556, 2015 Author(s): M. Majerova, B. T. Neilson, N. M. Schmadel, J. M. Wheaton, and C. J. Snow Beaver dams affect hydrologic processes, channel complexity, and stream temperature in part by inundating riparian areas, influencing groundwater–surface water interactions, and changing fluvial processes within stream systems. We explored the impacts of beaver dams on hydrologic and temperature regimes at different spatial and temporal scales within a mountain stream in northern Utah over a 3-year period spanning pre- and post-beaver colonization. Using continuous stream discharge, stream temperature, synoptic tracer experiments, and groundwater elevation measurements, we documented pre-beaver conditions in the first year of the study. In the second year, we captured the initial effects of three beaver dams, while the third year included the effects of ten dams. After beaver colonization, reach-scale (~ 750 m in length) discharge observations showed a shift from slightly losing to gaining. However, at the smaller sub-reach scale (ranging from 56 to 185 m in length), the discharge gains and losses increased in variability due to more complex flow pathways with beaver dams forcing overland flow, increasing surface and subsurface storage, and increasing groundwater elevations. At the reach scale, temperatures were found to increase by 0.38 °C (3.8 %), which in part is explained by a 230 % increase in mean reach residence time. At the smallest, beaver dam scale (including upstream ponded area, beaver dam structure, and immediate downstream section), there were notable increases in the thermal heterogeneity where warmer and cooler niches were created. Through the quantification of hydrologic and thermal changes at different spatial and temporal scales, we document increased variability during post-beaver colonization and highlight the need to understand the impacts of beaver dams on stream ecosystems and their potential role in stream restoration.

Description: Does the Budyko curve reflect a maximum power state of hydrological systems? A backward analysis Hydrology and Earth System Sciences Discussions, 12, 7821-7842, 2015 Author(s): M. Westhoff, E. Zehe, P. Archambeau, and B. Dewals Almost all catchments plot within a small envelope around the Budyko curve. This apparent behaviour suggests that organizing principles may play a role in the evolution of catchments. In this paper we applied the thermodynamic principle of maximum power as the organizing principle. In a top-down approach we derived mathematical formulations of the relation between relative wetness and gradients driving runoff and evaporation for a simple one-box model. We did this in such a way that when the conductances are optimized with the maximum power principle, the steady state behaviour of the model leads exactly to a point on the Budyko curve. Subsequently we derived gradients that, under constant forcing, resulted in a Budyko curve following the asymptotes closely. With these gradients we explored the sensitivity of dry spells and dynamics in actual evaporation. Despite the simplicity of the model, catchment observations compare reasonably well with the Budyko curves derived with dynamics in rainfall and evaporation. This indicates that the maximum power principle may be used (i) to derive the Budyko curve and (ii) to move away from the empiricism in free parameters present in many Budyko functions. Future work should focus on better representing the boundary conditions of real catchments and eventually adding more complexity to the model.

Description: Impacts of land use change and climate variations on annual inflow into Miyun Reservoir, Beijing, China Hydrology and Earth System Sciences Discussions, 12, 7785-7819, 2015 Author(s): J. K. Zheng, G. Sun, W. H. Li, X. X. Yu, C. Zhang, Y. B. Gong, and L. H. Tu Miyun reservoir, the only surface water source for Beijing city, has experienced water supply decline in recent decades. Previous studies suggest that both land use change and climate contributes the changes of water supply in this critical watershed. However, the specific causes of the decline in Miyun reservoir are debatable in a non-stationary climate in the past four decades. The central objective of this study was to quantify the separate and collective contributions of land use change and climate variability to the decreasing inflow into Miyun reservoir during 1961–2008. Different from previous studies, this work objectively identified breakpoints by analyzing the long-term historical hydrometeorology and land cover records. To effectively study the different impacts of the climate variation and land cover change during different sub-periods, annual water balance model (AWB), climate elasticity model (CEM), and rainfall–runoff model (RRM) were employed to conduct attribution analysis synthetically. We found a significant decrease in annual streamflow ( p 〈 0.01), a significant positive trend in annual potential evapotranspiration ( p 〈 0.01), and an insignificant negative trend in annual precipitation ( p 〉 0.1) during 1961–2008. Combined with historical records, we identified two breakpoints as in 1983 and 1999 for the period 1961–2008 by the sequential Mann–Kendall Test and Double Mass Curve. Climate variability alone did not explain the decrease in inflow to Miyun reservoir. Reduction of water yield was closely related to increase in evapotranspiration rates due to the expansion of forestlands and reduction in cropland and grassland, and was likely exacerbated by increased water consumption for domestic and industrial uses in the basin. Our study found that the contribution to the observed streamflow decline from land use change fell from 64–92 % during 1984–1999 to 36–58 % during 2000–2008, whereas the contribution from climate variation climbed from 8–36 % during the 1984–1999 to 42–64 % during 2000–2008. Model uncertainty analysis further demonstrated that climate warming played a dominant role in streamflow reduction in the 2000s. We conclude that future climate change and variability will further challenge the goal of water supply of Miyun reservoir to meet water demand. A comprehensive watershed management strategy needs to consider the climate variations besides vegetation management.

Description: Strange VLF bursts in northern Scandinavia: case study of the afternoon "mushroom-like" hiss on 8 December 2013 Annales Geophysicae, 33, 991-995, 2015 Author(s): J. Manninen, N. G. Kleimenova, A. Kozlovsky, I. A. Kornilov, L. I. Gromova, Y. V. Fedorenko, and T. Turunen We investigate a non-typical very low frequency (VLF) 1–4 kHz hiss representing a sequence of separated noise bursts with a strange "mushroom-like" shape in the frequency–time domain, each one lasting several minutes. These strange afternoon VLF emissions were recorded at Kannuslehto (KAN, ϕ = 67.74° N, λ = 26.27° E; L ∼ 5.5) in northern Finland during the late recovery phase of the small magnetic storm on 8 December 2013. The left-hand (LH) polarized 2–3 kHz "mushroom caps" were clearly separated from the right-hand (RH) polarized "mushroom stems" at the frequency of about 1.8–1.9 kHz, which could match the lower ionosphere waveguide cutoff (the first transverse resonance of the Earth–ionosphere cavity). We hypothesize that this VLF burst sequence could be a result of the modulation of the VLF hiss electron–cyclotron instability from the strong Pc5 geomagnetic pulsations observed simultaneously at ground-based stations as well as in the inner magnetosphere by the Time History of Events and Macroscale Interactions during Substorms mission probe (THEMIS-E; ThE). This assumption is confirmed by a similar modulation of the intensity of the energetic (1–10 keV) electrons simultaneously observed by the same ThE spacecraft. In addition, the data of the European Incoherent Scatter Scientific Association (EISCAT) radar at Tromsø show a similar quasi-periodicity in the ratio of the Hall-to-Pedersen conductance, which may be used as a proxy for the energetic particle precipitation enhancement. Our findings suggest that this strange mushroom-like shape of the considered VLF hiss could be a combined mutual effect of the magnetospheric ULF–VLF (ultra low frequency–very low frequency) wave interaction and the ionosphere waveguide propagation.

Description: Radiative and climate impacts of a large volcanic eruption during stratospheric sulfur geoengineering Atmospheric Chemistry and Physics Discussions, 15, 21837-21881, 2015 Author(s): A. Laakso, H. Kokkola, A.-I. Partanen, U. Niemeier, C. Timmreck, K. E. J. Lehtinen, H. Hakkarainen, and H. Korhonen Both explosive volcanic eruptions, which emit sulfur dioxide into the stratosphere, and stratospheric geoengineering via sulfur injections can potentially cool the climate by increasing the amount of scattering particles in the atmosphere. Here we employ a global aerosol-climate model and an earth system model to study the radiative and climate impacts of an erupting volcano during solar radiation management (SRM). According to our simulations, the radiative impacts of an eruption and SRM are not additive: in the simulated case of concurrent eruption and SRM, the peak increase in global forcing is about 40 % lower compared to a corresponding eruption into a clean background atmosphere. In addition, the recovery of the stratospheric sulfate burden and forcing was significantly faster in the concurrent case since the sulfate particles grew larger and thus sedimented faster from the stratosphere. In our simulation where we assumed that SRM would be stopped immediately after a volcano eruption, stopping SRM decreased the overall stratospheric aerosol load. For the same reasons, a volcanic eruption during SRM lead to only about 1/3 of the peak global ensemble-mean cooling compared to an eruption under unperturbed atmospheric conditions. Furthermore, the global cooling signal was seen only for 12 months after the eruption in the former scenario compared to over 40 months in the latter. In terms of the global precipitation rate, we obtain a 36 % smaller decrease in the first year after the eruption and again a clearly faster recovery in the concurrent eruption and SRM scenario. We also found that an explosive eruption could lead to significantly different regional climate responses depending on whether it takes place during geoengineering or into an unperturbed background atmosphere. Our results imply that observations from previous large eruptions, such as Mt Pinatubo in 1991, are not directly applicable when estimating the potential consequences of a volcanic eruption during stratospheric geoengineering.

Description: New fire diurnal cycle characterizations to improve fire radiative energy assessments made from MODIS observations Atmospheric Chemistry and Physics, 15, 8831-8846, 2015 Author(s): N. Andela, J. W. Kaiser, G. R. van der Werf, and M. J. Wooster Accurate near real time fire emissions estimates are required for air quality forecasts. To date, most approaches are based on satellite-derived estimates of fire radiative power (FRP), which can be converted to fire radiative energy (FRE) which is directly related to fire emissions. Uncertainties in these FRE estimates are often substantial. This is for a large part because the most often used low-Earth orbit satellite-based instruments such as the Moderate Resolution Imaging Spectroradiometer (MODIS) have a relatively poor sampling of the usually pronounced fire diurnal cycle. In this paper we explore the spatial variation of this fire diurnal cycle and its drivers using data from the geostationary Meteosat Spinning Enhanced Visible and Infrared Imager (SEVIRI). In addition, we sampled data from the SEVIRI instrument at MODIS detection opportunities to develop two approaches to estimate hourly FRE based on MODIS active fire detections. The first approach ignored the fire diurnal cycle, assuming persistent fire activity between two MODIS observations, while the second approach combined knowledge on the climatology of the fire diurnal cycle with active fire detections to estimate hourly FRE. The full SEVIRI time series, providing full coverage of the fire diurnal cycle, were used to evaluate the results. Our study period comprised of 3 years (2010–2012), and we focused on Africa and the Mediterranean basin to avoid the use of potentially lower quality SEVIRI data obtained at very far off-nadir view angles. We found that the fire diurnal cycle varies substantially over the study region, and depends on both fuel and weather conditions. For example, more "intense" fires characterized by a fire diurnal cycle with high peak fire activity, long duration over the day, and with nighttime fire activity are most common in areas of large fire size (i.e., large burned area per fire event). These areas are most prevalent in relatively arid regions. Ignoring the fire diurnal cycle generally resulted in an overestimation of FRE, while including information on the climatology of the fire diurnal cycle improved FRE estimates. The approach based on knowledge of the climatology of the fire diurnal cycle also improved distribution of FRE over the day, although only when aggregating model results to coarser spatial and/or temporal scale good correlation was found with the full SEVIRI hourly reference data set. We recommend the use of regionally varying fire diurnal cycle information within the Global Fire Assimilation System (GFAS) used in the Copernicus Atmosphere Monitoring Services, which will improve FRE estimates and may allow for further reconciliation of biomass burning emission estimates from different inventories.

Description: Validation of farm-scale methane emissions using nocturnal boundary layer budgets Atmospheric Chemistry and Physics Discussions, 15, 21765-21802, 2015 Author(s): J. Stieger, I. Bamberger, N. Buchmann, and W. Eugster This study provides the first experimental validation of Swiss agricultural methane emission estimates at the farm scale. We measured CH 4 concentrations at a Swiss farmstead during two intensive field campaigns in August 2011 and July 2012 to (1) quantify the source strength of livestock methane emissions using a tethered balloon system, and (2) to validate inventory emission estimates via nocturnal boundary layer (NBL) budgets. Field measurements were performed at a distance of 150 m from the nearest farm buildings with a tethered balloon system in combination with gradient measurements at eight heights on a 10 m tower to better resolve the near-surface concentrations. Vertical profiles of air temperature, relative humidity, CH 4 concentration, wind speed and wind direction showed that the NBL was strongly influenced by local transport processes and by the valley wind system. Methane concentrations showed a pronounced time course, with highest concentrations in the second half of the night. NBL budget flux estimates were obtained via a time–space kriging approach. Main uncertainties of NBL budget flux estimates were associated with instationary atmospheric conditions and the estimate of the inversion height z i (top of volume integration). The mean NBL budget fluxes of 1.60 ± 0.31 μg CH 4 m -2 s -1 (1.40 ± 0.50 and 1.66 ± 0.20 μg CH 4 m -2 s -1 in 2011 and 2012, respectively) were in good agreement with local inventory estimates based on current livestock number and default emission factors, with 1.29 ± 0.47 and 1.74 ± 0.63 μg CH 4 m -2 s -1 for 2011 and 2012, respectively. This indicates that emission factors used for the national inventory reports are adequate, and we conclude that the NBL budget approach is a useful tool to validate emission inventory estimates.

Description: New and improved infrared absorption cross sections for dichlorodifluoromethane (CFC-12) Atmospheric Measurement Techniques, 8, 3197-3207, 2015 Author(s): J. J. Harrison Despite its widespread commercial use throughout the twentieth century, primarily in the refrigeration industry, dichlorodifluoromethane (CFC-12) is now known to have the undesirable effect of depleting stratospheric ozone. As this long-lived molecule slowly degrades in the atmosphere, monitoring its vertical concentration profile using infrared sounders on satellite platforms crucially requires accurate laboratory spectroscopic data. This work describes new high-resolution infrared absorption cross sections of dichlorodifluoromethane over the spectral range 800–1270 cm −1 , determined from spectra recorded using a high-resolution Fourier transform spectrometer (Bruker IFS 125HR) and a 26 cm pathlength cell. Spectra of dichlorodifluoromethane/dry synthetic air mixtures were recorded at resolutions between 0.01 and 0.03 cm −1 (calculated as 0.9/MOPD; MOPD = maximum optical path difference) over a range of temperatures and pressures (7.5–761 Torr and 190–294 K) appropriate for atmospheric conditions. This new cross-section dataset improves upon the one currently available in the HITRAN and GEISA databases.

Description: Seasonality of ultrafine and sub-micron aerosols and the inferences on particle formation processes Atmospheric Chemistry and Physics Discussions, 15, 21803-21835, 2015 Author(s): H. C. Cheung, C. C.-K. Chou, M.-J. Chen, W.-R. Huang, S.-H. Huang, C.-Y. Tsai, and C. S.-L. Lee The aim of this study is to investigate the seasonal variations in the physicochemical properties of atmospheric ultrafine particles (UFPs, d ≤ 100nm) and submicron particles (PM 1 , d ≤ 1 μm) in an East-Asian urban area, which are hypothesized to be affected by the interchange of summer and winter monsoons. An observation experiment was conducted at the TARO, an urban aerosol station in Taipei, Taiwan, from October 2012 to August 2013. The measurements included the mass concentration and chemical composition of UFPs and PM 1 , as well as the particle number concentration (PNC) and size distribution (PSD) with size range of 4–736 nm. The results indicate that the mass concentration of PM 1 was elevated during cold seasons with peak level of 18.5 μg m -3 in spring, whereas the highest UFPs concentration was measured in summertime with a seasonal mean of 1.62 μg m -3 . Moreover, chemical analysis revealed that the UFPs and PM 1 were characterized by distinct composition; UFPs were composed mostly of organics, whereas ammonium and sulfate were the major constituents in PM 1 . The seasonal median of total PNCs ranged from 13.9 × 10 3 cm -3 in autumn to 19.4 × 10 3 cm -3 in spring. The PSD information retrieved from the corresponding PNC measurements indicates that the nucleation mode PNC ( N 4–25 ) peaked at 11.6 × 10 3 cm -3 in winter, whereas the Aitken mode ( N 25–100 ) and accumulation mode ( N 100–736 ) exhibited summer maxima at 6.0 × 10 3 and 3.1 × 10 3 cm -3 , respectively. The shift in PSD during summertime is attributed to the enhancement in the photochemical production of condensable organic matter that, in turn, contributes to the growth of aerosol particles in the atmosphere. In addition, remarkable photochemical production of particles was observed in spring and summer seasons, which was characterized with averaged particle growth and formation rates of 4.3 ± 0.8 nm h -1 and 1.6 ± 0.8 cm -3 s -1 , respectively. The prevalence of new particle formation (NPF) in summer is suggested as a result of seasonally enhanced photochemical oxidation of SO 2 , which contributes to the production of H 2 SO 4 , and low level of PM 10 ( d ≤ 10 μm) that serves as the condensation sink. Regarding the sources of aerosol particles, correlation analysis upon the PNCs against NO x revealed that the local vehicular exhaust was the dominant contributor of the UFPs throughout a year. On the contrary, the Asian pollution outbreaks can have significant influence in the PNC of accumulation mode particles during the seasons of winter monsoons. The results of this study underline the significance of secondary organic aerosols in the seasonal variations of UFPs and the influences of continental pollution outbreaks in the downwind areas of Asian outflows.

Description: The anthropogenic contribution to atmospheric black carbon concentrations in southern Africa: a WRF-Chem modeling study Atmospheric Chemistry and Physics, 15, 8809-8830, 2015 Author(s): F. Kuik, A. Lauer, J. P. Beukes, P. G. Van Zyl, M. Josipovic, V. Vakkari, L. Laakso, and G. T. Feig South Africa has one of the largest industrialized economies in Africa. Emissions of air pollutants are particularly high in the Johannesburg-Pretoria metropolitan area, the Mpumalanga Highveld and the Vaal Triangle, resulting in local air pollution. This study presents and evaluates a setup for conducting modeling experiments over southern Africa with the Weather Research and Forecasting model including chemistry and aerosols (WRF-Chem), and analyzes the contribution of anthropogenic emissions to the total black carbon (BC) concentrations from September to December 2010. The modeled BC concentrations are compared with measurements obtained at the Welgegund station situated ca. 100 km southwest of Johannesburg. An evaluation of WRF-Chem with observational data from ground-based measurement stations, radiosondes, and satellites shows that the meteorology is modeled mostly reasonably well, but precipitation amounts are widely overestimated and the onset of the wet season is modeled approximately 1 month too early in 2010. Modeled daily mean BC concentrations show a temporal correlation of 0.66 with measurements, but the total BC concentration is underestimated in the model by up to 50 %. Sensitivity studies with anthropogenic emissions of BC and co-emitted species turned off show that anthropogenic sources can contribute up to 100 % to BC concentrations in the industrialized and urban areas, and anthropogenic BC and co-emitted species together can contribute up to 60 % to PM 1 levels. Particularly the co-emitted species contribute significantly to the aerosol optical depth (AOD). Furthermore, in areas of large-scale biomass-burning atmospheric heating rates are increased through absorption by BC up to an altitude of about 600hPa.

Description: Calibration of 3-D wind measurements on a single-engine research aircraft Atmospheric Measurement Techniques, 8, 3177-3196, 2015 Author(s): C. Mallaun, A. Giez, and R. Baumann An innovative calibration method for the wind speed measurement using a boom-mounted Rosemount model 858 AJ air velocity probe is introduced. The method is demonstrated for a sensor system installed on a medium-size research aircraft which is used for measurements in the atmospheric boundary layer. The method encounters a series of coordinated flight manoeuvres to directly estimate the aerodynamic influences on the probe and to calculate the measurement uncertainties. The introduction of a differential Global Positioning System (DGPS) combined with a high-accuracy inertial reference system (IRS) has brought major advances to airborne measurement techniques. The exact determination of geometrical height allows the use of the pressure signal as an independent parameter. Furthermore, the exact height information and the stepwise calibration process lead to maximum accuracy. The results show a measurement uncertainty for the aerodynamic influence of the dynamic and static pressures of 0.1 hPa. The applied parametrisation does not require any height dependencies or time shifts. After extensive flight tests a correction for the flow angles (attack and sideslip angles) was found, which is necessary for a successful wind calculation. A new method is demonstrated to correct for the aerodynamic influence on the sideslip angle. For the three-dimensional (3-D) wind vector (with 100 Hz resolution) a novel error propagation scheme is tested, which determines the measurement uncertainties to be 0.3 m s −1 for the horizontal and 0.2 m s −1 for the vertical wind components.

Description: Biomass burning emissions and potential air quality impacts of volatile organic compounds and other trace gases from temperate fuels common in the United States Atmospheric Chemistry and Physics Discussions, 15, 21713-21763, 2015 Author(s): J. B. Gilman, B. M. Lerner, W. C. Kuster, P. D. Goldan, C. Warneke, P. R. Veres, J. M. Roberts, J. A. de Gouw, I. R. Burling, and R. J. Yokelson A comprehensive suite of instruments was used to quantify the emissions of over 200 organic gases, including methane and volatile organic compounds (VOCs), and 9 inorganic gases from 56 laboratory burns of 18 different biomass fuel types common in the southeastern, southwestern, or northern United States. A gas chromatograph-mass spectrometer (GC-MS) provided extensive chemical detail of discrete air samples collected during a laboratory burn and was complemented by real-time measurements of organic and inorganic species via an open-path Fourier transform infrared spectrometer (OP-FTIR) and 3 different chemical ionization-mass spectrometers. These measurements were conducted in February 2009 at the U.S. Department of Agriculture's Fire Sciences Laboratory in Missoula, Montana. The relative magnitude and composition of the gases emitted varied by individual fuel type and, more broadly, by the 3 geographic fuel regions being simulated. Emission ratios relative to carbon monoxide (CO) were used to characterize the composition of gases emitted by mass; reactivity with the hydroxyl radical, OH; and potential secondary organic aerosol (SOA) precursors for the 3 different US fuel regions presented here. VOCs contributed less than 0.78 ± 0.12 % of emissions by mole and less than 0.95 ± 0.07 % of emissions by mass (on average) due to the predominance of CO 2 , CO, CH 4 , and NO x emissions; however, VOCs contributed 70–90 (±16) % to OH reactivity and were the only measured gas-phase source of SOA precursors from combustion of biomass. Over 82 % of the VOC emissions by mole were unsaturated compounds including highly reactive alkenes and aromatics and photolabile oxygenated VOCs (OVOCs) such as formaldehyde. OVOCs contributed 57–68 % of the VOC mass emitted, 42–57 % of VOC-OH reactivity, and aromatic-OVOCs such as benzenediols, phenols, and benzaldehyde were the dominant potential SOA precursors. In addition, ambient air measurements of emissions from the Fourmile Canyon Fire that affected Boulder, Colorado in September 2010 allowed us to investigate biomass burning (BB) emissions in the presence of other VOC sources (i.e., urban and biogenic emissions) and identify several promising BB markers including benzofuran, 2-furaldehyde, 2-methylfuran, furan, and benzonitrile.

Description: Stable isotope study of a new chondrichthyan fauna (Kimmeridgian, Porrentruy, Swiss Jura): an unusual freshwater-influenced isotopic composition for the hybodont shark Asteracanthus Biogeosciences Discussions, 12, 12899-12921, 2015 Author(s): L. Leuzinger, L. Kocsis, J.-P. Billon-Bruyat, S. Spezzaferri, and T. Vennemann Chondrichthyan teeth (sharks, rays and chimaeras) are mineralised in isotopic equilibrium with the surrounding water, and parameters such as water temperature and salinity can be inferred from the oxygen isotopic composition (δ 18 O p ) of their bioapatite. We analysed a new chondrichthyan assemblage, as well as teeth from bony fish (Pycnodontiformes). All specimens are from Kimmeridgian coastal marine deposits of the Swiss Jura (vicinity of Porrentruy, Ajoie district, NW Switzerland). While the overall faunal composition and the isotopic composition of bony fish are consistent with marine conditions, unusually low δ 18 O p values were measured for the hybodont shark Asteracanthus . These values are also lower compared to previously published data from older European Jurassic localities. Additional analyses on material from Solothurn (Kimmeridgian, NW Switzerland) also have comparable, low- 18 O isotopic compositions for Asteracanthus . The data are hence interpreted to represent a so far unique, freshwater-influenced isotopic composition for this shark that is classically considered as a marine genus. While reproduction in freshwater or brackish realms is established for other hybodonts, a similar behaviour for Asteracanthus is proposed here. Regular excursions into lower salinity waters can be linked to the age of the deposits and correspond to an ecological adaptation, most likely driven by the Kimmeridgian transgression and by the competition of the primitive shark Asteracanthus with the rapidly diversifying neoselachians (modern sharks).

Description: Climate change impacts on net primary production (NPP) and export production (EP) regulated by increasing stratification and phytoplankton community structure in CMIP5 models Biogeosciences Discussions, 12, 12851-12897, 2015 Author(s): W. Fu, J. Randerson, and J. K. Moore We examine climate change impacts on net primary production (NPP) and export production (sinking particulate flux; EP) with simulations from nine Earth System Models (ESMs) performed in the framework of the fifth Coupled Model Inter-comparison Project (CMIP5). Global NPP and EP are reduced considerably by the end of the century for the intense warming scenario of Representative Concentration Pathway (RCP) 8.5. Relative to the 1990s, global NPP in the 2090s is reduced by 2.3–16 % and EP by 7–18 %. The models with the largest increases in stratification (and largest relative reductions in NPP and EP) also show the largest positive biases in stratification for the contemporary period, suggesting some potential overestimation of climate impacts on NPP and EP. All of the CMIP5 models show an increase in stratification in response to surface ocean warming and freshening that is accompanied by decreases in NPP, EP, and surface macronutrient concentrations. There is considerable variability across models in the absolute magnitude of these fluxes, surface nutrient concentrations, and their perturbations by climate change, indicating large model uncertainties. The negative response of NPP and EP to stratification increases reflects a bottom-up control, as nutrient flux to the euphotic zone declines. Models with dynamic phytoplankton community structure show larger declines in EP than in NPP. This is driven by phytoplankton community composition shifts, with a reduced percentage of NPP by large phytoplankton under RCP 8.5, as smaller phytoplankton are favored under the increasing nutrient stress. Thus, projections of the NPP response to climate change in the CMIP5 models are critically dependent on the simulated phytoplankton community structure, the efficiency of the biological pump, and the resulting (highly variable) levels of regenerated production. Community composition is represented relatively simply in the CMIP5 models, and should be expanded to better capture the spatial patterns and the changes in export efficiency that are necessary for predicting climate impacts on NPP.

Description: Examining the effects of anthropogenic emissions on isoprene-derived secondary organic aerosol formation during the 2013 Southern Oxidant and Aerosol Study (SOAS) at the Look Rock, Tennessee ground site Atmospheric Chemistry and Physics, 15, 8871-8888, 2015 Author(s): S. H. Budisulistiorini, X. Li, S. T. Bairai, J. Renfro, Y. Liu, Y. J. Liu, K. A. McKinney, S. T. Martin, V. F. McNeill, H. O. T. Pye, A. Nenes, M. E. Neff, E. A. Stone, S. Mueller, C. Knote, S. L. Shaw, Z. Zhang, A. Gold, and J. D. Surratt A suite of offline and real-time gas- and particle-phase measurements was deployed at Look Rock, Tennessee (TN), during the 2013 Southern Oxidant and Aerosol Study (SOAS) to examine the effects of anthropogenic emissions on isoprene-derived secondary organic aerosol (SOA) formation. High- and low-time-resolution PM 2.5 samples were collected for analysis of known tracer compounds in isoprene-derived SOA by gas chromatography/electron ionization-mass spectrometry (GC/EI-MS) and ultra performance liquid chromatography/diode array detection-electrospray ionization-high-resolution quadrupole time-of-flight mass spectrometry (UPLC/DAD-ESI-HR-QTOFMS). Source apportionment of the organic aerosol (OA) was determined by positive matrix factorization (PMF) analysis of mass spectrometric data acquired on an Aerodyne Aerosol Chemical Speciation Monitor (ACSM). Campaign average mass concentrations of the sum of quantified isoprene-derived SOA tracers contributed to ~ 9 % (up to 28 %) of the total OA mass, with isoprene-epoxydiol (IEPOX) chemistry accounting for ~ 97 % of the quantified tracers. PMF analysis resolved a factor with a profile similar to the IEPOX-OA factor resolved in an Atlanta study and was therefore designated IEPOX-OA. This factor was strongly correlated ( r 2 〉 0.7) with 2-methyltetrols, C 5 -alkene triols, IEPOX-derived organosulfates, and dimers of organosulfates, confirming the role of IEPOX chemistry as the source. On average, IEPOX-derived SOA tracer mass was ~ 26 % (up to 49 %) of the IEPOX-OA factor mass, which accounted for 32 % of the total OA. A low-volatility oxygenated organic aerosol (LV-OOA) and an oxidized factor with a profile similar to 91Fac observed in areas where emissions are biogenic-dominated were also resolved by PMF analysis, whereas no primary organic aerosol (POA) sources could be resolved. These findings were consistent with low levels of primary pollutants, such as nitric oxide (NO ~ 0.03 ppb), carbon monoxide (CO ~ 116 ppb), and black carbon (BC ~ 0.2 μg m −3 ). Particle-phase sulfate is fairly correlated ( r 2 ~ 0.3) with both methacrylic acid epoxide (MAE)/hydroxymethyl-methyl-α-lactone (HMML)- (henceforth called methacrolein (MACR)-derived SOA tracers) and IEPOX-derived SOA tracers, and more strongly correlated ( r 2 ~ 0.6) with the IEPOX-OA factor, in sum suggesting an important role of sulfate in isoprene SOA formation. Moderate correlation between the MACR-derived SOA tracer 2-methylglyceric acid with sum of reactive and reservoir nitrogen oxides (NO y ; r 2 = 0.38) and nitrate ( r 2 = 0.45) indicates the potential influence of anthropogenic emissions through long-range transport. Despite the lack of a clear association of IEPOX-OA with locally estimated aerosol acidity and liquid water content (LWC), box model calculations of IEPOX uptake using the simpleGAMMA model, accounting for the role of acidity and aerosol water, predicted the abundance of the IEPOX-derived SOA tracers 2-methyltetrols and the corresponding sulfates with good accuracy ( r 2 ~ 0.5 and ~ 0.7, respectively). The modeling and data combined suggest an anthropogenic influence on isoprene-derived SOA formation through acid-catalyzed heterogeneous chemistry of IEPOX in the southeastern US. However, it appears that this process was not limited by aerosol acidity or LWC at Look Rock during SOAS. Future studies should further explore the extent to which acidity and LWC as well as aerosol viscosity and morphology becomes a limiting factor of IEPOX-derived SOA, and their modulation by anthropogenic emissions.

Description: Real-time measurements of secondary organic aerosol formation and aging from ambient air in an oxidation flow reactor in the Los Angeles area Atmospheric Chemistry and Physics Discussions, 15, 21907-21958, 2015 Author(s): A. M. Ortega, P. L. Hayes, Z. Peng, B. B. Palm, W. Hu, D. A. Day, R. Li, M. J. Cubison, W. H. Brune, M. Graus, C. Warneke, J. B. Gilman, W. C. Kuster, J. A. de Gouw, and J. L. Jimenez Field studies in polluted areas over the last decade have observed large formation of secondary organic aerosol (SOA) that is often poorly captured by models. The study of SOA formation using ambient data is often confounded by the effects of advection, vertical mixing, emissions, and variable degrees of photochemical aging. An Oxidation Flow Reactor (OFR) was deployed to study SOA formation in real-time during the CalNex campaign in Pasadena, CA, in 2010. A high-resolution aerosol mass spectrometer (AMS) and a scanning mobility particle sizer (SMPS) alternated sampling ambient and reactor-aged air. The reactor produced OH concentrations up to 4 orders of magnitude higher than in ambient air, achieving equivalent atmospheric aging from hours up to several weeks in 3 min of processing. OH radical concentration was continuously stepped, obtaining measurements of real-time SOA formation and oxidation at multiple equivalent ages from 0.8 days–6.4 weeks. Enhancement of OA from aging showed a maximum net SOA production between 0.8–6 days of aging with net OA mass loss beyond 2 weeks. Reactor SOA mass peaked at night, in the absence of ambient photochemistry, and correlated with trimethylbenzene concentrations. Reactor SOA formation was inversely correlated with ambient SOA and O x , which along with the short-lived VOC correlation, indicates the importance of relatively reactive (τ OH ∼ 0.3 day) SOA precursors in the LA-Basin. Evolution of the elemental composition in the reactor was similar to trends observed in the atmosphere (O : C vs. H : C slope ∼ -0.65). Oxidation state of carbon (OS C ) in reactor SOA increased steeply with age and remained elevated (OS C ∼ 2) at the highest photochemical ages probed. The ratio of OA in the reactor output to excess CO (ΔCO, ambient CO above regional background) vs. photochemical age is similar to previous studies at low to moderate ages and also extends to higher ages where OA loss dominates. The mass added at low-to-intermediate ages is due primarily to condensation of oxidized species, not heterogeneous oxidation. The OA decrease at high photochemical ages is dominated by heterogeneous oxidation followed by fragmentation/evaporation. A comparison of urban SOA formation in this study with a similar study of vehicle SOA in a tunnel supports the dominance of vehicle emissions in urban SOA. Pre-2007 SOA models underpredict SOA formation by an order of magnitude, while a more recent model performs better but overpredicts at higher ages. These results demonstrate the value of the reactor as a tool for in situ evaluation of the SOA formation potential and OA evolution from ambient air.

Description: On the consistency of 2-D video disdrometers in measuring microphysical parameters of solid precipitation Atmospheric Measurement Techniques, 8, 3251-3261, 2015 Author(s): F. Bernauer, K. Hürkamp, W. Rühm, and J. Tschiersch Detailed characterization and classification of precipitation is an important task in atmospheric research. Line scanning 2-D video disdrometer devices are well established for rain observations. The two orthogonal views taken of each hydrometeor passing the sensitive area of the instrument qualify these devices especially for detailed characterization of nonsymmetric solid hydrometeors. However, in case of solid precipitation, problems related to the matching algorithm have to be considered and the user must be aware of the limited spatial resolution when size and shape descriptors are analyzed. Clarifying the potential of 2-D video disdrometers in deriving size, velocity and shape parameters from single recorded pictures is the aim of this work. The need of implementing a matching algorithm suitable for mixed- and solid-phase precipitation is highlighted as an essential step in data evaluation. For this purpose simple reproducible experiments with solid steel spheres and irregularly shaped Styrofoam particles are conducted. Self-consistency of shape parameter measurements is tested in 38 cases of real snowfall. As a result, it was found that reliable size and shape characterization with a relative standard deviation of less than 5 % is only possible for particles larger than 1 mm. For particles between 0.5 and 1.0 mm the relative standard deviation can grow up to 22 % for the volume, 17 % for size parameters and 14 % for shape descriptors. Testing the adapted matching algorithm with a reproducible experiment with Styrofoam particles, a mismatch probability of less than 3 % was found. For shape parameter measurements in case of real solid-phase precipitation, the 2-DVD shows self-consistent behavior.

Description: Schneefernerhaus as a mountain research station for clouds and turbulence Atmospheric Measurement Techniques, 8, 3209-3218, 2015 Author(s): S. Risius, H. Xu, F. Di Lorenzo, H. Xi, H. Siebert, R. A. Shaw, and E. Bodenschatz Cloud measurements are usually carried out with airborne campaigns, which are expensive and are limited by temporal duration and weather conditions. Ground-based measurements at high-altitude research stations therefore play a complementary role in cloud study. Using the meteorological data (wind speed, direction, temperature, humidity, visibility, etc.) collected by the German Weather Service (DWD) from 2000 to 2012 and turbulence measurements recorded by multiple ultrasonic sensors (sampled at 10 Hz) in 2010, we show that the Umweltforschungsstation Schneefernerhaus (UFS) located just below the peak of Zugspitze in the German Alps, at a height of 2650 m, is a well-suited station for cloud–turbulence research. The wind at UFS is dominantly in the east–west direction and nearly horizontal. During the summertime (July and August) the UFS is immersed in warm clouds about 25 % of the time. The clouds are either from convection originating in the valley in the east, or associated with synoptic-scale weather systems typically advected from the west. Air turbulence, as measured from the second- and third-order velocity structure functions that exhibit well-developed inertial ranges, possesses Taylor microscale Reynolds numbers up to 10 4 , with the most probable value at ~ 3000. In spite of the complex topography, the turbulence appears to be nearly as isotropic as many laboratory flows when evaluated on the "Lumley triangle".

Description: Technical Note: A simple calculation algorithm to separate high-resolution CH 4 flux measurements into ebullition and diffusion-derived components Biogeosciences Discussions, 12, 12923-12945, 2015 Author(s): M. Hoffmann, M. Schulz-Hanke, J. Garcia Alba, N. Jurisch, U. Hagemann, T. Sachs, M. Sommer, and J. Augustin Processes driving the production, transformation and transport of methane (CH 4 ) in wetland ecosystems are highly complex. Thus, serious challenges are constitutes in terms of the mechanistic process understanding, the identification of potential environmental drivers and the calculation of reliable CH 4 emission estimates. We present a simple calculation algorithm to separate open-water CH 4 fluxes measured with automatic chambers into diffusion- and ebullition-derived components, which helps facilitating the identification of underlying dynamics and potential environmental drivers. Flux separation is based on ebullition related sudden concentration changes during single measurements. A variable ebullition filter is applied, using the lower and upper quartile and the interquartile range (IQR). Automation of data processing is achieved by using an established R-script, adjusted for the purpose of CH 4 flux calculation. The algorithm was tested using flux measurement data (July to September 2013) from a former fen grassland site, converted into a shallow lake as a result of rewetting ebullition and diffusion contributed 46 and 55 %, respectively, to total CH 4 emissions, which is comparable to those previously reported by literature. Moreover, the separation algorithm revealed a concealed shift in the diurnal trend of diffusive fluxes throughout the measurement period.

Description: Plume-based analysis of vehicle fleet air pollutant emissions and the contribution from high emitters Atmospheric Measurement Techniques, 8, 3263-3275, 2015 Author(s): J. M. Wang, C.-H. Jeong, N. Zimmerman, R. M. Healy, D. K. Wang, F. Ke, and G. J. Evans An automated identification and integration method has been developed for in-use vehicle emissions under real-world conditions. This technique was applied to high-time-resolution air pollutant measurements of in-use vehicle emissions performed under real-world conditions at a near-road monitoring station in Toronto, Canada, during four seasons, through month-long campaigns in 2013–2014. Based on carbon dioxide measurements, over 100 000 vehicle-related plumes were automatically identified and fuel-based emission factors for nitrogen oxides; carbon monoxide; particle number; black carbon; benzene, toluene, ethylbenzene, and xylenes (BTEX); and methanol were determined for each plume. Thus the automated identification enabled the measurement of an unprecedented number of plumes and pollutants over an extended duration. Emission factors for volatile organic compounds were also measured roadside for the first time using a proton transfer reaction time-of-flight mass spectrometer; this instrument provided the time resolution required for the plume capture technique. Mean emission factors were characteristic of the light-duty gasoline-dominated vehicle fleet present at the measurement site, with mean black carbon and particle number emission factors of 35 mg kg fuel −1 and 7.5 × 10 14 # kg fuel −1 , respectively. The use of the plume-by-plume analysis enabled isolation of vehicle emissions, and the elucidation of co-emitted pollutants from similar vehicle types, variability of emissions across the fleet, and the relative contribution from heavy emitters. It was found that a small proportion of the fleet ( 〈 25 %) contributed significantly to total fleet emissions: 100, 100, 81, and 77 % for black carbon, carbon monoxide, BTEX, and particle number, respectively. Emission factors of a single pollutant may help classify a vehicle as a high emitter; however, regulatory strategies to more efficiently target multi-pollutant mixtures may be better developed by considering the co-emitted pollutants as well.

Description: High-resolution measurement of cloud microphysics and turbulence at a mountaintop station Atmospheric Measurement Techniques, 8, 3219-3228, 2015 Author(s): H. Siebert, R. A. Shaw, J. Ditas, T. Schmeissner, S. P. Malinowski, E. Bodenschatz, and H. Xu Mountain research stations are advantageous not only for long-term sampling of cloud properties but also for measurements that are prohibitively difficult to perform on airborne platforms due to the large true air speed or adverse factors such as weight and complexity of the equipment necessary. Some cloud–turbulence measurements, especially Lagrangian in nature, fall into this category. We report results from simultaneous, high-resolution and collocated measurements of cloud microphysical and turbulence properties during several warm cloud events at the Umweltforschungsstation Schneefernerhaus (UFS) on Zugspitze in the German Alps. The data gathered were found to be representative of observations made with similar instrumentation in free clouds. The observed turbulence shared all features known for high-Reynolds-number flows: it exhibited approximately Gaussian fluctuations for all three velocity components, a clearly defined inertial subrange following Kolmogorov scaling (power spectrum, and second- and third-order Eulerian structure functions), and highly intermittent velocity gradients, as well as approximately lognormal kinetic energy dissipation rates. The clouds were observed to have liquid water contents on the order of 1 g m −3 and size distributions typical of continental clouds, sometimes exhibiting long positive tails indicative of large drop production through turbulent mixing or coalescence growth. Dimensionless parameters relevant to cloud–turbulence interactions, the Stokes number and settling parameter are in the range typically observed in atmospheric clouds. Observed fluctuations in droplet number concentration and diameter suggest a preference for inhomogeneous mixing. Finally, enhanced variance in liquid water content fluctuations is observed at high frequencies, and the scale break occurs at a value consistent with the independently estimated phase relaxation time from microphysical measurements.

Description: An adsorption theory of heterogeneous nucleation of water vapour on nanoparticles Atmospheric Chemistry and Physics Discussions, 15, 21883-21906, 2015 Author(s): A. Laaksonen and J. Malila Heterogeneous nucleation of water vapour on insoluble nuclei is a phenomenon that can induce atmospheric water and ice cloud formation. However, modelling of the phenomenon is hampered by the fact that the predictive capability of the classical heterogeneous nucleation theory is rather poor. A reliable theoretical description of the influence of different types of water-insoluble nuclei in triggering the water condensation or ice deposition would help to decrease uncertainty in large scale model simulations. In this paper we extend a recently formulated adsorption theory of heterogeneous nucleation to be applicable to highly curved surfaces, and test the theory against laboratory data for water vapour nucleation on silica, titanium dioxide and silver oxide nanoparticles. We show that unlike the classical heterogeneous nucleation theory, the new theory is able to quantitatively predict the experimental results.

Description: High-resolution daily gridded datasets of air temperature and wind speed for Europe Earth System Science Data Discussions, 8, 649-702, 2015 Author(s): S. Brinckmann, S. Krähenmann, and P. Bissolli New high-resolution datasets for near surface daily air temperature (minimum, maximum and mean) and daily mean wind speed for Europe (the CORDEX domain) are provided for the period 2001–2010 for the purpose of regional model validation in the framework of DecReg, a sub-project of the German MiKlip project, which aims to develop decadal climate predictions. The main input data sources are hourly SYNOP observations, partly supplemented by station data from the ECA&D dataset ( http://www.ecad.eu ). These data are quality tested to eliminate erroneous data and various kinds of inhomogeneities. Grids in a resolution of 0.044° (5 km) are derived by spatial interpolation of these station data into the CORDEX area. For temperature interpolation a modified version of a regression kriging method developed by Krähenmann et al. (2011) is used. At first, predictor fields of altitude, continentality and zonal mean temperature are chosen for a regression applied to monthly station data. The residuals of the monthly regression and the deviations of the daily data from the monthly averages are interpolated using simple kriging in a second and third step. For wind speed a new method based on the concept used for temperature was developed, involving predictor fields of exposure, roughness length, coastal distance and ERA Interim reanalysis wind speed at 850 hPa. Interpolation uncertainty is estimated by means of the kriging variance and regression uncertainties. Furthermore, to assess the quality of the final daily grid data, cross validation is performed. Explained variance ranges from 70 to 90 % for monthly temperature and from 50 to 60 % for monthly wind speed. The resulting RMSE for the final daily grid data amounts to 1–2 °C and 1–1.5 m s −1 (depending on season and parameter) for daily temperature parameters and daily mean wind speed, respectively. The datasets presented in this article are published at http://dx.doi.org/10.5676/DWD_CDC/DECREG0110v1 .

Description: Storm-triggered landslides in the Peruvian Andes and implications for topography, carbon cycles, and biodiversity Earth Surface Dynamics Discussions, 3, 631-688, 2015 Author(s): K. E. Clark, A. J. West, R. G. Hilton, G. P. Asner, C. A. Quesada, M. R. Silman, S. S. Saatchi, W. Farfan-Rios, R. E. Martin, A. B. Horwath, K. Halladay, M. New, and Y. Malhi In this study, we assess the geomorphic role of a rare, large-magnitude landslide event and consider the effect of this event on mountain forest ecosystems and the erosion of organic carbon in an Andean river catchment. Proximal triggers such as large rain storms are known to cause large numbers of landslides, but the relative effects of such low-frequency, high-magnitude events are not well known in the context of more regular, smaller events. We develop a 25 year duration, annual-resolution landslide inventory by mapping landslide occurrence in the Kosñipata Valley, Peru, from 1988 to 2012 using Landsat, Quickbird and Worldview satellite images. Catchment-wide landslide rates were high, at 0.076 % yr −1 by area, indicating landslides may completely turn over hillslopes every ~ 1320 years and strip 28 tC km −2 yr −1 of soil (73 %) and vegetation (27 %). A single rain storm in March 2010 accounted for 27 % of all landslide area observed during the 25 year study and removed 26 % of the organic carbon that was stripped from hillslopes by all landslides during the study. An approximately linear magnitude–frequency relationship for annual landslide areas suggests that large storms contribute an equivalent landslide failure area to the sum of smaller frequency landslides events occurring over the same period. However, the spatial distribution of landslides associated with the 2010 storm is distinct. On the basis of precipitation statistics and landscape morphology, we hypothesize that spatial focusing of storm-triggered landslide erosion at lower elevations in the Kosñipata catchment may be characteristic of longer-term patterns. These patterns may have implications for the source and composition of sediments and organic material supplied to river systems of the Amazon basin, and, through focusing of regular ecological disturbance, for the species composition of forested ecosystems in the region.

Description: A Late Pleistocene sea level stack Climate of the Past Discussions, 11, 3699-3728, 2015 Author(s): R. M. Spratt and L. E. Lisiecki Late Pleistocene sea level has been reconstructed from ocean sediment core data using a wide variety of proxies and models. However, the accuracy of individual reconstructions is limited by measurement error, local variations in salinity and temperature, and assumptions particular to each technique. Here we present a sea level stack (average) which increases the signal-to-noise ratio of individual reconstructions. Specifically, we perform principal component analysis (PCA) on seven records from 0–430 ka and five records from 0–798 ka. The first principal component, which we use as the stack, describes ~80 % of the variance in the data and is similar using either five or seven records. After scaling the stack based on Holocene and Last Glacial Maximum (LGM) sea level estimates, the stack agrees to within 5 m with isostatically adjusted coral sea level estimates for Marine Isotope Stages 5e and 11 (125 and 400 ka, respectively). When we compare the sea level stack with the δ 18 O of benthic foraminifera, we find that sea level change accounts for about ~40 % of the total orbital-band variance in benthic δ 18 O, compared to a 65 % contribution during the LGM-to-Holocene transition. Additionally, the second and third principal components of our analyses reflect differences between proxy records associated with spatial variations in the δ 18 O of seawater.

Description: Secondary organic aerosol formation from photochemical aging of light-duty gasoline vehicle exhausts in a smog chamber Atmospheric Chemistry and Physics, 15, 9049-9062, 2015 Author(s): T. Liu, X. Wang, W. Deng, Q. Hu, X. Ding, Y. Zhang, Q. He, Z. Zhang, S. Lü, X. Bi, J. Chen, and J. Yu In China, a rapid increase in passenger vehicles has led to the growing concern of vehicle exhaust as an important source of anthropogenic secondary organic aerosol (SOA) in megacities hard hit by haze. In this study, the SOA formation of emissions from two idling light-duty gasoline vehicles (LDGVs) (Euro 1 and Euro 4) operated in China was investigated in a 30 m 3 smog chamber. Five photo-oxidation experiments were carried out at 25 °C with relative humidity at around 50 %. After aging at an OH exposure of 5 × 10 6 molecules cm −3 h, the formed SOA was 12–259 times as high as primary organic aerosol (POA). The SOA production factors (PF) were 0.001–0.044 g kg −1 fuel, comparable with those from the previous studies at comparable OH exposure. This quite lower OH exposure than that in typical atmospheric conditions might however lead to the underestimation of the SOA formation potential from LDGVs. Effective SOA yields in this study were well fit by a one-product gas-particle partitioning model but quite lower than those of a previous study investigating SOA formation from three idling passenger vehicles (Euro 2–4). Traditional single-ring aromatic precursors and naphthalene could explain 51–90 % of the formed SOA. Unspeciated species such as branched and cyclic alkanes might be the possible precursors for the unexplained SOA. A high-resolution time-of-flight aerosol mass spectrometer was used to characterize the chemical composition of SOA. The relationship between f 43 (ratio of m/z 43, mostly C 2 H 3 O + , to the total signal in mass spectrum) and f 44 (mostly CO 2 + ) of the gasoline vehicle exhaust SOA is similar to the ambient semi-volatile oxygenated organic aerosol (SV-OOA). We plot the O : C and H : C molar ratios of SOA in a Van Krevelen diagram. The slopes of ΔH : C / ΔO : C ranged from −0.59 to −0.36, suggesting that the oxidation chemistry in these experiments was a combination of carboxylic acid and alcohol/peroxide formation.

Description: Nitrogen export from a boreal stream network following forest harvesting: seasonal nitrate removal and conservative export of organic forms Biogeosciences Discussions, 12, 12061-12089, 2015 Author(s): J. Schelker, R. Sponseller, E. Ring, L. Högbom, S. Löfgren, and H. Laudon Boreal streams are under pressure from large scale disturbance by forestry. Recent scenarios predict an increase in forest production in Scandinavia to meet market demands and to mitigate higher anthropogenic CO 2 emissions. Increased fertilization and shorter forest rotations are anticipated which will likely enhance the pressure on boreal streams in the near future. Among the major environmental impacts of forest harvesting is the increased mobilization of inorganic nitrogen (N), primarily as nitrate (NO 3 - ) into surface waters. But whereas NO 3 - inputs to first-order streams have been previously described, their downstream fate and impact is not well understood. We evaluated the downstream fate of N inputs in a boreal landscape that has been altered by forest harvests over a 10 year period to estimate the effects of multiple clear-cuts on aquatic N export in a boreal stream network. Small streams showed substantial leaching of NO 3 - in response to harvests with concentrations increasing by ~ 15 fold. NO 3 - concentrations at two sampling stations further downstream in the network were strongly seasonal and increased significantly in response to harvesting at the medium size, but not at the larger stream. Nitrate removal efficiency, E r , calculated as the percentage of "forestry derived" NO 3 - that was retained within the landscape using a mass balance model was highest during the snow melt season followed by the growing season, but declined continuously throughout the dormant season. In contrast, export of organic N from the landscape indicated little removal and was essentially conservative. Overall, net removal of NO 3 - between 2008 and 2011 accounted for ~ 70 % of the total NO 3 - mass exported from harvested patches distributed across the landscape. These results highlight the capacity and limitation of N-limited terrestrial and aquatic ecosystems to buffer inorganic N mobilization that arises from multiple clear-cuts within meso-scale boreal watersheds.

Description: Uncertainties in global aerosols and climate effects due to biofuel emissions Atmospheric Chemistry and Physics, 15, 8577-8596, 2015 Author(s): J. K. Kodros, C. E. Scott, S. C. Farina, Y. H. Lee, C. L'Orange, J. Volckens, and J. R. Pierce Aerosol emissions from biofuel combustion impact both health and climate; however, while reducing emissions through improvements to combustion technologies will improve health, the net effect on climate is largely unconstrained. In this study, we examine sensitivities in global aerosol concentration, direct radiative climate effect, and cloud-albedo aerosol indirect climate effect to uncertainties in biofuel emission factors, optical mixing state, and model nucleation and background secondary organic aerosol (SOA). We use the Goddard Earth Observing System global chemical-transport model (GEOS-Chem) with TwO Moment Aerosol Sectional (TOMAS) microphysics. The emission factors include amount, composition, size, and hygroscopicity, as well as optical mixing-state properties. We also evaluate emissions from domestic coal use, which is not biofuel but is also frequently emitted from homes. We estimate the direct radiative effect assuming different mixing states (homogeneous, core-shell, and external) with and without absorptive organic aerosol (brown carbon). We find the global-mean direct radiative effect of biofuel emissions ranges from −0.02 to +0.06 W m −2 across all simulation/mixing-state combinations with regional effects in source regions ranging from −0.2 to +0.8 W m −2 . The global-mean cloud-albedo aerosol indirect effect (AIE) ranges from +0.01 to −0.02 W m −2 with regional effects in source regions ranging from −1.0 to −0.05 W m −2 . The direct radiative effect is strongly dependent on uncertainties in emissions mass, composition, emissions aerosol size distributions, and assumed optical mixing state, while the indirect effect is dependent on the emissions mass, emissions aerosol size distribution, and the choice of model nucleation and secondary organic aerosol schemes. The sign and magnitude of these effects have a strong regional dependence. We conclude that the climate effects of biofuel aerosols are largely unconstrained, and the overall sign of the aerosol effects is unclear due to uncertainties in model inputs. This uncertainty limits our ability to introduce mitigation strategies aimed at reducing biofuel black carbon emissions in order to counter warming effects from greenhouse gases. To better understand the climate impact of particle emissions from biofuel combustion, we recommend field/laboratory measurements to narrow constraints on (1) emissions mass, (2) emission size distribution, (3) mixing state, and (4) ratio of black carbon to organic aerosol.

Description: The Ice Selective Inlet: a novel technique for exclusive extraction of pristine ice crystals in mixed-phase clouds Atmospheric Measurement Techniques, 8, 3087-3106, 2015 Author(s): P. Kupiszewski, E. Weingartner, P. Vochezer, M. Schnaiter, A. Bigi, M. Gysel, B. Rosati, E. Toprak, S. Mertes, and U. Baltensperger Climate predictions are affected by high uncertainties partially due to an insufficient knowledge of aerosol–cloud interactions. One of the poorly understood processes is formation of mixed-phase clouds (MPCs) via heterogeneous ice nucleation. Field measurements of the atmospheric ice phase in MPCs are challenging due to the presence of much more numerous liquid droplets. The Ice Selective Inlet (ISI), presented in this paper, is a novel inlet designed to selectively sample pristine ice crystals in mixed-phase clouds and extract the ice residual particles contained within the crystals for physical and chemical characterization. Using a modular setup composed of a cyclone impactor, droplet evaporation unit and pumped counterflow virtual impactor (PCVI), the ISI segregates particles based on their inertia and phase, exclusively extracting small ice particles between 5 and 20 μm in diameter. The setup also includes optical particle spectrometers for analysis of the number size distribution and shape of the sampled hydrometeors. The novelty of the ISI is a droplet evaporation unit, which separates liquid droplets and ice crystals in the airborne state, thus avoiding physical impaction of the hydrometeors and limiting potential artefacts. The design and validation of the droplet evaporation unit is based on modelling studies of droplet evaporation rates and computational fluid dynamics simulations of gas and particle flows through the unit. Prior to deployment in the field, an inter-comparison of the optical particle size spectrometers and a characterization of the transmission efficiency of the PCVI was conducted in the laboratory. The ISI was subsequently deployed during the Cloud and Aerosol Characterization Experiment (CLACE) 2013 and 2014 – two extensive international field campaigns encompassing comprehensive measurements of cloud microphysics, as well as bulk aerosol, ice residual and ice nuclei properties. The campaigns provided an important opportunity for a proof of concept of the inlet design. In this work we present the setup of the ISI, including the modelling and laboratory characterization of its components, as well as field measurements demonstrating the ISI performance and validating the working principle of the inlet. Finally, measurements of biological aerosol during a Saharan dust event (SDE) are presented, showing a first indication of enrichment of bio-material in sub-2 μm ice residuals.

Description: Impacts on wave-driven harbour agitation due to climate change in Catalan ports Natural Hazards and Earth System Science, 15, 1695-1709, 2015 Author(s): J. P. Sierra, M. Casas-Prat, M. Virgili, C. Mösso, and A. Sánchez-Arcilla The objective of the present work is to analyse how changes in wave patterns due to the effect of climate change can affect harbour agitation (oscillations within the port due to wind waves). The study focuses on 13 harbours located on the Catalan coast (NW Mediterranean) using a methodology with general applicability. To obtain the patterns of agitation, a Boussinesq-type model is used, which is forced at the boundaries by present/future offshore wave conditions extracted from recently developed high-resolution wave projections in the NW Mediterranean. These wave projections were obtained with the SWAN model forced by present/future surface wind fields projected, respectively, by five different combinations of global and regional circulation models (GCMs and RCMs) for the A1B scenario. The results show a general slight reduction in the annual average agitation for most of the ports, except for the northernmost and southernmost areas of the region, where a slight increase is obtained. A seasonal analysis reveals that the tendency to decrease is accentuated in winter. However, the inter-model variability is large for both the winter and the annual analysis. Conversely, a general increase with a larger agreement among models is found during summer, which is the period with greater activity in most of the studied ports (marinas). A qualitative assessment of the factors of variability seems to indicate that the choice of GCM tends to affect the spatial pattern, whereas the choice of RCM induces a more homogeneous bias over the regional domain.

Description: Utilisation of CryoSat-2 SAR altimeter in operational ice charting The Cryosphere Discussions, 9, 4117-4145, 2015 Author(s): E. Rinne and M. Similä We present methods to utilise Cryosat-2 (CS-2) Synthetic Aperture (SAR) mode data in operational ice charting. We compare CS-2 data qualitatively to Synthetic Aperture Radar (SAR) mosaics over Barents and Kara seas. Furthermore, we compare the CS-2 to archived operational ice charts. We present distributions of four CS-2 waveform parameters for different ice types as presented in the ice charts. We go on to present an automatic classification method for CS-2 data which, after training with operational ice charts, is capable of determining open water from ice with a hit rate of 〉 90 %. The training data is dynamically updated every five days using the most recent 15 days CS-2 data and operative ice charts. This helps the adaption of the classifier to the evolving ice/snow conditions throughout winter. The classifier is also capable of detecting three different ice classes (thin and thick first year ice as well as old ice) with success rates good enough for the output to be usable to support operational ice charting. Finally, we present a near real time CS-2 product just plotting the waveform characteristics and conclude that even such a simple product is usable for some of the needs of ice charting.

Description: Inverse modeling of black carbon emissions over China using ensemble data assimilation Atmospheric Chemistry and Physics Discussions, 15, 20851-20879, 2015 Author(s): P. Wang, H. Wang, Y. Q. Wang, X. Y. Zhang, S. L. Gong, M. Xue, C. H. Zhou, H. L. Liu, X. Q. An, T. Niu, and Y. L. Chen Emissions inventories of black carbon (BC), which are traditionally constructed using a "bottom-up" approach based on activity data and emissions factors, are considered to contain a large level of uncertainty. In this paper, an ensemble optimal interpolation (EnOI) data assimilation technique is used to investigate the possibility of optimally recovering the spatially resolved emissions bias of BC. An inverse modeling system for emissions is established for an atmospheric chemistry aerosol model and two key problems related to ensemble data assimilation in the top-down emissions estimation are discussed: (1) how to obtain reasonable ensembles of prior emissions; and (2) establishing a scheme to localize the background-error matrix. An experiment involving a one month simulation cycle with EnOI inversion of BC emissions is performed for January 2008. The bias of the BC emissions intensity in China at each grid point is corrected by this inverse system. The inversed emission over China in January is 240.1 Gg, and annual emission is about 2750 Gg, which is over 1.8 times of bottom-up emission inventory. The results show that, even though only monthly mean BC measurements are employed to inverse the emissions, the accuracy of the daily model simulation improves. Using top-down emissions, the average root-mean-square error of simulated daily BC is decreased by nearly 30 %. These results are valuable and promising for a better understanding of aerosol emissions and distributions, as well as aerosol forecasting.

Description: Metrology of ground-based satellite validation: co-location mismatch and smoothing issues of total ozone comparisons Atmospheric Measurement Techniques Discussions, 8, 8023-8082, 2015 Author(s): T. Verhoelst, J. Granville, F. Hendrick, U. Köhler, C. Lerot, J.-P. Pommereau, A. Redondas, M. Van Roozendael, and J.-C. Lambert Comparisons with ground-based correlative measurements constitute a key component in the validation of satellite data on atmospheric composition. The error budget of these comparisons contains not only the measurement uncertainties but also several terms related to differences in sampling and smoothing of the inhomogeneous and variable atmospheric field. A versatile system for Observing System Simulation Experiments (OSSEs), named OSSSMOSE, is used here to quantify these terms. Based on the application of pragmatic observation operators onto high-resolution atmospheric fields, it allows a simulation of each individual measurement, and consequently also of the differences to be expected from spatial and temporal field variations between both measurements making up a comparison pair. As a topical case study, the system is used to evaluate the error budget of total ozone column (TOC) comparisons between on the one hand GOME-type direct fitting (GODFITv3) satellite retrievals from GOME/ERS2, SCIAMACHY/Envisat, and GOME-2/MetOp-A, and on the other hand direct-sun and zenith-sky reference measurements such as from Dobsons, Brewers, and zenith scattered light (ZSL-)DOAS instruments respectively. In particular, the focus is placed on the GODFITv3 reprocessed GOME-2A data record vs. the ground-based instruments contributing to the Network for the Detection of Atmospheric Composition Change (NDACC). The simulations are found to reproduce the actual measurements almost to within the measurement uncertainties, confirming that the OSSE approach and its technical implementation are appropriate. This work reveals that many features of the comparison spread and median difference can be understood as due to metrological differences, even when using strict co-location criteria. In particular, sampling difference errors exceed measurement uncertainties regularly at most mid- and high-latitude stations, with values up to 10 % and more in extreme cases. Smoothing difference errors only play a role in the comparisons with ZSL-DOAS instruments at high latitudes, especially in the presence of a polar vortex. At tropical latitudes, where TOC variability is lower, both types of errors remain below about 1 % and consequently do not contribute significantly to the comparison error budget. The detailed analysis of the comparison results, including now the metrological errors, suggests that the published random measurement uncertainties for GODFITv3 reprocessed satellite data are potentially overestimated, and adjustments are proposed here. This successful application of the OSSSMOSE sytem to close for the first time the error budget of TOC comparisons, bodes well for potential future applications, which are briefly touched upon.

Description: Relating seasonal dynamics of enhanced vegetation index to the recycling of water in two endorheic river basins in north-west China Hydrology and Earth System Sciences, 19, 3387-3403, 2015 Author(s): M. A. Matin and C. P.-A. Bourque This study associates the dynamics of enhanced vegetation index in lowland desert oases to the recycling of water in two endorheic (hydrologically closed) river basins in Gansu Province, north-west China, along a gradient of elevation zones and land cover types. Each river basin was subdivided into four elevation zones representative of (i) oasis plains and foothills, and (ii) low-, (iii) mid-, and (iv) high-mountain elevations. Comparison of monthly vegetation phenology with precipitation and snowmelt dynamics within the same basins over a 10-year period (2000–2009) suggested that the onset of the precipitation season (cumulative % precipitation 〉 7–8 %) in the mountains, typically in late April to early May, was triggered by the greening of vegetation and increased production of water vapour at the base of the mountains. Seasonal evolution of in-mountain precipitation correlated fairly well with the temporal variation in oasis-vegetation coverage and phenology characterised by monthly enhanced vegetation index, yielding coefficients of determination of 0.65 and 0.85 for the two basins. Convergent cross-mapping of related time series indicated bi-directional causality (feedback) between the two variables. Comparisons between same-zone monthly precipitation amounts and enhanced vegetation index provided weaker correlations. Start of the growing season in the oases was shown to coincide with favourable spring warming and discharge of meltwater from low- to mid-elevations of the Qilian Mountains (zones 1 and 2) in mid-to-late March. In terms of plant requirement for water, mid-seasonal development of oasis vegetation was seen to be controlled to a greater extent by the production of rain in the mountains. Comparison of water volumes associated with in-basin production of rainfall and snowmelt with that associated with evaporation seemed to suggest that about 90 % of the available liquid water (i.e. mostly in the form of direct rainfall and snowmelt in the mountains) was recycled locally.

Description: Investigating suspended sediment dynamics in contrasting agricultural catchments using ex situ turbidity-based suspended sediment monitoring Hydrology and Earth System Sciences, 19, 3349-3363, 2015 Author(s): S. C. Sherriff, J. S. Rowan, A. R. Melland, P. Jordan, O. Fenton, and D. Ó hUallacháin Soil erosion and suspended sediment (SS) pose risks to chemical and ecological water quality. Agricultural activities may accelerate erosional fluxes from bare, poached or compacted soils, and enhance connectivity through modified channels and artificial drainage networks. Storm-event fluxes dominate SS transport in agricultural catchments; therefore, high temporal-resolution monitoring approaches are required, but can be expensive and technically challenging. Here, the performance of in situ turbidity sensors, conventionally installed submerged at the river bankside, is compared with installations where river water is delivered to sensors ex situ, i.e. within instrument kiosks on the riverbank, at two experimental catchments (Grassland B and Arable B). The in situ and ex situ installations gave comparable results when calibrated against storm-period, depth-integrated SS data, with total loads at Grassland B estimated at 12 800 and 15 400 t, and 22 600 and 24 900 t at Arable B, respectively. The absence of spurious turbidity readings relating to bankside debris around the in situ sensor and its greater security make the ex situ sensor more robust. The ex situ approach was then used to characterise SS dynamics and fluxes in five intensively managed agricultural catchments in Ireland which feature a range of landscape characteristics and land use pressures. Average annual suspended sediment concentration (SSC) was below the Freshwater Fish Directive (78/659/EEC) guideline of 25 mg L −1 , and the continuous hourly record demonstrated that exceedance occurred less than 12 % of the observation year. Soil drainage class and proportion of arable land were key controls determining flux rates, but all catchments reported a high degree of inter-annual variability associated with variable precipitation patterns compared to the long-term average. Poorly drained soils had greater sensitivity to runoff and soil erosion, particularly in catchments with periods of bare soils. Well drained soils were less sensitive to erosion even on arable land; however, under extreme rainfall conditions, all bare soils remain a high sediment loss risk. Analysis of storm-period and seasonal dynamics (over the long term) using high-resolution monitoring would be beneficial to further explore the impact of landscape, climate and land use characteristics on SS export.

Description: Landscape heterogeneity drives contrasting concentration–discharge relationships in shale headwater catchments Hydrology and Earth System Sciences, 19, 3333-3347, 2015 Author(s): E. M. Herndon, A. L. Dere, P. L. Sullivan, D. Norris, B. Reynolds, and S. L. Brantley Solute concentrations in stream water vary with discharge in patterns that record complex feedbacks between hydrologic and biogeochemical processes. In a comparison of three shale-underlain headwater catchments located in Pennsylvania, USA (the forested Shale Hills Critical Zone Observatory), and Wales, UK (the peatland-dominated Upper Hafren and forest-dominated Upper Hore catchments in the Plynlimon forest), dissimilar concentration–discharge ( C – Q ) behaviors are best explained by contrasting landscape distributions of soil solution chemistry – especially dissolved organic carbon (DOC) – that have been established by patterns of vegetation and soil organic matter (SOM). Specifically, elements that are concentrated in organic-rich soils due to biotic cycling (Mn, Ca, K) or that form strong complexes with DOC (Fe, Al) are spatially heterogeneous in pore waters because organic matter is heterogeneously distributed across the catchments. These solutes exhibit non-chemostatic behavior in the streams, and solute concentrations either decrease (Shale Hills) or increase (Plynlimon) with increasing discharge. In contrast, solutes that are concentrated in soil minerals and form only weak complexes with DOC (Na, Mg, Si) are spatially homogeneous in pore waters across each catchment. These solutes are chemostatic in that their stream concentrations vary little with stream discharge, likely because these solutes are released quickly from exchange sites in the soils during rainfall events. Furthermore, concentration–discharge relationships of non-chemostatic solutes changed following tree harvest in the Upper Hore catchment in Plynlimon, while no changes were observed for chemostatic solutes, underscoring the role of vegetation in regulating the concentrations of certain elements in the stream. These results indicate that differences in the hydrologic connectivity of organic-rich soils to the stream drive differences in concentration behavior between catchments. As such, in catchments where SOM is dominantly in lowlands (e.g., Shale Hills), we infer that non-chemostatic elements associated with organic matter are released to the stream early during rainfall events, whereas in catchments where SOM is dominantly in uplands (e.g., Plynlimon), these non-chemostatic elements are released later during rainfall events. The distribution of SOM across the landscape is thus a key component for predictive models of solute transport in headwater catchments.

Description: Future changes in flash flood frequency and intensity of the Tha Di River (Thailand) based on rainfall–runoff modeling and advanced delta change scaling Hydrology and Earth System Sciences Discussions, 12, 7327-7352, 2015 Author(s): S. Hilgert, A. Wagner, and S. Fuchs As a consequence of climate change, extreme and flood-causing precipitation events are expected to increase in magnitude and frequency, especially in today's high-precipitation areas. During the north-east monsoon seasons, Nakhon Si Thammarat in southern Thailand is flash-flooded every 2.22 years on average. This study investigates frequency and intensity of harmful discharges of the Tha Di River regarding the IPCC emission scenarios A2 and B2. The regional climate model (RCM) PRECIS was transformed using the advanced delta change (ADC) method. The hydrologic response model HBV-Light was calibrated to the catchment and supplied with ADC-scaled daily precipitation and temperature data for 2010–2089. Under the A2 (B2) scenario, the flood threshold exceedance frequency on average increases by 133 % (decreases by 10 %), average flood intensity increases by 3 % (decreases by 2 %) and the annual top five discharge peaks intensities increase by 46 % (decrease by 5 %). Yearly precipitation sums increase by 30 % (10 %) towards the end of the century. The A2 scenario predicts a precipitation increase during the rainy season, which intensifies flood events; while increases projected exclusively for the dry season are not expected to cause floods. Retention volume demand of past events was calculated to be up to 12 × 10 6 m 3 . Flood risks are staying at high levels under the B2 scenario or increase dramatically under the A2 scenario. Results show that the RCM scaling process is inflicted with systematic biases but is crucial to investigate small, mountainous catchments. Improvement of scaling techniques should therefore accompany the development towards high-resolution climate models.

Description: A pan-African medium-range ensemble flood forecast system Hydrology and Earth System Sciences, 19, 3365-3385, 2015 Author(s): V. Thiemig, B. Bisselink, F. Pappenberger, and J. Thielen The African Flood Forecasting System (AFFS) is a probabilistic flood forecast system for medium- to large-scale African river basins, with lead times of up to 15 days. The key components are the hydrological model LISFLOOD, the African GIS database, the meteorological ensemble predictions by the ECMWF (European Centre for Medium-Ranged Weather Forecasts) and critical hydrological thresholds. In this paper, the predictive capability is investigated in a hindcast mode, by reproducing hydrological predictions for the year 2003 when important floods were observed. Results were verified by ground measurements of 36 sub-catchments as well as by reports of various flood archives. Results showed that AFFS detected around 70 % of the reported flood events correctly. In particular, the system showed good performance in predicting riverine flood events of long duration (〉 1 week) and large affected areas (〉 10 000 km 2 ) well in advance, whereas AFFS showed limitations for small-scale and short duration flood events. The case study for the flood event in March 2003 in the Sabi Basin (Zimbabwe) illustrated the good performance of AFFS in forecasting timing and severity of the floods, gave an example of the clear and concise output products, and showed that the system is capable of producing flood warnings even in ungauged river basins. Hence, from a technical perspective, AFFS shows a large potential as an operational pan-African flood forecasting system, although issues related to the practical implication will still need to be investigated.

Description: Gas adsorption and desorption effects on cylinders and their importance for long-term gas records Atmospheric Measurement Techniques Discussions, 8, 8083-8112, 2015 Author(s): M. C. Leuenberger, M. F. Schibig, and P. Nyfeler It is well known that gases adsorb on many surfaces, in particular metal surfaces. There are two main forms responsible for these effects (i) physisorption and (ii) chemisorption. Physisorption is associated with lower binding energies in the order of 1–10 kJ mol −1 compared to chemisorption ranging from 100 to 1000 kJ mol −1 . Furthermore, chemisorption forms only monolayers, contrasting physisorption that can form multilayer adsorption. The reverse process is called desorption and follows similar mathematical laws, however, it can be influenced by hysteresis effects. In the present experiment we investigated the adsorption/desorption phenomena on three steel and three aluminium cylinders containing compressed air in our laboratory and under controlled conditions in a climate chamber, respectively. We proved the pressure effect on physisorption for CO 2 , CH 4 and H 2 O by decanting one steel and two aluminium cylinders completely. The CO 2 results for both cylinders are in excellent agreement with the pressure dependence of a monolayer adsorption model. However, adsorption on aluminium ( 〈 0.05 and 0 ppm for CO 2 and H 2 O) was about 10 times less than on steel ( 〈 0.41 ppm and about 〈 2.5 ppm, respectively). The CO 2 amounts adsorbed (5.8 × 10 19 CO 2 molecules) corresponds to about the five-fold monolayer adsorption indicating that the effective surface exposed for adsorption is significantly larger than the geometric surface area. Adsorption/desorption effects were minimal for CH 4 and for CO. However, the latter dependence requires further attention since it was only studied on one aluminium cylinder with a very low mole fraction. In the climate chamber the cylinders were exposed to temperatures between −10 and +50 °C to determine the corresponding temperature coefficients of adsorption. Again, we found distinctly different values for CO 2 ranging from 0.0014 to 0.0184 ppm °C −1 for steel cylinders and −0.0002 to −0.0003 ppm °C −1 for aluminium cylinders. The reversed temperature dependence for aluminium cylinders point to significantly lower desorption energies than for steel cylinders and might at least partly be due to temperature and gas consumption induced pressure changes. Temperature coefficients for CH 4 , CO and H 2 O adsorption were, within their error bands, insignificant. These results do indicate the need for careful selection and usage of gas cylinders for high precision calibration purposes such as requested in trace gas applications.

Description: A two-layer canopy model with thermal inertia for an improved snowpack energy balance below needleleaf forest (model SNOWPACK, version 3.2.1, revision 741) Geoscientific Model Development, 8, 2379-2398, 2015 Author(s): I. Gouttevin, M. Lehning, T. Jonas, D. Gustafsson, and M. Mölder A new, two-layer canopy module with thermal inertia as part of the detailed snow model SNOWPACK (version 3.2.1) is presented and evaluated. As a by-product of these new developments, an exhaustive description of the canopy module of the SNOWPACK model is provided, thereby filling a gap in the existing literature. In its current form, the two-layer canopy module is suited for evergreen needleleaf forest, with or without snow cover. It is designed to reproduce the difference in thermal response between leafy and woody canopy elements, and their impact on the underlying snowpack or ground surface energy balance. Given the number of processes resolved, the SNOWPACK model with its enhanced canopy module constitutes a sophisticated physics-based modeling chain of the continuum going from atmosphere to soil through the canopy and snow. Comparisons of modeled sub-canopy thermal radiation to stand-scale observations at an Alpine site (Alptal, Switzerland) demonstrate improvements induced by the new canopy module. Both thermal heat mass and the two-layer canopy formulation contribute to reduce the daily amplitude of the modeled canopy temperature signal, in agreement with observations. Particularly striking is the attenuation of the nighttime drop in canopy temperature, which was a key model bias. We specifically show that a single-layered canopy model is unable to produce this limited temperature drop correctly. The impact of the new parameterizations on the modeled dynamics of the sub-canopy snowpack is analyzed. The new canopy module yields consistent results but the frequent occurrence of mixed-precipitation events at Alptal prevents a conclusive assessment of model performance against snow data. The new model is also successfully tested without specific tuning against measured tree temperature and biomass heat-storage fluxes at the boreal site of Norunda (Sweden). This provides an independent assessment of its physical consistency and stresses the robustness and transferability of the chosen parameterizations. The SNOWPACK code including the new canopy module, is available under Gnu General Public License (GPL) license and upon creation of an account at https://models.slf.ch/ .

Description: A comprehensive filtering scheme for high-resolution estimation of the water balance components from high-precision lysimeters Hydrology and Earth System Sciences, 19, 3405-3418, 2015 Author(s): M. Hannes, U. Wollschläger, F. Schrader, W. Durner, S. Gebler, T. Pütz, J. Fank, G. von Unold, and H.-J. Vogel Large weighing lysimeters are currently the most precise method to directly measure all components of the terrestrial water balance in parallel via the built-in weighing system. As lysimeters are exposed to several external forces such as management practices or wind influencing the weighing data, the calculated fluxes of precipitation and evapotranspiration can be altered considerably without having applied appropriate corrections to the raw data. Therefore, adequate filtering schemes for obtaining most accurate estimates of the water balance components are required. In this study, we use data from the TERENO (TERrestrial ENvironmental Observatories) SoilCan research site in Bad Lauchstädt to develop a comprehensive filtering procedure for high-precision lysimeter data, which is designed to deal with various kinds of possible errors starting from the elimination of large disturbances in the raw data resulting e.g., from management practices all the way to the reduction of noise caused e.g., by moderate wind. Furthermore, we analyze the influence of averaging times and thresholds required by some of the filtering steps on the calculated water balance and investigate the ability of two adaptive filtering methods (the adaptive window and adaptive threshold filter (AWAT filter; Peters et al., 2014), and a new synchro filter applicable to the data from a set of several lysimeters) to further reduce the filtering error. Finally, we take advantage of the data sets of all 18 lysimeters running in parallel at the Bad Lauchstädt site to evaluate the performance and accuracy of the proposed filtering scheme. For the tested time interval of 2 months, we show that the estimation of the water balance with high temporal resolution and good accuracy is possible. The filtering code can be downloaded from the journal website as Supplement to this publication.

Description: Assimilation of Antarctic velocity observations provides evidence for uncharted pinning points The Cryosphere, 9, 1427-1443, 2015 Author(s): J. J. Fürst, G. Durand, F. Gillet-Chaulet, N. Merino, L. Tavard, J. Mouginot, N. Gourmelen, and O. Gagliardini In ice flow modelling, the use of control methods to assimilate the dynamic and geometric state of an ice body has become common practice. These methods have primarily focussed on inverting for one of the two least known properties in glaciology, namely the basal friction coefficient or the ice viscosity parameter. Here, we present an approach to infer both properties simultaneously for the whole of the Antarctic ice sheet. After the assimilation, the root-mean-square deviation between modelled and observed surface velocities attains 8.7 m a −1 for the entire domain, with a slightly higher value of 14.0 m a −1 for the ice shelves. An exception in terms of the velocity mismatch is the Thwaites Glacier Ice Shelf, where the RMS value is almost 70 m a −1 . The reason is that the underlying Bedmap2 geometry ignores the presence of an ice rise, which exerts major control on the dynamics of the eastern part of the ice shelf. On these grounds, we suggest an approach to account for pinning points not included in Bedmap2 by locally allowing an optimisation of basal friction during the inversion. In this way, the velocity mismatch on the ice shelf of Thwaites Glacier is more than halved. A characteristic velocity mismatch pattern emerges for unaccounted pinning points close to the marine shelf front. This pattern is exploited to manually identify seven uncharted features around Antarctica that exert significant resistance to the shelf flow. Potential pinning points are detected on Fimbul, West, Shackleton, Nickerson and Venable ice shelves. As pinning points can provide substantial resistance to shelf flow, with considerable consequences if they became ungrounded in the future, the model community is in need of detailed bathymetry there. Our data assimilation points to some of these dynamically important features not present in Bedmap2 and implicitly quantifies their relevance.

Description: Long-term coastal-polynya dynamics in the Southern Weddell Sea from MODIS thermal-infrared imagery The Cryosphere Discussions, 9, 3959-3993, 2015 Author(s): S. Paul, S. Willmes, and G. Heinemann Based upon high-resolution thermal-infrared Moderate-Resolution Imaging Spectroradiometer (MODIS) satellite imagery in combination with ERA-Interim atmospheric reanalysis data, we derived long-term polynya parameters such as polynya area, thin-ice thickness distribution and ice-production rates from daily cloud-cover corrected thin-ice thickness composites. Our study is based on a thirteen year investigation period (2002–2014) for the austral winter (1 April to 30 September) in the Antarctic Southern Weddell Sea. The focus lies on coastal polynyas which are important hot spots for new-ice formation, bottom-water formation and heat/moisture release into the atmosphere. MODIS has the capability to resolve even very narrow coastal polynyas. Its major disadvantage is the sensor limitation due to cloud cover. We make use of a newly developed and adapted spatial feature reconstruction scheme to account for cloud-covered areas. We find the sea-ice areas in front of Ronne and Brunt Ice Shelf to be the most active with an annual average polynya area of 3018 ± 1298 and 3516 ± 1420 km 2 as well as an accumulated volume ice production of 31 ± 13 and 31 ± 12 km 3 , respectively. For the remaining four regions, estimates amount to 421 ± 294 km 2 and 4 ± 3 km 3 (Antarctic Peninsula), 1148 ± 432 km 2 and 12 ± 5 km 3 (Iceberg A23A), 901 ± 703 km 2 and 10 ± 8 km 3 (Filchner Ice Shelf) as well as 499 ± 277 km 2 and 5 ± 2 km 3 (Coats Land). Our findings are discussed in comparison to recent studies based on coupled sea-ice/ocean models and passive-microwave satellite imagery, each investigating different parts of the Southern Weddell Sea.

Description: Climatic controls and climate proxy potential of Lewis Glacier, Mt Kenya The Cryosphere Discussions, 9, 3887-3924, 2015 Author(s): R. Prinz, L. I. Nicholson, T. Mölg, W. Gurgiser, and G. Kaser The Lewis Glacier on Mt Kenya is one of the best studied tropical glaciers and has experienced considerable retreat since a maximum extent in the late 19th century (L19). From distributed mass and energy balance modelling, this study evaluates the current sensitivity of the surface mass and energy balance to climatic drivers, explores climate conditions under which the L19 maximum extent might have sustained, and discusses the potential for using the glacier retreat to quantify climate change. Multiyear meteorological measurements at 4828 m provide data for input, optimization and evaluation of a spatially distributed glacier mass balance model to quantify the exchanges of energy and mass at the glacier–atmosphere interface. Currently the glacier loses mass due to the imbalance between insufficient accumulation and enhanced melt, because radiative energy gains cannot be compensated by turbulent energy sinks. Exchanging model input data with synthetic climate scenarios, which were sampled from the meteorological measurements and account for coupled climatic variable perturbations, reveal that the current mass balance is most sensitive to changes in atmospheric moisture (via its impact on solid precipitation, cloudiness and surface albedo). Positive mass balances result from scenarios with an increase of annual (seasonal) accumulation of 30 % (100 %), compared to values observed today, without significant changes in air temperature required. Scenarios with lower air temperatures are drier and associated with lower accumulation and increased net radiation due to reduced cloudiness and albedo. If the scenarios currently producing positive mass balances are applied to the L19 extent, negative mass balances are the result, meaning that the conditions required to sustain the glacier in its L19 extent are not reflected in today's observations. Alternatively, a balanced mass budget for the L19 extent can be explained by changing model parameters that imply a distinctly different coupling between the glacier's local surface-air layer and its surrounding boundary-layer. This result underlines the difficulty of deriving paleoclimates for larger glacier extents on the basis of modern measurements of small glaciers.

Description: A Web-based spatial decision supporting system for land management and soil conservation Solid Earth, 6, 903-928, 2015 Author(s): F. Terribile, A. Agrillo, A. Bonfante, G. Buscemi, M. Colandrea, A. D'Antonio, R. De Mascellis, C. De Michele, G. Langella, P. Manna, L. Marotta, F. A. Mileti, L. Minieri, N. Orefice, S. Valentini, S. Vingiani, and A. Basile Today it is evident that there are many contrasting demands on our landscape (e.g. food security, more sustainable agriculture, higher income in rural areas, etc.) as well as many land degradation problems. It has been proved that providing operational answers to these demands and problems is extremely difficult. Here we aim to demonstrate that a spatial decision support system based on geospatial cyberinfrastructure (GCI) can address all of the above, so producing a smart system for supporting decision making for agriculture, forestry, and urban planning with respect to the landscape. In this paper, we discuss methods and results of a special kind of GCI architecture, one that is highly focused on land management and soil conservation. The system allows us to obtain dynamic, multidisciplinary, multiscale, and multifunctional answers to agriculture, forestry, and urban planning issues through the Web. The system has been applied to and tested in an area of about 20 000 ha in the south of Italy, within the framework of a European LIFE+ project (SOILCONSWEB). The paper reports – as a case study – results from two different applications dealing with agriculture (olive growth tool) and environmental protection (soil capability to protect groundwater). Developed with the help of end users, the system is starting to be adopted by local communities. The system indirectly explores a change of paradigm for soil and landscape scientists. Indeed, the potential benefit is shown of overcoming current disciplinary fragmentation over landscape issues by offering – through a smart Web-based system – truly integrated geospatial knowledge that may be directly and freely used by any end user ( www.landconsultingweb.eu ). This may help bridge the last very important divide between scientists working on the landscape and end users.

Description: Towards predictive data-driven simulations of wildfire spread – Part II: Ensemble Kalman Filter for the state estimation of a front-tracking simulator of wildfire spread Natural Hazards and Earth System Science, 15, 1721-1739, 2015 Author(s): M. C. Rochoux, C. Emery, S. Ricci, B. Cuenot, and A. Trouvé This paper is the second part in a series of two articles, which aims at presenting a data-driven modeling strategy for forecasting wildfire spread scenarios based on the assimilation of the observed fire front location and on the sequential correction of model parameters or model state. This model relies on an estimation of the local rate of fire spread (ROS) as a function of environmental conditions based on Rothermel's semi-empirical formulation, in order to propagate the fire front with an Eulerian front-tracking simulator. In Part I, a data assimilation (DA) system based on an ensemble Kalman filter (EnKF) was implemented to provide a spatially uniform correction of biomass fuel and wind parameters and thereby, produce an improved forecast of the wildfire behavior (addressing uncertainties in the input parameters of the ROS model only). In Part II, the objective of the EnKF algorithm is to sequentially update the two-dimensional coordinates of the markers along the discretized fire front, in order to provide a spatially distributed correction of the fire front location and thereby, a more reliable initial condition for further model time-integration (addressing all sources of uncertainties in the ROS model). The resulting prototype data-driven wildfire spread simulator is first evaluated in a series of verification tests using synthetically generated observations; tests include representative cases with spatially varying biomass properties and temporally varying wind conditions. In order to properly account for uncertainties during the EnKF update step and to accurately represent error correlations along the fireline, it is shown that members of the EnKF ensemble must be generated through variations in estimates of the fire's initial location as well as through variations in the parameters of the ROS model. The performance of the prototype simulator based on state estimation (SE) or parameter estimation (PE) is then evaluated by comparison with data taken from a reduced-scale controlled grassland fire experiment. Results indicate that data-driven simulations are capable of correcting inaccurate predictions of the fire front location and of subsequently providing an optimized forecast of the wildfire behavior at future lead times. The complementary benefits of both PE and SE approaches, in terms of analysis and forecast performance, are also emphasized. In particular, it is found that the size of the assimilation window must be specified adequately with the persistence of the model initial condition and/or with the temporal and spatial variability of the environmental conditions in order to track sudden changes in wildfire behavior. The present prototype data-driven forecast system is still at an early stage of development. In this regard, this preliminary investigation provides valuable information on how to combine observations with a fire spread model in an efficient way, as well as guidelines to design the future system evolution in order to meet the operational requirements of wildfire spread monitoring.

Description: High resolution satellite turbidity and sea surface temperature observations of river plume interactions during a significant flood event Ocean Science Discussions, 12, 1669-1692, 2015 Author(s): V. E. Brando, F. Braga, L. Zaggia, C. Giardino, M. Bresciani, D. Bellafiore, C. Ferrarin, F. Maicu, A. Benetazzo, D. Bonaldo, F. M. Falcieri, A. Coluccelli, A. Russo, and S. Carniel Sea surface temperature (SST) and turbidity ( T ) derived from Landsat-8 (L8) imagery were used to characterize river plumes in the Northern Adriatic Sea (NAS) during a significant flood event in November 2014. Sea surface salinity (SSS) from an operational coupled ocean-wave model supported the interpretation of the plumes interaction with the receiving waters and among them. There was a good agreement of the SSS, T , and SST fields at the sub-mesoscale and mesoscale delineation of the major river plumes. L8 30 m resolution enabled also the description of smaller plume structures. Sharp fronts in T and SST delimited each single river plume. The isotherms and turbidity isolines coupling varied among the plumes due to differences in particle loads and surface temperatures in the discharged waters. The different plumes reflectance spectra were related to the lithological fingerprint of the sediments in the river catchments.

Description: Sea level trend and variability around Peninsular Malaysia Ocean Science, 11, 617-628, 2015 Author(s): Q. H. Luu, P. Tkalich, and T. W. Tay Sea level rise due to climate change is non-uniform globally, necessitating regional estimates. Peninsular Malaysia is located in the middle of Southeast Asia, bounded from the west by the Malacca Strait, from the east by the South China Sea (SCS), and from the south by the Singapore Strait. The sea level along the peninsula may be influenced by various regional phenomena native to the adjacent parts of the Indian and Pacific oceans. To examine the variability and trend of sea level around the peninsula, tide gauge records and satellite altimetry are analyzed taking into account vertical land movements (VLMs). At annual scale, sea level anomalies (SLAs) around Peninsular Malaysia on the order of 5–25 cm are mainly monsoon driven. Sea levels at eastern and western coasts respond differently to the Asian monsoon: two peaks per year in the Malacca Strait due to South Asian–Indian monsoon; an annual cycle in the remaining region mostly due to the East Asian–western Pacific monsoon. At interannual scale, regional sea level variability in the range of ±6 cm is correlated with El Niño–Southern Oscillation (ENSO). SLAs in the Malacca Strait side are further correlated with the Indian Ocean Dipole (IOD) in the range of ±5 cm. Interannual regional sea level falls are associated with El Niño events and positive phases of IOD, whilst rises are correlated with La Niña episodes and negative values of the IOD index. At seasonal to interannual scales, we observe the separation of the sea level patterns in the Singapore Strait, between the Raffles Lighthouse and Tanjong Pagar tide stations, likely caused by a dynamic constriction in the narrowest part. During the observation period 1986–2013, average relative rates of sea level rise derived from tide gauges in Malacca Strait and along the east coast of the peninsula are 3.6±1.6 and 3.7±1.1 mm yr −1 , respectively. Correcting for respective VLMs (0.8±2.6 and 0.9±2.2 mm yr −1 ), their corresponding geocentric sea level rise rates are estimated at 4.4±3.1 and 4.6±2.5 mm yr −1 . The geocentric rates are about 25 % faster than those measured at tide gauges around the peninsula; however, the level of uncertainty associated with VLM data is relatively high. For the common period between 1993 and 2009, geocentric sea level rise values along the Malaysian coast are similar from tide gauge records and satellite altimetry (3.1 and 2.7 mm yr −1 , respectively), and arguably correspond to the global trend.

Description: A case study on generation mechanisms of a sporadic sodium layer above Tromsø (69.6° N) during a night of high auroral activity Annales Geophysicae, 33, 941-953, 2015 Author(s): T. Takahashi, S. Nozawa, T. T. Tsuda, Y. Ogawa, N. Saito, T. Hidemori, T. D. Kawahara, C. Hall, H. Fujiwara, N. Matuura, A. Brekke, M. Tsutsumi, S. Wada, T. Kawabata, S. Oyama, and R. Fujii We have quantitatively evaluated generation mechanisms of a sporadic sodium layer (SSL) based on observational data obtained by multiple instruments at a high-latitude station: Ramfjordmoen, Tromsø, Norway (69.6° N, 19.2° E). The sodium lidar observed an SSL at 21:18 UT on 22 January 2012. The SSL was observed for 18 min, with a maximum sodium density of about 1.9 × 10 10 m −3 at 93 km with a 1.1 km thickness. The European Incoherent Scatter (EISCAT) UHF radar observed a sporadic E layer (Es layer) above 90 km from 20:00 to 23:00 UT. After 20:00 UT, the Es layer gradually descended and reached 94 km at 21:18 UT when the SSL appeared at the same altitude. In this event, considering the abundance of sodium ions (10 % or less), the Es layer could provide only about 37 % or less of the sodium atoms to the SSL. We have investigated a temporal development of the normal sodium ion layer with a consideration of chemical reactions and the effect of the (southwestward) electric field using observational values of the neutral temperature, electron density, horizontal neutral wind, and electric field. This calculation has shown that those processes, including contributions of the Es layer, would provide about 88 % of sodium atoms of the SSL. The effects of meteor absorption and auroral particle sputtering appear to be less important. Therefore, we have concluded that the major source of the SSL was sodium ions in a normal sodium ion layer. Two processes – namely the downward transportation of sodium ions from a normal sodium ion layer due to the electric field and the additional supply of sodium ions from the Es layer under relatively high electron density conditions (i.e., in the Es layer) – played a major role in generating the SSL in this event. Furthermore, we have found that the SSL was located in a lower-temperature region and that the temperature inside the SSL did not show any remarkable temperature enhancements.

Description: Conceptual design of a measurement network of the global change Atmospheric Chemistry and Physics Discussions, 15, 21063-21093, 2015 Author(s): P. Hari, T. Petäjä, J. Bäck, V.-M. Kerminen, H. K. Lappalainen, T. Vihma, T. Laurila, Y. Viisanen, T. Vesala, and M. Kulmala The global environment is changing rapidly due to anthropogenic emissions and actions. Such activities modify aerosol and greenhouse gas concentrations in the atmosphere, leading to regional and global climate change and affecting e.g. food and fresh-water security, sustainable use of natural resources and even demography. Here we present a conceptual design of a global, hierarchical observation network that can provide tools and increased understanding to tackle the inter-connected environmental and societal challenges that we will face in the coming decades. The philosophy behind the conceptual design relies on physical conservation laws of mass, energy and momentum, as well as on concentration gradients that act as driving forces for the atmosphere-biosphere exchange. The network is composed of standard, flux/advanced and flagship stations, each of which having specific and identified tasks. Each ecosystem type on the globe has its own characteristic features that have to be taken into consideration. The hierarchical network as a whole is able to tackle problems related to large spatial scales, heterogeneity of ecosystems and their complexity. The most comprehensive observations are envisioned to occur in flag ship stations, with which the process-level understanding can be expanded to continental and global scales together with advanced data analysis, earth system modelling and satellite remote sensing. The denser network of the flux and standard stations allow application and up-scaling of the results obtained from flag ship stations to the global level.

Description: Sensitivities of UK PM 2.5 concentrations to emissions reductions Atmospheric Chemistry and Physics Discussions, 15, 20881-20910, 2015 Author(s): M. Vieno, M. R. Heal, M. L. Williams, E. J. Carnell, J. R. Stedman, and S. Reis The reduction of ambient concentrations of fine particulate matter (PM 2.5 ) is a key objective for air pollution control policies in the UK and elsewhere. Long-term exposure to PM 2.5 has been identified as a major contributor to adverse human health effects in epidemiological studies and underpins ambient PM 2.5 legislation. As a range of emission sources and atmospheric chemistry transport processes contribute to PM 2.5 concentrations, atmospheric chemistry transport models are an essential tool to assess emissions control effectiveness. The EMEP4UK atmospheric chemistry transport model was used to investigate the impact of reductions in UK anthropogenic emissions of primary PM 2.5 , NH 3 , NO x , SO x or non-methane VOC on surface concentrations of PM 2.5 in the UK for a recent year (2010) and for a future current legislation emission scenario (2030). In general, the sensitivity to UK mitigation is rather small. A 30 % reduction in UK emissions of any one of the above components yields (for the 2010 simulation) a maximum reduction in PM 2.5 in any given location of ~ 0.6 μg m −3 (equivalent to ~ 6 % of the modelled PM 2.5 ). On average across the UK, the sensitivity of PM 2.5 concentrations to a 30 % reduction in UK emissions of individual contributing components, for both the 2010 and 2030 CLE baselines, increases in the order NMVOC, NO x , SO x , NH 3 and primary PM 2.5 , but there are strong spatial differences in the PM 2.5 sensitivities across the UK. Consequently, the sensitivity of PM 2.5 to individual component emissions reductions varies between area and population weighting. Reductions in NH 3 have the greatest effect on area-weighted PM 2.5 . A full UK population weighting places greater emphasis on reductions of primary PM 2.5 emissions, which is simulated to be the most effective single-component control on PM 2.5 for the 2030 scenario. An important observation is that weighting corresponding to the Average Exposure Indicator metric (using data from the 45 model grids containing a monitor whose measurements are used to calculate the UK AEI) further increases the emphasis on the effectiveness of primary PM 2.5 emissions reductions (and of NO x emissions reductions) relative to the effectiveness of NH 3 emissions reductions. Reductions in primary PM 2.5 have the largest impact on the AEI in both 2010 and the 2030 CLE scenario. The summation of the modelled reductions to the UK PM 2.5 AEI from 30 % reductions in UK emissions of primary PM 2.5 , NH 3 , SO x , NO x and VOC totals 1.17 and 0.82 μg m −3 for the 2010 and 2030 CLE simulations, respectively.

Description: A dual, single detector relaxed eddy accumulation system for long-term measurement of mercury flux Atmospheric Measurement Techniques Discussions, 8, 8113-8156, 2015 Author(s): S. Osterwalder, J. Fritsche, M. B. Nilsson, C. Alewell, J. Sommar, G. Jocher, M. Schmutz, J. Rinne, and K. Bishop The fate of anthropogenic emissions of mercury (Hg) to the atmosphere is influenced by the exchange of elemental Hg with the earth surface. This exchange which holds the key to a better understanding of Hg cycling from local to global scales has been difficult to quantify. To advance and facilitate research about land–atmosphere Hg interactions, we developed a dual-intake, single analyzer Relaxed Eddy Accumulation (REA) system. REA is an established technique for measuring turbulent fluxes of trace gases and aerosol particles in the atmospheric surface layer. Accurate determination of gaseous elemental mercury (GEM) fluxes has proven difficult to technical challenges presented by extremely small concentration differences (typically 〈 0.5 ng m −3 ) between updrafts and downdrafts. To address this we present an advanced REA design that uses two inlets and two pair of gold cartridges for semi-continuous monitoring of GEM fluxes. They are then analyzed sequentially on the same detector while another pair of gold cartridges takes over the sample collection. We also added a reference gas module for repeated quality-control measurements. To demonstrate the system performance, we present results from field campaigns in two contrasting environments: an urban setting with a heterogeneous fetch and a boreal mire during snow-melt. The observed emission rates were 15 and 3 ng m −2 h −1 . We claim that this dual-inlet, single detector approach is a significant development of the REA system for ultra-trace gases and can help to advance our understanding of long-term land–atmosphere GEM exchange.

Description: Perceived temperature in the course of climate change: an analysis of global heat index from 1979 to 2013 Earth System Science Data, 7, 193-202, 2015 Author(s): D. Lee and T. Brenner The increase in global mean temperatures resulting from climate change has wide reaching consequences for the earth's ecosystems and other natural systems. Many studies have been devoted to evaluating the distribution and effects of these changes. We go a step further and propose the use of the heat index, a measure of the temperature as perceived by humans, to evaluate global changes. The heat index, which is computed from temperature and relative humidity, is more important than temperature for the health of humans and animals. Even in cases where the heat index does not reach dangerous levels from a health perspective, it has been shown to be an important factor in worker productivity and thus in economic productivity. We compute the heat index from dew point temperature and absolute temperature 2 m above ground from the ERA-Interim reanalysis data set for the years 1979–2013. The described data set provides global heat index aggregated to daily minima, means and maxima per day ( doi:10.1594/PANGAEA.841057 ). This paper examines these data, as well as showing aggregations to monthly and yearly values. Furthermore, the data are spatially aggregated to the level of countries after being weighted by population density in order to facilitate the analysis of its impact on human health and productivity. The resulting data deliver insights into the spatiotemporal development of near-ground heat index during the course of the past three decades. It is shown that the impact of changing heat index is unevenly distributed through space and time, affecting some areas differently than others. The data can serve as a basis for evaluating and understanding the evolution of heat index in the course of climate change, as well as its impact on human health and productivity.

Description: Multiscale evaluation of the standardized precipitation index as a groundwater drought indicator Hydrology and Earth System Sciences Discussions, 12, 7405-7436, 2015 Author(s): R. Kumar, J. L. Musuuza, A. F. Van Loon, A. J. Teuling, R. Barthel, J. Ten Broek, J. Mai, L. Samaniego, and S. Attinger The lack of comprehensive groundwater observations at regional and global scales has promoted the use of alternative proxies and indices to quantify and predict groundwater droughts. Among them, the Standardized Precipitation Index (SPI) is commonly used to characterize droughts in different compartments of the hydro-meteorological system. In this study, we explore the suitability of the SPI to characterize local and regional scale groundwater droughts using observations at more than 2000 groundwater wells in geologically different areas in Germany and the Netherlands. A multiscale evaluation of the SPI is performed using the station data and their corresponding 0.5° gridded estimates to analyze the local and regional behavior of groundwater droughts, respectively. The standardized anomalies in the groundwater heads (SGI) were correlated against SPIs obtained using different accumulation periods. The accumulation periods to achieve maximum correlation exhibited high spatial variability (ranges 3 to 36 months) at both scales, leading to the conclusion that an a priori selection of the accumulation period (for computing the SPI) would result in inadequate characterization of groundwater droughts. The application of the uniform accumulation periods over the entire domain significantly reduced the correlation between SPI and SGI (≈ 21–66 %) indicating the limited applicability of SPI as a proxy for groundwater droughts even at long accumulation times. Furthermore, the low scores of the hit rate (0.3–0.6) and high false alarm ratio (0.4–0.7) at the majority of the wells and grid cells demonstrated the low reliability of groundwater drought predictions using the SPI. The findings of this study highlight the pitfalls of using the SPI as a groundwater drought indicator at both local and regional scales, and stress the need for more groundwater observations and accounting for regional hydrogeological characteristics in groundwater drought monitoring.

Description: Comparing the Ensemble and Extended Kalman Filters for in situ soil moisture assimilation with contrasting soil conditions Hydrology and Earth System Sciences Discussions, 12, 7353-7403, 2015 Author(s): D. Fairbairn, A. L. Barbu, J.-F. Mahfouf, J.-C. Calvet, and E. Gelati Two data assimilation methods are compared for their ability to produce a deterministic soil moisture analysis on the Météo-France land surface model: (i) SEKF, a Simplified Extended Kalman Filter, which uses a climatological background-error covariance, (ii) EnSRF, the Ensemble Square Root Filter, which uses an ensemble background-error covariance and approximates random forcing errors stochastically. The accuracy of the deterministic analysis is measured on 12 sites with in situ observations and various soil textures in Southwest France (SMOSMANIA network). In the experiments with real observations, the two methods perform similarly and improve on the open loop. Both methods suffer from incorrect linear assumptions which are particularly degrading to the analysis during water-stressed conditions: the EnSRF by a dry bias and the SEKF by an over-sensitivity of the model Jacobian between the surface and the root zone layers. These problems are less severe for sandy soils than clay soils because sandy soils are less sensitive to perturbations in the initial conditions. A simple bias correction technique is tested on the EnSRF. Although this reduces the bias, it also suppresses the ensemble spread, which degrades the analysis performance. However, the EnSRF flow-dependent background-error covariance evidently captures seasonal variability in the soil moisture errors and should exploit planned improvements in the model physics. Synthetic experiments demonstrate that when there is only a random component in the precipitation forcing errors, the correct stochastic representation of these errors enables the EnSRF to perform better than the SEKF. But in the real experiments the same rainfall error specification does not improve the EnSRF analysis. It is likely that the actual rainfall errors are underestimated and that other sources of errors could limit the usefulness of this information. More comprehensive ways of representing the rainfall errors are suggested, which might improve the EnSRF performance.

Description: Patterns in atmospheric carbonaceous aerosols in China: emission estimates and observed concentrations Atmospheric Chemistry and Physics, 15, 8657-8678, 2015 Author(s): H. Cui, P. Mao, Y. Zhao, C. P. Nielsen, and J. Zhang China is experiencing severe carbonaceous aerosol pollution driven mainly by large emissions resulting from intensive use of solid fuels. To gain a better understanding of the levels and trends of carbonaceous aerosol emissions and the resulting ambient concentrations at the national scale, we update an emission inventory of anthropogenic organic carbon (OC) and elemental carbon (EC) and employ existing observational studies to analyze characteristics of these aerosols including temporal, spatial, and size distributions, and the levels and shares of secondary organic carbon (SOC) in total OC. We further use ground observations to test the levels and inter-annual trends of the calculated national and provincial emissions of carbonaceous aerosols, and propose possible improvements in emission estimation for the future. The national OC emissions are estimated to have increased 29 % from 2000 (2127 Gg) to 2012 (2749 Gg) and EC by 37 % (from 1356 to 1857 Gg). The residential, industrial, and transportation sectors contributed an estimated 74–78, 17–21, and 4–6 % of the total emissions of OC, respectively, and 49–55, 30–34, and 14–18 % of EC. Updated emission factors (EFs) based on the most recent local field measurements, particularly for biofuel stoves, led to considerably lower emissions of OC compared to previous inventories. Compiling observational data across the country, higher concentrations of OC and EC are found in northern and inland cities, while higher OC / EC ratios are found in southern sites, due to the joint effects of primary emissions and meteorology. Higher OC / EC ratios are estimated at rural and remote sites compared to urban ones, attributed to more emissions of OC from biofuel use, more biogenic emissions of volatile organic compound (VOC) precursors to SOC, and/or transport of aged aerosols. For most sites, higher concentrations of OC, EC, and SOC are observed in colder seasons, while SOC / OC is reduced, particularly at rural and remote sites, attributed partly to weaker atmospheric oxidation and SOC formation compared to summer. Enhanced SOC formation from oxidization and anthropogenic activities like biomass combustion is judged to have crucial effects on severe haze events characterized by high particle concentrations. Several observational studies indicate an increasing trend in ambient OC / EC (but not in OC or EC individually) from 2000 to 2010, confirming increased atmospheric oxidation of OC across the country. Combining the results of emission estimation and observations, the improvement over prior emission inventories is indicated by inter-annual comparisons and correlation analysis. It is also indicated, however, that the estimated growth in emissions might be faster than observed growth, and that some sources with high primary OC / EC, such as burning of biomass, are still underestimated. Further studies to determine changing EFs over time in the residential sector and to compare to other measurements, such as satellite observations, are thus suggested to improve understanding of the levels and trends of primary carbonaceous aerosol emissions in China.

Description: Morphological dynamics of an englacial channel Hydrology and Earth System Sciences Discussions, 12, 7615-7664, 2015 Author(s): G. Vatne and T. D. L. Irvine-Fynn Despite an interest in the hydraulic functioning of supraglacial and englacial channels over the last four decades, the processes and forms of such ice-bounded streams have remained poorly documented. Recent glaciological research has demonstrated the potential significance of so-called "cut and closure" streams, where englacial or subglacial flowpaths are created from the long-term incision of supraglacial channels. These flowpaths are reported to exhibit step-pool morphology, comprising knickpoints and/or knickzones, albeit exaggerated in dimensions in comparison to their supraglacial channel counterparts. However, little is known of the development of such channels' morphology. Here, we examine the spatial organization of step-pools and the upstream migration of steps, many of which form knickzones, with repeated surveys over a 10 year period in an englacial conduit in cold-based Austre Brøggerbreen, Svalbard. The observations show upstream knickpoint recession to be the dominant process for channel evolution. This is paralleled by an increase in average step height and conduit gradient over time. Characteristic channel reach types and step-riser forms are consistently observed in each of the morphological surveys reported. We suggest that the formation of steps has a hydrodynamic origin, where step-pool geometry is more efficient for energy dissipation than meanders, and that the englacial channel system is one in rapid transition rather than in dynamic equilibrium. The evolution and recession of knickzones reported here result in the formation of a 37 m moulin, suggesting over time the englacial channel may evolve towards a stable end-point characterised by a singular vertical descent to the local hydraulic base level. In light of this, our observations highlight the need to further examine the adjustment processes in cut-and-closure channels to better understand their coupling to supraglacial meltwater sources and role and potential significance in cold-based glacier hydrology and ice dynamics.

Description: Scenario based approach for multiple source Tsunami Hazard assessment for Sines, Portugal Natural Hazards and Earth System Sciences Discussions, 3, 4663-4693, 2015 Author(s): M. Wronna, R. Omira, and M. A. Baptista In this paper, we present a scenario-based approach for tsunami hazard assessment for the city and harbour of Sines – Portugal, one of the test-sites of project ASTARTE. Sines holds one of the most important deep-water ports which contains oil-bearing, petrochemical, liquid bulk, coal and container terminals. The port and its industrial infrastructures are facing the ocean southwest towards the main seismogenic sources. This work considers two different seismic zones: the Southwest Iberian Margin and the Gloria Fault. Within these two regions, we selected a total of six scenarios to assess the tsunami impact at the test site. The tsunami simulations are computed using NSWING a Non-linear Shallow Water Model With Nested Grids. In this study, the static effect of tides is analysed for three different tidal stages MLLW (mean lower low water), MSL (mean sea level) and MHHW (mean higher high water). For each scenario, inundation is described by maximum values of wave height, flow depth, drawback, runup and inundation distance. Synthetic waveforms are computed at virtual tide gauges at specific locations outside and inside the harbour. The final results describe the impact at Sines test site considering the single scenarios at mean sea level, the aggregate scenario and the influence of the tide on the aggregate scenario. The results confirm the composite of Horseshoe and Marques Pombal fault as the worst case scenario. It governs the aggregate scenario with about 60 % and inundates an area of 3.5 km 2 .

Description: A new high-precision and low-power GNSS receiver for long-term installations in remote areas Geoscientific Instrumentation, Methods and Data Systems Discussions, 5, 285-310, 2015 Author(s): D. H. Jones, C. Robinson, and G. H. Gudmundsson We have developed a new high-precision GNSS receiver specifically designed for long-term unattended deployments in remote areas. The receiver reports its status, and can be reprogrammed remotely, through an integrated satellite data link. It uses less power than commercially available alternatives while being equally, if not more, accurate. Data is saved locally on dual SD card slots for increased reliability. Deployments of a number those receivers in several different locations on the Antarctic Ice Sheet have shown them to robust and able to operate flawlessly at low temperatures down to −40 °C.

Description: The impact of near-surface soil moisture assimilation at subseasonal, seasonal, and inter-annual time scales Hydrology and Earth System Sciences Discussions, 12, 7971-8004, 2015 Author(s): C. Draper and R. Reichle Nine years of Advanced Microwave Scanning Radiometer – Earth Observing System (AMSR-E) soil moisture retrievals are assimilated into the Catchment land surface model at four locations in the US. The assimilation is evaluated using the unbiased Mean Square Error (ubMSE) relative to watershed-scale in situ observations, with the ubMSE separated into contributions from the subseasonal (SM short ), mean seasonal (SM seas ) and inter-annual (SM long ) soil moisture dynamics. For near-surface soil moisture, the average ubMSE for Catchment without assimilation was (1.8 × 10 −3 m 3 m −3 ) 2 , of which 19 % was in SM long , 26 % in SM seas , and 55 % in SM short . The AMSR-E assimilation significantly reduced the total ubMSE at every site, with an average reduction of 33 %. Of this ubMSE reduction, 37 % occurred in SM long , 24 % in SM seas , and 38 % in SM short . For root-zone soil moisture, in situ observations were available at one site only, and the near-surface and root-zone results were very similar at this site. These results suggest that, in addition to the well-reported improvements in SM short , assimilating a sufficiently long soil moisture data record can also improve the model representation of important long term events, such as droughts. The improved agreement between the modeled and in situ SM seas is harder to interpret, given that mean seasonal cycle errors are systematic, and systematic errors are not typically targeted by (bias-blind) data assimilation. Finally, the use of one year subsets of the AMSR-E and Catchment soil moisture for estimating the observation-bias correction (rescaling) parameters is investigated. It is concluded that when only one year of data is available, the associated uncertainty in the rescaling parameters should not greatly reduce the average benefit gained from data assimilation, but locally and in extreme years there is a risk of increased errors.

Description: Century-scale simulations of the response of the West Antarctic Ice Sheet to a warming climate The Cryosphere, 9, 1579-1600, 2015 Author(s): S. L. Cornford, D. F. Martin, A. J. Payne, E. G. Ng, A. M. Le Brocq, R. M. Gladstone, T. L. Edwards, S. R. Shannon, C. Agosta, M. R. van den Broeke, H. H. Hellmer, G. Krinner, S. R. M. Ligtenberg, R. Timmermann, and D. G. Vaughan We use the BISICLES adaptive mesh ice sheet model to carry out one, two, and three century simulations of the fast-flowing ice streams of the West Antarctic Ice Sheet, deploying sub-kilometer resolution around the grounding line since coarser resolution results in substantial underestimation of the response. Each of the simulations begins with a geometry and velocity close to present-day observations, and evolves according to variation in meteoric ice accumulation rates and oceanic ice shelf melt rates. Future changes in accumulation and melt rates range from no change, through anomalies computed by atmosphere and ocean models driven by the E1 and A1B emissions scenarios, to spatially uniform melt rate anomalies that remove most of the ice shelves over a few centuries. We find that variation in the resulting ice dynamics is dominated by the choice of initial conditions and ice shelf melt rate and mesh resolution, although ice accumulation affects the net change in volume above flotation to a similar degree. Given sufficient melt rates, we compute grounding line retreat over hundreds of kilometers in every major ice stream, but the ocean models do not predict such melt rates outside of the Amundsen Sea Embayment until after 2100. Within the Amundsen Sea Embayment the largest single source of variability is the onset of sustained retreat in Thwaites Glacier, which can triple the rate of eustatic sea level rise.

Description: Numerical study of upper hybrid to Z -mode leakage during electromagnetic pumping of groups of striations in the ionosphere Annales Geophysicae, 33, 1019-1030, 2015 Author(s): B. Eliasson and T. B. Leyser We investigate numerically the interaction between ionospheric magnetic field-aligned density striations and a left-hand circularly polarized ( L )-mode wave. The L -mode wave is scattered into upper hybrid (UH) waves which are partially trapped in the striations, but leak energy to electromagnetic waves in the Z -mode branch. For small-amplitude (1 %) striations, this loss mechanism leads to a significant reduction in amplitude of the UH waves. For several striations organized in a lattice, the leaking of Z -mode waves is compensated by influx of Z -mode radiation from neighboring striations, leading to an increased amplitude of the weakly trapped UH waves. For large-amplitude (10 %) striations the trapped UH waves rapidly increase in amplitude far beyond the threshold for parametric instabilities, and the Z -mode leakage is less important. The results have relevance for the growth of striations and the onset of UH and lower hybrid turbulence during electromagnetic high-frequency pumping of ionospheric plasma, which require large-amplitude UH waves.

Description: Observation of a new type of low-frequency waves at comet 67P/Churyumov-Gerasimenko Annales Geophysicae, 33, 1031-1036, 2015 Author(s): I. Richter, C. Koenders, H.-U. Auster, D. Frühauff, C. Götz, P. Heinisch, C. Perschke, U. Motschmann, B. Stoll, K. Altwegg, J. Burch, C. Carr, E. Cupido, A. Eriksson, P. Henri, R. Goldstein, J.-P. Lebreton, P. Mokashi, Z. Nemeth, H. Nilsson, M. Rubin, K. Szegö, B. T. Tsurutani, C. Vallat, M. Volwerk, and K.-H. Glassmeier We report on magnetic field measurements made in the innermost coma of 67P/Churyumov-Gerasimenko in its low-activity state. Quasi-coherent, large-amplitude (δ B / B ~ 1), compressional magnetic field oscillations at ~ 40 mHz dominate the immediate plasma environment of the nucleus. This differs from previously studied cometary interaction regions where waves at the cometary ion gyro-frequencies are the main feature. Thus classical pickup-ion-driven instabilities are unable to explain the observations. We propose a cross-field current instability associated with newborn cometary ion currents as a possible source mechanism.

Description: Hydroxy fatty acids in fresh snow samples from northern Japan: long-range atmospheric transport of Gram-negative bacteria by Asian winter monsoon Biogeosciences Discussions, 12, 13375-13397, 2015 Author(s): P. Tyagi, S. Yamamoto, and K. Kawamura Hydroxy fatty acids (FAs) in fresh snow from Sapporo, one of the heaviest snowfall regions in the world, have been studied to ascertain the airborne bacterial endotoxin concentrations and their biomass. The presence of β-hydroxy FAs (C 9 –C 28 ), constituents of Gram-negative bacteria (GNB), suggests long-range transport of soil microbes. Likewise, the occurrence of α- and ω-hydroxy FAs (C 9 –C 30 and C 9 –C 28 , respectively) in snow reveals their contribution from epicuticular waxes and soil microorganisms. Estimated endotoxin and GNB mass can aid in assessing their possible impacts on the diversity and functioning of aquatic and terrestrial ecosystems, as well as lethal effects on pedestrians through dispersal of microbes. Air mass back trajectories together with hydroxy FAs unveil their sources from Siberia, Russian Far East and North China by the Asian monsoon. This study highlights the role of fresh snow that reduces the human health risk of GNB and endotoxin by scavenging from the air.