ALBERT

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.

Description: A method for predicting the factor of safety of an infinite slope based on the depth ratio of the wetting front induced by rainfall infiltration Natural Hazards and Earth System Science, 15, 1835-1849, 2015 Author(s): B.-G. Chae, J.-H. Lee, H.-J. Park, and J. Choi Most landslides in Korea are classified as shallow landslides with an average depth of less than 2 m. These shallow landslides are associated with the advance of a wetting front in the unsaturated soil due to rainfall infiltration, which results in an increase in water content and a reduction in the matric suction in the soil. Therefore, this study presents a modified equation of infinite slope stability analysis based on the concept of the saturation depth ratio to analyze the slope stability change associated with the rainfall on a slope. A rainfall infiltration test in unsaturated soil was performed using a column to develop an understanding of the effect of the saturation depth ratio following rainfall infiltration. The results indicated that the rainfall infiltration velocity due to the increase in rainfall in the soil layer was faster when the rainfall intensity increased. In addition, the rainfall infiltration velocity tends to decrease with increases in the unit weight of soil. The proposed model was applied to assess its feasibility and to develop a regional landslide susceptibility map using a geographic information system (GIS). For that purpose, spatial databases for input parameters were constructed and landslide locations were obtained. In order to validate the proposed approach, the results of the proposed approach were compared with the landslide inventory using a ROC (receiver operating characteristics) graph. In addition, the results of the proposed approach were compared with the previous approach used: a steady-state hydrological model. Consequently, the approach proposed in this study displayed satisfactory performance in classifying landslide susceptibility and showed better performance than the steady-state approach.

Description: Regional trends and controlling factors of fatal landslides in Latin America and the Caribbean Natural Hazards and Earth System Science, 15, 1821-1833, 2015 Author(s): S. A. Sepúlveda and D. N. Petley A new data set of landslides that caused loss of life in Latin America and the Caribbean in the 10-year period from 2004 and 2013 inclusive has been compiled, providing new insight into the impact of landslides in this key part of the world. This data set indicates that in the 10-year period a total of 11 631 people lost their lives across the region in 611 landslides. The geographical distribution of the landslides is highly heterogeneous, with areas of high incidence in parts of the Caribbean (most notably Haiti), Central America, Colombia, and southeast Brazil. There is significant interannual variation in the number of landslides, with the El Niño/La Niña cycle emerging as a key control. Our analysis suggests that on a continental scale the mapped factors that best explain the observed distribution are topography, annual precipitation and population density. On a national basis we have compared the occurrence of fatality-inducing landslide occurrence with the production of locally authored research articles, demonstrating that there is a landslide research deficit in Latin America and the Caribbean. Understanding better the mechanisms, distribution causes and triggers of landslides in Latin America and the Caribbean must be an essential first step towards managing the hazard.

Description: Spatial scales of temperature and salinity variability estimated from Argo observations Ocean Science Discussions, 12, 1793-1814, 2015 Author(s): F. Ninove, P. Y. Le Traon, E. Remy, and S. Guinehut Argo observations from 2005 to 2013 are used to characterize spatial scales temperature and salinity variations from the surface down to 1500 m. Simulations are first performed to analyze the sensitivity of results to Argo sampling; they show that several years of Argo observations are required to estimate the spatial scales of ocean variability over 20° × 20° boxes. Spatial scales are then computed over several large scale areas. Zonal and meridional spatial scales ( Lx and Ly which are also zero crossing of covariance functions) vary as expected with latitudes. Scales are of about 100 km at high latitudes and more of 700 km in the Indian and Pacific equatorial/tropical regions. Zonal and meridional scales are similar: except in these tropical/equatorial regions where zonal scales are much larger (by a factor of 2 to 3) than meridional scales. Spatial scales are the largest close to the surface and have a general tendency for temperature to increase in deeper layers. There are significant differences between temperature and salinity scales, in particular, in the deep ocean. Results are consistent with previous studies based on sparse in-situ observations or satellite altimetry. They provide, however, for the first time a global description of temperature and salinity scales of variability and a characterization of their variations according to depths.

Description: Evaluation of an operational ocean model configuration at 1/12° spatial resolution for the Indonesian seas – Part 2: Biogeochemistry Geoscientific Model Development Discussions, 8, 6669-6706, 2015 Author(s): E. Gutknecht, G. Reffray, M. Gehlen, I. Triyulianti, D. Berlianty, and P. Gaspar In the framework of the INDESO (Infrastructure evelopment of Space Oceanography) project, an operational ocean forecasting system was developed to monitor the state of the Indonesian seas in terms of circulation, biogeochemistry and fisheries. This forecasting system combines a suite of numerical models connecting physical and biogeochemical variables to population dynamics of large marine predators (tunas). The physical/biogeochemical coupled component (INDO12BIO configuration) covers a large region extending from the western Pacific Ocean to the Eastern Indian Ocean at 1/12° resolution. The OPA/NEMO physical ocean model and the PISCES biogeochemical model are coupled in "on-line" mode without degradation in space and time. The operational global ocean forecasting system (1/4°) operated by Mercator Ocean provides the physical forcing while climatological open boundary conditions are prescribed for the biogeochemistry. This paper describes the skill assessment of the INDO12BIO configuration. Model skill is assessed by evaluating a reference hindcast simulation covering the last 8 years (2007–2014). Model results are compared to satellite, climatological and in situ observations. Diagnostics are performed on chlorophyll a , primary production, mesozooplankton, nutrients and oxygen. Model results reproduce the main characteristics of biogeochemical tracer distributions in space and time. The seasonal cycle of chlorophyll a is in phase with satellite observations. The northern and southern parts of the archipelago present a distinct seasonal cycle, with higher chlorophyll biomass in the southern (northern) part during SE (NW) monsoon. Nutrient and oxygen concentrations are correctly reproduced in terms of horizontal and vertical distributions. The biogeochemical content of water masses entering in the archipelago as well as the water mass transformation across the archipelago conserves realistic vertical distribution in Banda Sea and at the exit of the archipelago.

Description: A numerical study of the early stages of a tropical cyclogenesis in relation to the MJO Natural Hazards and Earth System Sciences Discussions, 3, 4919-4935, 2015 Author(s): J. Guerbette, M. Plu, C. Barthe, and J.-F. Mahfouf The role of an active phase of the Madden–Julian Oscillation (MJO) on the evolution of a mesoscale convective systems (MCS) leading to a tropical depression is investigated in the South-West Indian Ocean during the Dynamics of the Madden–Julian Oscillation (DYNAMO) field experiment, with a numerical limited-area atmospheric model. A mesoscale vortex is followed in the low-troposphere from the initiation of the active MJO phase. It is shown that the interaction of the vortex with the Equatorial jet associated with the MJO plays an important role on the vortex development. As the vortex encounters the southern part of the low-level jet, it undergoes intensification that is explained by the barotropic conversion of kinetic energy from the low-level jet to the vortex.

Description: A semi-analytical model for diffuse reflectance in marine and inland waters Ocean Science Discussions, 12, 1893-1912, 2015 Author(s): J. D. Pravin, P. Shanmugam, and Y.-H. Ahn A semi-analytical model for predicting diffuse reflectance of coastal and oceanic waters is developed based on the water-column optical properties and illumination conditions. Diffuse reflectance ( R ) is an apparent optical property that is related to the Gordon's parameter ( b b /(a+b b ) ) through a proportionality factor " f ". The conventional assumption of " f " as a constant (0.33) yields large errors in case of turbid and productive coastal waters and a predictive model based on this assumption is generally restricted to open-ocean waters (low chlorophyll case). In this paper, we have sorted the dependent factors that influence " f " values in the water column. Here, the parameter " f " is modeled as a function of wavelength, depth, inherent optical properties (IOPs) and illumination conditions. This work eliminates the spectral constants ( K Chl and K SS ) associated with our previous model and constrains the present model to be solely dependent on the IOPs and illumination conditions. Data used for parameterization and validation are obtained from in situ measurements in different waters within coastal environments. Validation shows good agreement between the model R and in situ R values with the overall mean relative error of less than a few percent. The model is valid for a wide range waters within coastal and open-ocean environments.

Description: Impact of variable sea-water conductivity on motional induction simulated with an OGCM Ocean Science Discussions, 12, 1869-1891, 2015 Author(s): C. Irrgang, J. Saynisch, and M. Thomas Carrying high concentrations of dissolved salt, ocean water is a good electrical conductor. As sea-water flows through the Earth's ambient geomagnetic field, electric fields are generated, which in turn induce secondary magnetic fields. In current models for oceanic induced magnetic fields, a realistic consideration of sea-water conductivity is often neglected and the effect on the variability of the oceanic induced magnetic field unknown. To model magnetic fields that are induced by non-tidal global ocean currents, an electromagnetic induction model is implemented into the Ocean Model for Circulation and Tides (OMCT). This provides the opportunity to not only model oceanic induced magnetic signals, but to assess the impact of oceanographic phenomena on the induction process. In this paper, the sensitivity of the induction process due to spatial and temporal variations in sea-water conductivity is investigated. It is shown that assuming an ocean-wide uniform conductivity is insufficient to accurately capture the temporal variability of the magnetic signal. Using instead a realistic global sea-water conductivity distribution increases the temporal variability of the magnetic field up to 45 %. Especially vertical gradients in sea-water conductivity prove to be a key factor for the variability of the oceanic induced magnetic field. However, temporal variations of sea-water conductivity only marginally affect the magnetic signal.

Description: A new model of the global biogeochemical cycle of carbonyl sulfide – Part 2: Use of carbonyl sulfide to constrain gross primary productivity in current vegetation models Atmospheric Chemistry and Physics, 15, 9285-9312, 2015 Author(s): T. Launois, P. Peylin, S. Belviso, and B. Poulter Clear analogies between carbonyl sulfide (OCS) and carbon dioxide (CO 2 ) diffusion pathways through leaves have been revealed by experimental studies, with plant uptake playing an important role for the atmospheric budget of both species. Here we use atmospheric OCS to evaluate the gross primary production (GPP) of three dynamic global vegetation models (Lund–Potsdam–Jena, LPJ; National Center for Atmospheric Research – Community Land Model 4, NCAR-CLM4; and Organising Carbon and Hydrology In Dynamic Ecosystems, ORCHIDEE). Vegetation uptake of OCS is modeled as a linear function of GPP and leaf relative uptake (LRU), the ratio of OCS to CO 2 deposition velocities of plants. New parameterizations for the non-photosynthetic sinks (oxic soils, atmospheric oxidation) and biogenic sources (oceans and anoxic soils) of OCS are also provided. Despite new large oceanic emissions, global OCS budgets created with each vegetation model show exceeding sinks by several hundred Gg S yr −1 . An inversion of the surface fluxes (optimization of a global scalar which accounts for flux uncertainties) led to balanced OCS global budgets, as atmospheric measurements suggest, mainly by drastic reduction (up to −50 %) in soil and vegetation uptakes. The amplitude of variations in atmospheric OCS mixing ratios is mainly dictated by the vegetation sink over the Northern Hemisphere. This allows for bias recognition in the GPP representations of the three selected models. The main bias patterns are (i) the terrestrial GPP of ORCHIDEE at high northern latitudes is currently overestimated, (ii) the seasonal variations of the GPP are out of phase in the NCAR-CLM4 model, showing a maximum carbon uptake too early in spring in the northernmost ecosystems, (iii) the overall amplitude of the seasonal variations of GPP in NCAR-CLM4 is too small, and (iv) for the LPJ model, the GPP is slightly out of phase for the northernmost ecosystems and the respiration fluxes might be too large in summer in the Northern Hemisphere. These results rely on the robustness of the OCS modeling framework and, in particular, the choice of the LRU values (assumed constant in time) and the parameterization of soil OCS uptake with small seasonal variations. Refined optimization with regional-scale and seasonally varying coefficients might help to test some of these hypothesis.

Description: Upper tropospheric water vapour variability at high latitudes – Part 1: Influence of the annular modes Atmospheric Chemistry and Physics Discussions, 15, 22291-22329, 2015 Author(s): C. E. Sioris, J. Zou, D. A. Plummer, C. D. Boone, C. T. McElroy, P. E. Sheese, O. Moeini, and P. F. Bernath Seasonal and monthly zonal medians of water vapour in the upper troposphere and lower stratosphere (UTLS) are calculated for both Atmospheric Chemistry Experiment (ACE) instruments for the northern and southern high-latitude regions (60–90 and 60–90° S). Chosen for the purpose of observing high-latitude processes, the ACE orbit provides sampling of both regions in eight of 12 months of the year, with coverage in all seasons. The ACE water vapour sensors, namely MAESTRO (Measurements of Aerosol Extinction in the Stratosphere and Troposphere Retrieved by Occultation) and the Fourier Transform Spectrometer (ACE-FTS) are currently the only satellite instruments that can probe from the lower stratosphere down to the mid-troposphere to study the vertical profile of the response of UTLS water vapour to the annular modes. The Arctic oscillation (AO), also known as the northern annular mode (NAM), explains 64 % ( r = −0.80) of the monthly variability in water vapour at northern high-latitudes observed by ACE-MAESTRO between 5 and 7 km using only winter months (January to March 2004–2013). Using a seasonal timestep and all seasons, 45 % of the variability is explained by the AO at 6.5 ± 0.5 km, similar to the 46 % value obtained for southern high latitudes at 7.5 ± 0.5 km explained by the Antarctic oscillation or southern annular mode (SAM). A large negative AO event in March 2013 produced the largest relative water vapour anomaly at 5.5 km (+70 %) over the ACE record. A similarly large event in the 2010 boreal winter, which was the largest negative AO event in the record (1950–2015), led to 〉 50 % increases in water vapour observed by MAESTRO and ACE-FTS at 7.5 km.

Description: Decadal evaluation of regional climate, air quality, and their interactions using WRF/Chem Version 3.6.1 Geoscientific Model Development Discussions, 8, 6707-6756, 2015 Author(s): K. Yahya, K. Wang, P. Campbell, T. Glotfelty, J. He, and Y. Zhang The Weather Research and Forecasting model with Chemistry (WRF/Chem) v3.6.1 with the Carbon Bond 2005 (CB05) gas-phase mechanism is evaluated for its first decadal application during 2001–2010 using the Representative Concentration Pathway (RCP 8.5) emissions to assess its capability and appropriateness for long-term climatological simulations. The initial and boundary conditions are downscaled from the modified Community Earth System Model/Community Atmosphere Model (CESM/CAM5) v1.2.2. The meteorological initial and boundary conditions are bias-corrected using the National Center for Environmental Protection's Final (FNL) Operational Global Analysis data. Climatological evaluations are carried out for meteorological, chemical, and aerosol-cloud-radiation variables against data from surface networks and satellite retrievals. The model performs very well for the 2 m temperature (T2) for the 10 year period with only a small cold bias of −0.3 °C. Biases in other meteorological variables including relative humidity at 2 m, wind speed at 10 m, and precipitation tend to be site- and season-specific; however, with the exception of T2, consistent annual biases exist for most of the years from 2001 to 2010. Ozone mixing ratios are slightly overpredicted at both urban and rural locations but underpredicted at rural locations. PM 2.5 concentrations are slightly overpredicted at rural sites, but slightly underpredicted at urban/suburban sites. In general, the model performs relatively well for chemical and meteorological variables, and not as well for aerosol-cloud-radiation variables. Cloud-aerosol variables including aerosol optical depth, cloud water path, cloud optical thickness, and cloud droplet number concentration are generally underpredicted on average across the continental US. Overpredictions of several cloud variables over eastern US result in underpredictions of radiation variables and overpredictions of shortwave and longwave cloud forcing which are important climate variables. While the current performance is deemed to be acceptable, improvements to the bias-correction method for CESM downscaling and the model parameterizations of cloud dynamics and thermodynamics, as well as aerosol-cloud interactions can potentially improve model performance for long-term climate simulations.

Description: Retrieving the paleoclimatic signal from the deeper part of the EPICA Dome C ice core The Cryosphere, 9, 1633-1648, 2015 Author(s): J.-L. Tison, M. de Angelis, G. Littot, E. Wolff, H. Fischer, M. Hansson, M. Bigler, R. Udisti, A. Wegner, J. Jouzel, B. Stenni, S. Johnsen, V. Masson-Delmotte, A. Landais, V. Lipenkov, L. Loulergue, J.-M. Barnola, J.-R. Petit, B. Delmonte, G. Dreyfus, D. Dahl-Jensen, G. Durand, B. Bereiter, A. Schilt, R. Spahni, K. Pol, R. Lorrain, R. Souchez, and D. Samyn An important share of paleoclimatic information is buried within the lowermost layers of deep ice cores. Because improving our records further back in time is one of the main challenges in the near future, it is essential to judge how deep these records remain unaltered, since the proximity of the bedrock is likely to interfere both with the recorded temporal sequence and the ice properties. In this paper, we present a multiparametric study (δD-δ 18 O ice , δ 18 O atm , total air content, CO 2 , CH 4 , N 2 O, dust, high-resolution chemistry, ice texture) of the bottom 60 m of the EPICA (European Project for Ice Coring in Antarctica) Dome C ice core from central Antarctica. These bottom layers were subdivided into two distinct facies: the lower 12 m showing visible solid inclusions (basal dispersed ice facies) and the upper 48 m, which we will refer to as the "basal clean ice facies". Some of the data are consistent with a pristine paleoclimatic signal, others show clear anomalies. It is demonstrated that neither large-scale bottom refreezing of subglacial water, nor mixing (be it internal or with a local basal end term from a previous/initial ice sheet configuration) can explain the observed bottom-ice properties. We focus on the high-resolution chemical profiles and on the available remote sensing data on the subglacial topography of the site to propose a mechanism by which relative stretching of the bottom-ice sheet layers is made possible, due to the progressively confining effect of subglacial valley sides. This stress field change, combined with bottom-ice temperature close to the pressure melting point, induces accelerated migration recrystallization, which results in spatial chemical sorting of the impurities, depending on their state (dissolved vs. solid) and if they are involved or not in salt formation. This chemical sorting effect is responsible for the progressive build-up of the visible solid aggregates that therefore mainly originate "from within", and not from incorporation processes of debris from the ice sheet's substrate. We further discuss how the proposed mechanism is compatible with the other ice properties described. We conclude that the paleoclimatic signal is only marginally affected in terms of global ice properties at the bottom of EPICA Dome C, but that the timescale was considerably distorted by mechanical stretching of MIS20 due to the increasing influence of the subglacial topography, a process that might have started well above the bottom ice. A clear paleoclimatic signal can therefore not be inferred from the deeper part of the EPICA Dome C ice core. Our work suggests that the existence of a flat monotonic ice–bedrock interface, extending for several times the ice thickness, would be a crucial factor in choosing a future "oldest ice" drilling location in Antarctica.

Description: Impacts of an unknown daytime HONO source on the mixing ratio and budget of HONO, and hydroxyl, hydroperoxyl, and organic peroxy radicals, in the coastal regions of China Atmospheric Chemistry and Physics, 15, 9381-9398, 2015 Author(s): Y. Tang, J. An, F. Wang, Y. Li, Y. Qu, Y. Chen, and J. Lin Many field experiments have found high nitrous acid (HONO) mixing ratios in both urban and rural areas during daytime, but these high daytime HONO mixing ratios cannot be explained well by gas-phase production, HONO emissions, and nighttime hydrolysis conversion of nitrogen dioxide (NO 2 ) on aerosols, suggesting that an unknown daytime HONO source ( P unknown ) could exist. The formula P unknown ≈ 19.60[NO 2 ] · J (NO 2 ) was obtained using observed data from 13 field experiments across the globe. The three additional HONO sources (i.e., the P unknown , nighttime hydrolysis conversion of NO 2 on aerosols, and HONO emissions) were coupled into the WRF-Chem model (Weather Research and Forecasting model coupled with Chemistry) to assess the P unknown impacts on the concentrations and budgets of HONO and peroxy (hydroxyl, hydroperoxyl, and organic peroxy) radicals (RO x ) (= OH + HO 2 + RO 2 ) in the coastal regions of China. Results indicated that the additional HONO sources produced a significant improvement in HONO and OH simulations, particularly in the daytime. High daytime average P unknown values were found in the coastal regions of China, with a maximum of 2.5 ppb h −1 in the Beijing–Tianjin–Hebei region. The P unknown produced a 60–250 % increase of OH, HO 2 , and RO 2 near the ground in the major cities of the coastal regions of China, and a 5–48 % increase of OH, HO 2 , and RO 2 in the daytime meridional-mean mixing ratios within 1000 m above the ground. When the three additional HONO sources were included, the photolysis of HONO was the second most important source in the OH production rate in Beijing, Shanghai, and Guangzhou before 10:00 LST with a maximum of 3.72 ppb h −1 and a corresponding P unknown contribution of 3.06 ppb h −1 in Beijing, whereas the reaction of HO 2 + NO (nitric oxide) was dominant after 10:00 LST with a maximum of 9.38 ppb h −1 and a corresponding P unknown contribution of 7.23 ppb h −1 in Beijing. The whole RO x cycle was accelerated by the three additional HONO sources, especially the P unknown . The daytime average OH production rate was enhanced by 0.67 due to the three additional HONO sources; [0.64], due to the P unknown , to 4.32 [3.86] ppb h −1 , via the reaction of HO 2 + NO, and by 0.49 [0.47] to 1.86 [1.86] ppb h −1 , via the photolysis of HONO. The OH daytime average loss rate was enhanced by 0.58 [0.55] to 2.03 [1.92] ppb h −1 , via the reaction of OH + NO 2 , and by 0.31 [0.28] to 1.78 [1.64] ppb h −1 , via the reaction of OH + CO (carbon monoxide) in Beijing, Shanghai, and Guangzhou. Similarly, the three additional HONO sources produced an increase of 0.31 [0.28] (with a corresponding P unknown contribution) to 1.78 [1.64] ppb h −1 , via the reaction of OH + CO, and 0.10 [0.09] to 0.63 [0.59] ppb h −1 , via the reaction of CH 3 O 2 (methylperoxy radical) + NO in the daytime average HO 2 production rate, and 0.67 [0.61] to 4.32 [4.27] ppb h −1 , via the reaction of HO 2 + NO in the daytime average HO 2 loss rate in Beijing, Shanghai, and Guangzhou. The above results suggest that the P unknown considerably enhanced the RO x concentrations and accelerated RO x cycles in the coastal regions of China, and could produce significant increases in concentrations of inorganic aerosols and secondary organic aerosols and further aggravate haze events in these regions.

Description: A modification to the standard ionospheric correction method used in GPS radio occultation Atmospheric Measurement Techniques, 8, 3385-3393, 2015 Author(s): S. B. Healy and I. D. Culverwell A modification to the standard bending-angle correction used in GPS radio occultation (GPS-RO) is proposed. The modified approach should reduce systematic residual ionospheric errors in GPS radio occultation climatologies. A new second-order term is introduced in order to account for a known source of systematic error, which is generally neglected. The new term has the form κ(a) × (α L1 (a)-α L2 (a)) 2 , where a is the impact parameter and (α L1 , α L2 ) are the L1 and L2 bending angles, respectively. The variable κ is a weak function of the impact parameter, a , but it does depend on a priori ionospheric information. The theoretical basis of the new term is examined. The sensitivity of κ to the assumed ionospheric parameters is investigated in one-dimensional simulations, and it is shown that κ ≃ 10–20 rad −1 . We note that the current implicit assumption is κ=0, and this is probably adequate for numerical weather prediction applications. However, the uncertainty in κ should be included in the uncertainty estimates for the geophysical climatologies produced from GPS-RO measurements. The limitations of the new ionospheric correction when applied to CHAMP (Challenging Minisatellite Payload) measurements are noted. These arise because of the assumption that the refractive index is unity at the satellite, made when deriving bending angles from the Doppler shift values.

Description: Nitrogen cycling in the subsurface biosphere: nitrate isotopes in porewaters underlying the oligotrophic North Atlantic Biogeosciences Discussions, 12, 13545-13591, 2015 Author(s): S. D. Wankel, C. Buchwald, W. Ziebis, C. B. Wenk, and M. F. Lehmann Nitrogen (N) is a key component of fundamental biomolecules. Hence, the cycling and availability of N is a central factor governing the extent of ecosystems across the Earth. In the organic-lean sediment porewaters underlying the oligotrophic ocean, where low levels of microbial activity persist despite limited organic matter delivery from overlying water, the extent and modes of nitrogen transformations have not been widely investigated. Here we use the N and oxygen (O) isotopic composition of porewater nitrate (NO 3 − ) from a site in the oligotrophic North Atlantic (IODP) to determine the extent and magnitude of microbial nitrate production (via nitrification) and consumption (via denitrification). We find that NO 3 − accumulates far above bottom seawater concentrations (∼ 21 μM) throughout the sediment column (up to ∼ 50 μM) down to the oceanic basement as deep as 90 mbsf, reflecting the predominance of aerobic nitrification/remineralization within the deep marine sediments. Large changes in the δ 15 N and δ 18 O of nitrate, however, reveal variable influence of nitrate respiration across the three sites. We use an inverse porewater diffusion–reaction model, constrained by the N and O isotope systematics of nitrification and denitrification and the porewater NO 3 − isotopic composition, to estimate rates of nitrification and denitrification throughout the sediment column. Results indicate variability of reaction rates across and within the three boreholes that are generally consistent with the differential distribution of dissolved oxygen at this site, though not necessarily with the canonical view of how redox thresholds separate nitrate regeneration from dissimilative consumption spatially. That is, we provide isotope evidence for expanded zones of co-ocurring nitrification and denitrification. The isotope biogeochemical modeling also yielded estimates for the δ 15 N and δ 18 O of newly produced nitrate (δ 15 N NTR and δ 18 O NTR ), as well as the isotope effect for denitrification ( 15 ϵ DNF ), parameters with high relevance to global ocean models of N cycling. Estimated values of δ 15 N NTR were generally lower than previously reported δ 15 N values for sinking PON in this region. We suggest that these values can be related to sedimentary N-fixation and remineralization of the newly fixed organic N. Values of δ 18 O NTR generally ranged between −2.8 and 0.0 ‰, consistent with recent estimates based on lab cultures of nitrifying bacteria. Notably, some δ 18 O NTR values were elevated, suggesting incorporation of 18 O-enriched dissolved oxygen during nitrification, and possibly indicating a tight coupling of NH 4 + and NO 2 − oxidation in this metabolically sluggish environment. Our findings indicate that the production of organic matter by in situ autotrophy (e.g., nitrification, nitrogen fixation) supply a large fraction of the biomass and organic substrate for heterotrophy in these sediments, supplementing the small organic matter pool derived from the overlying euphotic zone. This work sheds new light on an active nitrogen cycle operating, despite exceedingly low carbon inputs, in the deep sedimentary biosphere.

Description: Editorial Note "A novel Whole Air Sample Profiler (WASP) for the quantification of volatile organic compounds in the boundary layer" published in Atmos. Meas. Tech., 6, 2703–2712, 2013 Atmospheric Measurement Techniques, 8, 3405-3406, 2015 Author(s): T. Wagner, H. Harder, J. Joiner, P. Laj, and A. Richter No abstract available.

Description: Spring bloom onset in the Nordic Seas Biogeosciences Discussions, 12, 13631-13673, 2015 Author(s): A. Mignot, R. Ferrari, and K. A. Mork The North Atlantic spring bloom is a massive annual growth event of marine phytoplankton, tiny free-floating algae that form the base of the ocean's food web and generates a large fraction of the global primary production of organic matter. The conditions that trigger the onset of the spring bloom in the Nordic Seas, at the northern edge of the North Atlantic, are studied using in-situ data from five bio-optical floats released above the Arctic Circle. It is often assumed that spring blooms start as soon as phytoplankton cells daily irradiance is sufficiently abundant that division rates exceed losses. The bio-optical float data instead suggest the tantalizing hypothesis that Nordic Seas blooms start when the photoperiod, the number of daily light hours experienced by phytoplankton, exceeds a critical value, independently of division rates. This bloom behavior may be explained by realizing that photosynthesis is impossible during polar nights and phytoplankton enters in a dormant stage in winter, only to be awaken by a photoperiodic trigger. While the first accumulation of biomass recorded by the bio-optical floats is consistent with the photoperiod hypothesis, it is possible that some biomass accumulation started before the critical photoperiod but at levels too low to be detected by the fluorometers. Thus more precise observations are needed to test the photoperiod hypothesis.

Description: A fuzzy intelligent system for land consolidation – a case study in Shunde, China Solid Earth, 6, 997-1006, 2015 Author(s): J. Wang, A. Ge, Y. Hu, C. Li, and L. Wang Traditionally, potential evaluation methods for farmland consolidation have depended mainly on the experts' experiences, statistical computations or subjective adjustments. Some biases usually exist in the results. Thus, computer-aided technology has become essential. In this study, an intelligent evaluation system based on a fuzzy decision tree was established, and this system can deal with numerical data, discrete data and symbolic data. When the original land data are input, the level of potential of the agricultural land for development will be output by this new model. The provision of objective proof for decision-making by authorities in rural management is helpful. Agricultural land data characteristically comprise large volumes, complex varieties and more indexes. In land consolidation, it is very important to construct an effective index system. A group of indexes need to be selected for land consolidation. In this article, a fuzzy measure was adopted to accomplish the selection of specific features. A fuzzy integral based on a fuzzy measure is a type of fusion tool. The optimal solution with the fewest non-zero elements was obtained for the fuzzy measure by solving a fuzzy integral. This algorithm provides a quick and optimal way to identify the land-index system when preparing to conduct land consolidation. This new research was applied to Shunde's "Three Old" consolidation project which provides the data. Our estimation system was compared with a conventional evaluation system that is still accepted by the public. Our results prove to be consistent, and the new model is more automatic and intelligent. The results of this estimation system are significant for informing decision-making in land consolidation.

Description: Impact of debris cover on glacier ablation and atmosphere–glacier feedbacks in the Karakoram The Cryosphere, 9, 1617-1632, 2015 Author(s): E. Collier, F. Maussion, L. I. Nicholson, T. Mölg, W. W. Immerzeel, and A. B. G. Bush The Karakoram range of the Hindu-Kush Himalaya is characterized by both extensive glaciation and a widespread prevalence of surficial debris cover on the glaciers. Surface debris exerts a strong control on glacier surface-energy and mass fluxes and, by modifying surface boundary conditions, has the potential to alter atmosphere–glacier feedbacks. To date, the influence of debris on Karakoram glaciers has only been directly assessed by a small number of glaciological measurements over short periods. Here, we include supraglacial debris in a high-resolution, interactively coupled atmosphere–glacier modeling system. To investigate glaciological and meteorological changes that arise due to the presence of debris, we perform two simulations using the coupled model from 1 May to 1 October 2004: one that treats all glacier surfaces as debris-free and one that introduces a simplified specification for the debris thickness. The basin-averaged impact of debris is a reduction in ablation of ~ 14 %, although the difference exceeds 5 m w.e. on the lowest-altitude glacier tongues. The relatively modest reduction in basin-mean mass loss results in part from non-negligible sub-debris melt rates under thicker covers and from compensating increases in melt under thinner debris, and may help to explain the lack of distinct differences in recent elevation changes between clean and debris-covered ice. The presence of debris also strongly alters the surface boundary condition and thus heat exchanges with the atmosphere; near-surface meteorological fields at lower elevations and their vertical gradients; and the atmospheric boundary layer development. These findings are relevant for glacio-hydrological studies on debris-covered glaciers and contribute towards an improved understanding of glacier behavior in the Karakoram.

Description: Introduction: The Pan-Eurasian Experiment (PEEX) – multi-disciplinary, multi-scale and multi-component research and capacity building initiative Atmospheric Chemistry and Physics Discussions, 15, 22567-22596, 2015 Author(s): M. Kulmala, H. K. Lappalainen, T. Petäjä, T. Kurten, V.-M. Kerminen, Y. Viisanen, P. Hari, V. Bondur, N. Kasimov, V. Kotlyakov, G. Matvienko, A. Baklanov, H. D. Guo, A. Ding, H.-C. Hansson, and S. Zilitinkevich The Pan-Eurasian Experiment (PEEX) is a multi-disciplinary, multi-scale and multi-component research, research infrastructure and capacity building program. PEEX has originated from a bottom-up approach by the science communities, and is aiming at resolving the major uncertainties in Earth System Science and global sustainability issues concerning the Arctic and boreal Pan-Eurasian regions, as well as China. The vision of PEEX is to solve interlinked global grand challenges influencing human well-being and societies in northern Eurasia and China. Such challenges include climate change, air quality, biodiversity loss, urbanization, chemicalization, food and fresh water availability, energy production and use of natural resources by mining, industry, energy production and transport sectors. Our approach is integrative and supra-disciplinary, recognizing the important role of the Arctic and boreal ecosystems in the Earth system. The PEEX vision includes establishing and maintaining long-term, coherent and coordinated research activities as well as continuous, comprehensive research and educational infrastructures and related capacity building across the PEEX domain. In this paper we present the PEEX structure, summarize its motivation, objectives and future outlook.

Description: Particulate organic nitrates observed in an oil and natural gas production region during wintertime Atmospheric Chemistry and Physics, 15, 9313-9325, 2015 Author(s): L. Lee, P. J. Wooldridge, J. deGouw, S. S. Brown, T. S. Bates, P. K. Quinn, and R. C. Cohen Organic nitrates in both gas and condensed (aerosol) phases were measured during the Uintah Basin Winter Ozone Study from January to February in 2012. A high degree of correlation between total aerosol volume at diameters less than 500 nm and the particulate organic nitrate concentration indicates that organic nitrates are a consistent, if not dominant, fraction of fine aerosol mass. In contrast, a similar correlation with sub-2.5 μm aerosol volume is weaker. The C : N atomic ratio inferred from field measurements of PM 2.5 and particulate organic nitrate is 34 : 1. Calculations constrained by the observations indicate that both condensation of gas-phase nitrates and heterogeneous reactions of NO 3 / N 2 O 5 are responsible for introducing organic nitrate functionality into the aerosol and that the source molecules are alkanes. Extrapolating the results to urban aerosol suggests organic nitrate production from alkanes may be a major secondary organic aerosol source.

Description: The influences of mass loading and rapid dilution of secondary organic aerosol on particle volatility Atmospheric Chemistry and Physics, 15, 9327-9343, 2015 Author(s): K. R. Kolesar, C. Chen, D. Johnson, and C. D. Cappa The thermally induced evaporation of secondary organic aerosol (SOA) has been characterized for SOA formed from the dark ozonolysis of α-pinene at initial mass concentrations ranging from 1 to 800 μg m −3 . Temperature-dependent particle size distributions were measured using a thermodenuder and the resulting mass thermograms were compared between the SOA formed at the various SOA mass concentrations. Negligible differences were observed between the mass thermograms for SOA concentrations 〈 300 μg m −3 . At higher SOA concentrations, the observed mass thermograms indicated the SOA was actually slightly less volatile than the SOA at lower concentrations; this is likely an artifact due to either saturation of the gas phase or to recondensation during cooling. The thermograms observed when the SOA was formed at high concentrations (〉 380 μg m −3 ) and then rapidly isothermally diluted to low concentrations (1–20 μg m −3 ) were identical to those for the SOA that was initially formed at low concentrations. The experimental results were compared to a kinetic model that simulates particle evaporation upon heating in a thermodenuder for a given input volatility distribution and particle composition. Three cases were considered: (1) the SOA was composed of semi-volatile monomer species with a volatility distribution based on that derived previously from consideration of SOA growth experiments; (2) the initial SOA was composed almost entirely of non-volatile dimers that decompose upon heating into their semi-volatile monomer units, which can then evaporate; and (3) where a volatility distribution was derived by fitting the model to the observed mass thermograms. It was found that good agreement is obtained between model predictions and the observations when the particle composition is dominated by either compounds of low volatility or by dimers. These same models were used to simulate isothermal evaporation of the SOA and were found to be broadly consistent with literature observations that indicate that SOA evaporation occurs with multiple timescales. The use of the semi-volatile monomer volatility distribution fails to reproduce the observed evaporation. The presence of dimers and larger oligomers in secondary organic aerosol formed from products of the reaction of α-pinene and O 3 has been well established in laboratory studies. However, the timescale and relative importance of the formation of oligomers or low-volatility compounds in the growth and evaporation of SOA has been debated. This study provides further support that low-volatility compounds and oligomers are formed in α-pinene + O 3 in high abundances and suggests that their formation occurs rapidly upon particle formation.

Description: Possible influence of atmospheric circulations on winter hazy pollution in Beijing-Tianjin-Hebei region, northern China Atmospheric Chemistry and Physics Discussions, 15, 22493-22526, 2015 Author(s): Z. Zhang, X. Zhang, D. Gong, S.-J. Kim, R. Mao, and X. Zhao Using the daily records derived from the synoptic weather stations and the NCEP/NCAR and ERA-Interim reanalysis data, the variability of the winter hazy pollutions (indicated by the mean visibility and number of hazy days) in Beijing-Tianjin-Hebei (BTH) region during the period 1981 to 2015 and its relationship to the atmospheric circulations in middle-high latitude were analyzed in this study. The winter hazy pollution in BTH had distinct inter-annual and inter-decadal variabilities without a significant long-term trend. According to the spatial distribution of correlation coefficients, six atmospheric circulation indices ( I 1 to I 6 ) were defined from the key areas in sea level pressure (SLP), zonal and meridional winds at 850 hPa (U850, V850), geopotential height field at 500 hPa (H500), zonal wind at 200 hPa (U200), and air temperature at 200 hPa (T200), respectively. All of the six indices have significant and stable correlations with the winter visibility and number of hazy days in BTH. Both the visibility and number of hazy days can be estimated well by using the six indices and fitting and the cross-validation with leave- N -out method, respectively. The high level of the prediction statistics and the reasonable mechanism suggested that the winter hazy pollutions in BTH can be forecasted or estimated credibly based on the optimized atmospheric circulation indices. However, we also noted that the statistic estimation models would be largely influenced by the artificial control of a pollutant discharge. Thus it is helpful for government decision-making departments to take actions in advance in dealing with probably severe hazy pollutions in BTH indicated by the atmospheric circulation conditions.

Description: The radiative impact of desert dust on orographic rain in the Cevennes–Vivarais area: a case study from HyMeX Atmospheric Chemistry and Physics Discussions, 15, 22451-22492, 2015 Author(s): C. Flamant, J.-P. Chaboureau, P. Chazette, P. Di Girolamo, T. Bourrianne, J. Totems, and M. Cacciani The study is focused on Intensive Observation Period (IOP) 14 of the Hydrology Cycle in the Mediterranean Experiment first Special Observing Period (HyMeX SOP 1) that took place from 17 to 19 October and was dedicated to the study of orographic rain in the Cevennes Vivarais (CV) target area. During this IOP a dense dust plume originating from North Africa (Maghreb and Sahara) was observed to be transported over the Balearic Islands towards the south of France. The plume was characterized by an aerosol optical depth between 0.2 and 0.8 at 550 nm, highly variable in time and space over the Western Mediterranean basin. The impact of this dust plume, the biggest event observed during the 2 month long HyMeX SOP 1, on the precipitation over the CV area has been analyzed using high resolution simulations from the convection permitting mesoscale model Meso-NH validated against measurements obtained from numerous instruments deployed specifically during SOP 1 (ground-based/airborne water vapor and aerosol lidars, airborne microphysics probes) as well as space-borne aerosol products. The 4 day simulation reproduced realistically the temporal and spatial variability (incl. vertical distribution) of the dust. The dust radiative impact led to an average 0.6 K heating at the altitude of the dust layer in the CV area (and up to +3 K locally) and an average 100 J kg -1 increase of most unstable convective available potential energy (and up to +900 J kg -1 locally) with respect to a simulation without prescribed dust aerosols. The rainfall amounts and location were only marginally affected by the dust radiative effect, even after 4 days of simulation. The transient nature of this radiative effect in dynamical environments such as those found in the vicinity of heavy precipitation events in the Mediterranean is not sufficient to impact 24 h accumulated rainfall in the dusty simulation.

Description: Global validation of SCIAMACHY limb ozone data (versions 2.9 and 3.0, IUP Bremen) using ozonesonde measurements Atmospheric Measurement Techniques, 8, 3369-3383, 2015 Author(s): J. Jia, A. Rozanov, A. Ladstätter-Weißenmayer, and J. P. Burrows In this paper, the latest SCIAMACHY limb ozone scientific vertical profiles, namely the current V2.9 and the upcoming V3.0, are extensively compared with ozonesonde data from the World Ozone and Ultraviolet Radiation Data Centre (WOUDC) database. The comparisons are made on a global scale from 2003 to 2011, involving 61 sonde stations. The retrieval processors used to generate V2.9 and V3.0 data sets are briefly introduced. The comparisons are discussed in terms of vertical profiles and stratospheric partial columns. Our results indicate that the V2.9 ozone profile data between 20 and 30 km are in good agreement with ground-based measurements, with less than 5 % relative differences in the latitude range of 90° S–40° N (with the exception of the tropical Pacific region, where an overestimation of more than 10 % is observed), which corresponds to less than 5 DU partial-column differences. In the tropics the differences are within 3 %. However, this data set shows a significant underestimation northwards of 40° N (up to ~ 15 %). The newly developed V3.0 data set reduces this bias to below 10 % while maintaining a good agreement southwards of 40° N with slightly increased relative differences of up to 5 % in the tropics.

Description: Unusual biogenic calcite structures in two shallow lakes, James Ross Island, Antarctica Biogeosciences Discussions, 12, 13593-13629, 2015 Author(s): J. Elster, L. Nedbalová, R. Vodrážka, K. Láska, J. Haloda, and J. Komárek The floors of two shallow endorheic lakes, located on volcanic surfaces on James Ross Island, are covered with calcareous organosedimentary structures. Their biological and chemical composition, lake water characteristics, and seasonal variability of the thermal regime are introduced. The lakes are frozen down to the bottom eight-nine months per year and their water chemistry is characterized by low conductivity and neutral to slightly alkaline pH. The photosynthetic microbial mat is composed of filamentous cyanobacteria and microalgae that are considered to be Antarctic endemic species. The mucilaginous black biofilm is covered by green spots formed by a green microalga and the macroscopic structures are packed together with fine material. Thin sections consist of rock substrate, soft biofilm, calcite spicules and mineral grains originating from different sources. The morphology of the spicules is typical of calcium carbonate monocrystals having a layered structure and worn surface, which reflect growth and degradation processes. The spicules chemical composition and structure correspond to pure calcite. Lakes age, altitude, morphometry, geomorphological and hydrological stability, including low sedimentation rates, together with thermal regime predispose the existence of this community. We hypothesize that the precipitation of calcite is connected with the photosynthetic activity of the green microalgae that were not recorded in any other lake in the region. This study has shown that the unique community producing biogenic calcite spicules is quite different to any yet described.

Description: The GEWEX LandFlux project: evaluation of model evaporation using tower-based and globally-gridded forcing data Geoscientific Model Development Discussions, 8, 6809-6866, 2015 Author(s): M. F. McCabe, A. Ershadi, C. Jimenez, D. G. Miralles, D. Michel, and E. F. Wood Determining the spatial distribution and temporal development of evaporation at regional and global scales is required to improve our understanding of the coupled water and energy cycles and to better monitor any changes in observed trends and variability of linked hydrological processes. With recent international efforts guiding the development of long-term and globally distributed flux estimates, continued product assessments are required to inform upon the selection of suitable model structures and also to establish the appropriateness of these multi-model simulations for global application. In support of the objectives of the GEWEX LandFlux project, four commonly used evaporation models are evaluated against data from tower-based eddy-covariance observations, distributed across a range of biomes and climate zones. The selected schemes include the Surface Energy Balance System (SEBS) approach, the Priestley-Taylor Jet Propulsion Laboratory (PT-JPL) model, the Penman-Monteith based Mu model (PM-Mu) and the Global Land Evaporation: the Amsterdam Methodology (GLEAM). Here we seek to examine the fidelity of global evaporation simulations by examining the multi-model response to varying sources of forcing data. To do this, we perform parallel and collocated model simulations using tower-based data together with a global-scale grid-based forcing product. Through quantifying the multi-model response to high-quality tower data, a better understanding of the subsequent model response to coarse-scale globally gridded data that underlies the LandFlux product can be obtained, while also providing a relative evaluation and assessment of model performance. Using surface flux observations from forty-five globally distributed eddy-covariance stations as independent metrics of performance, the tower-based analysis indicated that PT-JPL provided the highest overally statistical performance (0.72; 61 W m −2 ; 0.65), followed closely by GLEAM (0.68; 64 W m −2 ; 0.62), with values in parenthesis representing the R 2 , RMSD and Nash-Sutcliffe Efficiency (NSE) and respectively. PM-Mu (0.51; 78 W m −2 ; 0.45) tended to underestimate fluxes, while SEBS (0.72; 101 W m −2 ; 0.24) overestimated values relative to observations. A focused analysis across specific biome types and climate zones showed considerable variability in the performance of all models, with no single model consistently able to outperform any other. Results also indicated that the global gridded data tended to reduce the performance for all of the studied models when compared to the tower data, likely a response to scale mismatch and issues related to forcing quality. Rather than relying on any single model simulation, the spatial and temporal variability at both the tower- and grid-scale highlighted the potential benefits of developing an ensemble or blended evaporation product for global scale LandFlux applications. Challenges related to the robust assessment of the LandFlux product are also discussed.

Description: Drought assessment in the Dongliao River basin: traditional approaches vs. generalized drought assessment index based on water resources systems Natural Hazards and Earth System Science, 15, 1889-1906, 2015 Author(s): B. S. Weng, D. H. Yan, H. Wang, J. H. Liu, Z. Y. Yang, T. L. Qin, and J. Yin Drought is firstly a resource issue, and with its development it evolves into a disaster issue. Drought events usually occur in a determinate but a random manner. Drought has become one of the major factors to affect sustainable socioeconomic development. In this paper, we propose the generalized drought assessment index (GDAI) based on water resources systems for assessing drought events. The GDAI considers water supply and water demand using a distributed hydrological model. We demonstrate the use of the proposed index in the Dongliao River basin in northeastern China. The results simulated by the GDAI are compared to observed drought disaster records in the Dongliao River basin. In addition, the temporal distribution of drought events and the spatial distribution of drought frequency from the GDAI are compared with the traditional approaches in general (i.e., standard precipitation index, Palmer drought severity index and rate of water deficit index). Then, generalized drought times, generalized drought duration, and generalized drought severity were calculated by theory of runs. Application of said runs at various drought levels (i.e., mild drought, moderate drought, severe drought, and extreme drought) during the period 1960–2010 shows that the centers of gravity of them all distribute in the middle reaches of Dongliao River basin, and change with time. The proposed methodology may help water managers in water-stressed regions to quantify the impact of drought, and consequently, to make decisions for coping with drought.

Description: Current model capabilities for simulating black carbon and sulfate concentrations in the Arctic atmosphere: a multi-model evaluation using a comprehensive measurement data set Atmospheric Chemistry and Physics, 15, 9413-9433, 2015 Author(s): S. Eckhardt, B. Quennehen, D. J. L. Olivié, T. K. Berntsen, R. Cherian, J. H. Christensen, W. Collins, S. Crepinsek, N. Daskalakis, M. Flanner, A. Herber, C. Heyes, Ø. Hodnebrog, L. Huang, M. Kanakidou, Z. Klimont, J. Langner, K. S. Law, M. T. Lund, R. Mahmood, A. Massling, S. Myriokefalitakis, I. E. Nielsen, J. K. Nøjgaard, J. Quaas, P. K. Quinn, J.-C. Raut, S. T. Rumbold, M. Schulz, S. Sharma, R. B. Skeie, H. Skov, T. Uttal, K. von Salzen, and A. Stohl The concentrations of sulfate, black carbon (BC) and other aerosols in the Arctic are characterized by high values in late winter and spring (so-called Arctic Haze) and low values in summer. Models have long been struggling to capture this seasonality and especially the high concentrations associated with Arctic Haze. In this study, we evaluate sulfate and BC concentrations from eleven different models driven with the same emission inventory against a comprehensive pan-Arctic measurement data set over a time period of 2 years (2008–2009). The set of models consisted of one Lagrangian particle dispersion model, four chemistry transport models (CTMs), one atmospheric chemistry-weather forecast model and five chemistry climate models (CCMs), of which two were nudged to meteorological analyses and three were running freely. The measurement data set consisted of surface measurements of equivalent BC (eBC) from five stations (Alert, Barrow, Pallas, Tiksi and Zeppelin), elemental carbon (EC) from Station Nord and Alert and aircraft measurements of refractory BC (rBC) from six different campaigns. We find that the models generally captured the measured eBC or rBC and sulfate concentrations quite well, compared to previous comparisons. However, the aerosol seasonality at the surface is still too weak in most models. Concentrations of eBC and sulfate averaged over three surface sites are underestimated in winter/spring in all but one model (model means for January–March underestimated by 59 and 37 % for BC and sulfate, respectively), whereas concentrations in summer are overestimated in the model mean (by 88 and 44 % for July–September), but with overestimates as well as underestimates present in individual models. The most pronounced eBC underestimates, not included in the above multi-site average, are found for the station Tiksi in Siberia where the measured annual mean eBC concentration is 3 times higher than the average annual mean for all other stations. This suggests an underestimate of BC sources in Russia in the emission inventory used. Based on the campaign data, biomass burning was identified as another cause of the modeling problems. For sulfate, very large differences were found in the model ensemble, with an apparent anti-correlation between modeled surface concentrations and total atmospheric columns. There is a strong correlation between observed sulfate and eBC concentrations with consistent sulfate/eBC slopes found for all Arctic stations, indicating that the sources contributing to sulfate and BC are similar throughout the Arctic and that the aerosols are internally mixed and undergo similar removal. However, only three models reproduced this finding, whereas sulfate and BC are weakly correlated in the other models. Overall, no class of models (e.g., CTMs, CCMs) performed better than the others and differences are independent of model resolution.

Description: Investigation of post-depositional processing of nitrate in East Antarctic snow: isotopic constraints on photolytic loss, re-oxidation, and source inputs Atmospheric Chemistry and Physics, 15, 9435-9453, 2015 Author(s): G. Shi, A. M. Buffen, M. G. Hastings, C. Li, H. Ma, Y. Li, B. Sun, C. An, and S. Jiang Snowpits along a traverse from coastal East Antarctica to the summit of the ice sheet (Dome Argus) are used to investigate the post-depositional processing of nitrate (NO 3 − ) in snow. Seven snowpits from sites with accumulation rates between 24 and 172 kg m −2 a −1 were sampled to depths of 150 to 300 cm. At sites from the continental interior (low accumulation, 〈 55 kg m −2 a −1 ), nitrate mass fraction is generally 〉 200 ng g −1 in surface snow and decreases quickly with depth to 〈 50 ng g −1 . Considerably increasing values of δ 15 N of nitrate are also observed (16–461 ‰ vs. air N 2 ), particularly in the top 20 cm, which is consistent with predicted fractionation constants for the photolysis of nitrate. The δ 18 O of nitrate (17–84 ‰ vs. VSMOW (Vienna Standard Mean Ocean Water)), on the other hand, decreases with increasing δ 15 N, suggestive of secondary formation of nitrate in situ (following photolysis) with a low δ 18 O source. Previous studies have suggested that δ 15 N and δ 18 O of nitrate at deeper snow depths should be predictable based upon an exponential change derived near the surface. At deeper depths sampled in this study, however, the relationship between nitrate mass fraction and δ 18 O changes, with increasing δ 18 O of nitrate observed between 100 and 200 cm. Predicting the impact of post-depositional loss, and therefore changes in the isotopes with depth, is highly sensitive to the depth interval over which an exponential change is assumed. In the snowpits collected closer to the coast (accumulation 〉 91 kg m −2 a −1 ), there are no obvious trends detected with depth and instead seasonality in nitrate mass fraction and isotopic composition is found. In comparison to the interior sites, the coastal pits are lower in δ 15 N (−15–71 ‰ vs. air N 2 ) and higher in δ 18 O of nitrate (53–111 ‰ vs. VSMOW). The relationships found amongst mass fraction, δ 15 N, δ 18 O and Δ 17 O (Δ 17 O = δ 17 O–0.52 × δ 18 O) of nitrate cannot be explained by local post-depositional processes alone, and are instead interpreted in the context of a primary atmospheric signal. Consistent with other Antarctic observational and modeling studies, the isotopic results are suggestive of an important influence of stratospheric ozone chemistry on nitrate formation during the cold season and a mix of tropospheric sources and chemistry during the warm season. Overall, the findings in this study speak to the sensitivity of nitrate isotopic composition to post-depositional processing and highlight the strength of combined use of the nitrogen and oxygen isotopes for a mechanistic understanding of this processing.

Description: NO x emission estimates during the 2014 Youth Olympic Games in Nanjing Atmospheric Chemistry and Physics, 15, 9399-9412, 2015 Author(s): J. Ding, R. J. van der A, B. Mijling, P. F. Levelt, and N. Hao The Nanjing Government applied temporary environmental regulations to guarantee good air quality during the Youth Olympic Games (YOG) in 2014. We study the effect of those regulations by applying the emission estimate algorithm DECSO (Daily Emission estimates Constrained by Satellite Observations) to measurements of the Ozone Monitoring Instrument (OMI). We improved DECSO by updating the chemical transport model CHIMERE from v2006 to v2013 and by adding an Observation minus Forecast (OmF) criterion to filter outlying satellite retrievals due to high aerosol concentrations. The comparison of model results with both ground and satellite observations indicates that CHIMERE v2013 is better performing than CHIMERE v2006. After filtering the satellite observations with high aerosol loads that were leading to large OmF values, unrealistic jumps in the emission estimates are removed. Despite the cloudy conditions during the YOG we could still see a decrease of tropospheric NO 2 column concentrations of about 32 % in the OMI observations when compared to the average NO 2 columns from 2005 to 2012. The results of the improved DECSO algorithm for NO x emissions show a reduction of at least 25 % during the YOG period and afterwards. This indicates that air quality regulations taken by the local government have an effect in reducing NO x emissions. The algorithm is also able to detect an emission reduction of 10 % during the Chinese Spring Festival. This study demonstrates the capacity of the DECSO algorithm to capture the change of NO x emissions on a monthly scale. We also show that the observed NO 2 columns and the derived emissions show different patterns that provide complimentary information. For example, the Nanjing smog episode in December 2013 led to a strong increase in NO 2 concentrations without an increase in NO x emissions. Furthermore, DECSO gives us important information on the non-trivial seasonal relation between NO x emissions and NO 2 concentrations on a local scale.

Description: Sensitivity of remotely-sensed trace gas concentrations to polarisation Atmospheric Measurement Techniques Discussions, 8, 8779-8816, 2015 Author(s): D. M. O'Brien, I. N. Polonsky, and J. B. Kumer Current and proposed space missions estimate column-averaged concentrations of trace gases (CO 2 , CH 4 and CO) from high resolution spectra of reflected sunlight in absorption bands of the gases. The radiance leaving the top of the atmosphere is partially polarised by both reflection at the surface and scattering within the atmosphere. Generally the polarisation state is unknown, and could degrade the accuracy of the concentration measurements. The sensitivity to polarisation is modelled for the proposed geoCARB instrument, which will include neither polarisers nor polarisation scramblers to select particular polarisation states from the incident radiation. The radiometric and polarimetric calibrations proposed for geoCARB are outlined, and a model is developed for the polarisation properties of the geoCARB spectrographs. This model depends principally upon the efficiencies of the gratings to polarisations parallel and perpendicular to the rulings of the gratings. Next an ensemble of polarised spectra is simulated for geoCARB observing targets in India, China and Australia from geostationary orbit at longitude 110° E. The spectra are analysed to recover the trace gas concentrations in two modes, the first denied access to the polarimetric calibration and the second with access. The retrieved concentrations using the calibration data are almost identical to those that would be obtained with polarisation scramblers, while the retrievals without calibration data contain outliers that do not meet the accuracies demanded by the mission.

Description: Observations of atmospheric mercury in China: a critical review Atmospheric Chemistry and Physics, 15, 9455-9476, 2015 Author(s): X. W. Fu, H. Zhang, B. Yu, X. Wang, C.-J. Lin, and X. B. Feng China presently contributes the largest amount of anthropogenic mercury (Hg) emission into the atmosphere in the world. Over the past decade, numerous studies have been conducted to characterize the concentration and forms of atmospheric Hg in China, which provide insights into the spatial and temporal distributions of atmospheric Hg through ground-based measurements at widely diverse geographical locations and during cruise and flight campaigns. In this paper, we present a comprehensive review of the state of understanding in atmospheric Hg in China. Gaseous elemental mercury (GEM) and particulate-bound mercury (PBM) measured at the remote sites in China are substantially elevated compared to the background values in the Northern Hemisphere. In Chinese urban areas, the highly elevated GEM, PBM and gaseous oxidized mercury (GOM) were mainly derived from local anthropogenic Hg emissions, whereas regional anthropogenic emissions and long-range transport from domestic source regions are the primary causes of the elevated GEM and PBM concentrations at remote sites. Using 7–9 years of continuous observations at a remote site and an urban site, a slight increase in atmospheric GEM (2.4–2.5 % yr −1 ) was identified (paired samples test: p 〈 0.01), which is in agreement with the increasing domestic anthropogenic emissions. Anthropogenic GEM emission quantity in China estimated through the observed GEM / CO concentration ratios ranged from 632 to 1138 t annually over the past decade, 2–3 times larger than published values using emission activity data. Modeling results and filed measurements show dry deposition is the predominant process for removing Hg from the atmosphere, 2.5–9.0 times larger than wet deposition, due to the elevated atmospheric GEM and PBM concentrations that facilitate dry deposition to terrestrial landscapes. Further studies to reconcile the observed and simulated Hg concentrations, to understand the impact of domestic emission reduction on Hg concentration and deposition and to delineate the role of Hg emission and deposition of China in the global Hg biogeochemical cycle, are needed.

Description: Exploiting the sensitivity of two satellite cloud height retrievals to cloud vertical distribution Atmospheric Measurement Techniques, 8, 3419-3431, 2015 Author(s): C. K. Carbajal Henken, L. Doppler, R. Lindstrot, R. Preusker, and J. Fischer This work presents a study on the sensitivity of two satellite cloud height retrievals to cloud vertical distribution. The difference in sensitivity is exploited by relating the difference in the retrieved cloud heights to cloud vertical extent. The two cloud height retrievals, performed within the Freie Universität Berlin AATSR MERIS Cloud (FAME-C) algorithm, are based on independent measurements and different retrieval techniques. First, cloud-top temperature (CTT) is retrieved from Advanced Along Track Scanning Radiometer (AATSR) measurements in the thermal infrared. Second, cloud-top pressure (CTP) is retrieved from Medium Resolution Imaging Spectrometer (MERIS) measurements in the oxygen-A absorption band and a nearby window channel. Both CTT and CTP are converted to cloud-top height (CTH) using atmospheric profiles from a numerical weather prediction model. First, a sensitivity study using radiative transfer simulations in the near-infrared and thermal infrared was performed to demonstrate, in a quantitative manner, the larger impact of the assumed cloud vertical extinction profile, described in terms of shape and vertical extent, on MERIS than on AATSR top-of-atmosphere measurements. Consequently, cloud vertical extinction profiles will have a larger influence on the MERIS than on the AATSR cloud height retrievals for most cloud types. Second, the difference in retrieved CTH (ΔCTH) from AATSR and MERIS are related to cloud vertical extent (CVE), as observed by ground-based lidar and radar at three ARM sites. To increase the impact of the cloud vertical extinction profile on the MERIS-CTP retrievals, single-layer and geometrically thin clouds are assumed in the forward model. Similarly to previous findings, the MERIS-CTP retrievals appear to be close to pressure levels in the middle of the cloud. Assuming a linear relationship, the ΔCTH multiplied by 2.5 gives an estimate on the CVE for single-layer clouds. The relationship is stronger for single-layer clouds than for multi-layer clouds. Due to large variations of cloud vertical extinction profiles occurring in nature, a quantitative estimate of the cloud vertical extent is accompanied with large uncertainties. Yet, estimates of the CVE provide an additional parameter, next to CTH, that can be obtained from passive imager measurements and can be used to further describe cloud vertical distribution, thus contributing to the characterization of a cloudy scene. To further demonstrate the plausibility of the approach, an estimate of the CVE was applied to a case study. In light of the follow-up mission Sentinel-3 with AATSR and MERIS like instruments, Sea and Land Surface Temperature Radiometer (SLSTR) and (Ocean and Land Colour Instrument) OLCI, respectively, for which the FAME-C algorithm can be easily adapted, a more accurate estimate of the CVE can be expected. OLCI will have three channels in the oxygen-A absorption band, possibly providing enhanced information on cloud vertical distributions.

Description: The electrical self-potential method is a non-intrusive snow-hydrological sensor The Cryosphere Discussions, 9, 4437-4457, 2015 Author(s): S. S. Thompson, B. Kulessa, R. L. H. Essery, and M. P. Lüthi Our ability to measure, quantify and assimilate hydrological properties and processes of snow in operational models is disproportionally poor compared to the significance of seasonal snowmelt as a global water resource and major risk factor in flood and avalanche forecasting. Encouraged by recent theoretical, modelling and laboratory work, we show here that the diurnal evolution of aerially-distributed self-potential magnitudes closely track those of bulk meltwater fluxes in melting in-situ snowpacks at Rhone and Jungfraujoch glaciers, Switzerland. Numerical modelling infers temporally-evolving liquid water contents in the snowpacks on successive days in close agreement with snow-pit measurements. Muting previous concerns, the governing physical and chemical properties of snow and meltwater became temporally invariant for modelling purposes. Because measurement procedure is straightforward and readily automated for continuous monitoring over significant spatial scales, we conclude that the self-potential geophysical method is a highly-promising non-intrusive snow-hydrological sensor for measurement practice, modelling and operational snow forecasting.

Description: Association of the pre-monsoon thermal field over north India and the western Tibetan Plateau with summer monsoon rainfall over India Annales Geophysicae, 33, 1051-1058, 2015 Author(s): S. D. Bansod, S. Fadnavis, and S. P. Ghanekar In this paper, interannual variability of tropospheric air temperatures over the Asian summer monsoon region during the pre-monsoon months is examined in relation to Indian summer monsoon rainfall (ISMR; June to September total rainfall). For this purpose, monthly grid-point temperatures in the entire troposphere over the Asian summer monsoon region and ISMR data for the period 1949–2012 have been used. Spatial correlation patterns are investigated between the temperature field in the lower tropospheric levels during May over the Asian summer monsoon region and ISMR. The results indicate a strong and significant northwest–southeast dipole structure in the spatial correlations over the Indian region, with highly significant positive (negative) correlations over the regions of north India and the western Tibetan Plateau region – region R1 (north Bay of Bengal: region R2). The observed dipole is seen significantly up to a level of 850 hPa and eventually disappears at 700 hPa. Thermal indices evaluated at 850 hPa level, based on average air temperatures over the north India and western Tibetan Plateau region (TI 1 ) and the north Bay of Bengal region (TI 2 ) during May, show a strong, significant relationship with the ISMR. The results are found to be consistent and robust, especially in the case of TI 1 during the period of analysis. A physical mechanism for the relationship between these indices and ISMR is proposed. Finally the composite annual cycle of tropospheric air temperature over R1 during flood/drought years of ISMR is examined. The study brings out the importance of the TI 1 in the prediction of flood/drought conditions over the Indian subcontinent.

Description: A wing pod-based millimeter wavelength airborne cloud radar Geoscientific Instrumentation, Methods and Data Systems, 4, 161-176, 2015 Author(s): J. Vivekanandan, S. Ellis, P. Tsai, E. Loew, W.-C. Lee, J. Emmett, M. Dixon, C. Burghart, and S. Rauenbuehler This paper describes a novel, airborne pod-based millimeter (mm) wavelength radar. Its frequency of operation is 94 GHz (3 mm wavelength). The radar has been designed to fly on the NCAR Gulfstream V HIAPER aircraft; however, it could be deployed on other similarly equipped aircraft. The pod-based configuration occupies minimum cabin space and maximizes scan coverage. The radar system is capable of collecting observations in a staring mode between zenith and nadir or in a scanning mode. Standard pulse-pair estimates of moments and raw time series of backscattered signals are recorded. The radar system design and characteristics as well as techniques for calibrating reflectivity and correcting Doppler velocity for aircraft attitude and motion are described. The radar can alternatively be deployed in a ground-based configuration, housed in the 20 ft shipping container it shares with the High Spectral Resolution Lidar (HSRL). The radar was tested both on the ground and in flight. Preliminary measurements of Doppler and polarization measurements were collected and examples are presented.

Description: A rainfall design method for spatial flood risk assessment: considering multiple flood sources Hydrology and Earth System Sciences Discussions, 12, 8005-8033, 2015 Author(s): X. Jiang and H. Tatano Information about the spatial distribution of flood risk is important for integrated urban flood risk management. Focusing on urban areas, spatial flood risk assessment must reflect all risk information derived from multiple flood sources: rivers, drainage, coastal flooding etc. that may affect the area. However, conventional flood risk assessment deals with each flood source independently, which leads to an underestimation of flood risk in the floodplain. Even in floodplains that have no risk from coastal flooding, flooding from river channels and inundation caused by insufficient drainage capacity should be considered simultaneously. For integrated flood risk management, it is necessary to establish a methodology to estimate flood risk distribution across a floodplain. In this paper, a rainfall design method for spatial flood risk assessment, which considers the joint effects of multiple flood sources, is proposed. The concept of critical rainfall duration determined by the concentration time of flooding is introduced to connect response characteristics of different flood sources with rainfall. A copula method is then adopted to capture the correlation of rainfall amount with different critical rainfall durations. Rainfall events are designed taking advantage of the copula structure of correlation and marginal distribution of rainfall amounts within different critical rainfall durations. A case study in the Otsu River Basin, Osaka prefecture, Japan was conducted to demonstrate this methodology.

Description: Tilt error in cryospheric surface radiation measurements at high latitudes: a model study The Cryosphere Discussions, 9, 4355-4376, 2015 Author(s): W. S. Bogren, J. F. Burkhart, and A. Kylling We have evaluated the magnitude and makeup of error in cryospheric radiation observations due to small sensor misalignment in in-situ measurements of solar irradiance. This error is examined through simulation of diffuse and direct irradiance arriving at a detector with a cosine-response foreoptic. Emphasis is placed on assessing total error over the solar shortwave spectrum from 250 to 4500 nm, as well as supporting investigation over other relevant shortwave spectral ranges. The total measurement error introduced by sensor tilt is dominated by the direct component. For a typical high latitude albedo measurement with a solar zenith angle of 60°, a sensor tilted by 1, 3, and 5° can respectively introduce up to 2.6, 7.7, and 12.8 % error into the measured irradiance and similar errors in the derived albedo. Depending on the daily range of solar azimuth and zenith angles, significant measurement error can persist also in integrated daily irradiance and albedo.

Description: A better understanding of hydroxyl radical photochemical sources in cloud waters collected at the puy de Dôme station – experimental versus modelled formation rates Atmospheric Chemistry and Physics, 15, 9191-9202, 2015 Author(s): A. Bianco, M. Passananti, H. Perroux, G. Voyard, C. Mouchel-Vallon, N. Chaumerliac, G. Mailhot, L. Deguillaume, and M. Brigante The oxidative capacity of the cloud aqueous phase is investigated during three field campaigns from 2013 to 2014 at the top of the puy de Dôme station (PUY) in France. A total of 41 cloud samples are collected and the corresponding air masses are classified as highly marine, marine and continental. Hydroxyl radical (HO•) formation rates (R HO• f ) are determined using a photochemical setup (xenon lamp that can reproduce the solar spectrum) and a chemical probe coupled with spectroscopic analysis that can trap all of the generated radicals for each sample. Using this method, the obtained values correspond to the total formation of HO• without its chemical sinks. These formation rates are correlated with the concentrations of the naturally occurring sources of HO•, including hydrogen peroxide, nitrite, nitrate and iron. The total hydroxyl radical formation rates are measured as ranging from approximately 2 × 10 −11 to 4 × 10 −10 M s −1 , and the hydroxyl radical quantum yield formation (Φ HO• ) is estimated between 10 −4 and 10 −2 . Experimental values are compared with modelled formation rates calculated by the model of multiphase cloud chemistry (M2C2), considering only the chemical sources of the hydroxyl radicals. The comparison between the experimental and the modelled results suggests that the photoreactivity of the iron species as a source of HO• is overestimated by the model, and H 2 O 2 photolysis represents the most important source of this radical (between 70 and 99 %) for the cloud water sampled at the PUY station (primarily marine and continental).

Description: Dust–air pollution dynamics over the eastern Mediterranean Atmospheric Chemistry and Physics, 15, 9173-9189, 2015 Author(s): M. Abdelkader, S. Metzger, R. E. Mamouri, M. Astitha, L. Barrie, Z. Levin, and J. Lelieveld Interactions of desert dust and air pollution over the eastern Mediterranean (EM) have been studied, focusing on two distinct dust transport events on 22 and 28 September 2011. The atmospheric chemistry–climate model EMAC has been used at about 50 km grid spacing, applying an online dust emission scheme and calcium as a proxy for dust reactivity. EMAC includes a detailed tropospheric chemistry mechanism, aerosol microphysics and thermodynamics schemes to describe dust "aging". The model is evaluated using ground-based observations for aerosol concentrations and aerosol optical depth (AOD) as well as satellite observations. Simulation results and back trajectory analysis show that the development of synoptic disturbances over the EM can enhance dust transport from the Sahara and Arabian deserts in frontal systems that also carry air pollution to the EM. The frontal systems are associated with precipitation that controls the dust removal. Our results show the importance of chemical aging of dust, which increases particle size, dust deposition and scavenging efficiency during transport, overall reducing the lifetime relative to non-aged dust particles. The relatively long travel periods of Saharan dust result in more sustained aging compared to Arabian dust. Sensitivity simulations indicate 3 times more dust deposition of aged relative to pristine dust, which significantly decreases the dust lifetime and loading.

Description: Accuracy of retrieving temperature and humidity profiles by ground-based microwave radiometry in truly complex terrain Atmospheric Measurement Techniques, 8, 3355-3367, 2015 Author(s): G. Massaro, I. Stiperski, B. Pospichal, and M. W. Rotach Within the Innsbruck Box project, a ground-based microwave radiometer (RPG-HATPRO) was operated in the Inn Valley (Austria), in very complex terrain, between September 2012 and May 2013 to obtain temperature and humidity vertical profiles of the full troposphere with a specific focus on the valley boundary layer. In order to assess its performance in a deep alpine valley, the profiles obtained by the radiometer with different retrieval algorithms based on different climatologies are compared to local radiosonde data. A retrieval that is improved with respect to the one provided by the manufacturer, based on better resolved data, shows a significantly smaller root mean square error (RMSE), both for the temperature and humidity profiles. The improvement is particularly substantial at the heights close to the mountaintop level and in the upper troposphere. Lower-level inversions, common in an alpine valley, are resolved to a satisfactory degree. On the other hand, upper-level inversions (above 1200 m) still pose a significant challenge for retrieval. For this purpose, specialized retrieval algorithms were developed by classifying the radiosonde climatologies into specialized categories according to different criteria (seasons, daytime, nighttime) and using additional regressors (e.g., measurements from mountain stations). The training and testing on the radiosonde data for these specialized categories suggests that a classification of profiles that reproduces meaningful physical characteristics can yield improved targeted specialized retrievals. A novel and very promising method of improving the profile retrieval in a mountainous region is adding further information in the retrieval, such as the surface temperature at fixed levels along a topographic slope or from nearby mountaintops.

Description: An inkling of the relation between the monofractality of temperatures and pressure anomalies Nonlinear Processes in Geophysics Discussions, 2, 1339-1353, 2015 Author(s): A. Deliège and S. Nicolay We use the discrete "wavelet transform microscope" to study the monofractal nature of surface air temperature signals of weather stations spread across Europe. This method reveals that the information obtained in this way is richer than previous works studying long range correlations in meteorological stations: the approach presented here allows to bind the Hölder exponents with the standard deviation of surface pressure anomalies, while such a link does not appear with methods previously carried out.

Description: Explanation of the values of Hack's drainage basin, river length scaling exponent Nonlinear Processes in Geophysics Discussions, 2, 1355-1361, 2015 Author(s): A. G. Hunt Percolation theory can be used to find water flow paths of least resistance. The application of percolation theory to drainage networks allows identification of the range of exponent values that describe the tortuosity of rivers in real river networks, which is then used to generate the observed scaling between drainage basin area and channel length, a relationship known as Hack's law. Such a theoretical basis for Hack's law allows interpretation of the range of exponent values based on an assessment of the heterogeneity of the substrate.

Description: Spectral diagonal ensemble Kalman filters Nonlinear Processes in Geophysics, 22, 485-497, 2015 Author(s): I. Kasanický, J. Mandel, and M. Vejmelka A new type of ensemble Kalman filter is developed, which is based on replacing the sample covariance in the analysis step by its diagonal in a spectral basis. It is proved that this technique improves the approximation of the covariance when the covariance itself is diagonal in the spectral basis, as is the case, e.g., for a second-order stationary random field and the Fourier basis. The method is extended by wavelets to the case when the state variables are random fields which are not spatially homogeneous. Efficient implementations by the fast Fourier transform (FFT) and discrete wavelet transform (DWT) are presented for several types of observations, including high-dimensional data given on a part of the domain, such as radar and satellite images. Computational experiments confirm that the method performs well on the Lorenz 96 problem and the shallow water equations with very small ensembles and over multiple analysis cycles.

Description: Design and validation of MEDRYS, a Mediterranean Sea reanalysis over 1992–2013 Ocean Science Discussions, 12, 1815-1867, 2015 Author(s): M. Hamon, J. Beuvier, S. Somot, J. M. Lellouche, E. Greiner, G. Jordà, M. N. Bouin, T. Arsouze, K. Béranger, F. Sevault, C. Dubois, M. Drevillon, and Y. Drillet The French research community on the Mediterranean Sea modelling and the French operational ocean forecasting center Mercator Océan have gathered their skill and expertise in physical oceanography, ocean modelling, atmospheric forcings and data assimilation, to carry out a MEDiterranean sea ReanalYsiS (MEDRYS) at high resolution for the period 1992–2013. The ocean model used is NEMOMED12, a Mediterranean configuration of NEMO with a 1/12° (∼ 7 km) horizontal resolution and 75 vertical z levels with partial steps. At the surface, it is forced by a new atmospheric forcing dataset (ALDERA), coming from a dynamical downscaling of the ERA-Interim atmospheric reanalysis by the regional climate model ALADIN-Climate with a 12 km horizontal and 3 h temporal resolutions. This configuration is used to carry a 34 year free simulation over the period 1979–2013 (NM12-FREE) which is the initial state of the reanalysis in October 1992. The first version of MEDRYS uses the existing Mercator Océan data assimilation system SAM that is based on a reduced-order Kalman filter with a 3-D multivariate modal decomposition of the forecast error. Altimeter data, satellite SST and temperature and salinity vertical profiles are jointly assimilated. This paper describes the configuration we used to perform the MEDRYS simulation. We then first validate the skills of the data assimilation system. It is shown that the data assimilation restores a good averaged temperature and salinity in intermediate layers compared to the free simulation. No particular biases are identified in the bottom layers. However, the reanalysis show slight positive biases of 0.02 psu and 0.15 °C above 150 m depth. In the validation stage, it is also shown that the assimilation allows to better reproduce water, heat and salt transports through the Strait of Gibraltar. Finally, the ability of the reanalysis to represent the sea surface high frequency variability is pointed out.

Description: Hydroelastic analysis of ice shelves under long wave excitation Natural Hazards and Earth System Science, 15, 1851-1857, 2015 Author(s): T. K. Papathanasiou, A. E. Karperaki, E. E. Theotokoglou, and K. A. Belibassakis The transient hydroelastic response of an ice shelf under long wave excitation is analysed by means of the finite element method. The simple model, presented in this work, is used for the simulation of the generated kinematic and stress fields in an ice shelf, when the latter interacts with a tsunami wave. The ice shelf, being of large length compared to its thickness, is modelled as an elastic Euler-Bernoulli beam, constrained at the grounding line. The hydrodynamic field is represented by the linearised shallow water equations. The numerical solution is based on the development of a special hydroelastic finite element for the system of governing of equations. Motivated by the 2011 Sulzberger Ice Shelf (SIS) calving event and its correlation with the Honshu Tsunami, the SIS stable configuration is studied. The extreme values of the bending moment distribution in both space and time are examined. Finally, the location of these extrema is investigated for different values of ice shelf thickness and tsunami wave length.

Description: A method for merging nadir-sounding climate records, with an application to the global-mean stratospheric temperature data sets from SSU and AMSU Atmospheric Chemistry and Physics, 15, 9271-9284, 2015 Author(s): C. McLandress, T. G. Shepherd, A. I. Jonsson, T. von Clarmann, and B. Funke A method is proposed for merging different nadir-sounding climate data records using measurements from high-resolution limb sounders to provide a transfer function between the different nadir measurements. The two nadir-sounding records need not be overlapping so long as the limb-sounding record bridges between them. The method is applied to global-mean stratospheric temperatures from the NOAA Climate Data Records based on the Stratospheric Sounding Unit (SSU) and the Advanced Microwave Sounding Unit-A (AMSU), extending the SSU record forward in time to yield a continuous data set from 1979 to present, and providing a simple framework for extending the SSU record into the future using AMSU. SSU and AMSU are bridged using temperature measurements from the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), which is of high enough vertical resolution to accurately represent the weighting functions of both SSU and AMSU. For this application, a purely statistical approach is not viable since the different nadir channels are not sufficiently linearly independent, statistically speaking. The near-global-mean linear temperature trends for extended SSU for 1980–2012 are −0.63 ± 0.13, −0.71 ± 0.15 and −0.80 ± 0.17 K decade −1 (95 % confidence) for channels 1, 2 and 3, respectively. The extended SSU temperature changes are in good agreement with those from the Microwave Limb Sounder (MLS) on the Aura satellite, with both exhibiting a cooling trend of ~ 0.6 ± 0.3 K decade −1 in the upper stratosphere from 2004 to 2012. The extended SSU record is found to be in agreement with high-top coupled atmosphere–ocean models over the 1980–2012 period, including the continued cooling over the first decade of the 21st century.

Description: On the potential for regolith control of fluvial terrace formation in semi-arid escarpments Earth Surface Dynamics Discussions, 3, 715-738, 2015 Author(s): K. P. Norton, F. Schlunegger, and C. Litty Cut-fill terraces occur throughout the western Andes where they have been associated with pluvial episodes on the Altiplano. The mechanism relating increased rainfall to sedimentation is however not well understood. Here, we apply a hillslope sediment model and reported cosmogenic nuclide concentrations in terraces to examine terrace formation in semi-arid escarpment environments. We focus on the Rio Pisco system in western Peru in order to determine probable hillslope processes and sediment transport conditions during phases of terrace formation. Specifically, we model steady state and transient hillslope responses to increased precipitation rates. The measured terrace distribution and reconstructed sediment loads measured for the Rio Pisco agree with the transient model predictions, suggesting strong climatic control on the cut-fill sequences in western Peru primarily through large variations in sediment load. Our model suggests that the ultimate control for these terraces is the availability of sediment on the hillslopes with hillslope stripping supplying large sediment loads early in wet periods. At the Rio Pisco, this is manifest as an approximately 4 × increase in erosion rates during pluvial periods. We suggest that this mechanism may also control terrace occurrence in other semi-arid escarpment settings.

Description: Post-Chernobyl surveys of radiocaesium in soil, vegetation, wildlife and fungi in Great Britain Earth System Science Data, 7, 215-221, 2015 Author(s): J. S. Chaplow, N. A. Beresford, and C. L. Barnett The data set "Post Chernobyl surveys of radiocaesium in soil, vegetation, wildlife and fungi in Great Britain" was developed to enable data collected by the Natural Environment Research Council after the Chernobyl accident to be made publicly available. Data for samples collected between May 1986 (immediately after Chernobyl) to spring 1997 are presented. Additional data to radiocaesium concentrations are presented where available. The data have value in trying to assess the contribution of new sources of radiocaesium in the environment, providing baseline data for future planned releases and to aid the development and testing of models. The data are freely available for non-commercial use under Open Government Licence terms and conditions. doi:10.5285/d0a6a8bf-68f0-4935-8b43-4e597c3bf251 . Supporting information to assist with the reuse of this data is available from the Environmental Information Data Centre (EIDC) ( http://eidc.ceh.ac.uk/ ).

Description: Three-dimensional finite element model of the drilling process used for fixation of Salter–Harris type-3 fractures by using a K-wire Mechanical Sciences, 6, 147-154, 2015 Author(s): A. Gok, K. Gok, and M. B. Bilgin In this study, the drilling process was performed with Kirschner wire (K-wire) for stabilization after reduction of Salter–Harris (SH) type-3 epiphyseal fractures of distal femur. The study was investigated both experimentally and numerically. The numerical analyses were performed with finite element method (FEM), using DEFORM-3D software. Some conditions such as friction, material model and load and boundary must be identified exactly while using FEM. At the same time, an analytic model and software were developed, which calculate the process parameters such as drilling power and thrust power, heat transfer coefficients and friction coefficient between tool–chip interface in order to identify the temperature distributions occurring in the K-wire and bone model (Keklikoǧlu Plastik San.) material during the drilling process. Experimental results and analysis results have been found as consistent with each other. The main cutting force, thrust force, bone model temperature and K-wire temperature were measured as 80° N, 120° N, 69 °C and 61 °C for 400 rpm in experimental studies. The main cutting force, thrust force, bone model temperature and K-wire temperature were measured as 65° N, 87° N, 91 °C and 82 °C for 800 rpm in experimental studies. The main cutting force, thrust force, bone model temperature and K-wire temperature were measured as 85° N, 127° N, 72 °C and 67 °C for 400 rpm in analysis studies. The main cutting force, thrust force, bone model temperature and K-wire temperature were measured as 69° N, 98° N, 83 °C and 76 °C for 800 rpm in analysis studies. A good consistency was obtained between experimental results and finite element analysis (FEA) results. This proved the validity of the software and finite element model. Thus, this model can be used reliably in such drilling processes.

Description: Pre-earthquake magnetic pulses Natural Hazards and Earth System Science, 15, 1873-1880, 2015 Author(s): J. Scoville, J. Heraud, and F. Freund A semiconductor model of rocks is shown to describe unipolar magnetic pulses, a phenomenon that has been observed prior to earthquakes. These pulses are suspected to be generated deep in the Earth's crust, in and around the hypocentral volume, days or even weeks before earthquakes. Their extremely long wavelength allows them to pass through kilometers of rock. Interestingly, when the sources of these pulses are triangulated, the locations coincide with the epicenters of future earthquakes. We couple a drift-diffusion semiconductor model to a magnetic field in order to describe the electromagnetic effects associated with electrical currents flowing within rocks. The resulting system of equations is solved numerically and it is seen that a volume of rock may act as a diode that produces transient currents when it switches bias. These unidirectional currents are expected to produce transient unipolar magnetic pulses similar in form, amplitude, and duration to those observed before earthquakes, and this suggests that the pulses could be the result of geophysical semiconductor processes.

Description: Ash iron mobilization through physicochemical processing in volcanic eruption plumes: a numerical modeling approach Atmospheric Chemistry and Physics, 15, 9361-9379, 2015 Author(s): G. A. Hoshyaripour, M. Hort, and B. Langmann It has been shown that volcanic ash fertilizes the Fe-limited areas of the surface ocean through releasing soluble iron. As ash iron is mostly insoluble upon the eruption, it is hypothesized that heterogeneous in-plume and in-cloud processing of the ash promote the iron solubilization. Direct evidences concerning such processes are, however, lacking. In this study, a 1-D numerical model is developed to simulate the physicochemical interactions of the gas–ash–aerosol in volcanic eruption plumes focusing on the iron mobilization processes at temperatures between 600 and 0 °C. Results show that sulfuric acid and water vapor condense at ~ 150 and ~ 50 °C on the ash surface, respectively. This liquid phase then efficiently scavenges the surrounding gases (〉 95 % of HCl, 3–20 % of SO 2 and 12–62 % of HF) forming an extremely acidic coating at the ash surface. The low pH conditions of the aqueous film promote acid-mediated dissolution of the Fe-bearing phases present in the ash material. We estimate that 0.1–33 % of the total iron available at the ash surface is dissolved in the aqueous phase before the freezing point is reached. The efficiency of dissolution is controlled by the halogen content of the erupted gas as well as the mineralogy of the iron at ash surface: elevated halogen concentrations and presence of Fe 2+ -carrying phases lead to the highest dissolution efficiency. Findings of this study are in agreement with the data obtained through leaching experiments.